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

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

воздействии

Л.Н. Панасюк, В. С. Тюрина, Ю.Ш. Чубка, А.У-Б. Пошев Академия строительства и архитектуры Донского государственного технического

университета, Ростов-на-Дону

Аннотация: В данной статье рассматривается механико-математическая модель деформирования и развития зон пластических деформаций в оползневых склонах при статическом и динамическом воздействии с учетом физической и геометрической нелинейности.

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

Одной из основных проблем современной строительной механики является разработка уточненных методов расчета искусственных сооружений на деформируемых основаниях [1-4]. Если для зданий и сооружений, находящихся или проектируемых на площадках со сравнительно «спокойным» рельефом зачастую используют упрощенные параметрические модели деформируемого основания, то для объектов, возводимых на площадках со сложным рельефом (особенно вблизи склонов) такие модели могут дать не только количественно, но и качественно неверные результаты. В свою очередь, использование уточненных механико-математических моделей позволит не только получить более достоверные результаты работы системы «сооружение - фундамент - деформируемое основание», но и корректно решить вопросы усиления основания. Проблема является достаточно актуальной, т.к. усиление основания в сложных геологических условиях может приводить к достаточно высоким затратам. Например, при строительстве на структурно-неустойчивых грунтах, стоимость работ нулевого цикла составляет 25-30% от общей сметной стоимости строительства.

Принципиальным в развитии расчетных схем стал учет нелинейных (геометрической, физической и конструктивной) свойств грунтов при различных динамических воздействиях [1,3]. Физическая нелинейность при описании математических моделей поведения грунтовой среды учитывается нелинейной зависимостью о~е. Геометрическая нелинейность обуславливается нелинейной зависимостью между перемещениями и деформациями, особенно проявляющейся при динамических воздействиях. Конструктивная нелинейность обусловлена тем, что в процессе деформирования грунтового основания происходит разрушение старых и возникновение новых структурных связей.

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

При решении задачи о распределении остаточных деформаций по глубине грунтового основания Л.Р.Ставницером показано [5], что характер диаграммы о~ е при динамических воздействиях является слишком сложным для того, чтобы было целесообразно аппроксимировать ее непрерывной кривой. Кроме того, свойства грунтов менее стабильны и однородны, чем для большинства искусственных сред, и поэтому нет необходимости точного подборе кривой, так как естественный разброс результатов может быть весьма значительным. В этих условиях представляется оправданным аппроксимировать зависимость моделью Прандтля с линейным упрочнением.

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

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

Известно, что процесс деформирования грунтов обуславливается воздействием всех трех инвариантов тензора напряжений. Поскольку грунт является дисперсной средой, интерес представляет собой изучение взаимного влияния первого и второго инвариантов тензора напряжений на развитие объемных деформаций. Экспериментальные исследования многих авторов (А.И.Боткина, Б.Н.Барщевского, Г.М.Ломизе, С.С.Вялова, Е.В.Виноградова, В.А.Иоселевича, Ю.К.Зарецкого, А.С. Строганованого), показали, что и формоизменение и объемные деформации грунтов зависят одновременно от величины всестороннего обжатия и девиатора напряжений. В [6] отмечается, что на основании анализа результатов экспериментов и их сопоставления в единообразной инвариантной форме сформулированы следующие общепринятые в нелинейной механике грунтов представления о деформируемости грунтовой среды при ее нагружении:

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

2. Объемная деформация зависит не только от шарового тензора напряжений, но и от девиатора напряжений.

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

Итерационный процесс Ньютона решения физически и геометрически нелинейной задачи (1) деформирования грунтового массива позволяет получить основные соотношения для формирования «секущей» и

«касательной» матриц жесткости при использовании метода конечных элементов.

Ата + р = 0 где а = О^Ли на (п+1) приближении дает:

Ат (ап + Зп Ае)+р = 0

<п = О0 п Аип

J -матрица Якоби, элементы которой вычисляются при

3 =

