CC BY
11
УДК 539.3
DOI: 10.33764/2618-981Х-2019-2-5-89-97
ХРУПКОЕ РАЗРУШЕНИЕ НЕСЖИМАЕМОГО МАТЕРИАЛА ПРИ ПЛОСКОЙ ДЕФОРМАЦИИ
Анвар Исмагилович Чанышев
Институт горного дела им. Н. А. Чинакала СО РАН, 630091, Россия, г. Новосибирск, Красный пр., 54, доктор физико-математических наук, главный научный сотрудник, тел. (383)335-97-50; Новосибирский государственный университет экономики и управления, 630099, Россия, г. Новосибирск, Каменская, 52, зав. кафедрой математики и естественных наук, e-mail: [email protected]
Ольга Евгеньевна Белоусова
Институт горного дела им. Н. А. Чинакала СО РАН, 630091, Россия, г. ный пр., 54, кандидат технических наук, старший научный сотрудник, e-mail: [email protected]
Ольга Анваровна Лукьяшко
Институт горного дела им. Н. А. Чинакала СО РАН, 630091, Россия, г. Новосибирск, Красный пр., 54, инженер, тел. (383)335-97-50, e-mail: [email protected]
Исследуется напряженно-деформируемое состояние массива пород при плоской деформации в случае его объемной несжимаемости и идеально-хрупком разрушении. Получены характеристики системы дифференциальных уравнений равновесия, соотношения на них. При этом использовалось условие соосности тензоров напряжений и деформаций. Сформулирована граничная задача для определения напряженно-деформированного состояния в зоне разрушения. В качестве примера рассмотрены уравнения идеально-хрупкого запредельного деформирования массива пород в виде прямоугольной пластины (целика) при однородном сжатии. Показано, что при слежении за смещением боковой границы пластины возможно предсказать наступление ее разрушения.
Ключевые слова: хрупкое разрушение, напряженное состояние, характеристики, соотношения на характеристиках, граничная задача.
BRITTLE DESTRUCTION OF INCOMPRESSIBLE MATERIAL DURING FLAT STRAIN
Anvar I. Chanyshev
Chinakal Institute of Mining SB RAS, 54, Krasny Prospect St., Novosibirsk, 630091, Russia, D. Sc. Chief Researcher; Novosibirsk State University of Economics and Business Administration, 52, Kamenskaya St., Novosibirsk, 630099, Russia, Head of Department of Mathematics and Natural Sciences, phone: (383)335-97-50, e-mail: [email protected]
Olga E. Belousova
Chinakal Institute of Mining SB RAS, 54, Krasny Prospect St., Novosibirsk, 630091, Russia, Ph. D., Senior Researcher, phone: (383)335-97-50, e-mail: e-mail: [email protected]
Olga A. Lukyashko
Chinakal Institute of Mining SB RAS, 54, Krasny Prospect St., Novosibirsk, 630091, Russia, Engineer, phone: (383)335-97-50, e-mail: e-mail: [email protected]
Новосибирск, Крас-тел. (383)335-97-50,
In the paper stress-strain behavior of solid during flat strain in case of its volumetric incompressible behavior and ideally brittle destruction is studied. Parameters of the system of differential equations of balance and its correlations are obtained. In this case, condition of stress and strain tensors axiality is used. Boundary problem for determination of stress-strain behavior at destruction zone is formulated. As example, equations of ideally brittle out-of-limit deformation of solid in form of rectangular plate (pillar) during uniform compression are considered. It is shown that when displacement of side border of the plate is observed, it is possible to predict its destructions.
Key words: brittle destruction, stress behavior, parameters, correlations of the parameters, boundary problem.
Введение
Диаграмма изменения напряжений с ростом деформаций характеризует изменение сопротивления материала деформированию [1-4]. Если напряжение возрастает, то это означает рост сопротивления, если напряжение падает, то сопротивление материала деформированию падает, то есть происходит его разрушение. Кривая запредельного деформирования может быть пологой. Если она представляет собой отрезок вертикальной прямой, то этот факт относится к идеально-хрупкому разрушению [5-11]. Плоская деформация идеально-хрупкого материала с объемной несжимаемостью исследовалась в [8]. Здесь использовалось условие несжимаемости среды и два уравнения совместности деформаций, содержащие производные от компоненты вектора поворота в направлении оси 2. Были получены две характеристики системы дифференциальных уравнений для перемещений и два соотношения на них, связывающие компоненту вектора поворота и угол, задающий главные оси тензора деформаций. В [8] не сформулирована граничная задача, и, кроме того, здесь нет уравнений для определения напряжений. В предлагаемой работе будут рассмотрены вопросы, которые не решались в [8]: будут сформулированы уравнения для определения напряжений и граничная задача для нахождения и напряжений, и деформаций.
Часть 1. Определение напряжений, деформаций при плоской деформации для объемно - несжимаемого идеально-хрупкого разрушаемого материала
Постановка задачи. Пусть рассматривается среда, для которой в прямоугольной декартовой системе координат хОуг выполняются уравнения:
= 0 % + = ^ С1)
Sx S y
-= Tcos2Q, sXy =Tsin2Q, r=Tp, (2)
где Г p - значение максимальной касательной деформации
г
г _ ^
8х 8у
+ 8
ХУ
при идеально хрупком разрушении на диаграмме Т = Т(Г), где Т - максимальное касательное напряжение (рис. 1).
Рис. 1. Диаграмма изменения максимального касательного напряжения Т с ростом максимальной касательной деформации (сдвига) Г
Полагаем, что в процессе деформирования и разрушения среды
tg 2О =
28
ху _
tg 20 =
2т
ху
(3)
х
у
х
у
где 0 - угол между первым главным направлением для тензора напряжений Та
и осью х, О - угол между первым главным направлением для тензора деформаций Т8 и осью х. Равенство (3) означает то, что главные оси тензоров Та
и Т8 совпадают во все время деформирования среды.
Кроме (1)—(3) предлагаем выполнение условий совместности деформаций:
дву двху да2
дх ду ду двх д8ху даг
0,
ду дх
+ ■
дх
= 0,
где со2 - компонента вектора поворота
®2 =т
г ду дил дх ду
(5)
и, у,ю- компоненты вектора перемещений. А также предполагаем выполнение уравнений равновесия:
'х. + = о,
дтху да
-ху
дх ду
+ ■
у
дх ду
= 0.
(6)
В [8] соотношения (1), (2) подставлялись в (4). В результате подстановки получались соотношения
2Г
Р
дП дП Б1п2П--еов2П-
2Г
Р
дх
в1п2ПдП + соб2п
ду
ду
ал
дх
дю2
ду
д®2 дх
(7)
Система (7) - гиперболического типа, ее характеристики
£=%.
оик оик
(8)
Они, как видно из (8), совпадают по направлению с направлениями главных осей тензора деформаций Т^ . Соотношения на характеристиках имеют вид [5]:
со2 -П = % со2 + П = г),
(9)
где константы.
Для решения задач необходимо знать граничные значения функций ю2 = ю2 (х, у), П = П(х, у). С этой целью будем считать, что на границе тела Ь заданы смещения
и = и (х, у), V = V (х, у), (10)
где координаты х,у е Ь. Вычисляя полные дифференциалы от и и у, имеем
ди 7 ди 7 7 ду 7 ду ,
аи = — ах л--ау, ау = — ах л--ау
дх ду дх ду
или
аи = £хОх + (ёху - ао2)ау, ау = (ёху + со2)ах + буйу
(11)
Подставляя в (11) выражения (1), (2), находим
du = Гр cos 2Q.dx + (Г sin 2Q- az ) dy, dv = (гp sin 2Q + wz )x - Гp cos2Qdy.
Поделив (12) на величину ds =
2 2 + (dy) , имеем соотношения:
u's = -Гp cos 2Q sin (р + Гp sin 2Qcos р - ®z cos р, v's =-Гp sin2Qsinр-Гp cos2Qcosр- ®z sinр,
(13)
где , - производные от смещений и, упо касательной к Ь, а р- угол между нормалью к границе Ь и осью х . Из (13) находим
cos
2(П- ) = (u's sin р - v's cos р) = u's cos (2Q-р( + v's sin (2Р-р) (14)
2( р)= , ®z = \ . (14)
^ cos2 (Q-р)
Г,
Формулы (14) позволяют определять граничные значения угла О и компоненты вектора поворота a)z. Дальше возникает задача: как определить напряженное состояние в области разрушения?
Для определения напряженного состояния имеем соотношения (3), то есть в = О, и общепринятые зависимости
<jx = < + Tcos20, <jy = <-Tcos20, Txv = Tsin20
ХУ
(15)
где угол в = в{х, у) из предыдущего анализа считается известной функцией координат х,у . Подставляя (15) в (6) при известном значении углав = в(х,у), получаем систему уравнений:
d< dT
dx dx
d< dT
dy dx
dT
cos20 +--sin 20 = 2T
dy
ат
sin 20--cos20 = -2T
dy
. 20d0 20d0 sin 20--cos20—
dx dy
20d0+ • 20d0 cos20--+ sin20—
dx dy
(16)
Здесь правые части (16) с точностью до Т считаются известными функциями. Сравнивая (16) с (7), видим, что правые части (16) соответственно равны следующим выражениям
daz T
daz T
dy Г
p
dx Г
p
Определяя характеристики (16), находим, что они совпадут с (8), где 0 = П. При этом соотношения на характеристиках принимают вид:
ау = tg0, аа+ат=— ах г
р
дю2 7 дю2 7 ах--— ау
ду
дх
— = -ctg0, аа-ат=— ах г
р
дю2 7 дю2 7 ■ах--— ау
(18)
ду
дх
Правые части (18) представляет собой векторное произведение векторов
graаa
г
^д®2 д®2 0 дх ду
и ая = (ах, ау, о)
Это означает, что при интегрировании соотношений на характеристиках, полученные интегралы будут зависеть от пути интегрирования, так как (18) через векторное произведение связаны с площадями параллелограммов, построенных на векторах graаaz, ая.
Возникает вопрос: как найти граничные значения функций а = а(х, у), Т = Т(х, у), входящих в (18)? Для этого имеем вектор напряжений Коши на границе тела Ь. Его координаты вдоль осей х,у имеют соответственно вид:
рП =(а + Т сов20)совр + Т в1п2081пр, рП = Тв1п20со8^ + (а-Тсо820)в1п^,
(19)
где р - угол между нормалью к Ь и осью х . Из (19) получаем граничные значения
у соор- рх с'пр рУ з1п(20- р(- рЩ соб(20 -р(
Т= рщ собр-рпБ1пр а= б1п2 (0-р)
б1П2 (0-р)
(20)
Таким образом, применяя (20) и используя (18), устанавливаем напряженное состояние в области разрушения.
Часть 2. Определение НДС пластины в состоянии запредельного деформирования при сжатии
Данная задача имеет отношение к устойчивости целиков при разработке полезных ископаемых. Она исследовалась во многих работах [12-20].
Как пример приведем уравнения и простейшее решение задачи о запредельном деформировании пластины при плоской деформации при однородном сжатии. Пусть в процессе нагружения выполняются условия:
8ху = тху = (21)
Это означает, что углы О = в = 0, то есть первые главные направления для тензоров Т^ и совпадают с направлением оси х. Далее, пусть в достигнутом напряженно-деформированном состоянии
£х = Гр, 8у = -Гр,
где Г р - предел прочности (см. рис. 1).
Отсюда находим смещения и, V при условиях 8ху = az = 0 :
и = Грх, V = -Гру.
Напряжение удовлетворяет уравнениям равновесия:
==0.
дх ду
Тогда
^у = / ( х ), ^х =р( у ),
где /, g - заданные функции координат.
Данная ситуация представлена на рис. 2.
Рис. 2. Пластина высоты H и ширины L сжата в вертикальном направлении усилием cjy = f (х). Характеристиками
системы дифференциальных уравнений служат линии х = const, y = const
Из предыдущих формул и рис. 2 не трудно заметить, что для контроля момента разрушения необходимо следить за изменением бокового смещения u .Как только оно достигнет критического значения u = Гр L, так тут же произойдет разрушение пластины.
Замечание
Из предыдущего анализа не видно: когда величина Г достигнет значения Гр ? Для решения этой задачи следует рассматривать предыдущие решения,
например, относящиеся к упругости (см. рис. 1). Для упругого состояния имеем
du = sxdx + ~ c°z )dy, dv = ^sxy + (Oz )dx + sydy .
Поделив это выражение наds = ij(dx)2 +(dy)2 , как и прежде, получаем
us = -% cos (p + (£xy -®z)sinp, vS =(sxy + ®z)cos p + Sy sinp. (22) Эти соотношения дополняются следующими:
+ Ox + °y Г T 2sxy 2Ty (23)
sx + sy =-^, Г = —, -=-¿—, (23)
' 2k 2JU Sx -Sy °x -°y
Pn = ax cos P + Txy sin p, = Txy cos P + Gy sin P, (24)
(рП, рУ - заданные величины на границе L). В результате имеем семь уравнений (22)-(24), для определения семи неизвестных sx, s , sy, а> , о , ст, т .
x xy y z x y xy
Система имеет единственное решение, из которого устанавливается величина Г.
Выводы
1. Определено напряженно-деформированное состояние в области запредельных деформаций при идеально хрупком разрушении объемно несжимаемой среды.
2. Рассмотрено решение задачи об однородном сжатии пластинки в условиях идеально хрупкого разрушения. Показано, что для контроля целика массива горных пород необходимо следить за изменением его смещения в поперечном направлении.
Работа выполнена при финансовой поддержке проекта РФФИ (№ 18-05-00757 А).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Курленя М.В., Афиногенов Ю.А., Жигалкин В.М., Опарин В.Н., Усольцева О.М., Ча-нышев А.И. Блочная феноменологическая механическая модель элемента деформируемой среды. Часть 2. Динамические эффекты// ФТПРПИ. - 1998. - № 5. - С. 34-49.
2. Чанышев А.И., Белоусова О.Е. Об одной интерпретации зональной дезинтеграции массива горных пород вокруг выработок // Физическая мезомеханика. - 2009. - Т. 12. -№ 1. - С. 89-99.
3. Феодосьев В.И. Сопротивление материалов: Учебник для вузов. - 10-е изд., перераб. и доп. - М.: Изд-во МГТУ им. Н. Э. Баумана, 1999. - Т. 2. - 592 с.
4. Старовойтов Э. И. Сопротивление материалов. - М.: Физматлит, 2008. -384с.
5. Никифоровский В.С., Шемякин Е.И. Динамическое разрушение твердых тел. -Новосибирск: Наука, 1979. - 272 с.
6. Шемякин Е.М. Задача о хрупком шарнире// Изд. АН.МТТ. - 1996. - № 2. -С. 138-144.
7. Леонов М.Я., Панасюк В.В. Развитие мельчайших трещин в твердом теле//Прикл. Механика. - 1959. - № 5. -С. 49 - 61.
8. Качанов Л.М. Основы механики разрушения. -М.: Наука, 1974. - 312 с.
9. Шемякин Е.М. О хрупком разрушении твердых тел (плоская деформация)//Известия Академии наук. Механика твердого тела. - 1997. -№ 2. -С. 145.
10. Черепанов Г.П. Механика хрупкого разрушения. - М.: Наука, Гл. ред. физ-мат литературы. - 1974. - 640 с.
11. Партон В.З., Морозов Е.М. Механика упругопластического разрушения 2-е изд., перераб. и доп. - М.: Наука, Гл. ред. физ-мат. Литературы.- 1985. - 504 с.
12. Борисов А. А. Механика горных пород и массивов. - М.: Недра. - 1980. - 360 с.
13. Проскуряков Н.М. Управление состоянием массива горных пород. - М.: Недра, 1991. - 368 с.
14. Шерман Д.И. К вопросу о напряженном состоянии междукамерных целиков. Упругая весомая среда, ослабленная двумя отверстиями эллиптической формы // Известия АН СССР. Отдел технических наук. - 1952. - № 6. - С. 840- 857.
15. Слесарев В. Д. Определение оптимальных размеров целиков различного назначения. - М.: Углетехиздат. - 1948. - 196 с.
16. Руппенейт К. В. Некоторые вопросы механики горных пород. - М.: Углетехиздат, 1954. - 384 с.
17. Нестеров М.П. Об инженерных методах расчета ленточных целиков // Горный журнал. -1968. - № 9.
18. Михлин С.Г. О напряжениях в породе над угольным пластом // Известия АН СССР. Отделение технических наук. - 1942. - № 7. - С. 13-28.
19. Динник А.Н. О давлении горных пород и расчет крепи круглой шахты // Инженерный работник. - 1925. - № 7. - С. 15-18.
20. Воробьев А.Н., Прутков С.Н. Метод расчета пролета камер, барьерных и междукамерных целиков на основе теории арочных систем // ГИАБ. - 2006. - № 10. - С. 201-205.
© А. И. Чанышев, О. Е. Белоусова, О. А. Лукьяшко, 2019
