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

книги / Проектирование и расчёт крепи капитальных выработок

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

величины средней нагрузки Рср при различных механических моделях массива приводятся, например, в работах [10, 12].

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

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

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

Однако полученное решение, как показывает приведенный в [12] модельный пример, в котором в качестве измеренной нормаль­ ной нагрузки задавались значения, найденные из решения кон­ тактной задачи теории упругости [41 ], хотя и дает возможность с хорошей степенью точности восстановить эпюру касательных на­ пряжений, тем не менее не решает вопроса расчета крепи по изме­ ренным нагрузкам. При детальном рассмотрении данного примера обнаруживается, что даже сравнительно небольшие отклонения касательных напряжений от воспроизводимой эпюры, полученной при решении контактной задачи, приводят к существенным иска­ жениям величин расчетных изгибающих моментов в крепи. Кроме того, указанное решение неустойчиво, так как весьма небольшие изменения в исходной эпюре измеренных нормальных напряже­ ний, как показывают расчеты, выполненные для ряда примеров, приводят в большинстве случаев к резкому изменению изгибающих моментов в крепи. Необходимо отметить также, что такой путь (равно как если бы имелись измерения нормальных и касательных контактных напряжений) вообще может дать возможность расчета

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

А,* = Ы2/Ыг и угла наклона главных осей а, а также корректирую­ щего множителя а* по результатам натурных измерений давления на крепь горных выработок и подземных сооружений, и в после­ дующем использовании найденных характеристик в качестве ис­ ходных данных для расчета крепи по методике, описанной в п. 6.1,

6. 2.

Указанные характеристики, найденные на основе натурных измерений, в отличие от нагрузок не зависят от параметров кон­ струкции и могут использоваться для расчета и проектирования другой, отличной от объекта измерений крепи в сходных горно­ геологических условиях при аналогичном способе возведения.

Такой путь расширяет возможности использования натурных

измерений при проектировании

крепи и дает наиболее близкие

к реальным результаты расчета

ее напряженного состояния, хотя

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

6.4.1. Постановка и решение обратной задачи об определении расчетных характеристик массива горных пород

Предлагаемая методика основана на аналитическом решении [41 ] плоской контактной задачи теории упругости для среды, ос­ лабленной отверстием произвольной формы (с одной осью симмет­ рии), подкрепленным кольцом из другого материала, имеющей начальное напряженное состояние с наклонными главными осями (расчетная схема дана на рис. 6.1). Задача ставится как обратная

и состоит в определении таких величин Nxа*, А,* = N а, при которых расчетная эпюра нормальных контактных напряжений ар

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

В нашем случае контактные напряжения в любой k-той точке измерения в силу линейности задачи могут быть представлены как суммы

<>„.=«• м ч '» > + w

+ r i w a

(6.89)

где о ?, о у , txy — вертикальные, горизонтальные и касательные начальные напряжения в ненарушенном массиве пород; ор1, о

ОрЦ — контактные напряжения в 6-той точке, полученные из ре­ шения задачи теории упругости, расчетная схема которой приве­ дена на рис. 6.1, соответственно в частных случаях:

1— N1cc* =

l,

Nacc* = 0, а = 0;

2—Nta* = О,

Naа * = 1 ,

а

= 0 ; 3— N1a* = l,

N2а * = — 1,

 

 

а = я/4.

 

На основании решения, алгоритм которого приведен в п. 6.2,

имеем

 

 

 

 

 

 

 

 

 

+

ар (1а) +

 

[ор (16) cos 2а +

 

 

 

+ ap(II) sin 2a],

 

 

(6.90)

 

откуда

 

 

 

 

 

 

 

 

o $ r= y - K ftO ^ +

Op^16)];

afk = Y

[ap* (Ia) —opk (16)];

 

 

 

 

 

 

 

 

(6.91)

 

 

 

 

<X$ =

<Tp*(II).

 

 

 

Подставляя

(6.91) в формулу (6.89) и приравнивая в п >

3 точ­

ках

измерений

расчетные

значения opk к

измеренным

сгр*

(k =

= 1,

. , /г),

придем

к

переопределенной системе п

линейных

