Основные положения теории колебаний осцилляторов

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

а. Аналитически гармонические колебания выражаются уравнением колебаний с использованием тригонометрических функций синуса или косинуса. Например:

x(t) = A cos(wt + j).

В соответствии с формулой Эйлера = cosj + i sinj (i = ) гармонические колебания аналитически можно представить и в комплексной форме. Под гармонической функцией в комплексном пространстве следует понимать экспоненту с мнимым показателем: x(t) = A

б. Гармонические колебания, например, напряжения в колебательном контуре U(t) = U₀ cos(wt + j), можно представить графически временной диаграммой, т.е. графиком U(t), где по горизонтальной оси отложено время t, по вертикальной оси колеблющаяся динамическая переменная U (рис. 1-2). Эта диаграмма может быть получена с помощью осциллографа, подключенного, например, к активному сопротивлению колебательного контура.

в.Если фазовые характеристики (фазовые соотношения) не интересны, то колебание удобно представлять графически спектральной диаграммой. На графике по горизонтальной оси откладывается частота w (или n = ), а по вертикали амплитуда динамической переменной (или энергия, пропорциональная квадрату амплитуды и частоты). В случае гармонического колебания x(t) = A cos(wt + j) спектр состоит из одной частоты w (рис. 1-3а). В реальной колебательной системе в ее спектре содержатся набор частот в определенном частотном диапазоне (рис.1-3б).

г. Гармонические колебания можно представить векторной диаграммой. Гармоническое колебание x(t) = A cos(wt + j) определено, если заданы амплитуда А, начальная фаза j и циклическая частота w. Эти характеристики на векторной диаграмме выражаются следующим образом. Выбирается полюс 0 и луч ОМ (рис.1-4). Относительно луча ОМ отсчитывается начальная фаза. Из полюса проводится вектор А, модуль которого ра вен амплитуде колебания, причем вектор Аповорачивается на угол j, равный начальной фазе колебания (положительный угол соответствует повороту против часовой стрелки, отрицательный – против). Мысленно считается, что вектор Авращается относительно полюса 0 с угловой скоростью w. Таким образом, векторная диаграмма отображает все характеристики уравнения колебаний.

д.Представление колебательного процесса с помощью фазовой траектории (диаграммы состояния). Этот способ рассмотрим на примере механической колебательной системы. Известно, что состояние и изменение состояния любой классической механической системы материальных точек с их известными массами полностью определено, если в данный момент времени одновременно заданы координаты и скорости всех материальных точек механической системы (в теоретическом курсе классической механики состояние системы задается координатами и импульсом). Пусть смещение материальной точки, прикрепленной к пружине, совершается по закону x(t) = A cos(wt + j), тогда скорость материальной точки изменяется по закону v_x = . Возведя в квадрат эти уравнения и принимая во внимание тождество cos² (wt + j) + sin² (wt + j) = 1, получим:

+ = 1 или x² + = A².

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

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

Фазовая траектория затухающих колебаний имеет вид спирали, скручивающейся к началу 0 (рис.1-6). Действительно, уравнение для смещения затухающих колебаний имеет вид:

x(t) = A cos wt

(для простоты за начальную фазу принято значение j = 0).

Уравнение скорости v_x = = - A (b cos wt + w sin wt). При малых затуханиях (b<<w) v_x = -Aw sin wt. Возведя в квадрат уравнения для смещения и скорости, получим:

x² + = A² .

Это и есть уравнение спирали, скручивающейся к началу 0.

1.1.3. Динамика гармонического осциллятора

Пример.

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

Решение

На ареометр массой действует сила тяжести и сила Архимеда. Если слегка притопить ареометр на глубину x, то сила Архимеда возрастет по модулю на величину F = rSxg . Понятно, что при всплытии ареометра выше положения равновесия сила тяжести будет на ту же величину F = rSxg больше силы Архимеда. Сила Архимеда направлена против направления смещения ареометра. Дополнительная сила Архимеда F_x = -rSxg играет роль возвращающей силы (напомним, в пружинном маятнике эту роль играет сила упругости). Второй закон Ньютона в проекции на ось 0Х имеет вид m = - rSxg. После соответствующих преобразований дифференциальное уравнение примет вид: или ,

