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

Литература / А.А. Борисевич. Строительная механика. Минск, 2009

.pdf
Скачиваний:
17
Добавлен:
23.01.2020
Размер:
10.86 Mб
Скачать

M Y = X TM a ,

(20.29)

что существенно уменьшает объем вычислений, когда матрица масс M диагональная.

Продифференцировав по времени выражение (20.25) и подставив

в него начальное время to = 0

и начальные скорости из начальных

условий (20.9), получим еще одну систему уравнений:

 

У

 

 

 

b = Z [X i^i cos(n )]ui = Z

X iyi = X V ,

 

 

 

 

(2°.30)

 

 

 

i=1

 

 

i=1

 

Т

 

 

 

 

 

 

 

 

 

где V - вектор новых переменных, компоненты которого равны:

 

 

 

 

 

 

 

 

Н

 

 

 

 

vi = uimi co s(n ).

 

Б

(20.31)

 

 

 

 

 

 

 

 

 

 

Вектор неизвестных V найдем

 

з с стемы совместных уравне­

ний (20.30) так же, как и (20.26), п

еобйазованных к виду с диаго­

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

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R V = X TRb

или M V = X TM b .

(20.32)

 

 

 

р

 

 

 

 

После вычислен я векоровY и

V параметрыut и rji, входя­

щие в формулу (20.25),топределяются покомпонентно из соотноше­

ний (20.27) и (20.31):

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

з

 

 

 

 

L.

(20.33)

 

ui = VУ2 + (yi , a i )2 ;tg (Г ) =

 

 

о

 

 

 

 

 

yi

 

 

 

 

 

 

 

 

 

 

 

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

п

 

 

 

 

 

 

 

 

 

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

ний, сведена к алгебраической задаче - полной обобщенной про­

е

 

 

 

 

 

 

 

 

 

 

блеме собственных значений для двух матриц. В настоящее время

разработаны эффективные, надежные алгоритмы и компьютерные

Рпрограммы для решения полной (и обобщенной) проблемы собст­

венных значений. Однако следует помнить следующее. Для реше-

611

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

собственных векторов, а также и для ряда рабочих массивов. Поэтому

даже современные компьютеры успешно решают полную проблему собст­

венных значений для систем с относительно небольшим числом степеней

свободы, порядка несколькихсотен или максимум несколькихтысяч.

 

 

 

 

 

 

Т

 

Примеры динамического расчета реальных сооружений показы­

вают, что в суммах вида (20.25) слагаемые, отвечающие высшим

 

 

 

 

 

Н

собственным частотам, получают значения, пренебрежимо малыеУ

по сравнению с вычисляемой суммой. Практически для получения

 

 

 

 

 

Б

 

искомого результата можно ограничиться вычислением суммы не­

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

ным частотам. Поэтому отпадает необходимость в решении полной

 

 

 

 

й

 

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

точно решить частичную проблему собственных значений, т. е. найти

 

 

 

Или

 

 

несколько низших собственных частот

соответствующих собствен­

ных форм (собственных векторов).

 

же на ти собственные час­

 

 

р

 

 

 

тоты и собственные формы, отвечающ е заданному диапазону час­

 

 

о

 

 

 

 

тот. Для этого применяют специальные ч сленные методы. Некото­

рые из них будут рассм трены ниже.

 

 

 

 

П р и м е р

20.1. Определи ь с бственные частоты и собственные

формы свободных колебаний неразрезной двухпролетной консоль­

ной балки постоянного сечения, несущей сосредоточенные массы

(рис. 20.2). Собственнойтраспределенной массой балки по сравне­

нию с сосредоточенными массами пренебречь.

 

 

 

При с вершенииипоперечных, изгибных колебаний балка, несущая

четыре с средзт ченные массы, имеет четыре степени свободы, направ­

ления к т рых б значены стрелками F\-F4 (рис. 20.2). Чтобы сокра­

тыре

 

 

 

 

 

тить вычисления,оиспользуем симметрию системы. Для этого введем

ч

гру

овые степени свободы (рис. 20.3): две прямосимметрич­

ныепо(ПС)

направлениям F\ и F2 и две кососимметричные (КС) по

направл ниям F3 и F4.

 

 

 

 

 