уравнений относительно трех неизвестных величин начальных на­ пряжений ах = а(х0)а*у оу = а|,0)а*, хху = т{хуа*у имеющей вид

°х [vpk(la) + oQk(16)] + оу [Op* (la) —apk (16)] + 2тxyo9k (II) = 2ар* (6.92)

(* = 1 ,2 ,

п).

Как известно, такую переопределенную систему можно свести к симметричной относительно диагонали системе трех уравнений вида

01,Л

+ ai,2°V + 01,з

= *ъ

 

02,1°^

+ 02,2^у

+

02,3^ х у

=

*2i

(6 • 93)

0 3 , Л

~+~ а З Л ° у

+

0 3 , 3 т х // =

* 3 ,

 

где коэффициенты матрицы и свободные члены выражаются фор­ мулами:

п

п

\о1к(1а)-с1к(1б)1

ai,i = Z К П 1а) + арЛ1б)12>

01,2= 02,1 = У

к=1

/с-1

п

01,з= 0з,1 = 2 У Орд^(II) [а9к (I#) -f-а9к(16)],

/с=1

 

fl2,2 =

Ё К * (la) — <V (I6)]2,

 

(6.94)

 

 

K=1

 

 

 

 

#2,3 — #3,2 — 2

^ Q K (II)

[СГрдс (I#)

S ACC1 6 )]»

#3,3

4 ^ Gp/c(II),

 

AC= 1

 

 

 

 

 

II

 

 

lb

 

 

 

bi = 2 2 CTpK[crp/c (la) + стрк(16)],

62 = 2 E

0p*

[огрл (la)—apx(16)],

/C=l

 

 

/C=l

 

 

 

&з = 4 £ огрл(Трк(П).

K=1

Решив систему (6.93) и определив неизвестные a*, ap и тху, можно найти величину главного напряжения

Ni<x — [crx + Оу + д/(сг*— Oyf -f 4т*р ],

(6.95)

отношение главных напряжений

N2

ах + Q y - /\J(ox - ° y f + ^ly

к*

(6.96)

Nl

2

^ + ар + д /(0 Л- а р)2 + 4т:'ху

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

игоризонтали

а =

где

р

 

 

при

гху

о ,— о„>0

1 OJ

.со

при

 

 

°х° у < 0

л/2 -f р

при

*ху< 0,

СГд;—СГу<0

ао

 

при

т ^ < 0 ,

Ох—Оу>0

0

 

 

при

II о

1

V о

л/4

 

 

при

%ху

0,

ох— Оу= 0

я/2

 

 

при

Т'ху=

 

ох— Оу<0

Зя/2

 

при

Т*!/<0,

1

«? II о

 

Р = —

arctg ■ 2 1чХ 1

 

 

 

 

2

 

>? 1

Q

 

 

 

 

 

«с

 

 

(6.97)

(6.98)

Входящие в формулы (6.94) значения ap/c (I а), арк (I б), арА (II) вычисляются по алгоритму расчета крепи в массиве, подвер­ женном действию тектонических сил, приведенному в п. 6.2 при

значениях углов

= 1,

, п). Дальнейший расчет крепи

производится по

указанному

выше алгоритму, приведенному в

п. 6.2, с использованием величин Nxa*, А,*, а в качестве исходных данных. Отметим, что полученное выше решение обратной задачи является единственным и точным, т. е. если в качестве измеренных

величин a*pft принять полученные в результате решения прямой

задачи при некоторых значениях Л^а*, Я*, а нормальные контакт­ ные напряжения, то после решения обратной задачи эти величины получатся теми же, которые задавались.

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

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

Исходными данными для

расчета

являются: заданные форма

и размеры поперечного сечения крепи;

Eltvx — модуль деформа­

ции и коэффициент Пуассона

материала крепи; Е 0, v0 — модуль

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

тактных напряжений

<Хр* (k = 1,

2,

, п).

Порядок обработки результатов измерений следующий.

1. Строится эпюра

измеренных

величин <Тр, значения которых

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

2. Строится конформное отображение внешности единичной окружности на внешность внутреннего контура поперечного се­

чения крепи и определяются коэффициенты а0, аъ

аъ ряда

отображающей функции, имеющей

вид:

 

г = < о ( У = £

а £ \~ у

(6.99)

v=0

 

3. Определяются координаты х , у точек внутреннего контура поперечного сечения крепи, соответствующих точкам единичной окружности при изменении угла 0 от 0 до 180° через 15° по фор­ мулам:

5

х = Y J avcos(l — v) 0, v=0

5

 

у = Y avsin( 1—v)0.

(6.100)

v=0

 

Эти точки — их 13 на половине периметра — наносятся на чер­ теж (заодно проверяется и точность конформного отображения), а на второй половине периметра отмечаются точки, симметричные с нанесенными относительно вертикальной оси. Всего получаются 24 точки, из которых проводятся нормали к внутреннему контуру

сечения крепи и с чертежа снимаются значения напряжений сгр,

которые

принимаются как

исходные

измеренные величины орк

(к = 1,

24) в точках,

лежащих

на пересечении нормалей с

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

 

9«= n ( ]

7 l)

(K=I>2’

24)-

(6Л01)

4.

Производится

расчет

по алгоритму,

приведенному

в

п. 6.2,

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

0Я,

определяемых формулой (6.101), вычисляются

значения арл (I а),

«V (I б), (V (И).

Составляется с использованием формул (6.94) система трех уравнений (6.93), после решения которой по формулам (6.95) —

(6.98) находятся величины Л^а*, X*, а. Указанные процедуры об­ работки данных измерений запрограммированы для ЕС ЭВМ на алгоритмическом языке ФОРТРАН (программа FOK6, автор про­ граммы А. Н. Козлов). Обработка результатов измерений вместе с последующим расчетом крепи при полученных исходных данных по алгоритму, изложенному в п. 6.2, составляет экспериментально­ аналитический метод расчета крепи произвольного поперечного сечения.

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

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

* Г

Ni Ч ~

Nj

(ci— 1)Н---- ——1 -(% —Зс3— 1) cos 2 (0 —a)J,

стр = а *

2

 

 

 

 

 

 

 

(6. 102)

 

где

 

 

 

 

 

 

Cl =•

Я.

 

Ся =

 

<*?(*!- О

 

1 — 2d (l - R \ )

sR\) (dR4 +

l) + 3tdR\ (R\ -

1)

 

(/ +

 

 

 

 

 

 

(6.103)

 

 

 

d R \ + l

 

 

 

 

 

CL\—------------ Co.

 

 

 

 

- ~1 , d =

* ~ P

, t=

1+XoP

,

£I (1+VQ)

 

 

4 ( 1 - v , )

 

4(1 - v , )

 

£o(l+Vx)

 

 

 

 

 

 

(6.104)

rx — внутренний,

r2 — наружный

радиус

поперечного сечения

крепи.

 

 

 

 

 

 

Тогда, согласно (6.62) и (6.66), имеем

 

Op (1а) = сх — 1,

сгр (16) = (ах—Зс3— 1) cos 20,

(6.105)

Op(II) = (ах—3с8— 1) sin 20.

 

Соответственно, учитывая формулы (6.94) и обозначив

 

А = сх— 1, В = ах—Зс3— 1,

(6.106)

 

■J

 

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

(6.93) относительно неизвестных оХУ оуу

тху в виде

 

01.1 = Z

И

+

Я cos 20Л)2,

а1Л= а%1 =

£

{А2 + В2cos2 20J ,

К=\

 

п

 

 

 

 

К=\

п

 

 

 

 

 

 

 

 

 

01.8= 08.1 =

2 Z

iA + В cos 20Д.) В sin 20к,

агл =

Z А — В cos 20х)2,

 

/С— 1

 

 

 

 

 

К — 1

 

 

 

 

п

 

 

 

 

п

(6.107)

 

 

 

 

 

 

 

 

02,з = 03.2 =

2 Z

(А — В cos 20J В sin 20Л,

a3l3 = 4 Z

В2sin2 20к,

п

 

ж=1

 

 

п

 

/с=1

 

 

 

 

 

 

 

 

bi = 2 Z

а1к (А + в cos 20J ,

fe2 =

2 X

tfp* {А— В cos20к),

к=1

 

 

 

п

 

А=|

 

 

 

 

 

 

кВ sin 2QK.

 

 

 

 

 

 

 

&з = 4 £

 

 

 

 

 

 

 

к—l

 

 

 

 

 

После решения

системы величины

N xa*, X*, а

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

по формулам

(6.95) — (6.98).

 

 

 

 

 

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

0р= Nia* ^-L±

(Cl_

1) -f l ~- k* (ax—3c8— l)cos2(0—a )j ,

ae внеш = Nxa*

[(4d— 1) ct — 1]----- X* x

X [3(4d.— 1)c3 + (1—4/) a3 — 1] cos2 (0—a)},

aeвнутр= Nxa* ((1 + X*) (2cxd -1)—2(1 -Л *) x

X (3 c 3d R t — t a ^ r 2) COS 2 (0— a )),

M = q9 внутр - Об внеш_ b (Гг_ Г1)2>

N

^6 внутр ~f~g6

внеш

(6.108)

2

Ь{Г2— Гг),

 

 

 

6=1 М.

6.4.4. Пример расчета

Ниже в качестве иллюстрации приводятся результаты расчета, произ­ веденного на основе измерений нормальных нагрузок на крепь околоствольного двора шахты «Красная Звезда» гор. 770 м.

Исходные данные:

Е0 = 27,7 ГПа, v0 = 0,33; Е г = 37,2 ГПа, v x = 0,15, толщина железобе­ тонной крепи Д = 0,2 м.

Рис. 6.41. Эпюры контактных напряжений и усилий в крепи:

а — измеренные контактные нормальные напряжения (сплошная линия), расчетные на­ пряжения (пунктирная линия); б — расчетные продольные силы; в — расчетные изгибающие моменты

Эпюра измеренных нормальных нагрузок о*, заимствованная из работы

[21], показана на рис. 6.41, а сплошной линией. Необходимые для расчета значения измеренных нагрузок приведены ниже.

 

градус

а*. МПа

0, градус

а*, МПа

0, градус

ст*, МПа

0, градус

<jp. МПа

 

0

 

—0,16

90

—0,10

180

—0,68

270

—0,30

 

15

 

—0,62

105

—0,60

195

—0,74

285

—0,23

 

30

 

—0,60

120

-0 ,3 0

210

—0,64

300

—0.30

 

45

 

—0,50

135

—0,10

225

—0,29

315

—0,39

 

60

 

—0,21

150

—0,11

240

—0,41

330

—0,35

 

75

 

—0,30

165

—0,33

255

—0,11

345

—0.30

=

В

результате

расчета по программе FOK6 получены значения: ох =

1,6

МПа,

оу =

1,37 МПа, тху =

— 0,26 МПа, N xа* =

1,78 МПа, А* =

=

0,67,

а

=

— 33°.

 

 

 

 

 

 

Расчетные эпюры контактных напряжений ор, изгибающих моментов М

и продольных

N в

крепи

приведены соответственно на рис. 6.41, а (пунк­

тирная линия),

6,41, в и 6,41, б. Отметим, что полученный угол наклона

главных осей а

= — 33° показывает, что большее

главное напряжение в не­

нарушенном массиве в данном случае направлено

почти строго по нормали

к напластованию.

 

Расчет набрызгбетонной крепи

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

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

До настоящего времени методы расчета набрызгбетонной крепи с учетом неровностей поверхности и полного контакта с породами отсутствовали. Не останавливаясь подробно на имеющихся пред­ ложениях и рекомендациях по расчету набрызгбетонной крепи, не учитывающих указанных особенностей ее работы, например [34 j, где крепь рассматривается как прямоугольная пластинка, работаю­ щая на изгиб под заданной равномерной нагрузкой [35], где ре­ комендуется определение толщины набрызгбетонной крепи стволов на основе формулы Ламе для толстостенной трубы при равномер­ ной внешней нагрузке [16], где крепь и массив с учетом зоны, упрочненной проникновением набрызгбетона в трещины или ан­ керами, рассматриваются как многослойное круговое кольцо. На формуле для толщины крепи, приведенной в работе [1], полученной также из решения осесимметричной задачи в условиях задан­ ных деформаций, или на эмпирических рекомендациях области применения этого типа крепи в зависимости от коэффициента кре­ пости пород по М. М. Протодьяконову, отметим работы [20], в ко­ торых оценивается влияние неровностей поверхности на несущую способность набрызгбетонной крепи, однако крепь рассматривается как круговая цилиндрическая оболочка (со срединной поверхностью, имеющей неровности, аппроксимируемые гипотрохоидой) под дейст­ вием равномерного внешнего давления, т. е. вне массива пород.

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

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

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

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

7.1. Постановка задачи

Рассматривается плоская контактная задача теории упругости для кольца произвольной формы (с одной осью симметрии), под­ крепляющего отверстие в линейно-деформируемой однородной изотропной среде из другого материала, имеющей начальное поле

Соседние файлы в папке книги