Стационарное решение дифференциального уравнения вынужденных колебаний

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

Определение вынужденных колебаний

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

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

Дифференциальное уравнение вынужденных колебаний

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

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

где $x$ — координата; $\delta $ — коэффициент затухания; $<\omega >_0$ — циклическая частота свободных незатухающих колебаний (если $\delta $=0, то $<\omega >_<0\ >$называют собственной частотой колебаний).

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

где $q$ — заряд; $\delta =\frac<2L>$ — коэффициент затухания; $<\omega >_0=\frac<1><\sqrt>$; $U=U_m<\cos \left(\omega t\right)\ >$ — внешнее переменное напряжение.

Уравнения (2) и (3) можно свести к линейному неоднородному дифференциальному уравнению вида:

где $s$ — колеблющийся параметр; $x_0=\frac$ если колебания механические ($x_0=\frac-\ в\ случае\ электрических\ колебаний$).

Решением уравнения (4) является сумма общего решения однородного уравнения и частного решения неоднородного уравнения. Однородное уравнение при этом имеет вид:

Его общее решение:

где $A_0$ — начальная амплитуда колебаний.

Частное решение уравнения (4) в представлено выражением:

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

Установившись, вынужденные колебания происходят с частотой $\omega $ и являются гармоническими. Амплитуда и фаза этих колебаний определяются равенствами (8) и (9), и они зависят от частоты $\omega $.

Резонанс вынужденных колебаний

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

Из выражения (8) видно, что амплитуда имеет максимум. Для нахождения резонансной частоты (частоты при которой $A=max$), следует найти максимум функции $A(\omega )$. Взяв производную $\frac$ и приравняв ее к нулю получим:

Равенство (10) справедливо при:

Получается, что резонансная частота ($<\omega >_r$) равна:

При $<\delta >^2\ll <\omega >^2_0$ резонансная частота совпадает с собственной частотой колебаний $<\omega >_0.$ Подставим вместо частоты правую часть выражения (11) в формулу (8), получим выражение для резонансной амплитуды вынужденных колебаний:

При небольшом затухании колебаний (если $<\delta >^2\ll <\omega >^2_0$) амплитуда при резонансе равна:

где $Q=\frac<<\omega >_0><2\delta >$ — добротность колебательной системы, величина, характеризующая резонансные свойства колебательной системы. С увеличением добротности увеличивается амплитуда резонанса.

Примеры задач с решением

Задание. Какова добротность колебательного контура, представленного на рис.1?

Решение. Добротность электрического колебательного контура найдем как:

При этом собственная частота колебаний в таком контуре равна:

коэффициент затухания находим как:

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

Ответ. $Q=10$

Задание. Пружинный маятник выполняет вынужденные колебания в вязком веществе. Масса груза на пружине равна $m$, коэффициент упругости пружины $k$. Коэффициент сопротивления среды равен $r$. Систему заставляет совершать колебания сила $F=<\cos \left(\omega t\right)(Н).\ \ \ >$Чему равна резонансная амплитуда заданных колебаний ($A_r$)?

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

где коэффициент затухания равен $\delta =\frac<2m>$. Из функции, которая задает вынуждающую силу:

мы видим, что амплитуда силы равна единице:

Собственная частота колебаний груза на пружине:

Амплитуда при резонансе таких колебаний равна:

Лекция №8. МЕХАНИЧЕСКИЕ КОЛЕБАНИЯ

5.6. Затухающие гармонические колебания.

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

где μ − постоянная, называемая коэффициентом сопротивления.

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

Учитывая , что a= $$d^2x\over dt^2$$ , а υ= $$dx\over dt$$ и разделив на массу m , получим

Применив обозначения $$ = ω_0<^2>$$ , $$ <μ\over m>= 2β$$ и $$ = f_0$$ получим

дифференциальное уравнение затухающих колебаний . Отметим, что ω₀ представляет собой ту частоту, с которой совершались бы свободные колебания системы в отсутствие сопротивления среды. Эта частота называется собственной частотой .

Для решения уравнения (5.6.4) сделаем подстановку

Проведем замену переменных

Подставим (5.6.5 и 5.6.6) в выражение (5.6.4)

Преобразуем , сократив на e -βt

Рассмотрим случай, когда сопротивление среды настолько мало, что ω₀ 2 -β 2 >0 есть величина положи мы можем ввести тельная, и обозначение ω₀ 2 -β 2 =ω 2 , после чего уравнение (5.6.8) примает вид

В случае большого сопротивления среды ω₀ 2 -β 2 , движение становится непериодическим.

