Теоретические основы динамики машин

         

Асимптотический метод расчета пластин


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

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

Рассмотрим применение этого метода на примере заделанной по контуру прямоугольной пластины размерами

, у которой, как и ранее, размер а соответствует оси
. Ограничимся расчетом симметричных относительно осей
 форм колебаний. В средней части пластины (начало координат располагается в центре тяжести пластины) принимаем

 

Вблизи границ

:

где

 - быстро изменяющаяся функция, позволяющая удовлетворить условиям закрепления.

Аналогично вблизи границ

:

Таким образом, общее выражение для

 имеет вид

           (322)

В средней части пластинки функции

 и
 пренебрежимо малы и поэтому первый член выражения (322) должен удовлетворять уравнению (316). Отсюда находим

Вблизи границ

 существенными являются первый и второй члены выражения (322). Учитывая, что первый член удовлетворяет уравнению (316), потребуем, чтобы и второй удовлетворял ему:

Выполняя дифференцирование, приходим к уравнению:

которое распадается на два уравнения:



Так как

, то затухающие решения имеет только первое из этих уравнений. Решение, затухающее с удалением от стороны
, имеет вид

где                          

В силу симметрии вблизи стороны

:

Аналогично вблизи стороны

:

где                           

и вблизи стороны

:

Рассмотрим граничные условия при

:

При вычислении

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




Для одновременного выполнения этих уравнений необходимо, чтобы определитель, составленный из коэффициентов при
 и
, равнялся нулю, что приводит к уравнению:

                             (323)

Аналогично условия при
 приводят к уравнению:

                                  (324)

Так как
 связаны с
, то трансцендентные уравнения (323) и (324) позволяют определить значения
 и
, а затем вычислить
 и частоты:



Рассмотрим, например, колебания квадратной пластинки с одинаковым числом узловых линий в направлениях
. В этом случае
 
 и уравнения (323) и (324) приводят к зависимости:

,

откуда



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



Достаточно хороший результат получается уже для низшей частоты:



Точное значение:



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


Бегущие волны в круглых пластинках


Рассмотренные выше собственные колебания круглых пластинок описываются уравнением

                                    (330)

Они соответствуют стоячим волнам на поверхности пластинки, при которых узловые диаметры неподвижны.

Наряду с (330) решением уравнения движения является также выражение

                                (331)

Но поскольку уравнение движения линейно, сумма и разность  (330) и (331) также являются его решениями:

Эти выражения представляют собой уравнения бегущих волн. Первое выражение соответствует вращению всей картины деформаций вокруг оси симметрии пластинки в направлении возрастания угла

 с угловой скоростью
. Второе выражение соответствует движению волны с той же скоростью в противоположном направлении.

Если имеются внешние нагрузки, вращающиеся по периферии пластинки со скоростью, близкой к скорости

 распространения собственных колебаний, то такие нагрузки вызовут большие резонансные колебания пластинки.

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

Критические скорости вращения диска

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

,                                               (332)

где

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



Частотное уравнение и собственные формы


Развёрнутая запись граничных условий приводит к однородным уравнениям относительно постоянных C1, C2, C3, C4.

Чтобы эти постоянные не равнялись нулю, должен равняться нулю определитель, составленный из коэффициентов системы; это приводит к частотному уравнению. При этих операциях выясняются соотношения между C1, C2, C3, C4, т.е. определяются собственные формы колебаний (с точностью до постоянного множителя).

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

Для балки с шарнирно-опёртыми концами согласно (203) имеем следующие граничные условия: X=0; X''=0 при x=0 и x=

. При помощи  (197)-(200) получим из первых двух условий: C1=C3=0. Два оставшихся условия можно записать в виде

Чтобы C2 и C4 не были равны нулю, необходимо равенство нулю определителя:

.

Таким образом, частотное уравнение имеет вид

.

Подставляя выражения T и U, получим

.

Так как

, то окончательно частотное уравнение записывается так:

.                                               (207)

Корни этого уравнения:

, (n=1,2,3,...).

Учитывая (196), получим

.                                                (208)

Перейдём к определению собственных форм. Из записанных выше однородных уравнений вытекает следующее соотношение  между постоянными C2 и C4:

.

Следовательно, (197) приобретает вид

или

Согласно (207), имеем

,                                                (209)

где

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



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


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

вызовет сложные колебания

.                                  (101)

Если частоты

и
близки одна к другой, то суммарное движение принимает своеобразный характер и называется биением (рис.43,a). Особенностью биений является периодическое изменение амплитуды колебаний. Для анализа этого явления преобразуем сначала решение (101).

   =

,

или

,                           (102)

где

Обозначая 

, (102) перепишем в виде

.

Вследствие того, что частоты

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

,

где      

- среднее значение частоты;

- медленно меняющаяся амплитуда колебаний;

 - медленно меняющаяся фаза.

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

TA=

.

Так как разность

 мала, то период ТА значительно больше периода колебаний

.

Биения могут возникнуть и при действии одной возмущающей силы

 вблизи резонанса, когда частота p близка к собственной частоте
. Из решения (96) видно, что в данном случае колебания состоят из двух гармоник с близкими частотами p и
. Преобразования, подобные выполненным  выше, приводят к выводу, что и здесь суммарные  колебания носят синусоидальный характер с переменной амплитудой. Однако этот процесс не является стационарным; свободные колебания вследствие затухания постепенно исчезают, и остаются только вынужденные колебания, так что биения прекращаются (рис.43,б).

Рис. 43



Действие гармонической силы


Случай, когда возмущающая сила изменяется по гармоническому закону

F=F0  sin pt,                                            (94)

где  F0 - амплитуда силы; p - её частота, является наиболее распространённым в практике.

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

.                                   (95)

Вычисляя интеграл, при

  находим                   

.                              (96)

Заменим

 и обозначим
 (прогиб, вызванный статически приложенной постоянной силой F0 ), тогда

.                                (97)

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

. Обычно первые колебания называют вынужденными, а вторые - свободными. Такая терминология является условной. Дело в том, что и вторые колебания вызваны действующей возмущающей силой, и их амплитуда зависит от этой силы; в этом смысле вторые колебания также являются вынужденными. Указанные наименования получили широкое распространение потому, что первое слагаемое имеет частоту возмущающей силы, а второе меняется с собственной частотой системы.

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

 .                                            (98)

Амплитуда вынужденных колебаний

                                                 (99)

отличается от прогиба

, подсчитанного в предположении статического действия силы F0. Отношение
 можно назвать динамическим коэффициентом

 .                                                (100)

Динамический коэффициент

 зависит только от отношения частот
. На рис.41,a приведена кривая зависимости
от отношения
.

Рис.  41

При малой частоте возмущающей силы динамический коэффициент близок к единице. С ростом частоты p динамический коэффициент быстро увеличивается и при

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

   Если частота p больше собственной частоты 
, то амплитуды становятся конечными; при
 динамический коэффициент становится меньше единицы, т.е. динамический эффект слабее соответствующего статического эффекта. При очень больших значениях отношения
 динамический коэффициент становится весьма малым. Это означает, что сила высокой частоты не вызывает ощутимых колебаний в упругой низкочастотной системе, которая как бы "не успевает" отзываться на быстрые изменения возмущающей силы.

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

,

где
- масса ротора; e- её эксцентриситет;
- угловая скорость.

В данном случае амплитуда возмущающей силы
пропорциональна квадрату
, и вместо решения (96) при
 следует принимать

.

Амплитуда стационарных колебаний при этом определяется выражением

,


в котором параметр системы
 не зависит от частоты p.

На рис.41,б представлено изменение амплитуды колебаний в зависимости от отношения
.

Как видно, при
 имеет место резонанс, а при p»
амплитуда стремится к значению
.

Остановимся подробнее на случае совпадения частот
(резонанс).



Рис. 42

При этом интеграл (95) принимает вид

.

После вычисления получим

.

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


Действие непериодической возмущающей силы


1. Действие линейно возрастающей силы (рис.38,а)

Используем выражение (91), полагая

График этого движения показан на рис.38,б.

Рис. 38

