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

Айзерман М.А. Классическая механика (1980)

.pdf
Скачиваний:
1498
Добавлен:
15.08.2013
Размер:
13 Мб
Скачать
☆

242

Г л . VI. ДВИЖЕНИЕ ВБЛИЗИ ПОЛОЖЕНИЯ РАВНОВЕСИЯ

т. е. опустим знак *, указывающий, что соответствующие коэффициенты получены при разложении в ряд функции Q*(q, q), так как это обстоятельство не будет играть далее существенной роли.

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

(55)

где Qj (t) —зависящие явно от времени части обэбдэнных сил. Предполагается, что в процесса движения они остаются малыми по модулю и не выводят систему из малой окрестности положения равновесия —к этому вопросу мы еще вернемся ниже.

В связи с тем, что система

уразнгний (55) является линейной,

а для линейных систем имеет место принцип суперпозиции, можно

рассмотреть движение системы

под действием какой-либо одной

силы из Qj(t) (/' = 1, ...,/г),

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

равны нулю. Определив порознь движения, возникающие под

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

затем

сложить.

 

Учитывая это обстоятельство, положим

 

Q.(Q = Qs(')=" - =Q»(O = O,

(56)

т.

е.

будем считать, что отлична от нуля только обобщенная сила

Q1

(t),

относящаяся к первой обобщенной координате, а все осталь-

ные обобщенные силы такого рода равны нулю.

 

Общее решение системы линейных дифференциальных уравне-

ний (55) складывается из двух слагаемых: первым является общее решение соответствующей однородной системы, получающейся из

(55) при Qj(t) = O ( / = 1 , .... п)

(обозначим его через q*), а вто-

рым — частное решение

соответствующей неоднородной

системы

(обозначим это частное решение через q**). Тогда

 

qj^qj

+ qj*

( / = 1 , . . . , л ) .

(57)

Что касается общего решения однородной системы q*, то оно находится по общим правилам интегрирования линейных однородных уравнений и в рассматриваемом случае движения вблизи положения асимптотически устойчивого равновесия заведомо стремится к положению равновесия при неограниченном возрастании времени t. В связи с этим движение q(t) стремится в пределе к движению q** (t), которое обусловлено наличием в правых частях уравнений зависящей явно от времени вынуждающей силы Qr (t).

Движения q*, которые бы возникали при отсутствии такой вынуждающей силы, называются свободными. Если этими движе-

§ 7. ВНЕШНЯЯ СИЛА, ЗАВИСЯЩАЯ ЯВНО ОТ ВРЕМЕНИ

243

ниями являются колебания, то их называют свободными колебаниями. Движение q**, которое возникает благодаря наличию вынуждающей силы, зависящей явно от времени, и к которому в пределе стремится суммарное движение, называют вынужденным движением (в случае колебаний говорят о вынужденных колебаниях).

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

1. Гармоническая вынуждающая сила. Частотная характеристика. Предположим, что обобщенная сила Q1 является гармонической функцией от t,

Q, (0 = A sin Q^.

(58)

Здесь А — амплитуда, а и —частота внешней вынуждающей силы Qx. Нам удобно далее вместо истинной обобщенной силы рассмат-

ривать комплексную функцию

Qj (/) = Aeiat = A (cos Qt + i sin Qf).

(59)

Интересующая нас действительная функция (58) получается как мнимая часть этого комплексного выражения:

0-

(60)

Рассмотрим уравнения (55), в которых Q1 (t) определяется выражением (59), a Qi = Qs= ... = Qn = 0. Нас интересует частное решение такой системы дифференциальных уравнений. Зная это решение и учитывая формулу (60), можно выделить в полученном решении мнимую часть и найти истинные вынужденные колебания.

Частное решение уравнений (55) с правой частью (59) будем искать в виде

qf = B^«

(/=1

п).

(61)

Подставив эти выражения

в наши уравнения и сократив на экспо-

ненциальный множитель, отличный от

нуля при любом t,

получим

J ] [a,n*Q)2 + M ^ )+ <7*]5ft = M

(/= 1, ..., п) (62)

А = 1

 

 

 

 

 

 

(здесь 6jy —обычный

символ

Кронекера: 6i/ = l при / = 1 и 61;

= О

при }Ф\.)

 

 

 

 

 

 

Если ввести

обозначения

 

 

 

 

 

dlk

(iQ) = ajk

(iQf

+ b,k (iQ) + c/k,

(63)

то уравнения (62) представятся в

виде

 

 

2

d}k

