Анализ устойчивости склона — это статический или динамический, аналитический или эмпирический метод оценки устойчивости склонов плотин из грунта и камня, насыпей, вырытых склонов и естественных склонов в почве и скале. Он выполняется для оценки безопасной конструкции искусственных или естественных склонов (например , насыпей , дорожных выемок , открытых горных работ , выемок, свалок и т. д.) и условий равновесия. [1] [2] Устойчивость склона — это сопротивление наклонной поверхности разрушению путем скольжения или обрушения. [3] Основными целями анализа устойчивости склона являются поиск опасных зон, исследование потенциальных механизмов разрушения, определение чувствительности склона к различным механизмам срабатывания, проектирование оптимальных склонов с учетом безопасности , надежности и экономичности , а также проектирование возможных мер по исправлению положения, например, барьеров и стабилизации . [1] [2]
Успешное проектирование склона требует геологической информации и характеристик участка, например, свойств почвы / скального массива, геометрии склона , условий грунтовых вод , чередования материалов из-за разломов , систем соединений или разрывов , перемещений и напряжений в соединениях, сейсмической активности и т. д. [4] [5] Наличие воды оказывает пагубное влияние на устойчивость склона. Давление воды, действующее в поровых пространствах, трещинах или других разрывах в материалах, из которых состоит склон карьера, снизит прочность этих материалов. [6] Выбор правильного метода анализа зависит как от условий участка, так и от потенциального режима разрушения, при этом тщательно учитываются различные сильные и слабые стороны и ограничения, присущие каждой методологии . [7]
До компьютерной эры анализ устойчивости выполнялся графически или с помощью ручного калькулятора. Сегодня у инженеров есть много возможностей использовать программное обеспечение для анализа , от простых методов предельного равновесия до подходов вычислительного предельного анализа (например, предельный анализ конечных элементов , оптимизация компоновки разрывов ) и сложных и замысловатых численных решений ( коды конечных / отдельных элементов). [1] Инженер должен полностью понимать ограничения каждого метода. Например, предельное равновесие является наиболее часто используемым и простым методом решения, но он может стать недостаточным, если склон разрушается из-за сложных механизмов (например, внутренняя деформация и хрупкое разрушение , прогрессирующая ползучесть , разжижение более слабых слоев почвы и т. д.). В этих случаях следует использовать более сложные методы численного моделирования . Кроме того, даже для очень простых склонов результаты, полученные с помощью типичных методов предельного равновесия, используемых в настоящее время (Бишоп, Спенсер и т. д.), могут значительно отличаться. Кроме того, сегодня все чаще используется концепция оценки риска . Оценка риска касается как последствий разрушения склона, так и вероятности разрушения (оба требуют понимания механизма разрушения). [8] [9]
Традиционные методы анализа устойчивости склонов можно разделить на три группы: кинематический анализ, анализ предельного равновесия и симуляторы камнепада . [10] Большинство компьютерных программ анализа устойчивости склонов основаны на концепции предельного равновесия для двух- или трехмерной модели. [11] [12] Двумерные сечения анализируются в предположении условий плоской деформации . Анализ устойчивости двумерной геометрии склонов с использованием простых аналитических подходов может дать важную информацию о первоначальном проектировании и оценке риска склонов.
Методы предельного равновесия исследуют равновесие грунтового массива, имеющего тенденцию скользить вниз под действием силы тяжести . Поступательное или вращательное движение рассматривается на предполагаемой или известной потенциальной поверхности скольжения под грунтовым или скальным массивом. [13] В проектировании скальных склонов методы могут быть весьма значимыми для простого разрушения блока вдоль отдельных разрывов. [10] Все эти методы основаны на сравнении сил , моментов или напряжений, сопротивляющихся перемещению массы, с теми, которые могут вызвать неустойчивое движение (возмущающие силы). Выходом анализа является коэффициент безопасности , определяемый как отношение прочности на сдвиг (или, альтернативно, эквивалентной меры сопротивления сдвигу или емкости) к напряжению сдвига (или другой эквивалентной мере), необходимому для равновесия. Если значение коэффициента безопасности меньше 1,0, склон неустойчив.
Все методы предельного равновесия предполагают, что прочность материалов на сдвиг вдоль потенциальной поверхности разрушения регулируется линейными ( Мора-Кулона ) или нелинейными соотношениями между прочностью на сдвиг и нормальным напряжением на поверхности разрушения. [13] Наиболее часто используемой вариацией является теория прочности на сдвиг Терцаги, которая гласит, что
где — прочность на сдвиг интерфейса, — эффективное напряжение ( — общее напряжение, нормальное к интерфейсу, — давление поровой воды на интерфейсе), — эффективный угол трения, — эффективное сцепление.
Методы срезов являются наиболее популярным методом предельного равновесия. В этом подходе почвенная масса дискретизируется на вертикальные срезы. [12] [14] Используется несколько версий метода. Эти вариации могут давать разные результаты (коэффициент безопасности) из-за разных предположений и граничных условий между срезами. [13] [15]
Местоположение интерфейса обычно неизвестно, но может быть найдено с помощью численных методов оптимизации. [16] Например, функциональный дизайн склона рассматривает критическую поверхность скольжения как место, где она имеет наименьшее значение коэффициента безопасности из диапазона возможных поверхностей. Широкий спектр программного обеспечения для устойчивости склона использует концепцию предельного равновесия с автоматическим определением критической поверхности скольжения.
Типичное программное обеспечение для оценки устойчивости склонов может анализировать устойчивость склонов слоистых грунтов, насыпей, выемок в земле и анкерных листовых конструкций . Также могут быть включены эффекты землетрясений, внешняя нагрузка , условия грунтовых вод, стабилизирующие силы (например, анкеры, геоармирование и т. д.).
Многие инструменты анализа устойчивости склонов используют различные версии методов срезов, такие как упрощенный метод Бишопа , обычный метод срезов ( метод шведского круга/Петтерсона/Феллениуса ), Спенсер , Сарма и т. д. Методы Сарма и Спенсера называются строгими, поскольку они удовлетворяют всем трем условиям равновесия: равновесию сил в горизонтальном и вертикальном направлении и условию равновесия моментов. Строгие методы могут давать более точные результаты, чем нестрогие методы. Упрощенный метод Бишопа или Феллениуса — это нестрогие методы, удовлетворяющие только некоторым условиям равновесия и делающие некоторые упрощающие предположения. [14] [15] Некоторые из этих подходов обсуждаются ниже.
Метод шведского круга скольжения предполагает, что угол трения почвы или скальной породы равен нулю, т. е . . Другими словами, когда угол трения считается равным нулю, эффективный член напряжения становится равным нулю, тем самым приравнивая прочность на сдвиг к параметру сцепления данной почвы. Метод шведского круга скольжения предполагает круговой интерфейс разрушения и анализирует параметры напряжения и прочности с использованием круговой геометрии и статики. Момент, вызванный внутренними движущими силами склона, сравнивается с моментом, вызванным силами, противодействующими разрушению склона. Если силы сопротивления больше движущих сил, склон считается устойчивым.
В методе срезов, также называемом OMS или методом Феллениуса, скользящая масса над поверхностью разрушения делится на ряд срезов. Силы, действующие на каждый срез, получаются путем рассмотрения механического (силового и моментного) равновесия для срезов. Каждый срез рассматривается сам по себе, а взаимодействия между срезами пренебрегаются, поскольку результирующие силы параллельны основанию каждого среза. Однако третий закон Ньютона не удовлетворяется этим методом, поскольку, как правило, результирующие слева и справа от среза не имеют одинаковой величины и не являются коллинеарными. [17]
Это позволяет выполнить простой расчет статического равновесия, учитывая только вес почвы, а также сдвиговые и нормальные напряжения вдоль плоскости разрушения. Для каждого слоя можно учитывать как угол трения, так и сцепление. В общем случае метода слоев силы, действующие на слой, показаны на рисунке ниже. Нормальные ( ) и сдвиговые ( ) силы между соседними слоями ограничивают каждый слой и делают задачу статически неопределимой, когда они включены в расчет.
Для обычного метода срезов результирующие вертикальные и горизонтальные силы равны
где представляет собой линейный фактор, который определяет увеличение горизонтальной силы с глубиной среза. Решение для дает
Далее, метод предполагает, что каждый срез может вращаться вокруг центра вращения и что баланс моментов вокруг этой точки также необходим для равновесия. Баланс моментов для всех срезов, взятых вместе, дает
где - индекс сечения, - плечи момента, а нагрузки на поверхности были проигнорированы. Уравнение момента можно использовать для решения для сдвигающих сил на границе раздела после подстановки выражения для нормальной силы:
Используя теорию прочности Терцаги и преобразуя напряжения в моменты, имеем
где - поровое давление. Коэффициент запаса прочности - это отношение максимального момента из теории Терцаги к расчетному моменту,
Метод модифицированного Бишопа [18] немного отличается от обычного метода срезов тем, что предполагается, что нормальные силы взаимодействия между соседними срезами коллинеарны, а результирующая сила сдвига между срезами равна нулю. Подход был предложен Аланом В. Бишопом из Имперского колледжа . Ограничение, введенное нормальными силами между срезами, делает задачу статически неопределимой. В результате для решения коэффициента безопасности приходится использовать итерационные методы. Было показано, что метод дает значения коэффициента безопасности в пределах нескольких процентов от «правильных» значений.
Коэффициент безопасности для равновесия момента в методе Бишопа можно выразить как
где
где, как и прежде, — индекс среза, — эффективное сцепление, — эффективный внутренний угол внутреннего трения, — ширина каждого среза, — вес каждого среза, — давление воды у основания каждого среза. Для решения необходимо использовать итерационный метод, поскольку фактор безопасности появляется как в левой, так и в правой части уравнения.
Метод Лоримера — это метод оценки устойчивости склонов в связных грунтах. Он отличается от метода Бишопа тем, что использует клотоидную поверхность скольжения вместо круга. Этот режим разрушения был определен экспериментально для учета эффектов цементации частиц. Метод был разработан в 1930-х годах Герхардтом Лоримером (20 декабря 1894 г. — 19 октября 1961 г.), учеником пионера геотехники Карла фон Терцаги .
Метод анализа Спенсера [19] требует компьютерной программы, способной выполнять циклические алгоритмы, но упрощает анализ устойчивости склона. Алгоритм Спенсера удовлетворяет всем равновесиям (горизонтальному, вертикальному и движущему моменту) на каждом срезе. Метод допускает неограниченные плоскости скольжения и, следовательно, может определять коэффициент безопасности вдоль любой поверхности скольжения. Жесткое равновесие и неограниченная поверхность скольжения приводят к более точным коэффициентам безопасности, чем, например, метод Бишопа или обычный метод срезов. [19]
Метод Сармы [20] , предложенный Сарадой К. Сармой из Имперского колледжа, представляет собой метод предельного равновесия, используемый для оценки устойчивости склонов в сейсмических условиях. Его также можно использовать для статических условий, если значение горизонтальной нагрузки принять за ноль. Метод позволяет анализировать широкий спектр обрушений склонов, поскольку он может учитывать многоклиновой механизм обрушения и, следовательно, не ограничивается плоскими или круглыми поверхностями обрушения. Он может предоставить информацию о коэффициенте безопасности или о критическом ускорении, необходимом для обрушения.
Предположения, сделанные рядом методов предельного равновесия, перечислены в таблице ниже. [21]
В таблице ниже показаны условия статического равновесия, которым удовлетворяют некоторые популярные методы предельного равновесия. [21]
Анализ устойчивости скального склона на основе методов предельного равновесия может учитывать следующие виды отказов:
Более строгий подход к анализу устойчивости склона — предельный анализ . В отличие от предельного равновесного анализа, который делает ad hoc, хотя часто и разумные предположения, предельный анализ основан на строгой теории пластичности. Это позволяет, среди прочего, вычислять верхние и нижние границы истинного коэффициента безопасности.
Программы, основанные на анализе пределов, включают:
Кинематический анализ изучает, какие режимы разрушения могут возникнуть в скальном массиве. Анализ требует детальной оценки структуры скального массива и геометрии существующих разрывов, способствующих неустойчивости блоков . [32] [33] Используется стереографическое представление ( стереосетей ) плоскостей и линий. [34] Стереосети полезны для анализа разрывных скальных блоков. [35] Программа DIPS позволяет визуализировать структурные данные с использованием стереосетей, определять кинематическую осуществимость скального массива и статистический анализ свойств разрывов. [32]
Анализ устойчивости скального склона может разрабатывать защитные меры вблизи или вокруг структур, находящихся под угрозой падения блоков. Симуляторы камнепада определяют пути перемещения и траектории нестабильных блоков, отделенных от поверхности скального склона. [36] Метод аналитического решения, описанный Хангром и Эвансом [37], предполагает, что скальный блок представляет собой точку с массой и скоростью, движущуюся по баллистической траектории с учетом потенциального контакта с поверхностью склона. Для расчета требуются два коэффициента восстановления, которые зависят от формы фрагмента, шероховатости поверхности склона, импульса и деформационных свойств, а также от вероятности определенных условий при данном ударе. [38]
Методы численного моделирования обеспечивают приблизительное решение проблем, которые в противном случае не могут быть решены обычными методами, например, сложная геометрия, анизотропия материала , нелинейное поведение, напряжения на месте. Численный анализ позволяет деформировать и разрушать материал, моделировать поровое давление , деформацию ползучести , динамическую нагрузку, оценивать эффекты изменений параметров и т. д. Однако численное моделирование ограничено некоторыми ограничениями. Например, входные параметры обычно не измеряются, и доступность этих данных, как правило, плохая. Пользователь также должен знать о граничных эффектах, ошибках сетки, аппаратной памяти и временных ограничениях. Численные методы, используемые для анализа устойчивости склона, можно разделить на три основные группы: континуальное , дискретное и гибридное моделирование. [39]
Моделирование континуума подходит для анализа склонов грунта, массивных нетронутых скальных пород или сильно сочлененных скальных массивов. Этот подход включает методы конечных разностей и конечных элементов , которые дискретизируют всю массу до конечного числа элементов с помощью сгенерированной сетки (рис. 3). В методе конечных разностей (FDM) решаются дифференциальные уравнения равновесия (т. е. соотношения деформация-смещение и напряжение-деформация). Метод конечных элементов (FEM) использует приближения к связности элементов, непрерывности смещений и напряжений между элементами. [40] Большинство числовых кодов позволяют моделировать дискретные трещины , например, плоскости напластования , разломы . Обычно доступно несколько конститутивных моделей, например , упругость , упругопластичность, размягчение деформаций, упруговязкопластичность и т. д. [39]
Дискретный подход полезен для скальных склонов, контролируемых поведением разрыва. Скальный массив рассматривается как совокупность отдельных взаимодействующих блоков, подвергающихся внешним нагрузкам и, как предполагается, подвергающихся движению со временем. Эта методология в совокупности называется методом дискретных элементов (DEM). Дискретное моделирование допускает скольжение между блоками или частицами. DEM основано на решении динамического уравнения равновесия для каждого блока повторно до тех пор, пока не будут удовлетворены граничные условия и законы контакта и движения . Дискретное моделирование относится к наиболее часто применяемому численному подходу к анализу скальных склонов, и существуют следующие вариации DEM: [39]
Подход отдельных элементов описывает механическое поведение как разрывов, так и твердого материала. Эта методология основана на законе силы-смещения (определяющем взаимодействие между деформируемыми блоками породы) и законе движения (определяющем смещения, вызванные в блоках неравновесными силами). Соединения рассматриваются как [граничные условия. Деформируемые блоки дискретизируются во внутренние элементы постоянной деформации. [39]
Дискретная программа UDEC [41] (Универсальный код отдельных элементов) подходит для скальных склонов с высокой трещиноватостью, подверженных статической или динамической нагрузке. Двумерный анализ механизма трансляционного разрушения позволяет моделировать большие смещения, моделировать деформацию или текучесть материала. [41] Трехмерный дискретный код 3DEC [42] содержит моделирование множественных пересекающихся разрывов и поэтому подходит для анализа неустойчивости клиньев или влияния опоры горных пород (например, анкерных болтов, тросов). [39]
В анализе разрывных деформаций (DDA) смещения неизвестны, и уравнения равновесия затем решаются аналогично методу конечных элементов . Каждая единица сетки конечного элемента представляет собой изолированный блок, ограниченный разрывами. Преимуществом этой методологии является возможность моделирования больших деформаций, движений жесткого тела, состояний сцепления или отказа между блоками горных пород. [39]
Прерывистый скальный массив может быть смоделирован с помощью методологии отдельных элементов в форме кода потока частиц , например, программы PFC2D/3D . [43] [44] Сферические частицы взаимодействуют посредством фрикционных скользящих контактов. Моделирование совместных ограниченных блоков может быть реализовано посредством заданных прочностей связи. Закон движения многократно применяется к каждой частице, а закон силы-смещения — к каждому контакту. Методология потока частиц позволяет моделировать зернистый поток, разрушение неповрежденной породы, переходные движения блоков, динамическую реакцию на взрыв или сейсмичность, деформацию между частицами, вызванную сдвиговыми или растягивающими силами. Эти коды также позволяют моделировать последующие процессы разрушения склона скалы, например, моделирование скалы [39]
Гибридные коды включают в себя объединение различных методологий для максимизации их ключевых преимуществ, например, анализ предельного равновесия в сочетании с конечноэлементным анализом потока грунтовых вод и анализом напряжений; связанный анализ потока частиц и конечно-разностный анализ; гидромеханически связанные методы конечных элементов и материальных точек для моделирования всего процесса оползней, вызванных ливнями. [45] Гибридные методы позволяют исследовать отказы склонов трубопроводов и влияние высокого давления грунтовых вод на отказ слабого склона скалы. Связанные конечно-отдельно-элементные коды обеспечивают моделирование как поведения неповрежденной породы, так и развития и поведения трещин.