Перемещения нарастают по сложному закону, представляющему собой сумму синусоиды и линейной функции. Дополнительное синусоидальное колебание тем существеннее, чем быстрее нарастает сила F, т.е. чем больше а.

Колебания подрессоренного груза при движении по неровной дороге (рис.39)

Пусть профиль дороги задан уравнением

=h(1-
),

где h - предел, к которому стремится высота профиля;

- параметр, характеризующий кривизну профиля.

Рис. 39

Обозначим  через V скорость движения груза массой m и примем начало отсчёта времени в мгновение, когда опорная точка проходит начало неровности. Тогда х=Vt и движение опорной точки по вертикали определяется законом

).

Дифференцируя, находим

.

На основании  (92) получим закон движения груза по вертикали:

Z=f(t)-

.

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

.

Интегрируя, находим

,

где

 определяется соотношением

.

Отсюда следует, что при весьма малой скорости параметр

  стремится к
, а разность z-f - к нулю. Наоборот, при весьма большой скорости (а также при весьма большом значении параметра
) параметр
стремится к нулю, и колебания приближённо описываются законом

2. Действие медленно изменяющихся сил

Рассмотрим полученное выше решение задачи о вынужденных колебаниях (91). Первое слагаемое представляет собой статическое отклонение, вызванное силой F(t). Второе слагаемое - это поправка к статическому отклонению, причём эта поправка зависит от скорости изменения силы

.

При малой скорости нарастания внешней нагрузки динамическая поправка к статическому решению относительно мала и нагружение практически можно рассматривать как статическое.

А.Н.Крылов дал оценку динамической поправки для общего случая возмущающей силы. Если кривая F(t) имеет один максимум (рис.40,а), то, обозначая максимальное значение

 через
(рис.40,б), имеем


.

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

.



Рис. 40

Если сила возрастает равномерно в течение времени t0, то
 и динамическая "добавка" составляет

.

Её относительная величина

.

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

3. Действие быстро исчезающих сил

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

Рассмотрим действие силы F, которая внезапно прикладывается в момент времени t=0, действует в течение некоторого промежутка времени а, а затем также внезапно исчезает. Можно показать, что если a<
, то максимальное отклонение системы достигается после исчезновения силы. В таком случае для t>à, согласно решению (89), имеем

                               
.                 (93)

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

.

Максимальное отклонение  в соответствии с (93):

.

Следовательно, динамический коэффициент:



Значения динамического коэффициента
при действии силы малой продолжительности приведены в табл. 5.

                                                                                                 Таблица 5



0

0,01

0,02

0,03

0,05

0,10

0,15

0,25

0,5



0

0,062

0,126

0,188

0,313

0,618

0,908

1,413

2,000

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



Покажем, что действие кратковременной силы приближённо можно оценить её импульсом. Для t>à решение имеет вид

,

или

=

=
.

Но так как отношение
  меньше отношения
, то
 есть малое число.

Тогда приближённо можно записать

.

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


Действие периодических импульсов


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

Рассмотрим какой-либо период Т, принимая начало отсчёта времени в конце действия предшествующего импульса (например, в момент времени

). Обозначим перемещение и скорость в момент времени
 через
 и
.

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

, т.е. описываются уравнением

                                      (104)

и, следовательно,

В конце этого периода, непосредственно перед следующим импульсом (мгновение

)

В результате действия очередного импульса скорость мгновенно изменится на величину

 (где S - значение импульса). Поэтому непосредственно после следующего импульса  (мгновение
)

Вследствие периодичности процесса эти величины должны быть равны

и

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

и
.

Закон движения (104) принимает вид

                                (105)

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

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

Дробь

 есть максимальное отклонение, вызванное одним импульсом, поэтому выражение

  

можно назвать коэффициентом повторности. Этот коэффициент характеризует возрастание влияния повторных импульсов. Из формулы для величины

 видно, что при совпадении частот или их кратности
 
 
 возникает резонанс. Максимально возможное значение коэффициента повторности равно 0,5.



Действие произвольной периодической


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

Так, на рис.44,а показан закон изменения крутящего момента, создаваемого четырёхтактным двигателем внутреннего сгорания. Другой пример – периодические "безмассовые" удары – показан на рис.44,б.

Силы (моменты) рассматриваемого вида имеют чётко выраженный период колебаний Т, но не описываются единым аналитическим выражением. В подобных случаях чаще всего пользуются разложением периодической нагрузки в ряд Фурье. При этом сила представляется в виде суммы гармонических составляющих, а затем определяется эффект, вызываемый каждой из составляющих; после этого полученные частные эффекты суммируются.

Периодическую силу F(t) можно представить в виде ряда Фурье:

где

 - основная частота возмущения.

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

и
вычисляются по формулам

               

             

.   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .  .   .   .   .   .

Рис. 44

Опираясь на решение (98), полученное для одной гармоники, находим

  

             (103)

Это решение состоит из постоянного слагаемого

, соответствующего среднему значению возмущающей силы, и ряда, соответствующего гармоническим колебаниям с частотами p, 2p, .... Если собственная частота совпадает с частотой какой-либо одной гармоники np (n=1,2,...), то соответствующее слагаемое в  (103) стремится к бесконечности. Следовательно, в общем случае периодической возмущающей силы резонанс наступает не только тогда, когда собственная частота
равна основной частоте возмущающей силы p, но и когда
 кратно p (в некоторых частных случаях в  (103) пропадают некоторые слагаемые, и резонанс наступает не при любой кратности).

Рассмотренный способ чётко выявляет условия наступления резонанса. Недостатком этого способа является сложность вычислений, необходимых для учёта большого числа слагаемых в (103). Так, возмущающую силу, показанную на рис.44,a, для достаточной точности необходимо заменить примерно десятью гармониками.



Флаттер крыла в воздушном потоке


Как уже говорилось выше, флаттером называются автоколебания тел в потоке газа или жидкости.

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

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

Существенное значение имеет положение центра тяжести крыла, т.е. той точки его хорды, приложение вертикальной силы в которой вызывает только вертикальное перемещение крыла, но не его поворот. К этой точке (точка 0 на рис.65) будем приводить действующие на крыло силы.

Рис. 65

Если обозначить вертикальное перемещение центра жесткости крыла через

, а изменение угла атаки крыла в процессе движения через
, то упругие сила и момент, приложенные в точке 0, будут равны соответственно
 и
, где
 и
 - коэффициенты жесткости.

Сила инерции и момент сил инерции относительно точки 0 составляют соответственно

                 (165)

где

 - расстояние от центра жесткости крыла до его центра массы;
 - масса крыла;
 - радиус инерции массы крыла относительно центральной оси.

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

                               (166)

где

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


Формулы (166) представляют собой грубое приближение, так как в них полностью игнорируется влияние движения крыла на обтекание. Более точное решение задачи показывает, что если крыло совершает, например, гармонические колебания с частотой
, то следует учитывать еще инерцию присоединенной массы воздуха и то обстоятельство, что изменение подъемной силы оказывается смещенным по фазе относительно изменения угла атаки.

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

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

                      (167)

где       
.

Решение системы (167) отыщем в виде, соответствующем гармоническим колебаниям:

                                          (168)

Подставляя (168) в (167), получим систему однородных алгебраических уравнений относительно
 и
:

              (169)

Приравнивая нулю определитель системы (169), получим частотное уравнение. Для того чтобы привести это уравнение к более простому виду, введем следующие обозначения:



- собственные частоты поступательных (изгибных) и крутильных колебаний крыла;



- относительная плотность крыла.

Тогда частотное уравнение можно представить в так:

       (170)            

При нулевой скорости потока V=0 это уравнение даёт два положительных значения W2, соответствующих  двум  собственным  частотам  системы.

С увеличением скорости потока возможно появление двух типов неустойчивости. Так, один из корней уравнения (170) может обратиться в нуль, что соответствует обращению в нуль свободного члена уравнения (170):

                                   (171)

Обращение в нуль частоты собственных колебаний системы свидетельствует о её статической неустойчивости. Действительно, возвращаясь в формуле (171) к первоначальным обозначениям, приведём её к виду





Если это соотношение выполняется, то при повороте крыла на угол а момент дополнительной подъёмной силы



