книги / Несущая способность конструкций в условиях теплосмен
..pdfосновных принципов механики (например, минимизацией потенциальной механической энергии системы, находящейся в равновесии).
Процедура формирования разрешающих соотношений рассматривается на примере решения задачи для твердого деформируемого тела, находящегося в плоско-деформированном состоянии.
Предполагается, что вершины i, j и k треугольного конечного элемента на плоскости имеют координаты {xi, yi}, {xj, yj}, {xk, yk} (рис. 13.3). Для такого конечного элемента можно определить три линейные пробные функции, или функции формы ϕi(x, y), ϕj(x, y) и ϕk(x, y). Например, функция ϕi(x, y) может быть записана в виде
ϕi (x, y ) = α i + βi x + γ i y,
где αi, βi, γi – константы, значениякоторых зависятот координатузлов:
αi = (xj yk − xk yj ) |
, βi = (yj − yk ) |
, γi = (xk − xj ) . |
Величина численно равна удвоенной площади треугольного |
||
конечного элемента ijk. |
Следует отметить, |
что > 0 лишь в том |
случае, когда нумерация вершин треугольника производится в направлении, противоположном движению часовой стрелки. Особенность сконструированной ϕi(x, y) функции заключается в том, что в своем узле75 эта функция равна 1, а в двух других обращается в 0.
ϕi
z |
k |
y |
|
i |
j |
x
Рис. 13.3. Кусочно-линейная пробная функция ϕi на двумерном конечном элементе
75 Иными словами, в соответствующей одноименной вершине «i» треугольного конечного элемента.
211
|
ρfxϕi |
|||
|
|
ρf |
ϕ |
|
|
|
y |
i |
|
{F} = |
ρfx |
ϕ j |
||
|
|
|
dΩ , px, py – компоненты вектора поверхностной на- |
|
Ω |
ρfyϕ j |
|||
|
ρf |
ϕ |
|
|
|
|
x |
|
k |
|
ρfyϕk |
|||
|
|
|
|
|
грузки; fx, fx – компоненты вектора массовых сил; ρ – плотность материала.
Согласно выражениям (13.1) можно получить соотношения между компонентами вектора перемещения u, v и тензора деформа-
ции εxx, εyy, εxy: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
εxx |
= |
∂u |
|
= |
|
∂ |
|
(uiϕi + u j ϕ j + uk ϕk ) = |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
∂x |
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
= u |
|
∂ϕ |
|
+ u |
|
∂ϕ j |
|
+ u |
|
|
∂ϕ |
|
= u |
β |
|
+ u |
|
β |
|
+ u |
|
β |
|
; |
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
i ∂x |
|
|
|
|
|
j |
|
|
|
|
|
|
|
|
|
k |
|
∂x |
|
|
i |
|
i |
|
|
j |
|
|
j |
k |
|
|
k |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
ε yy |
= |
∂v |
|
= |
|
∂ |
|
(vi ϕi + v j ϕ j + vk ϕk ) = |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
∂y |
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
= v |
∂ϕ |
i |
|
+ v |
|
|
|
∂ϕ j |
|
+ v |
|
|
∂ϕ |
k |
= v |
γ |
|
+ v |
γ |
|
|
+ v |
γ |
|
; |
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
i ∂y |
|
|
|
|
j |
|
|
|
|
|
|
|
|
|
k |
|
|
|
i |
|
i |
|
|
j |
|
|
j |
|
k |
|
|
k |
|
|
|
|
|
|
|
|
||||||||
ε |
|
= 1 |
|
∂u |
+ ∂v = |
|
|
∂ |
(u |
ϕ + u |
|
ϕ |
|
+ u |
ϕ |
|
)+ |
∂ |
(v |
ϕ + v |
ϕ |
|
+ v |
ϕ |
|
) = |
||||||||||||||||||||||||||||
xy |
|
|
|
|
|
|
j |
k |
|
|
j |
k |
||||||||||||||||||||||||||||||||||||||||||
|
2 |
∂y |
|
|
|
|
|
∂y |
|
|
|
|
|
i |
i |
|
|
|
|
|
j |
|
|
|
|
k |
|
|
|
|
∂x |
|
|
|
i |
i |
|
|
|
j |
|
|
k |
|
|
|||||||||
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
= u |
∂ϕ |
i |
+ u |
|
|
∂ϕ j |
|
|
|
+ u |
|
|
∂ϕ |
k |
|
+ v |
∂ϕ |
i |
+ v |
|
∂ϕ j |
+ v |
|
|
∂ϕ |
k |
|
= |
|
|
|
||||||||||||||||||||
|
|
|
i ∂y |
|
|
∂y |
|
|
|
|
|
|
|
∂y |
|
|
|
|
j ∂x |
|
|
|
|
∂x |
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
j |
|
|
|
|
|
|
|
k |
|
|
|
|
i ∂x |
|
|
|
|
|
k |
|
|
|
|
|
|
|
= ui γ i + u j γ j + uk γ k + viβi + v jβ j + vk βk .
Полученные выражения записываются в матричной форме:
{ε} = [A]{u},
213
|
|
|
|
|
|
|
|
|
|
|
|
|
ui |
|
ε |
|
|
β 0 β |
|
0 β |
|
0 |
|
v |
|
||||
|
|
|
|
i |
||||||||||
|
xx |
, |
|
i |
γi |
|
j |
γ j |
|
k |
|
, |
u j |
|
где {ε} = εxy |
[A] = |
0 |
0 |
0 |
γk |
{u} = |
. |
|||||||
ε |
|
γ β γ |
|
β γ β |
|
vj |
||||||||
xy |
|
i |
i |
|
j |
j |
|
k |
k |
|
|
|
||
|
|
|
|
uk |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
v |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
Ранее было показано, что закон Гука в случае плоскодеформированного состояния записывается в форме
σ |
|
|
E |
|
|
1− 2ν |
2ν |
0 |
|
ε |
|
|
||
|
|
xx |
|
|
|
|
|
|
|
|
xx |
|
||
{σ} = σyy |
= |
|
|
2ν |
1− 2ν |
0 |
ε yy |
= [D]{ε}, |
||||||
(1+ ν)(1− |
2ν) |
|||||||||||||
|
|
|
|
|
|
|
|
0 |
0 |
|
|
|
|
|
σxy |
|
|
|
|
|
1− 2ν |
εxy |
|
||||||
|
|
|
|
E |
1− 2ν |
2ν |
0 |
|
|
|
|
|
||
где [D] = |
|
|
|
|
2ν |
1− 2ν |
0 |
. |
|
|
|
|
||
|
|
|
|
|
|
|
|
|||||||
|
(1+ ν )(1− 2ν ) |
|
0 |
|
0 |
|
|
|
|
|
|
|||
|
|
|
1− 2ν |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
С учетом предыдущего соотношения получается
{σ} = [D]{ε} = [D][A]{u}.
Подстановка этого выражения в формулу (13.2) приводит к разрешающему соотношению метода конечных элементов
[A]T [D][A]{u}dΩ ={P} + {F}.
Ω
Матрица-столбец {u} содержит числовые значения компонент вектора перемещения вершин (узлов) треугольного конечного элемента и может быть вынесена из-под знака интеграла:
[B]{u} ={P} + {F}. |
(13.3) |
В выражении (13.3) использовано обозначение матрицы жесткости конечного элемента:
[B] = [A]T [D][A]dΩ.
Ω
214
Вид записи системы линейных алгебраических уравнений (13.3), сформированной для произвольного треугольного конечного элемента, моделирующего физические свойства элементарного материального фрагмента твердого деформируемого тела, практически совпадает с выражением (12.7), полученным для произвольного стержневого элемента: в обоих случаях [B] – матрицы жесткости
элементов, |
{u} – узловые перемещения конечных (стержневых) |
элементов, |
{P} + {F} – внешние нагрузки и массовые силы, дейст- |
вующие на фрагмент твердого деформируемого тела (балочного элемента).
13.3. Процедура ансамблирования
Выполняется процедура ансамблирования конечных элементов в единую композицию, аппроксимирующую всю область, занятую твердым деформируемым телом (рис. 13.4, 13.5).
а |
б |
Рис. 13.4. Процедура ансамблирования: аппроксимация параллелепипеда (а) с помощью набора тетраэдральных конечных элементов (б)
Формирование общей системы линейных алгебраических уравнений, позволяющей определить узловые перемещения конечных элементов, выполняется аналогично алгоритму, рассмотренному в примере реализации метода стержневых элементов. Система линейных алгебраических уравнений каждого конечного элемен-
215
та (13.3) расширяется до полной системы уравнений, содержащей все искомые узловые перемещения и все узловые нагрузки и массовые силы. Сложение всех расширенных матриц приводит к построению полной системы линейных алгебраических уравнений для всего рассматриваемого тела.
13.4. Силовые и кинематические граничные условия
Учет силовых граничных условий рассматривается на примере двух конечных элементов, имеющих общую сторону (см. рис. 13.5, а). При объединении этих конечных элементов (см. рис. 13.5, б) в единую композицию учитывается условие их механического взаимодействия76, поэтому при ансамблировании неизвестные внутренние реакции взаимодействия Fx , Fy , Fz и Fx′, Fy′, Fz′ будут исключены из
разрешающих соотношений и останутся только внешние нагрузки {P} и массовые силы {F} , действующие на рассматриваемое твер-
дое деформируемое тело.
|
wi |
|
wi |
|
wi |
|
vi |
|
|
||
ui |
ui |
vi |
ui |
vi |
|
|
|
||||
|
|
|
|||
|
|
Fz |
|
|
|
|
|
Fy |
F’y |
F’x |
|
|
|
|
|
||
|
Fx |
|
|
|
|
|
|
|
|
F’z |
|
|
|
а |
|
б |
|
Рис. 13.5. Конечные элементы, имеющие общую сторону (а), и ансамблирование конечных элементов в единую композицию (б)
76 Вновь используется аксиома статики, согласно которой два тела взаимодействуют с силами Fx, Fy и Fz, равными по величине и направленными в противоположные стороны (см. рис. 13.5, а).
216
Кинематические граничные условия (в перемещениях) задаются для тех вершин (узлов) конечных элементов, перемещения которых известны (например, величина осадки или сдвига) или определяются условиями симметрии, неподвижного закрепления и проч.
13.5. Реализация алгоритма метода конечных элементов
Определяется решение сформированной системы алгебраических уравнений для определения поля перемещения частиц деформируемого твердого тела.
Для решения полной системы линейных алгебраических уравнений и определения узловых компонент вектора перемещения используются прямые (Гаусса, квадратного корня) или итерационные методы (Якоби, Зейделя, простых итераций и др.).
Найденные узловые перемещения {u} вершин конечных эле-
ментов позволяют найти распределение компонент тензора деформации с помощью соотношения {ε} = [ A]{u} и далее – компоненты
тензора напряжения с использованием закона Гука {σ} = [D]{ε} во
всех конечных элементах, а значит, во всей области, занятой твердым деформируемым телом.
13.6. Анализ результатов численного решения задачи
Полученное численное решение необходимо подвергнуть критическому анализу с точки зрения достоверности полученных данных, соответствия результатов расчетов поведению реального материального объекта. Лучший способ проверки адекватности результатов вычислительного моделирования – это сопоставление результатов вычислительного моделирования с данными натурного эксперимента. Такое сравнение позволяет корректно определить погрешность вычислительной модели.
В случае невозможности проведения натурного эксперимента выполняют ряд вычислительных экспериментов на последовательности сгущающихся сеток конечных элементов. Например, первона-
217
чально выполняется вычислительный эксперимент на сравнительно грубой сетке, содержащей N конечных элементов.
Второй вычислительный эксперимент проводится с использованием конечно-элементной сетки, более точно аппроксимирующей исследуемую область и содержащей 2N конечных элементов. Очевидно, что во втором случае характерный размер конечного элемента уменьшится вдвое.
Третий вычислительный эксперимент проводится с использованием сетки, содержащей 4N конечных элементов с характерным размером, в четыре раза меньшим, чем первоначальный, и т.д.
При выполнении такой последовательности экспериментов изучаются изменения каких-то характеристик исследуемого объекта: потенциальной энергии, компонент тензоров деформации, напряжения, вектора перемещения в рассматриваемой области. Если при уменьшении характерного размера конечных элементов (увеличение числа конечных элементов) наблюдаемые параметры устанавливаются (стабилизируются), делается вывод о правильности получаемого решения.
Дополнительно может проводиться исследование устойчивости полученного решения по отношению к возмущению каких-то параметров. Например, возмущаемым параметром может быть выбран диаметр отверстия в рассматриваемом теле, наблюдаемым результатом – уровень концентрации напряжения вблизи этого отверстия. Если малое возмущение диаметра (его уменьшение или увеличение) приводит к малому изменению напряжения в зоне концентрации, решение считается устойчивым.
Дальнейшее подробное рассмотрение полученных результатов – распределение компонент вектора перемещения u и v, компонент εxx, εyy, εxy тензора деформации и σxx, σyy, σxy тензора напряжения, интенсивности деформации εi и интенсивности напряжения σi – позволяет выполнить полный анализ напряженно-деформированного состояния: выявить места концентрации деформаций и напряжений; определить области деталей узлов и конструкций, опасные с точки зрения различных теорий прочности; выявить нерационально спроекти-
218
рованные элементы конструкций, несущие чрезмерную или недостаточную нагрузку; оценить запас прочности детали, узла или агрегата при превышении допустимого значения нагрузки и проч.
13.7.Тестирование метода конечных элементов
Вкачестве тестовой задачи, как правило, используется известное точное или приближенное решение упрощенной задачи.
|
y |
3 |
y |
Fy2 |
2 |
|
3 |
|
|
4 |
4 |
|
Fx |
|
|
||
|
|
|
Fy1 |
F1 |
|
4 |
|
|
|
|
|
|
x |
|
Fy |
|
|
|
|
|
|
|
|
|
|
|
P |
|
x |
|
Fy3 |
|
|
4 |
|
Ми |
1 |
|
Fx3 |
|
|
|||
2 |
|
|
|
Fx |
|
|||
|
|
|
1 |
|
F5 |
Fx5 |
2 |
|
|
|
|
|
y |
x |
|||
|
|
|
|
|
|
|
|
|
|
а |
|
|
|
б |
|
|
|
Рис. 13.6. Геометрическая схема осадки длинной полосы (а)
исхема кинематических граничных условий (б)
Вкачестве тестовой задачи для оценки точности алгоритма метода конечных элементов, реализованного при использовании допущения о плоско-деформированном состоянии, рассматривается задача об осадке длинной стальной полосы с квадратным поперечным сечением
размером 2×2 м2, зажатой между двумя гладкими горизонтальными плитами (рис. 13.6, а). Каждая из плит осаживается на величину . Необходимо определить, на какую величину сместятся боковые стороны этой полосы в результате деформирования. Для упрощения принимается, что массовые силы отсутствуют, то есть {F} = {0} .
Поскольку форма поперечного сечения является симметричной, можно рассматривать лишь четверть исследуемой области (рис. 13.6, б). Кинематические граничные условия указаны на этом же рисунке. Поскольку плиты, деформирующие полосу, гладкие,
219