Файл: Васильев, С. П. Приближенные методы расчета сооружений на устойчивость и динамику с применением ЭВМ Наири (учебное пособие).pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 02.11.2024

Просмотров: 28

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Рис. 12

Основная система метода перемещений показана на рис. 12,6. При расчете симметричной формы потери устойчи­ вости выберем 6 параметров (ZrrZ6), учитывающих повороты узлов и 6 параметров (Z7-rZ]2), учитывающих изгибные де­ формации сжатых стоек рамы. Соответствующие эпюры от единичных Zr :-Z12 построены на рис. 12, в. По табл. 2 находим единичные реакции:

48

 

16£/

И2

PI

Гчч= ^

_

1 « Р / .

 

 

 

 

/

 

420

 

 

/

 

420

 

 

 

 

r 33 — /'бб

 

24£/

 

Ш Р1

 

40£7

1008

Я/;

 

 

 

 

 

/

~~ 420

 

 

/

 

420

 

 

 

 

.48£ /

448

 

 

 

99 ~

г 12,12'

409,6 £ /

4096

_Р_

г-,.-,

/

 

420 Pt ;

 

/3

~

420

/

 

 

 

 

 

 

 

 

 

 

 

 

 

5*«8

5;i, и

409,6£/

12288Я

 

 

 

 

 

=

 

Р

 

4201

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5 10, 10

429,6 £ /

 

, 24576Р

 

 

 

4 £ /

56

 

 

£

 

 

420/

 