(iQ) Bk = б1;-А

(/ = 1

я).

(64)

244

ГЛ VI. ДВИЖЕНИЕ ВБЛИЗИ ПОЛОЖЕНИЯ РАВНОВЕСИЯ

 

Уравнения (62) представляют собой систему линейных алгеб-

раических уравнений относительно неизвестных Bk\ коэффициентами этих уравнений являются комплексные величины, определяемые по формулам (63).

Решая систему (64) линейных алгебраических уравнений с ком-

плексными коэффициентами по правилу

Крамера, получаем

 

В * = - ^ М ,

 

(65)

где в знаменателе стоит

определитель матрицы из коэффициентов

системы алгебраических

уравнений (64)

 

 

 

dn(iQ)

... dln(lQ)

 

 

 

(66)

 

dal(iQ) ...

dnn(iQ)

Легко видеть, что если

в этом определителе

заменить t'Q на ком-

плексное переменное Я, то он совпадает с характеристическим полиномом (24). Но у асимптотически устойчивой системы все нули характеристического полинома расположены слева от мнимой оси. Поэтому в таком случае определитель (66) отличен от

нуля

при любом Q.

 

 

 

 

В

числителе в формуле (65) стоит определитель Alft — алгебраи-

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

первой строке и k-м столбце

элемента определителя Д.

 

 

 

 

Как определитель А, стоящий в

знаменателе выражений (65)

и не зависящий от индекса k, так и определители Alft,

стоящие

в числителе этих выражений, представляют

собой полиномы от Я

с комплексными коэффициентами. Поэтому

отношение

определи-

телей

в выражении (65) является дробно-рациональной функцией.

Вид этой функции зависит от k. Заметим здесь же,

что степень

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

степени

знаменателя.

Более того, легко видеть, что в невырожденных случаях степень

числителя заведомо меньше степени знаменателя.

'

Обозначим эту дробно-рациональную функцию через

Wlh](iQ):

^• = Wlk(i®) = \Wlk\eiarsWik.

(67)

Тогда искомое частное решение неоднородной системы дифференциальных уравнений (55) при учете обобщенной силы (59) можно представить в виде

q*= А | Wlk (/Q) | е' IQi+arg *i* «аП (k = 1, ..., л). (68)

Вспомним теперь, что мы заменили интересующую нас обобщенную силу (58) комплексным выражением (59) и что (§8) является мнимой частью выражения (59). Выделив в полученном

S Г ВНЕШНЯЯ СИЛА, ЗАВИСЯЩАЯ ЯВНО ОТ ВРЕМЕНИ

245

частном решении (68) мнимую часть, находим искомое частное

решение системы

(55) при условии (58), а именно

 

 

 

qk

— A Wlk(iQ)

sin[Qt+

arg Wlh(iQ)]

(k = 1,

... ,

п).

(69)

Сравнивая формулы (69), определяющие вынужденные движе-

ния,

возникшие

благодаря

действию

вынуждающей

силы

(58),

с выражением для этой силы, устанавливаем, что

в этом случае

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

бания той

же

частоты, но с иными

амплитудами

и со

сдвигом

фаз. Амплитуды и фазы вынужденных

колебаний полностью опре-

деляются введенной выше комплексной функцией Wlk(iQ),

и для

данной системы зависят поэтому только от частоты внешней силы Q.

 

Введенная выше функция Wlk(iQ),

играющая столь

важную

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

ной

характеристикой системы, или, как говорят иногда, ее ампли-

тудно-фазовой характеристикой. Индексы при Wlk

показывают,

что

эта дробно-рациональная

функ-

imW

 

ция мнимого переменного зависит, с

 

одной стороны,

от

того,

к

какой

 

 

 

координате относится внешняя гар-

 

 

 

моническая сила, —это обстоятельст-

 

 

 

во

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

 

 

ЛеИ/

(в

рассматриваемом

нами

случае,

 

 

когда сила действует на первую

ко-

 

 

 

ординату,

этот

индекс

равен

1),

 

 

 

и,

с другой стороны, от того, вы-

 

 

 

нужденные

колебания какой

коор-

 

 

 

динаты исследуются,--это обстоя-

Рис.

VI.12

 

тельство подчеркивается вторым ин-

 

 

 

дексом k.

В зависимости от того, на какую координату действует

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

Рассмотрим теперь несколько подробнее только что введенную

важную характеристику системы

—ее частотную характеристику.

В комплексной плоскости можно

