Энергетический подход к нахождению периода колебаний

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

Рассмотрим одномерное движение частицы массой m в поле, описываемом потенциальной энергией U(x). Динамическое уравнение движения (второй закон Ньютона) для частицы имеет вид

Энергетический подход к нахождению периода колебаний - student2.ru ,

где Энергетический подход к нахождению периода колебаний - student2.ru - ускорение частицы a (вторая производная от координаты по времени, а = Энергетический подход к нахождению периода колебаний - student2.ru ); F – сила, действующая на частицу, определяется видом потенциальной энергии и в данном случае

F = - Энергетический подход к нахождению периода колебаний - student2.ru (93)

Чтобы найти зависимость координаты от времени, надо проинтегрировать уравнение движения, которое является дифференциальным уравнением второго порядка. Однако в случае одномерного движения нет необходимости это делать. Так как потенциальная энергия U(x) не зависит от времени t, то энергия частицы должна сохраняться: E = const. Энергия частицы определяется суммой кинетической и потенциальной энергий

E = Энергетический подход к нахождению периода колебаний - student2.ru + U(x), (94)

где Энергетический подход к нахождению периода колебаний - student2.ru Энергетический подход к нахождению периода колебаний - student2.ru - скорость частицы (первая производная от координаты по времени):

Энергетический подход к нахождению периода колебаний - student2.ru .

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

Энергетический подход к нахождению периода колебаний - student2.ru + Энергетический подход к нахождению периода колебаний - student2.ru = Е

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

m x" + kx = 0. (95)

Для электрического колебательного контура без потерь (сопротивление R = 0) закон сохранения энергии можно представить в виде Энергетический подход к нахождению периода колебаний - student2.ru , i = Энергетический подход к нахождению периода колебаний - student2.ru ,

где первое слагаемое – энергия накопленная в магнитном поле катушки индуктивности, второе – энергия накопленная в электрическом поле конденсатора.

После дифференцирования по времени L q" + Энергетический подход к нахождению периода колебаний - student2.ru q = 0

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

Энергетический подход к нахождению периода колебаний - student2.ru + Энергетический подход к нахождению периода колебаний - student2.ru = Е, ω = Энергетический подход к нахождению периода колебаний - student2.ru ,

где первое слагаемое – кинетическая энергия поворачивающегося тела с моментом инерции J, второе – энергия упругой деформации закрученной нити, на которой подвешено данное тело ( r – коэффициент пропорциональности).

После дифференцирования J α" + r α = 0.

Для физического маятника закон сохранения энергии можно представить в виде

Энергетический подход к нахождению периода колебаний - student2.ru + mgh = E, h = l (1- cos α) = 2l sin2 Энергетический подход к нахождению периода колебаний - student2.ru ≈ l Энергетический подход к нахождению периода колебаний - student2.ru , ω = Энергетический подход к нахождению периода колебаний - student2.ru ;

где первое слагаемое - кинетическая энергия поворачивающегося тела с моментом инерции J, второе – потенциальная энергия тела, чей центр масс в результате поворота поднят на высоту h.

После дифференцирования

J α" + mgl α = 0.

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

Энергетический подход к нахождению периода колебаний - student2.ru . (96)

Откуда, разделяя переменные и интегрируя, получим

t - t0 = ± Энергетический подход к нахождению периода колебаний - student2.ru . (97)

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

U(x) < E.

Заметим, что подынтегральное выражение в (97) при этом всегда положительно.

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

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

E = U(x), (98)

определят границы движения. Они являются точками остановки, поскольку в них скорость обращается в нуль. Более того, можно ут-

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

Если допустимая область движения ограничена точками остановки с обеих сторон (между точками х1 и х2 на рис. 23.), то движе

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

U(x)
Энергетический подход к нахождению периода колебаний - student2.ru

Рис. 23. График зависимости потенциальной энергии частицы от координаты и точки заворота осциллятора

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

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

T(E) = 2 Энергетический подход к нахождению периода колебаний - student2.ru . (99)

Причем приделы интегрирования х1 и х2 являются корнями уравнения (94) при данном значении E. Период колебаний частицы является функцией ее энергии. Вид этой функциональной зависимости определяется видом потенциальной энергии U(x).

Рассмотрим случай малых колебаний, которые система совершает в области своего положения устойчивого равновесия. В положении устойчивого равновесия функция U(x) имеет минимум. Условимся координату х и потенциальную энергию U(x) отсчитывать от положения равновесия. Тогда U(0) = 0.

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