При прямосимметричных деформациях сечение балки над цен­

тральной опорой не поворачивается, а при кососимметричных де­

Рформациях в этом сечении изгибающий момент равен нулю. Следо­

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

612

таточно рассмотреть ее половину, введя в центральном сечении за­ щемляющую опору (рис. 20.4). Кососимметричные колебания бу­ дем изучать с помощью полубалки, полученной из исходной балки

введением шарнира в центральном сечении (рис. 20.5).

 

 

 

 

Частотные (вековые) уравнения для полубалок составим в об­

ратной форме (20.15).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Построим эпюры изгибающих моментов от единичных сил, при­

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т

ложенных в направлении степеней свободы полубалок (рис. 20.4,

20.5), и вычислим по формуле Мора путем «перемножения» эпюр

компоненты матриц податливости полубалок.

 

 

 

 

 

У

 

Для прямосимметричных колебаний имеем:

Б

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

4

 

 

 

 

 

 

a 3

 

6

 

6I

Н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D nc =

 

 

 

 

 

 

 

 

й

7

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S 22 _ = 1 2 E J - 6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

Fi

 

 

 

 

 

 

 

р

F

 

 

 

 

F 4

 

 

 

 

 

 

 

 

о

 

 

4m

 

 

 

 

m

 

 

 

 

 

a

 

 

 

a

co nst

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

E J =

 

 

 

 

 

 

 

 

 

 

 

 

 

2

т

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

о

и

 

 

 

Рис. 20.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b______________ u

 

 

п

з

 

F 3

 

 

 

 

 

F 3

 

 

 

 

 

 

 

е

 

 

 

 

f 2

 

 

 

 

 

f 2

 

 

Д

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р

 

F i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fi

 

 

 

m V

 

 

4

4m

 

 

 

T

 

4m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

;

 

 

 

V /////,

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 20.3

 

 

 

 

 

 

 

 

613

 

 

 

 

 

 

 

 

 

 

Н

У

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 20.4

 

 

Б

Т

 

 

 

 

 

 

 

й

 

 

 

 

 

 

 

 

 

и

 

M3

 

 

 

 

 

 

 

р

 

 

F4 = 1

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

M4

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 20.5

 

 

 

 

 

 

Соответственно для косос

мметричных колебаний получим:

 

 

 

 

1

 

1

 

1

 

 

1

 

 

 

 

 

 

 

4

 

3

 

 

 

 

и33

34

= я 3

-

 

 

 

 

DJKC =

 

 

3EJ

- 3

 

5

 

 

 

 

з

_^43

 

1 4

 

 

 

е

 

 

 

 

 

 

 

 

 

 

 

Р

Матрицыомасс полубалок примут вид:

 

 

 

 

 

п

 

"1

0"

 

 

 

 

"4 0"

 

 

 

 

 

m

 

;

 

 

= m 0 1

 

 

 

 

M ПС = m 0 4

M KC

 

 

Частотное уравнение (20.15) для первой полубалки, соответствую­ щей прямосимметричным колебаниям, можно представить в виде:

614

 

 

 

 

 

 

16 - р

 

-

24

= 0

 

 

 

 

 

 

 

 

 

 

 

 

- 6

 

28 - р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

12EI

 

12EI

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

р=

 

3

Л = ~

 

3.

 

 

 

 

 

 

 

 

 

 

ma

 

 

со ma

 

 

 

 

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Раскрыв полученный определитель и решив квадратное уравне­

ние относительно параметра р , найдем, что:

 

Б

У

 

 

Р,2 = 22 ± V180 ; или

 

 

 

 

 

 

 

 

р = 35,4; Р2

= 8,58.

 

 

Откуда находим выражения для круговых частот Нпрямосиммет­

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

 

и

 

 

 

 

о>1 =

12EJ

= 0,582

E J

 

 

12EJ

= 1,182. E J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ний) найдем из сис емы дн рдных уравнений (20.14), которая в

 

 

q\ma~

 

 

ma

 

 

0 2

=йqpma~3

 

ma 3

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

Соотношения между амплитудами перемещений масс при пря­

 

 

 

 

 

т

 

 

мы п ямосимметричных колеба­

