Об одном методе решения задач реологии в лессовых просадочных грунтах

Об одном методе решения задач реологии в лессовых просадочных

грунтах

И.Ю.Дежина

Донской государственный технический университет, Ростов-на-Дону

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

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

Значительная часть территории южной России, Украины, Средней Азии занята лессовыми просадочными грунтами. Эти грунты отличаются высоким содержанием пылеватых частиц, изменением плотности скелета в пределах рй =1,3-1,6 г/см2 , пористостью до 50% и более, колебаниями

влажности от 3-4 (в засушливых районах) до 15-20 процентов, имеют характерные минералогический состав и структуру.

Замачивание лессовых грунтов под действием собственного веса или дополнительной нагрузки от сооружения может привести к возникновению значительных неравномерных деформаций - просадок. Величина относительной просадочности при этом может колебаться в пределах ^ =0,01-0,16.

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

Наибольшие затруднения возникают при проектировании зданий и сооружений на лессовых просадочных грунтах II типа по просадочности. Даже использование свайных фундаментов, прорезающих всю просадочную толщу и заглубленных в непросадочные грунты, еще не гарантируют надежной работы основания. За счет разной скорости просадки окружающего грунта и сваи по ее боковой поверхности могут развиваться силы негативного трения по величине сопоставимые с несущей способностью сваи по материалу.

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

Сокращения объема полевых испытаний можно добиться за счет использования при проектировании более достоверных математических моделей основания, повышения полноты исследований в лабораторных условиях. Это будет способствовать снижению стоимости, сроков строительства, повышению надежности эксплуатации возводимых зданий и сооружений.

Исследования показывают, что зависимость между напряжениями и деформациями при испытании просадочных грунтов носит явно нелинейный характер. Согласно литературных данных [1,2] пластические деформации являются преобладающими и могут составлять до 80-90% всей деформации. Средняя величина ползучести для лессовидных суглинков по результатам компрессионных испытаний образцов грунта при давлении 0,1-0,3 МПа для грунтов различного литологического типа составляет до 30%. В тоже время

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

Автором ведутся исследования по определению напряженно -деформированного состояния лесовых просадочных оснований при неполном водонасыщении с учетом объемного, сдвигового пластического деформирования, а также ползучести скелета и разработана методика расчета НДС оснований для лессовых просадочных грунтов как упруговязкопластических сред [3,4]. Физико-механические параметры предложенной механико-математической модели стандартны и могут быть получены в полевых или лабораторных условиях. На основе МКЭ и шагового метода разработан алгоритм расчета и программное обеспечение в составе программного комплекса "DEFOR2" [5].

Теоретические решения доведены до практической реализации.

1.Определяющее уравнение теории пластичности наследственного типа

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

Согласно лабораторным и натурным наблюдениям зависимости относительной просадки от времени Ъ1=/(1) и б51=/()} носят явно

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

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

Т.о. развитие теории вязкоупругости и теории упругости наследственного типа, сдвиговой теории упруговязкопластичности (или теории пластичности наследственного типа) является актуальной задачей определения НДС совместно с явлениями влагопереноса в лессовых просадочных грунтах.

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

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

времени появляются деформации течения, зависящие от истории нагружения.

Используем представление нелинейного функционала в виде ряда

кратных интегралов.

£ = j/i

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

Обозначим наследственное напряжение

<Г =\1

Где L - ядро ползучести. Зависимость ( 1) преобразовывается в вид

£ - ai ~ + a2~ + аз~ +

Обращение (2) приводит к виду

~ -

