Математическое моделирование зонального разрушения массива вокруг подземных выработок и реверсивного деформирования образцов горных пород

Электронное периодическое издание «Вестник Дальневосточного государственного технического университета» 2011 год № 3/4

25.00.00 Науки о земле

УДК 622.248.34:517.9

Ци Чэнчжи, Цянь Циху, Ван Минян

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ЗОНАЛЬНОГО РАЗРУШЕНИЯ МАССИВА ВОКРУГ ПОДЗЕМНЫХ ВЫРАБОТОК И РЕВЕРСИВНОГО ДЕФОРМИРОВАНИЯ ОБРАЗЦОВ ГОРНЫХ ПОРОД

Ци Чэнчжи - д.ф.м. н., профессор, заместитель декана стройтельного и коммуникационного факультета

Цянь Циху - академик Китайской инженерной академии

Ван Минян - профессор Инженерного Института

Когда горное давление превышает предел прочности пород на одноосное сжатие, вокруг выработок возникает явление зональной дезинтеграции. Экспериментально установлено, что эффект смены знака приращения деформаций (реверсивное деформирование) наблюдается в

92

предварительно нарушенных образцах горных пород в предразрушающей области нагружения при сжатии. Эксперименты показывают, что процесс деформирования горных пород может рассматриваться как многостадийный процесс фазового перехода. В данной работе эти явления исследуются на основе теории размытого фазового перехода. На основе разложения Гинзбурга-Ландау по свободной энергии относительно праметра порядка получено управляющее уравнение относительно параметра порядка. Показаны и обоснованы общая природа эффекта реверсивного деформирования образцов горных пород при сильном сжатии и явления зональной дезинтеграции вокруг выработок на большой глубине.

Ключевые слова: зональная дезинтеграция, сильно сжатое состояние, эффект реверсивного деформирования, размытый фазовый переход.

Mathematical modeling of zone failure of the rock mass around the underground openings and reversive deformation of rock samples. Qi Chengzhi (Beijing University of Civil Engineering and Architecture, Beijing, China), Qian Qihu (Chinese Engineering Academy, Beijing, China), Wang Mingyang (Nanjing University of Science and Technology, Nanjing, China).

When rock mass stresses exceeds strength of rock on uniaxial compression, around the openings there is a phenomenon of zone failure. It is experimentally determined, that the effect of change of a sign on an increment of deformations (reversive deformation) is taken place in preliminary broken samples of rocks in predestroying

93

loading area at compression. Experiments show, that process of deformation of rocks can be considered as multistage process of phase transition. In the given work these phenomena are investigated on the basis of the theory of fuzzy phase transition. On the basis of Ginzburg-Landau failure on free energy relative to order parameter the operating equation concerning order parameter is constructed. The general nature of the effect of reversive deformation of rocks samples at high stresses and the phenomena of zone failure round the openings on the big depth are shown and proved.

Key words: zone failure, high stresses condition, effect of reversive deformation, fuzzy phase transition.

При большой глубине, когда горное давление превышает предел прочности на одноосное сжатие горных пород, вокруг выработок возникает явление зональной дезинтеграции. Аномальность этого явления заключается в локализации трещин отрыва на дискретных зонах массива, чередующихся с зонами относительно не нарушенных пород (рис. 1).

1

УН

Рис. 1. Зональная дезинтеграция вокруг выработок

SlSlfl-'

Рис. 2. Образец породы под сжатием (а) и кривые деформрования образца

породы

