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

книги из ГПНТБ / Булычев, Н. С. Расчет крепи капитальных горных выработок

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

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

При 6 — 1м , lL = 0,3901807? формулы для перемещений имеют следующий вид [9]:

Е1Ьп ^ 28,304526 ■ El

-5,576335

FR

1,151548/?;

К^НЗ

 

 

El

Я/614 = -70,660031 — 7 , 8 3 9 3 6 8 - ] - 0,332877/?;

 

 

 

E l

 

0,502896 FR

(14.10)

Е 1 \ ь = 40,595367 KMR3

ЕІ8Ы- EI8l7 = Я/ö,

 

 

 

 

 

E I

 

 

1

2Я/61В= 0,502896 A'M R3

1 0,502896 FR

ЕІ8Ш=211,013074

 

EI

13,550159 FR

0,783898/?;

KMR3

£ / 6« = -170,754995

E I

 

 

 

I

 

0,130060/?;

K M R 3 -1,111324

FR

EI8W = 33,687286 i S

r -

 

 

 

/

 

1-« « M T S - >

Е18„=Е18Л 2£/649= -

1,314136

E l

 

 

 

/

KMR3

 

1,314133 *

 

 

 

EI

0,405619

 

 

0,520240/?;

£/S 55 = £/Ö88= 224,506586 K ^ R 3

 

 

 

 

 

 

 

 

 

 

 

E l

Е18Ъ9= £ /ö 67= EI87&= EI8S9 = -129,348384 KMR3

4-0,202809 FR

0,130060/?;

 

 

ЕІ8Ъ7 = EI8es

£ /ö79= 35,001422

E l

,

Ä(i)i?3

EI866 = EI877 = 2EI899 =

189,505163

E I

 

RMR3

1

f 0,405619 ~

+ 0,520240/?;

658

^59 = ^69 = 0;

EIAlp = -10,989885

EI

p _

0,03188857/?3p -

5,531740 -E Rp;

£/A5p = -8,200079 K^R

p -

0,951021 -j- Rp;

 

 

 

EI

 

 

 

 

 

(14.11)

 

 

 

 

 

 

 

Rp;

EIAep = -0,865051 KMR - p — 0,966459

 

 

 

 

 

PT

 

 

 

j

EIA7p= ElABp = 2EIA9p = —0,966459 -p M R - P 0,966459 — Rp-

99

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

61к, (14.12)

где Ог — площадь эпюры изгибающих моментов состояния і.

7. Определение неизвестных моментов в узлах расчетной схемы путем решения системы канонических уравнений (14.2). Решение производится последовательным исключением неизвестных при по­ мощи алгоритма Гаусса.

Полная величина изгибающих моментов в точках 2 и 3 опреде­

ляется по формулам, смысл которых не требует пояснений:

 

М г = 0,876692717! + 0,123308Л74-0,020598Я2р;

 

М 3 = 0,525539л/! г 0,47446Ш 4-0,047511Я 2р.

(14‘1Л)

Эпюра изгибающих моментов показана на рис. 45.

8. Вычисление окончательных значений опорных реакций и нор­ мальных сил производится суммированием усилий в основной си­

стеме от действия нагрузки и найденных моментов:

 

 

 

 

 

 

Ri -- Rip

MkRik,

 

 

 

(14.14)

 

 

 

N i ^ N ip- r ^ M kN ik.

 

 

 

 

 

 

 

 

 

 

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