уравновешивается  упругим моментом 
.

Явление статической потери устойчивости крыла при достижении скоростью потока значения Vg называется дивергенцией.

Для крыльев самолётов, как правило, скорость дивергенции существенно превышает скорость полета и дивергенция не представляет реальной опасности.

Другой вид потери устойчивости – изгибно-крутильный флаттер - связан с тем, что частоты, определяемые из  (170), становятся комплексными числами. Если имеются сопряжённые комплексные частоты
, то соответствующие решения уравнений движения имеют множители

.

Экспоненциальные множители с действительными   положительными показателями неограниченно возрастают.

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

Итак, условием наступления флаттера является появление комплексных корней уравнения (170), что происходит при обращении в нуль (назовём это «граничным условием») его дискриминанта:

                          (172)

Из уравнения (172) легко вычислить скорость флаттера.

Проследим на числовом примере характер изменения частоты свободных колебаний крыла по уравнению (170) при увеличении скорости потока.

Допустим, что
;
;
;
;
.

Этим данным соответствуют скорость дивергенции, вычисленная по уравнению (171),
, и скорость флаттера, вычисленная по уравнению (172),
.

График  изменения частот колебаний системы в зависимости от скорости потока, построенный в соответствии с уравнением (170),показан на рис.66.



Рис. 66

При V=0 система имеет две частоты собственных колебаний, мало отличающиеся от частот чисто крутильных и чисто изгибных колебаний. С увеличением скорости потока частоты сближаются, и при скорости флаттера оказываются равными друг другу.

Наличие кратных собственных частот для консервативной системы не связано с какими-либо особенностями её поведения.Для неконсервативной системы, которую представляет собой крыло, находящееся в потоке воздуха, слияние двух частот ведёт к потере устойчивости движения. В процессе колебаний система начинает интенсивно потреблять энергию потока и амплитуды колебаний неограниченно возрастают. Механизм этого явления легко понять, если представить себе, что происходящие с одинаковой частотой крутильные и изгибные колебания крыла сдвинуты по фазе на
, так что, когда крыло движется вверх, его угол атаки (а значит, и подъёмная сила) больше, чем когда оно движется вниз. При этом за полный цикл подъёмная сила будет совершать положительную работу и энергия колебания будет непрерывно возрастать.


Формула Донкерлея


Так как метод Рэлея приводит к завышенному значению частоты колебаний, то весьма полезным является применение метода (формулы), дающего заниженное значение низшей частоты колебаний. Простейшей из такого рода формул является формула Донкерлея.

Рассмотрим какую-либо многомассовую систему, например балку (рис.74,а). Пусть на этом рисунке изображена точная форма собственных колебаний системы.

Рис. 74

Тогда точное значение собственной частоты системы выражается формулой

                                                      

.                                          (275)

Теперь рассмотрим ту же балку, но только с одной массой mi

(рис.74, б). В этом случае частота колебаний будет определяться по формуле

,

где dii - податливость балки при приложении силы в точке расположения массы mi .

С другой стороны, приближённое значение wi той же частоты можно определить по формуле Рэлея, считая, что форма колебаний совпадает с изображённой на рис.74,а

.                                                (276)

Здесь По и xi имеют те же значения, что и в (275).

Так как форма, изображённая на рис.74,a, не является точной формой колебаний одномассовой системы, то выполняется неравенство:

 .                                                  (277)

Сравнивая (275) и (276), находим

.                                                       (278)

Eсли в правой части полученного равенства заменить

 меньшими значениями
, то равенство превратится в неравенство:

.

Таким образом, приближенная формула Донкерлея всегда дает преуменьшенное значение частоты:

.                                              (279)

Рассчитав частоту одной и той же системы по методу Рэлея и формуле Донкерлея, получим вилку, в которой заключено истинное значение частоты колебаний.




Фрикционные автоколебания


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

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

Рис. 62

На ленте лежит груз массой

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

                                                 (161)

Для дальнейших рассуждений необходимо учесть, что сила трения R зависит от относительной скорости движения V (при заданном нормальном давлении груза, которое будем считать постоянным). Пусть характеристика трения имеет вид, показанный на рис.62,б, где значения V0 и R0 соответствуют состоянию равновесия груза.

Положим, что вследствие какого-либо возмущения в мгновение t=0 груз выведен из состояния покоя. Выясним характер движения, которое возникает после такого возмущения, считая, что скорость движения ленты в процессе колебаний груза не изменяется. Пусть в текущее мгновение

 смещение груза равно
, его скорость равна
, а скорость скольжения   

.

В это мгновение на груз действуют три силы: сила трения R, реакция пружины

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

Сила

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

,                                                (162)

где

 - тангенс угла наклона характеристики трения в точке с координатами
.

Уравнение движения груза:

.

Подставляя сюда  (162), получим

или

                             (163)

Сместим начало отсчета перемещений в точку

, т.е. введем новую переменную
.

Тогда (163)  имеет  вид

Согласно (161), сумма двух последних слагаемых равна нулю, поэтому получим уравнение


                              (164)

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

Однако при небольших значениях
 (нисходящий участок характеристики трения) величина
 становятся отрицательной. Если при этом
, то сумма
 остается положительной, и колебания будут затухающими. Если сумма
 обращается в нуль, т.е.
, то в уравнении (164) исчезает член, соответствующий затуханию, и становятся возможными колебания с постоянной амплитудой. Если же сумма
 отрицательна, то система обладает как бы «отрицательным затуханием», и колебания с течением времени будут возрастать за счет энергии, передаваемой движущейся лентой.

Природу «отрицательного затухания» при падающей характеристике трения можно выяснить исходя из следующих рассуждений. Когда груз движется вправо, т.е. по ходу движения ленты, относительная скорость скольжения уменьшается; вместе с этим сила трения увеличивается и ее приращение направлено вправо, т.е. также в сторону движения. В другом интервале движения, когда груз движется влево, приращение силы трения направлено также влево, т.е. опять в сторону движения. Такой характер изменения силы трения и является причиной возрастания колебаний.

Итак, для возрастания колебаний необходимо выполнение условия:
, что возможно только при достаточной крутизне падения характеристики трения. Обычно указанное условие выполняется лишь при малой скорости
.

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

К этому же классу явлений относятся автоколебания, возникающие при резании металлов на станках. Рассмотрим природу этих автоколебаний.


Со стороны заготовки на резец действует реакция
, которая может быть разложена на составляющие
 и
 (рис.63,в,г). Резец упруго закреплен и его конец может совершать колебания как в горизонтальном, так и в вертикальном направлениях. Для выявления возможности автоколебаний достаточно рассмотреть колебания только в горизонтальном направлении
 и учесть важный экспериментальный результат: горизонтальная составляющая
 зависит от скорости горизонтальных колебаний резца
.

На резец действуют три силы: сила упругости - су системы резец-суппорт (где с - коэффициент жесткости); реакция заготовки
 и сумма различных неупругих сопротивлений, которая может быть объединена в одно слагаемое вида
.



Рис. 63

Таким образом, уравнение движения системы резец-суппорт имеет вид



где
 - приведенная масса системы резец-суппорт.

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

Как установлено, неустойчивость состояния равновесия может быть обнаружена в предположении малости колебаний, т.е. при помощи линейного приближения. Однако, если отказаться от этого предположения и  проследить за дальнейшим течением процесса, то обнаруживается, что рост амплитуд постепенно замедляется и в конечном итоге амплитуда полностью стабилизируется. Этот процесс установившихся (стационарных) автоколебаний (рис.64) называется предельным циклом.

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



Рис. 64

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

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


Граничные условия


Для каждого конца стержня можно указать два граничных условия.

Свободный конец стержня (рис. 70,а). Нулю равны поперечная сила Q=EJX'''T и изгибающий момент M=EJX''T. Поэтому граничные условия имеют вид

X''=0; X'''=0 .                                      (202)

Рис. 70

Шарнирно-опёртый конец стержня (рис.70,б). Нулю равны прогиб y=XT и изгибающий момент M=EJX''T. Следовательно, граничные условия  таковы:

X=0 ; X''=0 .                                            (203)

Защемленный конец (рис.70,в). Нулю равны прогиб y=XT и угол поворота

. Граничные условия:

X=0; X'=0 .                                                       (204)

На конце стержня имеется точечный груз массы

 (рис.70,г). Его сила инерции
 может быть при помощи уравнения (194) записана так:
; она должна быть равна поперечной силе  Q=EJX'''T , поэтому граничные условия принимают вид

; X''=0 .                                        (205)

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

Упруго-опертый конец стержня (рис.70,д). Здесь изгибающий момент равен нулю, а поперечная сила Q=EJX'''T равна реакции опоры

 (Co-коэффициент жёсткости опоры).

Граничные условия:

X''=0 ;

                                    (206)

(знак минус принимается в случае, когда упругая опора является левой, и знак плюс, когда она является правой).



Изгибные колебания балок


Рассмотрим случай, когда возмущающая нагрузка задана в виде сосредоточенной силы

                                           (305)

или комбинации нескольких нагрузок того же вида с одинаковой частотой. Решение для прогибов будем искать в виде

                                     (306)

сводя задачу к определению формы колебаний (кривой амплитуд прогибов)

.

В случае

, подставляя в (192)  (306), получим

                                            (307)

Решение дифференциального уравнения (307) имеет вид

                                 (308)

где

 - функции Крылова (198), в которых вместо (196) нужно принять

Для определения постоянных

, входящих в общее решение (308), необходимо использовать граничные условия. Рассмотрим два случая, которые не освещались при расчете на свободные колебания.

1.Возмущающая сила

 приложена на конце балки. Поперечная сила в сечении должна быть равна этой силе

и граничное условие принимает вид

где знак «+» соответствует силе, приложенной к правому концу, знак «-» - силе, приложенной к левому концу. Кроме того,

.

2. Возмущающая сила

 приложена в промежуточном сечении балки.

В этом сечении должны выполняться четыре условия сопряжения:

где  а - абсцисса сечения, в котором приложена возмущающая сила; индексы «-» и «+» соответствуют сечениям, расположенным бесконечно близко слева и справа от сечения а.

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

.

Приведенные выше рассуждения представляют собой непосредственное решение задачи. Теперь рассмотрим другой способ - разложение решения в ряд по собственным функциям.

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

дифференциальное уравнение движения приобретает вид

                                    (309)

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


Как и выше, представим
 в виде ряда

                          (310)

Также в виде ряда будем искать решение для прогиба

                            (311)

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

                                        (312)

Эта формула совпадает по записи с (301), выведенной выше для продольных колебаний, но в  (312)
 представляет собой собственные формы задачи о свободных колебаниях балки («балочные функции»). Поэтому здесь также справедлива формула (303), относящаяся к случаю сосредоточенных возмущающих сил.

Учитывая, что каждое слагаемое ряда (310) вызывает движение, описываемое соответствующим слагаемым ряда (311), можно записать уравнение (309) в виде



Разделив обе части на
, получим



Левая часть этого равенства равна
, поэтому



Отсюда получим дифференциальное уравнение для




Общее решение этого уравнения имеет вид

                           (313)

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


Кинематическое возбуждение колебаний


К такому же стандартному уравнению можно привести задачу о вынужденных колебаниях, вызываемых кинематическим способом. Рассмотрим вновь одномассовую систему, но предположим, что причиной колебаний груза m являются колебания точки крепления пружины (рис. 37). Пусть закон движения этой точки задан в виде  f = f(t). Удлинение пружины в текущий момент времени равно x - f, а на груз действует сила упругости пружины - С(x - f), и дифференциальное уравнение движения

,

или

Рис. 37

Произведение Сf(t) можно считать приведенной возмущающей силой F(t), и тогда уравнение движения принимает стандартный вид (81).



Колебания оболочек


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

где

          (340)

          (341)

Величина

 представляет собой энергию растяжения оболочки,
 - энергию ее изгиба;
 - компоненты деформации срединной поверхности;
 - параметры изменения ее кривизны. Интегрирование в (340) и (341) выполняется по всей срединной поверхности
 оболочки. Величины
,
по известным формулам выражаются через компоненты амплитудного перемещения
 точек оболочки.

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

,

Частота колебаний определяется формулой Рэлея

                                            (342)

Числитель и знаменатель дроби (342) зависят от выбора функции перемещений

. При этом истинные формы собственных колебаний сообщают выражению (342) стационарные значения, а первая собственная форма - минимум. Обозначим характерный размер оболочки и характерную толщину через
 и
 соответственно. Тогда  (342)  можно записать так:

                              (343)

где

Величины

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

Если геометрия оболочки и условия ее закрепления это допускают, то наименьшие значения частот отвечают такому выбору функций

, при котором

Но требование

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



Колебания, перпендикулярные плоскости кольца


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

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

Рис. 72

Установим зависимость моментов от перемещений. Так как задача линейная, то рассмотрим сначала силовые факторы, связанные со смещением х3, а затем - с х4.

Если х3 постоянно по длине окружности, то кольцо смещается как жёсткое целое, и внутренние силы не возникнут. Если х3 изменяется в зависимости от центрального угла по линейному закону

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

,

где GJкр- крутильная жёсткость бруса.

Если при этом отлична от нуля и вторая производная

, то

меняется кривизна бруса и возникает изгибающий момент

,

где J1 - момент инерции сечения относительно центральной оси, лежащей в плоскости кривизны.

Найдём силовые факторы, связанные с поворотом х4. Если х4 постоянно, то происходит осесимметричный изгиб кольца, причём в его сечениях возникает изгибающий момент

.

При переменном по длине повороте х4 соседние сечения поворачиваются друг относительно друга и возникает крутящий момент

.

Суммируя силовые факторы, связанные с перемещениями х3 и х4, ,получаем

                                 (256)

Составим уравнение движения элемента Rd

бруса (рис.73).

 Будем пренебрегать инерцией поворота элемента вокруг своей оси.

Рис. 73

Условие динамического равновесия в направлении нормали к плоскости кольца приводит к уравнению:

.                                               (257)

Сумма моментов относительно нормали к оси элемента:

.                                     (258)

Сумма моментов относительно касательной к оси элемента:

.                                          (259)

Исключая поперечную силу из (257) и (258) и заменяя моменты в полученном уравнении и уравнении (259) их значениями (256), приходим к системе уравнений, в которую входят только перемещения х3 и х4:

            (260)

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

кольца, решение уравнений (260) можно представить в виде

x3 =Acoskj×coswt , x4 = Bcoskj×coswt .                    (261)

Подставляя значения (261) в уравнение движения (260), получим

                  (262)

Из равенства нулю определителя этой системы получим частотное уравнение, корни которого - собственные частоты - таковы:

                                         (263)

Наименьшая отличная от нуля частота соответствует k=2.



 Колебания пластин и оболочек


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

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



Колебания стержней переменного сечения


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

.                                        (238)

Уравнение крутильных колебаний (187) должно быть заменено уравнением

,                                       (239)

а уравнение поперечных колебаний (192) - уравнением

.                                    (240)

Уравнения (238)-(240) при помощи однотипных подстановок

;  
;  
 можно привести к обыкновенным дифференциальным уравнениям для функции

                                                 (241)

                                             (242)

                                           (243)

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

.

Уравнения (241)-(243) в отличие от уравнений, решённых выше, имеют переменные коэффициенты.

Замкнутую форму решений можно получить лишь в отдельных случаях, когда переменные

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



Колебания в плоскости кольца


Рассмотрим круговой брус малой кривизны постоянного сечения с радиусом R осевой линии (рис.71,а). Будем считать груз нерастяжимым. Перемещение центра тяжести поперечного сечения, зафиксированного угловой координатой

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

                                

.                                                        (244)

а                                                                            б

         

                           

Рис. 71

Угол поворота поперечного сечения бруса в процессе движения определяется формулой

.                                             (245)

Изменение кривизны бруса

 равно производной от
 по дуге:

.                                        (246)

Изгибающий момент в поперечном сечении кольца:

.                                     (247)

Теперь составим уравнение движения элемента

 бруса (рис.71,б).