мосимметричных колебаниях (ф

 

данном случае

примет

 

 

 

 

 

 

 

 

 

 

 

 

вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

з

- р

 

-

24

 

 

 

 

 

 

 

 

о

16

 

*1

= 0.

 

 

 

 

 

-

6

28 - р

X 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученная система однородных уравнений для найденных зна­

чений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р

р является вырожденной и эквивалентна одному уравнению:

п

 

 

(16 - р ) Х 1 -

24X 2 = 0,

 

 

 

 

 

 

 

 

 

откуда получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X 2

= (16-р)Х 1

 

 

 

 

 

 

 

 

 

 

 

=

 

 

24

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

615

Итак, первая собственная форма прямосимметричных свободных колебаний при р= р1 определяется соотношением:

X 2 = -0,809Xi.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У

Вторая собственная форма прямосимметричных свободных ко­

лебаний при р = р2 определяется соотношением:

 

Т

 

 

 

 

X 2 = 0,309X 1.

 

 

 

 

 

 

 

 

 

 

 

 

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

собственные

векторы

 

прямосимметричных

колебаний для полубалки примут вид:

 

Б

 

 

 

 

 

"

1

"

 

 

 

"

1

" Н

 

 

X (1) =

 

 

 

 

 

 

й

 

 

 

 

-

0,809

;

 

 

 

X (2) =

0,309

 

 

 

Для балки в целом будем иметь:

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

"

1

"

 

 

 

 

"

1

"

 

 

 

 

 

-

о

 

 

 

0,309

 

 

 

 

 

 

0,809

р; X (2) =

 

 

 

 

X (1) =

т

 

 

 

 

 

 

-

0,809

0,309

 

 

 

 

 

и

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

Частотное (вековое) уравнение (20.15) для второй полубалки, со­

о

 

 

 

 

 

 

 

 

 

 

 

 

 

ответствующей косос мметричным колебаниям, примет вид:

 

п

з

 

16 - р

 

- 3

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

Р

 

 

 

-1 2

 

5 - р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда найдем:

 

 

 

 

 

 

 

 

 

 

 

 

 

ер 4 = 10,5 ±-у]66,3 , или

 

 

р

= 18,64,р4 = 2,61.

 

 

При кососимметричных колебаниях принято, что:

616

 

 

 

 

 

 

 

3EJ

3EJ

 

 

 

 

 

 

 

 

 

 

 

<р = -

А -

~

ma

3 '

 

 

 

 

 

 

 

 

 

 

 

 

ma

 

о

 

 

 

 

 

 

Следовательно, круговые частоты кососимметричных свободных

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

 

 

 

 

 

 

 

У

 

3EJ

 

 

 

 

E J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т

о = V Р ma

3

= 0,401V ma

 

Л

 

 

3EJ

= 1,127

, EJ

 

 

= VP4mar3

 

ma

3

 

 

 

 

 

 

 

 

 

 

 

 

 

Н

 

 

Соответствующие собственные формы кососимметричных сво­

бодных колебаний найдем из соотношения:

Б

 

 

 

 

 

 

 

 

 

 

 

 

 

(16- p ) X 3 - 3 X 4 = 0,

 

 

 

 

 

откуда при р —р

имеем

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

р

 

 

 

 

 

 

 

 

 

 

 

 

 

X 4 = -0,880Xй3 ,

 

 

 

 

а при р —Р4

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

 

 

 

и

X 4 = 4,55X 3.

 

 

 

 

 

 

Для балки в целом собс венные формы кососимметричных коле­

 

 

з

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

"

0,880

"

 

 

 

- 4,55'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

х (3) =

 

- 1

;

X (4) =

-1

 

 

 

 

 

 

 

1

1

 

 

 

 

е

о

 

 

 

 

 

 

 

 

 

 

 

 

-

0,880

 

 

 

 

4,55

 

 

 

 

 

 

 

 

 

 

 

 

 

 

РВ динамике сооружений принято найденные собственные часто­ ты располагать в порядке возрастания (строго говоря, в порядке не­ убывания, так как среди найденных частот могут быть одинаковые)' ассмотренная симметричная двухпролетная балка с консолями, имеющая четыре степени свободы, характеризуется спектром из четырех круговых собственных частот и четырех соответствующих