[9,

261.

 

К сказанному необходимо добавить, что зона отлипания опре­

деляется

по

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

 

 

 

приближений,

при этом

необходимо всякий

 

 

 

раз повторять решение канонической системы

 

 

 

уравнений (14.2).

 

 

 

 

 

 

 

 

Метод

Метрогипротранса

был применен

 

 

 

Л. Б. Капланом [84]

для

расчета крепи

 

 

 

стволов,

которые

проходились

бурением в

 

 

 

Подмосковном

бассейне.

При

расчете было

 

 

 

принято,

что

нагрузки

на

крепь

ствола

 

 

 

глубиной

100 м и диаметром 4,5 м (в

свету)

Рис. 45. Эпюра изгибающих

складываются

из

равномерного

гидростати­

моментов в крени от дейст­

ческого давления

и давления пород,

часть

вия активной

веотикальной

нагрузки

 

которого

распределена равномерно, а часть

ние (рис. 46).

Расчет

представляет собой «одностороннее»

давле­

на «одностороннее» давление ведется с учетом

упругого отпора пород (рис. 47).

Методика С. А. Орлова. Метод Метрогипротранса получил разви­ тие в работах С. А. Орлова, который исследовал взаимодействие крепи, нагруженной активными нагрузками, с упругой средой. Расчетная схема крепи, как и в методе Метрогипротранса, при­ нимается в виде многоугольника, прикрепленного к массиву всеми вершинами, за исключением зоны отлипания А (рис. 48, а). За лиш­ ние неизвестные принимаются изгибающие моменты в верхней точке

100

и в местах опор. В точках опирания многоугольника на контур кру­ гового выреза в упругой плоскости при действии на него активной нагрузки возникают сосредоточенные реакции, которые G. А. Орлов условно называет упругими опорами [135].

{ШШШШШК

рг --ZOtcJm2— р-ШШИШПТТ|

 

- 60тс/нг \

[

'Рг -ZOtc/ m?

 

 

p,-Z0rc/n

 

 

Pl*Pr ,

 

 

=80тс/мг

I

I у

рг -ZQTcfn2

^шшпппУ

Рис. 46. Расчетные нагрузки на крепь ствола:

о — гидростатическое давление; б — давление пород; в — суммарные нагрузки

Лишние неизвестные определяются решением системы канони­ ческих уравнений (14.2), при этом в формулах для перемещений (14.3) изменяется последний член:

бik

MiMu

Д1 ,

'S ?

N iN k ,

(14.15)

~Ë1

 

 

г _ г 2 j

EF

 

 

 

 

где ut — радиальное смещение упругих опор.

 

Смещения определяются

сле­

 

иниипдиитимнмтінт

дующим образом. С. А. Орловым

 

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

 

 

 

контура кругового выреза в упру­

 

 

 

гой плоскости,

характеризуемой

 

 

 

центральным углом 2ß, приложена

 

 

 

единичная равномерная радиаль­

 

 

 

ная нагрузка.

Смещения

точек

 

 

 

контура в этом случае определя­

 

 

 

ются выражением в виде

триго­

 

 

 

нометрического ряда

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2пО X

 

 

 

 

 

 

X °’5+т 2

2 [А- p (А + 1)] +

1

X

 

 

 

к 21)

 

 

 

 

 

h=2

 

 

 

 

 

ЦЦЩЦЦЩЦЦЦШНННЦНШШ

 

 

 

 

 

 

X sin /cßcos АО

(14.16)

Рис. 47. Расчетная схема крепи ствола при

 

 

 

 

 

 

действии «односторонней» нагрузки

Для получения достаточной точности необходимо при расчете учиты­ вать не менее 16 членов ряда.

Эпюра

упругого

отпора пород принимается ступенчатой

(рис. 48, б),

середина

каждого участка соответствует узлу расчетной

101

схемы крепи. Смещения опор рассматриваются как суммарные от действия рассматриваемых нагрузок:

иі = 2

R iu Іѳ-anß,

(14.17)

П

 

М { = 1.

где Ri — реакции опор от действия единичных моментов

Методика С. II. Веркиной.

Существенное усовершенствование,

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

Рис. 48. Расчетная схема к методике С. А. Орлова:

а — расчетная схема крепи; б — эпюра упругого отпора пород

появляется возможность представить его в компактной форме и полностью механизировать [24].

Рассматривается произвольный ломаный стержень (рис. 49), одна из главных плоскостей инерции которого совпадает с плоскостью чертежа. Нагрузка приложена в узлах. Предполагается, что все кинематические и силовые факторы в начале координат (начальные

параметры)

имеют вполне определенные заданные значения: С/0,

Р 0, ф0, Х 0,

У 0, М 0. Перемещения точек оси зависят от изгиба и осе­

вого сжатия стержней. В узле 1 проекции перемещений на коорди­ натные оси составляют:

и 1 = и 0- (фсА + м о 2^ - + Qo е й г ) sin «1 — ^0 - щ - cos “ В

Vi = Vo + (% h + Mo l§ j^ + Qo-öWF[) cos ai ~ N o~gJr[ silia6 (14.18)

Ф! = фо + М 0 + Q0 2^ - .

где

Q0= Уо cos а х — Х 0sin а х + R0cos (у0—ocj; N 0 = Y 0sin + X 0cos a x — R 0sin (y0— a t);

= K0(V0cos у,, - U0sin Yo);

102

К 0 — жесткость опоры в нулевом узле.

При последовательном переходе от одного узла к другому, учи­ тывая условия равновесия и совместности деформаций, получены общие формулы метода начальных параметров для всех шести расчет­ ных величин (Un, Ѵп, <р„, Х п, У„, Мп) в произвольной точке крепи:

и п= и 0АѴЬ + Ѵот + Ф о ^ + X 0AWx + Y 0A^y + М 0АШі + A%';

(14.19)

........................................................................................4 ........................

Mn= UQAmu + V0A(ffv -j- %А%\ + Х 0А (мХ + Y 0A (ffy + M0A(mm + A%\

Для определения коэффициентов указанных линейных преобразо­ ваний получены рекурентные зависимости [24].

♦ui и тгп *

»01

Рис. 49. Расчетная схема к методу начальных пара-

Рис. 50. Расчетная

схема крепи

метров

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

 

симметричной относительно вер­

 

тикальной оси

нагрузке

При расчете крепи выработки круглого сечения рассматриваются две системы координат (рис. 50). Уравнения сопряжения двух систем в общей для них точке 4 имеют вид:

UА1- \ - =

V41 Ѵ42 = 0 ; ф 41 + ср42 = 0 ?

 

* 4 1 -^4 2 = 0;

У41 + г 42 = 0; м 41- м 42 = 0.

(14.20)

При симметричной относительно вертикальной оси нагрузке опре­ делению подлежат шесть начальных параметров: £/01, М 01, У 01,

^ 02» М 02> Y 02-

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

103

Метод С. С. Давыдова [54, 150] основан на предположении, что вокруг выработки образуется активный или упругий слой пород, за пределами которого массив рассматривается как жесткий (рис. 51). Под действием собственного веса активного слоя пород возникают вертикальные и горизонтальные составляющие объемных сил, кото­

рые уравновешивают часть распора крепи,

благодаря чему умень­

 

шаются горизонтальные

уп­

 

ругие

деформации

системы.

 

Толщина

упругого

слоя

Н г

 

(рис. 51) определяется по

со­

 

отношению

между

горизон­

 

тальными

 

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

в

 

нетронутом массиве и распо­

 

ром свода, действующим на

 

верх боковой стенки. Тол­

 

щина упругого слоя в осно­

 

вании

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

по

Рис. 51. Расчетная схема к методике С. С. Давы­

формуле

Н. А. Цытовича,

учитывающей

ширину

по­

дова:

1 — неподвижный слой; 2 — объемные силы; 3

дошвы стенки и соотношение

давление пород; 4 — упругий свод; 5 — жесткая

между нею

и

длиной выра­

или упругая стенка

ботки.

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

Боковые стенки рассматриваются как жесткие или как упругие,

в зависимости от величины

 

 

 

 

a = JT T - Т

^ ( ° ’2/*ст)3.

О ^ І )

где Е о,

р,0 — обобщенные модуль

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

поперечной де­

формации пород;

hCT — высота стенки.

 

При а

0,05

стенка считается жесткой.

составляющие

На боковые

стенки действуют

горизонтальные

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

Кроме объемных сил на боковые стенки действует упругий отпор пород, определяемый из расчетной схемы (рис. 52). Отпор модели­ руется стержнями; предусматривается по пять стержней вдоль боковой поверхности и подошвы стенки. Таким образом, в основной системе имеется 15 или 16 неизвестных (в зависимости от конструкции лотка), в том числе 5 усилий Х п действующих по боковой плоскости стенки; 5 усилий Y t по основанию, усилие X , заменяющее влияние

104

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

Системы уравнений для определения X t и Y t решены примени­ тельно к основным случаям практики, а результаты сведены в табл., вследствие чего расчет стенки сводится практически к статически определимой задаче.

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

Ap(hcosff)

Рис. 52. Расчетная схема стей­

Рис. 53. Расчетная схема крепи выработки круглого се­

ки выработки

чения но С. С.

Давыдову:

 

1 — сила, уравновешивающая

вертикальную нагрузку

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

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

на рис. 53

[54]. В качестве расчетных нагрузок приняты вертикаль­

ное горное

давление р х +

р 2 + Ар,

собственный вес крепи, гори­

зонтальные

составляющие

объемных

сил активного слоя пород р3

и упругий отпор пород р4. Ордината эпюры упругого отпора пород'р4 определяется из анализа перемещений точек крепи, лежащих на горизонтальном диаметре.

Задача рассматривается как дважды статически неопределимая. Расчет производится для каждого вида нагрузки (вертикальной и боковой) отдельно. Полные усилия в крепи находятся путем сум­ мирования усилий от каждого вида загружения.

1 0 5

Приближенные методики расчета

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

Рис. 54. Приближенные расчетные схемы крепи выработки круглого сечения (табл. 2,3)

совместности деформаций крепи и пород в точках, лежащих на гори­ зонтальном диаметре. В табл. 23, заимствованной из работ [89, 150], приведены характеристики некоторых предложений по приближен­ ному учету отпора пород (рис. 54).

Т а б л и ц а 23

Автор

Форма эпюры отпора пород

Направление отпора

А. М. Новиков К. А. Вахуркин В. С. Эристов

Л. И. Горелик

Г. К. Клейн

М. Шпанглер М. Друккер

Эллипс (рис. 54, а)

 

Горизонтальное

Квадратная парабола (рис. 54, б)

»

Подобна эпюре горизонтальных пе­

»

ремещений от вертикальной на­

 

грузки (рис. 54, б)

 

 

Подобна эпюре радиальных переме­

Радиальное

щений (рис. 54, в)

 

 

Подобна эпюре горизонтальных пе­

Горизонтальное

ремещений от вертикальной на­

 

грузки (рис. 54, г)

(рис. 54, д)

»

Квадратная парабола

Кубическая парабола

(рис. 54, е)

»

106

Методика О. Е. Бугаевой. Эта методика получила большое рас­ пространение, особенно для предварительного расчета крепи, так как дает достаточно близкую сходимость с методикой Метрогипротранса. Рекомендуя приближенные методики расчета, О. Е. Бугаева исхо­ дила из положения, что точность метода расчета крепи должна соответствовать степени точности исходных данных (нагрузок на крепь, характеристик отпора пород и др.).

О. Е. Бугаева предложила опре­ делять эпюру отпора пород по не­ скольким точкам, которые соеди­ няются плавной кривой [29, 79]. Разработаны приближенные мето­ дики расчета ряда конструкций крепи (замкнутых и незамкнутых). Наи­ большей известностью пользуется ме­ тодика расчета крепи выработки круглого сечения (рис. 55).

Размер зоны отлипания прибли­ женно принят постоянным, харак-

ji

теризуемым центральным углом — .

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

сечениях Ѳ

± — ,

двумя

на

55.

Расчетная

схема к методике

уровне горизонтального

диаметра

 

О. Е. Бугаевой

координатами

Ра =

‘иА и точкой

Очертание

эпюры отпора

с координатой рв = к и 'иВ — в подошве.

пород описывается выражениями:

 

 

 

 

- р . cos20

/ я

ѳ

- л

(14.22)

 

 

 

 

 

 

 

 

= 2 ■):

 

Р = Р А sin20 -f р

COS'

я

ѳ =

(14.23)

 

і2Ѳ

Т

 

 

 

 

 

 

 

 

 

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

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

а) от вертикальной нагрузки р:

М = pRR1 [Ат !-В + Сп(і т)\\

(14.24)

N = pR x[Dm г F Л Gn( 1 + m)],

где

т = 2

Ял

Я

 

 

1

0,06416-1-

Е І

 

Ь — ширина кольца крепи;

 

107

н

А

В

С

D

F

0

0,1628

0,0872

-0,00700

0,2122

—0,2122

л /4

-0,0250

0,0250

-0.00084

0,1500

0,3500

л /2

—0,1250

—0,1250

0,00825

0,0

1,0

Зл/4

0,0250

-0,0250

0.00025

-0,1500

0,9000

Л

0,0872

0,1628

-0,00837

—0,2122

0,7122

б) от собственного веса крепи g (тс/м2):

M = gR*(Al + B ln)-

 

 

 

N = gR {Сг+ Dpi).

 

 

 

(14.25)

 

 

 

 

 

 

 

 

 

Значения

входящих в

формулы

 

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

приведены

в табл. 24.

 

 

 

 

 

 

 

 

 

 

 

 

Аналогичные формулы получены и для случая гидростатического

давления на крепь [29,

79].

 

 

 

 

 

 

 

 

 

 

Интересная методика расчета предложена Г. А. Скобенниковым

(Труды ЛИИЖТ, вып.

241,

М. — Л-, «Транспорт»,

1965).

 

 

 

 

Методики расчета

крепи

стволов

 

 

 

 

 

Методика

немецких

инженеров.

В 50-х годах немецкие ученые

и инженеры (Ф. Мор,

Г.

Линк,

Ф.

Ванслебен), развивая

идеи

О.

Домке, разработали

методику

расчета крепи

стволов, которая

 

 

 

 

 

в дальнейшем нашла отраже­

 

 

 

 

 

ние

 

в официальных

рекомен­

 

 

 

 

 

дациях угольного объединения.

 

 

 

 

 

В

основу

расчета

положено

 

 

 

 

 

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

о

дополнитель­

 

 

 

 

 

ных

неравномерных

нагрузках

 

 

 

 

 

(рис. 56), действующих на крепь

 

 

 

 

 

помимо основных

равномерных

 

 

 

 

 

нагрузок

со

стороны пород и

 

 

 

 

 

воды. Общая нагрузка, согласно

 

 

 

 

 

Ф. Мору, описывается выраже­

 

 

 

 

 

нием

(13.5),

а

максимальная

 

 

 

 

 

ордината

дополнительной

на­

Рис.

51). Схема к

расчету крепи

по методике

грузки р 2 составляет 5—10%

 

 

Г. Линка

 

 

от основной равномерной на­

 

Механизм появления

 

 

грузки.

 

 

 

нагрузок

 

дополнительных неравномерных

Г. Линк иллюстрирует схемой, согласно которой неравномерность

нагрузок

вызывается неравнокомпонентностью

поля напряжений

в массиве

пород [236]. Г. Линк отмечает, что

деформации крепи

 

 

 

 

Т а б л II ц а 24

G

А 1

в,

С,

I)t

0.02100

0,3447

-0,02198

-0,1667

0.06592

0,01485

0,0334

-0,00267

0,3375

0,04661

0,00575

—0,3928

0,02589

1,5708

0,01804

0,01380

-0,0335

0,00067

1,9186

0.04220

0,02240

0,4405

-0,02620

1,7375

0,07010

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

В 1953 г. Ф. Ванслебен предложил определять степень неравно­ мерности нагрузок на крепь ствола, исходя из предположения о не­ возможности перемещений крепи в сторону массива пород [271]. При расчетной нагрузке на крепь, определяемой формулой (13.5), коэффициент неравномерности нагрузок

3El

0) - £ і

EF

'AEI

Pol*3

 

(14.26)

Ро

 

 

 

 

H

ЗроК

LA+'3

 

 

 

 

 

пород

 

тс/м2);

 

0

где p ’v = p v ~~; pv — давление

{ р ѵ ^ 0 , Ш

р

основная равномерно распределенная

нагрузка

(р 0 = р ' ѵ

+

рв);

Рв — Рв~]Р’ Яв — радиус

поверхности

приложения гидростатиче­

ского давления.

Для такого же условия (отсутствие перемещений крепи в сторону

массива пород) Г. Линк получил более простую формулу

[232]:

со

_________ Ив_________

(14.27)

 

 

 

( « в + 0 ,3 /0

( - і і г ~ 0,5)

 

При аппроксимации нагрузок на крепь выражением

 

 

р ---р0 *-р2cos20

(14.28)

указанное условие приводит к соотношению

 

 

Рг _

_91_

(14.29)

 

Po

F W

 

Отмечая, что в полученных выражениях степень неравномерности нагрузок на крепь зависит только от характеристик крепи и не зави­ сит от пород, Г. Линк считает это маловероятным. В работе [232]

108

109

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