РАЗВИТИЕ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ ДЛЯ РАСЧЕТА

advertisement
ПРОБЛЕМЫ ПРОЧНОСТИ И ПЛАСТИЧНОСТИ, вып. 73, 2011 г.
УДК 519.6
РАЗВИТИЕ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ ДЛЯ РАСЧЕТА
ВОДОНАСЫЩЕННОГО ОДНОРОДНОГО ОСНОВАНИЯ
Т.В. Мальцева, Т.Ю. Володина
Тюменский государственный университет
Рассматривается задача о нагружении водонасыщенного основания (тела)
при стабилизированном состоянии (в трехмерной постановке). Напряженнодеформированное состояние двухфазного основания описывается системой
дифференциальных уравнений типа Ламе, в которую входят младшие производные. Для решения математической задачи используется «ажурный» вариант метода конечных элементов, который приводит к понижению порядка системы линейных алгебраических уравнений при одинаковом количестве элементов. Построена новая матрица коэффициентов системы линейных алгебраических уравнений.
Ключевые слова: водонасыщенное основание, двухфазное тело, уравнения
типа Ламе, метод конечных элементов, «ажурная схема».
Введение
Строительство объектов нефтегазодобывающего комплекса и транспортировка углеводородного сырья в основном ведется на водонасыщенных глинистых грунтах, заторфованных, заболоченных территориях.
Ряд натурных и лабораторных экспериментов указывают на факт существования избыточного давления в поровой воде [1−3], когда кривая порового давления
σ l начиная с некоторого времени t stab имеет постоянное конечное значение σ ost
(рис.1), которое не учитывается моделями фильтрационной консолидации.
σl
σ0
σost
tstab
t
Рис. 1
Исследования, связанные с построением и развитием математической модели,
описывающей вклад поровой воды в напряженное и деформированное состояние
скелета грунта с помощью младших производных по координатам, а также получением некоторых аналитических решений классических задач типа Фламана и Буссинеска, приводятся в работах [4, 5].
150
С целью проведения численных расчетов напряженно-деформированного состояния водонасыщенного основания при действии различных нагрузок необходимо получить новый вариант метода конечных элементов (МКЭ), учитывающий остаточные поровые давления.
Постановка задачи
На части границы двухфазного тела (скелет грунта + поровая вода) действуют
внешние силы Q(х1, х2, х3), на части − нулевые кинематические условия. Влияние
избыточных остаточных поровых давлений описывается с помощью кинематической
модели [4], которая представляет собой систему трех дифференциальных уравнений эллиптического типа с положительными постоянными коэффициентами G, λ,
b, c относительно вектора перемещений скелета грунта u = (u1, u2, u3):

∂ 2u
∂u
∂θ
−  (G + λ )
+ G∆u i + b 2i + c i
∂xi
∂x i
∂x i

G=

 = Fi , i = 1, 2, 3,


