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

книги из ГПНТБ / Скуба, В. Н. Исследование устойчивости горных выработок в условиях многолетней мерзлоты

.pdf
Скачиваний:
8
Добавлен:
20.10.2023
Размер:
7.98 Mб
Скачать

Т а б л и ц а 13

Показатели

Размер трещин,

мм

0,5

2

4

 

Максимальная величина области со­

 

 

 

прикосновения блоков, мм:

 

3,5

 

расчетная ..................................

6,6

1,3

замеренная в модели . . .

5—7

3 ^ 6

1—0

Разрушающая нагрузка, кг:

 

 

0,5

расчетн ая ..................................

12

5,8

замеренная в модели . . .

12—11

4—6

1,0

в табл. 13. Там же указаны значения этих переменных величин, рассчитанные аналитически. Величины области соприкоснове­ ния блоков определены из выражения (IV.5), а разрушающая нагрузка, приложенная в центре, рассчитана по формуле, предложенной для таких случаев А. А. Борисовым (1948):

Р =

2<ТС Ж ^ б (h ~ ~ у)

(IV.6)

 

L

 

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

Зная величину соприкосновения блоков Сь , можно опреде­ лить стрелу прогиба арки (рис. 19), которая равна расстоянию между равнодействующими горизонтального распора в замке и пяте блока. Величина стрелы прогиба зависит от величины опускания замка, размера области соприкосновения блоков и характера распределения напряжений в шарнирах. Прини­ маем эпюру распределения напряжений в шарнирах для хруп­ ких пород близкой к треугольной, для пластических пород — к прямоугольной.

Тогда

 

h

= h t -----з-6 - у

,

(IV.7)

 

и

= К - С 6 - у .

(IV.8)

Горизонтальный

распор

составляет:

 

 

кино

Г * =

т

^~2с1

Г*

(IV.9)

инедо

 

8 v

6 — ~з- ~~~У)

 

ш

=

QL2

 

(IV.10)

ГуП

8 (h6 ~~

У)

Условия прочности арки:

 

 

tfx сж

qL2

 

(IV.11)

2Cj

 

4С

3

 

 

 

 

сж

__________

(IV.12)

(h6 ~

~ у)

 

 

Из уравнений (IV.11)

и (IV.12)

можно

определить макси­

мальное значение пролета, при котором возможно равновесие системы, т. е. предельный пролет незакрепленной кровли. При этом необходимо рассчитать ряд значений Съ для несколь­ ких пролетов и, подставив их в выражения (IV. 11) или (IV. 12), убедиться, какие из них удовлетворяют данному условию рав­ новесия.

По настоящей методике рассчитаны предельные пролеты незакрепленной кровли для нескольких моделей, отработанных в лаборатории горного давления Ленинградского горного ин­ ститута. Результаты расчетов хорошо согласуются с данными опытов. Так, предельный пролет для модели из песчано-пара­ финовой смеси (ах сж = 12 кг/см3, h = 1) равен 37 см при опуска­ нии среднего шарнира на 1,5—2,0 мм (ширина трещины соот­ ветствовала сжатию блока). При отработке этой модели пре­ дельный пролет незакрепленной кровли оказался равным

40—42 см.

По формуле (IV.5) рассчитаны величины области соприкос­ новения блоков, а по формуле (IV. 11) — величины предельных пролетов незакрепленной кровли для двух практических слу­ чаев. В первом случае кровля состояла из алевролита с толщи­ ной слоев 0,2 м, а во втором — из песчаника с толщиной слоев 1,0—1,5 м. Породы были разбиты трещинами, перпендикуляр­ ными напластованию, густотой от 0,2 до 3 на 1 пог. м и шириной до 0,2—0,5 см. Предельный пролет кровли из алевролита имеет величину примерно 1—2 м, в шахтных условиях такая кровля обрушается в выработках шириной немногим более 1 м. Пре­ дельный пролет кровли из песчаника составляет 7—8 м, в натуре такая кровля в подготовительных выработках при

обнажении

не обрушается и ее устойчивый пролет находится

в

пределах

5—6 м. Хорошая сходимость расчетных величин

с

данными

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

рекомендовать рассматриваемый

метод расчета предельных

обнажений

оттаивающих

пород

с открытыми трещинами для

решения практических

задач.

 

53

§ 2. Приближенные методы определения нагрузок на крепь при оттаивании интенсивно-трещиноватых пород

вокруг выработок

Впервые для определения нагрузок на крепь при оттаива­ нии пород А. Ф. Зильберборд (I960) предложил исходить из полного веса столба пород над крепью в пределах ореола оттаи­ вания:

 

Р —г LR TyСр .

(IV.13)

