Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Механика деформирования и разрушения горных пород

..pdf
Скачиваний:
4
Добавлен:
12.11.2023
Размер:
8.26 Mб
Скачать

Перейдем теперь к нахождению коэффициентов 1п и 1р. Система для их определения при горизонтальных пластических элементах ((3 = 0), записанная в матричном виде, выглядит так

 

-Тп

 

ТР

0

 

 

1„

=

Л''

12 sin 2а+1

, (3.43)

 

Pu

 

 

- Р р

 

 

*2~

 

 

tu

 

 

0

~ta

 

 

0

 

 

®1.

 

_ о 2

 

 

 

 

 

 

 

 

где 5,= >

S2

<*п'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решения

системы

ь I ' 1' И I

 

 

 

членов

 

(3.43) со

столбцами свободных

 

 

 

 

 

 

 

 

 

1

 

 

0

 

 

 

 

 

 

 

 

 

0

И

 

1

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

Тогда

искомые

интенсивности

равны

 

 

 

 

 

 

I,-= I <1»((5l -

J2 )/2 sin 2<х +1 )+ 1 f

52-

(3.44)

Определитель системы

(3.43)

записывается

 

 

 

 

 

Л = /с(ТпЛ,-Л,Тр) + /п Л>Тс-

(3.45)

Решения 1Ф с помощью определителя Д можно записать так

 

 

 

 

1(1>=Р

1 Д-1-

I (1)= P

t Д-1-

 

 

 

 

 

1 J

1

р *с

ч

12

1 п

1с ^

ч

 

 

 

 

 

1«з1» = Рр/пД -1;

1<2>= т„ Г Ж 1;

(3.46)

 

 

 

12:' = (х„ tc—тс t„ ) Д- и,

I (32)= тр t„Д- 1

 

Обозначим через

 

коэффициенты интенсивности, вычислен­

ные по формуле (3.41) через значения коэффициентов К1’ (/=1, 2).

Вводя

наибольшее

касательное напряжение Ï = (S1—S2)/2,

критерий

(3.44) роста

для вертикальных трещин

запишется

в виде

 

 

 

 

*{‘41 + т sin 2а) + К j2)s2 = Kf/xn,

(3.47)

где К\ — вязкость разрушения.

Выразим из условия (3.47) значение т касательного на­

пряжения:

 

т= ц+ Ъ 2 + п —.

(3-48)

где ц= —1/sin2сх; А = —K\2)/Kjl) 1/sin2а;

т| = 1/(АГ*1>sin2ос).

Связь (3.48) позволяет наити наибольшее касательное на­ пряжение в зависимости от напряжения s2, вязкости разрушения К\, коэффициентов ц, А., т|, являющихся функциями длины трещины 21, и вида функции а п.

Теперь перейдем к записи в явном виде значения параметров X, г) в зависимости от длины трещины.

Для этого заметим, что зависимость (3.41) может быть

представлена

в виде

 

 

*■ = 1п/п + 1р./р,

(3.49)

где /„=/„(/);

/ р =/„(/).

 

Тогда

 

 

* !1) = I 'Л /п + I (21 ,./р = teА -1(./„ Рр+/р Ра);

 

