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

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

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

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

Добавлен: 02.11.2024

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

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

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

Mu = 4(A„-A i2) = 4 а п—4A,2

M 12= 0 A7o= 0 jV)20

Общая матрица А принимает тогда вид

0

0

0

0 1

0

0

0

0

 

0

0

0

0

0

3

0

3 1 - 1 1

0 0

 

0 0 0

0

0

6

0

6 1 - 1 0 1 0 0 0 0 0

1

 

0

0,5

о 1

0

0

0

0

0

 

0

0

0

0

0 - 1

о 1 0 0 0 0 0 0 0 0

 

0

0 - 1

2 I 0 0 - з

3

 

0 0 0

0

0

0 - 1

4 I 0 0 - 3 0

 

3 0 0

0

ъ

 

0 - 1

8 I 0 0 - 3

0 - 6 9 0 0

0

 

0 - 1

12 |

0

0

- з

0 - 6 0 0

0

1

 

0

—0,5

0 1

0

0

0

0

 

0

0

0

0

0

6

0 - 6 |

0

0

4 0 0 0 0 - 4

0

- з

0 з 1 0 0 0 0 0 0

 

4 - 4

0

 

0

0

0 1

0

0

0

0

0

0

0

0

0

- 1

 

0

- 1 1

0

0 0 0 0 0 0 0

 

0

- 1

0

1 1

0 0 0 0 0

 

0 0 0

33


Здесь каждая из матриц, входящая в А, обведена пунктиром. Правило знаков при формировании матрицы М — обычное, а N>0 при растяжении.

Формируем матрицу L' на основании

 

 

 

 

 

 

 

 

 

 

 

1'м,

 

/.' v. v

 

 

 

 

 

 

 

 

 

 

 

I-м

 

Ls

 

 

 

 

 

 

 

где

L'м и А'л'— матрицы, строки которых состоят

из внут­

ренних усилии в основной системе от единичных сил,

прило­

женных по направлению Лг (рис. 8).

 

 

 

 

 

 

 

 

0

0

0

1

0

0

0

0

0

1

0

0

0

0

0

 

0

3

6

0

0

0 0 0 0 0 - 6 —3 0 —1

 

 

0

0

0

0,5

- 1

- 1

- 1

—1 -

1 - 0 ,5

0

0

0

'0

0

 

0

3

6

0

0

2

4

8 12

0 -Г)

3

0

- 1

1

 

0 —3 —6 0 0 - 1 —2 - 4 - ( i

0 0

0

0

0

0

 

0

0

^ 3

0

0

- 1 - 2 —4 - 6

0

0

0

0

0

0

 

0

0

0

0

0

- 1 - 2 —4 —6 0 0 0 0 0 0

 

0

0

0

0

0

0 _ -1 —3 - 5

0

0

0

0

0

0

 

0

0

0

0

0

0

0 - 2 - 4

0

0

0

0

0

0

 

0

0

0

0

0

0

0

0 - 2

0

0

0

0

0

0

 

0

0

0

0 0

- 1 _ 2 - 4 —6

0

3 0

0

0

0

 

0

0

0

0

0

- 1 - 2 - 4 - 6

0

6

3

0

0

0

Перемножив сформированные матрицы по программе, ука­

занной в Приложении, и приняв

 

 

 

 

 

 

 

 

1

Р12

 

 

 

1

 

 

 

 

 

 

 

 

 

А

------ (с учетом —

перед матрицей В), получим

 

■9б£/

 

J

 

24

 

 

 

 

 

 

 

 

 

>-шах = = 8 2 0 ,2 6 ;

п_ 0,0293 £ /

" к р (тт)

-

Соответствующая матрица-столбец V

0,010

0 ,6 8 7

1,000

0 ,9 3 4 V= 0 ,7 6 6

0 ,2 8 9

0 ,7 3 1

0,026

34


§ 5. Собственные колебания стержнёй

 

 

При динамическом расчете сооружений

важно

знать пе­

риодические колебания — собственные и вынужденные.

Соб­

ственные колебания наблюдаются при отсутствии

нагрузки,

вызвавшей их, а вынужденные — когда вызвавшая

их

сила

продолжает возмущать систему.

определяются час­

В динамическом расчете прежде всего

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

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

апримеры рассматривать с малым числом степеней свободы. Задача собственных колебаний, согласно принципу Далам-

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

// = wS m iVi >

где уг— перемещения масс, определяющие главную форму колебаний, то есть такую, при которой все массы колеблются с: одной частотой со. Эти перемещения определяются из систе­ мы уравнений

У\^ б п (со2mxyi )+ 612(со2т2у2) +...+ */2=621 (со2mitjx) +622 (<о2т2у2) +-..+

Уп = Кп iv-щуд + v (‘^woys) +■...

Обозначим в этих равенствах

V,

 

611

612

...

г‘\п

II

621

622

•••

г‘->п

.

1 и"тпуп)

 

ф»„ (c"'-w„y„)

(а)

\ п ^

1т пУп) ■

 

т1

0

 

пи

 

1

 

 

оз-

°п\

°л2 ••• г‘пп

тп

Уп

0

36


Записав систему уравнений (а) в матричной форме

 

(DM—KE) К=0 ,

(15)

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

В практических расчетах обычно отыскивается лишь на­ именьшая собственная частота___

ч - / ^

не»

I^ШОХ

исоответствующая основная форма колебаний после решения системы уравнений

(D M-).msxE ) Y = 0.