да

х

де

х

дт

_

де

Ь

да

ду

уг

м

Ь

дт ду

уг

Ь

' дЦ_ де 2

х

д 2и

ду де

* уг х

д 2и

де ду

х

х уг

м

Ь

д 2и

ду 2

уг

= И,

Подставив (2-4) в (1), с учетом геометрических зависимостей, получаем матричное уравнение (5).

АтН1 ААи = -р - АTD0пАип, е V

Уравнение (5) можно рассматривать как уравнение Эйлера для квадратичной аппроксимации (6) функционала Лагранжа в окрестности, достигнутого на итерации п, уровня напряженно-деформированного состояния.

При использовании теории пластического течения с упрочнением используются гипотезы, что среда изотропна, пластически сжимаема, приращение деформаций складываются из приращений упругой и пластических составляющих (7):

йе, _ йе^ + йеи

и и и

р

После ряда преобразований пластические деформации записываются как (8):

йе] _

где а1

* 3

—а

2 а

йа0 +

г У

1 * ' 3 ^

— д„а2 +----а3

V 3 У 2 2 а 3 у

да

а2 _

да

а2 _

да

а _

д^2 да

или йер _дфф йа ^дф^ _ м£р) + м£р) г

г, ^ 0 /-ч г V г, / 0 V г, / г

1 да да 1 1

г1 г1

Функции являются пластическими потенциалами. Для

появления пластических деформаций необходимо выполнение условий (9):

йа

^ йа > 0; а М- йа> 0;

да

да.

При численной реализации определяющие уравнения записывают [7,8]

так:

а _ ап +Н1 пЛе

В (10):

/-Ч л

Н _ Д +

9е

3е

АаТ _ [йап ...йт23} АеТ _ [йеп ... йу23}

]

0 0 0

0 0 0

0 0 0

0 0

0

СИММЕТРИЧНО

1 1 1 0 0 0

1 1 1 0 0 0

1 1 1 0 0 0

0 0 0 0 0 0

0 О 0 0 0 0

0 О 0 0 0 0

11

««в

«11в

аа

11т12

11 т13

«11