построить годограф функцииW

(рис. VI.12). Для этого надо, подставляя в выражение (67) зна-

чения Q, меняющиеся от 0 до + о о , подсчитывать порознь дей-

ствительные

и мнимые части этого выражения и по точкам строить

годограф. Начинаясь на

действительной оси (так как W при

Q = 0 равно

отношению

свободных членов полиномов Alft и Д),

246

ГЛ. VI ДВИЖЕНИЕ ВБЛИЗИ ПОЛОЖЕНИЯ РАВНОВЕСИЯ

 

этот годограф стягивается к нулю

(как уже было указано выше,

в выражении (67) степень знаменателя в невырожденных

случаях

выше степени числителя). Если

теперь

отдельно

рассмотреть

изменения модуля

и аргумента

вектора

W {19.) в

зависимости

от й, то получатся

характеристики, которые называются соответ-

 

 

 

 

ственно

амплитудной

 

 

 

 

и фазовой характеристи-

 

чно)!

 

 

ками системы (рис.VI. 13

 

 

 

 

и VI.14).

 

 

 

 

 

 

Зная частотную ха-

 

 

 

 

рактеристику

системы,

 

 

 

 

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

Рис.

VI.13.

Рис.

VI.14.

плитуду и сдвиг

фазы вынужденного

колебания

координаты

при любой частоте внешней силы, скажем Qt. Для того чтобы сделать это, надо по рис. VI.12 (либо по рис. VI.13 и VI.14

J2?

д

 

Рис. VI.15.

 

соответственно) найти при Q = Qx

модуль и аргумент

функции

W (Ш), затем умножить амплитуду вынуждающей силы на \W (iQJl

и сдвинуть фазу гармонической функции, выражающей

силу, на

arg W(iQ,).

 

 

Если амплитудная характеристика системы (рис.VI.13) при некотором значении Q = Q* имеет отчетливо выраженный пик,

 

 

 

 

7 ВНЕШНЯЯ СИЛА, ЗАВИСЯЩАЯ ЯВНО ОТ ВРЕМЕНИ

 

 

247

т. е. если

| W (iQ) \ при йя« Q* значительно больше, чем | W (iQ)

при остальных

 

значениях

Q, то при одной и той же

амплитуде

внешней

силы

Qx

 

амплитуда

отклика

резко возрастает,

когда

частота

внешней

силы приближается к значению

Q = Q*. Это

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

а частота Q* —резонансной часто-

той системы.

Разумеется,

амплитудная

характеристика системы

может

иметь

несколько

пиков

при

различных

значениях

Q*

(рис. VI. 15). В

этом случае говорят,

что система

резонирует

на

нескольких

частотах.

В

рассматриваемом

здесь

общем

случае

стационарной

системы

понятие

 

«резонанс»

является

не точным,

а скорее интуитивным, так как

 

оно связано с недостаточно чет-

ким

понятием

 

«отчетливо

выраженный

пик». Это понятие при-

обретает

точный

 

смысл

лишь в случае,

когда система

консерва-

тивна

(см. ниже).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В силу того что амплитудная

характеристика,

как

уже

было

указано

выше,

стягивается к нулю с ростом Q,

можно указать

такое значение «частоты среза» Q = Qcp, что при всех

Q >• Qcp и

при заданном

А

 

решение

q%* п о

 

 

 

 

 

 

 

 

модулю будет

меньше некоторого

 

 

 

 

 

 

 

 

наперед

заданного

числа

е > 0

 

 

 

 

 

 

 

 

(см. рис. VI.13 и VI.16, а). Задав-

 

 

 

 

 

 

 

 

шись

 

этим

числом

е и определив

 

 

 

 

 

 

 

 

таким образом значение Qcp, гово-

 

 

 

 

 

 

 

 

рят,

что рассматриваемая

система

 

 

 

 

 

 

 

 

в малых

колебаниях «пропускает»

 

 

 

 

 

 

 

 

лишь частоты

от

нуля до Q = Qcp,

 

 

 

 

 

 

 

 

а частоты,

большие

Qcp, «не про-

 

 

 

 

 

 

 

 

пускает».

 

В

этом

смысле меха-

 

 

 

 

 

 

 

 

нические системы

 

являются филь-

 

 

Рис. VI.16.

 

 

 

тром высоких частот.

Наоборот,

 

 

 

 

 

 

 

 

 

 

 

 

 

можно

таким

 

образом

подоб-

 

 

 

 

 

 

 

 

рать

 

коэффициенты

a,

b

и с

рассматриваемой

системы

урав-

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

представленный на

рис. VI. 16, б,

т. е.

чтобы модуль

ампли-

тудной характеристики был больше

наперед выбранного

числа е

лишь тогда, когда Q

лежит между

двумя

числами пср1

и Йс р 2 ,

и был бы меньше числа е вне этого диапазона. В таких

случаях

говорят, что система

является полосовым фильтром, пропускаю-

щим полосу частот, заключенную между й с р 1 и Qcp2-

 

Рассмотрим теперь консервативную систему. Пусть а>и о)2, ...