U(x) = U(0) + U'(0) · x + Энергетический подход к нахождению периода колебаний - student2.ru U"(0) x2 + … . (100)

Ввиду малости х остальными членами пренебрегаем. Они важны в нелинейных процессах.

Так как U(x) при х = 0 имеет минимум, то U'(0) = 0, а вторая производная U"(0) положительна. Кроме того, U(0) = 0. Тогда

U(x) = Энергетический подход к нахождению периода колебаний - student2.ru ;

где введено обозначение

k = U"(0).

Сила, действующая на частицу, в этом случае

F = - Энергетический подход к нахождению периода колебаний - student2.ru = - k x,

и второй закон Ньютона

m Энергетический подход к нахождению периода колебаний - student2.ru = - k x.

Отсюда находим x(t) в виде гармонического колебания, введя обозначение

Энергетический подход к нахождению периода колебаний - student2.ru .

Добротность

Добротность – очень емкое понятие в теории колебаний. Дело в том, что любая одномерная колебательная система всегда характеризуется двумя основными параметрами: w0 – собственной частотой и b - коэффициентом затухания. Следует отметить, что понятие коэффициента затухания вводится при условии, что сила трения, действующая на осциллятор, пропорционально скорости Fсопр ~ rV . Если закон сопротивления другой, как, например, при колебаниях с сухим трением, то введение b затруднительно. Пусть указанные два параметра ввести можно; тогда можно ввести и их отношения w0/b0 - мы получим безразмерный параметр, который и называется добротностью. Иногда в частных исследованиях вводят величины, отличающиеся от w0/b на некоторый постоянный множитель. Но это ничего не меняет.

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

х¢¢ + 2b х¢ + w02 х = 0.

Здесь на осциллятор действует только собственная возвращающая сила Fв= -w02x и сила трения Fт = -2bx’, х(t) – любая колеблющаяся физическая величина: смещение частицы от положения устойчивого равновесия, сила тока в колебательном контуре, смещение столбика газа в акустическом резонаторе и т.д.

Решением уравнения (110) при условии w0 > b является затухающая гармоника (рис. 26)

х(t) = Ае -bt cos (wt + j),

где w2 = w02 - b2, а величина а = Ае -bt может рассматриваться как переменная во времени амплитуда. Пусть t - время, за которое амплитуда уменьшается в е раз, т.е.

А/а = е = е-bt

за это время осциллятор успевает совершить Ne колебаний

Ne = t/Т = 1/(bТ); Т = 2p/w.

Введем параметр

Энергетический подход к нахождению периода колебаний - student2.ru ,

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

1. Рассмотрим убыль энергии осциллятора при затухающих колебаниях. Запас механической энергии

E = Энергетический подход к нахождению периода колебаний - student2.ru ( к – жесткость пружины)

Отсюда

E = Энергетический подход к нахождению периода колебаний - student2.ru (101)

(рис. 24).

здесь Е0 – начальный запас энергии в системе. Продифференцируем (101) по времени и определим скорость убыли энергии осциллятора:

Энергетический подход к нахождению периода колебаний - student2.ru . Энергетический подход к нахождению периода колебаний - student2.ru

Будем считать период колебаний Т достаточно малым, таким, что Т ~ dt, что оправдано для высокодобротных колебаний, тогда изменение энергии за период DЕ из (101):

Энергетический подход к нахождению периода колебаний - student2.ru Энергетический подход к нахождению периода колебаний - student2.ru .

Таким образом, добротность осциллятора в 2p раз больше убыли его энергии за период (рис. 24).

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

внешняя вынуждающая сила F = F0 cos wt будет гармонической. Уравнение вынужденных колебаний возьмем в виде

Энергетический подход к нахождению периода колебаний - student2.ru ,

где f = Энергетический подход к нахождению периода колебаний - student2.ru - амплитуда (приведенная) внешней (вынуждающей) силы, w - ее частота.

Решение этого уравнения в виде установившихся колебаний (без учета переходного процесса) имеет вид (86 - 87):

х(t)=Аcos(wt+j), A(w) = Энергетический подход к нахождению периода колебаний - student2.ru Энергетический подход к нахождению периода колебаний - student2.ru

Функция А(w) – амплитуда установившихся колебаний, она имеет максимум на резонансной частоте