Явление дискретных разрушений впереди лавы выработки впервые обнаружено G.R. Adams [9], и явление зональной дезинтеграции вокруг выработки на большой глубине Е.И. Шемякиным, М.В. Курленей, Ю.С. Кузнецовым, A.J. Jager и другими [2, 9, 15]. Это явление отличается от упруго-пластического решения данной проблемы, что и вызвало большой интерес у многих исследователей. Систематические in situ и экспериментальные исследования проведены Е.И. Шемякиным, М.В. Курленей и другими [5, 10-12, 15]. Эти исследователи и другие считают, что явление зональной дезинтеграции вызвано расклиниванием породы в зоне опорного давления вдоль максимального сжимающего напряжения, при малом боковом давлении [5, 11]. Эксперименты констатируют [5], что из-за существования слабой прослойки между блоками породы деформация не удовлетворяет условию совместности Сен-Венана, условие малости бокового давления всегда может выполнятся, и расклинивание осуществимо. В.Н. Одинцев, анализирует явление зональной дезинтеграции с точки зрения геологической структуры [6,

7], считает, что появление зон разрушенности связано с образованием новой устойчивой структуры. Она возникает в квазистатическом процессе деформирования из-за динамики разрушения. Эта новая структура сохраняет устойчивость при последующем нагружении определенного диапазона нагружения. Такой процесс повторяется, в результате чего возникает зональная дезинтеграция. Это явление аналогично пространственно-временной структуре в физико-химических процессах. Рассматривая соотношение между пространственно-временными параметрами, явление зональной дезинтеграции может быть интерпретировано аналогично соответствующим уравнениям квантовой механики. Эксперименты показывают, что процесс деформирования горных пород можно рассматривать как многостадийный размытый фазовый переход [1]. С точки зрения неравновесной термодинамики горные породы на большой глубине находятся в сильно неравновесном состоянии после экскавации выработки. По некоторым источникам, сильно неравновесное состояние приводит к инициированию процесса самоорганизации и образованию диссипативных структур [18], что необходимо учитывать при моделировании.

При исследовании зональной дезинтеграции В.В. Макаров открыл эффект смены знака приращения деформаций образцов горных пород при сильном сжатии [3, 4]. Он и другие исследователи при экспериментах на предварительно

поврежденных цилиндрических образцах на сжатие обнаружили, что, когда

*

осевая нагрузка а 2 превышает определенную величину а , с приращением нагрузки появляется смены знака приращения продольных и окружных

96

деформаций. Для образца, показанного в рис. 2 (а), типичная кривая деформирования образца породы показана на рис. 2(Ь) [3]. Измерение таких деформаций образцов свидетельствует о том, что существуют две аномалии: «отрицательная аномалия», и «положительная аномалия». Сущность положительной аномалии состоит в том, что одновременно с появлением отрицательного приращения деформаций в определенных участках образца в соседних участках появляется значительное положительное приращение деформаций, которое на 1.5-3 раза больше, чем среднее значение. Рис. 3 демонстрирует разнознаковое приращение деформаций по высоте образца при трехточечной схеме испытания.

Верхний торец

Ь (тт) /2 \ 1

70 /' •

35 -2 0 /I с 4 2 4

Нижний горец СХЮ3

Рис. 3. Разнознаковое приращение деформаций по высоте образца

Рис. 4. Периодическое распределение приращения объемных деформаций при трехточечной схеме испытания (кривая 1 соответствует объемных деформаций (кривая 2) по периметру диапазону изменения нагружения 136144 МРа, а кривая 2 образца (круг 1) - 144 - 150МРа) (абсцисса -приращение деформаций разделенное на приращение напряжения, т.е. податливость) .

При анализе объемных деформаций обнаружен периодический характер изменения приращения объемных деформаций по периметру образца (рис.4) [3]. Эксперниментальное и теоретическое исследования показывают, что процессы трещинообразования являются причиной смены знака приращения деформаций [3].

В физике теория фазового перехода второго порядка успешно объяснила сверхпроводимость, сверхтекучесть и явление периодической вихревой структуры в сверхпроводимости второго рода [14, 17, 19]. Если процесс деформирования горных пород является многостадийным размытым фазовым переходом, то его эволюция должна подчиняться законам фазового перехода. Поэтому в данной работе сделана попытка объяснить явление зональной

дезинтеграции с точки зрения фазового перехода.

Определение большой глубины с точки зрения классической теории пластичности

Обсудим напряженно-деформированное состояние горных пород вокруг подземных выработок.

Круглые подземные выработки без крепи подтверждаются внешним нагружением величиной <-уН, где у - плотность горной породы, H -глубина выработок. В упругом состоянии напряженно-деформированное состояние горных пород вокруг подземных выработок можно выразить следующими выражениями (здесь сжимаюшие напряжения предполагаются положительными):

Jr/V-J О + (о-ОО- qia VГ 2 )' J Z =СТОЭ

Т „2

s*-+J^t. sz -0

2^0 r

где f > <jv > < > sr > sv > sz - напряжения и деформации в полярной координате r and z ; < - гидростатическое давление в нетронутом горном массиве; — q внешнее давление, действующее на стенке подземных выработок;