..., «„ —ее собственные частоты (см. § 6 этой главы). Собственными колебаниями системы служат гармонические колебания с этими частотами, а это означает, что все корни характеристического уравнения консервативной системы —чисто мнимые и что они равны

±

±i

±

Рис. VI.17.

248

ГЛ VI ДВИЖЕНИЕ ВБЛИЗИ ПОЛОЖЕНИЯ РАВНОВЕСИЯ

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

(Я) =Ао П (К- i

©I)

(об этом уже шла речь выше), т. е. содержит лишь четные степени к. Ясно, что и алгебраическое дополнение любого элемента определителя А(X) обладает этим же свойством. Поэтому в случае консервативной системы частотные характеристики, определяемые формулами (67), являются не комплексными, а действительными функциями, и | W (iQ) | принимает бесконечные значения при

Q Q

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

Здесь понятие резонансной частоты имеет точный смысл: уконсервативной системы резонансной частотой называется значение Q, при котором амплитудная |MJ2)| характеристика системы имеет разрыв второго рода.

Из изложенного следует, что в случае, когда частота внешней силы приближается к любой из собственных частот консервативной системы, амплитуды вынужденных колебаний всех ее обобщенных координат неограниченно возрастают1).

Рассмотрим теперьфазовую характеристику консервативной системы. В связи с тем, что для консервативной систе-

мы Wlk(iQ) —действительная функция, аргумент ее может быть равен лишь нулю или кратен —л в зависимости от знака этой функции. Поэтому фазовой характеристикой консервативной системы служит кусочно-постоянная функция, равная 0 или крат-

! ) Мы не выясняем

здесь, по

какому закону

происходит

это нарастание

амплитуд

во

времени.

Читатель,

интересующийся

подробностям, найдет рх

в любом

курсе

теории колебаний (см.,например, Б у л г а к о в

Б. В. Колеба-

н и я . — М - Гостехиздат,

1954).

 

 

 

S 7 ВНЕШНЯЯ СИЛА, ЗАВИСЯЩЛЯ ЯВНО ОТ ВРЕМЕНИ

249

ная — я . Она имеет разрыв первого рода при всех значениях Q, равных собственным частотам системы

