Файл: Д’Анжело, Г. Линейные системы с переменными параметрами. Анализ и синтез-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