Full Text

Введение Консолидация - это процесс уплотнения пористого водонасыщенного грунта под нагрузкой при возможности оттока жидкости. Грунты являются естественным природным пористым материалом, который состоит из твердого скелета с порами, заполненными жидкостью. Твердая и жидкая фазы грунта взаимодействуют друг с другом. С одной стороны, увеличение порового давления приводит к росту сжимающих напряжений в скелете грунта и к деформированию скелета. С другой стороны, сжатие скелета грунта может вызвать уменьшение пористости и, как следствие, увеличение давления флюида. Эти связанные явления совместного упругого деформирования грунта и фильтрации жидкости были впервые описаны в теории линейной пороупругости М. Био [1], в которой соединены уравнение равновесия для твердого скелета и закон линейной фильтрации Дарси. Однако процесс деформирования флюидонасыщенных дисперсных грунтов под нагрузкой, как правило, нелинейный и сопровождается пластическими деформациями. В песчаных и глинистых грунтах под действием нагрузки деформации развиваются во времени с гидродинамическим запаздыванием и зависят от скорости оттока жидкости из грунта. Реакция скелета грунта на изменение порового давления может оказаться критичной, например, при ухудшении ровности дорожного покрытия и образовании колеи или при неравномерной осадке инженерных сооружений. Следовательно, разработка строгой и эффективной нелинейной связанной модели консолидации грунта, учитывающей упругопластическое поведение материала, сопровождающееся изменением пористости и проницаемости [2], весьма актуальна и имеет большое практическое значение. Большинство разработанных моделей консолидации базируются на линейной теории упругости или упругопластичности для малых деформаций [1; 3-5]. В настоящее время в литературе представлены обобщения классических уравнений консолидации с бесконечно малых на конечные деформации, например в работах [6-9]. Обобщения постановки классических уравнений консолидации на большие деформации в основном базируются на использовании определяющих соотношений в скоростях [10]. Альтернативная формулировка для упругопластичности с конечными деформациями базируется на мультипликативном разложении градиента деформации [6-9]. Задачи совместного деформирования грунта и движения жидкости могут решаться в несвязанной постановке (упругий режим фильтрации) и в связанной постановке. При решении задач об откачке нефти из глубоких скважин упругий режим фильтрации в отдельных случаях бывает вполне достаточен. Задачи о деформировании водонасыщенного грунта под нагрузкой (задачи консолидации), как правило, решаются в связанной постановке. Наиболее часто для решения подобных задач используется метод конечных элементов (МКЭ) [6-9, 11-13]. В последнее время взаимодействие между течением жидкости и скелетом грунта моделируется с использованием различных связанных схем. Уравнение фильтрации и уравнение равновесия среды могут решаться на каждом шаге по времени в виде общей системы [14; 15]. В этих работах метод решения называется полностью связанным. В другом варианте метода решения уравнение фильтрации и уравнение равновесия связываются последовательно в итерационном процессе типа метода Узавы. Сначала решается или задача фильтрации, или механическая задача, а потом решается другая задача, используя результат решения первой [12-17]. Если заданы достаточно жесткие критерии сходимости, итерационно связанное решение может быть таким же точным, как и полностью связанное, но количество требуемых итераций может быть чрезмерно большим [5]. Итерационно связанные последовательные методы могут использовать независимые решатели механики сплошных сред и средства моделирования пластового потока, развитые в нефтяной промышленности [17]. Численное решение уравнений трехмерной консолидации Био по-прежнему остается сложной задачей, несмотря на интенсивные исследования в этой области [5]. Во-первых, полностью связанный подход приводит к большим алгебраическим системам, которые могут быть довольно сложными для решения. При этом в реальных трехмерных приложениях количество неизвестных может достигать нескольких сотен тысяч. Во-вторых, матрица коэффициентов, полученная в результате численной дискретизации, может быть плохо обусловлена. Следовательно, для работы с полностью связанной схемой необходимы продвинутые решатели [18; 19]. При использовании итерационно связанной схемы сходимость итерационного процесса, зависящая от механических свойств жидкости и грунта, может быть медленной [16]. Третья проблема связана с различными формами неустойчивости при численном решении, которые проявляются в виде ложных осцилляций порового давления. Это может происходить, если сжимаемость жидкости достаточно мала по сравнению с объемной сжимаемостью скелета грунта [5]. Как результат, на начальном этапе скелет грунта ведет себя как несжимаемая среда. Тогда вначале скорость изменения объема скелета грунта близка нулю, что приводит к патологии конечных элементов (КЭ), известной как «запирание» (locking), которое обычно вызывает ложные осцилляции порового давления [5]. Справиться с этой проблемой можно, используя аппроксимации разного порядка для перемещений и порового давления. Известно, что поровое давление следует аппроксимировать линейными конечными элементами, а перемещения - квадратичными КЭ (см., например, [20]). В нашей работе формулируется физически и геометрически нелинейная связанная модель консолидации в скоростях, для решения которой используется итерационный метод Узавы. Для пространственной дискретизации седловой системы уравнений используются квадратичные серендиповы элементы для аппроксимации приращений перемещений и элементы трилинейного типа для аппроксимации изменения давления жидкости. Известно, что при решении связанных динамических задач консолидации различные исследователи часто используют метод ALE (Arbitrary Lagrangian - Eulerian) для исправления сильно искаженной лагранжевой сетки [21; 22], что приводит к перестроению сетки и, соответственно, к некоторой потере точности решения. Следует отметить, что в нашей работе метод ALE используется только для объединения Эйлерова и Лагранжева подходов и не связан с перестроением лагранжевой сетки. 1. Математическая постановка нелинейной связанной задачи консолидации и метод решения Дифференциальная постановка связанной задачи консолидации «в скоростях» включает три уравнения - равновесия, фильтрации и пористости [23-25]: (1) Система дифференциальных уравнений (1) сформулирована относительно дифференциалов вектора перемещений давления жидкости и пористости . Здесь - лагранжевы координаты твердого скелета, и - касательные модули соответственно в текущей и начальной конфигурациях, F - деформационный градиент, , - второй тензор эффективных напряжений Пиола - Кирхгофа, E - тензор деформаций Грина - Лагранжа, f - массовая сила, - эффективные напряжения Коши (осредненные напряжения при нулевом поровом давлении), - средняя плотность пористого водонасыщенного грунта, - тензор коэффициентов фильтрации, - параметр Био, - удельный вес жидкости, - сжимаемость жидкости, - пористость. Уравнения системы (1) выведены из закона фильтрации Дарси, уравнения состояния жидкости, уравнений неразрывности и уравнения равновесия [23; 25]. При выводе были сделаны следующие предположения и допущения. Процесс деформирования грунта изотермический. Поры грунта полностью заполнены флюидом. Жидкость ньютоновская и однородная, изменение плотности жидкости при изменении порового давления подчиняется закону баротропии . Удельный вес жидкости принимается постоянным при конкретной температуре и минерализации. Пористость зависит от перемещений в грунте: коэффициенты фильтрации зависят от пористости: В уравнениях неразрывности материал скелета грунта принимается несжимаемым по сравнению с сжимаемостью жидкости, так как плотность минеральных зерен rs, в отличие от плотности жидкости rf, слабо зависит от изменения давления [26]. Несжимаемость материала скелета не исключает объемную деформацию осредненного каркаса за счет перемещения частиц твердой фазы грунта или переупаковки зерен в дисперсных несвязных грунтах. Отметим, что в уравнении равновесия сжимаемость материала скелета учитывается. Уравнения фильтрации и пористости в (1) были переформулированы в движущейся системе координат, связанной с лагранжевыми координатами x твердого скелета согласно методу ALE (Arbitrary Lagrangian-Eulerian) [27] с использованием относительной скорости движения жидкости (скорости флюида относительно скелета грунта) [23; 25]. В задаче консолидации задавались следующие граничные условия. Часть границы неподвижна и непроницаема, остальная часть границы проницаема, и на ней действует поверхностная нагрузка: , (2) где n - внешняя нормаль к границе области , S - поверхностная сила. Для дальнейшего решения задачи консолидации численными методами была получена вариационная постановка в текущей конфигурации, а уравнение равновесия было линеаризовано с использованием техники дифференцирования по Гато [23; 25]. Системе (1) с граничными условиями (2) соответствует линеаризованная вариационная постановка задачи консолидации [25]: (3) Здесь S - поверхностная сила из граничных условий (2), w, q, h - пробные функции из подпространств Соболева, удовлетворяющих заданным граничным условиям 1-го рода: В качестве определяющих соотношений в (3) могут использоваться соотношения для гиперупругого материала при существовании упругого потенциала, а также теория пластического течения или деформационная теория пластичности. Сначала решаются первые два уравнения (равновесия и фильтрации) системы (3) в предположении постоянной пористости. Для решения используется обобщение неявной схемы, предложенной еще в [28; 29], с внутренними итерациями на каждом шаге по времени по методу Узавы [30]. Для реализации итерационного процесса на каждом шаге по времени применялся метод простой итерации с использованием предобусловливателя для ускорения сходимости. После решения этой задачи вычисляется пористость из третьего уравнения системы (3). Далее посчитанное значение пористости учитывается при решении уравнений равновесия и фильтрации на следующем временном шаге. Дискретизация линеаризованных вариационных уравнений равновесия и фильтрации по пространственным координатам осуществлялась с помощью метода конечных элементов (МКЭ). При дискретизации использовались 20-узловые серендиповы квадратичные элементы (Q2) для аппроксимации приращения перемещений скелета грунта и 8-узловые элементы трилинейного типа (Q1) для аппроксимации приращения давления воды [24; 25]. Такой выбор конечных элементов был обусловлен тем, что система уравнений равновесия и фильтрации является седловой, а при численном решении дискретизированной седловой системы возникает проблема корректности. Достаточным условием корректности этой системы является LBB-условие [31], выполнение которого зависит от выбранных при дискретизации конечных элементов для аппроксимации и . Как показывает опыт моделирования седловых задач [7; 32], аппроксимация квадратичными конечными элементами Q2 перемещений скелета грунта и трилинейными КЭ Q1 давления жидкости обеспечивает корректность численного решения задачи. А выполнение LBB-условия для элементов Q1 - Q1 зависит от соотношения шага по времени и шага сетки h: численное решение седловой задачи корректно при [33; 34]. Уменьшение и, как следствие, нарушение LBB-условия приводит к возникновению ложных осцилляций порового давления, которые появляются на некотором временном шаге и далее усиливаются с каждым шагом по времени, в результате чего решить седловую систему становится невозможно [24; 25]. 2. Пример численного решения задачи об упругопластическом деформировании водонасыщенного грунта Рассмотрим пример численного решения задачи о деформировании водонасыщенного песчаного грунта под действием перфорированного штампа. Для моделирования упругопластического деформирования грунта при кратковременных нагрузках предложено обобщение модели С.С. Григоряна [35] на большие деформации. Для формулировки соотношений теории пластичности использовался тензор Кирхгофа , как предлагается в работах [36-38]. Функция нагружения имеет вид: (4) Здесь - девиатор тензора напряжений, - шаровая часть тензора напряжений. Экспериментальные исследования песчаных грунтов [39] подтвердили справедливость условия пластичности f (4), а также позволили установить, что для песчаных грунтов в реальных условиях при функция F имеет вид где A и M - константы. В работе [35] для сдвиговых свойств выбрана модель идеальной пластичности, но в численных расчетах удобно использовать малое изотропное упрочнение для устойчивости вычислений, как показано в формуле (4). В критерии Кулона - Мора, который лучше всего описывает свойства слабо связных грунтов, параметры A и M из (4) выражаются формулами [40]: где - угол внутреннего трения грунта, c - сцепление. Рис. 1. Зависимость объемных напряжений и деформаций при упругопластическом деформировании грунта Fig. 1. Dependence of volumetric stresses and strains during elastoplastic deformation of the soil Сделаем допущение, что скорость объемной деформации близка к постоянной: Для этой скорости можно принять при нагрузке . Диаграмма для шаровых частей (рис. 1) также соответствует конкретной скорости нагружения. В настоящий момент в расчетах упругопластического деформирования грунта используются определяющие соотношения деформационного типа (5) в рамках предположения, что для деформаций деформационная теория сохраняет приемлемую точность, за исключением областей вблизи угловых точек, где имеют место большие деформации и сложное нагружение. (5) где - функция Ильюшина, - компоненты девиатора накопленной деформации Грина - Лагранжа, - интенсивность накопленной деформации, - мера деформаций Коши - Грина, - касательный модуль на диаграмме - предел текучести при - упругий модуль сдвига, - касательный модуль сдвига, соответствующий малому упрочнению и выбираемый как вычислительный параметр для улучшения сходимости, - тангенс угла внутреннего трения. По деформационной теории решена модельная задача о действии перфорированного штампа на область песчаного грунта кубической формы с ребром 1 м, штамп - квадратный со стороной 0,4 м. Расчеты выполнялись для деформаций 6, 8 и 10 %. На верхней границе области задавалось перемещение штампа равномерно на первых шести, восьми и десяти шагах по времени соответственно. Остальная часть верхней границы свободна от напряжений. Боковые и нижняя границы области непроницаемы и закреплены в направлении, перпендикулярном границе. Вся верхняя граница проницаема, через штамп также может просачиваться жидкость. Общее время - 47,5 с, число шагов по времени - 20, величина шага по времени - 2,375 с. Зависимость коэффициента фильтрации от пористости n (выраженной в долях единиц) задавалась по формуле Козени [41; 42]. Свойства грунта задавались по литературным источникам [43], как для среднезернистого песка: (начальный коэффициент фильтрации до деформирования), (начальная пористость), (рис. 1), (для функции нагружения (4)), Расчеты проводились в собственном программном коде. Результаты расчетов приведены на рис. 2 и 3. На рис. 2 показаны графики зависимости давления воды от времени в двух точках, находящихся на разной глубине - 0,1 м (рис. 2, а) и 0,4 м (рис. 2, б) под серединой штампа при общей деформации 6, 8 и 10 %. На рис. 3 показано распределение давления воды на 10-м шаге по времени, когда только произошла деформация грунта 10 % в результате опускания штампа. На рис. 3, а показано распределение давления во всей области, а на рис. 3, б - в вертикальной плоскости, проходящей через середину штампа. На рис. 4 приведено распределение порового давления в вертикальной плоскости на 11-м (рис. 4, а) и 15-м (рис. 4, б) шагах по времени. Можно заметить, что сразу после опускания штампа непосредственно в области под штампом формируется зона повышенного порового давления (см. рис. 3, б). На следующих шагах по времени (см. рис. 4), по мере просачивания воды через перфорированный штамп, давление воды под штампом постепенно падает, распределение порового давления в исследуемой области становится более равномерным. На 11-м шаге по времени величины порового давления снижаются примерно в 2 раза (см. рис. 4, а), на 15-м шаге по времени - становятся ниже на два порядка в результате оттока жидкости через верхнюю границу модели (см. рис. 4, б). а б Рис. 2. Зависимость давления воды от времени в точке под серединой штампа на глубине: а - 0,1 м; б - 0,4 м Fig. 2. Dependence of water pressure on time at a point under the middle of the stamp at a depth: а - 0.1 m; б - 0.4 m а б Рис. 3. Распределение давления воды на 10-м шаге по времени в области (а), в плоскости (б) (общая деформация 10 %) Fig. 3. Water pressure distribution at the 10th time step in the area (а), in the plane (б) (total deformation 10 %) а б Рис. 4. Распределение давления воды в плоскости на 11-м (а) и 15-м (б) шагах по времени (общая деформация 10 %) Fig. 4. Water pressure distribution in the plane at the 11th time step (а), at the 15th (б) time step (total deformation 10 %) Заключение В работе получена физически и геометрически нелинейная связанная формулировка задачи консолидации в скоростях. Для описания процесса деформирования твердой фазы грунта применялся лагранжев подход, для описания движения жидкости - метод ALE. С помощью метода ALE моделируется относительное движение жидкости через подвижную лагранжеву сетку для твердого каркаса. Для численного решения нелинейной задачи консолидации была получена линеаризованная вариационная постановка. Связывание уравнений деформирования скелета грунта и фильтрации жидкости осуществляется методом Узавы. Решена модельная задача об упругопластическом деформировании пористого водонасыщенного песчаного грунта под действием перфорированного штампа. Для моделирования определяющих соотношений упругопластического деформирования грунта используется обобщение модели С.С. Григоряна и деформационной теории пластичности на конечные деформации. В будущем при решении практических задач предполагается использовать теорию течения.

About the authors

S. V Sheshenin

Lomonosov Moscow State University

N. B Artamonova

Lomonosov Moscow State University

References

Statistics

Views

Abstract - 287

PDF (Russian) - 172