к \ 2) = I (2>УП+ 1( )/р=Д - 1[У п Тр /с - У р (тп /с - /п тс )] ,

(3.50)

где

I

 

 

 

*■'; Ур=/>0

-I

С учетом связей (3.50) коэффициенты А, и т| легко могут быть представлены в виде функции длины трещины.

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

ê„. = jje 0(-v, j’)dfi/j{dQ,

(3.51)

п

а

 

где Q— геометрическая область, соответствующая геометричес­ кому периоду задачи.

В системе координат (д-^О ;^), связанной с пластическим элементом, тензор средних деформаций системы пластических элементов

"(Х+1) 4G

—sin 2 (а —(3)

cos 2 (а —Р)

cos 2 (а —Р)

sin 2 (а —Р)

где

х = 3—4v;

v— коэффициент

Пуассона; G— модуль сдвига.

для

В системе координат (.vOj’ ) тензор средних

деформаций

системы

пластических

деформаций выглядит

так:

 

 

е = л(х+1)

- s in 2а

cos2а

(3.52)

 

 

4G

cos 2а

sin 2а

 

 

В этом случае в системе координат (дОу) деформации

системы трещин разрыва

определяются зависимостью

 

 

P - K / l X + D

1

0

 

.

4G

0 1

которая в системе

координат (л'Оу ) записывается:

;P _ 7t/(X+1)

 

COS2 p

 

cos Psin P

£ н = -

4G

cos Psin P

(3.53)

 

sin2 P

Аналогично для

трещин

сдвига

имеем

 

е = я/(х+ О

0

1

 

 

 

8G

1

О

и

 

 

 

 

 

„с

я/(х+1)

—sin2P

cos2p

8 “

8G

 

cos2P

(3.54)

 

sin2p

Деформации, определяемые упругим полем внешних напряжений

1+V

Е е

Е

va! —(1—v )a2

О

О

(3.55)

vcr2 —(1—v)a!

где Е — модуль упругости.

среды

 

Тогда средние деформации неоднородной

 

 

eij= InE"j+ Ip£f2+Iceij + efj-

 

(3.56)

Значения главных перемещений еь е2 для вертикальных

трещин

после преобразований записываются

так:

 

 

Ei = [(v—1 ) T +

( V 0,5)52 + к (1 —v) In sin 2<x] тп/<7;

 

82 =

[ V T + ( V - 0 , 5 ) 5 2 -

J I (1 — V ) Insin2a+n/(l - v ) I p] t n/G.

(3.57)

Учитывая реологические соотношёния (3.57), уравнения связи

Для деформаций преобразуем к виду

 

 

 

zl ={\il + 'klsl +T\l K\lxn)xJG\

 

 

 

Е2 = (Р2 + >-2*2 + П 2 ^ / О т„/<7,

 

(3.58)

где

 

 

 

 

|Д2= ц (v —1); p2 = vp; A,i=A.(v—1)+ я(1 —v)(M(11,sin2a+I (2))sin2a; X2 = À,v+v —0,5 —л (1 —v)(A.I<11,sin 2 a+ I(12))sin 2a+

+ ni (1 - v)I " >sin 2a+ 1 f );

r| j = r| [v—1+ 7 t l(11)(l —vsin2 2a)];

T|2 = Л [v—1+ л(1 —v) —л (1 —v)sin22 a(/I<21)—I <11,sin2a)].

Уравнение (3.56) в совокупности с уравнением (3.31) позволя­ ет получить замкнутые реологические соотношения дефор­ мирования неоднородных горных пород, связывающие про­ чностные, деформационные характеристики и структурные Параметры.

Рис. 3.26. Зависимости напряжения Астрах —а 3 от деформации ех д ля гранита

(а) и невыбросоопасного песчаника (б)

Параметры моделирования можно разделить на 3 группы.

Во-первых, это

силовые параметры — упругие константы

£ и V и предел упругости х. При моделировании они берутся

непосредственно

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

Во-вторых, это геометрические параметры: отношение длин сторон упруго-пластического элемента, углы Р и а (либо выражающиеся через последний угол отношением критической длины трещины к стороне элемента (см. рис. 3.25). Эти параметры определяют форму, но не размеры элемента. Учитывая нормированность соотношений (3.43), при модели­ ровании можно выбирать любые, например, единичные раз­ меры упруго-пластических элементов. Параметры подбираются так, чтобы получающиеся реологические кривые максимально соответствовали реальным.

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

K\ = K\JS/ R J R „

144

где Ru — размер упруго-пластического элемента, принятый в модели; Л,— реальный размер зерна.

Однако вследствие высокой степени неоднородности горных пород непосредственное экспериментальное определение К']3 за­ труднено, поэтому К\ подбирается исходя из соотношения (3.48) при 52= 0 таким образом, чтобы при одноосном сжатии пределы прочности реальной породы и модели совпадали.

В виде примера приведем построение модели неоднородной горной породы применительно к экспериментальным данным для гранита и невыбросоопасного песчаника. Деформационные характеристики материала среды для этих пород приняты соответственно следующими: £=0,55 • 105 МПа и v= 0,17; £=0,3 • 105 МПа и v= 0,12. Отношение длины трещины к длине структурного элемента в обоих случаях составляло 1:44, а Р= 0.

В качестве структурного элемента среды принят прямо­ угольник с соотношением сторон 2,35:1 для гранита и 2,9:1 для песчаника.

Сопоставление расчетных и экспериментальных значений Аа1= ст1 —ст2 в зависимости от et для различных а 2 приведено на рис. 3.26. Анализ приведенных данных показывает, что математическая модель неоднородной породы удовлетворитель­ но описывает деформационное разупрочнение горных пород.

3.7. СИСТЕМА УРАВНЕНИЙ ПЛАСТИЧЕСКОГО ДЕФОРМИРОВАНИЯ ПРИ РАЗУПРОЧНЕНИИ В УСЛОВИЯХ ПЛОСКОЙ ДЕФОРМАЦИИ

При решении различных прикладных задач горной геомеха­ ники значительное распространение получила задача о плоской деформации.

Вработе [32] приведены различные зависимости для аппроксимации экспериментальных данных между напряже­ ниями и деформациями при упрочнении горных пород, разупрочнении и в области остаточной прочности. Эти зависимости представлены набором нелинейных и кусочно­ линейных функций.

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

r = f ( a ) - M ( y- у .), (у^у.),

(3.59)

где Y = v/(ei - e y)2 + 4yi,,; а = 0,5(а,-|-а,.); у = е1- е 3;

М — мо­

дуль спада; т = / (а) — условие предельного состояния на грани­ це раздела упругой и пластически разупрочняющейся области.

Рассмотрим систему уравнений пластического деформирова­ ния при разупрочнении в условиях плоской деформации.

Компоненты напряжений ах, ау, тху в этой области должны удовлетворять условию предельного состояния с разупрочне­ нием (3.59) и уравнениям равновесия.

Результаты экспериментальных данных за деформированием пород показывают, что при разупрочнении их в запредельной области имеет место (см. рис. 3.5, 3.8) пластическая дилатансия. При исследовании системы ограничимся случаем линейной дилатансии, уравнение которой для плоской деформации записывается следующим образом: ._____________

ди/дх + dv/dy = А ^/(du/dx dv/ôy)2 + (du/dy+ dv/dx)2 + Ai, (3.60)

где и, v— компоненты перемещений; A, A i— постоянные. Для замыкания системы необходимо к уравнениям присо­

единить еще одно соотношение. Для этого используем рас­ пространенное в теории пластичности условие соосности деви-

аторов напряжений

и деформаций.

 

 

 

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

 

стх ау/(2хху)= (du/dx — dv/dy)/(du/dy + dv/dx).

(3.61)

Таким

образом,

для пяти неизвестных

ах,

ау,

хху, и,

v имеем

пять

уравнений.

виде

 

Компоненты

напряжений будем искать в

 

 

 

 

 

 

 

(3.62)

где т=/(сг) —М{у у,); а = 1/2 (стх + аД; р — угол

между

поло­

жительным направлением оси Ох и направлением максималь­ ного главного напряжения а 4.

Компоненты напряжений (3.62) удовлетворяют условию

(3.59). Тогда

система уравнений для нахождения неизвестных

ст, р, и, v состоит из уравнения (3.60), уравнений

 

da/dx+dx/dxcos 2Р+ 2т sin 2$d$/dx +dx/dy sin 2Р+

 

 

+ 2т cos 2§d$/dy = 0;

 

dx/dx sin 2P + 2т cos 2$d$/dx -I- dajdy — dx/dy cos 2p —

 

 

—2Tsin2Pc!P/^v = 0,

 

где

 

 

dxidx =f

(CT)dçs/dx — Mdy/dx; dx/dy = /' (a) da/dy — Mdy/dy;

 

у = J (du/dx — dv/dy)2 +(du/dy + dv/dx)2

(3.63)

и условия коаксиальности (3.61) тензоров напряжений и дефор­ маций, которое может быть записано так:

(du/dx dv/dy) sin 2P = (du/dy + dv/dx) cos 2p.

(3.64)

Рассмотрим частный случай приведенной выше системы уравнений. Условие пластичности разупрочняющейся среды примем в виде

T = 4tf-4A f(y-y,),

(3.65)

где К — пластическая постоянная.

Заметим, что условие (3.65) представляет собой условие пластичности Мизеса, обобщенное на разупрочнящуюся среду.

На первом этапе исследования системы дилатансию пород не будем учитывать, считая, что имеет место условие несжима­

емости.

виде:

Компоненты напряжений будем искать в

М = а + 4 [К —М(у —y.)]cos2P;

(3.66)

Gy)

 

хху= 4 [ К - М ( у - у.)] sin 2р.

 

Соотношения (3.66) тождественно удовлетворяют условию (3.65). Используя соотношения (3.66) уравнения равновесия, условие соосности (3.64) и условие несжимаемости, запишем замкнутую систему дифференциальных уравнений относительно неизвестных ст, и, v:

да

 

 

 

дх

 

 

 

 

, дI2и д2и

 

 

~ 4M |, à ? + ï ÿ l=0;

да .

ал \ f ^ П {du

<9iA~|

ô [ \ f d v ôu\

ду +ЦК+М,.)

 

 

 

I d2u

d2v 1

du/dx + dv/dv = 0. (3.67)

 

—4М( — + —у ) = 0;

 

дх~

dv

 

Система уравнений (3.67) является системой второго поряд­ ка. Можно показать, используя условие неопределенности вторых производных, что система (3.67) является эллиптической.

Заметим, что условие предельного состояния (3.65) для разупрочняющейся среды имеет такой же вид, как условие предельного состояния для вязко-пластической среды:

 

т = 4£+4ру,

(3.68)

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

вязкости; у— скорость

наибольшего

сдвига.

зрения для вязко-пластической среды

С .физической точки

с увеличением скорости у предельное напряжение (3.68) воз­ растает, а для пластически разупрочняющейся (3.65) с ростом

у напряжение т

убывает. Аналогом коэффициента вязкости

P для пластически разупрочняющейся среды является модуль

спада.

на приведенной выше аналогии, для задач

Основываясь

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

вид, аналогичный

(3.67), и приведена в работе Г А. Гениева

и М. И. Эстрина

(1972 г.).

3.8. ДЕФОРМАЦИОННОЕ РАЗУПРОЧНЕНИЕ ЛЕНТОЧНЫХ ЦЕЛИКОВ

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

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

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

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

В настоящее время методы расчета целиков в достаточной степени разработаны [40] при упругом и упруго-пластическом режимах их деформирования. Попытки учета запредельных свойств при расчете целиков были предприняты в работах [5, 19]. При этом вводился ряд допущений, в частности, предполагалось, что в запредельное состояние сразу полностью переходит любое вертикальное сечение или весь целик. Диа­ граммы напряжение—деформация для образцов и целиков приняты идентичными.

Исследовалось упруго-пластическое деформирование лен­ точных целиков с учетом их разупрочнения. Анализ напряжен- но-деформировг юго состояния производился с помощью вариационных методов и численных экспериментов на ЭВМ.

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

На ниспадающих участках зависимости интенсивность ка­ сательных напряжений Т — интенсивность деформаций сдвига Г —вопрос'о существовании экстремумов функционалов остает­

ся открытым.

задач

теории

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

потенциальная

Для

двумерных

энергия

системы записывается в виде

 

 

 

 

 

 

 

 

п(«.

tH

J

^ d5+ W s(/ ’)d S -/l,

 

 

(3-69)

 

 

 

 

 

(S) СТЛ

 

,5)

 

 

 

 

 

где 0— объемная деформация;

К — модуль объемного сжатия;

А — работа

внешних

сил;

g (Г) — работа

деформации

формы;

и, V— значения перемещений;

S — область, занимаемая

телом.

Необходимым условием того, что минимум функционала

(3.69)

достигается в точке (и0, i>°) является

 

 

 

 

 

 

 

5п(м°

г°) = 0;

82п (м°, г°)>0.

 

 

(3.70)

Условие. положительности

функционала может

быть

запи­

сано

в

виде

 

 

 

 

 

 

 

 

 

 

 

 

 

£ц(51)/С ?-М р(52)/С1>0,

 

 

(3.71)

где

E=m'mdT/dr;

dT/dr>0\

М = -тах б Г /б Г ;

dT/dr<0;

С) =тахГ;

С2 = т т Г ;

S,,

S2— соответственно

области

до­

предельного

и

 

запредельного

состояний;

p(Si),

p(S2) — площади этих

областей.

 

 

 

 

 

 

Заметим, что критерий (3.71) может быть обобщен на объемные задачи. При этом структура критерия сохраняется, следует только заменить значения площадей S, и S2 на соответствующие значения объемов Vx и V2.

Из выражения (3.71) следует, что минимум функционала существует до тех пор пока зона разупрочнения мала. Даль­

нейшее

увеличение нагрузки приводит к росту этой зоны, что

может

привести к нарушению необходимого условия. При

решении конкретных прикладных задач

это

часто связано

с неограниченным ростом перемещений,

что

свидетельствует

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

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

При исследовании напряженно-деформационного состояния целиков примем кусочно-линейную аппроксимацию зависимо­

сти

интенсивность касательных напряжений

Т — интенсивность

деформаций

сдвига Г,

состоящую из трех участков:

 

 

 

 

(GT

 

при Г < Г п

(3.72)

 

 

 

Т(Г) = <М( Г —Г„)+ Т„ при Г „ ^ Г ^ Г 0,

 

 

 

 

( Т0

 

при Г > Г 0

 

где

G,

М — соответственно

модули

сдвига

и спада; Тп, Гп,

Т0,

Г0— соответственно

значения интенсивностей касательных

напряжений

и

деформаций

сдвига

на пределах

прочности

и остаточной

прочности.

 

связь

между

объемной

 

На

первом

этапе

исследований

деформацией 0 и средним давлением

будем принимать в виде

 

 

0 = ЗЯа.

 

(3.73)

Работа

деформации

формы g (Г)

в энергии системы (3.69)

с учетом

связи (3.72)

записывается

 

 

g ( r ) = f r ( r ) d r =

 

 

 

О

 

 

 

при Г <Г„;

GГ 2/2

 

 

 

СГ 2/2 + [Т„+ М + {Г -Г„ )/2\{Г - Гп)

при Г „ ^ Г ^ Г 0- (3.74)

СГ 2/2 + (Т„+ Т0)(Гп — Г0)/2+Т0(Г—Г0) при Г0< Г

Деформирование однородных целиков

в режиме заданных

нагрузок. При проектировании целиков значительное распространение получила известная схема расчета Турнера — Шевякова [40].

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

ала

целика.

 

в

Методы расчета целиков, использующие известные нагрузки,

том

числе по

схеме Турнера — Шевякова, могут быть

отнесены

к методам

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

 

Наряду с этой группой методов получили развитие методы,

учитывающие взаимодействие целиков с вмещающими поро­ дами. Среди них в первую очередь следует отметить работы А. К. Черникова [40]. в которых рассматривается взаимодейст-