SIMULATION OF SOIL LIQUEFACTION DUE TO EARTHQUAKE LOADING USING THE UBC3D-PLM MODEL
Abstract and keywords
Abstract (English):
Liquefaction is a phenomenon in which the stiffness and strength of a soil are reduced as a result of seismic effect or other dynamic effects. Liquefaction was the basic reason of the big damages caused by many earthquakes around the world. The basic step in the processes of predicting the soil liquefaction is the modeling of soil behavior. At the present time, numerous soil models are presented. Nevertheless, only some of them can simulate this process. Model UBC3D-PLM is one of these models which can be used. In this paper, the possibilities of this model are considered by modeling on the PLAXIS software package the seismic impact on a building with its different heights. The actual data of Upland earthquake 1990 near Los Angeles city was used. Results of this simulation showed us the difference in the behavior of the soil mass under the impact of an earthquake compared with the elastic behavior, as well as showed us the necessary to use the UBC3D-PLM model to estimate the seismic impact.

Keywords:
Liquefaction, dynamic Impact, UBC3D-PLM, pore pressure, earthquake, seismic impact
Text
Text (PDF): Read Download

Введение

Землетрясения — это колебания участка земной коры, вызванные прохождением сейсмических волн, возникающих в результате кратковременного выделения большого количества энергии. Характерной особенностью воздействия землетрясений является их кратковременность, большая величина амплитуды смещения в момент воздействия и сравнительно малые амплитуды последующих колебаний. Волны, вызванные воздействием землетрясений, распространяясь по частицам грунта и находящейся в порах жидкости, могут привести к нарушению структуры массива грунта, а сравнительно длительный процесс к перекомпоновке частиц грунта под действием его собственного веса и внешней нагрузки. Характерной особенностью водонасыщенного несвязанного и слабосвязанного грунта является его возможность разжижаться. Несвязный водонасыщенный грунт временно переходит в состояние тяжелой вязкой жидкости, что может привезти к потере устойчивости зданий и сооружений в основании которых он находится. При это очевидно, что сооружения большой массы дают большие осадки и крены, а сравнительно легкие (заглубленные подземные части зданий, резервуары и т.д.) дают подъем. В разжиженном состоянии частицы грунта отделены и практически не имеют контактов. Уплотнение происходит вследствие взаимного смещения частиц (вниз) и поровой жидкости (вверх). Критерием оценки потенциального разжижения может служить критическое ускорение, возникающее в грунте при внешних воздействиях.

Количественная оценка устойчивости грунтовых массивов при землетрясении, вибрационных и других видах динамических воздействий связано с механическими свойствами грунтов. Для описания механических свойств грунтов при динамических воздействиях применяются различные модели. Они характеризуются совокупностью уравнений, определяющих поведение среды.

В рассматриваемом процессе всякой модели соответствует некоторая схематизация свойств реальных сред, полученных из лабораторных и полевых испытаний. Критерием применимости моделей является проверка соответствия результатов расчётов напряженно-деформированного состояния (НДС) оснований с наблюдаемыми деформациями. Основное различие между используемыми моделями грунтов заключается в характере принимаемых гипотез относительно вида физических уравнений, т.е. законов деформирования и условий предельного состояния грунтов при динамических воздействиях.

 

Модель UBC3D-PLM

Модель UBC3D-PLM основана на оригинальной модели UBCSAND и представляет собой упругопластическую модель, которая позволяет моделировать поведение разжижения песчаных и глинистых грунтов под воздействием сейсмической нагрузкой.

Оценка потенциала разжижения грунтового массива может быть выполнена путем проведения анализа динамического воздействия. Для проведения такого анализа возможно применение следующего алгоритма:

  • Определение геотехнической модели грунтового массива - распределение инженерно-геологических элементов в плане и по глубине, уровней грунтовых вод, граничных условий и механических свойств грунтов для описания его поведения при статическом и динамическом нагружении.
  • Определение сейсмического воздействия в соответствии с конкретным участком и вероятностей данного воздействия, как указано в ряде действующих нормативных документов.
  • Выполнение расчетов с использованием численной модели и дальнейшего анализа результатов.

