Оценка влияния выбора модели на результат недренированного расчета | Статья в журнале «Молодой ученый»

Отправьте статью сегодня! Журнал выйдет 6 апреля, печатный экземпляр отправим 10 апреля.

Опубликовать статью в журнале

Автор:

Рубрика: Архитектура, дизайн и строительство

Опубликовано в Молодой учёный №22 (469) июнь 2023 г.

Дата публикации: 01.06.2023

Статья просмотрена: 7 раз

Библиографическое описание:

Свербаев, А. Н. Оценка влияния выбора модели на результат недренированного расчета / А. Н. Свербаев. — Текст : непосредственный // Молодой ученый. — 2023. — № 22 (469). — С. 106-111. — URL: https://moluch.ru/archive/469/103512/ (дата обращения: 29.03.2024).



Для описания механического поведения грунта в современной инженерной практике используются сложные математические модели. При этом, их создание базируется, в основном, на результатах лабораторных испытаний (компрессионных, сдвиговых, трехосных и т. д.): конечным итогом считается набор завершающих уравнений состояния для единичного объема описываемого массива. При этом, сами уравнения, равно как и порядок определения компонент тензоров силовых и кинетических характеристик, обычно оставляется без должного внимания. В работе указаны следствия различий алгоритмов вычисления искомых величин наиболее популярных моделей грунта на 2-х характерных расчетных схемах. Автором приведены блок-схемы последовательности моделирования, подтвержденные аналитическим расчетом. Также даны рекомендации по выбору моделей для расчета в недренированной постановке.

Ключевые слова : нелинейные модели грунта, недренированное поведение.

In modern engineering practice, complex mathematical models are used to describe the mechanical behavior of soil. However, their creation is based mainly on the results of laboratory tests (compression, shear, triaxial, etc.): the result is considered to be a set of final equations of state for a unit volume of the described body. In this case, the equations themselves, as well as the order of determining the components of the force and kinetic characteristics tensors, are usually left without due attention. The paper points out the consequences of the differences in the algorithms for calculating the required values of the commonly used soil models on two most popular calculation schemes. The author gives block diagrams of the modeling sequence, confirmed by analytical calculations. In addition, there are recommendations on the choice of models for calculations in the undrained formulation.

Keywords : non-linear soil models, undrained behavior.

1. Введение

Для расчета оснований, сложенных слабыми связными грунтами, предлагается использовать разные модели грунта. Наиболее изученными и распространёнными являются модели, семейств Сam Clay (и её модификации, например, Soft Soil, Soft Soil Creep), Hardening Soil (и её модификации, например, Hardening Soil Small, Generalized Hardening Model). Однако не только данные модели позиционируются как универсально-нелинейные, пригодные для расчета слабых глин. В частности, в программном комплексе Plaxis, представлены еще несколько моделей (NGI-ADP, Sekiguchi-Ohta и т. д.)

Такое разнообразие, впрочем, приводит к закономерному вопросу о выборе наиболее качественной, равно как и о влиянии выбора модели на конечный результат расчета.

2. Методы

Для детального изучения моделей использована виртуальная среда Plaxis Soiltest, а также численный анализ Plaxis для простых, но в то же время принципиальных расчетных схем (рис. 1)

Общий вид расчётных схем

Рис. 1. Общий вид расчётных схем

Три основных модели грунта, рассматриваемые в данной работе, наглядно демонстрируют основные фрэймворки (англ. framework — основа, базис, имеется в виду теоретическая экспликация пластического поведения грунта) которые используются для обоснования методов геотехнических расчетов [1, 2, 4]:

— Модель слабого грунта Soft soil — эволюция шатровых моделей критического состояния типа Cam Clay

— Модель Норвежского геотехнического института NGI-ADP — опирается на эмпирический метод SHANSEP и линейную связь между индексами сжатия и пластичности

— Модель упрочняющегося грунта Hardening Soil — как и предыдущая модель, использует коэффициент переуплотнения (OCR) и, в добавок, апеллирует к теории устойчивого состояния, предполагая связь между скоростями деформаций и приращением напряжений по ветви первичного нагружения

— Также для получения контрольной группы значений, произведены расчеты с использованием идеально-упругопластической модели Кулона-Мора.

Одной из самых значительных проблем прикладной геотехники на момент написания статьи можно назвать зачаточное состояние методик интерпретации лабораторных испытаний в рамках предлагаемых математических моделей, в купе с отсутствием унифицированной системы обозначений используемых параметров [3, 6, 7, 8].

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

