1.КОНЕЧНЫЕ ЭЛЕМЕНТЫ
ИРАЗРЕШАЮЩИЕ УРАВНЕНИЯ
При расчете конструкции МКЭ рассматриваемое тело делится на конечное число малых частей – конечных элементов (КЭ), связанных между собой в определенном числе точек – узлах. Если предположить, что перемещения в любой точке конечного элемента можно выразить в зависимости от перемещений в узлах, то задача сводится к вычислению перемещений в узлах. Перемещения узлов области и на контуре определяются системой уравнений, которые представляют собой условия
равновесия узлов, а также условия на контуре. Эти уравнения могут быть |
||
|
|
И |
сформированы несколькими путями, в том числе на основе принципа ми- |
||
нимума потенциальной энергии. |
|
|
1.2. Дискретизация (разбиение) расчетной схемы |
||
|
А |
|
Сущность дискретизации (разбиения) тела по МКЭ состоит в |
||
следующем: |
|
|
-рассматриваемое тело с помощьюДвоображаемых линий или |
||
поверхностей разбивается на элементы конечных размеров. Эти |
||
элементы называются конечными элементами, а их совокупность – |
|
сеткой конечных элементов; |
|
|
С |
-конечные элементыбсоединены в конечном числе точек, кото- |
|
рые находятся на контуре элемента и называются узловыми точками |
|
или узлами; |
и |
-состояние каждого конечного элемента (перемещение, деформация, напряжение и т.д.) описывается с помощью интерполирующих функций и конечного числа параметров в узлах (степеней свободы), представляющих основные неизвестные МКЭ (рис. 1.1).
Число узловых перемещений и соответствующих им узловых сил равно числу степеней свободы. Числом степеней свободы называют количество независимых параметров, определяющих положение в пространстве узловой точки конечного элемента и всей дискретной расчетной схемы конструкции. В соответствии с этим различают количество степеней свободы одного узла mу, конечного элемента mэ и всей конструкции m. Если вся расчетная схема имеет n узлов, то общее число степеней свободы определяется по формуле
6
m = n my. |
(1.1) |
Рис. 1.1. Плоская расчетная схема (а) и ее конечно-элементное представление (б)
четыре прямоугольных элемента. Каждая узловаяИточка на плоскости имеет две степени свободы. В качестве степеней свободы произвольной точки k (рис. 1.2,б) рассматриваются линейные
Степени свободы имеют определенный физический смысл.
Они представляют собой узловые перемещения.
На рис. 1.2 представлен пример плоской области, разбитой на
перемещения этой точки ик и vк соответственно в направлении осей X |
||||
|
|
|
|
Д |
и Y. Таким образом, количество степеней свободы одного узла равно |
||||
двум: my = 2. |
|
|
А |
|
|
|
|
||
|
|
б |
|
|
|
и |
|
|
|
|
С |
|
|
|
Рис.1.2. Конечно-элементная дискретизация области
Так как один прямоугольный конечный элемент содержит четыре узла, то количество степеней свободы такого элемента равно восьми: тэ = 8. Соседние конечные элементы могут иметь общие узлы.
Полное количество узлов расчетной схемы, представленной на рис. 1.2,а, равно девяти: n = 9. Общее количество степеней свободы всей системы т = n ту = 9×2 = 18.
7
В методе перемещений основную систему получают путем наложения дополнительных связей на систему заданную. Причем связи накладывают по направлению искомых узловых перемещений или степеней свободы. Для произвольного узла к в дополнительно наложенных связях возникают реактивные силы Rxk и Ryk, которые по существу являются узловыми силами взаимодействия между конечными элементами. Число узловых сил совпадает с числом узловых перемещений и равняется числу степеней свободы.
Конечные элементы (КЭ) можно классифицировать по следующим признакам:
- одномерные – стержневые КЭ, используемые для расчета как плоских, так и пространственных стержневых конструкций (табл.
1.1);
|
|
|
А |
|
|
|
|
|
|
Таблица 1.1 |
||||||
|
Параметры одномерных конечныхИэлементов |
|
|
|
|
|
||||||||||
Тип |
Стержневые элементы |
Узловые |
|
|
|
|
Узловые силы |
|||||||||
Д |
|
|
||||||||||||||
|
|
|
перемещения |
|
|
|
|
|
|
|
|
|||||
|
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
б{∆}= |
∆ |
|
|
|
|
F |
|
|
||||||
І |
С |
∆1 |
|
|
|
|
{F}= F1 |
|
||||||||
|
|
|
|
|
2 |
|
|
|
|
|
2 |
|||||
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
∆ |
|
|
w |
|
F |
|
|
|
||||
|
|
|
|
1 |
|
|
|
|
1 |
|
|
|
z1 |
|
||
ІІ |
|
|
∆2 |
|
= |
|
θ1 |
|
M1 |
|
||||||
|
|
{∆}= |
|
|
|
|
|
|
{F}= |
|
|
|
|
|
||
|
|
|
∆3 |
|
|
w2 |
|
M 2 |
|
|||||||
|
|
|
∆ |
|
|
θ |
2 |
|
F |
|
|
|
||||
|
|
|
|
4 |
|
|
|
|
|
|
|
z2 |
|
|||
|
|
|
∆1 |
|
|
u1 |
|
Fx1 |
|
|||||||
|
|
|
∆ |
|
|
w |
|
F |
|
|
|
|||||
|
|
|
|
2 |
|
|
|
|
1 |
|
|
|
z1 |
|
||
|
|
|
∆ |
|
= |
|
θ |
|
|
M |
|
|
|
|||
ІІI |
|
|
{∆}= |
3 |
|
|
u |
1 |
|
{F}= |
F |
1 |
|
|||
|
|
|
|
∆ |
|
|
|
2 |
|
|
|
|
|
|||
|
|
|
4 |
|
|
|
x2 |
|||||||||
|
|
|
∆ |
|
|
w |
|
F |
|
|
|
|||||
|
|
|
|
5 |
|
|
|
|
2 |
|
|
|
z2 |
|
||
|
|
|
∆6 |
|
θ2 |
M 2 |
||||||||||
|
|
|
8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
- двумерные, треугольной и четырехугольной формы, используемые для решения плоской задачи теории упругости и расчета пластин и оболочек (табл.1.2);
Таблица 1.2
Параметры двумерных конечных элементов
Тип |
Плоские элементы |
|
|
Узловые |
|
|
|
|
|
|
|
Узловые силы |
||||||||||
|
перемещения |
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
∆1 |
|
u1 |
|
|
|
Fx1 |
|
||||||||||
|
|
|
|
|
∆ |
|
v |
|
|
|
F |
|
|
|||||||||
|
|
|
|
|
|
|
2 |
|
|
|
1 |
|
|
|
|
y1 |
|
|||||
|
|
|
|
|
|
∆ |
|
|
|
|
|
|
|
|
F |
|
|
|||||
ІV |
|
|
|
{∆}= |
|
|
|
u |
|
|
|
{F}= |
|
|||||||||
|
|
|
|
3 |
= |
v |
2 |
|
|
x2 |
|
|||||||||||
|
|
|
|
|
|
∆ |
|
|
2 |
|
|
|
|
F |
|
|
||||||
|
|
|
|
|
|
4 |
|
|
|
|
y2 |
|||||||||||
|
|
|
|
|
∆ |
|
u |
|
|
|
F |
|
|
|||||||||
|
|
|
|
|
|
|
5 |
|
|
|
|
3 |
|
|
|
|
x3 |
|
||||
|
|
|
|
|
|
Д |
|
|
|
F |
|
|
||||||||||
|
|
|
|
|
|
∆ |
|
|
v |
3 |
|
|
|
|
||||||||
|
|
|
|
|
|
6 |
|
|
|
|
|
y3 |
|
|||||||||
|
|
|
|
|
|
∆1 |
|
|
u1 |
|
Fx1 |
|||||||||||
|
|
|
|
А |
|
И |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
{∆} |
∆2 |
|
v1 |
|
|
Fy1 |
|
|||||||||||
V |
|
б |
|
3 |
|
|
|
|
|
2 |
|
{F} |
|
x2 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Fy4 |
||||||
|
|
|
|
|
∆8 |
|
v4 |
|
||||||||||||||
|
и |
|
∆1 |
|
|
|
w1 |
|
|
|
Fz1 |
|
||||||||||
|
|
|
∆2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
θx1 |
|
M x1 |
||||||||||||||
VІ |
|
|
∆ |
|
|
= |
θ |
|
|
|
{F}= |
M |
|
|
||||||||
С |
|
|
{∆}= |
|
3 |
|
|
|
|
|
y1 |
|
|
F |
y1 |
|||||||
|
|
|
|
|
∆ |
|
|
|
|
w |
|
|
|
|
|
|
||||||
|
|
|
|
|
4 |
|
|
|
|
2 |
|
|
z2 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
∆ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
θ |
|
|
|
|
M |
|
|
||||||
|
|
|
|
|
|
12 |
|
|
|
|
y4 |
|
|
|
y4 |
|||||||
- трёхмерные, объёмные, или пространственные, используемые для расчета сложных пространственных конструкций и массивных тел (табл. 1.3).
9
Таблица 1.3
Параметры трехмерных конечных элементов
Тип |
Объемные элементы |
|
Узловые |
|
|
|
|
Узловые силы |
||||||
|
перемещения |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
∆1 |
|
|
u1 |
|
Fx1 |
|
||||
|
|
|
|
|
∆ |
|
|
v |
|
|
F |
|
||
|
|
|
|
|
2 |
|
|
|
1 |
|
|
y1 |
|
|
|
|
|
|
|
∆3 |
|
= |
|
|
|
|
|
|
|
VIІ |
|
|
{∆}= |
|
w1 |
|
{F}= Fz1 |
|
||||||
|
|
|
|
|
∆ |
|
|
|
u |
2 |
|
|
F |
|
|
|
|
|
4 |
|
|
x2 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∆12 |
|
|
w4 |
|
Fz4 |
|
||||
|
|
|
|
∆1 |
|
|
u1 |
|
Fx1 |
|
||||
|
|
|
|
|
∆ |
|
|
v |
|
F |
|
|||
|
|
|
|
|
2 |
|
|
|
|
1 |
|
|
y1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
VIІІ |
|
|
|
Д |
|
|
|
|
{F}= Fz1 |
|
||||
|
|
{∆}= |
∆3 |
|
= w1 |
|||||||||
|
|
|
|
|
∆ |
|
|
|
u |
|
|
|
F |
|
|
|
|
|
|
4 |
|
|
|
|
2 |
|
|
x2 |
|
|
|
|
|
|
|
И |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
∆ |
|
|
w |
|
F |
|
||||
|
|
|
|
|
24 |
|
|
|
|
8 |
|
|
z8 |
|
|
|
б |
|
|
|
|
|
|
|
|
|
|
|
|
|
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.4. С стемыАкоординат |
|
|
|
|
|
|
|||||||
|
С |
|
|
|
|
|
|
|
|
|
|
|
|
|
В методе конечных элементов используют две системы координат: локальную (местную) глобальную (общую) (рис. 1.3).
Локальные координаты применяют для описания отдельных элементов, глобальные – для описания всей конструкции.
Рис.1.3. Системы координатных осей
10