- модуль сдвига горных пород; a - радиус круглой выработки; r- текущая координата.

В этом случае имеются три взаимноперпедикулярные максимальные сдвигающие напряжения:

T =

ap-°r

; T =-

a -a

z r

; t2 =

a (p a z

В качестве простого условия пластичности мы принимаем критерий Треска, который изображен на рис. 5, где предел текучести обозначен Т0, и предел текучести при одноосном сжатии а0. С ростом глубины главные сдвигающие напряжения растут постепенно. Когда Т = Т0 , пластические деформации впервые появляются на стенке выработок. В этом случае мы имеем

ap-ar _ a0

= ~Г~ = a-= T0 =т

т.е.

ao = T0 = ao/2

Рис. 5. Зависимость напряжения от деформации для идеальной

пластичности

2

2

2

Согласно теории пластичности С.А. Христиановича и Е.И. Шемякина, пластичность считается неполной, если толко T достигает предела текучести [8].

С дальнейшим ростом глубины условие пластичности выполняется и для т

T = = ^ = т0 = a .

1 2 2 0 2

В этом случае кроме плоскости r ,p, линии скольжения также появляются в

плоскости г, г. Такая пластичность называется полной, согласно концепции С.А. Христиановича и Е.И. Шемякина [8]. Горные породы разбиваются на блоки, что делает более ярким проявление геодинамических событий (например, горных ударов). При полной пластичности из (6) получим

^оо= 2Г0 =ст0 .

Поэтому с точки зрения появления динамических событий большую глубину можно считать глубиной, при которой гидростатическое давление достигает предела прочности при одноосном сжатии.

В качестве управляющего параметра возьмем максимальную деформацию сдвига

Г = ^г.

А в качестве модели горных пород принимаем модель идеальной пластичности, согласно которой при достижении сдвигающим напряжением Т предела текучести Т горные породы переходят в пластическое состояние. Но эти породы разрушатся только при достижении сдвигом Г предельного значения Г0, т.е. при

г = г .

Для хрупких пород Г=Г, где Г - предел упругой деформации.

Когда горные породы вокруг выработок переходят в пластическое состояние , напряжения в горных породах выражаются следующим образом

^ = 2То1п г 2То

<р

а \ а

Радиус зоны пластической деформации определяется формулой

2 ст»-То

- I = е а

т.е.радиус пластической зоны быстро растет с ростом глубины. В случае, когда = 2Т0 =а0 , из уравнения (11) получим радиус пластической зоны: с = 1.65а.

На основании упругого изменения объема и гипотезы Кармана перемещения

в зоне необратимой деформации определяются следующим образом [13]:

и =

Т011 -Ч , О

О

г 1п— + -

г г

где Э - постоянная, которая определяется из условий на границе г = с с упругой зоной.

На границе г = с имеем:

Г = р(р-рг =Ге = Т„/ О.

Сдвиг и компоненты деформаций определяются формулами

г=Г

-1 + 2у + 2(1 -1

р = Г

а г Г е

(1 - 21)1 1п а -11 + (1 -у)-

р = Г

р г Г е

(1 - 2у)1п а -(1 -у)-г I г

Когда максимальная деформация сдвига достигает своего критического значения г = г0 , стенка выработок начинает разрушаться, в этом случае радиус пластической зоны может определяться из

2

Г0

V Ге

Л

+1 - 2у

Из (16) когда г0/ге = 10иу = 0.3, получим - «2.73а. В этом случае имеем

СТоэ * 3Т0 — 1.5^0

2

с

г

2

2

2

1

Деформация на стенке выработок определяется как:

г =г

а е

зд 10

-1 + 2у + 2(1 -у)е То

Особенности деформирования горных пород при внешнем нагружении

Эксперименты на образцах горных пород показывают, что напряжение в породе, достигнув максимума, начинает падать при продолжающемся росте деформаций из-за возникновения локальных поверхностей скольжения [1]. В этих поверхностях в процессе дальнейшего деформирования постепенно подготавливаются микроразрушения. Можно связывать пластическую деформацию с образованием плоскостей скольжения, которые близки к площадкам действия максимальных касательных напряжений. Поэтому кривые зависимости максимального касательного напряжения Т от максимальной сдвиговой деформации г содержат основную информацию о деформировании и разрушении горных пород.

На рис. 6 показаны паспортные зависимости максимального касательного напряжения Т от максимальной сдвиговой деформации г (кривая 1) и зависимости относительного изменения объема в от максимальной сдвиговой деформации г (кривая 2).

В точке те деформирование переходит от упругого участка ( г<ге,т < те) в участок упрочнения, последний участок переходит в точке (г = г0, т = т ) в участок разупрочнения. После достижения критического сдвига Т наступает

разрушение. Объемная деформация уменьшается до предела упругости, после чего постепенно увеличивается. в «0 в точке ( т0 , г0 ). При дальнейшем увеличении деформаций объемная деформация положительна.

На рис. 7 показана кинетика образования микротрещин [18], гдеж/жтах -отношение текущего количества трещин n к максимальному числу трещин ытах в момент разрушения.

о в

к 2

"7

г / 1 1Тс

го г г г с

Рис. 6. Кривые т

0.2 0.4 0.6 0.8 1-0 Рис. 7. Зависимость от

# = г/гс

Соотношение можно рассматривать как параметр повреждения с . В области упругих деформаций с растет незаметно. При переходе пластическое состояние с начинает быстро расти. Помимо чисто сдвигового процесса, существенный вклад в деформацию материала дает разрыхление, сопровождающееся потерей связности и наличием поверхностей отрыва. Кроме того, на этапе упрочнения следует ожидать появление локализованных областей, в которых частицы материала измельчаются. В образовавшейся структуре зарождаются ростки новой, более сложной структуры, переход к которой совершается при смене типа устойчивости. Наступает новое состояние, при

0

в

котором в материале при дальнейшем росте деформаций происходит падение напряжения. Параметр повреждения в запредельной области быстро возрастает и при появлении магистральных трещин стремится к единице.

Процесс деформирования горных пород можно рассматривать как многостадийный размытый фазовой переход. На рис. 7 видны характерные для размытого фазового перехода особенности. Переход материала в запредельную область сопровождается сменой типа устойчивости, обусловленной сильными структурными флуктуациями. Поэтому наличие на графике точки перегиба, соответствующей пределу прочности материала, вполне предсказуемо.

Такой вывод позволяет нам исследовать зональную дезинтеграцию с помощью теории фазового перехода.

Параметр порядка запредельного деформирования

и управляющее уравнение

Для описания изменения запредельного состояния горных пород

целесообразно ввести параметр порядка:

. <19>

1 С 1 е

В горных породах существует структурная иерархия. На каждом уровне имеются структурные ослабленные поверхности - здесь главным образом концентрируются деформации и разрушение горных пород. На каждом уровне раскрытие таких структурных поверхностей пропорционально характерному размеру данного уровня. Естественно предположить, что чем меньше раскрытие структурных поверхностей, т.е. чем меньше характерный размер

данного уровня, тем выше здесь механические характеристики горных пород. При внешнем нагружении разрушение сначала происходит на уровне с большим размером. С увеличением нагрузки деформирование и разрушение постепенно охватывают уровни с меньшим размером.

В физике самой распространенной симметрией является однородность и изотропность пространства. В механике деформируемых твердых тел среда с изотропией обладает свойством симметрии, механическое поведение которой описывается двумя параметрами упругости. При разрушении среды изотропия исчезает, симметрия нарушается. При увеличении внешней нагрузки разрушение происходит в последовательно уменьшающихся структурных уровнях. При первом разрушении формируются блоки самого большого размера. В этом случае глобальная симметрия геосреды нарушается. Но вновь сформировавшиеся блоки можно рассматривать как однородные и изотропные, в этом случае симметрия локализируется в блоках. С дальнейшим ростом нагрузки размеры последовательно сформировавшихся блоков, следовательно, и размеры локальной симметрии постепенно уменьшаются. Поэтому разрушение горных пород можно рассматривать как процесс последовательного нарушения симметрии среды и процесс последовательной локализации симметрии.

Нарушение симметрии сопровождается возникновением необратимой деформации в горной породе, которое сопровождается возникновением диссипации и образованием диссипативной структуры [18]. Деформируемое

тело можно рассматривать как статически консервативную систему, которая

106

обладает общим потенциалом н.

Предполагая минимум н при равновесии, можно получить уравнение равновесия в виде

УН = о.

Можно принять термодинамический потенциал н в виде Н = Уо + Уу-1 Уу2 +1УУ +•••

Согласно особенности размытого фазового перехода мы имеем У > о, у > о.

Функция (21) должна удовлетворить условию равновесия при у =0, т.е.

У - Уу У У +--Уо = 0.

ёщ

Из (24) получим у =0; не теряя общности, можно полагать у = 0

При изучении изменения системы при нагружении от 0 до критического значения можно оставить члены щ до четвертого порядка:

H = -1 УЩ +1 УЩ.

Когда горная порода находится в предельном состоянии, согласно особенностям критических явлений, корреляционная длина увеличивается, т.е. между частицами породы появилось дальнодействующее взаимодействие. Поэтому необходимо ввести градиентный член в (6) по аналогии с разложением Гинзбурга - Ландау [1]

н = -1 ущщ2 +1у4щ4 +1 c

где C- коэффициент, и с >0.

Поэтому общий потенциал системы имеет вид

н = -1 г2г2 +1+1С

А полная свободная энергия есть

Ф = 1 +1 +1С (У у)2

О

¿о.

При устойчивом состоянии система выбирает значение параметров, которое делает потенциал минимальным, поэтому функция должна принять минимум. Варьируя ф относительно у, получим

дФ = | [- У2уду + У3у2ду + СУ у-У(ду)}ю

О

= |ду[- Уу+У3у2 ]ао+|сууу

о о

Из-за

У(Уу-ду) = У2у-ду + Уу-У(ду) мы имеем

Уу-У

Применяя для первого члена в правой стороне (29) метод интегрирования по

частям и теорему Гаусса - Остроградского, получим

| СУу-У(ду)ёО =

\ е

= Г ду (СУ уУШБ - Г СУ у - 8усЮ

О

где первый член в правой стороне (30) является поверхностным интегрированием,, dS - элемент поверхности, а ¿о - объемный элемент. Подставляя (30) в (28), получим

дФ = Г ду [- У у у + У у у3 ¿О + Г СУ уУ(ду)аЮ =

о о

= Г ду [- у2у + у4у3 - сУ 2у]ю + Гсду[

у\

о ¿о

Так как (31) справедливо для любой величины у и для любой точки в теле, мы имеем

У у-У у3 + СУ У = о

с граничным условием

Если в уравнении (32) только сохранить линейный член, то получим следующее уравнение:

у 2у+Бу = о

с граничным условием (33), где Б = У2/С.

Моделирование реверсивного деформирования образцов горных пород

Согласно условиям испытания [3], торцовые поверхности тщательно полированы на шлифовальных кругах с добавлением абразивных порошков. Если знаками г,в,х обозначены радиальная, окружная и аксиальная координаты, то граничные условия можно написать следующим образом [13]:

аг\ п и=-&>Тгг\ к и= о гч. = о ?сгг| = о ' тгп\ = о ' тА = о .

2\г=о,п =о,Л 2а> 2=о,п г\у=а га> г=а гг\г=а

Напряженное состояние, удовлетворяющее уравнениям равновесия и граничным условиям, равно

^ = -<У,&г =°в= Тгв = ТК =Т2в= о .

Когда нагрузка о2 меньше критического значения ст* , поведение материала подчиняется закону Гука, и ненулевые компоненты деформации определяются следующим образом:

а а а

— ,£в = У — £г = У — ,

Е в Е г Е

где Е — модуль Юнга; у - коэффициент Пуасона.

Два максимальных сдвиговых напряжения равны

т = = а = т23 = .

2 2 23 2

А соответствующие две максимальные сдвиговые деформации определяются следующим образом:

Т = ев-е2 =

Е

Из уравнений (4) , ( 5 ) видно, что, два максимальных сдвиговых напряжений равны между собой, как и соответствующие деформации. С ростом нагрузки они достигают критических значений одновременно. Поэтому использование единого параметра порядка у для описания запредельного поведения горных образцов является правомерным.

В цилиндрической координате при осесимметричном напряженно-деформационном состоянии (34) и (33) выглядят

д2у у 1 д 2у д 2у о Л ^ллл

—Ьт + — + ^—Ьт + —Ьт + БАу = 0

дг2 гдг г2 дв2 дт2

Для решения уравнения (40) используем метод разделения переменных [19]. Выражаем параметр порядка в форме произведения трех функций: у(г,в,г)=я(г)-©(в)- 2(т.), после подстановки чего в уравнение (40) получим 3 независимые уравнения

^ + ж = 0

¿г2

2 / x

dR + ™ + | D2-Я-41Л = 0,

V r2 у

dr 2 rdr

+ л = 0.

d 20 de2

где i, л — коэффициенты, которые не зависят от координат r,e,z. Решение уравнения (44) имеет вид

©

Используем условие периодичности ©(о)=©(2^) и ©'(о)=©'(2^), получим

л = m2, m = 0,1,2,-••

Поэтому уравнение (45) приобретает вид

то

0

m=0

где Am, Bm - коэффициенты.

Решение уравнения (42) имеет следующий вид:

Z

где коэффициенты P, Q определяются из граничных условий Z'(о) = 0, Z '(h)=0

P=0, A = ,и = 0,1,2,-

Поэтому уравнение (21) становится

то

Z

и=0

где Qn — коэффициенты.

Для уравнения (43 ) введем новые переменные: р = л/Б2-Я • г = кг, где к = л/Б2 -Я )' ,у(р) = Я(г)' получим

д.. ( 2\

д ym_ + dym |

m

y.m = 0 ,

др2 рдр

где ут является соответствующей параметру т величиной у(р). Уравнение (51) является уравнением Бесселя т-го порядка, решение которого имеет следующий вид:

Ут

где Jm и Ыт являются функциями Бесселя и Неймана т-го порядка. Из условия ограниченности величины у

Общее решение уравнения (51) у(р) является суперпозицией всех решений уравнений (25) по т.

то то

у

т=0 т=0

Согласно граничным условиям на поверхности образца:

у'(р)р=р = 0 при р = 4в2 -я • а = ра

получим

-} юр

=0

y'iPa )= Z °m\jm(Pa )= Z Cm1

m=0 m=0

Jm-l

P

Таким образом, можно определить ра, а затем D.

Поэтому общее решение уравнения (40) можно представить в следующем виде

^ п=да ftZTT

у(г,0,z) = Z[cml(Ат sinm0 + 5m cosm6)jт (kr)]• ZQn COS—-

m=0 n=0 h

да n да vf7TT

= £[

m=0 n=0 h

где A'm, B'm - коэффициенты.

Для описания дилатансии горных пород можно принять пропорциональное соотношение между объемной деформацией и сдвигом следующим образом [20]:

eV = Лу ,

где Л - коэффициент пропорциональности.

Для фиксированной вертикальной координаты z на поверхности образца r = a, если принять вертикальную ось на рис. 3 как начало отсчета угловой координаты, можно принять 5 гармоник для аппроксимации распределения объемной деформации по периметру образцу.

ev = Л у = -

Некоторые результаты расчета по формуле (57) показаны в таблице 1 и изображены на рис. 6 кривой 3. Как видно из рис. 3, аппроксимация достаточно хорошая.

Таблица 1

Некоторые результаты расчета

в(°) 0 30 45 90 120 150 180 210 240

270 315 330 345 360

eV (10-3) -1415 -394 257 399 500 208 -320 8 316

187 107 -500 -1150 -1415

Рис. 8. Сравнение результатов расчета по (58) (кривая 3) с

экспериментальными данными

Рассмотрим аппроксимацию распределения деформации по высоте образца. Мы приняли относительный сдвиг как параметр порядка, а экспериментальные данные являются аксиальной деформацией е2 и окружной деформацией ев . Для аппроксимации экспериментальных данных нужно делать некоторое упрощение. Из рис. 2 (Ь) видно, что в определенной точке на поверхности образца приращения вертикальной деформации Ав2 и окружной деформации А^в имеют противоположные знаки, но примерно пропорциональны по величине. Поэтому можно принять следующее соотношение:

А^в = — лАе2,

где 7 — коэффициент пропорциональности. Из-за того, что АГ = Авв - Аа2, мы имеем

Аб2 =кАГ ,

где к — коэффициент пропорциональности.

Для фиксированной угловой координаты в , поскольку датчики расположены

114

в 3 точках по высоте, то согласно решению (56) для аппроксимации кривых 1 и 2 на рис. 3 достаточно сохранить 3 члена:

(-) (тт тт -п тт 2-п

—= \ И + H cos — + И cos-

Аа I 01 h 2 h

•10

-3

После расчета получим следующее выражение для аппроксимации кривой 1

А^Л-) ( -п 2-п \ т. -у ' = 1 2.63 + 0.173cos--0.57cos-10-3.

Actz ^ h h

и изображено кривой 3 на рис. 9

Для кривой 2 мы имеем следующее выражение для аппроксимации

= \ 1.733 +1.617 • cos П- 2.267 • cos ^ • 10-3

- 2.267 • cos — 1-10-3 h h

и изображено как кривая 4 на рис. 9.

верхний торец

h (mm) 4 \ уЪ

70 2 W

35 -2 0 2 4

нижний торец

CxlCf;

Рис. 9 Сравнение результатов по уравнению

Из результатов физического моделирования эффекта смены знака приращения деформаций видно, что подход теории размытого фазового перехода достаточно хорошо описывает закритическую деформационную аномальность пород.

Моделирование зонального разрушения массива вокруг подземных

Для круглой выработки при осесимметричном

напряженно-деформационном состоянии (34) выглядит

^ д2ш дш^

I ¡Г

V

+ Б 2ш = 0.

дг гдг

У

Уравнение (64) является однородным уравнением Бесселя, решение которого имеет вид [19]

ш = СЛ

где 70 и N - функция Бесселя нулевого порядка и функция Неймана нулевого порядка соответственно (цилиндрические функции); а С и С2 -подгоночные параметры.

В полярной системе координат граничное условие (33) имеет вид

-ш/-Гг=в = 0.

Далее используем следующие итерационные соотношения

а [70

, '/-,1^/ II— 1 17 11 У г,

из (66) получим

-Ш-Аг=а = -Б[С171 (Ба)+ С2N1 (Ба)] = 0 , То есть С2 = -С 7 (Ба )/N (Ба).

Сдвиговая деформация не имеет положительного или отрицательного знаков, поэтому мы можем принять абсолютную величину параметра ш

w=

J 0 D )-S N0 D )"

Из у=(г-г. V(rc -ге) вытекает

г = ^

Выбор параметра с осуществляется таким образом, что на краю выработки решение (68) совпадает с упругопластическим решением. Параметр D зависит от радиуса выработки и подбирается таким образом, чтобы период функции у был равен расстоянию от края выработки до первой зоны дезинтеграции.

Принимая Da = 3.2 (i.e. D = 3.2/а ), v = 0.3 , гс/ге = 10 , = 3Т0 = 1.5 сг0 , по условию

г(а) =

Q

J 0 D )- N0 D )

получим С «1.862 о .

Далее из (68) определяем распределение деформации в массиве вокруг выработок. Результаты расчета показаны на рис. 5, где решение (68) обозначено кривой 1, а решение по идеальной пластической модели изображено кривой 2. Если в качестве первого приближения считаем, что область появления зональной дезинтеграции ограничена областью необратимых деформаций, то из рис. 5 видно, что кроме зоны разрушения вблизи стенки выработки имеется еще одна зона дезинтеграции. Этот результат совпадает с экспериментальными данными [10].

Такой результат может показаться абсурдным с первого взгляда, потому что слаборазрушенная зона имеет деформацию, равную величине деформации при

117

пределе упругости. Но эксперименты действительно показывают, что участки массива, расположенные между зонами дезинтеграции, нарушены значительно слабо, и трещины в них расположены хаотически. При разборке моделей создается впечатление, что трещины на этих участках возникли позднее, в результате перемещения их в сторону выработки [4].

Необходимо подчеркнуть, что образование узора деформаций является популярным в процессе деформирования материалов. Экспериментальные наблюдения, представленные в технической литературе, убедительно показывают, что в макроскопически однородных областях материала при поверхностном нагружении развиваются неоднородные и часто высоколокализованные поля деформаций, что противоречит прогнозированию локальных теорий [21]. Поэтому образование зональной дезинтеграции является не случайным, это реализация общего закона деформирования материалов. Совместное использование классических упругопластических теорий и общих законов физики позволяет понять явление зональной дезинтеграции более глубоко.

Л г/г.

2

г

Рис. 11. Пластическая деформация сдвига в породе вокруг выработки

Таким образом, в работе показана единая природа явлений зональной дезинтеграции и реверсивного деформирования образцов горных пород.

Финансовая поддержка: национальный научный фонд Китая (^РС)^о.50825403), Пекинский научный фонд (N0.^200810016007).

СПИСОК ЛИТЕРАТУРЫ

1. Адигамов Н.С., Рудаев Я.И. Уравнение состояния, учитывающие разупрочнение материала // ФТПРПИ. Физико-технические проблемы разработки полезных ископаемых, 1996. № 4. С. 24-32.

2. Глушихин Ф.П., Кузнецов Г.Н., Шклярский М.Ф. и др. Моделирование в геомеханике. М.: Недра, 1991. 240 с.

3. Гузев М.А., Макаров В.В. Деформирование и разрушение сильно сжатых горных пород вокруг выработок. Владивосток: Дальнаука, 2007. 232 с.

4. Гузев М.А., Макаров В.В., Ушаков А.А. Моделирование поведения сжатых

горных образцов в предразрущающей области // ФТПРПИ. 2006. № 5. С. 3-13.

5. Курленя М.В., Опарин В.Н. и др. О расклинивающем эффекте зон опорного давления // ФТПРПИ. 1995. № 3. С. 7-17.

6. Одинцев В.Н. О механизме зональной дезинтеграции массива горных пород вблизи глубоких выработок // ФТПРПИ. 1994. № 4. С. 10-19.

7. Одинцев В.Н. Отрывное разрушение массива горных пород. М.: ИПКОН РАН1996. 166 с.

8. Христианович С.А., Шемякин Е.И. К теории идеальной пластичности // Механика Твердого Тела. 1967. № 4. С. 86- 97.

9. Шемякин Е.И., Фисенко Г.Л. Курленя М.В., Опарин В.Н. и др. Зональной дезинтеграции горной породы вокруг подземных выработок. Ч. 1. Данные натурных наблюдений // ФТПРПИ. 1986. № 3. С. 3-15.

10. Шемякин Е.И., Фисенко Г.Л., Курленя М.В., Опарин В.Н. и др. Зональной дезинтеграции горной породы вокруг подземноых выработок. Ч. II: Разрушение горных пород на моделях из эквивалентных материалов. ФТПРПИ // 1986. № 4. С. 3-13.

11. Шемякин Е.И, Фисенко Г.Л., Курленя М.В., Опарин В.Н. и др. Зональная

119

дезинтеграция горных пород вокруг подземных выработок. Ч. III. Теоретическое представление // ФТПРПИ. 1987. № 1. С. 3-8.

12. Шемякин Е.И, Фисенко Г.Л., Курленя М.В., Опарин В.Н. Зональная дезинтеграция горных пород вокруг подземных выработок. Ч. IV. Практические приложения // ФТПРПИ.

1989. № 4. С. 3-9.

13. Шемякин Е.И. Две задачи механики горных пород, связанные с освоением глубоких месторождений руды и угля. ФТПРПИ. 1975. № 6. С. 29-45.

14. Abrikosov A.A. Nobel lecture: Type II superconductor and the vortex lattice // Rev. of Modern Physics. 2004. V. 76. P. 975-979.

15. Adams G.R., Jager A .J. Petroscopic observations of rock fracturing ahead of stop faces in deep-level gold mines // J. of The South Africa Institute of Mining and Metallurgy. 1980. V. 80. N 6. P. 204-209.

16. Asmar N.H. Partial differential equation with Fourier series and boundary value problems. Beijing: China Machine Press, 2005. 252 p.

17. Ginzburg V.L. Nobel lecture: On superconductivity and superfluidity as well as on the "physical minimum" at the beginning of the XXI century // Rev. of Modern Physics. 2004. V. 76. P. 981-998.

18. Haken H. Advanced synergetics. Berlin: Springer-Verlag, 1983. 384 p.

19. Landau L.D., Lifshits E.M. Statistical physics / Physical Institute. Beijing: Beijing World Publishing Corporation, 1999. 252 p.

20. Nikolaevsky V. N. Mechanics of porous and fractured media. Singapore: World Science,

1990. 284 p.

21. Valanis K.C. A gradient thermodynamic theory of self-organization // Acta Mechanica. 1998. N 127. P. 1-23