Исходными данными для работы послужили результаты трёхосных испытаний слабого связного грунта по КД схеме, выполненных в соответствии с ГОСТ 12248.3–2020.

Далее, с помощью утилиты Parameter Optimisation в программе Plaxis Soiltest, были собраны пакеты исходных данных, гарантирующие максимальное соответствие модели материала реальным испытаниям (рис. 2). Исходные данные представлены в таблице 1.

Таблица 1

Hardening Soil

Soft Soil

NGI - ADP

Mohr - Coulomb

R inter

0,67

-

e

0,38

-

γ

15,8

кН/м 3

γ unsat

16,0

кН/м 3

E 50

9612

кПа

λ*

0,038

-

G ur /s UA

56

-

E’

5240

кПа

E oed

5240

кПа

κ*

0,015

-

γ fC

10

%

E ur

50417

кПа

M

0,8968

-

γ fE

28

%

m

0,7723

-

γ fDSS

16,9

%

c’

29,7

кПа

c’

29,5

кПа

s UA

72

кПа

c’

29,4

кПа

ϕ

14,63

°

ϕ

14,23

°

s CTX /s UA

0,99

-

ϕ

14,23

°

ψ

0

°

ψ

0

°

s UP /s UA

0,5

-

ψ

0

°

R f

0,7811

-

τ 0 /s UA

0,01

-

K 0

0,7543

-

K 0

s UDSS /s UA

0,8

-

ν ur

0,2146

-

ν

0,3

-

ν

0,32

-

ν’

0,3

Графики результатов реальных и виртуальных трехосных испытаний (КД)

Рис. 2. Графики результатов реальных и виртуальных трехосных испытаний (КД)

Исходя из предположения об относительной несжимаемости воды и гипотезы Терцаги [4, 5], очевидно, что наличие избыточного порового давления для РС 1 вызовет снижение устойчивости основания (рис. 3а), а, следовательно, и меньший коэффициент безопасности (M safety ). И, напротив, для РС 2 избыточное поровое давление будет «догружать» массив основания, уменьшая тем самым нагрузку на конструкцию ограждения (рис. 3б).

Принципиально ожидаемый вид для разных условий расчета

Рис. 3. Принципиально ожидаемый вид для разных условий расчета

а) графика «нагрузка — осадка» основания насыпи б)графиков распределения изгибающих моментов в ограждающей конструкции по глубине

3. Результаты и обсуждения

Характерные результаты расчетов приведены в таблицах 2, 3

Таблица 2

Коэффициенты устойчивости

Условия\Модель

Hardening Soil

Soft Soil

NGI-ADP

Mohr-Coulomb

P = 100 кПа

Undrained

1,545

1,331

1,414

1,792

Drained

2,470

2,465

-

2,471

P = 150 кПа

Undrained

1,126

1,348

1,073

1,250

Drained

1,832

1,822

-

1,826

Таблица 3

Максимальные абсолютные значения изгибающих моментов, кНм

Условия\Модель

Hardening Soil

Soft Soil

NGI - ADP

Mohr - Coulomb

H = 6 м; H plate = 15 м

Undrained

0,21

- 53,15

20,88

- 0,06

11,87

-8,14

30,74

- 0,24

Drained

0,87

- 27,04

32,77

- 7,77

-

26,31

- 41,73

Заметим имеющееся совпадение результатов дренированного расчета насыпи и разброс, как результатов недренированного расчета для РС-1, так результатов для РС-2, в том числе и по виду самой эпюры изгибающих моментов (рис. 4).

Различия в эпюрах изгиб. моментов, в зависимости от модели и учета/игнорирования избыточного порового давления

Рис. 4. Различия в эпюрах изгиб. моментов, в зависимости от модели и учета/игнорирования избыточного порового давления

Основная причина наблюдаемого — во внутренних связях компонент деформаций и напряжений и порядке их получения [6, 7].

Для идеально-упругопластической модели Кулона-Мора фаза неупругих деформаций вырождается в точку достижения критерия разрушения. До этого момента связь «напряжения-деформации» линейна. Это объясняет, почему коэффициенты устойчивости у неё относительно высоки, а наличие порового давления снижает устойчивость насыпи.

