Файл: Д’Анжело, Г. Линейные системы с переменными параметрами. Анализ и синтез-1.pdf

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

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

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

Добавлен: 23.10.2024

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

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

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

Натурные исследования колебания зданий и сооружений [1, 5, 1,12, 119, 122, 139, 140] показывают, что величина коэффициента у в большинстве случаев не превосходит 0,2. При этом коэффициент уменьшения частоты затухающих колебаний

w = V 1 - - V = 1 - 1° ’0 0 5 -

 

Так как в наихудшем

случае

ошибка

не превышает 0,5%,

в

дальнейшем можно принимать

= ш для всех форм колебаний.

Из формулы (11.25)

следует, что

критическое значение

у,

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

будет равно 2. Эта величина также не

зависит от

частоты ко­

лебаний.

 

 

что система

Из формул (11.21), (11.22), (11.25), (11.26) следует,

уравнений (11.14) с матрицей рассеяния

(11.27)

 

_

_і_

 

 

М х + ч(МК)2х + К х = 0

(11.28)

удовлетворяет всем требованиям, сформулированным в § 1 настоя­ щей главы, за исключением свойства устойчивости системы, кото­ рое этими формулами непосредственно не обусловливается. Для доказательства устойчивости можно сослаться на то, что системы уравнений (11.21) или (11.28) в этом смысле не отличаются от часто используемых уравнений колебаний механических систем с вязким трение;'', по гипотезе Фохта, устойчивость которых извест­ на. Независимо ~»т этого факта можно привести доказательство устойчивости на основании формулы (П.23). Передаточная функ­ ция для і-й главной координаты

ф*00

S2 “j“

1

2

имеет полюсы в левой полуплоскости, следовательно, система (II.22) устойчива. Отсюда следует устойчивость системы (11.28), которая получается из (П.22) линейным преобразованием коор­ динат.

При других значениях отношения — , отличных от — = у , решение системы (11.21) будет приводить к различным зависи­

мостям декремента от частоты. При

< -і- декремент колеба-

тп

1

ний будет убывающей, а при— > -^---- возрастающей функцией

частоты высших форм колебаний. Принципиально это дает воз­

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

аппроксимировать

различные заданные (монотонные) зависимости декремента высших форм колебаний от частоты. Выше было показано, что натурные на­

59



блюдения этой зависимости дают большой разброс результатов в- ту и другую сторону и предпосылка о постоянстве декремента для всех форм колебаний была принята в порядке первого приближения, а также для того, чтобы получить результаты, согласованные с вы­ водами хорошо обоснованной гистерезисной теории. Однако экспе­ риментальной базой гистерезисной теории являются испытания объектов другого, более простого вида и вполне возможно, чта дальнейшие исследования приведут к заключению о переменности декремента для высших форм колебаний зданий и сооружений. Систему уравнений (II.1) молено приближенно привести в соответ­ ствие с различными экспериментальными результатами путем вы­ бора различных матриц рассеяния. На основании формул (11.15) и (11.20) общий вид таких матриц будет

 

 

п—т т

 

С = СіМ + С2К + СЪМ ~ К Ч.

(11.29)

Как показано выше, эти матрицы оператором [У?] преобразуют­

ся к диагональному виду.

 

 

 

Рассмотрим в качестве примера матрицу

 

 

_1_ _2_

 

 

 

Сі = сМ 3 К 3 .

 

Система уравнений

(11.21) будет иметь

вид

 

__

2__

 

 

 

іс + с [ М ~ 1 к ) ъ X +

М

~ х К х = 0.

(11.30)

Преобразование координат с оператором [R] приводит

систему

(II.30) к главным координатам.

 

 

 

Уравнение для і-й координаты будет

 

 

 

_4_

 

 

 

 

а{ + счу3 ©г -}-

срг = 0,

 

частота затухающих колебаний —

 

 

 

По сравнению

с формулой (11.12) отношение

возрастает

значительно медленнее. Критическая

частота, при которой перио­

дическое движение в случае данного

с становится невозможным,

составляет сог>-

,

что

значительно

превосходит критическую

частоту по формуле

(11.12).

 

Декремент

колебаний

^принимая'

^ u>£j

Ьі = с™і3

возрастает с частотой, но значительно медленнее, чем согласно вы­ ражению (11.13).

60


Если матрицу рассеяния примем как

С2= с М 3К 3 ,

то получим:

2_

3

Sj = ежа 3

Эти результаты сходны с формулами (П.8) и (II.9), но изменения частоты и декремента происходят значительно медленнее. На рис. 28 показаны зависимости декремента от частоты для матриц и С2 (кривые IV и V соответственно).

§4. С вободны е и вы нуж денны е колебания

Вдальнейшем будем рассматривать систему уравнений (11.28). Если имеется в виду аналитическое решение с помощью разложе­ ния в ряд по формам свободных колебаний, то поступаем следую­ щим образом. Находим обычными методами частоты и формы ко­ лебаний системы без затухания (у= 0); главные координаты будут

иметь вид

(О = aLsin

Главные координаты системы (11.28) получаем умножением

где ги — элементы матрицы собственных векторов. Столбцами

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

бом

по заданным начальным условиям.

задачи о

Аналогично можно искать аналитическое решение

вынужденных колебаниях

 

 

М'х + ч (M K )2 х + К х = f ( t ),

(11.32)

где /

(£) — столбец внешних нагрузок—составляет

 

 

/ ( * ) = [ / і ( * ) . Л ( * ) . ..........> Г п Щ -

 

61

При сейсмических воздействиях внешние нагрузки имеют вид

f ( t ) = [ — m^w (t), т2 w (£),...] = — М w (t ),

где w (£) — ускорение основания сооружения. Уравнение в главных координатах запишется как

У + т 1Р®+

? = R M w (t),

где R — транспонированная матрица собственных векторов мат­

рицы М~1 К. Собственные векторы предполагаются нормирован­ ными по массе:

г=і

При нулевых начальных условиях решение для і-рі коорди­ наты будет (принимается шг = я>г )

Л

і

- 1ШІ ,,

■>

 

о,

f

(11.33)

— “^7

j w(t)e

2

sin шг (t — z) dr,

о

лV

°i= 2 i mk rki

*=i

Координаты x k определяются по формуле

Л

 

2 г« ч*і (*)•

(п-34)

г=і

 

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

баний.

Интеграл (II.33) не может быть найден в аналитическом виде, так как функция w(t) задается в табличной или графической форме. Поэтому для определения реакций на воздействие по за­ кону акселерограммы применяются другие методы решения, не связанные с разложением решения в ряды по формам собствен­ ных колебаний. К их числу относятся некоторые численные мето­ ды, методы теории линейных систем, применяемые при вычислении вероятностных характеристик реакций, и методы моделирования на различных устройствах непрерывного действия. Именно для этих целей и составлены системы уравнений (11.28) и (11.32).

Уравнение затухающих колебаний системы с одной степенью свободы на основании (11.28) будет

т х + іѴ'm k x - \ - k x = 0.

(11.35)

62


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

По форме уравнение (11.35) сходно с левой частью уравнения

т х -\— — х + k x = F(.eia>ft ,

которое применяется для описания циклических процессов сис­ тем с гистерезисным затуханием [43]. Существенная разница за­

ключается в том, что в последнем уравнении сила

неупругого'

сопротивления зависит от частоты внешнего воздействия

, тог­

да как в уравнении (11.35) она зависит от частоты

собственных

колебаний^|//7г k= -£—'J и не связана с внешним воздействием. Эта

обстоятельство позволяет распространить уравнение (11.35) и на непериодические процессы.

Рассмотрим динамические характеристики системы, описывае­

мой уравнением (11.35).

Импульсная переходная функция. Импульсной переходной функцией, или весовой функцией системы называется ее реакция на воздействие мгновенного единичного импульса. Единичный

импульс сообщает массе т скорость г)0 = — . Весовая функция

определяется путем решения уравнения (11.35) с начальными ус­ ловиями

.х(0) = 0 ; * ( 0) = -^ -.

Будем, как и выше, считать приближенно, что частота затухающих колебаний равна частоте колебаний без затухания. Последнюю будем обозначать а>о-

Решение уравнения (11.35) —

х — ае

sin (ш01 — а).

Подставив в это решение начальные условия, найдем весовую функцию для перемещений:

 

 

А. (г) = -------е 2

smcu0zf.

(11.36)

 

 

 

г 4 '

пш0

и

 

Весовая функция для

ускорений —

 

 

Kit)

=

rf2M O

=

I 1-j X")-~ ß

2ш°‘ sin (m0t+ ocj),

 

 

dt2

 

m

 

 

63