Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

bmp

.pdf
Скачиваний:
39
Добавлен:
05.02.2016
Размер:
36.09 Mб
Скачать

272

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

{ 4

5

T

j

 

5

 

6}j .

 

 

6

 

 

 

 

 

 

j

 

 

 

Повний вектор переміщень скінченного елемента e в локальній системі координат матиме вигляд:

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

T

.

 

 

 

 

 

e

 

i

 

j e

 

1

2

3

4

5

6

e

 

Так само можна записати вектор кінцевих реакцій стержня в локальній системі координат:

 

 

 

 

 

 

T

 

 

 

 

 

 

T

 

 

 

e

 

re

ri

rj

r1

r2

r3

r4

r5

r6

e .

Для стержня, який приєднується до вузлів шарнірно (стержень ферми), вектори кінцевих переміщень і реакцій в локальній системі координат матимуть вигляд:

 

 

 

 

 

 

T

 

 

 

 

 

 

δ

T

 

 

 

 

e

 

 

 

 

e

 

 

i

 

 

j

 

 

1

2

3

 

4

e

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

T

 

 

 

 

 

 

 

e

 

 

 

 

 

re

ri

rj

r1

 

r2

r3

r4

e

 

 

Кінцеві переміщення і реакції можуть бути представлені у вигляді компонентів, які являють собою проекції відповідних величин на осі глобальної і локальної систем координат. Позначення і додатні напрями зазначених компонентів для стержня жорстко прикріпленого до вузлів представ-

лено на рис.14.9,а,б, а для стержня, який приєднується до вузлів шарнірно на рис.14.9,в,г.

Нумерація цих величин також суворо фіксована і аналогічна нумерації в локальній системі координат.

Рис.14.9

Кінцеві переміщення і кінцеві реакції в глобальній системі координат можуть бути представлені у вигляді векторів:

для стержня жорстко прикріпленого до вузлів

273

e i

 

j T 1

 

2

3

4

 

6 Te ;

 

 

 

 

e

 

 

 

 

 

 

 

 

 

T

r1

 

 

 

 

 

T

re

ri

rj e

r2

r3

r4

r5

r6 e .

для стержня прикріпленого до вузлів шарнірно

 

 

 

 

T

1

 

 

T

e

i

 

j e

2

4

4 e

 

 

 

T

r1 r2

 

T

 

re

ri

rj e

r3

r4 e

 

Рис.14.10

Крім кінцевих переміщень, на стержневий скінченний елемент можуть діяти рівномірно розподілені навантаження, які орієнтуються за осями або локальної (рис.14.10,а) або глобальної (рис.14.10,б) системи координат. Зв’язок між цими двома системами навантажень встановлюється співвідношеннями

q

 

qxly

 

 

cos

qylx

 

sin ;

 

 

 

 

 

 

 

 

 

x

 

 

 

l

 

 

l

(14.10)

 

 

 

 

 

 

q

 

qxly

sin

qylx

cos

 

 

 

 

 

 

 

 

y

 

 

 

l

 

 

l

 

 

 

 

 

 

 

 

або навпаки

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qx

qxl

 

cos

qyl

sin ;

 

 

 

 

 

 

 

 

ly

 

ly

(14.11)

 

 

 

qxl

 

 

 

qyly

 

 

qy

sin

 

cos .

 

 

 

 

 

 

 

 

 

lx

 

 

lx

 

Так само за локальними (рис.14.11,а) або за глобальними координатними осями (рис.14.11,б) орієнтуються вантажні реакції, які зумовлюються розподіленими вздовж стержня навантаженнями.

274

а

 

б

Рис.14.11

Нумерація вантажних реакцій також суворо визначена.

Вантажні реакції в локальній і в глобальній системі координат також можуть бути представлені у вигляді векторів:

 

p

 

p

 

p

 

p

 

p

 

p

T

,

pe

1

2

3

4

5

6 e

 

 

p1

p2

p3

p4

p5

 

 

T

 

pe

 

p6 e .

 

Для стержня, що має на кінцях шарніри

 

p

 

p

 

p

 

p

 

T

pe

1

2

3

4 e

 

p1

p2

 

p3

 

 

T

 

pe

 

 

p4 e

 