Энергетический подход к нахождению периода колебаний - student2.ru .

На этой частоте амплитуда колебаний (резонансная) имеет значение

Энергетический подход к нахождению периода колебаний - student2.ru . (102)

Эти соотношения уточняют приведенные ранее.

Из (102) находим, что на частоте близкой к нулю

A(0) = Энергетический подход к нахождению периода колебаний - student2.ru

(рис. 27)

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

K = Энергетический подход к нахождению периода колебаний - student2.ru (103)

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

3. Рассмотрим сдвиг фаз между вынуждающей силой и установившимися колебаниями, задаваемый формулой (87). Преобразуем ее к виду

Энергетический подход к нахождению периода колебаний - student2.ru

Отсюда при

Энергетический подход к нахождению периода колебаний - student2.ru

имеет место

Энергетический подход к нахождению периода колебаний - student2.ru Энергетический подход к нахождению периода колебаний - student2.ru

Таким образом, добротность равна котангенсу сдвига фаз на некоторой частоте w = 0, 618033988w0 (между вынуждающей силой и установившимися колебаниями)

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

A= f0/( Энергетический подход к нахождению периода колебаний - student2.ru )= (1/ k)Арез =(1/ k)f0/2β Энергетический подход к нахождению периода колебаний - student2.ru (104)

Отсюда

20 – ω2)2 + 4β2 ω2 = 4 β2k220 - β2).

Возведем это выражение в квадрат и получим биквадратное уравнение относительно ω :

ω4 + ω2 (4 β2 - 2 ω20 ) + ω40 - 4 β2k2 ω20 +4 β4k2 = 0.

Пусть корни этого уравнения ω12 и ω22 . Для дальнейшего анализа применим теорему Виетта. Имеем

ω12 + ω22 = - (4 β2 - 2 ω20 )

ω12 * ω22 = ω4 0 - 4 β2k2 ω20

Членом, содержащим β4, мы пренебрегли ввиду его малости.

Составим выражение, содержащее разность частот ω1 и ω2 :

1 + ω2 )2 = ω12 + ω22 – 2 ω1* ω2 =

= - 4 β2 + 2 ω20 - 2 Энергетический подход к нахождению периода колебаний - student2.ru .

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

Энергетический подход к нахождению периода колебаний - student2.ru

Окончательно имеем:

Энергетический подход к нахождению периода колебаний - student2.ru . (105)

Таким образом, относительная ширина частотной полости (по уровню 1/k) Dw/w0 обратно пропорциональна добротности осциллятора. Чем выше добротность, тем острее резонансная кривая (рис. 24).

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

Энергетический подход к нахождению периода колебаний - student2.ru

а вынужденные - на частоте вынуждающей силы w близкой к w0

Энергетический подход к нахождению периода колебаний - student2.ru

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

Энергетический подход к нахождению периода колебаний - student2.ru (106)

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

Энергетический подход к нахождению периода колебаний - student2.ru Энергетический подход к нахождению периода колебаний - student2.ru Энергетический подход к нахождению периода колебаний - student2.ru (107)

Найдем скорость v (t) и ускорение a(t), продифференцировав (106)

Энергетический подход к нахождению периода колебаний - student2.ru

Энергетический подход к нахождению периода колебаний - student2.ru

Подставив сюда начальные условия (107), найдем амплитуды A и В (постоянные интегрирования) и получим

B = -A; Энергетический подход к нахождению периода колебаний - student2.ru

Тогда общее решение (80) примет вид

Энергетический подход к нахождению периода колебаний - student2.ru

Это выражение представляет собой биения.

Устремим теперь частоту вынуждающей силы w к частоте собственных колебаний осциллятора w0

Энергетический подход к нахождению периода колебаний - student2.ru

(рис. 24).

Таким образом, амплитуда колебаний при раскачке осциллятора в нашем случае происходит по линейному закону в зависимости от времени. Но раскачка заканчивается, когда амплитуда достигает значения, задаваемого формулой (75). Пусть раскачка продолжается в течение времени t, тогда

Энергетический подход к нахождению периода колебаний - student2.ru

отсюда

Энергетический подход к нахождению периода колебаний - student2.ru (108)

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

Порядки величин добротностей некоторых осцилляторов даны в табл. 3.

Таблица 3

Добротность осцилляторов

Механические (резина) 10 на 100 Гц; 5 на 2*103 Гц; 3 на 104 Гц;