где - квадрат собственной частоты колебаний ареометра в жидкости. Решением уравнения является уравнение гармонических колебаний x = X_max cos (w₀t + j), X_max – амплитуда колебаний ареометра.

Динамические уравнения колебаний пружинного маятника и ареометра и их решения имеют одинаковый вид. Роль возвращающей силы в пружинном маятнике играет сила упругости F_x = - kx. Для ареометра возвращающей силой является сила Архимеда, которая формально имеет тот же вид, что и сила упругости: F_x = -rSgx = - bx, где b = rSg. В этой связи возвращающая сила вида F = - bx получила название квазиупругойсилы, т.е. силы, формально схожей с упругой силой (приставка «квази» – от латинского quasi – якобы, как будто, вроде).

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

Рис.1-9

Математический маятник (рис. 1-8). Дифференциальное уравнение, описывающее динамику колебаний математического маятника, имеет вид:

J = - mgla,

где J = ml ² – момент инерции маятника относительно оси 0. Собственная частота маятника , период T = = . Решение уравнения:

a = a_max cos (w₀t + j).

Колебания материальной точки m, прикрепленной к двум струнам (рис. 1-9). Дифференциальное уравнение имеет вид: m = - . Собственная частота системы , период T₀ = = .

Решение уравнения: x = X_max cos (w₀t + j).

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

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

F = = - kx.

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

1.1.4. Динамика затухающих колебаний

Пусть пружинный маятник находится в вязкой среде (рис. 1-10). На материальную точку в этом случае действует сила сопротивления среды F_c = - mv_x = - m , где m - коэффициент сопротивления вязкой среды. 2-ой закон Ньютона в проекции на ось 0Х примет вид: , или . Обозначим и учтем, что - квадрат собственной частоты колебательной системы. Тогда 2-ой закон Ньютона примет вид:

. (5)

Уравнение (5) - линейное однородное дифференциальное уравнение второго порядка с постоянными коэффициентами.

 

 

Рис. 1-10

 

Решение дифференциального уравнения (5) ищется, как известно, с помощью подстановки экспоненциальной функции x(t) = , где l = const, е - основание натурального логарифма. Напомним, функцию записывают также в форме exp (lt), т.е. º exp (lt). Дважды продифференцируем x(t) = и результаты подставим в (5): l² + 2 bl + w₀² = 0.

Полученное алгебраическое уравнение, соответствующее дифференциальному уравнению (5), называется характеристическим уравнением. Корнями этого уравнения являются: l_1,2 = - b ± . При малых затуханиях (при w₀ >> b) подкоренное выражение будет отрицательным, следовательно, корень является мнимой величиной. Представим подкоренное выражение в виде (iw)², где: i = - мнимая единица, w = - вещественная величина. Таким образом, решениями характеристического уравнения являются корни:

l₁ = - b + iw ; l₂ = - b - iw ,

Условие w₀ >> b означает, что сопротивление среды не слишком большое и колебательная система до полного затухания успевает многократно пройти через положение устойчивого равновесия. Корни характеристического уравнения l₁ и l₂ - это те значения l, при которых подстановка x(t) = удовлетворяет дифференциальному уравнению . Решение дифференциального уравнения (5) при данных корнях характеристического уравнения l₁ и l₂ принимает вид:

x(t)= (C₁ + C₂ ).

Функция x(t) описывает колебание, следовательно, эта функция должна быть вещественной. Условие вещественности означает равенство комплексно сопряженных величин:

C₁ + C₂ = C₁^* + C₂^* .

Следовательно, C₁ = C₂^* и C₂ = C₁^*. В показательной форме постоянные комплексные коэффициенты C₁ и C₂ можно представить, например, в виде: C₁ = ; C₂ = . Подставим эти коэффициенты в решение (6):

x(t)= ( + ). (6)

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

= cosa + i sina ; = cosa - i sina, то cosa = ( + ) и sina = ( - ), следовательно: x(t)= ( + ) = A cos (wt + a). (7)

График затухающих колебаний (7) приведен на рис. 1-11.

 

 

Рис. 1-11.

Частота затухающих колебаний определяется соотношением w = (частота затухающих колебаний w не равна собственной частоте w₀ колебательной системы), период - T = , где коэффициент b носит название коэффициента затухания. Приведенное решение верно при условии w₀ > b. При w₀>> b частота и период затухающих колебаний колебательной системы почти совпадают с собственной частотой и периодом этой колебательной системы: w » w₀ и T » .

