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

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

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

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

Добавлен: 23.10.2024

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

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

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

где

7

tg<4 =

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

* w = ^

в 2

sinto0 ( t - x ) d * .

о

При единичном гармоническом внешнем воздействии с часто­ той р

уравнение

-ро X + Шц X = eipt ,

(11.37)

X

левая часть которого

получается

из (11.35) делением

на т, при

нулевых начальных условиях имеет решение

 

 

-*•

2

еірі

 

 

J I .

 

 

 

“о —

Р + Ч шйР

 

Частотная характеристика системы —

 

Ф(ір),

 

1

(11.38)

передаточная функция—

“о — Р°"+ г>0Р

 

 

1

 

Ф(з) =

 

(11.39)

s"-f- s

+ ы0

 

 

 

Полюсы передаточной функции

Sl,2-----Н Г ± * “<>] / ' 1 _ 'Т

лежат в левой полуплоскости.

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

Резонанс. При резонансе амплитуда колебаний, вызванных •единичным внешним воздействием, равна

а = _j__ Пшо

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

64


■ * + ü>o(l +it)x=elPl,

( 11.40;

однако соответствующие им частотные характеристики различны. Форма резонансной кривой определяется величиной модуля час­ тотной характеристики (11.38)

\Ф{ір)\=

1

V

(ши - р-уі + f шір -

Для уравнения (11.40)

 

(ір) I =

|/~ (_“иР І У +7' “и

Сдвиг фаз между внешним воздействием и реакцией системы в первом случае

о = arctg

(Im Ф{ір)

= arctg j -

 

Re Ф {ip)

во втором случае

 

 

шІ - р - /

 

 

 

cp =

arctg

ü)- — p~

 

 

 

При низких значениях у разница в форме резонансной кривой и величине фазы реакции мала и обычными измерениями не может быть обнаружена.

По гипотезе вязкого трения (гипотеза Фохта), выражаемой уравнением

х -]—2s

CDX ~ 01

частотная характеристика будет

 

Ф (ір) =

1______

“и — р'2+2івр

 

Это выражение будет совпадать с (II.38), если определять коэф­ фициент е по заданному декременту колебаний:

£оцо _ 7шо 2тс — ~~2~‘

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

Диссипативная функция. Матрице рассеяния энергии

С = т М * К *

соответствует диссипативная функция

П

 

 

 

Vт

ік

X.

(11.41)

г,а— 1

 

 

■5-248

65


1

 

 

 

где &>Ä2 ' — элемент матрицы К~.

 

 

 

Это выражение можно записать так:

 

 

1

j __

_

_

Z7 = -1- т х* Ж 2 /с 2 х = 4 - т

 

W x ,

где X* — вектор-строка, а х — вектор-столбец.

 

__ і_

j_

 

 

Элементы матрицы IF = Al 2 ЛГ2 имеют размерность частоты.. Для системы с одной степенью свободы

с-

 

1

 

-2

г

= ~Y

т и)0 X .

Следовательно, скорость

рассеяния

энергии пропорциональ­

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

Диссипативная

функция по гипотезе Фохта равна

 

 

 

 

 

1

П

 

 

 

 

С.

V

 

 

 

 

F = —

Z

СІ^І <?*•

 

 

 

 

 

 

і,

k= 1

 

 

Здесь

коэффициенты матрицы

рассеяния сосчитаются

постоян­

ными, зависящими

от свойств

материала. Поэтому скорость

рас­

сеяния

энергии и силы сопротивления непосредственно

не

зави­

сит от механических характеристик системы, таких как вид опорных закреплений, распределение масс и жесткостей. Функ­ ция F постоянно встречается в литературе, где излагается прин­ ципиальная сторона вопроса, однако не имеется никаких данных о методах определения и величине элементов матрицы рассеяния ||сіИ| применительно к строительным конструкциям и сооруже­

ниям. Этот вопрос полностью разрешен при применении матриц вида (11.27)! В сложных матрицах вида (11.29) соотношение между коэффициентами ск выбирается произвольно, с таким расчетом, чтобы получить удовлетворительную аппроксимацию заданного закона зависимости декремента колебаний от частоты высших форм. Величины коэффициентов ck определяются в соответствии

с величиной декремента основного тома колебаний. Тем самым метод вычисления всех элементов матрицы рассеяния становится достаточно определенным.

Рассмотрим пример составления и решения уравнений сво­ бодных колебаний системы с двумя степенями свободы. Объектом, исследования примем абстрактную схему с параметрами

66


2 - 1

Матрица жесткости соответствует схеме со сдвиговыми деформа­ циями при одинаковой жесткости обоих этажей. Коэффициент у принимаем 0,06, что приблизительно соответствует 6=0,2. Ввиду того что матрица жесткости симметрична, извлечение квадрат­ ного корня приводит к системе трех квадратных уравнений

^

-Хо\

/

-^*2 )