Акустические резонаторы – 50;

Радиотехнические контуры с частотой ~ 1 МГц – несколько сотен;

Медные резонаторы на СВЧ f >1МГц – 3*104;

Пьезоэлектрические кристаллы (кварц) – 5*105 ;

Колебания ядер атомов в эффекте Мессбауэра – 1010;

Лазерные колебания (резонатор Фабри-Перо) – 5*106 ;

Пульсары (нейтронные замагниченные звезды) – 7,5*1012

(Период 0,033с, замедление вращения и частоты следования импульсов 36,52 нс/сутки = 4,23*10-13. Период увеличивается в е раз за 2500 лет).

Энергетический подход к нахождению периода колебаний - student2.ru

Рис. 24. Добротность как характеристика различных сторон колебательных процессов

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

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

Рассмотренные затухающие колебания характеризуются двумя параметрами β и ω (или ω0). В ряде случаев весьма информативным оказывается их отношение ω/β, называемое добротностью. Аналогично вводится логарифмический декремент затухания. В различных разделах физики и техники определение добротности свое и на это следует обращать внимание. Характеристика различных аспектов колебательных процессов с помощью понятия добротности представлена в таблице 3.

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

Таблица 4

Основные колебательные уравнения и их решения

Гармонические колебания
х 11 + ω02х = 0 x(t) = Acos(ω0 t - j0)
Гармонические колебания около смещенного положения устойчивого равновесия. Реализуются, когда на колеблющуюся материальную точку действует постоянная внешняя сила
х 11 + ω02х = fc x(t) = х0 + Acos(ω0 t - j0)
Гармонические колебания при силе сопротивления пропорциональной скорости
х11 + 2βх1 + ω02х = 0 (β < ω0) x(t) = A0e-βt cos(ωt + j0). ω2= ω20 – β2
Гармонические колебания при силе сопротивления пропорциональной скорости и при действии постоянной силы
х11 + 2βх1 + ω02х = fc (β < ω0) x(t) = х0 + A0e-βt cos(ωt + j0). ω2= ω20 – β2
Гармонические колебания при постоянно действующей внешней силе, которая линейно зависит от времени t.
х11 + 2βх1 + ω02х = a + bt x(t) = A0e-βt cos(ωt + j0)+α+γt, где ω2= ω20 – β2, γ = b/ ω02, α = a/ ω02 - 2 β b/ ω04
Колебания под действием внешней гармонической силы. Решение в виде установившихся колебаний
х11 + 2βх1 + ω02х = f0 cosωt x(t) = А cos(ωt - j), где A= f0/( Энергетический подход к нахождению периода колебаний - student2.ru ), tg j = 2βω/(ω02 – ω2).
х11 + ω02х = f0 sinωt x(t) = Энергетический подход к нахождению периода колебаний - student2.ru
Колебания под действием внешней силы, представимой в виде ряда Фурье. Решение в виде установившихся колебаний
х11 + 2βх1 + ω02х = = Энергетический подход к нахождению периода колебаний - student2.ru f0i cosωi t x(t) = Энергетический подход к нахождению периода колебаний - student2.ru Аi cos(ωi t - ji ), где Ai= f0i /( Энергетический подход к нахождению периода колебаний - student2.ru ), tg ji = 2βωi /(ω02 – ωi 2).
Вынужденные колебания в нелинейной системе
х11 + ω02х +sx3 = f0 cosωt x(t) ≈ Энергетический подход к нахождению периода колебаний - student2.ru ( ω02A + Энергетический подход к нахождению периода колебаний - student2.ru sA3 - f0 ) cosωt + Энергетический подход к нахождению периода колебаний - student2.ru cos3ωt. ( Энергетический подход к нахождению периода колебаний - student2.ru ω02)

Фазовые траектории

Характер движения материальной точки определяется вторым законом Ньютона, т.е. дифференциальным уравнением 2-го порядка, и начальными значениями координат и скорости (всего 6 величин в самом общем случае). Большинство задач на движение можно свести к движению вдоль одной прямой (например, взяв проекции соответствующих векторов), положим, вдоль оси х. Любое движение осциллятора, при котором меняется одна единственная физическая величина, для которой записано дифференциальное уравнение движения, есть одномерное движение.