Модель слабого грунта показывает завышенные значения для всех условий расчета. Это связано с ошибкой, заложенной в самой модели: приоритет отдан объемным деформациям. Однако, в условиях отсутствия (или чрезвычайно медленной) фильтрации объемные деформации относительно невелики, что ведёт к мнимой минимизации перемещений (рис. 5)

Алгоритмы вычисления компонент моделей грунта для моделей Soft soil (а) и Hardening soil (б)

Рис. 5. Алгоритмы вычисления компонент моделей грунта для моделей Soft soil (а) и Hardening soil (б)

Теоретически, модель двойного упрочнения должна иметь преимущество пред Soft Soil, т. к. обладает дополнительными независимыми уравнениями состояния. Однако наследие шатровых моделей в виде ассоциированного закона течения приводит к вычислению пластических деформаций сдвига в том числе и на основе полученных пластических объемных деформаций, определяющихся из шатра-аналога (изолиний равных объемных деформаций) моделей типа «Cam Clay».

Подтверждение приведенных алгоритмов проведено в среде автоматизированного проектирования MathCAD. Сравнение аналитических графиков с результатами Plaxis Soiltest приведено на рисунке 6.

Сходимость аналитического и программного графиков e-q для Soft soil (а) и Hardening soil (б)

Рис. 6. Сходимость аналитического и программного графиков e-q для Soft soil (а) и Hardening soil (б)

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

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

4. Выводы

  • Применение современных моделей ограничено условиями получения для них входных параметров. В частности, в описании и технической документации по моделям грунта отсутствуют методики получения параметров на основе упоминающийся лабораторных испытаний.
  • Разница в последовательности вычисления компонент тензоров напряжений и деформаций даёт глобальную разницу в результатах расчёта. Алгоритм вычисления должен иметь физическое обоснование, а не математическое удобство. Сходимость результатов моделирования лабораторных испытаний не гарантирует точности при анализе конкретной расчетной схемы.
  • Для более консервативных результатов рекомендуется модели:

Hardening soil для н/дрен. расчета котлована

Soft soil или Hardening soil для дрен. расчета котлована

NGI - ADP или Soft soil — для н/дрен. расчета насыпи

Mohr - Coulomb — для дрен. расчета насыпи

5. Заключение

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

Литература:

  1. Шашкин А. Г. Критический анализ наиболее распространенных нелинейных моделей работы грунта // Инженерная геология. М.: 010. № 3. с. 9–37.
  2. Мирный А. Ю. Области применения современных механических моделей грунтов / А. Ю. Мирный, А. З. Тер-Мартиросян // Геотехника. — 2017.: No 1. — С. 20–26
  3. Lade, P. V. Overview of constitutive models for soils. // Soil Constitutive Models: Evaluation, Selection and Calibration, pp 1–34. 2005
  4. Yin, Z. Y. Practice of Constitutive Modelling for Saturated Soils / Z. Y. Yin, P. Y. Hicher, Y. F. Jin. — 1. — Singapore: Springer, 2020. — 400 c.
  5. Hai-Sui, Y. Plasticity and Geotechnics / Yu Hai-Sui. — 13. NY: Springer, 2006. — 522 c.
  6. Plaxis Material Models Manual 2D
  7. Свербаев, А. Н. Основные математические компоненты моделей грунта, сформулированных в рамках механики сплошной среды / А. Н. Свербаев. —// Серия «Строительство»: сборник статей магистрантов и аспирантов. Выпуск 6. В 2-х т. Т. 1 / Санкт-Петербургский государственный архитектурно-строительный университет. — Санкт-Петербург: СПбГАСУ, 2023. — 524 с. — С. 480–497.
  8. Строкова, Л. А. Научно-методические основы численного прогноза деформирования грунтовых оснований: специальность 25.00.08: диссертация на соискание ученой степени доктора геолого-минералогических наук / Строкова Людмила Александровна; ФГБОУ ВПО «Национальный исследовательский томский политехнический университет». — Томск, 2011. — 265 c. с ил.
  9. Практическая проблема использования прочности по теории плотности-влажности вместо недренированного сопротивления сдвигу. — Текст: электронный // База знаний PLAXIS: [сайт]. — URL: https://plaxisrussia.notion.site/aefab79ebc8f4d4ea5129c5a98a0322c (дата обращения: 31.05.2023).
Основные термины (генерируются автоматически): NGI-ADP, модель, ADP, NGI, избыточное поровое давление, модель грунта, CTX, OCR, идеально-упругопластическая модель Кулона-Мора, слабый грунт.