Упругое поведение

Модель UBC3D-PLM включает нелинейный изотропный закон упругого поведения, который определяется упругим модулем объемной деформации K и модулем упругого сдвига G, которые показаны в следующих уравнениях

 

Где  являются входными параметрами модели UBC3D-PLM и представляют собой модуль объемной и сдвиговой деформации соответственно;   опорное давление;   - параметры, определяющие изменение жесткости.

Абсолютно упругое поведение  определяется моделью во время процесса разгрузки. Пластическое поведение материала моделируется при достижении напряжений поверхности текучести вплоть до снижения напряжений ниже поверхности, где он снова работает упруго.

 

Пластическое поведение

Пластические деформации образуются при первичном и вторичном нагружении и догружении. Они определяются первичной и вторичной поверхностью текучести. Первичное нагружение реализуется, когда величина напряжений превышает природное напряженное состояние для нормально уплотненных грунтов и наибольшее значение напряжений для переуплотненных.

После этого материла работает в рамках первичной поверхности текучести на основе изотропного упрочнения. По мере увеличения напряжений, поверхность текучести трансформируется, расширяя область упругого поведения материала. При разгрузке материал снова начинает работать упруго. Во время повторного нагружения в случае, когда напряжения не превышают предыдущие значения - пластические деформации материала меньше по сравнению с первичной нагрузкой. Такое поведение обуславливается правилом кинематического упрочнения для вторичной поверхности текучести, которая активируется при вторичной нагрузке. Обе поверхности определяются функцией Мора-Кулона, следующим образом:

Где - максимальные и минимальные главные напряжения,   - пиковые и мобилизованные углы трения соответственно, c  – сцепление.

Функция пластического поведения основана на Drucker Prager и формулируется как:

Где q  - девиаторное напряжение;   - мобилизованный угол дилатации.

Мобилизованный угол дилатации определяется функцией от мобилизованного угла трения   и угла трения при постоянном объеме   следующим образом:

Правило упрочнения имеет вид:

Где   - модуль сдвига; np  - показатель модуля пластического сдвига; Rf  - коэффициент потери прочности;  dl  - приращение пластических сдвиговых деформаций.

При первичном нагружении , что является входным параметром. Однако при вторичной нагружении он заменяется обновленным модулем сдвига:

Где   - определяет изменение напряжений от нагружения и разгрузки;   - коэффициент от 0,5 до 1,0, корректирующий правило уплотнения,    - коэффициент уплотнения, определенный по кривой.

Как правило, при проектировании в сейсмоопасных районах, когда определение потенциала разжижения является целью расчета, методами определения параметров модели UBC3D-PLM являются лабораторные испытания на динамическом трехосном приборе или динамическом приборе простого сдвига. Однако во многих случаях доступны только данные из статических трехосных (CD TxC) или полевых испытаний (SPT).

 

Моделирование сейсмических воздействий

В рамках выполняемого исследования была составлена модель многоэтажного жилого дома (3-9-15 этажей) в программном комплексе PLAXIS. Высота этажа - 3m. Уровень подземных вод находится на уровне поверхности земли (рис. 1). Параметры грунтов основания приведены в таблице:

 

 

Таблица 1

Входные параметры, используемые в расчетах

 

Рис. 1. Расчетная схема модели

 

Запись ускорений землетрясения при магнитуде 5,40 (в открытой шкале Рихтера) и пиковым значением 239.87 cm/sec2 показана на рис. 2.

 

Рис. 2. Динамическая нагрузка: запись ускорений землетрясения

 

 

Рис. 3. Зависимость горизонтального перемещения в верхней точке здания от времени сейсмического воздействия. а. Модель Linear elastic

 

Рис. 4. Зависимость горизонтального перемещения в нижней точке здания от времени сейсмического воздействия. а. Модель Linear elastic, б. Модель UBC3D- PLM

 

 

Рис. 5. Зависимость порового давления в нижней точке здания от времени сейсмического воздействия. а. Модель Linear elastic, б. Модель UBC3D- PLM

 