Цілком очевидно, що в разі відсутності розподілених навантажень ці вектори є нульовими. При наявності розподілених навантажень значення елементів векторів вантажних реакцій залежать від граничних умов по кінцям стержня. Так, для скінченного елементу першого типу (рис.14.3,а), тобто для елемента, що має жорсткі вузли на обох кінцях

 

 

 

 

 

q l

 

 

 

 

 

 

 

 

 

 

 

q l

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

x

y

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

qyl

 

 

 

 

 

 

 

 

 

 

 

qylx

 

 

 

p

1

 

2

 

 

 

 

 

p1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

2

 

 

 

 

p

2

 

 

 

 

2

 

 

 

 

 

2

 

 

qyl

 

 

 

 

 

 

 

 

 

 

qyly

 

 

 

p

3

 

12

 

 

 

 

p3

 

 

12

 

(14.12)

 

 

;

 

 

 

pe

p

 

 

 

q l

 

 

 

pe

p

 

 

 

 

q

l

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

x

 

 

 

 

 

4

 

x

y

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

p

5

 

 

qyl

 

 

 

 

 

p5

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

p

6

 

 

 

qylx

 

 

 

 

 

6 e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

q l2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

qyly

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12

e

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

Матриці жорсткості і вектори вантажних реакцій для скінченних елементів інших типів наведено в додатках 1.

275

Між кінцевими характеристиками в глобальній та в локальній системах координат може бути встановлений формальний зв’язок:

 

 

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

 

 

 

 

 

 

 

 

 

5 1

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) визначає вектор вузлових переміщень

 

277

K 1F .

(14.21)

Таким чином стають відомими переміщення всіх вузлів дискретної моделі. Внаслідок нерозривності деформацій кінці стержнів, які приєднуються до вузлів, матимуть такі ж самі переміщення, що й вузли. Отже для кожного стержня дискретної моделі можна побудувати вектор кінцевих переміщень у глобальній системі координат:

e { i

 

 

 

1

 

 

 

xi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

yi

 

 

T

3

 

 

 

i

(14.22)

 

j

}e

 

 

 

 

 

 

.

 

 

4

 

 

 

xj

 

 

 

 

5

 

 

 

yj

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

e

 

Сумарні кінцеві реакції стержня в локальній системі координат (рис.14.12,а)

 

 

T

 

 

 

 

 

T

se

{si

sj }e

{s1

s 2

s 3

s 4

s 5

s 6}e .

складаються з кінцевих реакцій ri , що зумовлені переміщеннями вузлів дискретної моделі (14.22),

і кінцевими силами pi , які викликано дією розподіленого на стержні зовнішнього навантаження:

se re pe .

(14.23)

Рис.14.12

Сумарні кінцеві реакції можуть бути обчислені за формулою:

 

e

 

 

se he

pe .

(14.24)

Вигляд матриці he залежить від граничних умов на кінцях стержня. Так, для стержня, що має на обох кінцях жорсткі вузли

 

 

 

 

 

 

 

fC

 

 

 

 

 

12i

S

 

 

 

 

 

 

 

 

 

 

 

 

l

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6i

S

 

h

e

 

 

 

 

l

 

 

 

 

 

 

fC

 

 

 

 

 

 

 

 

12i

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6i

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

278

 

 

 

 

fS

 

 

 

0

 

 

fC

 

 

fS

 

 

0

 

 

 

12i

C

 

 

6i

12i

S

 

12i

C

 

 

6i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l2

 

 

 

l

 

 

l2

 

l2

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6i

 

C

 

4i

 

6i

 

S

 

 

6i

C

 

2i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

0

 

 

 

 

 

l

 

 

 

 

 

l

 

 

 

 

0

 

 

(14.25),

 

 

 

fS

 

 

 

 

 

 

 

fC

 

 

 

 

fS

 

 

 

 

 

 

 

12i

 

 

6i

 

12i

 

12i

 

 

 

6i

 

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

2

 

 

 

l

 

 

l2

 

 

 

l2

 

 

l

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6i

 

C

 

2i

 

6i

 

S

 

 

6i

C

 

4i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

l

 

 

 

 

 

 

 

e

 

а для стержня з шарнірами на кінцях

 

 

 

fC

fS

fC

fS

 

 

 

 