Введем ортогональную систему координат x и v(x) . Плоскость, заданную этими координатами назовём фазовой. С течением времени материальная точка в своем движении по фазовой плоскости очерчивает фазовую траекторию v(x). Значения x и v в заданный момент времени определяют изображающую точку. Изображающие точки с координатами f(x, v) = 0 и v= 0 называются особыми точками. Особым точкам соответствуют состояния равновесия системы.

Рассмотрим простой осциллятор и установим основные свойства фазовых плоскостей и фазовых траекторий на его примере. Смещение осциллятора задаётся в виде

x(t) = А cos(ω0 t + α),

Скорость осциллятора найдена дифференцированием x(t)

vx(t) = dx/dt = -A ω0 sin (ω0t + α)

Исключим время из выражений x(t) и vx(t) . Для этого выразим cos x и sin x, возведем в квадрат и сложим. Получим

Энергетический подход к нахождению периода колебаний - student2.ru (109)

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

2.1.18 а. Свойства фазовых траекторий

1. Семейство фазовых траекторий, отличающихся начальными условиями, т.е. значениями х0,, v0 , являются не двух-, а однопараметрическим, т.к. определяет не х и v поразнь, а их комбинацию, образующую полную энергию Е осциллятора (А = 1)

Е = Энергетический подход к нахождению периода колебаний - student2.ru (v2 + Энергетический подход к нахождению периода колебаний - student2.ru х2) . (110)

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

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

4. Установим, что движение изображающей точки М (х, v) происходит вдоль фазовой траектории с постоянной секторной скоростью. Понятие секторной (секториальной) скорости (рис. 25) введено в § 49 (т.1) / 1/. Определим ее в векторном виде по правилу

Энергетический подход к нахождению периода колебаний - student2.ru Энергетический подход к нахождению периода колебаний - student2.ru Энергетический подход к нахождению периода колебаний - student2.ru

Здесь Энергетический подход к нахождению периода колебаний - student2.ru - радиус-вектор движущейся точки, Энергетический подход к нахождению периода колебаний - student2.ru - его приращение за время dt (см. рис. 25). Из этого рисунка видно, что векторное произведение Энергетический подход к нахождению периода колебаний - student2.ru есть площадь dS треугольника, образованного векторами Энергетический подход к нахождению периода колебаний - student2.ru и Энергетический подход к нахождению периода колебаний - student2.ru , т.к. dS = Энергетический подход к нахождению периода колебаний - student2.ru sin Энергетический подход к нахождению периода колебаний - student2.ru , где θ – угол между векторами Энергетический подход к нахождению периода колебаний - student2.ru и Энергетический подход к нахождению периода колебаний - student2.ru .

Энергетический подход к нахождению периода колебаний - student2.ru

Рис. 25. К определению секторной скорости

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

Окончательно,

Энергетический подход к нахождению периода колебаний - student2.ru . (111)

Представим радиус – вектор Энергетический подход к нахождению периода колебаний - student2.ru изображающей точки в координатах v,х в виде Энергетический подход к нахождению периода колебаний - student2.ru = х Энергетический подход к нахождению периода колебаний - student2.ru + v Энергетический подход к нахождению периода колебаний - student2.ru , где Энергетический подход к нахождению периода колебаний - student2.ru и Энергетический подход к нахождению периода колебаний - student2.ru - орты (единичные векторы) по осям х и v соответственно. Тогда Энергетический подход к нахождению периода колебаний - student2.ru 1(t) = х1 Энергетический подход к нахождению периода колебаний - student2.ru + v1 Энергетический подход к нахождению периода колебаний - student2.ru = v Энергетический подход к нахождению периода колебаний - student2.ru + v1 Энергетический подход к нахождению периода колебаний - student2.ru . Подставим Энергетический подход к нахождению периода колебаний - student2.ru и Энергетический подход к нахождению периода колебаний - student2.ru 1 в (111). Получим

Энергетический подход к нахождению периода колебаний - student2.ru . (112)

Если х изменяется по гармоническому закону x = А cos(ω0 t + α), vx = -A ω0 sin (ω0t + α), v1 x = d2 x/dt = -A ω20 cos (ω0t + α), то после подстановки в (112) будем иметь

Энергетический подход к нахождению периода колебаний - student2.ru . Выражение в скобках есть ни что иное, как полная механическая энергия материальной точки (110).

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

2.1.18 б.Фазовые траектории затухающих колебаний и релаксаций

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

х11 + 2βх1 + ω02х = 0.

