DETERMINATION OF THE VALUE OF THE STABILITY MARGIN COEFFICIENT BASED ON THE ANALYSIS OF THE STRESSED STATE OF THE SOIL MASS IN THE FRAMEWORK OF THE LINEAR THEORY OF ELASTICITY
- Authors: Bogomolova O.A1, Zhidelev A.V2
- Affiliations:
- Volgograd State Technical University
- Branch «Vzlet» of Moscow Aviation Institute (National Research University)
- Issue: Vol 13, No 1 (2022)
- Pages: 58-72
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/CG/article/view/2001
- DOI: https://doi.org/10.15593/2224-9826/2022.1.05
- Cite item
Abstract
The introduction lists the main disadvantages of the traditionally used methods for calculating the stability of slopes and slopes and presents a calculation method in which, for the first time, in our opinion, an attempt was made to construct the most dangerous "slip line" based on the analysis of the stressed state of the soil massif. Further, in the first chapter of the article, three possible approaches to solving the problem of determining the angle of inclination of the most probable shear site at the studied point of the soil massif, based on the analysis of its stress state, are presented and analyzed. The second chapter presents the results of calculating the value of the coefficient of stability of the soil slope by L.G. Fisenko (one of the traditional calculation methods) for two types of sliding line: circular cylindrical (classical case) and constructed by the method of N.S. Nikitin taking into account the operating stresses; and by V.K. Tsvetkov - A.N. Bogomolov. All calculations are performed for one object, the geometric parameters of which and the numerical values of the physico-mechanical properties of the composing soil are given in the article cited below by S.N. Nikitin. As a result, it turned out that the first two coefficients of the margin of stability are respectively equal to K 1 = 1.6 and K 2 = 1.75, and the value of the third coefficient of the margin is equal to K 3 = 2.18. Comparing this value with the values of the coefficients of the margin of stability calculated by the method of S.N. Nikitin and the method of L.G. Fisenko, we see that it is 26 and 20 % higher, respectively, while K 2 > K 1 is 9.4 %. Based on the analysis of the calculation results for the example considered in this paper, it can be stated that methods based on the analysis of the stressed state of the soil mass give higher values of the stability margin coefficients than traditionally used. Thus, the value of the slope stability margin coefficient calculated by the method of V.K. Tsvetkov - A.N. Bogomolov turned out to be equal to K 3 = 2.18. Comparing this value with the values of the stability margin coefficients calculated by the method of S.N. Nikitin and the method of L.G. Fisenko, which are respectively equal to K 1 = 1.75 and K 2 = 1.6, we see that it is 20 and 26 % more than them, respectively. At the same time, K 1 > K 2 by 9.4 % - this indicates that taking into account the stresses during the construction of the sliding line makes its shape different from the circular cylindrical one, which affects the value of the stability margin coefficient.
Full Text
Введение Расчет величины коэффициента запаса устойчивости грунтовых массивов (откосов, склонов, оснований фундаментов) является одной из важнейших задач геомеханики. Как в свое время отметили Н.А. Цытович и З.Г. Тер-Мартиросян [1], решение этой задачи может быть получено известными методами механики грунтов, которые подразделены ими на три группы. К первой группе отнесены методы, основанные на предположении о том, что положение и форма поверхности разрушения известны [2; 3]. Методы второй группы позволяют определить величину критических нагрузок и геометрических параметров грунтовых массивов методами теории предельного равновесия [4; 5]. И, наконец, методы третьей группы опираются на анализ напряженного состояния грунтового массива, базирующегося на решениях задач линейной теории упругости и упругопластических задач [6; 7]. Рис. 1. Расчетные схемы методов первой группы: а - метода А.Л. Можеветинова; метода Г.М. Шахунянца б - при круглоциллиндрической и в - полилинейной форме следа поверхности разрушения (скольжения) Fig. 1. Calculation diagrams of the first group methods: а - A.L. Mozhevetinov's method; G.M. Shahunyants' method б - for circular-cylindrical and в - polylinear fracture surface trace (sliding) forms Известно, что результаты расчета устойчивости грунтового массива разными методами порой сильно отличаются. В этой связи представляется справедливым объяснение, которое дает этому факту проф. И.В. Федоров: «Характер происходящих деформаций грунта при обрушении земляных масс существенно зависит от свойств грунта. Невозможно создать такую расчетную модель, которая одновременно бы отражала все многообразие физико-механических свойств. Любая расчетная модель - это лишь идеализация реальных свойств того или иного материала. Расчетная модель должна выбираться таким образом, чтобы в основных чертах удовлетворять запросы практики» [8]. а б Рис. 2. Расчетные схемы методов второй группы: определение максимального значения интенсивности внешней нагрузки: а - при заданном очертании откоса, б - при заданной интенсивности внешнего воздействия определить контур равноустойчивого откоса Fig. 2. The calculation diagrams of the second group methods: determining the maximum value of the external load intensity: а - for a given slope outline, б - for a given external load intensity determine the contour of an equally stable slope Подавляющее большинство расчетных методов первой и второй группы (см. рис. 1 и 2) не в полной мере и весьма приближенно отвечают физической сущности задачи об устойчивости грунтового массива. По этому поводу отметим некоторые наиболее существенные их недостатки. 1. Форма и положение поверхности разрушения считаются заранее известными. Не устанавливается их непосредственная связь с физико-механическими свойствами грунтов, геометрией и геологическим строением массива грунта, поверхностными нагрузками. 2. Призмы разрушения (обрушения или выпора) грунта расчленяются на отдельные блоки, силы взаимодействия и трение между которыми не всегда учитываются. При рассмотрении равновесия блоков составляются одно или два уравнения статики (при плоской задаче их должно быть три), что часто приводит к получению результата K < 1, противоречащему условию предельного равновесия. 3. При определении суммарных удерживающих и сдвигающих сил операция сложения векторных величин выполняется как для скалярных (метод алгебраического сложения сил). 4. При вычислении величины коэффициента запаса устойчивости практически всегда учитывается лишь одна вертикальная, и то весьма приблизительно, составляющая напряжений (вес блоков). 5. При расчете устойчивости откоса игнорируется наличие его подошвы (рис. 3, а), т.е. не учитывается тот факт, что в месте его перехода в подошву возникают зоны концентрации напряжений, могущие повлечь за собой образование областей предельного состоянии грунта (областей пластических деформаций). 6. В расчетах устойчивости не учитывается величина коэффициента бокового давления грунта xо. Исключить эти и другие недостатки при построении следа наиболее вероятной поверхности разрушения СНВПР (в современной научной терминологии закрепилось название «линия скольжения», которая есть не что иное, как след поверхности разрушения или поверхности скольжения на вертикальном сечении откоса или основания сооружения) позволяют методы, основанные на анализе напряженного состояния (НС) грунтового массива. По всей видимости, впервые (работ иного автора, посвященных данному вопросу, мы не нашли в научной литературе) метод построения СНВПР по напряжениям был предложен С.Н. Никитиным [9] (рис. 3, а-г). Согласно версии С.Н. Никитина, «зарождение поверхности скольжения в однородных породах связано с развитием площадок скольжения, возникающих под углом 45º- к направлению наибольших главных напряжений s1 (где j - угол внутреннего трения пород) (рис. 3, б). Направления s1 в массиве откоса определяются по напряжениям , и , вычисляемым по следующим формулам: (1) где γ - объемный вес, a - коэффициент бокового распора (давления) горных пород, рассчитываемый через коэффициент Пуассона. (2) . (3) Рис. 3. Схема откоса (а); последовательность построения следа наиболее вероятной поверхности разрушения методом [9] (б); СНВПР, построенные методом [9], при различных значениях коэффициента бокового давления грунта (в); СНВПР, построенные методом [9] (1) и методом [10] (2) при величине коэффициента бокового давления xо = 0,5 (г) (цитируется по работе [9]) Fig. 3. Schematic of the slope (a); sequence of construction of the most probable fracture surface trace by method [9] (б); SNVPR constructed by method [9] at different values of the soil lateral pressure coefficient (в); SNVPR constructed by method [9] (1) and method [10] (2) at value of lateral pressure coefficient xо = 0.5 (г) (cited from work [9]) Угол β находим из соотношения (4) где - угол откоса карьера (рис. 3, а). Угол между касательной к следу НВПР в рассматриваемой точке грунтового массива и направлением максимального нормального напряжения определяется выражением , (5) где - острый угол» [9]. Сам процесс построения следа наиболее вероятной поверхности разрушения подробно описан в работе [9]. На рис. 3, в приведены следы наиболее вероятных поверхностей разрушения в породном откосе, построенные автором работы [9] при различных значениях коэффициента бокового давления (распора) горной породы ξо (в дальнейшем будем придерживаться этого принятого в механике грунтов обозначения величины коэффициента бокового давления). Из рисунка видно, что положение наиболее вероятных поверхностей разрушения весьма существенно зависит от этой величины. Справедливо будет предположить, что и соответствующие величины коэффициентов запаса устойчивости K также будут зависеть от ξо. К сожалению, данные об этом в работе [9] отсутствуют. Сделав большой шаг вперед, предложив метод построения потенциально наиболее опасной поверхности разрушения на основе анализа НС грунтового массива, автор не пошел дальше, а использовал для определения величины коэффициента запаса устойчивости K одну из разновидностей методов первой группы [10]. Определение угла ориентации наиболее вероятной площадки разрушения на основе анализа напряженного состояния грунтового массива В работах В.К. Цветкова [11] и А.Н. Богомолова [12] приведены методы расчета устойчивости грунтовых массивов, основанных на анализе напряженного состояния грунтовых откосов и оснований сооружений. Напряженное состояние исследуемых объектов определяется либо методом конечных элементов (МКЭ), либо на основе аналитических решений задач линейной теории упругости методами теории функций комплексного переменного (МТФКП). Идея построения СНВПР авторами [11; 12] заимствована из работы [9], а вот углы ориентации наиболее вероятной площадки разрушения определяются совершенно другим образом. Ниже приведены два способа определения угла наклона наиболее вероятной площадки разрушения, которые подробно описаны в нашей работе [15]. Предложение В.К. Цветкова Величина коэффициента запаса устойчивости в точке грунтового массива Кt определяется выражением (6) [5], полученным на основании условия прочности Кулона [13; 14]: . (6) Если заменить напряжения sn и tn в формуле (6) напряжениями sy; sx; txy, получим выражение , (7) где sx; sz и txz - безразмерные (в долях γH) компоненты напряжения в точке грунтового массива; sсв = С (γHtgj)-1 - приведенное давление связности; H - высота откоса или глубина заложения фундамента. Очевидно, что условиями минимизации величины Kt по параметру α будут (8) Выполнение первого из условий (8) сводится к решению квадратного уравнения, корни которого определяются по формуле , (9) где (10) (11) Выбирать из полученных двух углов следует тот, при подстановке значения которого во второе условие (8) будет выполняться неравенство (12) Оказалось, что величина угла a, при котором значение коэффициента устойчивости в рассматриваемой точке Kt = Kt min соответствует знаку «+» перед радикалом в формуле (9). Из формул (7)-(11) следует, что угол aэкстр является функцией напряженного состояния, физико-механических свойств и характерного размера грунтового массива H. На рис. 4 приведены графические зависимости величины коэффициента устойчивости Kt в точке грунтового массива по наклонной площадке, угол наклона которой изменяется от 0 до 90°. Величины безразмерных (в долях γН) напряжений в точке приоткосной области определены методом конечных элементов и равны sx = 1,3; sy = 5,6; txy = 2, а физико-механические свойства грунта таковы, что давление связности sсв = 1,4, а угол внутреннего трения j = 28,6°. Рис. 4. Графическая зависимость вида Kt = f (a) Fig. 4. Graphical relationship Kt = f (a) Из рис. 4 видно, что экстремальное значение Kt соответствует величине aэкстр = 42°, полученной на основании выражений (8)-(12). Величина коэффициента запаса устойчивости в рассматриваемой точке, вычисленная по формуле (7), равна Кt = 1,08. Предложение А. Н. Богомолова Рассмотрим выражение , (13) где fост - остаточное сопротивление сдвигу [12]. Выразив нормальное и касательное напряжения, действующие по наклонной площадке, через компоненты напряжения sy; sx; txy и угол ее наклона a, подставив их в выражение (8), а затем применив к полученному выражению условия (8), найдем величину угла наклона площадки, которому соответствует минимальное значение величины остаточного сопротивления сдвигу fост min. При вычислениях используем те же значения напряжений и физико-механических свойств грунта, что указаны в предыдущем примере. На рис. 5 приведены кривые зависимостей величины коэффициента запаса устойчивости Kt, остаточного сопротивления сдвигу fост, удерживающих fуд и сдвигающих fсд сил, вычисленных в исследуемой точке приоткосной области. При этом величина угла ориентации площадки наиболее вероятного сдвига изменяется в интервале α Î[0°; 90°]. Рис. 5. Графические зависимости вида Kt = f (aэкстр); fост = f (aэкстр); fуд = f (aэкстр) и fсд = f (aэкстр) для рассматриваемой точки приоткосной области Fig. 5. Graphical relationships of the form Kt = f (aэкстр); fост = f (aэкстр); fуд = f (aэкстр) and fсд = f (aэкстр) for the considered near-surface point Из рис. 5, иллюстрирующего результаты проведенных вычислений, видно, что минимум значения остаточного сопротивления сдвигу fост соответствует площадке, ориентированной к оси абсцисс под углом aэкстр = 31о. Данное значение на 11o, или на 24 %, меньше значения угла aэкстр, определенного для тех же условий выше. Тем не менее значения соответствующих коэффициентов запаса устойчивости Kt для первого (Кt = 1,08) и второго (Кt = 1,044) примеров, вычисленных по формуле (7), отличаются менее чем на 4 %. Последнее обстоятельство объясняется обоснованным в работах [11; 12] положением о том, что в приоткосных областях и активных зонах фундаментов имеются равноустойчивые области, характеризующиеся тем, что величины коэффициентов запаса устойчивости, вычисленные для множества поверхностей разрушения, расположенных внутри этих зон, либо совпадают, либо отличатся друг от друга весьма незначительно. Таким образом, при потере устойчивости грунтовым массивом в предельное состояние переходит некая прослойка грунта, внутри которой располагается наиболее вероятная поверхность разрушения. Предложение, реализованное в работе [9] Известно, что нормальное sn и касательное τn напряжения, действующие в направлениях главных осей, могут быть выражены через соответствующие главные напряжения следующим образом: (14) (15) На рис. 6 в качестве примера представлена геометрическая интерпретация соотношения (14) в системе координат (; ) и соответствующие круги Мора. Рис. 6. Круги Мора а - для предельного и б - допредельного состояния в точке грунтового массива Fig. 6. Mohr circles а - for the limit and б - sub-limit conditions at a point in the soil mass Величина коэффициента запаса устойчивости в точке грунтового массива определяется отношением . (16) В предельном состоянии эта величина равна единице и растет по мере снижения нагрузки. При Кt ³ 1 нагрузка меньше критической для всех углов a. Поэтому величину Кt можно называть коэффициентом запаса устойчивости грунта в точке. Линия, в каждой точке которой выполняется равенство Kt = 1, также является границей области предельного состояния грунта. Числитель выражения (16) определен прочностными свойствами грунта, а знаменатель определяется фактическими напряжениями, действующими в рассматриваемой точке и физико-механическими свойствами грунтового массива. В допредельном состоянии , так как PK < BK.По мере роста напряжений (нагрузки) угол увеличивается. При достижении предельного напряженного состояния (рис. 3, а), выполняются условия PK = BK и вытекающие из выражения (15). Отрезок КВ указывает направление максимального касательного напряжения в момент достижения предельного напряженного состояния. Из рис. 6, б видно, что в предельном состоянии ÐАКВ = ÐAKN = соответственно. Обозначив через угол, составляемый направлением и направлением максимального касательного напряжения , получим, что , или . (17) Угол a не зависит от того, при каких значениях главных напряжений s1 и s2 наступает предельное состояние. Автор работы [11] считает, что выражение, определяющее величину угла a, справедливо не только для критической, но и для любой, меньшей по величине, нагрузки. Это означает, что для площадки, ориентированной под углом a, разность между удерживающей и сдвигающей силами будет заведомо меньше, чем для любой другой площадки, угол ориентации которой b ¹ a. Угол отклонения направления действия напряжения от оси x расчетной схемы равен углу , определяемому по формуле . (18) Таким образом, угол наклона вектора максимального касательного напряжения tmax, а значит, и площадки наиболее вероятного сдвига относительно оси OX расчетной схемы будет равен . (19) Вычисленное по формуле (19) значение угла используется, например, в работах [11; 12] для построения наиболее вероятной линии скольжения. Анализ этого выражения показывает, что величина угла ориентации площадки наиболее вероятного сдвига является функцией напряженного состояния и угла внутреннего трения. В результате выполненных вычислений по формуле (19) при численных значениях напряжений в исследуемой точке, которые использованы в двух предыдущих примерах, определены два значения угла ориентации площадки наиболее вероятного сдвига: aэкстр1 = 11° и aэкстр2 = 74°. Сравнивая полученные значения углов ориентации наиболее вероятной площадки сдвига с численными значениями соответствующих углов, полученных в первых двух случаях, можно видеть, что они отличаются в 2-3 раза. Кроме того из графиков, приведенных на рис. 4 и 5, видно, что при углах aэкстр1 = 11° и aэкстр2 = 74° численные значения коэффициентов запаса устойчивости Кt в рассматриваемой точке грунтового массива приблизительно равны 3 и 2,8, что существенно больше, чем при aэкстр = 42° (1,08 - первый случай) и aэкстр = 31° (1,044 - второй случай). Вычисление величины коэффициента запаса устойчивости откоса Вычисление величины коэффициента запаса устойчивости методом Л.Г. Фисенко [10] Выше было отмечено, что автор работы [9] предложил новый по сравнению с ранее известными метод построения следа наиболее опасной поверхности разрушения, который отличается от метода Л.Г. Фисенко, изложенного в [10, с. 300-304]. Сама же процедура вычисления величины коэффициента запаса устойчивости осталась такой же, как в [4, с. 296-300]. Анализируя способы построения СНВПР, изложенные в работах [9] и [10], можно сделать вывод, что они существенно разнятся, а следовательно, положение соответствующих СНВПР в грунтовом откосе должно быть разным. Это подтверждается рис. 3, г (цитируется по работе [9]), на котором приведены следы круглоцилиндрической поверхности разрушения, построенной по методу Л.Г. Фисенко [4] (пунктирная кривая 2), и наиболее вероятной поверхности разрушения (кривая 1), построенной методом [3] для однородного откоса с углом β = 40о и высотой Н = 200 м, сложенного грунтом со следующими физико-механическими характеристиками: γ = 2 т/м3, С = 50 т/м2, φ = 30о и ξо = 0,5 (γ, С, φ и ξо - соответственно объемный вес, удельное сцепление, угол внутреннего трения и коэффициент бокового давления грунта; размерность автора работы [9]). Как отмечено выше, вычисление величины коэффициента запаса устойчивости в работе [9] проведено методом [10], являющимся одной из разновидностей методов первой группы. Это значит, что поверхность разрушения принимается круглоцилиндрической, призма разрушения разбивается на вертикальные отсеки, вес каждого из которых принимается приблизительно равным вертикальной составляющей напряжений σn; напряжение τn определяется через σn по паспорту пород, а величина коэффициента запаса устойчивости откоса вычисляется по формуле [10] , (20) где , KL - сумма сил трения и сцепления, удерживающих откос; A - другие силы, удерживающие откос; - сумма касательных сил, сдвигающих откос и действующих по расчетной поверхности; B - другие силы, сдвигающие откос (обозначения и терминология работы [10]). Проведя вычисление коэффициентов запаса устойчивости для этих СНВПР методом [10], автор работы [9] установил, что их величины соответственно равны K1 = 1,75 и K2 = 1,6, причем K1 больше K2 на 9,4 %. Вычисление величины коэффициента запаса устойчивости методом В.К. Цветкова - А.Н. Богомолова Для сопоставления результатов, полученных в предыдущем параграфе, с результатами, получаемые методом Цветкова - Богомолова, проведены вычисления величины коэффициента запаса устойчивости того же объекта при помощи компьютерных программ [16; 17], в которых для определения напряжений используются решения основных граничных задач теории упругости методами теории функций комплексного переменного [12] (однородные откосы) и метод конечных элементов (неоднородные откосы). Отметим, что формализованный в этих программах метод позволяет учесть все три компоненты напряжения при построении следа наиболее вероятной поверхности разрушения (формулы (7)-(13)) и вычисления величины коэффициента запаса устойчивости, которая отыскивается по формуле , (21) где Fуд и Fсд - удерживающие и сдвигающие силы в точках СНВПР, определяемые соответственно числителем и знаменателем формулы (7), а координаты начальной точки СНВПР приняты такими же, как в работе [9], что делает одинаковыми условия всех расчетов. На рис. 7, а-в приведены изображения изолиний компонент напряжений, возникающих в исследуемом объекте от собственного веса пород при величине коэффициента бокового давления xо = 0,5. Из этих рисунков видно, что в месте перехода откоса в подошву наблюдаются концентрации всех трех компонент напряжения, что уже отмечено нами выше. Положение СНВПР также несколько отличается: на рис. 7, г-е ясно виден небольшой ниспадающий ее участок. В местах концентрации напряжений возникают области пластических деформаций (ОПД), которые определяются и по условию прочности Кулона [13; 14] (рис. 7, д), и при использовании решения смешанной задачи теории упругости и теории пластичности грунта [12; 18] (рис. 7, е). Однако эти области весьма малы, СНВПР их не пересекает, и они не оказывают какого-либо влияния на устойчивость откоса. Кроме того, на рис. 7, г показана ориентация площадок наиболее вероятного сдвига (ПНВС), которые построены на основе использования формул (9)-(12). В результате вычислений установлено, что величина коэффициента запаса устойчивости откоса K получилась равной 2,18. Сравнивая эту величину со значениями коэффициентов запаса устойчивости, вычисленных методом С.Н. Никитина [9] и методом Л.Г. Фисенко [10], видим, что она больше их соответственно на 20 и 26 %. Выводы Опираясь на анализ результатов расчетов для рассмотренного в работе примера, можно констатировать, что методы, основанные на анализе напряженного состояния грунтового массива, дают более высокие значения коэффициентов запаса устойчивости, чем традиционно применяемые. Так, величина коэффициента запаса устойчивости откоса K3, вычисленная методом В.К. Цветкова - А.Н. Богомолова, получилась равной 2,18. Сравнивая эту величину со значениями коэффициентов запаса устойчивости, вычисленных методом С.Н. Никитина [9] и методом Л.Г. Фисенко [10], которые соответственно равны K1 = 1,75 и K2 = 1,6, видим, что она больше их соответственно на 20 и 26 %. В это же время K1 > K2 на 9,4 % - это указывает на то, что учет напряжений при построении линии скольжения Рис. 7. Картины изолиний безразмерных (в долях γh) вертикальных σz (а), горизонтальных σx (б) и касательных τxz (в) напряжений, следов наиболее вероятных поверхностей разрушения и ПНВС - площадок наиболее вероятного сдвига (г); областей пластических деформаций, построенных в исследуемом откосе исходя из условия прочности Кулона [13; 14] и на основе решения смешанной задачи (с. з.) теории упругости и теории пластичности грунта [15], соответственно (д) и (е) Fig. 7. Plots of dimensionless (in fractions of γh) vertical σz (a), horizontal σx (б) and tangential τxz (в) stresses, traces of most probable fracture surfaces and most probable shear areas (г); areas of plastic deformation constructed in the studied slope based on Coulomb strength condition [13; 14] and on solution of mixed problem of elasticity theory and soil plasticity theory [15], respectively (д) and (е) делает ее форму отличной от круглоцилиндрической, что сказывается на величине коэффициента запаса устойчивости. Повышенные значения расчетной величины коэффициента запаса устойчивости позволяют в случае необходимости увеличивать углы заложения откосов грунтов сооружений, что дает существенный экономический эффект в связи с резким уменьшением объемов земляных работ.About the authors
O. A Bogomolova
Volgograd State Technical University
A. V Zhidelev
Branch «Vzlet» of Moscow Aviation Institute (National Research University)
References
- Цытович Н.А., Тер-Мартиросян З.Г. Основы прикладной геомеханики в строительстве. - М.: Высшая школа, 1981. - 320 с.
- Fellenius W. Calculation of the stability of earth dams // Transactions of II Congress on Large Dams. - Washington DC, 1933. - Vol. 4. - P. 445-462.
- Tscytbatarioff G. Foundations, retaining and earth structures. - N.Y.: McCraw-Hill Book Company, 1973.
- Соколовский В.В. Статика сыпучей среды. - Изд. 3-е, перераб. и доп. - М.: Физ.-мат. лит., 1960. - 121 с.
- Голушкевич С.С. Статика предельных состояний грунтовых масс. - М.-Л.: Гостехтеориздат, 1957. - 288 с.
- Wang F.D., Sun M.C., Ropchan D.M.Compyter program for pit slope stability analysis bee the finite element stress analysis and limiting equilibrium method // RJ 7685. - Burin of Mints, 1972.
- Desai C.S., Lightner J.G. Mixed fenite element procedure for Soil-Structure iteraction and construction seguences // Inter. J. for Numerical Mhetods in Engineering. - 1985. - Vol. 21, № 5.
- Федоров И.В. Методы расчета устойчивости откосов и склонов. - М.: Госстройиздат, 1962. - 202 с.
- Никитин С.Н. Построение ожидаемой поверхности скольжения по напряжения в бортах карьеров // Уголь. - 1962. - № 1. - С. 36-38.
- Фисенко Л.Г. Устойчивость бортов карьеров и отвалов. - М.: Недра, 1965. - 378 с.
- Цветков В.К. Расчет рациональных параметров горных выработок. - М.: Недра, 1993. - 251 с.
- Богомолов А.Н. Расчет несущей способности оснований сооружений и устойчивости грунтовых массивов в упругопластической постановке. - Пермь: ПГТУ, 1996. - 150 с.
- Coulomb C.A. Application des riles de maximus et minimis a quelques problemes de statique relatifs a L`architecture // Memories de savants strangers de L`Academlie des sciences de Paris, 1773. - 233 p.
- Coulomb C.A. Essay sure one application des regales des maxims et minims a quelques problems statique relatives, a la architecture // Memo. Acad. Roy. Div. Sav. - 1776. - Vol. 7. - P. 343-387.
- Богомолова О.А., Жиделев А.В. Determination of an inclination angle of the most probable destruction site at the ground point // Construction and Geotechnics. - 2020. - Т. 11, № 1. - C. 20-29. doi: 10.15593/2294-9826/2020.1.02
- FEA: св-во о гос. рег. программы для ЭВМ № 2015617889 / А.Н. Богомолов [и др.]. Зарег. 23 июля 2015 г.
- Устойчивость (напряженно-деформированное состояние): св-во о гос. рег. программы для ЭВМ № 2009613499 / А.Н. Богомолов [и др.]. Зарег. 30 июня 2009 г.
- Богомолов А.Н., Богомолова О.А. Сопоставление результатов физических и численных экспериментов по определению несущей способности однородного основания штампа // Основания, фундаменты и механика грунтов. - 2015. - № 6. - С. 7-11.
Statistics
Views
Abstract - 179
PDF (Russian) - 298
Refbacks
- There are currently no refbacks.