CC BY
46
УДК 539.421.2
МОДЕЛИРОВАНИЕ РОСТА ТРЕЩИНЫ В ЗЕРНИСТОМ МАТЕРИАЛЕ
Е. О. Каракулина
Оренбургский государственный университет,
460018, Оренбург, пр-т Победы, 13.
E-mail: elok2004@yandex.ru
Разработана модель, описывающая процесс разрушения зернистых материалов. Установлены закономерности формирования и развития трещин в дискретных средах. Представлены результаты численного 'расчёта коэффициентов интенсивности напряжений первого рода с учётом упругих характеристик межзё-ренного заполнителя.
Ключевые слова: трещина, коэффициент интенсивности напряжений первого рода, зернистый материал.
Введение. Горные породы представляют собой хрупкие, неоднородные зернистые материалы, которые отличаются по форме, размеру зёрен, а также прочности связывающего вещества. Одна из первых работ по моделированию горных пород в виде зернистого материала принадлежит Троллопу. Элемент зернистой структуры представлен им в виде произвольного набора зёрен круглой формы.
В упрощённой модели Троллопа участки контакта между зёрнами изображают в виде отрезков [1]. При этом зерно материала представляет собой гексагональный элемент. Границы зёрен занимают большую объёмную долю структуры материала и существенно влияют на процесс его разрушения при механической нагрузке. Одним из параметров, характеризующих упругое состояние тела, является модуль Юнга. Горные породы с разной степенью уплотнённости имеют разные модули упругости. Наибольшее влияние на модуль Юнга оказывают минералогический состав, структура, условия залегания породы и другое. Однако критериев роста трещин, учитывающих упругие характеристики материала, до настоящего времени получено не было.
Целью данной работы является разработка математической модели роста трещин в зернистом материале и установление зависимости коэффициента интенсивности напряжений первого рода от значения модуля Юнга межзё-ренного заполнителя.
1. Постановка задачи и метод решения. Гассматривался модельный зернистый материал, содержащий плоскую трещину АВ длиной 21 = 0,0025 м. Структура материала была представлена в виде сетки гексагональных зёрен (рис. 1). Диаметр зерна исследуемого материала d = 1 мм. Контакты между зёрнами заполнены линейно-упругим материалом с упругими характеристиками, отличающимися от характеристик зёрен. Упругие характеристики при моделировании полагались следующими: для зерна модуль Юнга -Езер = 3 • 104 МПа; для межзёренного заполнителя Езап = 3-10; 3 • 102; 3 • 103; 3 • 104; 3 • 105; 3 • 106 МПа. Коэффициент Пуассона для зерна и заполнителя
Елена Олеговна Каракулина, старший преподаватель, каф. математического анализа.
V = 0,3. Ширина заполнителя между зёрнами породы /г = 0,1 мм. Жёсткость заполнителя определялась следующими равенствами:
Кп = Езап/1г, К3 = С/1г.
Численный расчёт выполнялся методом граничных элементов [2]. Граничные интегральные уравнения этого метода в процессе дискретизации границ представляются в виде следующей системы линейных уравнений:
N
3 = 1 N
(1)
&п — т І)ці I ^^[^ns^jDsj I BnnijDnj},
3 = 1
где г = 1,2,..., ТУ; N — количество всех граничных элементов; <т*, игп — касательные и нормальные напряжения на границе тела соответственно; Кщ, К3г — нормальная и касательная жёсткость заполнителя г-того элемента соответственно; Ду = ип(х, 0_) —ип(х, 0+), Dsj = щ(х, 0_) —ив(х, 0+) — компоненты разрывов смещений ^’-того отрезка трещины в нормальном и касательном направлениях соответственно; ип(х, 0_), ип(х, 0+)—нормальные смещения нижнего и верхнего берегов трещины; и3(х, 0_), и3(х, 0+)—касательные смещения нижнего И верхнего берегов трещины; В33у, Вапц, Впзу, Вппц — коэффициенты влияния.
При численном расчёте использовались следующие краевые условия на трещине АВ: ав = 0, ап = <т.°, где а3, (тп — соответственно касательные и нормальные напряжения; <т.°—задаваемое нормальное напряжение. Для вычисления коэффициента интенсивности напряжений первого рода использовали критерий [3]:
1<1 = -\/тг • _Л • Пт , (2)
2 2(1 — и) х уі у/1 — х’
где 21— длина трещины АВ, м; А/ — коэффициент интенсивности напряжений первого рода, Н/м3/2; О = Е/(2(1 + и)) — модуль сдвига, Н/м2.
2. Результаты моделирования. В результате решения системы линейных уравнений (1) были найдены значения величины Оп^. На рис. 2 показаны результаты численного решения задачи для сплошной и дискретной среды с различными упругими характеристиками межзёренного заполнителя. Так
Рис. 2. Величина раскрытия трещины в сплошном (штриховая линия) и зернистом (сплошные линии) материалах: 1 — Езап = 3 • 10Б МПа, 2 —Езап = 3 • 104 МПа, 3 —Езап = 3 • 103 МПа
как нормальные смещения верхнего и нижнего берегов трещины равны по величине, то на рисунке показана величина смещения только нижнего берега трещины, где у = 0 — линия симметрии. Из рисунка видно, что уменьшение упругих характеристик межзёренного заполнителя приводит к увеличению разрывов смещений и к отличию формы трещины от эллиптического вида.
Для раскрытия берегов трещины в сплошной среде используют представление
Оп = -2^с и\л/12 -х2.
Для разрывов смещений в зернистом материале были найдены аппроксимирующие функции вида Бп = /(ж) • л/1 — х, где /(ж) = ао + а\Х + агж2 + + • • • + а,кХк при к = 6. Коэффициенты сц (г = 0,1,... , 6) найдены методом наименьших квадратов. В таблице показана относительная погрешность значений аппроксимирующей функции в сравнении с численным расчётом Оп, которая подтверждает достоверность и точность полученных результатов.
После подстановки Ип в критерий (2) вычислялись значения коэффициентов интенсивности напряжений первого рода. В результате расчётов была установлена зависимость коэффициента интенсивности напряжений первого рода от величины модуля Юнга межзёренного заполнителя:
К/ = рл/тй ■ у , (3)
\ -^зап /
где р— нормальные напряжения, приложенные к берегам трещины, Па;
V (= 1,212 . 1еЕ (+ 0,658 .1е‘ (§=) - 4,427 ■ 1е3 +
V -С^зя/гт / V -С^зя/гт / V -С^зя/гт / V -С^зя/гт /
+0,418 ■ 182 (0) + 4,569 • 1к (1У + 1.76.
На рис. 3 показана зависимость коэффициента интенсивности напряжений первого рода от длины трещины при различных значения модуля Юнга
межзеренного заполнителя при фиксированном значении р = 1 КПа в формуле (3). Полученные зависимости показывают, что значение коэффициента интенсивности напряжений первого рода увеличивается при уменьшении модуля упругости межзеренного заполнителя.
К!, Н/м3/2
Рис. 3. Зависимость коэффициента интенсивности напряжений первого рода от длины трещины в сплошном (штриховая линия) и зернистом (сплошные линии) материалах: 1 — Езап = 3 • 104 МПа, 2 — Езап = 3 • 103 МПа, 3 — Е3ап = 3 • 102 МПа
3. Выводы. Анализируя полученные результаты и зависимости, можно сделать следующие выводы:
- уменьшение упругих характеристик межзёренного заполнителя приводит к увеличению разрывов смещений берегов трещины и к отличию её формы от эллиптического вида;
- прочностные характеристики материала прямо пропорциональны значениям модуля Юнга межзёренного заполнителя;
- полученное аналитическое выражение для нахождения коэффициента интенсивности напряжений первого рода позволяет прогнозировать рост трещин в дискретных средах с различным модулем Юнга межзёренного заполнителя.
Относительная погрешность вычисления разрывов смещений
Язап = 30 МПа Ез&п = 3 • 106 МПа
В„ ■ 104,м В„ ■ 104,м 5фп) Вп ■ 106,м В„ ■ 106,м 5фп)
-1,03 -1,04 0,77% -0,84 -0,84 0,24%
-1,04 -0,99 4,58% -1,19 -1,19 0,54%
-1,05 -1,16 10,76% -1,41 -1,41 0,16%
-1,05 -0,93 11,59% -1,53 -1,53 0,33%
-0,66 -0,68 4,33% -1,59 -1,58 0,47%
-0,66 -0,73 11,25% -1,59 -1,59 0,07%
-1,05 -0,96 8,19% -1,53 -1,53 0,02%
-1,05 -1,09 3,91% -1,41 -1,4 0,54%
-1,04 -1,03 0,84% -1,19 -1,19 0,06%
-1,03 -1,03 0,05% -0,84 -0,84 0,21%
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Герике Б. Л., Полкунов Ю. Г., Герике П. Б. Математические модели циклического разрушения крепких горных пород дисковым инструментом. Кемерово: Кузбассвузиз-дат, 2001. 171 с. [Gerike B.L., Polkunov Yu. G., Gerike P. В. Matematicheskie modeli tsiklicheskogo razrusheniya krepkikh gornykh porod diskovym instrumentom. Kemerovo: Kuzbassvuzizdat, 2001. 171 pp.]
2. Crouch S. L., Starfield A. M. Boundary element methods in solid mechanics: with applications in rock mechanics and geological engineering. London - Boston - Sydney: George Allen & Unwin, 1983. 322 pp.; русск. пер.: Крауч С., Старфилд А. Методы граничных элементов в механике твердого тела. М.: Мир, 1987. 328 с.
3. Линьков А. М. Комплексный метод граничных интегральных уравнений теории упругости. СПб.: Наука, 1999. 382 с.; англ. пер.: Linkov А. М. Boundary integral equations in elasticity theory / Solid mechanics and its applications. Vol. 99. Dordrecht, Boston: Kluwer Academic Publishers, 2002. 268 pp.
Поступила в редакцию 20/XI/2010; в окончательном варианте — 12/11/2011.
MSC: 74R10; 74S15
SIMULATION OF CRACK GROWTH IN A GRANULAR MATERIAL
E. O. Karakulina
Orenburg State University,
13, pr-t Pobedy, Orenburg, 460018, Russia.
E-mail: elok2004@yandex.ru
A model describing the process of destruction of granular materials was developed. The regularities of formation and development of cracks in discrete environment were determined. The results of numerical calculation of stress intensity factors of the first kind, taking into account the elastic characteristics of grain boundary core are presented.
Key words: crack, stress intensity factor of the first kind, granular material.
Original article submitted 20/XI/2010; revision submitted 12/11/2011.
Elena O. Karakulina, Lecturer, Dept, of Mathematical Analysis.