0

0

0

0

 

 

h

e

 

.

(14.26)

 

fC

fS

fC

fS

 

 

 

 

 

0

0

0

0

 

 

 

 

 

 

 

Тут позначено: S sin , C cos .

Нумерація сумарних кінцевих реакцій для стержня, що має на кінцях шарнірні вузли представлена на рис.14.12,б.

Приклавши сумарні кінцеві реакції до стержня (рис.14.12), можемо обчислити внутрішні зусилля в будь-якому його перерізі розглядаючи рівновагу однієї з його частин. Наприклад, у перерізі з абсцисою x (рис.14.12,а), маємо:

M x M x s 2 x s 3 q y x

2

;

лів.

 

 

 

2

 

лів.

 

 

Qx y

s 2 q y x;

(14.27)

лів.

 

 

Nx x s1 q x x.

 

14.2. Приклад розрахунку плоскої рами

Необхідно виконати статичний розрахунок плоскої рами (рис.14.13). На раму діють зосере-

джені сили

P1 10кН і P2 20кН , а також рівномірно розподілені навантаження з інтенсивністю

q1 2,4кН м

і q2 3,0кН м.

279

Рис.14.13

Стояки мають жорсткості на згинальні та поздовжні деформації відповідно EI 1, EA 100 , а

ригелі – вдвічі більшу: EI 2, EA 200 .

Перехід до дискретної (скінченно-елементної) моделі

Для переходу від заданої розрахункової схеми до дискретної (скінченно-елементної) моделі розіб’ємо схему на окремі стержневі елементи, які поєднуються у вузлах. За вузли приймаємо точки сполучення окремих стержнів, опорні точи, а також точку прикладення зовнішньої зосередже-

ної сили P1 . Пронумеруємо вузли в довільному порядку. Початок глобальної системи візьмемо у вузлі 7. Дискретна модель наведена на рис.14.14.

Фізико-геометричні характеристики скінченних елементів

Для кожного стержня дискретної моделі необхідно обчислити довжину (l), нахил осей локальної системи координат (sin , cos ), по-

гонні жорсткості (i,f), а також визначити тип його граничних умов (рис.14.14). Результати обчислень занесені до табл.14.1.

Рис.14.14

Таблиця 14.1

 

 

Початок

Кінець

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

СЕ

СЕ

 

 

 

 

EI

EA

i

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xi

 

yi

xj

 

y j

 

 

 

 

 

 

 

1-2

1

0

 

3

4

 

3

4

0

1

2

200

0,5

50

2-3

2

4

 

3

8

 

3

4

0

1

2

200

0,5

50

4-5

3

0

 

6

8

 

6

8

0

1

2

200

0,25

25

5-6

1

8

 

6

12

 

3

5

-0,6

0,8

2

200

0,4

40

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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.

281

Згідно до формули (14.3) вектор вузлових навантажень складається з двох векторів. Перший з них – це вектор зовнішніх зосереджених сил і моментів, що безпосередньо діють на вузли:

 

{F

F

F

F

F

F

F

F

F

F

F

F

F

F

F

F

F

}T .

F

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

 

Рис.14.16

Порівнюючи схему завантаження дискретної моделі (рис.14.14) зі схемою компонентів вузлових навантажень (рис.14.16), можна записати:

F5 10;

F13 20.

Усі інші компоненти дорівнюватимуть нулю, тому:

F

{0 0 0 0

10 0 0 0 0 0 0 0

20 0 0 0 0}T .

 

Вектор Q реактивних сил, що передаються на вузли від розподілених на стержні наванта-

жень, також має 17 компонентів:

Q {Q1 Q2 Q3 Q4 Q5 Q6 Q7 Q8 Q9 Q10 Q11 Q12 Q13 Q14

Q15 Q16 Q17}T .

F

 

Схема дії компонентів цього вектора аналогічна компонентам вектора

(рис.14.17). Для обчи-

слення компонентів вектора Q розглянемо завантажені стержні окремо.

 

Рис.14.17

Стержень 1-7 (рис.14.18,а). Скінченний елемент має з обох боків жорсткі вузли і відноситься до першого типу. Навантаження в локальній системі координат (рис.14.18,б):

qx 0;

q y 2,4.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]