Проведем небольшой анализ решения (7).

Роль амплитуды в (7) играет величина A , где А - амплитуда в момент времени t = 0. Амплитуда уменьшается по экспоненциальному закону (знак минус), следовательно, затухающие колебания не являются строго периодическими. Однако определенная периодичность проявляется, например, в повторяемости прохождения материальной точкой положения равновесия.

Из выражения для переменной амплитуды затухающих колебаний A определяется феноменологический смысл коэффициента затухания b. Показатель степени является всегда безразмерной величиной, следовательно, b в измеряется в обратных секундах. Таким образом, коэффициент затухания может быть представлен в виде b = , где t - некоторое характерное время колебательной системы. За промежуток времени t = t амплитуда уменьшается в e » 2,7 раз (A = A = ). Итак, коэффициент затухания b равен обратной величине промежутка времени t, за которое амплитуда уменьшается в e » 2,7 раз. Заметим: уменьшение амплитуды в 2,7 раза (почти в 3 раза) на практике воспринимается как почти полное затухание колебаний. Время t часто называют временем релаксациизатухающих колебаний.

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

l = ln = ln = bT.

Само отношение амплитуд = называется декрементом затухания.

За время релаксации t колебательная система совершает N_e = колебаний, где T – период колебания (время одного колебания). Из определения N_e следует, что N_e = = или l N_e = 1. Число N_e разных колебательных систем могут существенно различаться. Например, в обычных механических колебательных системах N_e может достигать нескольких десятков, в электрических контурах - нескольких сотен, у кварцевого резонатора может достигать значений N_e~10⁶. Величина N_e определяет добротность колебательной системы Q, которая выражается формулой Q = = = p N_e, или Q » .Чем выше N_e, тем добротнее колебательная система.

 

1.1.5. Динамика вынужденных колебаний. Импеданс колебательной

Системы

Вследствие потери энергии свободные колебания осциллятора являются затухающими. Для генерации незатухающих колебаний необходимо эти потери компенсировать. Например, затухание колебания качели можно компенсировать, если в нужный момент времени подталкивать качели. В этом примере внешняя вынуждающая сила является импульсной силой. Рассмотрим вынужденные колебания пружинного маятника под действием непрерывной гармонической вынуждающей силы F_x(t) = F₀ coswt, где w - частота внешней вынуждающей силы (рис. 1-12).

На колеблющееся тело маятника действуют сила сопротивления среды F_c = -m v_x = - m (примером является сила Стокса F = 6phRv) и сила упругости пружины F_уп. = - kx. Второй закон Ньютона примет вид: , или в компактной записи

, (6)

где ; . Уравнение (6) является линейным неоднородным дифференциальным уравнением второго порядка с постоянными коэффициентами 1, 2b, и .

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

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

Имеется некоторое неоднородное дифференциальное уравнение: . Здесь С, D - вещественные постоянные коэффициенты; f(t) и y(t) - вещественные функции; - мнимая единица. Функция z(t) является комплексной: z(t) = x(t) + iy(t), где x(t) и y(t) - вещественные функции. Комплексной функцией является и правая часть уравнения . После подстановки в это уравнение функцию z(t) = x (t) + iy (t)получим:

+ i( ) = .

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

= f (t) и = .

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

1) к правой части решаемого уравнения = f(t) прибавляется мнимая функция . При этом получается уравнение в комплексной форме . Решением уравнения является комплексная функция z(t) = x (t) + iy (t);

2) Из этого решения выделяется действительная часть x(t), которая и является решением уравнения = f(t). Выбор функции диктуется видом функции f(t) в соответствии с формулой Эйлера: если f(t) - функция косинуса, то в качестве принимается функция синуса.

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

(где и ) или по формуле Эйлера

. (7)

Искомой комплексной функцией является z(t) = x (t) + iy (t). В качестве решения выберем функцию z(t) = А , где А - комплексная амплитуда. Определим А, для чего подставим z(t) в уравнение (7):

[ - w²A + i2bwA + A] = .

Т.к. данное уравнение выполняется при всех t, то А =