Уравнение (3) - основной закон одномерной нелинейной теории наследственности [4]. Здесь (p(s) - функция, определяющая мгновенную диаграмму деформирования.

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

Для приближенных решений сформулированной краевой задачи использован метод линеаризации, представление потенциальной энергии в виде ряда Тейлора, обобщенный постулат Druckera [6].

0

2. Реологические модели упруговязкопластических сред

В основе определяющих уравнений - теория течения наследственного типа с учетом объемного, сдвигового пластического деформирования,

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

= йоеу - йоу , йоу = ду + ds<¡j

Где - приращение напряжений, которое дополняет истинное

приращение с1&у до величины приращения напряжений линейно-упругой задачи.

Обозначим по тексту знаком ~ наследственные средние объемные деформации и напряжения, знаком А - наследственные напряжения и деформации формоизменения.

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

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

dGii = 3Ко ду ds о + 20о dë

(

2еУ

\

^ дуР" 2 + —Рз 3 3si

dSi

1 2ёУ

- 8уР, +— Р" 2 3 3si

dsо'

Р1 = 9(ко-Кк) . Р3 = 3 (о, - 6к) Р'2 = -Р"2 = -~2

У

При наличии аппроксимации экспериментальных данных в виде

функции со = со (¿~о ) и а = (£о^1) коэффиЦиенты & могут быть

представлены через физические характеристики среды

дао л _ да г А _?да0 ~ - да г

к к = а Ск = таг, £ = 3 дО0, ~2

3део 3д£г де г 3део

Где К0, О0 - начальные модули объемной и сдвиговой деформации, кк, Ск -наследственные касательные модули объемной и сдвиговой деформации; ^ , ~2 - наследственные модули дилатансии; йект(йа }) - приращение тензора деформаций(напяжений); е г - интенсивность деформаций; е} -компоненты девиатора деформаций; §} -дельта Кронекера.

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

йОу,} + йр .г = о, еУ,

йе1} = — (йиг,} + йи},г), е

йаг} = Нукшйе кт + Й}ктйёкш, е К^, йО}П} - ы = 0, е

йиг - duгs = 0, S2,

Уравнения связи между наследственными и истинными деформациями и экспериментальные зависимости среды

йе } = й[(^е], йёу = й[(1-^2)ег}],

аг=аг

Где

Ьк ~Оо)

Н - Бо + 3К к - К о) ' А + — • Б'2 + — Б" 2 +-—- • Б

9£г 9£г 3£2

~ _ дао л _да, , _3дао ~ _да, Кк - ' Ьк - , - , ё2 - ~ 3д£о 3д£1 д£1 д£о

При решении задач шаговым методом уравнение (4) линеаризуется.

3. Конечноэлементные алгоритмы решения задач упруговязкопластичности

Применяем итерационные методы и МКЭ. На произвольном конечном элементе вектор-функция модифицированных перемещений

и -рГ4 , еуэл ,

Где (р - матрица координатных функций, q - вектор узловых модифицированных перемещений.

Модифицированные деформации

£ - А •

После подстановки (5) - (6) в вариационное уравнение Лагранжа имеем

5П\-у£ •айУ -\у5иГрсУ -\8&Г,

§4Г

(а" + Н• ФЬ£Т)йУ -УррйУ -^

- о ,

Из(7)следует

рп А "+1 " ÍO\

кк -Р - г",

Где наследственная касательная матрица жесткости

кк ФГН" Ф йУ ,

у

вектор узловых сил

р = \фр ёУ + \р gsdS

V s

вектор узловых сил, обусловленных начальным напряжением

г" = |ф"<пёУ ,

V

Для системы конечных элементов уравнения имеют вид

К"-АС1 = р"+1 - Я",

Где Ккп,Р,Я - общие матрицы жесткости системы конечных элементов.

При организации итерационного процесса при использовании шагового метода зависимость <-£1 представляется в виде изохронных кривых. Временной интервал ^=т-А . Начальные модули объемной, сдвиговой деформации и модули дилатансии принимаются равные конечным модулям предыдущего шага. Внутри временного интервала А1 расчет ведется на полную нагрузку Р и gs или интервал АР и Аgs .

Основное соотношение (8 ) представим в виде

к о ?о + кт-1 -(1-Я*)-Адт = Арт, еуэл

Интеграл Я* -Ацт в соответствии с методом квадратурных формул(формулы трапеции ):

Я*адт = |Ад(в) ёв = Я (тА! )Ад0 — + я\т-п)А( ] -

о 2 п=1

Ад" А!+Я(о)Адт ^ ,

Тогда вязкопластическая составляющая

к т-1Адт

1-Я(о) А!

-_Арт+—Я (тА! -1 Адо +

А! кГ1 Ъ Я\т-п)А! ] -Ад"