(17)

Как и при расчете на устойчивость перемещения сосредото­ ченных масс Yt определяются с точностью до постоянного множителя.

Таким образом, задача о собственных колебаниях статиче­ ски определимого сооружения свелась к решению уравнений (15), (16) и (17). Входящая в эти уравнения матрица D яв­ ляется матрицей единичных перемещений определяемых по известному соотношению

D = //B L ,

где матрица L формируется па основании эпюр, построенных от единичных сил, приложенных в точках приложения сосре­ доточенных масс.

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

доточенными массами (рис. 9, о). Приведем

балку к безраз­

мерному виду, приняв I

 

т!_

-характерный размер, М= —----ха-

рактерная масса, где

 

3

т -- масса, приходящаяся на единицу

длины балки, EI — характерная жесткость (рис. 9,6). Расчет­

ные участки и сечения показаны ,на рис. 9, в.

Матрица подат­

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

в ,

 

1 2

 

1

L . . I .

2

1

В I V - 6 - 6 1 1

 

2

I

36

1 2

Bi I

в ,

1

2

1

 

1

2

1

3-6

1

2

=

18

1

о

 

 

37


Матрицу податливости В формируем с учетом сдвижки

2

1

2

 

0

1

6

2

 

 

2

8

1

 

 

2

6

0

 

 

1

2

Матрицу влияния моментов формируем на основе эпюр рис. 9, г

0 0 0

5

3

1

О

3

О

9

1

3

5

0

0

0

Матрица сосредоточенных масс

1

38

Введем в оперативную память ЭВЛ\ «Напри» матрицы /., В и М . Проделав над ними операции по программе, изложен­ ной в Приложении, получим /.тах=720 и после решения урав­ нения (17) получим ординаты соответствующей формы коле­

бания (приняв l/i =

l)

 

 

 

 

 

 

 

У=

 

 

 

 

Т о г д а н а и м е н ь ш а я с о б с т в е н н а я ч а с т о т а

 

(с у ч е т о м с к а л я ­

ров, В ы н е с е н н ы х за

з н а к и м а т р и ц L ,

В )

 

 

Ж'Е!

/

36

9,86

 

/'/,/

!

ml

 

Г-Ш

 

720 I '

 

 

 

/

V

т

В дальнейшем,

приведяJсооружение к безразмерному ви­

ду, будем наименьшую собственную

частоту

определять как

0)min

 

Е!

 

 

(18)

РМ /-шах

 

 

 

 

 

сооружений для

При расчете статически неопределимых

определения частот и форм колебаний можно

воспользовать­

ся уравнением (12), но матрица D должна формироваться для статически неопределимой системы по формуле

D = L 0B M ,

где L 0 — матрица влияния моментов от действия единичных сил в основной системе метода сил, приложенных в точках приложения сосредоточенных масс;

М — матрица окончательных моментов* определяемая по формуле метода сил

М = М р - L ( L ’B L ) - xL ' B M p ,

зная, что

M

p = L 0,

 

M = L

0 - L ( L ' B L ) - l L ' B L 0 .

Матрица М

может быть получена и методом перемещений

по формуле

М = М Л — L R ~ XR P ,

здесь матрицы М р , L , R и R p формируются для кинемати­ чески определимой основной системы.

Однако для многих статически неопределимых задач полу­ чение матрицы D для конструкции, загруженной инерционны­ ми силами, может оказаться весьма трудоемкой задачей. Воз­ можный путь упрощения этой задачи указан в работе [4].

39