В рамках представленных результатов показано, что горизонтальные перемещения в верхней точке здания не превышают 5 см (рис. 3а), а в нижней не превышают 3,5 см (рис. 4а) в линейно-упругой модели. При этом, результаты расчетов на основе модели UBC3D-PLM показали, что происходит обрушение здания (рис. 3б, 4б).

Результаты анализа расчетов по модели UBC3D-PLM  показали (рис. 5б), что грунт полностью теряет несущую способность в результате высокого порового давления, то есть происходит разжижение, вызванное сейсмической нагрузкой через несколько секунд после ее начала, что, в свою очередь, приводит к обрушению здания, что не зафиксировано при расчетах по модели Elastic linear (рис. 5а).

 

Заключение

При воздействии динамической или циклической нагрузкой модель UBC3D-PLM позволяет учесть пластические деформации при мобилизации прочности материала (сдвиговое упрочнение). Демпфирование Рэлея может быть учтено для моделирования характеристик демпфирования грунта.

Модель UBC3D-PLM разработана для прогноза динамического поведения несвязных грунтов, в особенности моделировании избыточного порового давления при недренированном поведении и разжижении.

UBC3D-PLM является достаточно простой в использовании моделью, т.к. большинство параметров грунта может быть получено в рамках полевых испытаний (SPT).

Результаты моделирования показали разницу в поведении под воздействием землетрясения в модели UBC3D-PLM по сравнению с упругим поведением в модели Linear elastic. Кроме того, результаты расчетов избыточного порового давления по модели Linear elastic в случае полностью водонасыщенных грунтов не корректны в связи с чем рекомендуется использование модели UBC3D-PLM для оценки сейсмического воздействия.

References

1. GOST R., 2015. 56353-2015 Soils. Laboratory methods for determining dynamic properties of fine-grained soil.", Standartinform, Moscow. (in Russian)

2. Jaro M. N., 2013. Interaction of two-layer base with foundations under static and dynamic loads (Doctoral dissertation, Moscow state of civil engineering). (in Russian)

3. Ivanov P. L., 1962 Liquefaction of sandy soils, Publishing House of State Energy. (in Russian)

4. Krasnikov, N. D., 1970. Dynamic properties of soils and their determination methods. L.: Strojizdat, 240. (in Russian)

5. Ter-Martirosyan A. Z., 2010. Interaction of the foundations with the base soil under cyclic and vibration impacts, considering the rheological properties of soils (Doctoral dissertation, Moscow state of civil engineering). (in Russian)

6. Ter- Martirosyan Z. G., Ter- Martirosyan A. Z., Sobolev E. S., 2014. Creep and vibro-creep of sandy soils. Engineering survey, No. 5-6, pp. 24-28. (in Russian)

7. Ter- Martirosyan Z. G., Ter- Martirosyan A. Z., Sobolev E. S., 2015. Experimental and theoretical basis for determining the mechanical properties of sandy soils under dynamic effects. Publishing House of Moscow State of Civil Engineering, Moscow. (in Russian)

8. Ukhov S. B., Semenov V. V., Znamenski V. V., Ter- Martirosyan Z. G., Chernyshev S. N., 2007. Soil mechanics-Bases and foundations, , Publishing House of " High school". (in Russian)

9. Cytovich N. A., 1963. Soil mechanics, State publishing house of literature on construction, architecture and building. materials. (in Russian)

10. ASTM D4767-11., 2011. Standard Test Method for Consolidated Undrained Triaxial Compression Test for Cohesive Soils. ASTM International, West Conshohocken, PA.

11. EN 1998-1S., 2004. Eurocode 8: Design of structures for earthquake resistance - Part 1: General rules, seismic actions and rules for buildings [Authority: The European Union Per Regulation 305/2011, Directive 98/34/EC, Directive 2004/18/EC].

12. Laera A., Brinkgreve R. B. J., 2015. Site response analysis and liquefaction evaluation. Available in the Plaxis Knowledge Base website.

13. Petalas A., Galavi V., 2013. Plaxis Liquefaction Model UBC3DPLM. PLAXIS Report.

14. Tsegaye A. 2010. Plaxis liquefaction model. external report. PLAXIS knowledge base.


Login or Create
* Forgot password?