п=1

Окончательно для конечного элемента

л , 1 т—1

кт—1Адт=Арт+ ит—1'п Ад1 + 2— кт—1т—п

п=1

Алгоритм решения следующий:

1.Вычисляется условно-мгновенные компоненты НДС для начального момента 1=0. Решается линейно-упругая задача.

К 0 Ад'0 = АР'0

Где начальная матрица жесткости

К 0 = | ФTD0ФdV

V

Определяются начальные деформации и напряжения

£ = Ф д , и = DоS\

2.Решается упругопластическая задача

КП Ад'п = АР'п

З.Организуется итерационный процесс упруговязкопластической задачи при АР=0.

т—1

кт—10Адт=Арт+ кт—1,п Ад1 + 2 кт—1>т—п - Адп,

п=1

После решения системы уравнений с использованием (9) относительно

л т

Ад определяем

т т—1, * т

д = д + Ад

Значения модифицированных перемещений и деформаций организуем с использованием формул

лт - {л и* {л 0* \ т

д = Ц—ЯтЯ > £ =11—Ят£

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

использования выражения (1-Я*) - определение методом выравнивания кривых из экспериментальных данных.

Расчет повторяют до тех пор, пока КтА1. Число итераций т назначают таким образом, чтобы достигнуть стабилизации решения.

Для использования методики в практических целях необходимо проводить серию экспериментов в лабораторных условиях на компрессионных, сдвиговых, стабилометрических приборах. Определение объемных и сдвиговых модулей деформации с учетом фактора времени сводится к получению математической зависимости деформаций от времени и определению вида и параметров ядер ползучести и реалаксации при различных значениях влажности и уплотняющих давлениях [7,8]. Использование достоверных математических моделей основания с учетом деформаций ползучести при проектировании на лессовых просадочных грунтах в условиях неполного водонасыщения способствует повышению надежности эксплуатации возводимых строительных объектов [9,10,11,12].

Выводы

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

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

3.Алгоритм решения упруговязкопластических задач построен на использовании итерационных процессов. Вначале определяются условные

мгновенные компоненты НДС для начального временного интервала. При появлении пластических деформаций решается упруго вязкопластическая задача.

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

Литература

1. Вялов С.С. Реологические основы механики грунтов. 1978 изд. М.: Высш.школа, 1978. 447 с.

2. Зарецкий Ю.К. Вязкопластичность грунтов и расчеты сооружений. М.: Стройиздат, 1988. 352 с.

3. Васильков Г.В., Дежина И.Ю., Орта-Ранхель А. О реологических моделях упруговязкопластических сред // Изв.вузов Сев-Кавк.регион. Естественные науки. 1998. №3. С. 21-36.

4. Васильков Г.В., Дежина И.Ю., Орта Ранхель. Механико-математическая модель пористой упруговязкопластической среды при неполном влагонасыщении // Труды VII школы-семинара "Современные проблемы механики грунтов и охраны геологической среды. Ростов-на-Дону: изд.Рост.госуд.ун-та и НИИ механики и прикладной математики, 1998. С. 12-13.

5. Свидетельство о государственной регистрации программы для ЭВМ N2015619078. Определение напряженно-деформированного состояния оснований и фундаментов по деформациям с использованием нелинейных моделей DEFOR2.

6. Работнов Ю.Н. Элементы наследственной механики твердых тел. М.: Наука, 1977. 384 с.

7. Эюбов А.А. Ползучесть лессовых просадочных грунтов при компрессии // Основания, фундаменты и механика грунтов. 1986. №3. С. 29-31.

8. Месчан С.Р. Реологические процессы в глинистых грунтах. Ереван: Айастан, 1992. 392. с.

9. Чмшкян А.В. Совершенствование методов расчета просадочных деформаций // Инженерный вестник Дона, 2012, №4 (часть2) URL: ivdon.ru/magazine/achive/n4p2y2012/1256.

10. Филипов А.М. Аналитическое численное решение плоской контактной задачи с учетом ползучести материала// Инженерный вестник Дона, 2015, №1 URL: ivdon.ru/magazine/achive/nly2016/3510 .

11. Gioda, G. and A. Cividini, 1996. Numerical Methods for the Analysis of Tunnel Performance in squeezing Rocks. Rock Ме^ап^ and Rock Engineering (issue 29(4)), Austria:Springer-Veriag, pp: 171-193.

12. Yin, Z.Y., C.S. Chang, M. Karstunen and P.Y. Hicher, 2010. An anisotropic elasticviscoplastic model for soft clays. International Journal of Solids and Structures, 47, pp: 665-677.

References

1. Vyalov S.S. Reologicheskie osnovy mekhaniki gruntov [Rheological Bases of the Soil Mechanics]. 1978 izd. M.: Vyssh.shkola, 1978. 447 p.

2. Zaretskiy Yu.K. Vyazkoplastichnost' gruntov i raschety sooruzheniy [Visco-plastics of Soils and Calculations of Constructions]. M.: Stroyizdat, 1988. 352 p.

3. Vasilkov G.V., Dezhina I.Y., Orta-Rankhel А. Northern-Caucasus region. Estestvennie Nauki. 1998. №3. pp. 21-36.

4. Vasil'kov G.V., Dezhina I.Yu., Orta Rankhel'. Trudy VII shkoly-seminara "Sovremennye problemy mekhaniki gruntov i okhrany geologicheskoy sredy. Rostov-na-Donu: izd.Rost.gosud.un-ta i NII mekhaniki i prikladnoy matematiki, 1998. pp. 12-13.

5. Svidetel'stvo o gosudarstvennoy registratsii programmy dlya EVM N2015619078. Opredelenie napryazhenno-deformirovannogo sostoyaniya osnovaniy i fundamentov po deformatsiyam s ispol'zovaniem nelineynykh modeley DEFOR2 [Determination of Tense and Deformed State of Grounds and Foundations According to Deformations Using the Non-Linear Models DEFOR2]. 6 Rabotnov Yu.N. Elementy nasledstvennoy mekhaniki tverdykh tel. [Elements of Hereditary Mechanics of Solid Bodies]. M.: Nauka, 1977. 384p.

7. Eyubov A.A. Polzuchest'. Osnovaniya, fundamenty i mekhanika gruntov. 1986. №3. pp. 29-31.

8. Meschan S.R. Reologicheskie protsessy v glinistykh gruntakh [Rheological Processes in Clay Soils]. Erevan: Ayastan, 1992. 392. p.

9. Chmshkyan A.V. Inzenernyj vestnik Dona (Rus), 2012, №4 (part 2) URL: ivdon.ru/magazine/achive/n4p2y2012/1256.

10. Filipov A.M. Inzenernyj vestnik Dona (Rus), 2015, №1. URL: ivdon.ru/magazine/achive/nly2016/3510.

11. Gioda, G. and A. Cividini, 1996. Rock Mechanics and Rock Engineering (issue 29(4)), Austria:Springer-Veriag, pp. 171-193.

12. Yin, Z.Y., C.S. Chang, M. Karstunen and P.Y. Hicher, 2010. International Journal of Solids and Structures, 47, pp. 665-677.