ГЕОМЕХАНИКА
УДК 622.833.5
В.М. Волков, А.Е. Клыков ВЛИЯНИЕ ДИЛАТАНСИИ ПОРОД НА СМЕЩЕНИЯ В ВЫРАБОТКЕ
Вокруг выработки радиуса Я в гидростатическом поле с составляющей нетронутого массива д при прочности пород на одноосное сжатие осж < 2д формируется зона неупругих деформаций (ЗНД) радиусом Яр , точки которой на границе с упругой зоной находятся в предельном состоянии. В этом состоянии связь между главными напряжениями для прямолинейной огибающей предельных кругов Мора имеет вид:
Ог С(Гд + Осж , (1)
где осж - прочность на одноосное сжатие пород, а с = (1 + sinp)/(1 - sinp).
В упругой зоне напряжения связаны зависимостью иг + ид= 2, которая представлена прямой 3 на рис.1. На большом удалении от контура выработки напряжения равны д , что соответствует точке А . Уравнение (1) предельного состояния описывается прямой 1, точка пересечения которой с прямой 3 соответствует границе ЗНД с упругой зоной. Эта граница формирует вокруг выработки поверхность радиусом Яр. Внутри нее с уменьшением координаты г от радиуса Яр, т.е. с приближением к контуру выработки, напряжения уменьшаются. Если на диаграммах объёмного сжатия осадочных пород при предельном состоянии есть горизонтальный участок, то вблизи границы ЗНД будет наблюдаться зона предельных напряжений, соответствующая участку МК на прямой 1. С этой зоной соседствует зона, напряжения в которой соответствуют участкам спада на диаграммах. Эти напряжения меньше предельных и для них точки располагаются на некотором участке КЬ ниже прямой 1. В непосредственной близости от контура выработки породы находятся в запредельном состоянии, и если выработка не закреплена, то на контуре напряжение равно остаточной прочности ст0 при одноосном сжатии. Этому напряжению соответствует точка N на рис.1, а участок ЫЬ запредельного состояния описывает зависимость остаточной прочности пород в приконтурной зоне выработки от минимального напряжения. До настоящего времени вид этой зависимости не установлен, однако эксперименты показывают [1] возрастание максимального напряжения от минимального и во всех известных решениях задачи о смещениях в выработке с учетом дилатансии угол
внутреннего трения пород на стадии остаточной прочности принимается равным р для предельного состояния. В этом случае участок ЫЬ параллелен прямой 1.
Если при решении задач о напряженном состоянии окрестности выработки принимать во внимание участки ЫЬ, ЬК и КМ, то это сопряжено с определенными трудностями, заключающимися в необходимости определения положения точек К и Ь. Задача значительно упрощается при замене этих участков отрезком МЫ прямой. Это упрощение представляется приемлемым для практических задач. Кроме того, оно позволяет рассмотреть многие аспекты, неразрешимые в том случае, если учитывать полные диаграммы объёмного сжатия, аппроксимируя их каким-либо образом.
Тогда связь между главными напряжениями для ЗНД примет вид:
Ф2 = кФг + Ф0. (2)
Точка М не является фиксированной, поскольку представляет собой пересечение прямой 1, характеризующей прочность пород, и прямой 3, положение которой зависит от напряжений нетронутого массива. С ростом составляющей д при одинаковых свойствах пород точка М вместе с прямой 3 перемещается вверх и наклон прямой МЫ уменьшается, в связи с чем уменьшится и величина коэффициента к.
Таким образом, коэффициент к в (2) яв-
СТе
1 к X
2ц / А
с\ /і
Чж N ч
СТо :
Рис.1. Графическое представление напряженного
состояния пород за контуром выработки
ляется параметром, характеризующим состояние массива в окрестности выработки.
Решение задачи с физическим уравнением (2) для ЗНД при осесимметричном распределении напряжений вокруг выработки радиусом Я = 1 с реакцией крепи р приводит к следующим результатам [2]:
_ _ [(к- 1)р + ор]гк-1 -а0 ; г к-1 ’ к[(к - 1)р + ао]гк -1 - а0
аа = ■
Яр _
к-1
(Ц-а 0)(к-1)/(к + 1) + а 0
1/(к-1)
(3)
(к-1)Р + а о
к = [2яс + <УсЖ - с>о(с + 1)] / (2я - осж); ог (Яр) = 2 - Осж)/(с +1).
Смещение пород внутрь выработки обусловлено перемещением границы упругой зоны и дилатансией пород в ЗНД. Радиус ЗНД и напряжение на её границе можно определить по формулам (3), а затем найти упругую составляющую смещения контура выработки [2]:
л2
У
Я 2
/Я _^~[2(1-/и2)- - (1 + н)аг(Яр)/я] Е
Я
Я
2
(4)
где Е и ц - соответственно модуль упругости и коэффициент Пуассона породы в окрестности выработки.
С переходом в предельное и запредельное состояние происходит увеличение объема пород в ЗНД, дилатансия, что приводит к дополнительным смещениям контура выработки.
На основании экспериментальных данных [1] с учетом граничных условий можно принять относительную объемную деформацию гу в ЗНД связанной с минимальным напряжением зависимостью:
0 —г (Яр )—г
^ = ^ (Я ) ^------, (5)
—г (Яр ) + -г
где £у° - объемная деформация в запредельном состоянии при одноосном сжатии.
Увеличение объема ЗНД связано со смещением ид контура выработки, которое определится интегрированием:
I
(к - 1)ог (Яр )-[(к - 1)р + оо]г - + оо (к - 1)ог (Яр ) + [(к - 1)р + оо ]ГЫ - оо
г&.
(6)
Для вычисления этого интеграла можно, например, применить квадратуру Симпсона.
Расчеты показали, что дилатансионная составляющая смещения (6) на порядок больше упругой, рассчитываемой по формуле (4).
Полученные формулы реально отражают влияние прочностных свойств горных пород, начальных напряжений в массиве и сопротивления крепи на смещения в выработке. Уменьшение прочности на сжатие на 40% приводит к росту смещений в 3 раза, а увеличение составляющей массива на 33% - в 3 - 4 раза. Увеличение сопротивления крепи значительно снижает смещения при низких значениях остаточной прочности и сопротивления крепи и становится неэффективным при отношении этого сопротивления к составляющей массива р/д >0,006.
Полученные формулы позволяют определить смещения в выработке при наличии в ней начального свободного закрепного пространства и проанализировать его влияние.
Примем размер закрепного пространства на момент установки крепи равным А и выраженным в долях от радиуса Яп выработки в проходке. По мере проведения выработки и последующего развития горных работ давление д в окрестности выработки возрастает, растут и смещения пород. При этом до появления контакта между крепью и окружающими породами отпора со стороны крепи нет и на этой стадии р = 0. Данные натурных наблюдений [3] показывают, что в условиях развившихся деформаций породы на контуре выработки находятся в руинном состоянии. Примем на этой стадии развития смещений остаточную прочность при одноосном сжатии —о = 0.
В момент возникновения контакта между крепью и боковыми породами перемещение пород внутрь выработки ид = А и это соответствует определенному радиусу Яр' ЗНД. Принимая р = 0 и —о = 0, из формулы (6) получим:
Яр/Яп =41 + 2 А/Япе 0. (7)
При отсутствии отпора со стороны крепи и остаточной прочности —о = 0 напряжения —г и —д в зоне запредельного состояния согласно формулам (3) равны нулю, поэтому породы внутри этой зоны до появления контакта с крепью находятся в ненапряженном состоянии. Поскольку при этом предполагается, что остаточная прочность —о = 0, то тогда и коэффициент сцепления здесь к0 = 0. Это характерно для сыпучей среды и, рассматривая действие крепи на неё, можно воспользоваться уравнениями статики сыпучей среды, принимая ог активной, а од пассивной составляющими.
После ликвидации зазора А в работу вступает крепь и по кон-туру радиусом Я оказывает давление р , равное её сопротив-лению, на окружающие породы. В результате этого на границе предварительно сформировавшейся запредельной зоны возникает давление р', которое при постоянном отпоре р сохраняется с ростом смещений: р' = р/(Яр/Я)(с - 1)/с.
ид _
После перекрытия закрепного пространства деформирующимися породами с дальнейшим ростом давления д зона запредельного состояния увеличивается и её радиус в зависимости от давления р' из формулы (3) Яр = Яр' [2д/р '(к + 1]1/(к -1)
Перемещение ид' границы радиусом Яр' , возникающее при действии отпора р со стороны крепи, запишется из формулы (5):
и;=,: яр Яг <к -- (Я: ? -(к Г' ^г,г,
1(к -1)—r (Rp) + (k-1) p гк
(8)
где линейные размеры под интегралом выражаются в долях от Rp .
Это приводит к дополнительному перемещению пород внутрь выработки с момента их контакта с крепью, которое равно её смещению: ud = ud' Rp' /R..
Результаты расчета смещений крепи со свободным закрепным пространством согласуются с практическими данными: оставление свободного закрепного пространства размером 5-10 см приводит к росту деформаций крепи в 2-З раза, а увеличение пространства до 20-25 см - в 4-5 раз.
Полученные формулы можно применить для расчета смещений в выработке с забученным за-крепным пространством, которое при заполнении размельченной породой может рассматриваться как сыпучая среда с углом внутреннего трения p . Можно также оценить эффективность взрывной забутовки закрепного пространства.
Согласно этим расчетам эффект воздействия забутовки закрепного пространства является существенным: смещения при тщательно выполненном забучивании снижаются более чем в три раза. Качество забутовки значительно влияет на смещения. Обеспечение её плотности на 80% приводит к росту смещений в 1,5-2 раза в сравнении с плотной забутовкой.
Взрывное забучивание закрепного пространства приводит к снижению смещений примерно в 1,4 - 1,5 раза.
Результаты расчета смещений в окрестности выработок с учетом дилатансии пород по приведенным выше формулам соответствуют реальным величинам, однако не отражают некоторые особенности проявления горного давления. Так, при одинаковых свойствах окружающих пород смещения кровли, почвы и боков горизонтальной выработки, закрепленной кольцевой крепью с одинаковым за-крепным пространством, по расчетам одинаковы. В действительности смещение кровли превышают смещения боков и почвы. Кроме того, при отсутствии отпора со стороны крепи и остаточной прочности пород, равной нулю, смещения в выработке по расчету стремятся к бесконечности. Однако в реальных условиях смещение незакрепленной почвы,
состояние которой может быть близким к руинному, имеет конечную величину. Такое расхождение объясняется тем, что в расчетах не учитывается сила тяжести пород, находящихся в ЗНД.
На элементарный объем вертикальной оси симметрии выработки действуют главные напряжения — и —в и объемный вес у. Решая уравнение равновесия с учетом этих составляющих, аналогично (3) для кровли получим:
— = - —</(к -1) + [—0/(к - 1) + р]г(к - 1 + уЯ[г -
г(к - 1)]/(к - 2);
Я(- 1)[—0/(к - 1) + р - уЯ/(к - 2)] + ЯруЯ/(к - 2)
- [2ц(к - 1) + 2—0]/(к - 1) = 0. (9)
Определение Яр из последнего уравнения
осуществляется приближенным способом.
В качестве примера были произведены расчеты для следующих условий.
Выработка пройдена на глубине Н = 1000 м сечением в проходке £ = 13,85м2, закреплена
замкнутой крепью с сопротивлением р = 50 кПа =
0,05 МПа. Характеристики пород: —сж = 30 МПа; р= 30°; —0 = 0; £у° = 0,1; у = 25 кН/м3.
Давление пород в нетронутом массиве на глубине расположения выработки: д = уН = 25 МПа. Радиус равновеликой выработки кругового очертания: Я =4£7я = 2,1 м . Тогда уЯ = 0,0525 МПа. Параметры, имеющие размерность давления, в долях от составляющей д: р = 0,002; —сж = 1,2; уЯ = 0,0021. Характеристики состояния массива: к = 9; —г(Яр) = 0,2.
Согласно расчетам смещения контура выработки составляют: в кровле - 0,147, в почве - 0, 132, в боку - 0, 139 м. Сближение кровли и почвы равно 0,279, а боков - 0,278 м.
Таким образом, смещения кровли, почвы и боков выработки, закрепленной замкнутой крепью, установленной вплотную к окружающим породам, различаются незначительно.
Аналогичные расчеты были произведены для тех же условий при наличии свободного закреп-ного пространства А = 10 см. Смещение кровли увеличивается до 0,61 м, т.е. более чем в 4 раза, а боков - до 0,342 м, т.е. почти в 2,5 раза, и если в случае крепления без зазоров разница в смещениях кровли и боков несущественна, то при зазоре 10см смещение кровли почти в два раза превышает смещение боков.
Расчет напряжений и размера ЗНД почвы выработки с учетом силы тяжести пород производится по формулам (9) с противоположным знаком перед у. Эти формулы позволяют рассчитать смещение незакрепленной почвы. При этом оно имеет конечную величину и в том случае, когда остаточная прочность пород при одноосном сжатии равна нулю.
Согласно этим расчетам наиболее заметное влияние на подъем почвы в выработке оказы-
вает угол внутреннего трения. Уменьшение его от 30 до 20° сопряжено с увеличением смещения почвы более чем в два раза.
«Пучащие» породы зачастую требуют дополнительных мероприятий по уменьшению смещений почвы, например, анкерования. При расчете параметров анкерования задаются длиной анкеров. Поскольку концы анкеров необходимо располагать в упругой зоне, то радиус Яр ЗНД оказывается определенным и из формулы (9) находят сопротивление крепи р, соответствующее этому радиусу, после чего находят плотность установки анкеров.
СПИСОК ЛИТЕРАТУРЫ
1. Ставрогин А.Н., Протосеня А.Г. Механика деформирования и разрушения горных пород. - М.: Недра, 1992.
2. Волков В.М., Клыков А.Е., Широколобое Г.В. Влияние прочностных свойств горных пород на устойчивость выработок. - В сб. «Научно-технические проблемы подземной разработки месторождений».
- Кемерово: КузПИ, 1991.
3. Заславский Ю.З., Дружко Е.Б. Новые виды крепи горных выработок. - М.: Недра, 1989.
□ Авторы статьи:
Волков Клыков
Владимир Матвеевич Александр Ефимович
- канд. физ.-мат. наук, доц. каф. - канд. техн. наук, доц. каф. сопро-
прикладной математики тивления материалов
Расчеты показали, что применение анкеров для укрепления почвы может снизить пучение более чем в 5 раз. Плотность установки анкеров может быть уменьшена с увеличением их длины, однако при этом необходимо учитывать неизбежность возрастания смещений. При больших значениях допускаемого смещения почвы применение анкеров становится нерациональным вследствие требуемой большой их длины, которая на практике не превышает 2 м. В этом случае возможно применение замкнутой крепи, обладающей большой вертикальной податливостью.