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

Строительная механика

.pdf
Скачиваний:
0
Добавлен:
29.11.2025
Размер:
9.21 Mб
Скачать

Ф орм а конечного элем ента (К Э ) предопределяется особенностям и рассчиты ваемого объекта (систем ы , конструкции). В стерж невы х сис­ тем ах за К Э приним аю т стерж ень, как правило, постоянной ж есткости на растяж ение- сж атие и и згиб , или группу взаим освязанны х стерж ­ н ей . Д ля плоских и тонкостенны х континуальны х систем наиболее часто использую тся треугольны е или прям оугольны е (в общ ем слу ­ чае, четы рехугольны е) К Э , при реш ен и и трехм ерны х задач - объем ­ ные К Э в виде тетраэдра или параллелепипеда. В ы бор формы и разм е­ ров К Э оказы вает сущ ественное влияние на результаты расчета, кото­ ры е, естественно, долж ны позволять правильно оценивать напряж ен­ но- деф орм ированное состояние исходной систем ы . П редставление исследуем ой систем ы как достаточно больш ого набора КЭ значитель­ но увеличивает размерность задачи , ведет к повы ш ению точности рас­ чета, но связано со значительны м объемом вы числений. Ч итатели , ин ­ тересую щ иеся вопросам и оценки погреш ности дискретизации, м огут найти соответствую щ ие реком ендации в научной литературе.

Точки, в которых соединяются К Э , называю т узлами. Различаю т узлы жесткие и ш арнирные. В жестком узле предполагается наличие связей,

обеспечиваю щ их неразрывность линейных и угловых перемещ ений К Э ,

примы каю щ их к этому узлу . Связи ш арнирного узла позволяю т сохра­ нить неразрывность линейны х перемещ ений. Узловые перемещ ения и соответствую щ ие им узловые силы принимаю тся за обобщ енны е.

И дея м етода состоит в то м , чтобы описать напряж енно-деф ор­

мированное состояние К Э через обобщ енны е перем ещ ения Z узлов и установить их связь с действую щ ей на систему нагрузкой . Д ля реали­ зации этой идеи необходим о получить м атрицу ж есткости К Э .

Так как ф ункция п ерем ещ ен и й исходн ой систем ы н еи звестн а, то ее н еобходим о зад авать . Е сли в м етоде Р и тц а п редп олагалось, что базисны е ф ункции оп ределяю тся одним вы раж ен и ем на всей облас­ ти си стем ы , то в М К Э реали зуется альтерн ати вн ы й п о д х о д . О н за ­ клю чается в то м , что на каж дом К Э н еизвестны е ф ун кц и и п ерем е­

щ ен и й зам ен яю тся ап п рокси м и рую щ и м и и х таки м

о бр азо м , чтобы

п ерем ещ ен и я всех точек эл ем ен та бы ли вы раж ены

через у зл о вы е .

Д ля одн ом ерны х эл ем ен то в, с учетом зам ечан и я о п остоянстве ж е ­ сткостей , ф ункция п ерем ещ ен и й является то ч н о й (см . раздел 16.7),

для дву м ер н ы х и тр ех м ер н ы х К Э эти ф ункции зап и сы ваю тся при ­

бли ж ен н о , наиболее часто

- в виде

п о л и н о м о в . П одбор

и х п р ед ­

ставляет собой достаточн о

слож н ую

зад ач у . О т удачн ого

реш ен и я

501

ее сущ ественно зави си т точн ость окон чательн ы х

результатов. С

п ом ощ ью

апп рокси м и рую щ и х ф ункций на

основе

вари ац и он н ы х

при н ц и п ов строи тельн ой м ехан и ки реш ается

одна из осн овн ы х за ­

дач М К Э -

определение м атри ц ж есткостей кон ечны х элем ентов.

Так как

каж дый стержень в составе исследуемой

системы имеет

свою ориентацию , то вначале строят матрицы жесткостей в местной системе координат, а затем, при переходе от местной системы к общей, преобразую т их. М атрицу ж есткости всей системы получаю т соответст­ вую щ им объединением матриц ж есткостей отдельных элементов.

Разреш аю щ и е уравн ен и я М К Э зап и сы ваю тся в виде:

 

 

 

 

 

 

R Z + R F

— 0 ,

 

(16.24)

где

R F

-

векто р “ гр у зо в ы х ” р еак ц и й , р ав н ы й вектору у зл о вы х

 

 

нагрузок, взятом у с обратны м знаком .

 

 

П олн ая

нагрузка

на у зел определяется

как сум м а

нагрузок

от

п ри м ы каю щ и х к у зл у элем ентов. Т ак