Q 1 = (D1, Й 2 = ( О 2 , . . . , Q n = (Un.

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

aq+cq— A sin Q/.

Свободные колебания системы

q* = С sin (<o/-+-ф)

происходят с собственной частотой

<л = У с/а,

а вынужденные колебания описываются функцией

В этом простейшем случае

1

1

1

W (f Q ) ••

с — afi2

a(co2 — Q 2 ) '

" а ( Ш ) 2 + с

амплитудная

характеристика

имеет

вид

 

 

 

 

 

 

\W(iQ)

1

 

1

 

 

 

 

 

 

 

со

2 — Q2

 

 

 

 

 

 

 

 

 

 

 

 

 

а фазовая характеристика такова:

 

 

 

 

 

 

 

 

 

 

О при Q<co,

 

 

 

 

 

 

 

— я

при Q > со.

 

 

 

Эти три характеристики показаны на

рис. VI. 18.

 

 

 

Читателю

рекомендуется

самому

найти

явное выражение для

вынужденных

колебаний

консервативной

системы

с

п

степенями

свободы,

придерживаясь

следующего

плана.

 

 

 

1° Найти

матрицу преобразования

обобщенных

сил

Qly Q2, ...

. . . , Qn

(для

системы координат qx,

q2, ... , qn)

в

обобщенные

силы Qx,

0г>

•••. ®л (для главных

координат Вх,

6Я, . . . , бл). Для

этого надо приравнять выражения элементарной работы через обобщенные силы Q и в , представить в этом равенстве q как функции 0 при помощи преобразования (45) и изменить порядок суммирования. Читатель установит тогда, что искомая матрица преобразования Q в G получается транспонированием амплитудной матрицы v(ft.

250

ГЛ VI ДВИЖЕНИЕ ВБЛИЗИ ПОЛОЖЕНИЯ РАВНОВЕСИЯ

2° После того как обобщенные силы для главных координат найдены, вынужденные колебания каждой из этих координат 6у определяются по формуле (70).

k\mW

W(iQ)

\W(fO)\ щЩп)

 

 

со

 

 

 

 

 

•V* I

1

 

 

 

 

 

Рис. VI.18.

 

 

 

 

3°

Вынужденные колебания координат q} находятся

по фор-

муле

(45).

 

 

 

 

 

2. Периодическая, но не гармоническая вынуждающая сила.

Рассмотрим теперь

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

нату

действует не

гармоническая, а периодическая обобщенная

сила

с периодом Т,

заданная функцией Q*(t),

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

 

 

условию

 

 

 

Qj(t)

 

Ц\ (t) = Щ (t -)- / ),

 

(/1)

 

 

например

функцией,

график

 

 

которой

изображен

на

 

 

рис.VI. 19.

 

 

 

 

 

Мы

будем предполагать^-

 

 

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

 

 

Q* (t)

представима

 

рядом

 

 

Фурье

 

 

 

 

 

 

Qt(t)=

 

 

 

 

Рис. VI. 19.

 

 

 

(72)

где амплитуды гармоник Ak и соответствующие сдвиги фаз срА определяются по обычным правилам разложения периодической функции в ряд Фурье, a Q = 2я/Т —круговая частота. Теперь внешняя сила, действующая на первую координату, представлена как сумма гармонических колебаний. В силу линейности системы

§ 7 ВНЕШНЯЯ СИЛА, ЗАВИСЯЩАЯ ЯВНО ОТ ВРЕМЕНИ

25!

и действующего поэтому принципа суперпозиции каждая из этих гармонических сил вызывает независимое вынужденное колебание, а общее вынужденное колебание, возникающее под действием такой периодической силы, получается суммированием этих независимых колебаний. Для определения каждого из вынужденных колебаний, которое возникает в том случае, когда внешняя сила представляется не всем рядом (72), а лишь какой-либо одной гармоникой, например k-й, можно воспользоваться полученной выше формулой (69) —надо лишь заменить всюду Q на Ш. Поэтому вынужденное колебание /-й координаты <7/, которое возникает под действием периодической силы, действующей на первую координату <7Х и выражающуюся рядом (72), может быть представлено в виде

 

 

 

qf* =

 

 

 

 

 

 

sin [kQt + arg Wy (ikQ)].

(73)

Из этой формулы видно, что вынужденные колебания, возни-

кающие в системе

под

действием внешней силы (72),

полностью

определяются

частотной

характеристикой

системы так

же,

как и

в случае, когда

рассматривалась гар-

 

imW

 

 

моническая

сила.

Но теперь

на ча-

 

 

 

стотной характеристике надо рассмат-

 

 

 

 

ривать не только точку, соответствую-

 

 

 

 

щую

частоте

Q,

но

и все точки, со-

 

 

 

 

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

 

частотам

Ш

(k = 0,

 

 

 

ReW

1, ...). Отмечая эти точки

на частот-

 

 

 

 

 

 

 

ной

характеристике

(рис. VI.20) и

 

 

Л

 

вспоминая о наличии полосы про-

 

 

 

пускания,

благодаря

чему практиче-

 

 

 

 

ски

оказывается

 

необходимым

рас-

 

Рис. VI.20.

 

смотреть

лишь

конечное (и

обычно

 

 

 

 

 

 

небольшое)

число

таких

 

точек,

 

 

 

 

мы можем для

каждой

из этих точек

определить модуль

частот-

ной

характеристики

и ее аргумент и, подставив их в формулу (73),

найти

вынужденное

колебание. Этот ряд можно изобразить графи-

чески, откладывая

в точках

0,

Q, 2Q,

... оси Q значения ампли-

туд

гармоник

Ak

 

и соответствующих сдвигов фаз <pk

(рис. VI.21).

Такой

график

называется

линейчатым

спектром воздействия.

Аналогично возникающее в результате вынужденное движение также представимо рядом Фурье и изображается своим линейчатым спектром. Частотная характеристика W (iQ) в этом случае играет роль оператора, преобразующего линейчатый спектр возмущающей силы в линейчатый спектр вынужденного движения.

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

Соседние файлы в предмете Теоретическая механика