Помимо перечисленных сил, на элемент действует также сила инерции:

,

где

масса единицы длины бруса.

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

.                                   (248)

Равенство нулю суммы проекций всех сил на направление касательной приводит к уравнению:

.                                      (249)

Уравнение моментов имеет вид

.                                                 (250)

Исключим из (248) и (249) нормальную силу N, а поперечную силу Q заменим её значением из (250):

.                         (251)

Подставляя сюда значение M из (247), получим уравнение движения в перемещениях

, и, наконец, исключая один из компонентов перемещения, с помощью условия нерастяжимости (244) придём к уравнению, в которое входит единственная переменная
:

.                (252)

Решение уравнения движения (252) будем искать в виде

;  
.

При этом для

 получается обыкновенное дифференциальное уравнение

,               (253)

.

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

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

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

;  
.                                (254)

Подставляя (254) в (253), устанавливаем, что последнее удовлетворяется тождественно, если

.                                (255)

Формула (255) определяет частоты собственных колебаний кольца в своей плоскости. Значению
 соответствует нулевая частота, так как при
 формулы (254) описывают смещение кольца как жёсткого тела.


Кратные и нулевые корни частного уравнения


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

Рассмотрим эти случаи на примере системы с двумя степенями свободы.

Из соответствующего такой системе частотного уравнения (39) следует, что при выполнении равенства

                  (49)

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

                                                   (50)

один из корней частотного уравнения обращается в нуль.

В качестве примера рассмотрим свободные колебания плоской системы с двумя степенями свободы (рис.23).

Рис. 23

Обозначим через

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

 

После вычисления соответствующих производных уравнения

Лагранжа записываются так:

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

;  
,

тогда полученные дифференциальные уравнения принимают более простую форму

                                             (51)

   Следовательно, инерционные коэффициенты и обобщённые коэффициенты жёсткости в этих уравнениях

 

и условие (49) выполняется,  значит, рассматриваемая система имеет две одинаковые собст­венные частоты колебаний. Для выяснения этого можно было и не привлекать условие (49), так как из уравнений (51) непосредственно следует

.

Вследствие независимости уравнений (51) постоянные интегрирования одного уравнения не связаны с постоянными интегрирования другого уравнения

 

Для определения постоянных

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

Рассмотрим пример системы с одной нулевой собственной частотой (рис.24).

Рис. 24

Обозначим жёсткость вала на кручение через С; моменты инерции дисков относительно про­дольной оси системы - через

и
.



Крутильные колебания валов


Рассмотрим крутильные колебания многомассовой системы (рис. 26), которая является обще­принятой эквивалентной схемой для расчёта крутильных колебаний коленчатых валов. Колен­чатый вал приводится к эквивалентной схеме путём следующих замен: момент инерции заме­няющего диска относительно оси вала должен быть равен моменту инерции колена относи­тельно той же оси, при этом учитывается присоединённая масса шатуна; жёсткость на кручение участка заменяющего вала должна быть равна жёсткости на кручение соответствующего участка коленчатого вала.

Рис. 26

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

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

Обозначим через

 моменты инерции масс дисков относительно продольной оси вала;
- коэффициенты жёсткости участков при кручении;
- углы по­воротов дисков вокруг продольной оси вала (рис. 26,а).

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

на первом участке

;

на втором участке

;

на

- м участке

Уравнения движения удобнее всего составлять прямым способом (рис. 26,б).

                    (65)

Число этих уравнений

совпадает с числом дисков, т.е. с числом степеней свободы сис­темы.

Одним из решений системы (65) является

,                                          (66)

описывающее равномерное вращение вала и дисков как жёсткого целого.


Кроме того, возможно решение, описывающее упругие колебания системы

                                                 (67)

Подставляя (67) в (65), получим

           (68)

Система уравнений (68) содержит
 неизвестных: n амплитуд и частоту колебаний
.

Если преобразовать систему (68) и рассматривать её как однородную систему линейных ал­гебраических уравнений относительно неизвестных амплитуд
, то ненулевое её решение можно получить, как это неоднократно делалось выше, из условия равенства нулю оп­ределителя. Раскрывая определитель, получим частотное уравнение.

Для иллюстрации сказанного запишем систему (68) в преобразованном виде при




Отсюда получим частотный определитель

                       (69)

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

.                   (70)

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

При
 решение частотного уравнения представляет значительные трудности. Но цепная структура уравнений (68) позволяет упростить определение собственных частот при помощи метода последовательных приближений (метода остатков). Суть метода состоит в следующем. Принимая
 и  задаваясь ориентировочным значением
, из первого уравнения системы (68) находят амплитуду
; из второго уравнения системы можно определить амплитуду
, из третьего уравнения - амплитуду
 и, наконец, из предпоследнего уравнения - амплитуду
. Если в последнее уравнение системы (68) подставить вычисленные значения
 и
, то оно, вообще говоря, не будет удовлетворяться вследствие произвольности исходного значения
 (которое, по сути, является первым приближением). Полученное значение левой части (остаток) характеризует меру неточности первого приближения
 и одновременно показывает, в какую сторону нужно изменить расчётное значение
 во втором приближении.

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


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

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

,                  (71)

которые получаются из уравнений (68) после сложения первых i уравнений системы. Соотно­шение (71) выражает равенство крутящего момента в сечении i-го участка вала (левая часть) сумме моментов сил инерции всех расположенных слева дисков (правая часть).

Задаваясь значением
 и принимая
, находим из соотношения (71) для
:

.

Далее из того же соотношения для
:

.

Общая формула имеет вид

.

Процесс продолжается таким образом до
- го уравнения. После определения из него
 можно переходить к последнему уравнению и вычислять его левую часть. Этот результат дол­жен быть равен нулю, так как если сложить все уравнения типа (71), то должно получиться

.

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

После нескольких расчётов такого типа (при разных значениях
) можно построить кривую зависимости остатка R от
 (рис.27).



Рис. 27

Точки пересечения кривой с осью абсцисс соответствуют истинным значениям частот.

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

При записи решения (67) предполагалось, что колебания являются одночастотными, т.е. для любого диска описываются одной гармоникой

.

Существование спектра частот
 требует обобщения решения (67) и записи его в виде

  
,

где первый индекс у амплитуды означает номер диска, а второй индекс - номер соответствую­щей частоты.

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


добавить слагаемое вида (66), и тогда общее решение уравнений движения (65)

 (72)

Уравнения (72) содержат 2n неизвестных:
 неизвестных амплитуд колебаний первого диска
;
 неизвестных начальных фаз
; угловое смеще­ние
 и угловую скорость
. Амплитуды колебаний всех остальных дисков
 определяются через амплитуды составляющих колебаний первого диска
; отношения
 зависят от но­мера частоты k и определяют соответствующие формы колебаний.

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

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


Маятник с колеблющейся точкой подвеса


Рассмотрим маятник (рис.61,а). Если точка подвеса неподвижна, то единственной силой, создающей момент относительно точки подвеса, является вес груза

.

 Уравнение малых колебаний маятника имеет вид

Если же точка подвеса колеблется вдоль оси y по закону

то при составлении уравнения моментов нужно учесть переносную силу инерции

 момент которой составляет
, и тогда уравнение колебаний маятника запишется так:

или

                                          

                                    (160)

Рис. 61

Уравнение (160) можно привести к стандартному виду - уравнению Матье.

 Для этого положим:

Теперь из диаграммы Айнса-Стретта видно, что параметр а не зависит от амплитуды колебаний точки подвеса и сколь бы малой ни была амплитуда

, неустойчивость нижнего положения маятника наступает вблизи значений
 т.е. при

Рассмотрим вопрос об устойчивости верхнего положения маятника (рис.61,б). При неподвижной опоре это положение, конечно, неустойчиво; однако вибрации основания могут придать ему устойчивость. Для получения уравнения движения в данном случае достаточно изменить знак перед членом, содержащим ускорение

 в уравнении (160); соответственно параметр
 становится отрицательным:

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

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



Метод Граммеля


