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

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

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

Рис. 7.1. Очертание набрызгбетонной крепи в двух сечениях полевого штрека гор. 710 м

Рис. 7.2. Очертание

набрызгбетонной

крепи

в двух сечениях

квершлага

гор. 710 м

 

 

 

 

 

напряжений, обусловленное

собственным

весом

 

ai0) =

— уН а,

o f' =

ХуНа, rfj = 0,

(7.1)

где Н — глубина заложения выработки; А, — коэффициент боко­ вого давления в ненарушенном массиве пород; у — вес единицы объема пород; а* — корректирующий множитель, вводимый для приближенного учета влияния отставания крепления от проведе­ ния выработки, неупругого деформирования пород и пр., опреде­ ляемый одним из способов, описанных в п. 6.4.

Некруговое в общем случае кольцо имеет неровности, аппрок­ симируемые гипотрохоидальной кривой, параметры которой — число неровностей п г и средняя амплитуда 6 отклонений от про­ ектного контура — могут задаваться на основе статистической об­ работки результатов натурных измерений или согласно [25], при-

Рис.

7.3.

Очертание

набрызгбетонной

крепи

в двух

сече­

 

ниях

штрека гор. 800

м

 

 

 

 

 

 

чем

 

где

R — средний

радиус

выработки. В

качестве

примеров

на рис.

7.1, 7.2

и

7.3

приведены

формы

реальных

контуров

поперечных сечений

выработок,

закрепленных

набрызг­

бетонной

крепью,

предоставленные Криворожским

филиалом

ВНИИОМШСа, которым выполнены соответствующие измерения (для каждой выработки представлено по два сечения, обозначенные на рисунках а и б).

Число выступов и впадин (соответственно) /гвыст и я Впад, сред­ ние величины отклонений от проектного контура б и толщины набрызгбетона во впадинах А занесены в табл. 7.1.

Из табл. 7.1 видно, что в реальных выработках условие n xb< .R, как правило, удовлетворяется. Расчетная схема рассматриваемой контактной задачи приведена на рис. 7.4.

Кольцо, моделирующее крепь, материал которого имеет дефор­ мационные характеристики E 1 w v 1 (соответственно модуль дефор­ мации и коэффициент Пуассона), деформируется совместно с ок­ ружающей средой, моделирующей массив пород и имеющей харак­

теристики £ 0, v0, т. е. на линии

контакта

L соблюдаются условия

 

 

Таблица 7.1

 

 

№ сечений

 

пвпад

б, см

Л.

 

 

 

 

" V

 

7.1,

а

11

10

6,5

2 - 5

7.1,

б

9

6

5,3

2—5

7.2,

а

11

8

8,2

3—8

7.2,

б

10

11

10,2

2—5

7.3,

а

10

9

14

2—6

Рис. 7.4. Расчетная схема на­ 7.3,

б

11

9

11

2 - 5

брызгбетонной крепи

 

 

 

 

 

непрерывности векторов напряжений и смещений. Внутренний контур кольца L x свободен от действия внешних сил.

Граничные условия поставленной задачи имеют вид:

м (1)

к

М (1)

 

К

М ___

К

М ___

К

на

г .

п 0 ч

Ор

—Op,

TpQ

 

—Тр0,

ир —tip,

IIQ

UQ

 

(7.^)

 

 

<Jp

= 0

,

Тр9 = 0

на

Lu

 

 

 

 

где ар, тре — нормальные

и касательные

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

и&

смещения по нормали и касательной

к контуру

L;

индексами «м»

и «к» обозначены факторы, относящиеся к среде (массиву) и кольцу (крепи).

Полные напряжения

в массиве а“ (1),

а“ (|>,

могут

быть

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

ai0),

 

действовавших

в массиве до

проведения

выработки и выражае­

мых формулами

(7.1), и дополнительных

напряжений ст“, а",

тм

обусловленных

наличием в

массиве

ослабляющей

 

т*</>

его выработки

=а<°>+I Gx,

 

(0) ,

м

г“ <1)

= Т.

+ т“„

(7.3)

<^(,) = а;У

 

1ху

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

(7 .2) вместо а£(|),

Тр0(1) фигурируют дополнительные напряжения