/*12

/*35=

420 Р1\

 

 

 

 

 

 

Г27

_

448 п

r^ ~

 

 

8£/

 

 

 

Г17

 

42()

r^ — f56------

 

 

 

 

\Е 1

 

168

п ,

Гой-

 

 

 

1344Р

 

 

Гла

-|- "

/ / I

-548= 55, И =

420/

 

 

24

/ т 420

 

 

 

 

/ 39 — 559 — 56, 12

 

448Я

54, Ю = —2688Р

 

 

 

 

 

 

 

420/

 

 

■420/

 

 

Остальные r i k =

0. Обозначив а =

Р12

 

 

 

и учитывая, что г ш =

= гк1 , сформируем матрицы ^ и S, необходимые для ввода з • оперативную память ЭВМ,

Матрица R принимает вид

16

4

 

4

 

 

4

40

8

4

 

 

8

24

0

 

 

4

0

40

8

 

 

 

 

8

48

8

 

 

 

 

8

24

409,6

409,6

409,6

409,6

409,6

409,6

Матрица поправок S формируется с учётом того, что в формуле (19) эта матрица вычитается из R и ее элементы бе­ рутся со знаками, противоположными данным табл. 2.

Введя эти матрицы в память

ЭВМ «Наири»,

получим

1

E I

amin = г— =0,0139 и соответственно

Якр(тт) =5,83 — ;

^тах

*

 

49


 

112

- 5 6

0

0

0

0

448

0

0

0

0

0

 

- -56

448

0

—168

0

0 -448

1344

0

0

0

0

 

0

0 112

0

0

0

0

0

448

0

0

0

 

0 -1 6 8

0

1008

0

0

0 -1344

0

2688

0

0

 

0

0

0

0

448

0

0

0

448

0

1344

0

5 =

0

0

0

0

0

112

0

0

0

0

0

448

418

-4 4 8

0

0

0

0

4096

0

0

0

0

0

 

0

1344

0 -

-1344

0

0

0

12288

0

0

0

0

 

0

0 448

0

448

0

0

0

4096

0

0

0

 

0

0

0

2688

0

0

0

0

0 24576

0

0

 

0

0

0

0

1344

0

0

0

0

0

12288

0

 

0

0

0

0

0 448

0

0

0

0

0

4096

^ 77 А W

^ Т,

L i

i s

S Г 777 7777) 77. 7

* ГГ, 7 77

 

£Г J *7 77 г

Рис. 13



Для обратно симметричной формы потери устойчивости схема рамы, нагрузки, основная система и единичные эпюры приведены на рис. 13. Проведя решение аналогично предыду-

EI

щему, получим amin =0,00629 и P Kp(min) =2,64 — . Как и сле­

довало ожидать, при обратно симметричной форме потерн устойчивости для данной рамы критическая сила получилась меньше, что свидетельствует о ‘большей вероятности именно такой формы потери устойчивости. Как показано в работе [81,

если пренебречь влиянием изгибных

деформаций стоек, то

есть принять Zio=Zii = ... = Zi5= 0, то

получим amin =0,0063,

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

чится только линейными перемещениями узлов

(Z\,

Z2 и Z3

не равны нулю), то погрешность при определении

аШт

дости­

гает 15%, что недопустимо.

 

 

§ 7. Пример расчета рамы на динамику

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

Определить низшую частоту собственных колебаний рамы, показанной на рис. 14, а. Наименьшую частоту колебаний для данной симметричной рамы можно ожидать только при обрат­ но симметричной форме колебаний, поэтому основную систему и обобщенные координаты выбираем с учетом этой особенно­ сти (рис. 14, б, в). Пользуясь табл. 3, строим эпюры Мх и М2 (рис. 14, г) и находим единичные реакции:

.,

\

. ^ « 7 4E I

■—

т

йз

-|-

4 £ 7 - 3

 

 

3 m « > W - 2 »

 

2 £ / - 3

тi

21

h

 

105

——

105 - З3

 

 

{

 

 

 

 

 

 

 

 

1

. _,,ж5

 

 

 

 

3т Д-/;3• 2 :!

\

J

26/:/

 

2 2 т о>2Л3

 

 

 

 

 

 

~1”

 

140-33

 

=

h _

 

 

945

 

 

 

 

Г12---- 2 (\

6 Е1

11 т со3 лз

 

1 2 £ /

1 1 т о,2/;3

 

 

/!2

210

 

 

Л

 

+-

210

 

 

 

 

 

 

 

 

13т о,2 h

 

 

Q

2 h

 

 

24Е1

96

т ш2h .

 

 

 

 

 

 

 

36

 

 

от

иг — ==

 

hз

35

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

51


я)

Jm

к)

 

 

 

-------- 7

 

т '

ОсноЗная

 

 

 

система

 

 

E h const

ГГТУ

гг

 

*■

 

Ш к

 

 

2 h

Здесь слагаемое Зт со2 — учитывает реакцию во введенном

опорном стержне от силы инерции, возникающей при колебанин ригеля в горизонтальном направлении.

п

^

 

„ ' •

гп

получаем

уравнение собст-

Введя

обозначение

 

 

веппой частоты

 

 

—124 0,105Я.

I

Q

 

 

26—0,023?.

 

 

 

-12+ 0,105л

 

24—2,743?,

|

 

$2

Наименьший корень этого уравнения Опт =6,95, откуда наинизшая частота

 

 

 

 

u)min

2>64 1 Г ЕГ

 

 

 

 

 

 

 

 

 

 

к2

г

гп

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как и при решении задач устойчивости,

запишем полные

уравнения собственных колебаний:

 

 

 

 

 

 

 

 

(26—0,023a)Z1+ (—12+0,105X)Z2 = 0 ;

 

 

 

(-12+0,105?.)Zi+ (24-2,7431)Z2=0.

 

<а)

В матричной форме

 

(Д—со2F)Z = 0,

 

 

 

 

(22)

где

 

 

 

 

 

 

 

26

- 1 2

 

F= 1 0,023

-0,1051

 

 

[Z,

 

,

,

Z =

 

-12

24

 

 

 

! -0,105

 

2,7431

 

 

' z 2

F — матрица динамических поправок к реакциям

в связях от

стержцей, подверженных колебаниям (см. табл. 3).

уравне­

Для приведения уравнения (22) к виду векового

ния

\А—A£ j =0 умножим его слева на R~l

и

разделим на

(—и2), получим

(Д -1/7- — F)Z=0,

 

 

 

 

(23)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О)^

 

 

 

 

 

 

 

где

A = R~1F и

 

 

Решив

уравнение

(21),

определим

Amin

и тогда, с учетом вынесенных за знаки матриц скаляров,

получим

 

 

1

 

 

 

 

 

 

 

 

 

(24)

 

 

 

 

 

 

Лтах

 

 

 

 

 

 

 

Ч’пип

I 2 V

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

Таким образом, решение задачи об определении наинизшей

частоты собственных колебаний рам сводится

к формирова­

нию матриц R и F, решению с ними векового уравнения (23),

а затем к определению по формуле (24)

ujmin .

 

 

 

 

Пример*. Определить паинизшую частоту собственных ко­

лебаний для рамы,

показанной на рис.

15, а.

Массы стержней

постоянны и распределены по их длине равномерно.

Так как

рама и распределение масс симметрично, то необходимо рас­ смотреть симметричную й обратно симметричную форму коле­ баний рамы. Рассмотрим подробнее симметричную форму ко­ лебаний. Охарактеризуем эту' форму колебаний двенадцатью обобщенными координатами ('Zi+Z6 учитывают повороты уз­ лов, Z7+Z!2— изгибные деформации стоек рамы). Основная

* Пример-взят из работы [81-

53