Принимая форму колебаний подобной статическим прогибам системы от некоторой подходящей нагрузки, можно существенно увеличить точность расчёта за счёт исключения операции дифференцирования. Еще большая точность достигается в методе Граммеля, в котором дифференцирование заменяется интегрированием. Последовательность операций здесь такова:

1. Задают форму колебаний и вычисляют максимальную кинетическую энергию движения:

.

2. Определяют максимальные силы инерции:

.

3. Определяют внутренние силы в элементах системы, вызываемые нагрузками Fi.

4. По внутренним силам вычисляют максимальную потенциальную энергию деформации П0.

5. Из равенства Кmax = П0 определяют частоту колебаний.

Применим метод Граммеля для вычисления частоты колебаний консольной балки. Принимая x = (z/

)2, находим

.

Интенсивность сил инерции:

.

Поперечная сила в сечении:

.

Изгибающий момент:

.

Потенциальная  энергия деформации:

.

Приравнивая Kmax = По, находим

,

что отличается от точного решения на 0,42 %.



Метод последовательных приближений


Докажем, что обычный процесс последовательных приближений приводит к первой собственной форме колебаний. Основой процесса является сравнение двух кривых аn и аn+1, из которых вторая получается как линия прогибов, вызванных нагрузкой man; при этом приближенное значение квадрата частоты определяется по формуле

.                                             (280)

Подобно выражению (265) представим исходную кривую

 в виде ряда

                           (281)

Тогда нагрузка, соответствующая прогибам

, такова:

                              (282)

Рассмотрим одно из слагаемых этой нагрузки -

. От нагрузки
 прогибы будут
, поэтому от нагрузки
 прогибы будут в
 раз больше, т.е. составят
. Следовательно, кривая прогибов от суммарной нагрузки определяется рядом

                     (283)

который отличается от ряда (282) тем, что каждый член ряда разделен на квадрат соответствующей частоты. Так как

 то кривая
 ближе к
, чем исходная кривая
; члены ряда, содержащие
,
 и искажающие основную форму
, представлены в ряде (283) слабее, чем в ряде (281). Продолжая процесс дальше, получим для

 кривой

.                (284)

Как видно, при

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

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

Однако несколько видоизменяя метод, можно добиться того, что в результате последовательных приближений «очистится» не первая, а именно вторая собственная форма колебаний. Этот прием нашел практическое применение при расчете изгибных колебаний крыльев самолетов и лопаток турбин.

Прием основан на устранении формы

 из исходной функции
.
Допустим, что в разложении (281) отсутствует слагаемое, соответствующее первой форме, тогда оно не сможет возникнуть при всех последующих операциях, и ряд (284) принимает вид



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

                                    (285)

где
 - «подходящая» функция;
 - предварительно найденная первая собственная форма.

Коэффициент
 следует принять таким, чтобы форма
 была ортогональна первой собственной форме
:

         


Подставляя сюда (285), получим



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

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

                                  (286)

где  коэффициент
 также определяется условием ортогональности

функций
 и
:



которое после подстановки  (286) дает



Затем следует определить кривую
 от нагрузки
, вновь исправить ее по формуле:

  и т.д.

В таком процессе последовательных приближений ортогонализация сопровождает каждый шаг выкладок и, непрерывно вытесняя «примесь» первой формы, приведет ко второй собственной форме и второй частоте, которая, подобно (280), определится формулой



Таким же образом при помощи сопровождающей ортогонализации можно определить третью собственную форму и третью частоту и т.д.


Метод Ритца


Зададимся несколькими функциями f1(x) , f2(x) ,... , fn(x), каждая из которых удовлетворяет геометрическим граничным условиям задачи, и образум функцию f(x) как сумму

f(x) = C1f1(x) +C2f2(x)+...+Cnfn(x).                       (272)

Если эту функцию подставить в формулу Рэлея

,                                      (273)

то результат будет зависеть от конкретного выбора коэффициентов С1 , С2 , …, Сn.

Метод Ритца основан на простой идее: коэффициенты С1 , С2 , …,  Сn должны быть выбраны так, чтобы вычисление по (273) дало наименьшее значение для w2. Из теоремы Рэлея вытекает, что такой выбор будет наилучшим (при данной системе функций fi).

Условия минимума w2 имеют вид

,                 (i = 1 , 2,..., n),

т.е.

.

Разделив это уравнение на интеграл

 и учитывая (273), получим

             (i = 1 , 2,..., n).            (274)

Уравнения (274) однородны и линейны относительно С1 , С2 , …, Сn и их число равно числу членов выражения (272). Приравнивая нулю определитель, составленный из коэффициентов при С1 , С2 , …, Сn, получим частотное уравнение. Это уравнение не только дает хорошее приближение для низшей частоты, но также определяет (хотя и с меньшей точностью) значения высших частот; при этом можно будет вычислить столько частот, сколько слагаемых принято в выражении (272).

Метод Ритца, как и метод Рэлея, позволяет решить задачу в случаях разрывных функций EJ и m и когда эти функции представлены различными аналитическими выражениями на различных участках.

Иногда та же идея используется в иной форме. Например, при исследовании поперечных колебаний турбинных лопаток задаются функцией f(x) = axs  (начало координат в закрепленном конце). Применяя затем формулу Рэлея (273), получают частоту в виде зависимости от показателя степени s. Затем при помощи числовых расчетов определяют значение s, которому отвечает наименьшая частота. Это позволяет достаточно надежно определить как форму, так и частоту колебаний первого тона.

Пример 20. Определить методом Ритца низшую собственную частоту поперечных колебаний консоли переменного сечения, имеющей толщину, равную единице; высота изменяется по линейному закону:

;
;
 

(

- длина консоли).



Непосредственное решение


Для применения первого способа необходимо предварительно разложить периодические возмущающие моменты в ряды Фурье. После этого уравнения (139) решаются несколько раз - отдельно для каждой гармоники возмущения. Это приводит к ряду однотипных частных задач, каждая из которых требует анализа действия возмущающих моментов одинаковой частоты

:

                                             (140)

При этом стационарные колебания будут происходить с частотой возмущения:

                                           (141)

Подставляя (140) и (141) в (139), получим систему алгебраических уравнений

                (142)

Решая эту систему, находим амплитуды вынужденных колебаний

, а затем и крутящие моменты в сечениях вала:
 на первом участке;
на втором и т.д.



Обратный способ


  Отделяем грузы и рассматриваем упругий безмассовый скелет системы под действием кинетических реакций - сил инерции

 и
(рис. 22,в). В этой схеме первая пружина нагружена силой
, а вторая  - силой
. Перемещение
 конца первой

 пружины, равное её удлинению, можно записать в виде

Перемещение правого конца второй пружины

равно сумме удлинений обеих пружин:

Из этих соотношений получим

 

Таким образом, совпали формы записей дифференциальных уравнений движения по основ­ному (уравнения Лагранжа) и прямому способам, а уравнения, полученные обратным способом, отличаются от них по форме. Это связано с тем, что при нашем выборе обобщённых координат кинетическая энергия имеет каноническую форму:

,

т.е. не содержит произведений скоростей

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

,

то уравнения Лагранжа совпали бы с уравнениями, полученными обратным способом.             

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

 при
, а при составлении по обратному способу
 при
.

Таким образом, пользуясь прямым способом, приходим в общем случае к системе:

            
,                              (33)

а применяя обратный способ - к системе:

            
.                             (34)

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

 называются нормальными, или главными. При этом

,  

и уравнения Лагранжа принимают вид

         
                                  (35)

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



Общее решение стандартного уравнения


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

Идея метода состоит в том, что частное решение уравнения (81) ищется в виде

                                         (82)

соответствующем решению однородного уравнения, но здесь величины С1 и С2 следует считать не постоянными, а переменными. В результате задача определения  функции x(t) заменяется задачей определения двух функций - C1(t) и C2(t). Так как для этого имеется только одно уравнение (82), то функции C1(t) и C2(t) можно связать еще одной произвольной зависимостью.

Составим выражение скорости:

и свяжем

 и
 соотношением

                                   (83)

тогда скорость запишется в форме

а ускорение

            (84)

Подставляя (82) и (84) в (81), получим

                             (85)

