CC BY
32
УДК 624.121:550.3
А.П.ГОСПОДАРИКОВ, д-р техн. наук, профессор, kafmat@spmi. ru М.В.МАКСИМЕНКО, аспирант, [email protected],
Национальный минерально-сырьевой университет «Горный», Санкт-Петербург
A.P.GOSPODARIKOV, Dr. in eng. sc., professor, [email protected] M.V.MAKSIMENKO,post-graduate student, [email protected], National Mineral Resources University (Mining University), Saint Petersburg
ОБ ОДНОМ ПОДХОДЕ К ИССЛЕДОВАНИЮ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ МАССИВА ГОРНЫХ ПОРОД С УЧЕТОМ НЕЛИНЕЙНОГО ХАРАКТЕРА ПРОЦЕССА
ИХ ДЕФОРМИРОВАНИЯ
Рассмотрена нелинейная краевая задача об исследовании напряженно-деформированного состояния (НДС) массива горных пород вокруг одиночной горизонтальной протяженной выработки. Предложен алгоритм определения основных параметров НДС горного массива, включающий комплекс вычислительных методов.
Ключевые слова: физическая нелинейность, метод линеаризации Ньютона -Канторовича - Рафсона.
THE METHOD OF ANALYSIS STRESS-AND-STRAIN STATE OF A ROCK MASSIF WITH REGARD TO NONLINEAR CHARACTER
OF STRAINING PROCESS
The nonlinear problem of analysis of the stress-and-strain state of rock massif is discussed. The numerical algorithm analysis of the stress-strain state of a massif, which includes a set of computational methods, is suggested.
Key words: physical nonlinearity, Newton - Kontorovich - Raphson linearization method.
В связи с переходом подземных горных работ в основных угольных бассейнах страны на более глубокие горизонты значительно ухудшаются горно-геологические условия разработки пластов. На глубоких шахтах существенно увеличиваются расходы на ремонт, перекрепление подготовительных выработок, усложняются мероприятия по борьбе с опасными проявлениями горного давления, что в конечном итоге приводит к значительному росту смещений горных пород вокруг выработок [2]. Поэтому традиционные методы охраны объектов от вредных влияний горных работ в современных условиях все чаще становятся или неэффективными, или вовсе неприемлемыми. С практи-
60
ческой точки зрения очень важно иметь представление о реальном распределении компонент тензоров напряжений и деформаций, а также вектора перемещений в массивах пород при наличии горных выработок различного назначения и очертания, целиков, с учетом выработанного пространства и т.д.
Получение аналитических решений в таких случаях ограничено, как правило, выбором простой физической модели среды (сплошная, изотропная и однородная) в рамках применимости линейного закона Гу-ка и формой выработки (круглая). Математическое моделирование неоднородного породного массива, ослабленного одной или несколькими подземными выработками
ISSN 0135-3500. Записки Горного института. Т.205
сложного очертания, а также учет нелинейного характера процесса деформирования горных пород предполагают использование эффективных численных методов.
Проведение выработки с физической точки зрения представляет собой образование некоторой пространственной полости в массиве горных пород. Анализ геометрических параметров рассматриваемого объекта приводит к тому, что граничные поверхности последнего представляют собой многосвязные области, учет которых существенно усложняет процедуры нахождения решения рассматриваемой задачи. Установлено, что для выработок, располагающихся на достаточно большой глубине Н, погрешность от замены реальных граничных условий (земная поверхность) на задачу с односвязным контуром (отверстие в тяжелой плоскости) быстро убывает с увеличением глубины залегания последних [1].
В работе принята физическая модель нелинейно-упругой среды (нелинейная связь между компонентами тензоров напряжений и деформаций). Начальное поле напряжений в массиве горных пород формируется под действием гравитационных сил.
Отметим также, что для случая одиночной протяженной горизонтальной выработки, у которой длина во много раз превышает два других характерных размера (высоту и ширину), пространственная задача по определению основных параметров НДС массива горных пород может быть сведена к плоской.
Учитывая, что фактическое влияние выработки на породный массив является локальным и что глубина заложения выработки Н превышает десять радиусов выработки и более, собственный вес горных пород имитируется приложением распределенной нагрузки интенсивностью уН к границам исследуемой невесомой области.
Таким образом, моделирование НДС массива горных пород, включающих выработку, в рамках плоской деформации сводится к следующему: рассматривается система разрешающих дифференциальных уравнений в частных производных, описывающих физически нелинейную деформацию массива, дополненную надлежащими
граничными условиями. Рассматриваемую двумерную нелинейную краевую задачу запишем в единой векторной форме, удобной для изложения разработанного алгоритма численного решения. Введем в рассмотрение вектор-функции:
с =
а,
Vх чу )
; в =
V8 чу )
; и =
(и\
V У )
Нелинейная краевая двухмерная задача содержит восемь уравнений: два уравнения равновесия, три уравнения Коши (соотношения между компонентами тензора деформаций и компонентами вектора перемещений), а также три уравнения, устанавливающие соотношения между компонентами тензоров напряжений и деформаций (в рамках принятого нелинейного физического закона). Тогда матричная запись разрешающих уравнений примет вид
Ат с = 0; в = Аи; с = Св,
- дифференциальный оператор,
матрица С определяется из
нелинейных физических соотношений между компонентами тензоров напряжений и деформаций на основе кривой деформирования о - е. В частности, в работе используется закон, предложенный Бахом, и связь между напряжениями и деформациями представлена в виде степенной зависимости оi - Вг ва, где В г (' = 1,2,3...) - коэффициенты деформирования, Па; аг < 1 получены экспериментальным путем [1].
Граничные условия, заданные в перемещениях на контуре выработки L и на границе области Г, имеют вид
где А -
дх 0
А = 0 д ду
д д
кдУ дх
где иь =
V )
V УГ)
ч
в
и = иь , и =иГ
; иг =
Итак, исследование НДС массива горных пород, включающего выработку, сводится к решению системы нелинейных дифференциальных уравнений (1) в частных производных при одновременном удовлетворении граничных условий (2) на контуре выработки и на границе рассматриваемой области.
В области П = X х Y,X е [0;a], Y е [0;Ь] ищутся вектор-функции основных
V (х, y) =
С,
ху
u
V V j fr Л
и
дополнительных
Г (x, y) =
Vy xv j
неизвестных, которые всюду
в области О удовлетворяют системе т дифференциальных уравнений в частных производных
IV = £(V,¥,х,у),х е (0,а),у е (0,Ь), (3) где А - дифференциальный оператор
(д_ дх
0
A =
- 0 д 0
ду д
ду дх
— — 0
0 0 0
д дх
0
0 0 0
ду д
g (V ,W, х, у) =
д_ _ ду дх
( 0 ^ 0 S х
s„
Vх ху j
и ассоциированной с (3) системе р трансцендентных уравнений [3]
f (V ,W, х, у) = 0,
(4)
где f (V, W, х, у) =
4s"1 -С х^ Bsa/ -Су
VB3У" -хху j
Граничные условия задаются некоторой системой уравнений общего вида
D (Vl V) = 0.
(5)
Отметим, что в рассматриваемой постановке исходной задачи т = 5, р = 3.
Рассмотрим нелинейное операторное уравнение, соответствующее системам (3)-(5) Т (ю) = 0. Под ю понимается вектор-функция
Ю(х, у) = [V (х, у), ¥ (х, у)]т , компоненты которой представляют собой непрерывные вместе с первыми частными производными ограниченные функции, заданные в области О.
Оператор Т(ю), символизирующий дифференциальный оператор А (3) и левые части уравнений (4) и (5), является дифференцируемым по Фреше нелинейным оператором в банаховом пространстве Е.
Используя метод линеаризации Ньютона - Канторовича - Рафсона, представим операторное уравнение в виде
T (ю)ю + T2 (ю) = 0; T (ю) = тю (ю), T2 (ю) = T (ю) - тю (ю)ю,
(6)
где ТЮ - производная Фреше нелинейного оператора Т(ю).
Вводя индексы итерации п,п - 1 в уравнение метода линеаризации (6), приходим к соотношению при каждом п = 1,2,...:
Т; (ю п-1 )ю п + Т(юп-1) - Т; (ю п-1 )юп-1 = 0.
Применяя метод линеаризации Ньютона - Канторовича - Рафсона к краевой задаче (3)-(5), получим систему
AV = G1(V ,W, х, y)V + G2(V ,W, х, у )W + + g1(V ,W, х, у);
Fj (V, W, х, у )V + F2 (V, W, х, y)W + + fj(V ,W, х, у) = 0;
D (Vl V) = 0.
х
s
62 -
ISSN 0135-3500. Записки Горного института. Т.205
Отметим, что, как правило, граничные условия, записанные в перемещениях, на контуре выработки L и на границе области Г описываются системой линейных алгебраических уравнений вида (5).
Окончательно итерационный процесс метода линеаризации Ньютона - Канторовича - Рафсона запишется в виде
AV" = G1"-1V" + G2"-1W" + g"-1;
F1"-1V" + F2"-1W" + J"'1 = 0; D (VL V) = 0.
(7)
В системе (7) G"-1, G2"
рл-1 ТГп~1 — М ' 2
матрицы Якоби от вектор-функций £ (V ,Ж, х, у) и / (V ,Ж, х, у) по тому векторному аргументу, который стоит множителем при данной матрице. Размерности этих матриц однозначно определяются размерностью искомых и заданных вектор-функций, присутствующих в уравнениях рассматриваемой краевой задачи.
Здесь и далее использованы сокращенные обозначения вида:
G1 = G1(V"- 1, W"-1 , x, y);
G2 = G2(V"- -\W"- 1, x, y);
F"-1 = F1(V"- 1, W"-1 , x, y);
F"-1 = F1(V"- 1, W"-1 , x, y);
g"-1 = g1(V"-1 l,W"-1 , x, y);
T-"-1 J1 = J1(V"-1 ,W"-1 , x, y).
Векторы g" 1, f" 1, согласно заданным системам уравнений (3)-(4), определяются по формулам
jg"-1 = g"-1 -G" -1V"-1 -G2"- W"-1; I f"-1 = f 1 - F"- 1V"-1 - F"-1W"-1.
Таким образом, последовательность ш„ на каждом шаге итерации " определяется как решение линейной краевой задачи, а с учетом необходимого для численной реализации перехода к соответствующей дискретной задаче (например, на основе метода конечных элементов) - из решения системы линейных алгебраических уравнений.
ЛИТЕРАТУРА
1. Вялов С.С. Реологические основы механики грунтов: Учеб. пособие. М., 1978. 447 с.
2. Геомеханика: Учебник для вузов / И.В.Баклашов, Б.А.Картозия, А.Н.Шашенко, В.Н.Борисов. М., 2004. Т.2. 249 с.
3. Господариков А.П. Метод расчета нелинейных задач механики горных пород при подземной разработке пластовых месторождений. СПб, 1999. 129 с.
REFERENCE
1. Vyalov S.S. Rheological fundamentals of soil mechanics: Training manual for building high schools. Moscow, 1978. 447 с.
2. Geomechanics: Textbook for universities / I.V.Baklashov, B.A.Kartozia, A.N.Shashenko, V.N.Borisov. Moscow, 2004. Vol.2. 249 с.
3. Gospodarikov A.P. Method of calculation of nonlinear problems of rock mechanics in underground mining of bedded deposits. Saint Petersburg, 1999. 129 с.