к ак

в н е у зл о в ая

н агр у зк а

з а ­

м ен яется эк в и в ал ен тн о й у зл о во й

по

н ап р авл ен и я м Z , то векто р

р ек ц и й

R f

[ ,

R 2F ,...,R n F ]

— - F , г д е F — [Fb

F 2,...,F n ] -

векто р у зл о в ы х н агр у зо к .

П о сл е р еш е н и я си с те м ы у р а в н е н и й (1 6 .2 4 ) ст а н о в я т с я и з ­

в е с тн ы м и п ер е м е щ е н и я Z у зл о в в о б щ ей си с те м е к о о р д и н ат . Д л я в ы ч и сл е н и я у с и л и й в к о н еч н о м э л ем ен т е у д о б н о вн ач ал е

н ай т и ве к то р п е р е м е щ е н и й Z ' в м е с т н о й с и с те м е к о о р д и н ат , а

за тем с

п о м о щ ь ю м атр и ц ы ж е с т к о с т и R '

о п р ед ел и ть к о н ц ев ы е

р еак ц и и . Ф о р м у л ы

д л я

с о о т в е т с т в у ю щ и х

п р е о б р а зо в а н и й п р и ­

ве д е н ы

в р а зд е л а х

16.11

и 16.12.

 

Т акая ф орм а расчета соответствует варианту М К Э ”в п ерем ещ е­ н и ях” . О на является наиболее распространенной .

В озм ож ен и другой п одход к реш ен и ю задачи по М К Э . Н ап ря­ ж ен н о -деф орм и рован н ое состояние К Э н еобходим о описать кон еч ­ ны м набором о боб щ ен н ы х узл о вы х сил, а затем у стан ови ть и х связь с нагрузкой. Т акая ф орм а расчета со о тветствует М К Э “в у си ли ях” .

502

16.11. Матрица жесткости стержня в местной системе координат

Существует несколько способов получения матриц жесткостей отдельных стержней. Одним из наиболее простых является способ, основанный на известных положениях метода перемещений.

Конец стержня, примыкающий к жесткому узлу, имеет три сте­ пени свободы: линейные перемещения по горизонтальному и вер­ тикальному направлениям и угол поворота. Соответствующими этим смещениям силовыми факторами будут концевые реактивные силы R1, R2 , R4 , R5 и моменты R3, R6 . (Перемещения, реакции,

матрица жесткости и ее элементы в местной системе координат обо­ значаются буквами со штрихами). Матрица жесткости (матрица еди­

ничных

реакций)

преобразует

вектор

перемещений

Z ' —[ ,

Z 2, Z 3,

Z 4, Z 5,

Z 6 ]

в вектор концевых реакций, то есть

имеет место соотношение:

 

 

 

 

 

 

 

R1"

'1,1

r12

•••

r16

"Z1"

 

 

 

R2

>21

r22

•••

r26

Z2

R' Z '

 

 

R6 _

Г61

r62

•••

r66_

Z

 

 

 

 

 

 

_ 6_

 

 

Положительные направления реакций

Rj

соответствуют поло­

жительным направлениям Z j .

 

 

 

 

 

Элементы матрицы R '

суть реакции в связях, вызываемые еди­

ничными смещениями Z j —1 (рис. 16.12).

 

 

 

В первом столбце записаны значения реакций от Z1 —1, во вто­

ром - от

Z 2 —1 и т. д. Следовательно, для вычисления элементов

матрицы R ' можно использовать данные из таблицы, применяемой в методе перемещений для определения реакций в опорных закреп­ лениях стержня постоянного сечения.

503

t 7*21 = 0

^ 5 1 = 0

Рис. 16.12

Матрица R ' имеет вид:

504

 

EA

 

EA

 

 

 

~ T

 

~ T

 

 

 

12EJ

6EJ

12EJ

6EJ

 

 

/ 3

/ 2

/ 3

/ 2

 

 

6EJ

4E J

6EJ

2EJ

 

R'

/ 2

/

/ 2

/

(16.25)

EA

 

EA

 

 

 

 

 

 

/

6EJ

/

6EJ

 

 

12EJ

12EJ

 

 

/ 3

/ 2

/ 3

/ 2

 

 

6EJ

2EJ

6EJ

4E J

 

 

/ 2

/

/ 2

/

 

Для стержней с другими условиями закрепления элементы мат­ рицы R' вычисляются аналогично.

Рассмотрим способы, основанные на использовании аппрокси­ мирующих функций перемещений. Необходимый порядок расчета по одному из них, например, для защемленного по концам стержня, представляется следующим.