Из  (83) и (85) можно найти производные

 и

.

Интегрируя, получим

                           (86)

где B1 и B2 - постоянные величины.

Подставляя (86) в ( 82), получим общее решение уравнения (81)

      (87)

или, внося

 и
 под знаки интегралов и объединяя их,

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

            (88)

Значения постоянных B1 и B2 можно определить только после того, как указаны начальные условия движения. Если

 при t=0, то из  (87) и  (88) найдем

Тогда решение принимает вид

Здесь первые два слагаемых описывают свободные колебания, вызванные начальными возмущениями x0 и V0, а третье слагаемое характеризует вынужденные колебания, вызванные действием возмущающей силы F(t).

В случае нулевых начальных условий, когда движение начинается при x0 = 0 и V0 = 0,

                                    (89)

В некоторых случаях удобнее использовать другую форму решения, которую получим, интегрируя по частям решение (89).

Положим

Тогда

Заменяя

, по формуле интегрирования по частям получим

     (90)

Если в начальный момент времени F(0)=0, то решение принимает вид

                              (91)

где

 - переменное «статическое» перемещение, вычисляемое в предположении, что силы инерции отсутствуют.

Применим полученные результаты к случаю кинематического возбуждения колебаний (рис.37). Полагая F(t)=Cf(t), основное решение (89) запишем в виде

Аналогично вместо формулы (90) при F(0)=0 получим

                               (92)



Общие положения и примеры автоколебательных систем


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

Способность системы так регулировать отбор энергии от источника, чтобы появлялись автоколебания, определяется некоторыми ее свойствами. Этими свойствами являются динамические характеристики - собственные частоты и формы колебаний и коэффициенты затухания, определяющие «динамическую индивидуальность» системы.

Приведем несколько примеров возникновения автоколебаний. Во многих случаях самовозбуждение колебаний возникает вследствие взаимодействия системы с потоком жидкости, доставляющим энергию, необходимую для поддержания колебаний системы. Так, если деревянный стержень полукруглого сечения, подвешенный на пружинах, обдувать потоком воздуха, то стержень будет совершать колебания в вертикальной плоскости, перпендикулярной к направлению потока, причем с течением времени эти колебания становятся весьма интенсивными. Такой стержень полукруглого сечения хорошо моделирует условия реальной проблемы - «галопирование» линий электропередачи. При некоторых метеорологических условиях (главным образом при обледенении) провода линий электропередачи, протянутые между пилонами, могут колебаться с весьма большими амплитудами и низкими частотами.

Колебания стержня полукруглого сечения называют колебаниями системы с одной степенью свободы. Это означает, что движение совершается по одной форме, и соответствующее перемещение системы «стержень - пружины» можно охарактеризовать одной переменной величиной, например вертикальным отклонением центра тяжести стержня от своего среднего положения.
Возможен такой тип автоколебаний, которые могут существовать лишь в случае возбуждения более чем одной степени свободы. В этом случае все положения, которые принимает система в процессе колебаний, могут быть описаны лишь при помощи более чем одной переменной величины.

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

Флаттер может возникать также у лопастей вертолетов и лопаток турбомашин. На практике явление флаттера служит объектом обширных исследований; условие отсутствия флаттера часто является одним из основных требований при проектировании.

Возникновение флаттера связано с выбором «формы флаттера», для которой выполняются определенные соотношения между амплитудами и фазами колебаний, соответствующих различным степеням свободы. Условие возникновения флаттера зависит от скорости потока, а также от плотности и температуры воздуха. Предположим, что скорость изменяется. От скорости потока зависит значение энергии, получаемой системой за один цикл колебаний, и значение энергии, рассеиваемой за цикл колебаний вследствие внутреннего и аэродинамического демпфирования. Когда отношение этих значений энергии становится равным единице, в системе могут установиться колебания постоянной амплитуды; соответствующая скорость самолета называется критической скоростью флаттера.


Каждой из возможных форм флаттера соответствует своя критическая скорость, и все расчеты флаттера проводятся с целью удостовериться, что наименьшая из критических скоростей с достаточным запасом превышает максимально возможную скорость полета самолета.

Существуют три пути борьбы с флаттером.

Первый способ основан на таком изменении характеристик системы, при котором достигается независимость колебаний, соответствующих различным степеням свободы, причем демпфирование всех этих различных форм колебаний положительно. Так, можно добиться того, чтобы поворот оси профиля крыла относительно продольной оси сечения слабо зависел от вертикального перемещения оси. Для этого нужно, чтобы ось занимала определенное положение, а распределение массы по сечению профиля удовлетворяло определенному условию. Другой способ заключается в увеличении собственных частот конструкции за счет увеличения отношений «жесткость/масса» отдельных ее частей. Этот способ основан на том, что энергия, получаемая системой при флаттере за один цикл колебаний, почти не зависит от частоты, тогда как энергия, рассеиваемая за один цикл, пропорциональна частоте.

Два рассмотренных метода обычно используются в практике самолетостроения. Третий метод борьбы с флаттером, вызываемым аэродинамическими силами, не всегда эффективен, но часто с успехом используется для устранения других видов автоколебаний. Этот метод заключается в демпфировании системы. Увеличивая трение в системе, склонной к флаттеру (т.е. увеличивая энергию, рассеиваемую за один цикл колебаний заданной амплитуды) можно, как правило, повысить критическую скорость системы.


Общие сведения о свободных колебаниях


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

Известно, что в ряде случаев тело, получившее некоторое начальное возмущение, после удаления причины этого возмущения продолжает совершать колебания. Эти свободные колебания играют важнейшую роль не только в плане проверки системы на резонанс, т.е. на выявление совпадения одной из собственных частот колебаний с частотой действующих на систему постоянно вибрационных нагрузок. Дело в том, что поведение системы при свободных колебаниях характеризует её "динамическую индивидуальность", которая определяет поведение системы при всех других условиях.

Рассмотрим некоторые примеры возникновения свободных колебаний. После того как по  струне рояля ударяет один из молоточков, струна некоторое время продолжает самостоятельно совершать колебания - свободные колебания. Это становится возможным, во-первых, потому, что струна имеет массу и при движении накапливает кинетическую энергию, а во-вторых, потому, что при отклонении от положения равновесия струна накапливает потенциальную энергию .    

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

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

Идеальным объектом исследования свободных колебаний может служить подвешенная за один конец велосипедная цепь. Пусть свободно висящая цепь первоначально находится в состоянии покоя. Свободные колебания можно вызвать, если отклонить цепь каким-либо образом, а затем отпустить или резко её ударить (но так, чтобы боковое перемещение любой точки цепи было малым по сравнению с длиной цепи).
При этом можно наблюдать следующее:

1. Развитие движения во времени зависит от того, как оно началось.

2.     Движение постепенно затухает.

3.     При своем движении цепь не имеет какой-либо определённой формы; с течением времени форма цепи изменяется, однако в конце движения колебания часто характеризуются более или менее отчетливой формой.

4.     Невозможно указать "частоту" колебаний, но с течением времени движение может принять определённую частоту.

   Кажущийся хаотический характер колебаний можно упорядочить, если надлежащим образом задать начальные условия движения. В этом случае удаётся получить колебания цепи с неизменной формой и определённой частотой. Проще всего этого можно добиться следующим  образом. Подвесим цепь к кулисному механизму (рис. 9).



Рис. 9

При вращении кривошипа точка подвеса перемещается в горизонтальном направлении по синусоидальному (т.е. гармоническому) закону. Если при этом регулировать скорость двигателя, приводящего механизм в движение, то можно изменять частоту колебаний точки подвеса. При очень низких частотах цепь будет просто перемещаться из стороны в сторону, сохраняя более или менее вертикальное положение. Однако при возрастании частоты наступает момент, когда цепь начинает интенсивно раскачиваться с частотой, равной частоте возбуждения (рис. 10,а). Эти колебания не являются свободными, так как цепь находится под действием вибрационной нагрузки. Но если внезапно остановить двигатель, т.е. удалить внешнюю нагрузку, то последующие колебания будут свободными и при этом совершенно не такими, как при случайно заданных начальных условиях. Эти колебания будут затухать, но  всё время будут сохраняться форма и частота колебаний цепи. Иначе говоря, характер движения с течением времени остаётся неизменным. Можно показать, что существует целый ряд частот начального возбуждения, для каждой из которых характер движения будет иным, но также неизменным во времени.


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

 рис. 10,б.

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