(?ПП в

32 ааа

Аллвнп

££ 1 _

^22^13

^ллйп

л л И^П ___ а

33^13

ааг::а

13

«1 я

1а=^а

СИММЕТРИЧНО

0 аи аа -аза

ага _ аи 0 «¡¡а — аа — е12 а1а —егэ

аЕ2 — а11 0 -а1! —а12 -аза

е12 е12 0 0 0

е13 Й1Э 0 0 0

Й2Э е21 0 0 0

Нелинейные уравнения движения по методу конечных элементов имеют вид [9-11]:

(2М + МРС + 0,5 М в2 К, ^ )дп+1 = (2М + АфС + 0,5 М в ККп )дп + (2М + М(в- 1)С + 0,5Д*2в2КК,п -0,5Д/2рККп)+

М2

^ =- ^ + 2(дп+1 - ^), 5 = Д<7

Для численного решения (11) построена абсолютно устойчивая схема (12) прямого интегрирования уравнения движения.

С 1 Л

I 1 1 Ядп1 = Я - - Ксп д„ +

V в

к=7Кём+05вК"

5.+1 =-э, + 2(4„, -Я„), 5 =Д'1

Литература

1. Работнов Ю.Н. Механика деформируемого твердого тела. М.: Наука, 1979. 650 с.

2. Lee W. Abramson, Thomas S. Lee Sunil Sharma, Glenn M. Boyce Slope stability and stabilization methods. Second Edition. New York: Wiley Pages, 2002. 736 p.

3. Ржаницын А.Р. Расчет сооружений с учетом пластических свойств материала. М.: Госстройиздат. 1954. 129 с

4. Batht K.-J. Finite Element Procedures. K.-J. Batht .New Jersey: Prentice Hall, 1996. pp. 10-12.

5. Ставницер Л.Р. Деформации оснований сооружений от ударных нагрузок. М.: Стройиздат, 1969. 350 с.

6. Бугров А.К., Нарбут Р.М., Сипидин В.П Исследование грунтов в условиях трехосного сжатия. М.: Стройиздат, 1987. 185 с.

7. Ананьев И.В. Некоторые методы решения динамических задач строительной механики: дис. ... д-р. техн наук: 05.23.17. Ростов-на-Дону, 1993. 356 с.

8. Ананьев И.В., Васильков Г.В., Хадисов М.К. Экспериментальное исследование ударного уплотнения лессовых грунтов ненарушенной структуры // Изв.Вузов. Строительство. 1992. С. 115-117.

9. Чубка Ю.Ш., Тюрина В. С., Панасюк Л.Н. Решение задач в постановке нелинейной наследственности// Инженерный вестник Дона, 2016, №3 URL: ivdon.ru/ru/magazine/ archive/n3y2016/3742

10. Тюрина В.С., Чубка Ю.Ш., Панасюк Л.Н. Моделирование работы подкрепленных оползневых склонов методом конечных элементов// Инженерный вестник Дона, 2016, №3 URL: ivdon.ru/ru/magazine/archive/n3y2016/3733

11. Панасюк Л.Н Прямые методы решения нестационарных задач теории сооружений: дис. ... д-р. техн. наук: 05.23.17. Ростов-на-Дону, 1996. 389 с.

References

1. Rabotnov Ju.N. Mehanika deformiruemogo tverdogo tela. [Mechanics of deformable solid]. M.: Nauka, 1979. 650 p.

2. Lee W. Abramson, Thomas S. Lee Sunil Sharma, Glenn M. Boyce Slope stability and stabilization methods. Second Edition. New York: Wiley Pages, 2002. 736 p.

3. Rzhanicyn A.R. Raschet sooruzhenij s uchetom plasticheskih svojstv materiala. M.: Gosstrojizdat. 1954, 129 p. [Calculation of structures with consideration of plastic properties of the material].

4. Batht K.-J. Finite Element Procedures. K.-J. Batht .New Jersey: Prentice Hall, 1996. pp. 10-12.

5. Stavnicer L.R. Deformacii osnovanij sooruzhenij ot udarnyh nagruzok. M.: Strojizdat, 1969. 350 p. [Deformation of the bases of structures from shock loads].

6. Bugrov A.K., Narbut R.M., Sipidin V.P Issledovanie gruntov v usloviyah trekhosnogo szhatiya. M.: Strojizdat, 1987. 185 p. [The investigation of soils under conditions of triaxial compression].

7. Anan'ev I.V. Nekotorye metody reshenija dinamicheskih zadach stroitel'noj mehaniki: dis. ... d-r. tehn nauk: 05.23.17. Rostov-na-Donu, 1993. 356 s. [Some methods of solving dynamic problems of building mechanics].

8. Anan'ev I.V., Vasil'kov G.V., Hadisov M.K. Eksperimental'noe issledovanie udarnogo uplotneniya lessovyh gruntov nenarushennoj struktury. Izv.Vuzov.Stroitel'stvo. 1992. pp. 115-117. [Experimental investigation of impact compaction of loess soils of undisturbed structure].

9. Chubka Ju.Sh., Tyurina V.S., Panasjuk L.N. Inzenernyj vestnik Dona (Rus), 2016, №3. URL: ivdon.ru/ru/magazine/archive/n3y2016/3742

10. Chubka Ju.Sh., Tyurina V.S., Panasjuk L.N. Inzenernyj vestnik Dona (Rus), 2016, №3. URL: ivdon.ru/ru/magazine/archive/n3y2016/3733

11. Panasjuk L.N Prjamye metody reshenija nestacionarnyh zadach teorii sooruzhenij: dis. d-r. tehn. nauk: 05.23.17. Rostov-na-Donu, 1996. 389 p. [Direct methods for solving nonstationary problems of the theory of structures].