CC BY
55
© П.С. Шпаков, В.Н. Долгоносов, Ю.Л. Юнаков, М.В. Шпакова, 2015
УДК 622.1: 622.271
П.С. Шпаков, В.Н. Долгоносов, Ю.Л. Юнаков, М.В. Шпакова
АНАЛИТИЧЕСКИЙ СПОСОБ РАСЧЕТА УСТОЙЧИВОСТИ ОТКОСА НА СЛАБОМ ОСНОВАНИИ НЕОГРАНИЧЕННОЙ МОЩНОСТИ
При разработке месторождений полезных ископаемых открытым способом и отсыпке вскрышных пород во внешние и внутренние отвалы часто возникает задача определения устойчивости откосов, расположенных на слабом основании большой мощности. Существующие способы решения данной задачи зачастую графические, обладают рядом недостатков, которые устраняются в предложенном аналитическом способе расчета. Результаты наблюдений за деформациями карьерных откосов на ряде месторождений Казахстана, выполненных авторами, позволили создать модель деформирования откосов на слабом основании и выявить форму и местоположение потенциальной поверхности скольжения. Показано, что потенциальная поверхность скольжения в откосах, расположенных на слабом основании большой мощности, может аппроксимироваться дугами трех окружностей разных радиусов. На основе положений теории предельного равновесия разработан аналитический способ расчета параметров предельного откоса. Реализация решения задачи на ЭВМ значительно облегчает расчеты и анализ по обоснованию параметров устойчивых откосов, что дает способу заметные преимущества по сравнению с существующими графическими решениями. Для оценки точности разработанного способа расчета были выполнены поверочные расчеты параметров предельного откоса на участках оползней бортов Тургайского бокситового рудоуправления предлагаемым аналитическим способом и численно аналитическим способом профессора П.С. Шпакова. Сравнительные результаты решения аналитическим и численно-аналитическим способами показали, что результаты расчетов не сильно отличаются друг от друга (в пределах 3%).
Ключевые слова: отвал, коэффициент запаса, основание отвала, мощность отвала, расчетная схема, предельная высота, генеральный угол откоса, предельный откос.
При разработке месторождений полезных ископаемых открытым способом часто возникает задача определения параметров откосов, расположенных на слабом основании большой мощности. Так, например, при складировании вскрышных пород на территориях пересыхающих водоемов (отвалы Качарского карьера) либо при залегании в основании отвалов мощных линз слабых пластичных глин, а также в случаях, когда в нижней части стратиграфического разреза борта карьера присутствуют мощные горизонты слабых песча-но-глинистых пород, которые отличаются более низкими прочностными характеристиками по сравнению с породами вышележащих горизонтов (Качарский и Тургайские карьеры).
Характерной особенностью механизма деформирования прибортового массива является наличие выпора - выдавливания пород слабого основания у подошвы откоса (рис. 1).
В случае залегания в основании откоса слабого слоя большой мощности поверхность скольжения представляет собой сложную линию переменной кривизны, состоящую из трех участков - дуг окружностей (рис. 1).
Рассмотрим общую схему построения поверхности скольжения.
1) Первая дуга (участок СМ) начинается с точки С на глубине Н90 от верхней площадки уступа под углом ц = (45° - р / 2) к вертикали или е = (45° + р / 2) к горизонтали. В точке М дуга пересекает верхнюю границу слабого слоя под углом входа 9. Формулы для вычисления параметров Н90 и 9 приведены в работе [1]. Координаты точки входа М определяем как результат пересечения прямых ОМ и МС
(1)
(2)
х = Ус - х^д[45° + (8 + 6 - ц) / 2] ; ум = хм ■ tg8 хм ~
tg8-tg[45° + (8 + 6-ц)/2]
где 8 - угол наклона слабого основания.
Определим координаты центра первой окружности - точки 01. Уравнения радиусов в точках М и С
у = ум -(х - хм) • <^д(8 + 6) ; у = Ус -(х - Хс) • <^де .
Координаты центра окружности 01 определяем как результат пересечения данных радиусов (прямых С01 и МО 1)
х = Ум - Ус - Хс^де + х„<^д(8 + 6) ; 0 ^д(8 + 6) -е*дв ; (
У о = Ум - (х0 - хм) • с*д(8 + 6) или у0 = ус - (х0 - хс) • Сде . Центральный угол первой дуги: Р2 = 90° - (8 + 6 + ц) .
2) На границе со слабым слоем происходит преломление поверхности скольжения под углом
Д6 = (р-р')/2. (4)
Рис. 1. Схема откоса на слабом основании большой мощности
В слабом слое на участках MN и ND может быть построено бесконечное множество поверхностей.
На втором участке MN поверхность скольжения аппроксимируется семейством дуг окружностей радиусами R2. (от R2min до R2max) с центрами в точках O2i. В точке наибольшего заглубления поверхности скольжения N. (в которой касательная направлена горизонтально, а нормаль - вертикально), вторая дуга переходит в третье семейство NA радиусами R3. (от R3min до R3max) и центрами в точках O3i. Ввиду преломления поверхности скольжения в точке М, угол между касательной к дуге MN в этой точке и контактом
0! = 0 - Д0 . (5)
Поверхность скольжения на третьем участке имеет обратный уклон и выходит в основание откоса под углом
ц = 45°-р'/2 . (6)
Минимальная граница семейства поверхностей в слабом слое соответствует выходу поверхности в основание откоса (с минимальным заглублением hmin) и нулевой ширине выпора, а максимальная граница соответствует выходу поверхности в точке D, максимальным ширине L и глубине h выпора, когда
1 1' 1 max ^ max 1 '
точка сопряжения второго и третьего участков поверхности скольжения находится под нижней бровкой откоса (s = 0).
3) Определим координаты центров второго семейства окружностей O2i. Уравнение радиуса второй окружности в точке М:
у = yM-(х - xM) •cig(8 + 01). (7)
Второй радиус представляет собой вертикаль х = s = const, х е [0; а]. При x = 0 поверхность скольжения выходит на расстоянии Lmax от основания откоса, а при x = a - в основание откоса. На прямых, которым принадлежат вторые радиусы (O2.N.), также будут находиться центры дуг окружностей третьего участка O3..
Построение начинаем с внешней поверхности, которая выходит на расстоянии от основания откоса. Подставляя в уравнение (7) значение x = 0, получим координаты точки 02max [х = 0; у = yM + xM • ctg(8 + 01)] .
Определим максимальную глубину поверхности скольжения под основанием откоса
hmax = R2max - Уо2т« = R2max - [Ум + XM ' ctg(8 + 01)] . (8)
Центральный угол дуг второго семейства: Р2 = (8 + 01) .
4) Вычисляем параметры третьего семейства дуг. Центр дуги максимального радиуса находится на прямой O2maxN.
Центральный угол дуг третьего семейства: Р3 = ц1. Ордината максимального радиуса третьего семейства
cos Ц
У0 3 = hmax^-. (9)
0 3 max 1 - cos ц1
Максимальная ширина выпора у основания
Lmax = Уо зтх ^ tgЦ1 - S . (10)
5) Определим общую методику построения промежуточных кривых скольжения на втором и третьем участках. Построения начинаем со второго участка. Задаемся величиной смещения вертикального радиуса в сторону откоса x = a и подставляем в уравнение (7). В результате получим промежуточные значения координат точек
°2, [х, = а; У, = Ум - (а - хм)' + ©а)] . (11)
Максимальное заглубление промежуточной поверхности скольжения
Кпах , = И2 , - У0 2, = И2, - [Ум - (а - ХМ) ' С^(8 + 61)] . (12)
Центры промежуточных дуг третьего семейства находится на вертикальных прямых O2N. Ордината центра дуги окружности
. еоэ ц
У = п,-1—
0 3■ ' 1 - еоэ ц1 . (13)
Ширина выпора у основания
Ц = уо 3, • . (14)
6) Методика интегрирования сдвигающих, нормальных и удерживающих сил, действующих по линиям скольжения АMND подробным образом описана в работе [2], отдельные исследования по этой проблемы детально изучались в работах [3-11].
Определим высоту и вес элементарного столба пород, расположенного выше поверхности скольжения, для чего призму разделяем на три части, соответствующие трем дугам поверхности скольжения.
Первый участок, соответствующий дуге CM, делим на два участка относительно проекции верхней бровки откоса на поверхности скольжения (т. Е). Участок под откосом (МЕ):
П(х) = х-^а-Уо . (15)
Участок под горизонтальной площадкой (ЕС):
П( х) = Н - Уо ^^И^-^Х-ХО)2 . (16)
На участке объемный вес пород постоянный и может быть вынесен за пределы интеграла, поэтому возможно использовать полученные выше решения.
Определим сдвигающие силы, действующие на призму возможного обрушения на первом участке, путем интегрирования элементарных сдвигающих (касательных) сил по поверхности скольжения
Т1 сдв = 1 4Тсдв = 1 4Тсдв + 1 4Тсдв (17)
мс МЕ ЕС
Результат интегрирования аналогичен решениям, полученным в работе [2], формулы - (2.37) и (2.38).
Определим суммарные нормальные силы, действующие по первому участку поверхности скольжения, путем интегрирования элементарных нормальных сил по дуге
N1 =1 =1 +1 dN
МС МЕ ЕС . (18)
Результат интегрирования соответствует решениям (2.42) и (2.43) [2].
Второй участок соответствует дуге MN. Общая высота породного столба на этом участке может быть определена по формуле
Ь(х) = х•да-У0 + у/Я2 -(х- х0)2 . (19)
Так как откос и его основание сложены различными породами с разным объемным весом, то вес элементарного столба
¿Р( х) = [ ух • tgа + Уг(-У0 + х0)^)№ , (20)
где у и у1 - объемный вес пород откоса и основания соответственно. Сдвигающие силы, действующие на втором участке
т г ,-р
2 сдв ~ 1 сдв
им . (21)
Интегрирование выполнено на основе решения (2.37) [2], которое с учетом условия (20) примет вид
1 ¿Тсдв = 1 1 (ух •да+ 7^-У0 + (х - х0 )dх =
им сдв Я ме ^ '
1 хм
= Я 1 [ух%а- х (YХоtgа + Y1 У0) + 71х0У0 -
хи
хм ---_
- 2Я ¡ТЯмх-хо7d [ я2 - (х - х0 г
2Я
хи
Вычисляем два полученных интеграла а)
.. хм ,-
К Я2 - (х - х
¡уЯ-х-х^б[Я2 -(х-х0)2] = -[Я2 -(х-х0)2]3/2|
2Я * 0 [ 0 ] 3Я^
хи
= 3Я([Я2 -(хи -х0)2]3/2-[Я2 -(хм - х0)2]3/2)
б)
1 хм
я | [У x2tgа - х(7 Хоtga + 71У0) + У1 х0У0 =
Я хИ
Х3 Х2
у tga —---^ (У Хоtga + 71У0) + У1 х • Х0У0
(22)
1 Х3 Х2 = Я [у tga^м —м (1 Хоtga + У1У0) + У1 ХмХ0У0 ] -
1 Х3 Х2 - Я [У tga^Г -"2^ (У Хоtga + 71 Уо) + 71 х^Уо ]
Нормальные силы, силы, действующие на втором участке И2 =| ¿И
им . (23)
Интегрирование выполнено на основе решения (2.42) [2], которое с учетом условия (20) примет вид
1 г ' ■ ' '21 Ш2 - (х - ^ )2 ¿х =
| ¿П = ш | (у х • tga-YlУ0 + у^ Я2 - (х - х0 )2 ^
NM ^М
Гя2 Г . 71 • уо
М I - . М I
хм
Ш ¡(2ххо - х2 - Хо2 + Я2 ) ¿х.
(24)
Результат интегрирования соответствует решению (2.42) [2] с учетом объемного веса пород основания уг
Третий участок (зона выпора основания) соответствует дуге ОД. В общем случае третий участок необходимо делить на два интервала точкой М, которая находится под нижней бровкой.
Тогда, при х > 0 высота породного столба определяется по формуле (19), а интегрирование сдвигающих и нормальных сил осуществляется аналогично второму участку по формулам (22) и (24) от нижнего предела х = 0 до верхнего х = а. При х < 0 высота породного столба определяется выражением (16) при Н = 0. Интегрирование сдвигающих и нормальных сил по дуге осуществляется по решениям (2.38) и (2.43) [2] для пород основания с объемным весом уг Разработан алгоритм и программа для решения данной задачи. Определим параметры предельного откоса на слабом основании большой мощности. Угол откоса а = 28°, прочностные характеристики пород массива к = 10 тн/м2, р = 18° и слабого слоя к' = 3,5 тн/м2, р' = 12°. Объемный вес пород массива у = 2,0 тн/м3 и основания у' = 1,8 тн/м3.
N
Рис. 2. Схема к аналитическому решению
Результаты расчета предельных параметров откоса
Метод расчета Параметры предельных откосов, м
Н пр r L h R! R2 R3
Аналитический 53,5 17,6 10,1 12,4 147,5 106,2 55,8
Численно-аналитический 56,4 17,6 9,0 12,6 157,8 107,9 56,5
Схема откоса с радиусами предполагаемой поверхности скольжения представлена на рис. 2.
Сравним результаты решения аналитическим и численно-аналитическим способами (таблица).
Приведенные в таблице результаты свидетельствуют о достаточно высокой сходимости (в пределах 3%) результатов расчета двумя способами.
_ СПИСОК ЛИТЕРАТУРЫ
1. Попов И.И., Шпаков П.С., Поклад Г.Г. Устойчивость породных отвалов. - Алма-Ата: Наука, 1987. - 224 с.
2. Долгоносов В.Н., Шпаков П.С., Низаметдинов Ф.К. и др. Аналитические способы расчета устойчивости карьерных откосов. - Караганда: Издательство ТОО «Санат-полиграфия», 2009. - 332 с.
3. Фисенко Г.Л. Устойчивость бортов карьеров и отвалов. - М.: Недра, 1965. - 378 с.
4. Методические указания по определению углов наклона бортов, откосов уступов и отвалов строящихся и эксплуатируемых карьеров. - Л.: ВНИМИ, 1972. - 163с,
5. Правила обеспечения устойчивости откосов на угольных разрезах. - СПб., 1998. - 208 с.
6. Шпаков П.С., Долгоносов В.Н., Шпакова А.П. Расчет параметров предельного откоса на слабом наклонном основании численно - аналитическим способом // Горный информационно-аналитический бюллетень. - 2003. - № 9. - С. 40-41.
7. Шпаков П.С., Долгоносов В.Н., Старастина О.В., Шпакова М.В. Сравнительный анализ способов расчета предельных параметров откоса на слабом наклонном основании // Горный информационно-аналитический бюллетень. - 2010 - № 4. - С. 302-306.
8. Попов В.Н., Шпаков П.С., Юнаков Ю.Л. Управление устойчивостью карьерных откосов. Учебник для вузов. - М.: Изд-во МГГУ, изд-во «Горная книга», 2008. - 683 с.: ил.
9. Шпаков П.С., Юнаков Ю.Л. Расчет устойчивости карьерных откосов. Учебное пособие. - Красноярск, 2006. - 202 с.: ил
10. Шпаков П.С., Поклад Г.Г. Расчет параметров предельного откоса на слабом основании большой мощности // Известия вузов. Горный журнал. - 1987. - № 4. - С. 35-39. 1ГЛЧ
КОРОТКО ОБ АВТОРАХ_
Шпаков Петр Сергеевич - доктор технических наук, профессор, e-mail: SPSP01@rambler.ru, e-mail: SPSP01@mail.ru,
Муромский институт (филиал) Владимирского государственного университета им. А.Г. и Н.Г. Столетовых,
Долгоносов Виктор Николаевич - доктор технических наук, профессор, e-mail: vdolgonosov@hotmail.ru,
Карагандинский государственный технический университет,
Юнаков Юрий Леонидович1 - кандидат технических наук, профессор,
зав. кафедрой, e-mail: yunakov11@rambler.ru,
Шпакова Мария Владимировна1 - кандидат технических наук,
ассистент кафедры, e-mail: SPSP01@rambler.ru,
1 Сибирский Федеральный университет,
Институт горного дела, геологии и геотехнологий.
UDC 622.1: 622.271
ANALYTICAL METHOD OF CALCULATION OF STABILITY SLOPES AT WEAK BASES UNLIMITED CAPACITY
Shpakov P.S., Doctor of Technical Sciences, Professor, e-mail: SPSP01@rambler.ru,
e-mail: SPSP01@mail.ru, Murom Institute (Division),
Vladimir State University named after Stoletovs, Murom, Russia,
Dolgonosov V.N., Doctor of Technical Sciences, Professor, e-mail: vdolgonosov@hotmail.ru, Karaganda State Technical University, Karaganda, Kazakhstan, Yunakov Yu.L.1, Candidate of Technical Sciences, Professor, Head of Chair, e-mail: yunakov11@rambler.ru,
Shpakova M.V.1, Candidate of Technical Sciences, Assistant of Chair, e-mail: SPSP01@rambler.ru, 1 Siberian Federal University, Institute of Mining, Geology and Geotechnology, Russia.
In the development of mineral deposits by open pit and dumping overburden in the external and internal piles often arises the problem of determining the stability of slopes, located on low basis high power. The existing methods for solving this task, often graphic, have a number of disadvantages which are eliminated in the proposed analytical method of calculation.
The results of observations of the pit slope deformations on the number of fields in Kazakhstan, carried out by the authors, has allowed to create a model of deformation of slopes on a weak base and identify the shape and location of the potential surface of sliding.
It is shown that the potential slip surface in slopes located on a weak base big capacity, can approximiately three arcs of circles of different radii. On the basis of the provisions of the limit equilibrium theory to develop an analytical method of calculating the parameters of the marginal slope. The implementation of the decision of tasks on the computer greatly facilitates the calculations and analysis of the justification of the parameters of stable slopes, giving way noticeable advantages over existing graphics solutions.
To assess the accuracy of the developed method of calculation were performed confirmatory calculations of the parameters of the limiting slope at the sites of landslides sides Turgay bauxite ore proposed analytical method and numerical analytical method of Professor P.S. Shpakov.
Comparative results of analytical solution and numerical-analytical methods showed that the results of the calculations are not much different from each other (within 3%)
Key words: blade, the safety factor, the base of the blade, the blade power, the design scheme, the maximum height, the general slope angle, slope limit.
REFERENCES
1. Popov I.I., Shpakov P.S., Poklad G.G. Ustoychivost' porodnykh otvalov (Stability of rocktion dumps), Alma-Ata, Nauka, 1987, 224 p.
2. Dolgonosov V.N., Shpakov P.S., Nizametdinov F.K. Analiticheskie sposoby rascheta ustoychivosti kar'ernykh otkosov (Analytical methods for calculating the stability of pit slopes), Karaganda, Izdatel'stvo TOO «Sanat-poligrafiya», 2009, 332 p.
3. Fisenko G.L. Ustoychivost' bortov karerov i otvalov (Stability pit and waste dumps), Moscow, Nedra, 1965, 378 p.
4. Metodicheskie ukazaniya po opredeleniyu uglov naklona bortov, otkosov ustupov i otvalov stroyash-chikhsya i ekspluatiruemykh karerov (Methodological guidelines for determining the angular Inclination on board otkosov ledges and piles constructed and operated karerov), Leningrad, VNIMI, 1972. 163 p.
5. Pravila obespecheniya ustoychivosti otkosov na ugol'nykh razrezakh (Rules to ensure the stability of slopes in the coal insert Zech), Saint-Petersburg, 1998, 208 p.
6. Shpakov P.S., Dolgonosov V.N., Shpakova A.P. Gornyy informatsionno-analiticheskiy byulleten'. 2003, no 9, pp. 40-41.
7. Shpakov P.S., Dolgonosov V.N., Starastina O.V., Shpakova M.V. Gornyy informatsionno-analitiches-kiy byulleten'. 2010, no 4, pp. 302-306.
8. Popov V.N., Shpakov P.S., Yunakov Yu.L. Upravlenie ustoychivost'yu kar'ernykh otkosov. Uchebnik dlya vuzov (Management of stable careereast slopes. Textbook for high schools), Moscow, Izd-vo MGGU, izd-vo «Gornaya kniga», 2008, 683 p.
9. Shpakov P.S., Yunakov Yu.L. Raschet ustoychivosti kar'ernykh otkosov. Uchebnoe posobie (Calculation of the stability of quarry slopes. Textbook), Krasnoyarsk, 2006, 202 p.
10. Shpakov P.S., Poklad G.G. Izvestiya vuzov. Gornyy zhurnal. 1987, no 4, pp. 35-39.