Таким образом, цепь обладает рядом форм свободных колебаний, которые называются собственными формами. Каждой собственной форме соответствуют определённая частота и скорость затухания колебаний.

Собственные частоты системы (разумеется, любой системы, а не только подвешенной цепи), её собственные формы и скорости затухания являются индивидуальными характеристиками системы; они не связаны ни с какими внешними воздействиями.



Рис. 10

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

Собственные формы колебаний системы обладают следующими свойствами: любая возможная конфигурация системы может быть представлена в виде суперпозиции (наложения) некоторого числа таких конфигураций, каждая из которых соответствует одной собственной форме колебаний. Таким образом, если задать системе некоторую статическую деформацию, а затем освободить систему, то возникнут свободные колебания по всем соответствующим формам; каждое из таких колебаний будет происходить с соответствующей собственной частотой независимо от остальных движений. Любая комбинация двух собственных форм колебаний также является собственной формой.

Теперь перейдём к рассмотрению такой особенности свободных колебаний, как их затухание.

Эффект затухания колебаний объясняется наличием трения; иногда его называют демпфированием. Звук колокола слышен в течение длительного времени после удара, так как нет значительных сил трения, которые привели бы к рассеянию механической энергии за счёт её перехода в тепловую энергию, а рассеяние энергии за счёт излучения звуковых волн происходит весьма медленно.


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

Рассеяние энергии имеет место в любой колебательной системе. Известно, например, что при вибрациях самолёта часть энергии рассеивается в панелях обшивки за счёт трения в заклёпочных соединениях. Значительным демпфированием должны обладать конструкции зданий, что очень важно с точки зрения поведения здания при землетрясении.

Иногда, если это особенно желательно, можно искусственно вводить трение; так, например, на автомобилях устанавливают демпферы колебаний. Прибор, стрелка которого безостановочно колеблется около положения, соответствующего истинному показанию, доставляет большое неудобство в работе. Поэтому, чтобы стрелка пришла в это положение достаточно быстро, вводится демпфирирование. Чрезмерно сильное трение является столь же вредным, как и недостаточное трение, поскольку в этом случае стрелка перемещается в положение отсчёта слишком медленно.

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

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

2.           Материалы с высоким уровнем рассеяния энергии. При ударе по "колоколу", изготовленному из специального сплава меди и марганца, вместо звона слышится глухой стук. В амортизирующих опорах часто используют резину; это связано с её высокими демпфирующими характеристиками.


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

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

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

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

6.           Пенопластовые или резиновые прокладки. Например, электрические лампочки, упакованные в такие прокладки, можно без всякого риска бросать с большой высоты на твёрдый пол.

Таким образом, существуют два вида демпфирования: искусственно вводимое и связанное с естественными силами трения. Если искусственно вводимое трение чаще всего допускает теоретическую оценку, то естественное трение, как правило, не поддается расчёту и должно определяться экспериментально.

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


Общий случай


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

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

Решение этой основной задачи чаще всего ведут одним из двух способов: непосредственного решения или разложения по собственным формам колебаний. Наибольшее распространение в расчётной практике получил второй способ. Помимо этого, иногда используется способ разложения по собственным формам колебаний при сохранении заданного вида периодических нагрузок, т.е. без разложения их на гармонические составляющие.

Особенности каждого из этих способов рассмотрим на примере простейшей двухмассовой системы (рис.48).

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

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

Рис. 48

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

и
:

                             (122)

Решение этой системы, как и решение одного уравнения, состоит из двух частей: решения соответствующей однородной системы и частного решения неоднородной системы (122).

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


Примем частное решение в виде

.                           (123)

После подстановки (123) в (122) получим два уравнения с двумя неизвестными амплитудами
и
:

                   (124)

Решая систему уравнений (124), находим

                           (125)

Знаменатели выражений для
и
 совпадают с левой частью частотного уравнения (42), если заменить в нём букву
 буквой
. Следовательно, если частота возмущения
 совпадает с любой из двух собственных частот
 или
, то знаменатели формул (125) обратятся в нуль, а амплитуды
 и
 станут бесконечно большими (резонанс).

При
 формулы (125) определяют статические отклонения обеих масс,  вызванные силами
 и
:



При
 решения системы (125) стремятся к нулю. Зависимость амплитуды
 от частоты показана на рис.49. Этот график построен для случая:
. В этом случае число резонансов равно двум, что соответствует числу  степеней свободы системы и числу ее собственных частот.

 При помощи (125) можно найти форму вынужденных колебаний, определяемую отношением
 к
. В общем случае эта форма не совпадает ни с одной из собственных форм колебаний и только при резонансах форма вынужденных колебаний совпадает с формой свободных колебаний. Рассмотрим действие двух сил одинаковой частоты, но сдвинутых по фазе:
 и
.



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



а затем решать две задачи: действуют только «синусные» составляющие
 и
; действуют только «косинусные» составляющие
 и
.

В первой задаче получаются уравнения:



а во второй (когда частное решение имеет вид:
):



Решая каждую из этих задач и суммируя результаты, получим решение исходной задачи.

Если возмущающие силы имеют полигармоническую структуру



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

;
 ; 
;

;
 ; 
,

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


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


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

                                 (210)

и использовать свойство ортогональности собственных форм:

   
.

Общее решение (201) запишем так:

.                           (211)

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

  

.                        (212)

Подставляя в правые части уравнений (211) и (212)

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

получим

        
.

Умножая эти выражения на

 и интегрируя по всей длине, имеем

                          (213)

Бесконечные суммы в правых частях исчезли вследствие свойства ортогональности. Из (213) следуют формулы для постоянных

и

                                      (214)

Теперь эти результаты нужно подставить в решение (211).

Снова подчеркнём, что выбор масштаба собственных форм несущественен. Если, например, в выражении собственной формы (209) принять вместо

 величину в
раз большую, то (214) дадут результаты в
раз меньшие; после подстановки в решение (211) эти различия компенсируют друг друга. Тем не менее часто пользуются нормированными собственными функциями, выбирая их масштаб таким, чтобы знаменатели выражений (214) равнялись единице, что упрощает выражения
 и
.




Определение форм и частот колебаний


Для круглой пластины в уравнениях (317) для амплитудной функции

 следует перейти к полярным координатам
. В этих координатах оператор Лапласа имеет вид

Таким образом, уравнение (317) в полярных координатах принимает форму

                         (325)

Решения этих уравнений, соответствующие колебаниям пластины с n

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

После подстановки этого выражения приходим к уравнениям:

                           (326)

                            (327)

Решениями уравнения (326) являются бесселевы функции порядка

 первого
 и второго
 рода; решениями уравнения (327) - модифицированные бесселевы функции
,
. Таким образом, общее выражение амплитудной функции с
 узловыми диаметрами таково:

          (328)

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

 Для сплошной пластинки в выражении (328) равны нулю коэффициенты
 и
 при функциях, стремящихся к бесконечности при
 Граничные условия на внешнем контуре пластинки образуют в этом случае однородную систему уравнений относительно
 и
. Частотное уравнение получается путем приравнивания нулю определителя системы.

В качестве примера рассмотрим колебания свободной сплошной круглой пластинки. В этом случае на контуре должны выполняться условия:

Изгибающий момент определяется формулой

Поперечная сила:

 

Крутящий момент:

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

                     (329)

Учитывая, что

 является решением уравнения
, а
 - уравнения
, находим

При подстановке в (329) вместо

 его выражения

учтем правила дифференцирования функций Бесселя:

В результате приходим к уравнениям

Здесь аргументом всех бесселевых функций является величина

, где
 - радиус пластинки.

Значения

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

 

Если ограничиться формами колебаний без узловых окружностей, то значениям

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



При
 наименьший корень этого уравнения


 и соответствующая частота собственных колебаний



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



Частотное уравнение