bmp
.pdf275
Між кінцевими характеристиками в глобальній та в локальній системах координат може бути встановлений формальний зв’язок:
|
|
e Te e , |
|
|
(14.13) |
|||
|
|
re Tere , |
|
|
(14.14) |
|||
|
|
pe Tepe |
, |
|
(14.15) |
|||
де |
|
|
|
|
|
|
|
|
cos |
sin |
0 |
|
0 |
0 |
0 |
||
|
||||||||
|
|
cos |
|
|
|
|
|
|
sin |
0 |
|
0 |
0 |
0 |
|||
|
0 |
0 |
1 |
|
0 |
0 |
0 |
|
Te |
|
|
|
|
|
|
|
. |
0 |
0 |
0 |
|
cos |
sin |
|
||
|
|
0 |
||||||
|
0 |
0 |
0 |
|
sin |
cos |
0 |
|
|
0 |
0 |
0 |
|
0 |
0 |
|
|
|
|
1 |
||||||
|
|
|
|
|
|
|
|
e |
|
|
|
|
|
|
|
Квадратна матриця Te називається матрицею перетворення стержневого елемента e.
Для стержня с шарнірами на кінцях матриця перетворення набуває вигляду:
|
cos |
sin |
|
0 |
0 |
|
|
|
|
||||||
|
|
sin |
cos |
|
0 |
0 |
|
T |
|
|
|
||||
0 |
0 |
|
cos |
|
|||
e |
|
|
sin |
||||
|
|
0 |
0 |
|
sin |
|
|
|
|
|
cos |
||||
|
|
|
|
|
Матриці жорсткості стержневих скінченних елементів
Між кінцевими реакціями і кінцевими переміщеннями існує зв’язок, який у глобальній системі координат може бути записаний у формі:
re ke e , |
(14.16) |
де ke – матриця жорсткості елемента e в глобальній системі координат. Для стержневого скінче-
ного елемента першого типу матриця має вигляд:
|
k11 |
k12 |
k13 |
k14 |
k15 |
k16 |
|
|
|||||||
|
|
|
|
k22 |
k23 |
k24 |
k25 |
k26 |
|
|
|||||
|
k21 |
|
|
||||||||||||
ke |
k |
31 |
k |
32 |
k |
33 |
k |
34 |
k |
35 |
k |
36 |
|
(14.17) |
|
|
|
|
|
|
|
|
|
||||||||
|
k41 |
k42 |
k43 |
k44 |
k45 |
k46 |
|
|
|||||||
|
k |
51 |
k |
52 |
k |
53 |
k |
54 |
k |
55 |
k |
56 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
k |
|
k |
|
k |
|
k |
|
k |
|
k |
|
|
|
|
|
61 |
62 |
63 |
64 |
65 |
66 |
e |
|
||||||
|
|
|
|
|
|
|
|
276
Елементи матриці жорсткості – це реакції на кінцях, тобто кінцеві реакції стержневого скінченного елемента, які зумовлюються його кінцевими переміщеннями. Перший індекс характеризує номер кінцевої реакції, а другий – номер кінцевого переміщення, яке цю реакцію викликало. Наприклад, елемент k25 – це кінцева реакція r2 (див.рис.14.9,б) від дії примусового переміщення
(див.рис.14.9,а) за умови, що всі інші кінцеві переміщення дорівнюють нулю. Величини ре-
акцій залежать від типу скінченного елемента. Так, для скінченного стержневого елемента, що має жорсткі вузли з обох сторін (рис.14.3,а), елементи матриці жорсткості мають значення:
k |
|
k |
|
k |
|
|
k |
|
|
f cos2 |
|
|
12i |
|
sin2 ; |
|
|
|
|
|
||||||||
|
44 |
|
|
41 |
|
|
|
|
|
|
|
|||||||||||||||||
11 |
|
|
14 |
|
|
|
|
|
|
|
|
l2 |
|
|
|
|
|
|
|
|
|
|
||||||
k |
22 |
k |
k |
25 |
k |
52 |
f sin2 |
|
|
12i |
cos2 ; |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
55 |
|
|
|
|
|
|
|
|
|
|
|
l2 |
|
|
|
|
|
|
|
|
|
|
||||
k33 k66 4i; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
12i |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(14.18) |
|||
k |
|
k |
|
k |
|
|
k |
|
k |
k |
|
k |
|
|
|
k |
|
f |
|
|
sin cos ; |
|||||||
|
|
|
|
|
|
|
|
|
|
l2 |
||||||||||||||||||
12 |
|
21 |
|
45 |
|
|
54 |
|
|
|
15 |
|
51 |
|
|
|
24 |
|
42 |
|
|
|
|
k13 k31 k16 k61 k46 k64 k34 k43 6li sin ; k35 k53 k56 k65 k23 k32 k26 k62 6li cos ; k36 k63 2i.
Тут позначено:
i |
EI |
; |
f |
EA |
, |
(14.19) |
|
l |
|
|
l |
|
|
де l – довжина, EI – жорсткість стержня на згин, EA – жорсткість стержня на поздовжні деформації.
Для стержневого скінченого елементу четвертого типу (шарнірне приєднання з обох сторін) матриця жорсткості в глобальній системі координат набуває вигляду:
|
|
k |
k |
|
k |
k |
|
|
|
|
|
|
|
||||||
|
|
11 |
12 |
|
13 |
14 |
|
|
|
k |
e |
|
k21 |
k22 |
|
k23 |
k24 |
|
|
|
k32 |
|
k33 |
k34 |
|||||
|
k31 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
k41 |
k42 |
|
k43 |
k44 |
e |
|
|
|
|
|
|
(14.20)
|
f cos2 |
|
f sin cos |
|
f sin2 |
f sin cos |
||
|
|
|
||||||
|
|
|
f sin2 |
|
f sin cos |
f sin2 |
|
|
|
f sin cos |
|
|
|||||
|
|
|
|
|
|
|
|
|
f sin2 |
|
f sin cos |
|
f cos2 |
f sin cos |
|
||
|
|
|||||||
|
|
|
f sin2 |
|
f sin cos |
f sin2 |
|
|
|
f sin cos |
|
|
|||||
|
|
|
|
|
|
|
e |
|
|
|
|
|
|
|
Визначення зусиль у стержнях
Розв’язок рівняння рівноваги (14.8) визначає вектор вузлових переміщень
280
1-7 |
1 |
0 |
3 |
0 |
0 |
3 |
-1 |
0 |
1 |
100 |
0,333 |
33,33 |
1-4 |
2 |
0 |
3 |
0 |
6 |
3 |
1 |
0 |
1 |
100 |
0,333 |
33,33 |
3-8 |
1 |
8 |
3 |
8 |
0 |
3 |
-1 |
0 |
1 |
100 |
0,333 |
33,33 |
3-5 |
1 |
8 |
3 |
8 |
6 |
3 |
1 |
0 |
1 |
100 |
0,333 |
33,33 |
6-9 |
1 |
12 |
3 |
12 |
0 |
3 |
-1 |
0 |
1 |
100 |
0,333 |
33,33 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Усі жорсткі вузли дискретної моделі мають три ступені вільності. Вузол 4, до якого всі елементи примикають шарнірно, має два ступеня вільності. Схема можливих вузлових переміщень наведена на рис.14.15.
Рис.14.15
Вектор вузлових переміщень має 26 компонентів:
1 |
2 |
3 |
4 |
5 |
6 |
24 |
25 |
26 T . |
Переміщення опорних вузлів дорівнюють нулю:
18 19 20 21 22 23 24 25 26 0 .
Виключивши нульові компоненти з вектора вузлових переміщень маємо остаточно:
1 |
2 |
3 |
4 |
5 |
6 |
15 |
16 |
17 T . |
Таким чином дискретна модель має 17 невідомих вузлових переміщень і отже 17 ступенів вільності.
Побудова вектора вузлових навантажень
Вектор вузлових навантажень має компоненти, які відповідають компонентам вектора вузлових переміщень:
F F1 F2 |
F3 F4 F5 F6 F15 F16 |
F17 T . |
Схема дії зовнішніх вузлових навантажень наведена рис.14.16.