CALCULATION OF THE «BASE - EMBANKMENT» SYSTEM STABILITY
- Authors: Bogomolova O.A1, Zhidelev A.V2
- Affiliations:
- Volgograd State Technical University
- Moscow Aviation Institute (National Research University)
- Issue: Vol 12, No 4 (2021)
- Pages: 19-36
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/CG/article/view/1936
- DOI: https://doi.org/10.15593/2224-9826/2021.4.02
- Cite item
Abstract
The article describes a method of calculating the stability of fill slopes, including the construction of the most probable sliding surface and the determination of plastic deformations areas with elastic and elastoplastic (mixed task) distribution of stresses in the soil massif. The comparative stability calculations results of the "embankment - weak foundation" systems are presented, the data for which were taken from scientific articles published by other authors. According to these calculation results, the stress and displacement fields arising in the "embankment - weak base" system, the position and shape of the areas of plastic deformations constructed both on the basis of the Coulomb - Mohr strength condition and on the basis of solving a mixed problem of the theory of elasticity and the theory of plasticity of soils the article has comparison of numerical values of the stability coefficients and the safety coefficients of a weak. It is established that the replacement of the soil mass of the embankment with an equivalent load, as is often done in engineering calculations, significantly transforms the stress and displacement fields in the soil mass of the base, the position and shape of the areas of plastic deformation and, ultimately, entails obtaining unreliable results. The article shows the results of calculations performed using the FEA computer program developed by one of the authors and the Plaxis software package are in satisfactory agreement: the stress and displacement isofields practically coincide, plastic areas have similar dimensions, shape and position in the soil mass, which allows us to recommend the FEA computer program for practical use in calculating the stability of slopes and slopes, as well as the "embankment - weak foundation" system.
Full Text
Введение В настоящей работе приведены результаты вычисления устойчивости системы «основание - насыпь», выполненные при помощи двух подходов: на основании использования упругопластической модели Кулона - Мора [1; 2] и решения смешанной задачи теории упругости и теории пластичности грунтов [3; 4]. Для использования этих подходов необходимо знать четыре параметра грунта: два деформационных - модуль деформации Е0 и коэффициент Пуассона ν (или коэффициент бокового давления x0), и два прочностных - угол внутреннего трения φ и удельное сцепление с. Напряжения в исследуемом грунтовом массиве определяются аналитическими или численными методами. В работах [5-8] приведены результаты расчетов устойчивости системы «основание - насыпь», полученные независимо от нас другими авторами. В настоящей статье представлены результаты расчета тех же объектов на основе методики, изложенной в [7; 8], и выполнено сравнение наших результатов и полученных независимо от нас. Предлагаемый метод расчета устойчивости грунтового массива Считаем, что напряженное состояние исследуемого грунтового массива определено одним из численных (метод конечных расчетов (МКЭ)) или аналитических методов, например, методами теории функций комплексного переменного. Запишем условие прочности Кулона - Мора [1; 2] в виде, предложенным А. Како [8], введя предварительно некую функцию напряженного состояния и физико-механических свойств грунта К в рассматриваемой точке грунтового массива (1) где - приведенное нормальное напряжение [6]; sсв - давление связности. При K = 1 выражение (1) принимает вид условия почности Кулона - Мора. Рассмотрим равновесие бесконечно малой призмы единичной высоты (рис. 1, а), расположенной в некоторой точке грунтового массива, и составим уравнения равновесия проекций всех сил на направления (2) Выразим нормальную sn и касательную tn компоненты напряжения, действующего по наклонной площадке, ориентированной под углом a, через составляющие sz; sx; txz полного напряжения в рассматриваемой точке (3) Подставляя выражения (3) в формулу (1) и выделяя величину К, получим выражение для функции K, которую будем называть коэффициентом запаса устойчивости в точке грунтового массива: . (4) Числитель и знаменатель полученной дроби равны соответственно удерживающей Fуд и сдвигающей Fсд силам, действующим в рассматриваемой точке по площадке сдвига, ориентированной под углом α (см. рис. 1, а). Отметим, что в работах [6; 9] величина, эквивалентная К, называется коэффициентом стабильности. При K = 1 и при выполнении всех преобразований в обратном порядке выражение (4) вырождается в условие прочности Кулона - Мора. Рис. 1. Расчетные схемы для составления уравнений равновесия бесконечно малой треугольной призмы (а); решения смешанной задачи (б-г): ОУД - область упругих деформаций; ОПД - область пластических деформаций Fig. 1. Calculation schemes for composing the equilibrium equations of an infinitesimal triangular prism (а); solving a mixed problem (б-г) Задача об определении угла ориентации наиболее вероятной площадки сдвига решается путем исследования функции (4) на экстремум согласно тривиальной процедуре, заключающейся в выполнении условий (5) Выполнение первого из условий (5) сводится к решению квадратного уравнения, корни которого определяются выражением , (6) где ; (7) (8) Выбирать из полученных двух углов следует тот, при подстановке значения которого во второе условие (5) будет выполняться неравенство . (9) В результате вычислений установлено, что величина угла a, при котором K = Kmin, соответствует знаку «+» перед радикалом в формуле (6). Из формул (4), (6)-(9) следует, что угол a является функцией напряженного состояния и физико-механических свойств (в том числе и коэффициента бокового давления xо) слагающих грунтов. Отметим, что в выражения (2)-(9) входят безразмерные напряжения (в долях γH, где γ - удельный вес грунта; H - характерный геометрический размер исследуемой области, например глубина зложения фундамента или высота откоса). Построение наиболее вероятной линии скольжения (НВЛС), точнее говоря, следа наиболее вероятной поверхности скольжения, проводится согласно методике, приведенной в работе [10] и усовершенствованной в работах [3; 11; 13]. Согласно этой методике при построении НВЛС в каждой ее точке выполняется условие Kst = Kst min. Если S - дуговая координата на НВЛС длиной l, то величина коэффициента запаса устойчивости грунтового массива, естественно, при неизменном положении наиболее вероятной линии скольжения определяется выражением (10) где Fуд и Fсд - удерживающие и сдвигающие силы, действующие в каждой точке НВЛС, численные значения которых определяются соответственно числителем и знаменателем фомулы (4). Построение областей пластических деформаций Использование условия прочности Кулона - Мора Если задача об определении напряженного состояния грунтового массива решается в рамках линейной теории упругости, то в подавляющем большинстве случаев для отыскания положения, формы и размера областей пластических деформаций (областей сдвига) используют условие прочности Кулона - Мора, записанное в той или иной форме. Выше показано, что выражение (4) при K = 1 представляет собой условие прочности Кулона - Мора, в которое входит угол ориентации площадки наиболее вероятного сдвига, определяемый на основе выражений (6)-(9). Таким образом, для определения областей пластических деформаций на рассматриваемом этапе нагружения следует определить точки грунтового массива, в которых выражение (4) принимает численное значение, равное 1. Соединяя эти точки плавной кривой, получим очертание области пластических деформаций. Решение смешанной задачи теории упругости и теории пластичности грунта В работах [3; 4] показано, что, если некоторая точка грунтового массива находится в области пластических деформаций, то наряжения в этой точке могут быть определены по следующим формулам: (11) где ; (см. рис. 1, б); угол α определяется выражениями (6)-(9). Для определения границ между областями упругих и пластических деформаций воспользуемся очевидным условием непрерывности поля напряжений: нормальные и касательные к границе напряжения в каждой ее точке должны быть одинаковыми по обеим ее сторонам (см. рис. 1, в, г). Это условие будет определять как напряженное состояние самого грунта, так и форму областей пластических деформаций [3; 5]. Пусть кривая АА¢ есть граница между упругой и пластическими областями, и точка В лежит на этой кривой. Вырежем мысленно в окрестности точки В бесконечно малую прямоугольную призму единичной высоты ДЕКМ, расположенную таким образом, что точка В лежит на ее диагонали, а сама диагональ ДК является касательной к кривой АА¢ в точке В (рис. 1, в). Разделим прямоугольную призму ДЕКМ на две треугольные призмы ДЕК и ДМК таким образом, что первая из них будет находиться в области упругих деформаций, а вторая - в пластической области. Уравнения равновесия выполняются во всем объеме грунтового массива. Выразим нормальные и касательные составляющие полного напряжения, действующего по наклонным граням призм ДЕК и ДМК через компоненты полного напряжения sz; sx; txz и s¢z, s¢x; t¢xz, используя тригонометрические функции двойного аргумента (12) (13) Из условия неразрывности поля напряжений следует, что . (14) Подставляя в (14) соответствующие значения из выражений (12) и (13) и имея в виду, что sz = s¢z, s¢x; t¢xz определяются формулами (11), получим выражения, описывающие границу упругих и пластических областей: (15) где q - угол между касательной к границе АА¢ в точке В и положительным направлением оси ОХ; jq - «граничное» значение угла внутреннего трения грунта, при котором данная точка массива «переходит» на данном этапе нагружения в область пластических деформаций; b и l - см. комментарий к формулам (11); sz; sx; txz - компоненты напряжения, найденные из аналитического или численного решения соответствующей задачи теории упругости для исследуемой области на рассматриваемом этапе нагружения. Аппарат для проведения вычислений Все вычисления, выполненные при подготовке настоящей статьи, проведены при помощи компьютерной программы [14], разработанной при непосредственном участии авторов. Для определения напряжений в исследуемых областях в этих программах формализованы метод конечных элементов и аналитическое решение второй основной краевой задачи теории упругости для весомой изотропной полуплоскости с криволинейной границей [3; 15], полученное на основе использования методов теории функций комплексного переменного. Примеры расчета устойчивости системы «основание - насыпь» Расчет устойчивости системы «водонасыщенное основание - насыпь» в стабилизированном состоянии Для проведения этого расчета выбран пример из работы [5], в которой, в частности, рассматривается устойчивость системы «водонасыщенное основание - дамба» в стабилизированном и нестабилизированном состоянии, т.е. с учетом и без учета фильтрационных процессов. При решении этой задачи об устойчивости системы «основание - дамба» в стабилизированном состоянии авторами [5] приняты следующие характеристики грунта: для дамбы - γ = 16 кН/м3; φ = 20°; с = 20 кПа; Е = 20 МПа; ν = 0,3, что соответствует величине коэффициента бокового давления xо = 0,43; для основания - γ = 20 кН/м3; φ = 17°; с = 17 кПа; Е = 10 МПа; ν = 0,4, что соответствует величине xо = 0,66. Расчетная область основания дамбы высотой 10 м, шириной по верху 20 м и по основанию 40 м имеет размеры 40 ´ 120 м, что исключает влияние на напряженно-деформированное состояние (НДС) граничных условий. Расчеты выполнены методом конечных элементов при помощи программы Plaxis, в которой формализован МКЭ. Авторы цитируемой статьи рассматривают проблемы взаимодействия дамбы и ее однородного водонасыщенного основания в рамках плоской задачи линейной теории упругости для анализа локальной и общей устойчивости основания и дамбы, что важно при выборе экономически эффективных конструкций дамб. Этот выбор зачастую сводится к оценке напряженного состояния основания под воздействием эквивалентной трапециевидной нагрузки [15]. Авторами [5] показано, что результаты расчетов с учетом и без учета жесткости тела дамбы существенно отличаются друг от друга. Критериями оценки этих различий могут, по мнению авторов [5], служить картины изолинии напряжений σx, σy, σxy, τmax и безразмерной величины назваемой авторами относительной прочностью (отметим, что точки грунтового массива, в которых выполняется условие принадлежат границе «кулоновской» пластической области). Кроме того, авторами [5] введены следующие обозначения: τ - действующее касательное напряжение; - предельное значение сопротивления сдвигу. Поэтому во избежание путаницы в тексте будем придерживаться этих обозначений. Нами выполнены все расчеты, результаты которых приведены в работе [5], и дополнительно построены картины изолиний полей перемещений и напряжений, возникающих в системе «насыпь - дамба» при учете и неучете жесткости последней. Отметим, что эпюры распределенной нагрузки при замене тела дамбы соответствующими силами (при неучете жесткости дамбы) показаны на рисунках красным цветом, тело дамбы окрашено голубым, а ее основание желтым цветом. Проанализируем эти изображения. На рис. 2 приведены картины изолиний безразмерных (в долях γдНд) вертикальных σz; горизонтальных σx и касательных τzx напряжений в системе «дамба - основание» без учета (а; в; д) и с учетом жесткости дамбы (б; г; е) соответственно. Сравнивая изолинии соответствующих полей напряжений, делаем вывод о том, что численные значения соответствующих напряжений в точках основания с одинаковыми координатами существенно отличаются друг от друга. Причем отличие численных значений вертикальных σz и горизонтальных σx напряжений может достигать 20-25 %, а касательных τzx напряжений 300 % и более. Аналогичный результат получен З.Г. Тер-Мартиросяном и соавторами, которые утверждают, что, «как и следовало ожидать, в случаях учета и неучета жесткости дамбы получается существенная разница в НДС ее основания» [5]. В качестве еще одного подтверждения полученного результата приведем картины изополей напряжений τmax (рис. 3). Рассматривая их, делаем вывод, аналогичный вышеприведенному. а б в г д е Рис. 2. Изолинии безразмерных (в долях γдНд) вертикальных σz; горизонтальных σx и касательных τzx напряжений в системе «дамба - основание» без учета (а; в; д) и с учетом жесткости дамбы (б; г; е) соответственно Fig. 2. Isolines of dimensionless stresses (in fractions of γдНд): vertical - σz; horizontal - σx and tangential - τzx in the “dam - base” system without (а; в; д) and taking into account (б; г; е) the rigidity of the dam, respectively Изучая поля перемещений, возникающих в основании дамбы, мы получили результат, несколько отличающийся от результатов исследования, в котором сказано, что «расчеты показывают, что максимальная осадка основания в центре при его толщине 40 м составляет 32 см при отсутствии пригрузки. В случае замены дамбы эквивалентной нагрузкой осадка также составила 32 см - следовательно, учет жесткости не оказывает влияние на осадку основания» [5]. Наши же расчеты показали (рис. 4), что вертикальное перемещение (осадка) дамбы в той же точке (точке пересечения основания дамбы с вертикальной осью ее симметрии) при условии, что отношение модулей деформации материала дамбы и основания Eод/Eо осн = 2, как это определено в работе [5], равно 19,5 см. Если Eод/Eосн = 1, то величина осадки равна 9,7 см. Если тело дамбы заменено эквивалентной вертикальной нагрузкой, то величина осадки той же точки равна всего 1,1 см, т.е. соответствующие напряжения отличаются по величине более чем на порядок. Из рис. 3 видно, что существенные отличия имеет и распределение горизонтальных перемещений. Таким образом, мы не можем согласиться с выводами авторов работы [5] о том, что «учет жесткости дамбы не оказывает влияние на осадку основания». а б в г Рис. 3. Изополя τmax в основании и в теле дамбы с учетом (а, в) и без учета (б, г) жесткости дамбы, полученные авторами работы [5] (а, б) и нами (в, г) при помощи компьютерной программы [14] Fig. 3. Isofields τmax in the base and in the dam body, taking (а, в) and without taking into account (б, г) the dam rigidity, obtained by the authors of [5] (а, б) and by proposed method (в, г) using a computer program [14] а б в г Рис. 4. Изолинии безразмерных (в долях Н) вертикальных (а; в) и горизонтальных (а; г) упругих перемещений системы «основание - дамба» с учетом (а; б) и без учета (в; г) жесткости дамбы Fig. 4. Isolines of dimensionless (in fractions of H) vertical (а; в) and horizontal (а; г) elastic displacements of the “base - dam” system taking (а; б) and without taking into account (в; г) the dam rigidity Следующим изученным нами вопросом является вопрос о локальной устойчивости основания. Для решения этого вопроса авторами работы [5] и нами построены области пластических деформаций (области сдвигов, области предельного равновесия, как они названы в [5]) с учетом и без учета жесткости дамбы. Как сказано в работе [5], «в первом случае области предельного равновесия находятся под краями дамбы, что может привести к локальной потере устойчивости ее откоса. Во втором случае область предельного напряженного состояния локализуется под центром дамбы на определенной глубине и не выходит на поверхность» (рис. 5, а, б). Нами получен совершенно аналогичный результат, графическая интерпретация которого приведена на рис. 5, в, г. Отметим лишь одно отличие: в нашем случае пластические области, возникающие под краями дамбы (учитывается ее жесткость), имеют направление развития не под дамбу, а в противоположные стороны. Этот результат также подтверждает предположение о неэквивалентности рассматриваемых расчетных схем. а б в г Рис. 5. Изополя τrel = τ/τ* в основании и в теле дамбы с учетом (а) и без учета (б) жесткости тела дамбы, полученные авторами работы [5]; области пластических деформаций в основании и теле дамбы, ограниченные линиями τrel = τ/τ* = 1, построенные нами (в, г) при помощи компьютерной программы [14] с учетом (в) и без учета (г) жесткости тела дамбы Fig. 5. Isofields τrel = τ/τ* in the base and in the body of the dam, taking (а) and without taking into account (б) the rigidity of the body of the dam, obtained by the authors of the work [5]; areas of plastic deformations in the base and body of the dam, bounded by the lines τrel = τ/τ* = 1, constructed by authors (в, г) using a computer program [14] taking (в) and without taking into account (г) the dam body rigidity На рис. 6, а представлено изополе относительной прочности τrel = τ/τ* в начальный момент t0 = 0 и область пластических деформаций в основании и теле дамбы, ограниченная линиями τrel = τ/τ* = 1, построенная нами (рис. 6, б) при помощи компьютерной программы [14]. Оба изображения соответствуют результатам расчета, выполненного с учетом жесткости тела дамбы. Сравнивая эти рисунки, определяем их практическую тождественность. а б Рис. 6. Изополе относительной прочности τге1 = τ/τ* в водонасыщенном основании в момент времени t0 = 0, полученное авторами работы [5] (а) и область пластических деформаций в основании и теле дамбы, ограниченная линиями τrel = τ/τ* = 1, построенная нами (б) при помощи компьютерной программы [14] с учетом жесткости тела дамбы Fig. 6. The relative strength isofield τге1 = τ/τ* in a water-saturated base at time t0 = 0, obtained by the authors of [5] (а) and the plastic deformations area in the base and body of the dam, bounded by the lines τrel = τ/τ* = 1, constructed by us (б) using a computer program [14] taking into account the dam body rigidity Последним этапом проведенных нами вычислений являются расчеты по определению коэффициента запаса устойчивости основания дамбы с учетом ее жесткости и без такового. В работе [5] эти коэффициенты запаса устойчивости определены равными соответственно т.е. отличаются друг от друга в 1,408 раза. К сожалению, в работе [5] этот результат никак не проиллюстрирован. На рис. 7, а, в приведены изображения следов наиболее вероятных поверхностей выпора грунта основания из-под дамбы, построенные с учетом и без учета ее жесткости, и «кулоновские» области пластических деформаций, построенные на основе выражений (4)-(9), как это описано выше. Вычисления, выполненные нами, дают следующие значения соответствующих коэффициентов запаса устойчивости - Как видно, численные значения этих коэффициентов отличаются друг от друга в 1,503 раза. Сравнивая численные значения этих коэффициентов с соответствующими значениями, приведенными в работе [5], видим, что величина больше величины на 7,1 %, а величина K больше величины на 12,9 %. Как видно, и в этом случае полученные нами и приведенные в работе [5] результаты очень близки. На рис. 7, б, г в качестве примера изображены области пластических дефрмаций (ОПД) в оснвании дамбы при условиях, описанных выше, но на основе приближенного смешанного решения задачи теории упругости и теории пластичности грунта с использованием формул (11), (15). Из рисунков видно, что эти области существенно отличаются от соответствующих «кулоновских областей». Отметим, что величина (соответствует рис. 7, б) на 18 % больше величины (соответствует рис. 6, а) из-за того, что длина участка находящегося внутри ОПД (рис. 6, б) существенно меньше длины участка CE, также находящегося внутри соответствующей пластической области (рис. 7, а). Рис. 7. Области пластических деформаций и следы наиболее вероятных поверхностей выпора, построенные с учетом (а; б) и без учета (в; г) жесткости дамбы на основании условия прочности Кулона - Мора (а; в) и решения смешанной задачи теории упругости и теории пластичности грунта (б; г) Fig. 7. Plastic deformations areas and the most probable outburst surfaces traces, constructed taking into account (а; б) and without taking into account (в; г) the dam rigidity based on of the Coulomb-Mohr strength condition (а; в) and solving a mixed problem of the theory of elasticity and the theory of soil plasticity (б; г) Расчет величины коэффициента стабильности системы «слабое основание - насыпь» Перейдем теперь к анализу результатов, приведенных в работе [6], где рассматривается пример расчета величины коэффициента стабильности системы «слабое основание - насыпь», который выбран из книги проф. Г.М. Шахунянца как наиболее подробно расписанный [9]. Высота насыпи составляет 30 м, крутизна ее откосов равна 1:2, а ширина основной площадки равна 10 м. Грунт основания имеет следующие характеристики: удельный вес 20 кН/м3, угол внутреннего трения φ = 12°, удельное сцепление с = 20 кН/м2. При расчете насыпь заменяется эквивалентной нагрузкой, эпюра которой имеет форму трапеции, а ее максимальная интенсивность равна Р = 600 кН/м2 [6]. Результаты вычислений, выполненные проф. Г.М. Шахунянцем [9], представлены на рис. 8, а; они указывают, что минимальное значение коэффициента стабильности Кстаб = 1,17. Область, оконтуренная изолинией со значением Кстаб = 1,25, располагается по оси насыпи на глубине от 15 до 30 м. Результаты численного моделирования такой схемы авторами работы [6] (для сопоставления насыпь была заменена трапецеидальной нагрузкой) дают весьма похожие результаты (рис. 8, б) и немного завышенный коэффициент стабильности (минимальное значение Кстаб = 1,25), что объясняется невозможностью полного сопоставления двух расчетных методик (численное моделирование использует модели грунтов с параметрами, не указанными в примере). Область минимальных значений располагается на глубине от 5 до 30 м. а б Рис. 8. Результаты расчета величины коэффициента стабильности для рассматриваемого примера: полученные проф. Г.М. Шахунянцем [9] (а) и авторами работы [6] (б) Fig. 8. The calculating results of the stability coefficient value for example under consideration: obtained by Prof. G.M. Shakhunyants [9] (а) and by actual article authors [6] (б) Аналогичные расчеты выполнены и нами при помощи компьютерной программы FEA [14] при величине коэффициента бокового давления грунта основания xо = 0,75, графичекая интерпретация которых представлена на рис. 8. Стоит отметить, что расчет коэффициента стабильности при приложении эквивалентной трапецеидальной нагрузки дает некорректные результаты, поскольку не учитывает, из какого грунта отсыпается земляное полотно. Сравнивая изополя Кстаб, приведенные на рис. 8, б и 9, можно говорить об их практической идентичности. Рис. 9. Изополя коэффициента стабильности, полученные при помощи компьютерной программы FEA [14] при величине коэффициента бокового давления грунта основания xо = 0,75 Fig. 9. The stability coefficient isofields obtained using the FEA computer program [14] at the value of the lateral pressure coefficient of the base soil xо = 0,75 Однако если отступить от допущений, принятых в инженерном расчете, и задать насыпь не трапецеидальной нагрузкой, а супесчаным грунтом с характеристиками: удельный вес 20 кН/м3, угол внутреннего трения φ = 24°, удельное сцепление c = 19 кН/м2, то картина распределения коэффициента стабильности в основании будет совершенно иная (рис. 10). На рис. 11 приведены изображения следов наиболее вероятных поверхностей скольжения (НВПС) для рассматриваемого примера, построенные авторами работы [6] (рис. 11, а) и нами (рис. 11, б, в) при помощи компьютерной программы FEA, из которых видно, что а б Рис. 10. Изополя Кстаб по данным работы [6] (а) и построенные при помощи компьютерной программы FEA для измененных условий работы [6] (б) Fig. 10. Isofields of Кстаб according to article [6] (а) and constructed using the FEA computer program for modified working conditions [6] (б) а б в Рис. 11. Следы наиболее вероятных поверхностей скольжения а - по данным работы [6], и б - построенные при помощи компьютерной программы FEA [14] для условий работы [6] при упругом решении, и в - при упругопластическом распределении напряжений (смешанная задача) в грунтовом массиве Fig. 11. The most probable sliding surfaces traces а - according to the data of acticle [6] and б - constructed with the help of the FEA computer program [14] for working conditions [6] with elastic solution and в - elastic-plastic stress distribution (mixed problem) in the soil massif форма следов НВПС практически одинакова в обоих случаях, но их положение в грунтовом массиве несколько отличается - следы НВПС, построенные нами, расположены несколько ближе к контуру насыпи, чем это получено авторами работы [6]. Отличаются и численные значения коэффициентов запаса устойчивости. В нашем случае при упругом распределении напряжений в теле насыпи величина коэффициента запаса устойчивости равна К = 1,34, т.е. отличается от величины Куст на 33 %. Если задача об устойчивости насыпи решается в смешанной постановке, то величина Ксм = 1,08, т.е. практически совпадает с величиной Куст, приведенной в работе [6]. Это объясняется тем, что в данном случае след НВПС практически совпадает (см. рис. 11, в) с границей упругих и пластических областей, на которой выполняется условие К = 1. Расчет насыпи на слоистом глинистом основании Анализ устойчивости [7] проведен для двух конструкций земляного полотна: это конструкция с бермами и однородным откосом заложения 1:1,75 (см. соответственно правую и левую части рис. 12, а). Моделирование в программе Plaxis 2D проведено с использованием опции stage construction, что позволило задавать промежутки времени возведения каждого отдельного фрагмента (уровня) насыпи. Графическая интерпретация результатов моделирования приведена на рис. 12, б, в. а б Рис. 12. Две конструкции земляной насыпи: с бермами и однородным откосом заложения 1:1,75 (соответственно правая и левая части а), конечно-элементная схема для работы в оболочке программы FEA (б) Fig. 12. Two earth embankment structures: with berms and a uniform slope of 1:1.75 (respectively, the right and left parts of a), a finite element scheme for operation in the FEA program shell (б) Оценка устойчивости насыпи в программном комплексе ZSoil выполнена с использованием моделей грунта Мора - Кулона и Друкера - Прагера. Численные значения величин коэффициентов запаса устойчивости сведены в обоих случаях в таблицу. Численные значения расчетных коэффициентов запаса устойчивости The numerical values of calculated stability margin coefficients Высота насыпи и заложение откоса Величина коэффициента запаса устойчивости Bishop Fellinius Шахунянц Plaxis ZSoil GEO 5 FEM FEA [13] M-K M-K D-P M-K D-P 15,72 (1:1,75) 1,94 1,70 1,76 1,69 1,72 1,72 1,88 2,06 1,58 15,72 (бермы) 2,31 2,05 2,06 1,98 2,10 1,95 2,26 2,28 1,72 На рис. 12, б приведена конечно-элементная схема для работы в оболочке программы FEA, которая составлена нами на основе схемы, приведенной в работе [7], и цитируемой нами на рис. 12, а. На этом же рисунке приведены численные значения физико-механических свойств грунтов слоев основания (слева на схеме) и параметры самой конечно-элементной схемы (справа). На рис. 13 приведена часть результатов компьютерного моделирования устойчивости насыпи на слоистом глинистом основании, полученных авторами работы [7] и настоящей статьи. Сравнивая положение и форму следов НВПС, можно сказать, что они подобны. а б в Рис. 13. Результаты моделирования конструкции в программе FEA и Plaxis 2D: области пластических деформаций и следы наиболее вероятных поверхностей скольжения, построенные нами (а), и изополя сдвигающих перемещений для левой (б) и правой (в) части насыпи (цитируется по работе [7]) Fig. 13. Structural modeling results in the FEA and Plaxis 2D program: areas of plastic deformations and the most probable sliding surfaces traces constructed by authors (а), and isofields of shifting displacements for the left (б) and right (в) parts of the embankment (cited in [7]) Численные значения коэффициентов запаса устойчивости откосов насыпи, полученные при помощи программы FEA, наиболее близки величинам коэффициентов запаса, полученным на основе использования программных комплексов Plaxis 2D и Zsoil. Отличия составляют соответственно 6,6 % и 8,2 % для откоса с заложением 1:1,75, а для откоса с бермами - 13,1 % и 15,1 % (см. таблицу). Выводы 1. В результате проведения сопоставительных расчетов подтверждено, что замена грунтового массива насыпи эквивалентной нагрузкой при расчеие устойчивости ее основания ведет к получению заведомо недостоверных результатов. 2. Результаты, получаемые при помощи программного комплекса Plaxis и компьютерной программы FEA для одних и тех же объектов, совпадают с достаточной для практики степенью точности, что позволяет рекомендовать последнюю для решения задач об устойчивости системы «насыпь - слабое основание».About the authors
O. A Bogomolova
Volgograd State Technical University
A. V Zhidelev
Moscow Aviation Institute (National Research University)
References
- 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. Essai sur une application des regles des maximis et minimis a quelquels problemesde statique relatifs, a la architecture // Mem. Acad. Roy. Div. Sav. - 1776. - Vol. 7. - P. 343-387.
- Богомолов А.Н. Расчет несущей способности оснований сооружений и устойчивости грунтовых массивов в упругопластической постановке. - Пермь: ПГТУ, 1996. - 150 с.
- Богомолов А.Н., Богомолова О.А. Сопоставление результатов физических и численных экспериментов по определению несущей способности однородного основания штампа // Основания, фундаменты и механика грунтов. - 2015. - № 6. - С. 7-11.
- Тер-Мартиросян З.Г., Сидоров В.В., Олодо Т.Д. Напряженно-деформированное состояние дамбы и ее основания с учетом их взаимодействия // Инженерная геология. - 2011. - № 6. - С. 30-34.
- Вавринюк Т.С., Федоренко Е.В. К вопросу об определении допустимой нагрузки земляного полотна на слабые основания // Путь и путевое хозяйство. - 2015. - № 10. - С. 10-14.
- Лунев А.А., Сиротюк В.В., Безделов Н.С. Сопоставление методов оценки устойчивости высокой насыпи из золошлаковой смеси // Вестник СибАДИ. - 2016. - Вып. 5 (51). - С. 106-113.
- Caquot A. Equilibrium des Massifs a frottement Interne. - Gauthier Villars, Paris, France, 1934. - P. 1-91.
- Шахунянц Г.М. Земляное полотно железных дорог. - М.: Трансжелдориздат, 1953. - 827 с.
- Никитин С.Н. Построение ожидаемой поверхности скольжения по напряжениям в бортах карьера // Уголь. - 1962. - Вып. 1. - С. 36-38.
- Цветков В.К. Расчет устойчивости откосов и склонов. - Волгоград: Нижне-Волжское кн. изд-во, 1979. - 238 с.
- Богомолова О.А., Жиделев А.В. Determination of an inclination angle of the most probable destruction site at the ground point // Construction and Geotechnics. - 2020. - Т. 11, № 1. - С. 20-29.
- Горбунов-Посадов М.И. Устойчивость фундаментов на песчаном основании. - М.: Госстройиздат, 1962. - 96 с.
- FEA / А.Н. Богомолов и др. // Свидетельство о государственной регистрации программы для ЭВМ № 2015617889 от 23 июля 2015 г.
- Богомолов А.Н., Ушаков А.Н. Методы теории функций комплексного переменного в задачах геомеханики. - Волгоград: ВолгГАСУ: Изд-во ВГСПУ «Перемена», 2014. - 227 с.
- Тер-Мартиросян З.Г. Механика грунтов. - М.: АСВ, 2009. - 550 с.
Statistics
Views
Abstract - 964
PDF (Russian) - 164
Refbacks
- There are currently no refbacks.