<Jp, Тро, смещения

же рассматриваются только дополнительные.

Тогда граничные условия для отыскания дополнительных на­

пряжений

имеют вид:

 

 

 

 

 

 

м

к

ГО)

м

к

(0)

.,м

к

, м

, к ¥тл г

сГр—

 

—^р »

тр9 — тр0

—тр0 »

^р —^р>

 

—tiQ на L;

 

 

 

ар =

0,

Трв= 0 на Lx.

 

(7.4)

В данной

постановке задачи кольцо

принято замкнутым, т. е.

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

Если вокруг выработки имеется слой пород, ослабленный бу­ ровзрывными работами, или слой пород, закрепленный анкерами, то приближенный учет наличия этих слоев (из-за отсутствия в на­ стоящее время решения контактной задачи для многослойного не­ кругового кольца) может быть осуществлен путем введения так на­ зываемого приведенного или эффективного модуля деформации массива [34], определяемого по формулам (6.8), (6.9). Может быть также приближенно учтено влияние ползучести пород на основе метода переменных модулей, как это предложено в работе [12] и описано в §6.1.

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

(Тр = dp Pit

Up=Upf UQ = UQ

на L;

cfp = 0 ,

Тр0 = О на Li.

(7.5)

Как будет показано далее, алгоритм расчета крепи на равно­ мерное давление подземных вод является частным случаем алго­ ритма расчета на действие собственного веса пород при значении А, = 1.

7 .2 . Решение контактной задачи

Пусть г = © (£) — рациональная функция, осуществляющая конформное отображение внешности окружности радиуса 1 на внешность области, занимаемой кольцом и средой таким обра­ зом чтобы единичная окружность Г перешла в линию контакта L.

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

т

 

Z = = £ a £ r v + ап+11Тя,

(7.6)

V=0

 

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

крепи,

имеющего п х = п + 1 неровностей с амплитудой

(отклоне­

нием от проектного контура) б =

г. Первое слагаемое функции

(7.6),

содержащее ряд с коэффициентами ау (v = 1

, m), осу­

ществляет отображение внешности единичной окружности на внеш­ ность внутреннего проектного (гладкого) контура сечения крепи; для достижения достаточной в практических целях точности можно

ограничиться пг =

5.

Второе слагаемое отражает

наличие неров­

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

аппрокси­

мируются гипотрохоидой.

 

 

Таким образом, функция (7.6) имеет вид

 

 

 

 

п

 

(7.7)

 

 

V=0

 

 

 

 

 

Оу= 0

при

v = m - f 1, m-j-2,

п.

(7.8)

Второй этап построения требуемого отображения состоит в пре­ образовании, при котором единичная окружность Г переходит во

внешний контур сечения L, а окружность радиуса /?х< 1 — во внутренний контур L v Для этого сначала определяется радиус

/?!> 1 окружности, которая при отображении (7.7) преобразуется во внешний контур поперечного сечения крепи. Если обозначить через расстояние от начала координат до верхней точки внутрен­

него контура

поперечного

сечения крепи, а через

А — толщину

крепи, то для

определения

R\ получим

уравнение

(п + 1)-й сте­

пени

 

п+ 1 ~

 

 

 

 

 

 

 

(7.9)

 

 

£

5 у - у+1 = 0,

где

 

v = 0

 

 

 

 

ач при

 

 

 

 

~

|

v = 0,

2,

n + 1,

(7.10)

 

|

a i— (4 + А + б)

при

v = l .

 

 

В качестве R* берется единственный, больший единицы, дейст­ вительный корень уравнения (7.9). Затем осуществляется преобра­ зование, при котором окружность единичного радиуса Г в плоско­

сти переменного | переходит в окружность радиуса R\ в плоскости При этом окружности единичного радиуса будет соответствовать

внешний контур, а окружности радиуса R x ~ —!—< 1 — внутрен-

ний контур поперечного сечения крепи. Отображающая функция (7.7) принимает вид

2= (,)©=*>

( g + i : ^ r vy

(7.11)

 

V

V=1

/

 

R==~ 1T ’ <7v= 0v+

f l+l

(v = l,

2,

n). (7.12)

В отображающей функции (7.11) отброшен член q0, влияющий лишь на положение начала координат.

Вследствие симметрии контуров относительно вертикальной оси Ох все коэффициенты отображающей функции вещественны.

Введем комплексные потенциалы (2), ф< (г) (i = 0, 1), регу­ лярные в соответствующих областях S0 (бесконечная область, ог­ раниченная контуром L) и S i (область, занимаемая кольцом, огра­ ниченная контурами L и Lx) и обращающиеся в нуль на бесконеч­ ности, связанные с компонентами напряжений известными форму­ лами Колосова—Мусхелишвили. Тогда поставленная задача тео­ рии упругости сводится к краевой задаче теории аналитических функций комплексного переменного при граничных условиях:

Фо (t)+t<?o (0 + фо (0 =

_____ _____

t

= Фх М+ <ф'| (0+ Фх (0-

i I (*п0) + д а ) ds + c на L;

 

to

— фо(0

1

[ftpi (0 + ^о (0] = — <Pi (О

Цо

М'О

 

Pi

 

 

4Pi (0 +

1 (0 + Фг (0 = OL

Здесь

[<ф1 (t) + 'Фх (О].

Pi

(7.13)

на Ьг.

 

^

= Н 7 7 Г Г Т ’

х / —3 —4vf (i =

0, 1),

(7.14)

 

 

 

Z (1 -f- Vi)

 

 

 

Ei,

v,- — модули

деформации

и коэффициенты

Пуассона

материа­

лов

сред Si (i

=

0, 1).

 

 

 

Определим входящее в правую часть первого уравнения (7.13)

выражение

tj (АГ^0) Н-

c/s на L. Из

(7.1) следует,

что

 

Х „0)( = а *1cos (л,

х) + T(°JxCOS (л, у) = — уНа

cos (л,

х),

Kn0>= txj cos (л,

х) + o f cos (л, у) =

—ЯуНа

cos (л,

г/). (7.15)

Тогда, учитывая,

что

 

 

 

 

 

 

 

 

 

cos (л,

л:) =

as

,

соэ(л,

у) = ---- (7.16)

 

 

 

 

 

 

 

 

 

as

 

 

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д а

+ д а ) ds =

— уЯа* (dy — хЯ dx).

 

(7.17)

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

i

J(X<0) + iY(„0)) d s=

—/уЯа* J (dy—t7.dx) =

 

= -iy H a * [ ( * L +

 

2

+ i J L _ i J L - i x _Ё*__а_Ё*. + Я -ЁХ- -

r

J \ 2

 

 

2

2

2

 

 

2

2

 

to

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

~ я “f " ) = — 1УН а * J [ ~ Y

(dY+ ldx) + y

(dy—tdx) +

 

 

 

 

 

to

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

+ у

(dy —■idx) —

 

(dy + idx)] = — iyHa* J

 

 

(dy —idx) -f-

 

 

 

 

 

 

 

t

to

 

 

 

 

 

1 — A, (dy + idx)J =

 

 

 

 

 

 

 

— ytfa* ^

1^ ^

(dx + idy) —

 

2

 

 

 

 

 

 

 

 

 

 

 

 

\ - X

(dx—idy)] = -

у Ha* J

dz—

 

dzj =

(7-18)

Таким образом, первое граничное условие (7.13) на контуре принимает вид

Фо (0 + *фо (0 + Ы О = Ф1 (0 + |*ф! (О + ФЛО + УНа*

1

1

(7.19)

 

Записывая граничные условия (7.13) с учетом (7.19) в преобра­ зованной области и произведя над ними операцию комплексного сопряжения, имеем на Г

 

 

(i)...

 

 

 

 

R

 

 

 

й

' 1

 

 

 

 

 

й ' 1

 

ф° 0 ~ )+

тЧоч

Фо (а) + th, (а) = Фх

V о )

,

U

;

ф!(а) +

W

 

 

 

 

 

со'(а)

 

 

+ ФеИ

+ УНа* [- L ± i-

со

___ L z A

СОи ].

(7.20)

Ф °(~ )'

 

 

(т) .

 

= —

Ф1 ( — ) -

 

L

'

фо (Ст) -I--4>,(<т)

щ>

\

о /

 

(а)

 

 

 

 

М-1

\

а /

 

 

 

 

_1_

5J' 1il

ф1И + ФхИ

 

 

(7.21)

 

 

 

 

Mi

со' (а)

 

 

 

 

 

 

 

 

Ф1

 

 

со' (RiO)

Ф1 С ^И Ф х (Яха) = 0,

(7.22)

 

 

 

 

 

 

 

 

 

 

 

где а =

е10 — точка единичной окружности Г.

 

 

 

Решение

данной

краевой

задачи может

быть

получено, как

и в работе

[41 ], на основе развития метода М. П. Шереметьева, но

с тем отличием,

что в работе [41] окончательный алгоритм и про­

грамма

расчета для ЭВМ построены для

гладкого

проектного кон­

тура (п = 4) здесь

же

решение,

алгоритм

и программа даются

в общем виде при произвольном числе п членов ряда отображающей функции (7.11); кроме того, ввиду отличия граничных условий (7.20) — (7.22) от рассмотренных в [41 ], здесь система уравнений будет иметь другие свободные члены; изменяются также и формулы для напряжений.

Отметим, что примененный в настоящей работе путь решения задачи отличается от метода М. П. Шереметьева и основан на пред-

. Ч т )

Ч ^ )

в виде рядов, для опре­

ставлении выражении —

 

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

кованных в работах [41, 43] и приведения их к виду, изложенному в § 6.2.

Согласно работы [43], умножим граничное условие (7.20) на Xi/pi и вычтем из него условие (7.21), затем умножим (7.20) на l/p.x и сложим с (7.21). После деления всех членов полученных равенств

на выражение 1/р х граничные условия примут вид

 

 

 

 

 

соR .

 

 

 

 

 

 

 

- (

1

 

 

 

 

 

s<p0( ± ) + I _

со' (о)

фо (ff) + Фо (а)

 

 

- /

1

 

 

 

 

 

 

 

 

(0( - )

Ф1 (О) + % (О) + А х(а) на Г;

 

(7.23)

 

° j

 

со' (а)

 

 

 

 

 

 

 

 

% ( - y ) + d

0) ( т )

 

 

= 9x( - L ) +

,42(a),

(7.24)

-

■ фо(ог) + ^о(сг)

 

(O'(а)

 

 

 

 

 

 

 

 

 

- (

*i

 

 

 

 

 

 

^ У а ) ^

(— )

ф ! ( а д + %

( а д = о ,

(7.25)

<D'(*ior)

где

 

 

 

 

 

 

 

 

 

j __ _1_+_^оР_

^

_ 1— Р

 

о _

Е1 (1

Уо)

(7.26)

1+ х,

 

1+Xi

р “

£o(l+Vi)

 

 

 

1=1— d,

 

s = l — t,

 

 

 

а функции A x (а), A 2 (а), выраженные в долях величины

у Н а *

1+Xi

имеют вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аг (а) = хх

 

 

 

 

 

® (°)]«

 

 

 

 

 

 

 

 

 

 

(7.27)

Комплексные потенциалы <р0 (£), ф0 (£)> характеризующие на­ пряженное состояние массива, голоморфны вне единичной окруж­ ности Г, включая бесконечно удаленную точку, т. е. они имеют при достаточно больших £ разложения

Фо (£) = Е а £ - \

ь (6)= Е bvr v .

(7.28)

v=0

v*=0

 

Не меняя напряженного состояния, можно положить а0 — 0. Комплексные потенциалы <рх (£), фх (£), определяющие напря­ женное состояние крепи, голоморфны в области кольца < р < 1

и могут быть представлены в виде сумм

Ф1 © = P i (6) + Рг © , Фх (1) = Qi © + Q2 © ,

150

где Р г (£), Qi(b) — функции, голоморфные внутри единичной ок­ ружности Г и, следовательно, разлагающиеся в ряды

P i (I) = £ с Х ,

Q1(6) = £ d , r ,

(7.30)

v=0

v=0

 

a P 2 (£). Qz (£)—функции, голоморфные вне окружности радиуса Rlt включая бесконечно удаленную точку, выражаемые разложениями

Р г(I) = £ a ~ v.

Q2(S) = £ f vl~ v-

(7.31)

V=1

V=1

 

В силу геометрической и силовой симметрии относительно оси коэффициенты разложения искомых функций в ряды (7.28), (7.30) и (7.31) вещественны.

Умножим равенство (7.24) на ядро Коши —?----- —— и проин-

а — £

тегрируем его с учетом представлений (7.28) — (7.31) по контуру единичной окружности дважды, считая точку £ последовательно расположенной вне и внутри Г.

Результаты интегрирования отдельных выражений, входящих в 7.24, полученные на основании использования свойств интегралов

типа Коши, приведены в работе (41 ].

Остается проинтегрировать

по окружности Г умноженное на

ядро Коши

выражение А 2 (а),

стоящее в правой части граничного условия (7.24).

 

Учитывая,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ШИ

= R (of +

X

4k<J~*) .

 

(7.32)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.33)

из выражения

(7.27)

имеем (в долях

у Н а *

 

>

 

R (l+ x i)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ni

J

а — £

 

 

 

 

 

 

1+ X-S-1 +

- L= ^ Z < 7 * |- ‘

при I вне

Г

 

 

 

 

2

 

 

 

2

£=х

 

 

 

(7.34)

 

 

1

 

п

 

чАк

1— Я. |

 

 

 

 

 

 

 

при |

внутри Г

 

 

 

 

Z

 

 

 

 

 

2

 

k=\

 

 

 

 

 

 

 

 

После интегрирования условия (7.24) получим

 

 

Г

i f i )

 

 

 

 

 

 

 

 

 

- P J - L ) + cо +

d - V

\

(O'

Л

' Уо И

a — £

(1) +db0=

J

(a)

 

 

 

 

 

 

 

\ l

/

 

г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ - Ц

^ -

 

£

ч А ~ к -

2

 

l

вне Г

(7.35)

 

 

 

 

2

 

 

k=\

 

 

 

 

 

 

Чт)+

 

 

 

-/

1

'

, 4

 

 

 

 

d

 

 

 

(

a )

 

d a

■№„-Р.(у) +

2 n i

J

 

со'(a)

фо (<*)

0 —5

 

 

 

г

 

 

 

 

 

 

 

 

 

+ * - - Ц

±

- б +

- ^ £ «

* 1

 

1—% £. £ внутри

Г (7.36)

 

 

 

 

2

k=i

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

- /

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ч т ) " ‘Ф .(у) + i

 

 

)

©(тЪ

 

- c . -

 

 

- 7

^

-

1 » W - ^ - + A

 

 

 

 

 

5

внутри Г

 

 

 

 

 

 

 

- Ц ^ - Е и б Ч - Ц

1^

I,

(7.37)

 

 

 

 

2

fc=i

 

 

 

 

 

 

* » -!= -) 4

^Г * w

+6"+T Pl( f ) '-f

+

 

Г

 

 

 

 

 

 

 

 

 

 

 

5 вне

Г

 

 

 

 

 

 

 

 

 

 

 

+

d

L

 

2

 

 

2

 

*=i

J

(7.38)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- /

1

 

 

 

 

Представляем

отношение

CiO

i l

в

виде

 

 

 

 

 

 

 

 

(O'(О)

 

 

 

 

 

5(т)

=

Е

М Ч

Е

 

 

(7.39)

 

 

 

.

 

 

 

 

 

со' (a)

 

 

А=о

 

 

ft=l

 

 

где коэффициенты hk и /г_* определяются по рекуррентным форму­

лам, полученным делением

многочлена

на

многочлен:

 

 

л-<*+1)

 

 

 

 

 

 

(7.40)

hk —Як+ S*+1

2

V7V^V+A+I

(k =

n>

0),

 

 

v=l

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

h-k =

2

 

л+1

(^=

1>

2,

°°).

(7.41)

где

V=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f 1

при

6 = 1

 

|

1

при

k<_n,

(7.42)

\ 0

при

k Ф 1

— ( 0

при

k > n.

 

Тогда, учитывая вытекающую из (7.28) формулу

 

 

 

Фо (<т) =

YJ vava - v- ‘ ,

 

 

V=1

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