1. В линейно деформируемом стержне продольные и попереч­ ные перемещения сечений стержней не взаимосвязаны. Поэтому функции, описывающие характер изменения перемещений по длине стержня, для них будут различными. Перемещения сечений стерж­ ня вдоль его оси в соответствии с дифференциальным уравнением

N E A и' будем аппроксимировать линейной функцией:

 

и(х) —ai + «4 х .

(16.26)

Изогнутая ось стержня при отсутствии по его длине распределен­ ной нагрузки описывается кривой третьего порядка, что является

следствием дифференциального уравнения v (lV) —0 . Следователь­ но, аппроксимирующий полином третьей степени позволяет точно задать функцию перемещений стержня.

505

П усть:

 

 

 

 

v (х ) = + a 3х + a 5 х

2

3

(16.27)

+ a 6х .

 

В вы раж ениях (16.26) и (16.27)

a f -

неизвестные параметры; число

и х равн о чи слу степ ен ей свободы .

 

 

 

Ф ун кц и и и (х ) и v ( х ) н азы ваю т ф ун кц и ям и формы .

 

У гол п оворота сечен и я стерж н я оп редели тся зн ачен ием

первой

производной:

 

 

 

 

d v

,

о

2

(16.28)

---- =a3+2a5х + 3a6х

.

2. И сп ользуя зави си м ости (16.26), (16.27) и (16.28), п редставим

"d v nT

вектор п ерем ещ ен и й 5 = и, v,

----

в следую щ ем виде:

 

 

 

 

 

d х

 

 

 

 

 

 

5

= L a ,

 

 

(16.29)

где

L

-

м атри ц а коэф ф ициентов:

 

 

 

 

 

 

1

х

 

 

 

 

 

 

L =

х

х2

х 3

(16.30)

 

 

 

1

 

 

 

 

1

3х2

 

 

a

-

вектор н еи звестн ы х парам етров:

 

 

 

 

 

a = [ a a 2 a 3 a 4 a 5 a6 ] .

 

 

3.

Д ля кон ц евы х то ч ек (х = 0, х = l) стерж н я с пом ощ ью вы р аж е­

н и й (16.29) и (16.30) получим :

 

 

 

 

506

Z1

1

 

 

a1

 

 

 

z 2

1

 

 

a 2

 

 

 

z 3

 

1

 

a 3

z ;

1

l

 

a ;

z 5

1

l

l 2

a5

z 6

l 3

 

1

2l

a 6

 

 

3l2

или

 

 

 

 

 

 

Z ' = H a ,

 

(16.31)

где H - матрица связи.

;. Из (16.31) следует, что:

 

 

 

 

(16.32)

при этом

 

 

 

 

1

 

 

 

 

1

1

 

 

H -

= -1/1

 

 

 

1/1

 

 

-3/12

-2/ 1

3/12

-1/ 1

 

2/13

1/12

-2/13

1/12

5.

Вектор перемещений 5 с

помощью зависимостей (16.29)

и (16.32) представляется в виде:

 

 

 

 

5 = L H -

Z '.

(16.33)

Выполнив перемножение матриц, получим:

 

507

 

 

 

 

,

 

х Л

 

 

х

1

 

 

 

 

 

и 11— l JZ1+ ~ jZ ; ,

 

 

 

3 х 2

 

2 х 3 Л

 

 

