- •ВВЕДЕНИЕ
- •1.1. Идея метода перемещений
- •1.2. Дискретизация (разбиение) расчетной схемы
- •1.3. Типы конечных элементов
- •Вопросы для самопроверки
- •2. ВЫВОД РАЗРЕШАЮЩИХ УРАВНЕНИЙ МКЭ
- •2.1. Аппроксимация перемещений на элементе
- •2.2. Общие требования к аппроксимирующим функциям перемещений и рекомендации по их применению
- •2.3. Механический смысл коэффициентов жесткости
- •2.4. Определение коэффициентов жесткости из принципа возможных перемещений
- •2.5. Приведение нагрузки на элементе к узловой
- •2.6. Вывод разрешающих уравнений МКЭ из принципа минимума полной потенциальной энергии
- •Вопросы для самопроверки
- •3. ПОСТРОЕНИЕ ОБЩЕЙ СИСТЕМЫ УРАВНЕНИЙ
- •3.1. Глобальная матрица жесткости
- •3.1.1. Плоская задача
- •3.1.2. Изгиб балки
- •3.2. Вектор свободных членов
- •3.3. Учет граничных условий
- •3.4. Определение перемещений и усилий в балке
- •Вопросы для самопроверки
- •4. ПРИМЕРЫ РАСЧЕТА БАЛОК
- •4.1. Общая схема алгоритма расчета МКЭ
- •4.2. Расчет статически определимой балки МКЭ
- •4.2.1. Расчет балки, представленной одним конечным элементом
- •4.2.2. Расчет балки, разбитой на два конечных элемента
- •4.3. Расчет статически неопределимой балки МКЭ
- •4.4. Задания к самостоятельной работе
- •Библиографический список
4.2.2.Расчет балки, разбитой на два конечных элемента
1.Рассмотрим задачу, приведенную на рис. 4.1, разбив балку на два конечных элемента с двумя степенями свободы в узле и с заданным направлением глобальных осей (рис. 4.3).
Рис. 4.3. Конечные элементы балки
2. Локальные матрицы жесткости [k1] и [k2 ] определим, подста-
вив численные значения l1 = |
|
l2 = 5 м в выражение (2.45): |
|
|
||||||||||
|
|
|
0,096 |
0,24 |
−0,096 |
0,24 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
И |
|
|
|
|
|
[k1 |
|
|
|
0,24 |
0,8 |
|
−0,24 |
0,4 |
|||||
|
] = [k 2 ] = EJ |
|
|
|
. |
|||||||||
|
|
|
|
− |
0,096 |
−0,24 |
0,096 |
−0,24 |
|
|||||
|
|
|
|
|
|
|
|
Д |
|
|
|
|||
|
|
|
|
|
0,24 |
0,4 |
|
−0,24 |
0,8 |
|||||
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
А |
|
|
|
|
|
||
|
3. Взаимосвязь между локальной и глобальной нумерацией ко- |
|||||||||||||
нечных элементов, узлов бстепеней свободы представим в табличной |
||||||||||||||
форме. |
|
и |
|
|
|
|
|
|
Таблица 4.1 |
|||||
|
|
|
|
|
|
|
|
|
||||||
|
|
Взаимосвязь локальной и глобальной нумераций |
|
|
||||||||||
|
|
С |
|
|
|
|
|
|||||||
|
Номер |
|
Номера узлов |
|
|
Номера степеней свободы |
||||||||
|
элемента |
|
локальные |
|
|
|
глобальные |
|
локальные |
|
глобальные |
|
||
|
I |
|
1 |
|
|
|
1 |
|
|
1,2 |
|
|
1,2 |
|
|
|
2 |
|
|
|
2 |
|
|
3,4 |
|
|
3,4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
II |
|
1 |
|
|
|
2 |
|
|
1,2 |
|
|
3,4 |
|
|
|
2 |
|
|
|
3 |
|
|
3,4 |
|
|
5,6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
4. Глобальная матрица жесткости формируется путем объединения общих степеней свободы соседних элементов. Глобальная матрица жесткости имеет размерность 6×6 и имеет вид
53
|
|
0,096 |
0,24 |
−0,096 |
0,24 |
0 |
0 |
|
||
|
|
0,24 |
0,8 |
−0,24 |
0,40 |
0 |
0 |
|
||
[K]= EJ |
|
|
||||||||
−0,096 |
−0,24 |
0,096 |
+0,096 |
0,024 |
−0,024 |
−0,24 |
0,24 |
. |
||
|
|
0,24 |
0,4 |
0,24 |
−0,24 |
0,8 |
+0,8 |
−0,24 |
0,24 |
|
|
|
|
||||||||
|
|
0 |
0 |
−0,096 |
−0,24 |
0,096 |
−0,024 |
|||
|
|
0 |
0 |
0,24 |
0,4 |
−0,024 |
0,8 |
|
||
|
|
|
5. Для построения локальных векторов узловых сил используем выражение (4.2), в котором произведем замену l на l1 = l2 = 5 м:
|
2,5 |
|
|
2,083 |
|
{FI } = {FII } = |
. |
|
|
2,5 |
|
Д |
|
|
−2,083 |
||
|
|
|
А |
|
|
6. Формирование глобального вектораИсвободных членов осу- |
ществляется по принципу объединения компонентов локальных векторов соседних конечных элементов, соответствующих общим степеням свободы. Глобальный вектор свободных членов {P } имеет раз-
мерность 1×6 и состоит из локальных векторов первого и второго
элементов {FI }и {FII }: |
бq l |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
||
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
С |
|
|
|
|
q l12 |
|
|
|
|
|
2,5 |
|
||||||
|
|
q l |
|
|
|
12q l |
|
|
|
|
|
2,083 |
|
|||||
|
|
1 |
2 |
|
|
|
||||||||||||
|
|
|
|
+ |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
5,0 |
|
|||||
{P } = |
|
|
2 |
|
|
|
|
|
2 |
|
|
= |
||||||
|
|
|
|
|
|
|
|
|
|
0,0 |
. |
|||||||
|
|
− |
q l12 |
|
q l 22 |
|
|
|
|
|||||||||
|
|
12 |
+ |
12 |
|
|
|
2,5 |
|
|||||||||
|
|
|
|
|
||||||||||||||
|
|
|
|
|
q l2 |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− 2,083 |
||
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
||
|
|
|
− |
|
q l |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
12 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
54 |
|
|
|
|
|
|
7. Определим номера степеней свободы в глобальной системе, в направлении которых невозможно перемещение конструкции. Это номера 1 и 5. В глобальной матрице жесткости столбцы и строки с этими номерами вычёркиваются:
|
|
|
0,096 |
|
0,24 |
−0,096 |
|
0,24 |
|
0 |
|
0 |
|
||||||||||
|
|
|
0,24 |
|
0,8 |
|
−0,24 |
|
0,40 |
|
0 |
|
0 |
|
|||||||||
[K ]= EJ |
|
|
|
|
|
|
|
||||||||||||||||
−0,096 |
−0,24 |
0,192 |
|
|
0 |
|
|
−0,24 |
|
0,24 |
. |
||||||||||||
|
|
|
0,24 |
|
0,4 |
|
|
0 |
|
|
|
1,6 |
|
−0,24 |
|
0,4 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
0 |
|
|
0 |
−0,096 |
−0,24 |
|
0,096 |
|
−0,024 |
|||||||||||
|
|
|
0 |
|
|
0 |
|
0,24 |
|
0,4 |
|
−0,024 |
|
0,8 |
|
||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
И |
|
|
|
|||||
8. Система разрешающих уравнений (4.3) принимает вид |
|
||||||||||||||||||||||
|
|
0,8 |
−0,24 |
0,4 |
|
|
0 |
|
∆2 |
|
|
|
|
2,083 |
|
|
|
||||||
|
−0,24 |
0,192 |
0 |
|
0,24 |
∆ |
|
|
|
|
5 |
|
|
|
|||||||||
EJ |
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
− |
|
|
|
|
= 0. |
|
|
0,4 |
0 |
|
1,6 |
|
0,4 |
|
∆ |
|
|
|
0 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
б |
|
|
|
4 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
0 |
0,24 |
|
0,4 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
0,04Д∆ |
|
− 2,083 |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
|
|
|
|
||
В результате решения этойАсистемы получаем значения узловых |
|||||||||||||||||||||||
перемещений, составляющ х вектор {∆}: |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
∆1 |
|
|
w1 |
|
|
|
|
|
|
|
0 |
|
|
|
|
|||
|
|
|
|
и∆ |
|
θ |
1 |
|
|
|
|
|
41,665 EJ |
|
|
|
|||||||
|
|
{∆} = |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
∆3 |
|
= |
w2 |
|
= |
|
|
130,20 EJ |
. |
|
|||||||||||
|
|
|
|
|
∆ |
|
|
|
θ |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
С4 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
∆ |
|
|
w |
3 |
|
|
|
|
|
|
|
0 |
|
|
|
|
|||
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
− 41,665 EJ |
|
|
|
|||||
|
|
|
|
|
∆ |
|
|
θ3 |
|
|
|
|
|
|
|
||||||||
|
|
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Для построения эпюры прогибов воспользуемся формулой (2.2), в которую подставим полученные значения узловых перемещений:
55
|
|
|
|
|
|
|
|
|
|
|
|
|
41,665 |
|
|
|
|
|
130,2 |
|
|
|
|
|
|
|||||
w |
(x ) = N |
(x ) 0 |
+ N |
|
(x ) |
|
|
|
|
|
+ N |
|
(x ) |
|
|
+ N |
|
(x ) 0 |
; |
|||||||||||
|
|
|
EJ |
|
|
EJ |
|
|||||||||||||||||||||||
1 1 |
|
1 |
|
1 |
|
|
|
2 |
1 |
|
|
|
3 |
|
1 |
|
|
|
|
|
4 |
1 |
|
|||||||
w |
(x |
|
|
|
|
(x |
|
) |
130,2 |
+ N |
|
|
(x |
|
) 0 + N |
|
(x |
|
) 0 − N |
|
(x |
|
) |
41,665 |
||||||
|
) = N |
|
|
|
EJ |
|
|
|
|
|
|
|
|
EJ |
, |
|||||||||||||||
2 |
|
2 |
|
1 |
|
2 |
|
|
|
|
2 |
|
2 |
|
|
3 |
|
2 |
|
|
|
4 |
|
2 |
|
|
где w1(х1) – прогиб на первом участке (0 < x1 < l1); w2(х2) – прогиб на втором участке (0 < x2 < l2).
Рис. 4.4. Эпюра прогибов w EI
|
|
Д |
|
||
9. Из курса сопротивления материалов известна зависимость |
|||||
А |
MИ |
|
|||
d 2 w |
|
(4.5) |
|||
d x |
2 |
= − |
EJ |
. |
|
б |
|
|
|
|
|
Для построения эпюры моментов воспользуемся формулой (3.23): |
M = − EJ w′′(x) = − EJ [N1′′(x) ∆1 + N2′′(x) ∆2 + N3′′(x) ∆3 + N4′′(x) ∆4 ],
где N”i(x) – вторые про зводные функций формы, определяемые по |
|
формулам (2.22). |
и |
Учитывая найденные значения узловых перемещений, запишем |
выражения для изгибающих моментов в пределах первого и второго |
||||||||||||||||||||||||||||||
|
|
|
|
С |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
конечных элементов в следующем виде: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
M |
|
|
|
|
′′ |
|
|
+ N |
′′ |
(x ) |
41,665 |
+ N |
|
′′ |
(x ) |
130,2 |
+ N |
′′ |
(x )0 |
|
||||||||||
1 |
(x ) = −ЕJ N (x )0 |
2 |
|
|
EJ |
|
3 |
EJ |
4 |
; |
||||||||||||||||||||
|
1 |
|
1 1 |
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
1 |
|
||||||||||||
M |
|
(x |
|
|
′′ |
|
) |
130,2 |
+ N |
′′ |
(x |
|
)0 |
+ N |
′′ |
(x |
|
)0 − N |
′′ |
(x |
|
|
) |
41,665 |
||||||
2 |
2 |
) = −ЕJ N (x |
2 |
EJ |
|
2 |
2 |
3 |
2 |
4 |
2 |
EJ |
. |
|||||||||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Так как вторые производные координатных функций являются функциями линейными, то и выражения для изгибающих моментов в
56
пределах каждого элемента являются линейными функциями координаты х. Поэтому для построения эпюры моментов на каждом элементе достаточно иметь значения изгибающих моментов в узловых точках:
Рис. 4.5. Эпюра изгибающих моментов М∆, кН м
Изображенная на рис. 4.5 эпюра М∆ является узловой, так как
она является результатом узловых воздействий, то есть перемещений
изгибающих моментов. Эпюра Мр от распределеннойИ нагрузки строится в основной системе метода перемещений (рис. 4.6).
∆1,…,∆6. Поскольку балочный элемент испытывает кроме узловых
воздействий еще и воздействие равномерно распределенной нагруз-
ки, то это должно быть учтено при построении окончательной эпюры
Для получения суммарной эпюры моментов ординаты эпюры
М∆ от узловой нагрузки суммируются с ординатами эпюры Мр от ме- |
|||||
стной нагрузки: |
|
|
|
Д |
|
|
M = M ∆ + M р. |
(4.5) |
|||
|
|
||||
|
|
|
А |
|
|
|
|
б |
|
|
|
|
и |
|
|
|
|
С |
|
|
|
|
Рис. 4.6. Эпюра Мр от распределенной нагрузки
Окончательная эпюра М представлена на рис. 4.7.
57