Полагая Энергетический подход к нахождению периода колебаний - student2.ru , находим Энергетический подход к нахождению периода колебаний - student2.ru = - 2 βy - ω02х. Введем новые переменные

u = ω01x, v = y +βx, ω201 = ω02 – β2,

и будем считать, что u и v – прямоугольные координаты.

Тогда из ω201 = ω02 – β2 получим

y2 + 2 βxy + ω02х = (y +βx)2+ ω201x2 = u2+v2,

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

u2+v2 = С Энергетический подход к нахождению периода колебаний - student2.ru ,

(113)

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

ρ2 = С Энергетический подход к нахождению периода колебаний - student2.ru , или ρ = С Энергетический подход к нахождению периода колебаний - student2.ru . (114)

Мы получили уравнение спирали, навивающейся на начало координат.

Если теперь перейти к исходным координатам у',х, т.е. v,x, то в силу линейности преобразования u = ω01x, v = y +βx, интегральные кривые качественно не изменятся и будут также представлять собой спирали (рис. 26). Спирали закручиваются при затухающих колебаниях и разворачиваются при нарастающих.

Фазовые траектории затухающих колебаний, релаксаций и биений представлены на рис. 26 – 28.

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

Энергетический подход к нахождению периода колебаний - student2.ru . (115)

Это уравнение выражает второй закон Ньютона для точки, на которую действует возвращающая сила сложного вида и силы сопротивления, сложным образом зависящие от скорости Энергетический подход к нахождению периода колебаний - student2.ru . Уравнение (115) может быть сведено к системе двух дифференциальных уравнений

Энергетический подход к нахождению периода колебаний - student2.ru = у, Энергетический подход к нахождению периода колебаний - student2.ru = f (x,y) (116)

Решением этой системы является пара функций x(t), y(t), задающих в параметрическом виде фазовую траекторию. Параметр t - это, ко-

нечно, время. Характерной особенностью системы (116) является то, что в нее время явно не входит.

Энергетический подход к нахождению периода колебаний - student2.ru Энергетический подход к нахождению периода колебаний - student2.ru

а б

Рис.26. Фазовые траектории затухающих колебаний в коор-

динатах u, v и х¢, х

. Энергетический подход к нахождению периода колебаний - student2.ru

Рис. 27. Фазовая траектория релаксационных движений

Энергетический подход к нахождению периода колебаний - student2.ru

Рис. 28. Фазовая траектория биений

Если окажется, что в некоторой точке M(x0, y0) фазовой плоскости выполняются равенства

y = 0, f(x,y) = 0 ,(117)

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

При анализе движения на фазовой плоскости стремятся установить:

1) наличие особых точек;

2) различия в расположение фазовых траекторий вблизи особых точек;

3) устойчивость и неустойчивость особых точек, т.е. остается ли с течением времени материальная точка вблизи особой точки фазовой траектории или удаляется от нее;

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

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

Спектры колебаний

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

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

Рассмотрим простейшее сложное колебание – гармонику, модулированную по амплитуде:

x = А cos(ωt - j ) = а(1 + к cosn t) cos(ωt - j )= а cos(ωt - j )+ Энергетический подход к нахождению периода колебаний - student2.ru cos((n+ω)t - j ) + Энергетический подход к нахождению периода колебаний - student2.ru cos((n- ω)t - j ) (118)

Это сложное колебание мы представили в виде набора трех гармоник с несущей частотой ω и двумя боковыми частотами ω+n и ω -n. Зависимость амплитуд гармоник от частоты называется амплитудным спектром. Для данного случая он представлен на рис. 29. Начальные фазы у всех наших гармоник одинаковы, т.е. фазы не зависят от частоты гармоник. Если в более сложных случаях такая зависимость появляется, то говорят о фазовом спектре. Например, если гармоническое колебание получено в результате сложения двух одночастотных колебаний со своими амплитудами и фазами, то можно говорить об амплитудном и фазовом спектрах. Заметим, что поскольку гармоническое колебание можно разложить на два (и больше) таких же бесконечным числом способов, то говорить о спектрах отдельных гармоник бессмысленно.

Рассмотрим сложение одночастотных колебаний с одинаковыми амплитудами и равномерно (по линейному закону) нарастающими фазами. Исходные формулы имеют вид (а = ωt, h = φ)

cosa + cos (a + h) +…..+ cos[a + (n-1)h] = Энергетический подход к нахождению периода колебаний - student2.ru ;