(

 

о

 

2

3 Л

 

 

 

 

 

2 х

 

х

v = 1—-

 

 

/3

 

 

 

 

х

 

 

 

Z 3 +

 

l

 

 

У

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 3 х 2

 

2 х3 Л

5 +

( х 2 х

 

 

3 Л

 

+

2

 

3

z

 

 

г

+ ?

 

 

l

 

 

 

 

 

 

d v

 

6х

 

6х

Z 2 +

 

 

; х

 

+

3 х 2

 

^

+

 

 

 

1 ——

 

z 3 +

l 3

 

2

 

 

 

l

 

 

V

l

 

 

 

 

V

 

 

l 2

 

*

 

*

 

 

 

 

 

 

*

 

6 х

6х

Z 51 +

 

 

2х

+

2

 

+ ------

 

 

—-----

 

l 2

Z 6.

 

l 2

 

l 3

 

V

 

l

 

 

 

 

 

 

 

 

 

 

*

У

 

V

 

 

*

 

 

 

 

 

В ы р а ж е н и я д л я и и v с о в п а д а ю т с з а п и с а н н ы м и р а н е е в р а з д е л е 16.7.

Т аким образом , с пом ощ ью (16.26) и ап п рокси м и рую щ его п о л и ­

н ом а (16.27) п олучен ы точн ы е ф ункции, п озволяю щ и е вы числить

гори зон тальн ое перем ещ ен и е, п рогиб и у го л п оворота лю бого сече­

н и я стерж ня.

 

 

1

1

1

 

1

1

 

1

 

 

 

 

 

5

стерж н я A B

6. Узловым перемещ ениям Z \, Z 2 , Z

3 и Z ; ,

Z

, Z g

соответствую т реакц и и

1

1

R 3 и

1

1

 

 

1

1

R j,

R2 ,

R; ,

R 5 ,

 

R , позволяю щ ие

най ти у си л и я N

, Q и

M

в кон ц евы х сечениях. Д ля определения

и х и сп ользуем ди ф ф ерен ц и альн ы е зависим ости:

 

 

 

N

= du

 

Q

d

v M

d

v

 

 

 

E A

d х

 

E J

d х 3 ,E J

d х 2 .

 

Д и ф

ф ерен ц и руя в вы раж ен и и (16.29) п ервую и третью строки по

одном у

р азу , а вторую - три ж ды , н ай д ем ком п он ен ты вектора

508

к =

 

 

, с п ом ощ ью которого

определяю тся уси -

d х

d х 3 ’

d х 2

 

 

 

 

лия N , Q и M

в пром еж уточн ы х сечен и ях стерж ня:

 

 

 

 

 

к = B H _1 Z ,

 

(16.34)

где

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

B =

 

 

 

6

 

 

 

 

 

 

2

 

 

Д ля определения усилий в

концевы х сечениях стерж ня образуем

матрицу

B а в , первы е три строки которой соответствую т матрице

B

при х = 0 , а остальн ы е - п р и

х = l .

 

 

 

 

 

 

 

1

6

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

2

6/

 

 

Т огда

вектор

у си л и й

S = [ N h , Q h , M h

, N _^ , Q ^ , M K

]

м ож но вы числить с пом ощ ью вы раж ения:

 

 

 

 

 

 

S A B =

B а в H — Z

 

(16.35)

где D - ди агон альн ая м атри ц а ж есткостей:

 

 

 

 

 

D

= d ia g [EA , E J , E J , E A , E J , E J ].

 

 

7.

Н ап равлен и я п олож и тельн ы х у си л и й N , Q и

M в кон ц евы х

сечен и ях

стерж н ей (вектор

S а в ) и п олож и тельн ы х

р еакц и й

R { ,

R 2 , R 3 ,

R 2 , R 5 ,

R 6 (их н ап равлен и я совп ад аю т с н ап равлен и ям и

509

ком п он ен тов вектора S ) не совпадаю т. В заи м освязь м еж ду ни м и

м ож но у стан ови ть с п ом ощ ью м атрицы соответстви я зн ако в у си л и й по вы раж ению :

R ' = И S A B ■.

(16.36)

где И = d ia g [—1,1, —1,1, —1,1].

 

П одставляя в (16.36) вы раж ен и е (16.35), получим :

 

R ' = И D B A B H _1 Z 1.

(16.37)

И з (1 6 .3 7 ) сл е д у е т, ч то м а тр и ц а р е а к ц и й о п р е д е л я е т с я по в ы ­

р аж ен и ю :

 

 

R ' = И

D B АВ H 1

(16.38)

В рассм атриваемом случае

м атрица реакций

(она ж е является и

м атрицей ж есткости) им еет вид, показанны й ранее в формуле (16.25).

В торой сп особ п олучен и я м атри ц ы ж естко сти стерж ня, о сн ован ­

н ы й н а и сп ользован и и п ри н ц и п а Л агранж а, со сто и т в следую щ ем .

П ри известном у р авн ен и и и зогн утой оси (см. р аздел

16.7) стер ж ­

н я м атрицу

ж естко сти его

м ож но п олучи ть из у сл о ви я

стац и он ар ­

н ости п олн ой энергии . П окаж ем это реш ение.

 

 

Запиш ем

вы раж ение п олн ой эн ерги и дл я защ ем лен н ого по к о н ­

цам стерж ня, загруж ен н ого р асп р ед ел ен н о й нагрузкой:

 

 

Э = 1

E J v ft2

- + ■

2 (х ) v =

 

 

 

2

 

 

 

 

2

 

 

 

= I

E J ^ d 2 ( Z 2 /2 + Z 3 /3 + Z 5 /5 + Z 6 / 6 )

 

 

2

 

2

 

 

(16.39)

 

 

 

 

 

+

E A f d ( 1 /1 + Z ; / ; )

 

 

 

2

V

J

 

 

 

 

) ( Z 2 ./2 + Z 3 / 3 + Z 5 / 5 + Z 6 У б)] ,

510