Решение уравнения (5.6.8) можно записать в виде

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

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

и амплитудой, изменяющейся по закону

На рисунке показан график данной функции. Пунктирными линиями показаны пределы, в которых находится смещение колеблющейся точки. Верхняя из пунктирных кривых дает график функции A(t) , причем величина A₀ представляет собой амплитуду в начальный момент времени. Начальное смещение зависит от A₀ и также от начальной фазы φ , т.е. x₀=A₀cosφ .

5.7. Коэффициент затухания и логарифмический декремент затухания.

Отношение значений амплитуд, соответствующих моментам времени, отличающимся на период, равно

и называется декрементом затухания .

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

Скорость затухания колебаний определяется величиной называем коэффициентом затухания $$β=<μ\over 2m>$$ .

Найдем время, называемое временем релаксации τ , за которое амплитуда уменьшается в e раз

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

За время релаксации τ система успевает совершить $$N_e=<τ\over T>$$ колебаний

Следовательно, $$δ=<1\over N_e>$$ логарифмический декремент затухания обратно пропорционален по величине числу колебаний, за которые амплитуда колебаний уменьшается в e раз.

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

которая называется добротностью колебательной системы.

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

5.8. Вынужденные колебания.

До сих пор мы рассматривали свободные колебания, когда выведенная из положения равновесия система совершает колебания будучи предоставленной самой себе. Рассмотрим колебательную систему, которая подвергается действию внешней силы, изменяющейся по гармоническому закону F=F₀cosωt . Колебания, возникающие под действием внешней периодически изменяющейся силы, называются вынужденными колебаниями . В этом случае уравнение второго закона Ньютона имеет вид

Учитывая , что a= $$d^2x\over dt^2$$ , а υ= $$dx\over dt$$ и разделив на массу m , получим

Применив обозначения $$ = ω_0<^2>$$ , $$ <μ\over m>= 2β$$ и $$ = f_0$$ получим

дифференциальное уравнение вынужденных колебаний.

Будем искать решение уравнения (5.8.3) в виде

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

Подставим (5.8.4) и (5.8.5) в уравнение (5.8.3)

Чтобы уравнение (69) обратилось в тождество необходимо, чтобы коэффициенты при cosωt и sinωt были равны нулю.

Из выражения (71) получаем

Возведем в квадрат уравнения (70) и сложим

Подставив полученные выражения (71) и (73) в выражение (64) получим уравнение вынужденных колебаний

5.9. Резонанс.

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

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

Найдем резонансную частоту. Амплитуда вынужденных колебаний будет max, когда выражение $$(ω_0<^2>-ω^2)^2 + 4β^2ω^2$$ в уравнении $$A=-ω^2)^2 + 4β^2ω^2>>$$ (5.8.13) будет минимальным.

Продифференцируем это выражение по ω и приравняем к нулю

Полученное уравнение имеет три решения: ω=0 и ω=± $$\sqrt <ω_0<^2>-2β^2>$$ . 2 . Первое решение соответствует максимуму знаменателя. Из остальных двух решений отрицательное не имеет физического смысла (частота не может быть отрицательной). Таким образом, резонансная циклическая частота

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

Из последнего уравнения (5.9.3) следует, что при отсутствии сопротивления среды амплитуда при резонансе обращалась бы в бесконечность, а резонансная частота, согласно (5.9.2), при тех же условиях (при β=0 ), совпадала бы с собственной частотой колебаний системы ω₀

Зависимость амплитуды вынужденных колебаний от частоты вынуждающей силы показана графически на рис. 5.9.1. В соответствии с (5.9.2) и (5.9.3), чем меньше параметр β , тем выше и правее лежит максимум данной кривой. Изображенная на рис. 5.9.1 совокупность графиков функций (5.8.13), соответствующих различным значениям параметра β , называется резонансными кривыми .

При стремлении ω к нулю все кривые приходят к одному и тому же, отличному от нуля, предельному значению, равному f₀ω₀ 2 . Это значение представляет собой смещение из положения равновесия, которое получает система под действием постоянной силы величины F₀

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

Наконец, отметим, что чем меньше β , тем сильнее изменяется с частотой амплитуда вблизи резонанса, тем «острее» получается максимум. При малом затухании (т. е. β ) амплитуда при резонансе приближенно равна A_pes≈f₀/2βω₀ . Разделим это выражение на смещение x₀ из положения равновесия под действием постоянной силы F₀ , равное x₀=f₀/ω_p 2 . В результате получим

где δ = βТ – логарифмический декремент затухания (5.7.2); Q – добротность колебательной системы (5.7.6).

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