sina + sin(a+h) +….+ sin[a + (n-1)h]= Энергетический подход к нахождению периода колебаний - student2.ru . (119)

Здесь амплитуды гармоник стоящих слева одинаковы, а их частоты равноотстоят друг от друга. Следовательно, спектр функций, стоящих справа в (119) представляет собой «гребенку» (см. рис. 30).

Другой пример – гребенка с амплитудами, нарастающими в арифметической прогрессии:

cosa+2cos2a+…..+n cos na = Энергетический подход к нахождению периода колебаний - student2.ru ;

sina + 2sin2a +….+ n sin na = Энергетический подход к нахождению периода колебаний - student2.ru . (120)

Фазовый спектр также может представлять гребенку (рис. 31).

Примеры разложения процессов в ряды Фурье представлены в таблице в /3/. Приведем несколько примеров.

Для четной периодической функции (рис. 32)

f(ωt) = Энергетический подход к нахождению периода колебаний - student2.ru (121)

Для нечетной периодической функции (рис. 33)

f(ωt) = Энергетический подход к нахождению периода колебаний - student2.ru (122)

Для периодической функции, не обладающей свойством четности и заданной на промежутке – π < β t < π функцией x(t) = exp(β t) (рис. 34), имеем ряд Фурье вида

f(ωt)= Энергетический подход к нахождению периода колебаний - student2.ru , (123)

где k – натуральное число.

Спектр затухающей гармоники f(t) = Энергетический подход к нахождению периода колебаний - student2.ru sin(ω1t) является сплошным и имеет вид:

Энергетический подход к нахождению периода колебаний - student2.ru ; (124)

Энергетический подход к нахождению периода колебаний - student2.ru . (125)

Иногда определяют следующие спектральные функции:

- амплитудный спектр

Энергетический подход к нахождению периода колебаний - student2.ru ; (126)

- фазовый спектр

Энергетический подход к нахождению периода колебаний - student2.ru .(127)

Все эти спектральные функции для f(t) = A Энергетический подход к нахождению периода колебаний - student2.ru sin(ω1t) при β = 0,01 с-1, ω1 = 10 рад/с представлены на рис. 37 - 40.

Энергетический подход к нахождению периода колебаний - student2.ru

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

Энергетический подход к нахождению периода колебаний - student2.ru

Рис. 30. Дискретный амплитудный спектр «гребенка» (119)

Энергетический подход к нахождению периода колебаний - student2.ru

Рис.31. Фазовый спектр «гребенки»

Энергетический подход к нахождению периода колебаний - student2.ru

Рис. 32. Дискретный амплитудный спектр (по формуле (120))

Энергетический подход к нахождению периода колебаний - student2.ru

Рис. 33. Равномерный фазовый спектр

Энергетический подход к нахождению периода колебаний - student2.ru

Рис. 34. Четная функция, описывающая колебательный процесс

Энергетический подход к нахождению периода колебаний - student2.ru

Рис. 35. Нечетная функция, описывающая колебательный процесс

Энергетический подход к нахождению периода колебаний - student2.ru

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

Энергетический подход к нахождению периода колебаний - student2.ru

Рис. 37. Косинусный спектр затухающего колебания

Энергетический подход к нахождению периода колебаний - student2.ru

Рис.38. Синусный спектр затухающего колебания

Энергетический подход к нахождению периода колебаний - student2.ru

Рис.39. Фазовый спектр затухающего колебания

 
  Энергетический подход к нахождению периода колебаний - student2.ru

Рис. 40. Амплитудный спектр затухающего колебания

Энергетический подход к нахождению периода колебаний - student2.ru

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

При суммировании конечного (или бесконечного) числа гармоник в узком спектральном диапазоне получаются функции, представляющие пакеты. Их часто называют волновыми. Дело в том, что волна, приходя в данную точку пространства, возбуждает в ней колебания. Если в эту же точку в данный момент времени приходят другие волны, то колебания, создаваемые ими, складываются тем или иным образом. Ранее (72 - 73) мы показали, что колебания разных частот можно сложить как колебания одной частоты, но с разными фазами. Наиболее практически важным является случай суммирования гармоник с дискретным набором равноотстоящих частот сосредоточенных в определенном интервале (рис. 41).

 
2.1.20. Нелинейные колебания

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

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

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

В математическом виде имеем систему дифференциальных уравнений, описывающую указанные процессы:

Энергетический подход к нахождению периода колебаний - student2.ru = k1x – kxy; Энергетический подход к нахождению периода колебаний - student2.ru = k1xy – k2y (128)

Здесь коэффициенты k – постоянны, а х и y - функции времени t.

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

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

Энергетический подход к нахождению периода колебаний - student2.ru ,

и при этом (128) превращается в алгебраическую систему

k1x0 – kx0 y0 = 0 ,

k1x0 y0 – k2y0 = 0.

Отсюда находим

у0 = Энергетический подход к нахождению периода колебаний - student2.ru , х0 = Энергетический подход к нахождению периода колебаний - student2.ru .

Теперь ищем решение системы (128) в виде

х = х0 + α (t), у = у0 + β(t) . (129)

Подставляем (129) в (128), имеем:

Энергетический подход к нахождению периода колебаний - student2.ru = k1α – kx0β - k y0α – kαβ,

Энергетический подход к нахождению периода колебаний - student2.ru = k1x0β+ k1 y0α + k1αβ – k2β.

Членами kαβ и k1αβ можно пренебречь, т.к. они содержат произведение двух малых величин, х0 и у0 уже найдены. В результате получаем линеаризованную систему:

Энергетический подход к нахождению периода колебаний - student2.ru = - Энергетический подход к нахождению периода колебаний - student2.ru ; Энергетический подход к нахождению периода колебаний - student2.ru = Энергетический подход к нахождению периода колебаний - student2.ru . (130)

Продифференцируем первое уравнение по t и подставим Энергетический подход к нахождению периода колебаний - student2.ru из второго. Получим

Энергетический подход к нахождению периода колебаний - student2.ru = - k1k2α . (131)

Продифференцируем второе уравнение из (130) по t и подставим Энергетический подход к нахождению периода колебаний - student2.ru из первого. Получим

Энергетический подход к нахождению периода колебаний - student2.ru = - k1k2β . (132)

Уравнения (131) и (132) – стандартные уравнения колебаний. Частота колебания численности жертв и хищников одинакова и равна ω0 = Энергетический подход к нахождению периода колебаний - student2.ru . С учетом начальных условий получим, что существует сдвиг фаз между гармоническими функциями х(t) и у(t), описывающими численность жертв и численность хищников:

х(t) = х0 + α (t) = х0 + Acos (ω0t + φ1),

у(t) = у0 + β(t)= у0 + Bcos (ω0t + φ2).

Эти функции описывают гармонические колебательные процессы около смещенного положения равновесия.

Содержание

Введение
Требования к выполнению курсового проекта (работы)
  1.1 Организация курсового проектирования (технология выполнения). Структура и содержание курсового проекта
  1.2 Общие требования к пояснительной записке
  1.3 Требования, предъявляемые к графической части
Методические указания по выполнению курсового проекта (работы)
  2.1 Исходные положения. Теория и типовые задачи в пояснительной записке
  2.1.1. Гармонические колебания
  2.1.2. Математический маятник
  2.1.3. Пружинный маятник
  2.1.4. Комбинированные осцилляторы
  2.1.5. Колебания в электрических цепях
  2.1.6. Колебания в электростатическом поле
  2.1.7. Колебания в магнитном поле
  2.1.8. Вынужденные колебания. Резонанс
  2.1.9. Рекомендации по решению задач
  2.1.10. Сложение колебаний
  2.1.11. Затухающие колебания
  2.1.12. Релаксационные колебания
  2.1.13. Резонанс смещений, скоростей, ускорений
  2.1.14 Ангармонический осциллятор
  2.1.15 Модулированные колебания
  2.1.16 Энергетический подход
  2.1.17 Добротность
  2.1.18 Фазовые траектории
  2.1.19 Спектры колебаний
  2.1.20. Нелинейные колебания
  2.1.21. Параметрические колебания. Параметрический резонанс
  2.1.22. Автоколебания
  2.1.23. Режимы работы осциллятора при подводе энергии
  2.1.24. Колебания со многими степенями свободы
  2.1.25 Колебания связанных систем
  Вспомогательные материалы
  3.1. Перечень тем курсовых работ
  3.2. Указания к выполнению курсовых работ
  3.3. Пример выполнения курсовой работы
  3.4. Приборы и стенды
  3.5. Заключительные замечания
  3.6. Рекомендуемая литература
  Приложения
  Содержание

Наши рекомендации