(1)
νES
ES
E
E
, λ=
, b = l2 , c = l , θ = divu,
2(1 + ν)
(1 + ν)(1 − 2ν)
ℵh
ℵ
с граничными условиями:
ui
S1
= 0, ~tij u j
S2
= Qi ( x1 , x2 , x3 ),
(2)
~
tij = λni ∂ j + (G + bδij ) n j ∂ i + G δ ij n k ∂ k ,
где ∂ j = ∂ / ∂x j , δij − символ Кронекера, ES − модуль деформации, ν − коэффициент
Пуассона скелета грунта, El − механическая постоянная поровой воды, ℵ − безразмерная величина ( 0 < ℵ < 1) , определяемая из эксперимента, h − высота сжимаемой толщи, Fi − объемные силы.
Согласно принципу Лагранжа скалярно умножим уравнения (1), записанные в
операторном виде, на вектор истинных перемещений u :
− ((A + B + C) u, u ) = (F, u ),
(3)
где A = (G + λ )grad div + G∆ − оператор Ламе, B = b( ∂ 2 / ∂x12 , ∂ 2 / ∂x 22 , ∂ 2 / ∂x32 ),
C = c(∂ / ∂x1 , ∂ / ∂x2 , ∂ / ∂x3 ) − операторы, описывающие влияние поровой воды на
скелет грунта. В [4] установлена связь между скалярными произведениями и полной энергией деформации, на основании которой строится матрица коэффициентов системы линейных алгебраических уравнений (СЛАУ).
«Ажурная» схема метода конечных элементов
Для численной реализации модели воспользуемся «ажурной» схемой метода
конечных элементов, изложенной в работе [6] для решения задач теории упругости.
В качестве конечного элемента выберем центральный тетраэдр lmnp параллелепипеда (один в центре и четыре по краям) (рис. 2). Такая схема называется «ажурной», так как в расчете участвует только один тетраэдр, что приводит к существенному (в несколько раз) снижению порядка СЛАУ.
3
На основании линейной аппроксимации u i = α i0 + ∑ j=1 α ij x j (i = 1, 2, 3) вектор искомых перемещений u представим в виде линейной комбинации узловых
151
перемещений
{δ} = (u1m , u2m , u3m , u1n , u2n , u3n , u1l , u2l , u3l , u1p , u2p , u3p ) :
ui =
1
[(d m + q m x1 + ym x2 + t m x3 )u im + ( d n + q n x1 + y n x 2 + t n x3 )u in +
∆
+ ( d l + ql x1 + y l x2 + t l x3 )u il + ( d p + q p x1 + y p x2 + t p x3 ) uip ] ( i = 1, 2, 3), (4)
где
d m = x1n x 2l x3p − x1n x2p x3l − x1l x 2n x3p + x1l x 2p x3n + x1p x 2n x3l − x1p x2l x3n ,
qm = x2l x3p − x2p x3l − x2n x3p + x2p x3n + x2n x3l − x2l x3n ,
ym = x1n x3p − x1n x3l − x1l x3p + x1l x3n + x1p x3l − x1p x3n ,
tm = x1n x l2 − x1n x2p − x1l x 2n + x1l x2p + x1p x 2n − x1p x2l ,
1
1
∆=
1
1
x1m
x1n
x1l
x1p
x2m
x2n
x2l
x2p
x3m
x3n
.
x3l
x3p
Верхний индекс указывает на номер узла тетраэдрального элемента. Коэффициенты d, q, y, t находятся путем круговой перестановки индексов.
p
m
x1
x2
l
n
x3
а)
б)
Рис. 2
Из уравнений Коши
∂u j
1  ∂u
ε ij =  i +

2  ∂x j ∂x i




относительные деформации внутри конечного элемента выражаются через искомые узловые перемещения {δ}:
{ε} = [ N ]{δ},
где
152
 qm

 0
 0
[N ] = 
 ym

 tm
 0

0
ym
0
qm
0
tm
0
0
tm
0
qm
ym
qn
0
0
yn
tn
0
0
yn
0
qn
0
tn
0
0
tn
0
yl
yn
ql
0
0
yl
tl
0
0
yl
0
ql
0
tl
0
0
tl
0
ql
yl
qp
0
0
yp
tp
0
0
yt
0
qp
0
tp
0 

0 
tp 

0 

qp 
y p 
− геометрическая матрица, соответствующая положению тетраэдра.
Построение СЛАУ
Для получения матрицы, отвечающей сумме операторов (А + В), согласно принципу Лагранжа используется равенство работ внешних {δ}T {F } (T − операция
транспонирования) и внутренних сил.
Удельная работа внутренних сил, отвечающих скелету грунта, с учетом влияния поровой воды через оператор В равна {ε}T {σ} , где {ε}T = {δ}T [ N ]T . Напряжения и деформации для скелета грунта связаны законом {σ} = [ D]{ε} , где D −
матрица механических характеристик водонасыщенного грунта для тетраэдрального элемента.
От удельной работы перейдем к работе внутренних сил в пределах объема конечного элемента, в результате получим:
{δ}T ∫ [ N ]T [ D][ N ] dV {δ} = {δ}T [k s ]{δ}, [k s ] = [ N ]T [ D][ N ]∆,
V
s
где [ k ] − матрица для скелета грунта с учетом поровой воды, зависящая от механических и геометрических характеристик элемента.
Построим матрицу, соответствующую третьему операторному слагаемому
(−Cu, u) в системе уравнений (3):
{δ}T [k s ]{δ} + (−Cu, u) = {δ}Т {F },
(5)
 ε11 


∂ ui
( − Cu , u ) = ∫ c
u i dV = с ∫ ( u1 u2 u3 )  ε 22  dV.
∂x i
ε 
V
v
 33 
Перемещения u запишем в матричном виде согласно формулам (4):
{u} = [ M *]{δ},
(6)
(7)
где
 ξ*m 0
0 ξ*n 0 0 ξ*l 0 0 ξ*p 0

 0 ξ*m 0
0 ξ*n 0 0 ξ*l 0 0 ξ*p

0
0 ξ*m 0 0 ξ*n 0 0 ξ*l 0
0
[M *] = 
 0
0
0
0 0 0 0 0 0 0
0