Ключевые слова

нелинейные модели грунта, недренированное поведение

Похожие статьи

Основные составляющие геомеханической модели резервуара

В данной статье приведены основные оперируемые величины, составляющие геомеханическую модель.

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

Модифицированные критерии Писаренко — Лебедева и Кулона...

Рис. 2. Глубинная колея, обусловленная пластическим деформированием грунтов земляного полотна и дискретных

– для расчета на сдвиг трехпараметрический критерий КулонаМора [4].

11. Александрова Н. П., Семенова Т. В., Долгих Г. В. Совершенствование моделей расчета главных

Долгих Г. В. Расчет грунтов земляного полотна по критерию безопасных давлений.

Прогнозирование процесса разработки с использованием...

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

Граничное состояние между упругой и неупругой деформацией определяется методом Кулона-Мора.

Упругопластическая модель учитывает и упругие деформации.

Модель деформирования впроцессе падения пластового давления.

Выбор варьируемых параметров при адаптации... | Молодой ученый

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

Математическая модель процесса топливоподачи системой...

Библиографическое описание: Субботенко, Д. И. Математическая модель процесса топливоподачи системой Common Rail с

- давление в аккумуляторе в начале каждого расчетного режима считается постоянным и равным остаточному давлению в системе

Модели блочной среды для исследования колебательных...

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

Моделирование коэффициентов активности компонентов системы...

Составы равновесных фаз для системы ацетон (1) –пропанол (2) при общем давлении Ри температуре 25º C.

Наибольшее распространение обрела модель «локальных составов». Первое уравнение, образованное на этой модели, было выведено в 1964 г. Вильсоном.

Особенности многовариантного геологического моделирования на...

Геологическое моделирование начинается с создания концептуальной геологической модели, которая служит основой модели фаций, исходя из знания региональных палеоусловий осадконакопления [3].

Рис. 2. Модель осадконакопления пласта Ю.

Моделирование ударных волн сжатия в пористых металлах

При построении модели используется реологический метод. На рис. 1 приведена реологическая схема пористой среды [1].

Модель описывает два варианта решения: - схлопывание пор происходит на этапе упругой деформации, пластичность проявляется после уплотнения среды

Адаптация гидродинамической модели месторождения N на...

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

Похожие статьи

Основные составляющие геомеханической модели резервуара

В данной статье приведены основные оперируемые величины, составляющие геомеханическую модель.

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

Модифицированные критерии Писаренко — Лебедева и Кулона...

Рис. 2. Глубинная колея, обусловленная пластическим деформированием грунтов земляного полотна и дискретных

– для расчета на сдвиг трехпараметрический критерий КулонаМора [4].

11. Александрова Н. П., Семенова Т. В., Долгих Г. В. Совершенствование моделей расчета главных

Долгих Г. В. Расчет грунтов земляного полотна по критерию безопасных давлений.

Прогнозирование процесса разработки с использованием...

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

Граничное состояние между упругой и неупругой деформацией определяется методом Кулона-Мора.

Упругопластическая модель учитывает и упругие деформации.

Модель деформирования впроцессе падения пластового давления.

Выбор варьируемых параметров при адаптации... | Молодой ученый

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

Математическая модель процесса топливоподачи системой...

Библиографическое описание: Субботенко, Д. И. Математическая модель процесса топливоподачи системой Common Rail с

- давление в аккумуляторе в начале каждого расчетного режима считается постоянным и равным остаточному давлению в системе

Модели блочной среды для исследования колебательных...

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

Моделирование коэффициентов активности компонентов системы...

Составы равновесных фаз для системы ацетон (1) –пропанол (2) при общем давлении Ри температуре 25º C.

Наибольшее распространение обрела модель «локальных составов». Первое уравнение, образованное на этой модели, было выведено в 1964 г. Вильсоном.

Особенности многовариантного геологического моделирования на...

Геологическое моделирование начинается с создания концептуальной геологической модели, которая служит основой модели фаций, исходя из знания региональных палеоусловий осадконакопления [3].

Рис. 2. Модель осадконакопления пласта Ю.

Моделирование ударных волн сжатия в пористых металлах

При построении модели используется реологический метод. На рис. 1 приведена реологическая схема пористой среды [1].

Модель описывает два варианта решения: - схлопывание пор происходит на этапе упругой деформации, пластичность проявляется после уплотнения среды

Адаптация гидродинамической модели месторождения N на...

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

Задать вопрос