617

 

 

 

 

 

 

 

 

 

 

 

 

У

 

 

 

 

 

 

 

 

 

 

 

Т

 

 

 

 

 

 

 

 

 

 

Н

 

 

 

 

 

 

 

 

 

 

 

Б

 

 

 

 

 

 

 

 

 

 

 

й

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

р

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

з

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

е

 

 

 

 

 

 

 

 

 

 

 

Р

 

 

 

 

 

 

 

 

 

 

 

 

20.3. Действие вибрационной нагрузки при отсутствии сил сопротивления

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

 

 

 

 

 

 

Н

описывать неоднородными дифференциальными уравнениями дви­У

жения (20.5), составленными в прямой форме. После подстановки в

 

 

 

 

 

 

Б

 

их правую часть вектора внешних вибрационных сил (18.2)Туравне­

ния движения примут вид:

 

 

 

 

 

 

 

M Z + R Z —F sin(# t).

 

(20.34)

 

 

 

 

 

и

 

 

 

Как известно из теории дифференциальных уравнений, общее

решение неоднородных линейных д фференциальных уравнений

можно представить как сумму общего йешения соответствующих

 

 

 

 

о

 

 

 

однородных уравнений вида (20.4) частного решения неоднород­

ных уравнений, в данн м случае (20.34). Общее решение однород­

 

 

 

 

т

 

 

 

ных уравнений движения, с держащеепостоянные интегрирования,

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

ния деформируемой с с емы. В реальных условиях из-за неизбеж­

 

 

 

з

 

 

 

 

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

ния с течением времени затухают, и деформируемая система про­

 

 

совершать

 

 

 

 

должает

 

чисто вынужденные колебания с частотой вы­

нуждающей нагру ки в так называемом установившемся режиме.

 

п

 

 

 

 

 

 

При к лебаниях в установившемся режиме перемещения всех

е

 

 

 

 

 

 

узлов деф рмируемой системы, усилия во всех ее элементах с тече­

ни м вр мени изменяются по гармоническому закону с той же час­

Р

 

 

 

 

 

 

 

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

Определение экстремальных усилий и экстремальных переме­ щений при действии динамических нагрузок составляет предмет

619

второй основной задачи динамики сооружений, известной под на­ званием динамический расчет сооружений.

Итак, будем искать частное решение системы дифференциальных

уравнений движения при вибрационном возбуждении в виде одночас­

тотных вынужденных колебаний, когда неизвестные динамические

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У

перемещения Z (t) изменяются по гармоническому закону с частотой

 

 

 

 

 

 

 

 

 

 

 

Т

вибрационной нагрузки 0 , но с неизвестными амплитудами Z

:

 

 

 

 

 

Z (t) = Z sin(0t).

 

Н

 

 

 

 

 

 

 

 

 

 

 

Б

 

 

 

 

Дифференцируя (20.35) дважды по времени, получим выражение

для вычисления ускорений:

 

 

 

й

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(20.36)

 

 

 

Z (t) = - 0 2Z sin(0

1).

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

Подставив динамические перемещен я (20.35) и ускорения (20.36)

 

 

 

 

 

 

р

 

я (20.34),получим после

в дифференциальные уравнения дв жен

сокращения на sin(0 t)систему л нейных алгебраических

 

уравне­

 

 

 

 

о

 

 

 

 

 

 

 

 

ний относительно амплитуд динамических перемещений:

 

 

 

 

 

 

т

 

 

 

 

 

 

 

(20.37)

 

 

 

(R

- 0 2М )Z = F .

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

з

 

 

 

 

 

 

 

 

 

 

 

 

Матрицу коэфф ц ен ов сис емы алгебраических уравнений (20.37)

обозначают как

 

 

 

 

 

 

 

 

 

 

 

 

 

е

о

 

 

R D = R - 0 2М

 

 

 

 

(20.38)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и называют матрицей динамической жесткости (точнее, внешней ди­

намич ской жесткости) деформируемой системы. Система алгебраи­

Р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ч скихпуравнений относительно амплитуд динамических перемеще­

ний, выраженная через матрицу динамической жесткости, примет вид:

 

 

 

 

 

R D Z = F .

 

 

 

(20.39)

620