0
0
0 0 0 0 0 0 0
0
 0
 0
0
0
0 0 0 0 0 0 0
0

0 

0 

ξ*p 
,
0 

0 
0 
153
здесь ξ*k = d k + qk x1 + y k x2 + tk x3 , k = l, m, n, p. Отметим, что нулевые элементы в
[M*] добавлены, чтобы выполнить требование одинаковой размерности матричных слагаемых (5).
Продолжим преобразование выражения (5):
( −Cu, u) = ∫ {δ}T [ M * ]T [ Dl ]{ε}dV = {δ}T ∫ [ M * ]T [ Dl ] [ N ]{δ}dV =
V
V
T
= {δ}
∫ [M
* T
] dV [ Dl ] [ N ]{δ}, [ Dl ] = с[ E ],
(8)
V
где [E] − единичная матрица. В формулы перемещений (7) входят аргументы х1, х2,
х3, поэтому подынтегральное выражение содержит произведения xidV, i = 1, 2, 3.
Интегралы с точностью до константы вычислим c помощью координат центра тяжести тетраэдра:
∫ xi dV = xicV , xic =
( x1p − x1n ) 2 + ( x 2p − x 2n ) 2 + ( x3p − x3n ) 2
2
V
.
В формуле (8) объемный интеграл заменим на произведение [M]{δ}, тогда с
учетом координат центра тяжести тетраэдра окончательно получим:
( −Cu, u) = {δ}T [ k l ]{δ},
где [ k l ] = [ M ]T [ D l ][ N ].
Выражение (5) после сокращения на {δ}T примет вид:
([ N ]T [D][ N ] + [M ]T [Dl ][N ]){δ} = {F } или ([k s ] + [k l ]){δ} = {F }.
Сумму матриц [ k s ] + [ k l ] назовем матрицей жесткости для тетраэдрального
двухфазного элемента.
Заключение
На основе матриц жесткости строится глобальная матрица коэффициентов
СЛАУ. Она имеет блочную структуру. Аналогичная матрица жесткости была построена для случая плоского конечного элемента (треугольника, прямоугольника) в
работе [7].
Литература
1. Амарян, Л.С. Свойства слабых грунтов и методы их изучения / Л.С. Амарян. − М.:
Недра, 1990. − 220 с.
2. Зехниев, Ф.Ф. Стабилизация оснований с плоскими вертикальными песчаными дренами: автореф. дис… канд. техн. наук: 05.23.02 / Зехниев Фархад Фархадович. − М., 1988. −
25 с.
3. Набоков, А.В. Исследование напряженно-деформированного состояния основания из
водонасыщенной глины: дис… канд. техн. наук: 05.23.02: защищена 05.07.2004: утв.
08.10.2004 / Набоков Александр Валерьевич. − Тюмень, 2004. − 142 с.
4. Мальцева, Т.В. О разрешимости смешанной задачи, отвечающей обобщенному
оператору Ламе / Т.В. Мальцева // Вестник ТюмГУ. − 2006. − №5. − С. 230−234.
5. Мальцева, Т.В. Моделирование двухфазного тела с учетом несущей способности
жидкой фазы / Т.В. Мальцева, Е.Р. Трефилина // Математическое моделирование. − 2004. −
Т.16, №11. − С. 47−60.
154
6. Чекмарев, Д.Т. Об одном классе двухмерных схем МКЭ / Д.Т. Чекмарев, К.М. Гладильщикова // Проблемы прочности и пластичности: Межвуз. сб. / Нижегород. ун-т. − 2006.
− Вып. 68. − С. 236−241.
7. Мальцева, Т.В. Сопоставление матриц жесткости при расчете двухфазной полуплоскости / Т.В. Мальцева, Т.В. Салтанова // Вестник ТюмГУ. − 2007. − №5. − C. 25−33.
[21.09.2011]
DEVELOPMENT OF THE FINITE ELEMENT METHOD FOR CALCULATION
OF WATER-SATURATED HOMOGENEOUS BASE
T.V. Maltseva, T.Yu. Volodina
A problem of loading of water-saturated base (body) in a stabilized state (a spatial case) is considered.
The stressed-strained state of the two-phase base is described by Lame's differential equations
which include the lower order derivatives. To solve the mathematical problem a rare mesh version
of the finite element method is used that leads to the decrease in the order of a system of linear
algebraic equations with the same number of elements. A new coefficient matrix for the system of
linear algebraic equations is constructed.
Keywords: water-saturated base, two-phase body, Lame's equations, finite element method, rare
mesh scheme.
155
Download