і 2

'—1

\Х2

x j

\ х 2

x j

\ —\

1(

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

хх х 2+ х 2 х 3= — 1; JCj+Jt* = 1.

Решая уравнения, находим

 

 

 

К 2 =

1

 

 

 

 

 

 

 

 

 

 

 

5

Ѵ 1

 

 

 

 

 

 

 

 

 

у

 

 

 

 

 

Знак

перед

матрицей должен быть положительным, так как глав-

ные

члены

матрицы

-L _L

 

 

 

 

 

 

 

 

 

М 2 К 2 положительны. Система уравнений

(11.28)

будет иметь вид

 

 

 

 

 

 

 

 

 

 

 

х і Н—

Хо) + 2 x t х 2= 0;

 

 

,

У 5

 

 

 

 

 

 

 

 

 

 

 

 

0,06

(2-х:, — X J

N

+ х

2Хі =

0.

 

 

х 2+

- у —

 

Оператор системы —

 

 

 

 

 

 

 

 

 

 

 

 

 

,

, 0,18

. 0

0,06

 

 

.

 

 

 

р- + -^=г-р +2

/ 5

 

р —1

 

 

L =

 

ѵт

.

 

 

 

 

 

 

 

 

 

0,06

 

 

0,12

р +1

 

 

 

 

7=ГР-1

 

 

■ ’

5

 

 

 

 

У 5

 

 

 

V

 

 

 

уравнение частот —

 

 

 

 

 

 

 

 

 

 

Х2 +

-^ ü -X + 2

 

АЩ-Х + 1

 

 

0,06

 

 

 

 

 

—_ /. +1 = 0 . (11.42)

 

V 5

 

 

К 5

^

 

 

 

У ь

 

 

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

 

 

Х4 +

0,134X3-}- 3,003X2+

0,134/, +

1 = 0 .

67


Вместо решения этого уравнения можно рассмотреть более прос­ тую задачу о незатухающих колебаниях:

 

 

 

х [ -f- 2 X у— -C-j — О,

 

 

 

Уравнение частот

 

Хо -f*

л'о —'-VI

- 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ш'1— 3 ш- 4- 1 =

0

 

 

 

имеет решение

 

(о, =0,618;

ш, =

1,618.

 

 

 

 

 

 

 

 

 

Формы колебаний —

 

 

 

 

 

 

 

 

 

 

 

?!

=

(1,

1,618) и 7 , = (1,-0,618).

 

 

 

Согласно (11.24)

решение

уравнения

частот

затухающих ко­

лебаний будет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_

т

4- і

ш.

 

л

л і

 

 

 

 

 

 

 

 

 

 

 

 

V

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Яй -

° '06 ХД 61-

± 0,618 і =

-

0,0185 ± 0,6181;

 

 

 

 

 

Із4 =

- 0,0485 ±

1,618 і.

 

 

 

Подстановка полученных четырех значений л в уравнение

(11.42)

показывает, что

они

действительно

удовлетворяют

уравнению

(с точностью до третьего десятичного знака).

форм

колеба­

Декременты

колебаний

 

одинаковы

для обеих

ний:

0,0185 X 6,28

 

 

 

 

 

 

 

 

 

°і

=

0,188 =

3,14 X 0,06;

 

 

 

 

0,618

 

 

 

о2

0,0485 X 6,26

=

0,188.

 

 

 

 

 

 

 

 

1,618

 

 

 

 

 

 

 

 

 

 

Матрица собственных

 

векторов —

 

 

 

 

 

 

 

 

П

 

 

 

 

 

 

і у

 

 

 

 

 

 

 

 

 

U.618 -0 ,6 1 8 /’

 

 

 

обратная матрица —

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ Г 1=

 

 

0,618

1

 

 

 

 

 

 

2-236 I 1,618

- 1

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычислим

преобразование

подобия матрицы К-

3,618

 

 

3 - 1 \

/

 

1

 

1

 

\

1 /1,382

 

- 1

2

\ 1,618

0,618/ = ~уг 5~ \2,236

- 2 , 236

68