УДК 622.271.332.624.13
Жабко Андрей Викторович
доцент кафедры маркшейдерского дела,
кандидат технических наук,
Уральский государственный
горный университет
620144, г. Екатеринбург, Куйбышева 30,
Тел.+79122192411,
е-таП: [email protected].
ПОВЕРХНОСТИ СКОЛЬЖЕНИЯ И ПРЕДЕЛЬНЫЕ ПАРАМЕТРЫ НЕОДНОРОДНЫХ И АНИЗОТРОПНЫХ ОТКОСОВ
Zhabko Andrey Victorovich
associate professor, mine surveying department, candidate of engineering sciences The Urals State Mining University, 620144, Yekaterinburg, 30, Kuibyshev st. Tel. +79122192411, е-mail: [email protected]
SURFACES OF SLIDING AND LIMIT PARAMETERS OF THE NON-UNIFORM AND ANISOTROPIC SLOPES
Аннотация:
В статье предлагается механический принцип для расчета устойчивости откосов. На его основе выстраивается теория расчета их устойчивости. Выводятся основные уравнения условия равновесия призмы смещения. Получены дифференциальные уравнения, определяющие геометрию потенциальной поверхности скольжения. Дается решение задачи на основе методов вариационного исчисления. Рассмотрены вопросы преломления поверхности скольжения в анизотропных горных массивах.
Ключевые слова: откос, условие равновесия, устойчивость, экстремум функционала, дифференциальное уравнение, поверхность скольжения, анизотропия, угол излома.
Abstract:
The article proposes a mechanical principle to calculate slopes' stability. The theory of calculation of their stability is based on it. The basic equations of conditions of a displacement prism equilibrium are derived. Differential equations that determine the geometry of a potential slip surface are received. The solution of the problem in terms of the methods of the variations calculus is given. The questions of slip surface refraction in anisotropic rock mass are discussed.
Key words: slope, equilibrium condition, stability, composite function extremum, differential equation, slip surface, anisotropy, torsion angle.
Рассмотрим участок горного массива (призму смещения), который может смещаться по некоторой произвольной поверхности скольжения. Разобьем призму смещения на отсеки и предположим, что функция поверхности скольжения известна. Тогда система (призма смещения), состоящая из п отсеков имеет 2п степеней свободы, то есть каждый отсек может независимо перемещаться по площадке скольжения и вращаться относительно некоторой точки. Однако заметим, что вращение отдельных отсеков не влияет на устойчивость призмы смещения в целом, а приводит лишь к ее деформированию1. Поэтому п степеней свободы, обусловленных вращением отсеков, можно исключить из рассмотрения. Если при смещении призмы отсеки взаимодействуют между собой, это свидетельствует о системе с удерживающими связями. Такой характер нарушения устойчивости будет иметь место тогда, когда нижерасположенные отсеки смогут перемещаться лишь под воздействием вышерасположенных. При выполнении последнего условия система будет иметь одну степень свободы, в качестве обобщенной координаты можно выбрать, например, горизонтальное смещение, одинаковое для всех отсеков.
Условие равновесия призмы смещения в общем виде определяется уравнением
[1 - 3]:
\ [ - у )( -1)- C(1 + У'2 )+ (' + ] + (( - Eo)- /(( - ^) = 0, (1) где 7 - объемный вес горных пород; у, у - функции линий откоса и поверхности скольжения, соответственно; у' - производная функции поверхности скольжения;
1 Здесь речь идет о незначительных деформациях, которые не приведут к состоянию системы, резко отличной от той, для которой составлялось условие равновесия.
24
аХ+(£,-Е>)-/(Т - Т>) = 0. (2)
1
ТО,Е0,Т1,Е1 - внешние касательные и нормальные реакции на вертикальных гранях призмы смещения, соответственно, слева и справа; / = tgф - коэффициент внутреннего трения (тангенс угла внутреннего трения); С - сцепление массива горных пород; Е', Т' - соответственно, производные функций нормальной и касательной составляющих межблоковой реакции по аргументу.
При углах наклона поверхности скольжения, превышающих угол внутреннего трения (( > ф), условие равновесия призмы смещения определяется уравнением
"т(Г - У) - /)-С(1 + у2)"
_ 1 + у'2 _
А для участков с углами наклона поверхности скольжения (( < ф) условие равновесия имеет вид:
- у К»'-. /)-С ( + у'2 ) + ((-Ео)-/ (Т - То ) = 0. (3)
1 + /у \
Поверхность скольжения в массиве может зарождаться в однородных массивах вследствие действия гравитационных сил, может быть предопределена геологическим строением массива (тектоническими нарушениями, трещиноватостью различных рангов, слоистостью, сланцеватостью) или являться следствием их комбинации. Кроме того, из уравнений (2, 3) видно, что призма смещения может находиться в равновесии как в свободном состоянии, так и под действием внешних сил.
Определим пределы применимости функционала (2). Дело в том, что для присутствия межблоковых связей необходимо, чтобы вышерасположенный (вдоль поверхности скольжения) отсек имел большую энергию, чем нижерасположенный. Для выполнения этого требования необходимо
> о,
ёх
где ¥ - подынтегральное выражение функционала (2).
Считая высоты смежных отсеков равными, получаем следующее неравенство:
= /(1 + 2/у - у'2 )> о, (4)
ах
где у' - вторая производная функции поверхности скольжения по х.
Учитывая, что функция поверхности скольжения вогнутая, определяем максимально возможное значение производной функции поверхности скольжения, определяемое функционалом (2). Таким образом, согласно неравенству (4), максимальный угол наклона зарождающейся поверхности скольжения (в предельном равновесии) не может
превышать величины < ^П + . Однако если поверхность скольжения предопределена геологическим строением массива, а угол поверхности ослабления превышает предельный, то в этом случае межблоковые связи (реакции) не совершают работу по перемещению. Это же относится к выпуклой поверхности скольжения. Поэтому при использовании функционала (1) условие равновесия будет иметь вид:
/ [у - У )(У - /)-С (1+у'2)]+(( - Ео)-/(Т - То ) = о. (5)
Таким образом, функционалами (2, 3, 5) описываются условия равновесия по произвольной поверхности скольжения.
Рассмотрим анизотропный откос. Под анизотропией в общем случае будем понимать систему поверхностей ослабления с характеристиками паспорта прочности С' и /' = tgф/, как правило, меньшими, чем показатели прочности массива С и ф . Рассмот-
рим вначале произвольный плоский откос, разрушение которого произойдет по поверхности скольжения, полностью совпадающей с плоской поверхностью ослабления (рис. 1). Очевидно, из постановки задачи имеем условие а > в . Если в < ф', то разрушение откоса произойти не может, поэтому будем исходить из неравенства в > ф'.
Поместим начало прямоугольной системы координат в точку выхода поверхности ослабления в откос, тогда условие равновесия будет иметь вид [1, 2]:
Hctgа
|
у(( - tgpx )(в- Г)
с '
Hctgа+а
йх + |
Нс%а
у(Я - tgвx )(в- /') 1 + tg2p
1 + tg2p
откуда предельная высота откоса выразится зависимостью [1, 2]:
^ tgа + Л/ tgаtgв = И tgа- tgв 1
Н = И
-С'
йх = 0.
(6)
Рис. 1 - Схема к расчету предельных параметров анизотропного откоса
Формула (6) рекомендуется Правилами [4] и получена из условий геометрической статики. Однако, как показано выше, она является следствием более общего уравнения. Заметим, что в случае, когда в>п /4 + ф' /2, в качестве условия равновесия необходимо использовать функционал (5), который также приводит к формуле (6).
Рассмотрим теперь более общий случай потери устойчивости анизотропного откоса и методику определения его предельных параметров (рис. 2). В этом случае поверхность скольжения может частично совпадать с поверхностями анизотропии и проходить вкрест ее, подобно поверхности в однородном откосе. В этом случае условие равновесия, а также угол излома 9 существенно зависят от координат точек излома, угла падения поверхностей ослабления, а также от физико-механических свойств пород вдоль плоскостей анизотропии и вкрест ее.
Ограничимся рассмотрением следующей расчетной схемы. Пусть в массиве имеются плоские поверхности ослабления, падающие в сторону выработанного пространства под углом в > ф' к оси х. В этом случае поверхность скольжения может частично совпадать с поверхностями ослабления.
Рис. 2 - Поверхность скольжения в анизотропном откосе
Необходимо определить условие в точке пересечения этих поверхностей с криволинейной частью поверхности скольжения, построенной для изотропной части массива на участке Ф > ф. Предположим, что найдется такой параметр п > 0, одинаковый для всех отсеков, и, если на него разделить величину сцепления (или умножить объемный вес), условие предельного равновесия будет выполняться в каждой точке поверхности скольжения. Поместим начало системы координат в точку пересечения поверхности скольжения с откосом и, используя функционал (2), получим следующую задачу вариационного исчисления:
) ( - гХУ - /)-41+у2)+1 ^ (7)
X 1 + у'2 X 1 + «ё2р V 7
X ^ хо
где g - постоянная; X = С / уп, X' = С' / уп - заведомо неизвестные параметры, определяющие предельную высоту откоса (очевидно, X '=(с ' / с )х).
Условие в точке излома поверхности скольжения имеет вид:
Е +(Ф' - /1X у, = Г2 +(Ф' - /2 )2 у2
2/2 '
(8)
где Е, Е2 - подынтегральные выражения слагаемых (7); Еу,, Е2у2 - частные производные подынтегральных выражений по производной функции поверхности скольжения; Ф' - производная функции, по которой перемещается точка разрыва (поверхность ослабления).
Условие (8) запишется так:
«ё¥- «ёФ «ёР- «ёф' + («ёв «ё„Х + 2«ёФ«ёУ- «ёУ = Х-Х'
^1ё^-^ёГ+(«ёв-«ёу) (1+«ё2^)2 =~,
(9)
где у - угол наклона поверхности скольжения к оси х (к горизонту) в точке преломления; к - глубина точки излома (расстояние по вертикали от поверхности откоса или бермы до точки излома).
Рассмотрим задачу о преломлении поверхности скольжения вследствие ее перехода в литологический слой с иными механическими характеристиками. Аналогом поставленной задачи является задача о преломлении луча света на границе сред с разными оптическими плотностями в постановке принципа Ферма (1660 г.). Примем у > ф - угол наклона поверхности скольжения к оси х до преломления; в - угол наклона контакта
литологических слоев к оси х; ю > ф - угол наклона поверхности скольжения к оси х после преломления (слой с параметром А'); С', ф' - механические характеристики лито-логического слоя, в который переходит поверхность скольжения. Составляя выражения, подобные (7) и (8), с использованием функционала (2), получим уравнение:
;§ф - (8ю-^+(гёР- (8¥)1+2;§ф1ё¥ -2ё ¥-(р- 18Ю)1+218Ф' 18°)-218 ю=, (¡о)
1 + 1ё ^ 1 + ю (1 + 1ё2у) ( + 1ё2ю
Анализируя условие (10), замечаем, что в случае С = С', ф = ф', преломления поверхности скольжения не происходит, то есть у = ю. Кроме того, очевидно, что поверхность скольжения не может после преломления вернуться в первоначальный литологи-ческий слой, то есть предельным значением угла ю является угол в. В этом случае скольжение произойдет по параллельной или совпадающей с контактом поверхности. Положив в уравнении (10) ю = в, получим условие (9).
Рассмотрим пример расчета неоднородного откоса. Требуется определить предельную высоту откоса Н , состоящего из двух горизонтально залегающих слоев ( в = 0° ), верхний слой имеет мощность, равную 1 / 3Н, при следующих условиях (нумерация слоев снизу): а = 40°;у = у 1 = 2,5т/м3;С1 = 2т/м2;С = 10т/м2;ф1 = 25°;ф = 30°.
Предположим, что в точке излома поверхности скольжения выполняются условия у > ф, ю > ф1, а абсцисса точки излома х0 находится под откосом. Тогда уравнения поверхности скольжения в нижнем слое имеют вид [3]:
Лк/ -1)'2 + 2(( + /)) +1 - к/ () (п а + .
кх-у1 =А^ —^—^^-= Ат]1 (р), - 1ё|---- I^У ^ 1ёф
* /у{2 - 2/2у1 + к - / + к/2 \ 4 2 ) *
о 1 + к/ (1 + У22 )2 о ( ) ' +
кх - У2 = А и ■ '3 ( - Л п ,-7 = П (Р) У2 > Ш
1 + /2 2У23 -(( + 3/)) + 2к/у2 + к - / где п(р)- некоторая функция параметра р = у'.
Произвольная постоянная для первого уравнения получена из условия трансверсальности на пересечении поверхности скольжения с откосом, а для второго - из условия равенства производных в точке у[ = у2 = 1ёф.
Уравнение поверхности скольжения под откосом после преломления имеет вид
[3]:
кх -У3 = + ^) '3 (, + 3 (1)+'У3+ ),, '+ . , = ( + 4 >13(Р).
2у3 -(к + 3/ )У3 + 2к/у3 + к - /
Для определения углов поверхности скольжения по обе стороны от точки излома,
абсциссы и глубины точки излома к, а также параметра + составляем условия:
2 2 -С
- 1ёф 1ёю - 1ёф^ + (в - 1ёу)1 + 21ёф1ё^- 1ё ¥ - (( - 1ёю) 1 + 21ёф11ёю - 1ё ю = I С
1 + 1ё2¥ 1 + 1ё2ю (1 + 1ё>) (1 + 1ё2ю)2 к
_(1 + 1ё2ю)2_
к ( + ^) 21ё3ю -(( + 3/ )2ю+ 2/ёю + к - /
к = А 1±/_(1 + 1ё2¥)2_
1 + /2 21ё3у - (к + 3/)1ё> + 2к/1ёУ + к - / к + 2 Н/3
х0 =
к 28
«ёф
а+ф п 2 4
К (р К+ТК (р)
к - р
ёр + |
«ёф
к - р
ёр = х0
где п'(р)- производная функции по параметру р .
Уравнение поверхности скольжения под бермой имеет вид [2]:
С (1+/42)
Н - /4
У 2/43 - Э/У2 - /
Производную поверхности скольжения на границе участков откоса и бермы Ь определяем из условия отсутствия ее излома:
(( + ))
(1+ь2 )2
= С (1+Ь2)
2Ь3-(к + 3/ )Ь2 + 2к/Ь + к - /1 у 2Ь3 - 3/Ь2 - /1 ' И, наконец, система уравнений, требующая выполнения условия предельного равновесия и осуществляющая связь X и Н :
1+/р 1+/р
а+ф п |
2 4 I
К(р)
к - р
«ёУ ёр + |
«ёф
Мь (р)-с'
1+р
(р )
к - р
ёр +
+
т(р - /) 1 + р2
( + ё )П3 (р )- С
«ёЮ
? К(р)
«ё| п+ф
(X + ё )л3 () +С2 1 V Т (р - /)
а+ф
2 4 I
л к - р
п 1 *
ёр+«7 ММ ёр + Ь к - р
«ёф
«ёЮ
к -ёр +- Ь
к - р у Ь
( + )3 (р) р=Н к - р к
1+р
П (р)-1
24
(р )
ёр = 0
Данных уравнений достаточно для определения всех неизвестных.
Рассмотрим пример расчета анизотропного откоса. Требуется определить предельную высоту откоса Н , имеющего поверхность ослабления, выходящую в его основание под углом в = 20°. Физико-механические характеристики пород и поверхности
ослабления: у = 2,5 т/м3; С' = 5 т/м2; С = 10 т/м2; ф' = 15°; ф = 30°, угол откоса а = 40°.
Предположим, что угол наклона криволинейной части поверхности скольжения в точке преломления поверхности скольжения у > ф и точка преломления находятся под откосом, тогда ее форма описывается уравнениями [2, 3]:
кх - /2 = (X + ё)
(1 + /22 )2
2 / 23 -(к + 3/))22 + 2/2 + к - /
= (X + ё)2 (р), /2 > «ёф, откос
Н - /3 =
С
(1+/32 )2
г 2/33 - 3//32 - / г
СП3 (р),
«ёф < /3 ^ «ё1 п + ф |, берма
Тангенс угла наклона поверхности скольжения в точке стыка откоса и бермы Ь определится из уравнения:
«ё
«ё
«ё
(1+ъ2 )2
= С (1 + ъ2 )2
( + й )2Ъ3-(к + 3 / )ъ2 + 2к/ъ + к - / у 2Ъ3 - 3 /ъ2 - / ' Условие в точке стыка имеет вид (9):
1ё¥- 1ёф 1ёР- 1ёф' + (.ёв + 21ёф1ё¥- 1ё2¥ = А-А'
"Г^ 1ё¥) (1+1ё2^)2 =~к~
Высота точки излома выражается уравнениями:
(1 + 1ё2¥)2_
к = (( + й)
21ё> - (к + 3/))ё2у + + к - /
Абсцисса точки излома:
Н
н - Г к •>
к = (к - 1ёв)х (А + й ) (р)
1ёУ
к - р
йр = х0
где п'2 (р)- производная функции по параметру р .
Система уравнений для определения А + й и Н :
х0 у(( - 1ёвх)ёв - 1ёф') - С'(1 + 1ё2в)
(1 + 1ё2 в)
йх +
1ёё¥
р - / 1 + р2
( + й ))п2 (()- С
(( + й ) '2 (р)
к - р
> +
+ -
Г п ф
с 2 чп+т
р - / 1 + р
2 П Ы- 1
П (р )
йр = 0,
I п ф
I 4 2
+ г рk±йM.íр) рйр + £л; (р »+* ]=Н.
•>к - р -1 V V 142)
1ё¥
У
У
42
Данных уравнений достаточно для определения всех неизвестных.
На рис. 3 показаны потенциальные поверхности скольжения и предельные параметры неоднородного и анизотропного откосов для рассмотренных выше примеров.
В заключение отметим, что уравнениями (9, 10) не исчерпываются все условия в точках излома, которые могут иметь место при решении задач для неоднородных и анизотропных откосов. Подобные условия находятся аналогично с применением уравнения (8) и использованием соответствующих функционалов (2, 3, 5), описывающих условия равновесия для конкретных углов наклона поверхностей скольжения и ослабления.
а = 40°
а = 40°
а
б
Рис. 3 - Примеры расчета неоднородного (а) и анизотропного (б) откосов
У
Литература
1. Жабко А.В. Расчет устойчивости откосов / А.В. Жабко // Маркшейдерия и недропользование. - 2011. - № 5. - С. 37 - 39.
2. Жабко А.В. Теория расчета устойчивости оснований и откосов / А.В. Жабко // Известия Уральского государственного горного университета. - 2011. - № 25 - 26. - С. 59 - 65.
3. Жабко А.В. Предельные параметры плоских однородных откосов / А.В .Жабко // Изв. вузов. Горный журнал. - 2012. - № 6. - С. 22 - 25.
4. Правила обеспечения устойчивости откосов на угольных разрезах. - Спб., 1998. - 208 с. (Минтопэнерго РФ. РАН. Гос. НИИ горн. геомех. и маркшейд. дела - Межотраслевой науч. центр ВНИМИ).