Справедливо

отмечая недостаток

формулы (IV. 13),

Ю. Д. Дядькин (1968) предлагает при расчете нагрузок на крепь учитывать пригрузку от мерзлых пород. Однако при рассмотре­ нии схемы нагружения крепи автор исходит из частного слу­ чая, когда между обрушившимися на крепь породами (талыми, а затем и мерзлыми) есть свободная полость (щель), образую­ щаяся вследствие оттока воды в выработку, деформации крепи и ее вдавливания в почву. В этом случае учитывается[пригрузка только от обрушенных мерзлых пород в пределах свода об­ рушения, а влияние массива пород не принимается во внимание.

При таком допущении интересно следующее предположе­ ние: при оттаивании пород контур выработки будет охваты­ ваться все расширяющейся зоной пластических деформаций, так как упругопластические деформации беспрерывно распрост­ раняются от границы оттаивания к мерзлому массиву. Поэтому при оценке нагрузок на крепь необходимо исходить из веса пород в пределах ореола оттаивания и разрушенных пород в пределах пластической зоны за контурами оттаивания, па­ раметры которой могут быть определены по К. В. Руппенейту (1957). Расчетная схема для нашего случая представлена на рис. 20. Тогда

^= V i ‘5i (t) + v2S2(t).

(VI. 14)

Для выработки с радиусом R 0 границей пластических де­ формаций при радиусе оттаивания R T будет эллипс с полуо­ сями ап и Ьп:

ап = С(1 + Р); 6а = с(1 - Р ) ,

(IV. 15)

где

(IV.16)

\'у гНт

(IV . 17)

54

Рис. 20. Расчетная схема для опреде

Рис. 21. Упрощенная схема

ления нагрузок на крепь при оттаива­

нагружения крепи при от­

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

таивании

горных

пород.

вокруг выработки.

 

 

 

Учитывая,

что

 

 

 

=

- г й ) ’ v = 4 - ( 1 - ^ 1 . *.2 =

г = - р

(IV. 18)

получаем

 

 

 

 

 

Г Ь1

( З о |- 1) + 1

 

(IV.19)

 

Р = ^ r(Y i+ Y 2'

 

В общем случае при оттаивании пород нагружение крепи идет по схеме, упрощенный вариант которой показан на рис. 21.

Нагрузка на крепь,

отнесенная к 1 м выработки,

составляет

Р =

2ау [кНт+ (1 — к) R т],

(IV.20)

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

при О ^Ге минус 10°; к =

0,08 |Ге|—°»8е °*31^Ге1

Подставляя

значение к и R T (гл. III, § 2),

определим

 

Р =2 аг[0,08|Г еГ 0’8е~°'3|Те'Яг +

(1 -

0,08) 1Ге|—0 8е~0-3|'Ге' ~ Х

L

 

 

 

Г е|

х ( 1 _

е“ Щ

г)]-

(IV.21)

Зависимость (IV.21) характеризует максимально возможные нагрузки, развивающиеся при оттаивании пород в горных вы­ работках.

55

РуГП/М

 

 

 

 

 

 

 

 

 

80-

 

 

5

Рис.

22.

Зависимость

 

 

^

нагрузок

на крепь

от

40-

 

 

1

при различных

ореолах

 

 

 

оттаивания пород

(1— 6

 

 

 

 

соответственноUUlbCltlBBHHU 1,1,; 2,2;U

О

100

200 300

400 Нп м

3,1;

4,5;

6,3;

7,5

м).

 

 

 

 

 

 

На рис. 22 приведены результаты расчетов нагрузок на

крепь,

выполненных на

ЭВМ М-220,

в зависимости от hд при

различных значениях R w и постоянных а= 1,5

м, у =2,3

т/м3

и Т е= —3° С.

расчетов

показывают,

что

в большинстве слу­

Результаты

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

§ 3. Аналитический метод определения перемещений контура выработки

и границы пластических деформаций при оттаивании пород

Устойчивость выработок при оттаивании пород зависит от глубины оттаивания пород и положения границы пластической зоны вокруг выработок. Допущение, что пластическая зона вокруг выработок охватывает ореол оттаивания пород, позво­ лило нам (Дубина, Красовицкий, Скуба, Первенцев, 1974) по­ дойти к аналитическому решению задачи определения переме­ щений контура выработки и границы пластического деформиро­ вания. Исходя из методов механики сплошной среды и решая уравнение плоской деформации, найдем напряженно-деформи­ рованное состояние пород вокруг выработки при условии пла­ стичности в виде прямолинейной огибающей кругов Мора.

Для безразмерной координаты r= R /R 0 любой точки массива (рис. 23) уравнения равновесия при плоской деформации в об­ ласти /<><[оо имеют следующий вид:

(ГУ.22)

дтгв, , 1 doet 2тг01 _

дг + г ’ дд1 +

г “ 0

56

 

 

 

III *в

а

 

- ( ■ - 1 - 1 - 1 J -

 

м з д Г А .

 

 

 

----

 

 

 

d

/

 

/

/ С

 

 

 

1

I

\?А

0

 

L

— (

\

f e w

Я )

)

—»-(

Область ///

-

 

 

__1

 

 

 

^ т т т т н ' " T i n t 1

 

 

V

811

 

 

 

 

 

 

 

Рис. 23. Схема к рас­ чету перемещений на контуре выработки и движения границы пла­ стических деформаций.

Решение системы (IV.22) при граничных условиях на конту­ ре выработки L q

ol = Р ■

= о при г = -'*0 = 1

(IV.23)

возможно при введении в нее функции напряжений Ф таким образом, чтобы выраженные через нее компоненты напряжения после подстановки в (IV.22) удовлетворяли системе (IV.22) тождественно:

Функция напряжений Ф в талой зоне (область 7, ограни­ ченная контурами L 0 и LT) должна на границе LTудовлетворять также условиям непрерывности напряжений

I _ II

I

_ 11ш

I

II

JlT

=

(IV.25)

Op — Ог ,

О01

— И0!,

*гвх=

Тг01 При —

 

 

 

 

 

 

-“о

 

 

где oj, g^, тг01, компоненты напряжения в талой зоне (об­

ласть

/); о*1,

о01г, т£01 — компоненты

напряжения в

области

II, ограниченной контурами LT и L.

 

 

Так

как в

областях I

и I I породы находятся в состоянии

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

то величина

Ф должна удовлетворять условию

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

 

 

 

 

(а01 — (тг) + т>01 = sin 2p (а01 + or +

2/cctg р)2.

(IV .26)

57

Предполагая, что Ф не зависит от 0Х, имеем Tret=0. Тогда условие (IV.25) после подстановки (IV.24) дает уравнение для определения Ф:

ф „ _ 1 -f- sin р ^_Ф_'

2к ctg р

1 — sin р г

1 — sin р

Отсюда получаем общий интеграл для Ф:

Ф =

Рг2

Аг«+®

+ £,

2 (1 - б)

1 + 6

где

 

 

 

 

1 + sin р

ctg р

 

1 — sin р ’

Р = 1 — sin р *

(IV.27)

(IV.28)

Подставляя Ф в уравнения (IV.24) и учитывая граничные условия (IV.23), находим'постоянные А г= А 2 и В = В Х для об­

ласти I. Постоянную В = В г можно приравнять нулю,

ибо она

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

не

влияет. Тогда + 1=P+/cctgp1,

где

рх

и кх

— параметры

Мора для

области /.

Если

обозначить сг[ =

2 = овр

получим

следующие

выражения

для компонент на­

пряжения в области I:

 

 

 

 

 

 

 

 

 

о[ +

A: ctg Pl = ( - 1){ 1

1 +

sinpl

(Р + кг ctg Pl) r«s

(IV.29)

 

 

 

 

(1 — г)г —sm pL

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_

2 sin P!

 

 

 

 

 

 

 

 

 

 

 

1

1 —sin pT

 

 

 

 

 

В

области

I I

функция

напряжений

Ф

выражается

той

же формулой (IV.28), а постоянная

А = А 2 находится

из

гра­

ничных

условий

(IV. 23),

 

которые являются

условиями

непрерывности

напряжений

при

переходе

через

границу

зоны Ьт.

Тогда

 

 

 

 

 

 

 

 

 

 

 

А 2 =

[k2ctg р2 — A^ctg pi +

+ кх ctg рх) г?1] г^~ Ьг\

(IV.30)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6,

1 +

sin Р2

 

 

 

 

 

 

 

 

 

 

 

1 — sin р2 '

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

58

Таким образом, компоненты напряжений в области I I опре­ деляются из уравнения

ОV + /с2 Ctg р2 = ( — !) *

А 2

i — 1 +

sin р2 гД,

(IV.31)

(1 — i)1

— sin р2

 

 

 

где

 

 

 

 

_

2 sin р2

 

 

^2

1 — sin р2*

 

 

Исходя из метода К. В. Руппенейта (1954), определяем пе­ ремещения на границе выработки и границу пластического де­ формирования. При этом для определения напряжений в упру­ гой зоне (в области III) необходимо решить уравнение упруго­ пластической границы L.

В области I I I компоненты напряжения должны удовлетво­

рять условию совместности деформаций,

которое при xret = О

запишется в виде

 

 

 

 

- | r [rsi?(ff1.- O 0 I)] =

О,

(IY.32)

откуда

 

 

 

 

R

=

С3

 

 

 

 

 

г2 (°г ~

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

ди

R .

 

и

 

 

R .

 

 

 

дг

4Gc

а0*)’

~ Г

4GC(а01 ~

^

 

откуда выражение перемещения

 

 

 

 

 

 

 

 

“ = - 4 § Т -

 

 

 

 

(1V.33)

В области I I I

функция напряжения Ф

должна удовлетво­

рять уравнению совместности

 

 

 

 

 

 

 

___ д_

j _

e а2_\

 

 

 

дФ

 

j _

а2Ф \

• г ' дг

"*■ г

* 302 J ^

дг*

1

г

дг

+

г

* 502 J ” 0

 

 

 

 

 

 

 

 

 

(IV .34)

59

и условиям на бесконечности

С7°°)

 

 

т^ 4 = Х'усрН г sinr0i. (IV.35)

М (Яз — A, cosrGj yCptfr;

CTeJ

 

 

 

Условия будут удовлетворены, если задать компоненты

напряжения

в форме

 

Or =

У с р # г

[я,3 (1 — а г ~ 2) V

(1 — 2 Ь г~ 2 + c r - 4) cos 2 0 j;

cre4 =

уср# г [Яз (1

ar~2) -f X' (1 -f- cr~k) cos 2 0 х] ;

^re, =

Уср#Д' (1 +

br~ 2 — сг~4) sin 20х.

 

 

 

 

(IV.36)

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

L. Уравнение

определяется

в виде

 

 

 

rL = r0 + kr1(Q1)t

(IV.37)

где г и

не известны.

 

 

Предполагается X' — малый параметр, поэтому справедливы

следующие равенства:

 

 

 

 

Гь 2 =

Го 2 [1

2Xr0 \ (6 j)];

 

 

гГ ‘ =

г5-4 [1-4Х 'г5"Ч (в1)];

(IV.38)

 

rL =

ro [1 +

аХго 1г1 (©!)].

 

Подставляя (IV.38) в (IV.29) и (IV.36) и уравнивая при ма­ лом параметре А, получаем следующие уравнения для опре­ деления а, в, с, т*0, 71(0 4 :

1 4~ Ьго2 — сго4 = 0;

4 2a 2ro2~Vi (0 i) 4 - усрЯ г [( 1 2Ъг^24- c r ^ ) cos 2 0 х —

— 2X3a r^3r 1(0Х)] = 0;

у42г? 2 к ctg р2 — уСрЯДз (1 а г^2) = 0 ;

Л (1 4 - sin р2) г? 2 (1 — sin р2) к2ctg р2 — уср# Д з (1 +

4 ~ar^2) (1 — sin р2) = 0 ;

А 2(1 -f sin р2) a 2r£2“ Vi (04 — [(1 4- сг^к) cos20i —

2X3a r^3r1 (0 j) (1 — sin yp2) ycp# r = 0 .

(IV .39)

60

Отсюда

JX t

1 — sin p2

 

r0

 

-------- J~

(^зТср^г + k2cl& Рг)>

 

 

 

" sin p2

 

 

 

V cp

(Л з7 ср # Г + k2 Ctg p2) ;

(IV.40)

 

 

ЭУсрЯрГо cos 20!

 

(0i) =

 

 

 

4 sin

HT + k2ctg f

 

b =

2r20;

c =

3rJ.

 

Таким

образом, уравнение границы L

 

r L

=

r 0

' 4 j______ 3^Ycptfrcos 20!

(IV.41)

 

4 sin p2 (Яд7 срЯг + k2ctg p2)

 

 

 

 

 

при г=гь, т. e. на упругопластической границе R = 1, тогда

: -

а 34,г£‘*+2),

(IV.42)

или приближенно

 

 

с3 = -

a 2A / ^ +i).

(IV .43)

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

(a2+2)

(IV.44)

“ “4G7

 

Для проверки полученных теоретических результатов на ЭВМ «НаириС» был произведен машинный расчет перемеще­ ний по формуле (IV.44) и

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

чи:

уср = 3 т/м3;

 

U =0,9;

 

X'

=0,1;

р ^ Б 0;

р2=30°;

 

кх = 117,2

кГ/см2;

/с2=

 

=792 кГ /см2. Для # г=500м

 

построены Гь при Гт—4 (кри

рис ^ Изменение положения пла-

вая

и .при гх

1,о (кри-

стических деформаций при различ-

вая В 2). Для Н т~

700 м пост-

ных начальных условиях.

61

Соседние файлы в папке книги из ГПНТБ