Информатика -продвинутый курс

       

КОЛЕБАНИЯ МАТЕМАТИЧЕСКОГО МАЯТНИКА


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

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

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

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

Рис.7.17. Колебания математического маятника

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

(7.27)

На груз действуют две силы: сила тяжести

 и упругая сила натяжения стержня
. При выводе уравнения движения достаточно учесть лишь компоненту силы
  , направленную по касательной к дуге: F = mg sin ?, направлена она в сторону уменьшения ?. Сила
 перпендикулярна к касательной и вклада в это уравнение не дает.
Уравнение движения примет вид

(7.28)

Обычно в курсе физики ограничиваются исследованием малых колебаний. Если |?|<< 1, то уравнение (7.28) можно считать эквивалентным (так как sin ?



? ?; здесь и далее используется радианная мера углов) уравнению



Решение его элементарно:



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



то



или, как часто записывают,



где ? - так называемая, начальная фаза; А - амплитуда колебания; А и ?

легко выразить через начальные условия ?0

и v0.



Движение, происходящее по закону (7.29), называют гармоническим колебательным движением. Слово «гармонический» связывают с простой тригонометрической функцией (синусом или косинусом); так, гармоническим является и движение A sin (?t + ?), к которому также можно свести (7.29) (оно отличается лишь сдвигом фазы на ?/2).

Для изучения колебаний с большой амплитудой следует обратиться к уравнению (7.28), которое заведомо не интегрируется в элементарных функциях. Обезразмерим его, взяв за характерный масштаб времени период малого колебания. Если ? = t/T, то

(7.30)

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

Сведем (7.30) к системе двух уравнений первого порядка:

(7.31)

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

(7.32)

В безразмерных переменных x и ?

(7.33)

Как и при моделировании движения небесных тел, сохранение ? в ходе интегрирования - прекрасный критерий для изучения устойчивости метода, выбора шага и т.д. На рис. 7.18 представлен график зависимости ?(?) для ?0

= ?/2 и v0 = 0 (сплошная линия). На первый взгляд, это косинусоида (7.29), но, во-первых, это не так (зрительным впечатлениям в таких случаях доверять особо не следует), а, во-вторых, у этого движения период отнюдь не определяется формулой, следующей из решения задачи о малых колебаниях.


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



Рис. 7.18. Графики зависимости ?(?) для ?0 = ?/2 и v0 = 0 (сплошная линия) и гармонического движения с той же амплитудой ?/2 (пунктирная линия)

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

Вернемся снова к разговору о периодическом, но не гармоническом движении. Период колебаний в рассмотренном примере приблизительно равен 1,18 (определено в численном эксперименте). Уравнение гармонического движения с периодом Т и амплитудой A



(в нашем конкретном случае A = ?/2, T ? 1,18, ? =

0). В табл. 7.5 сведены результаты численного решения уравнений (7.31) (вторая строка) и табулирования функции при A

= ?/2, T ? 1,18, ? = 0 (третья строка) на промежутке времени, чуть большем периода. Хотя различия и невелики, но видно, что движение не является гармоническим.

Таблица 7.5

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

t

0,0

0,1

0,2

0,3

0.4

0,5

0,6

?реал

1,5708

1.3737

0,7971

-0,0437

-0,8688

-1,4104

-1,5689

?гарм

1,5708

1,3533

0.7611

-0,0418

-0.8332

-1,3938

-1.5686

?(t)

1,5710

1.3737

0.7938

-0,0473

-0,8696

-1.4077

-1.5631

t

0.7

0,8

0,9

1,0

1,1

1,2

1,3

?реал

-1,3331

-0,7228

0,1308

0,9374

1,4434

1,5632

1,2889

?гарм

-1,3090

-0,6870

0,1253

0,9028

1,4304

1,5619

1,2609

?(t)

-1,3299

-0,7216

0,1297

0,9371

1,4448

1,5631

1,2869

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


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

Гармоническое разложение функции f(t)

с периодом Т в общем случае имеет вид



причем число гармоник-слагаемых формально бесконечно велико. Те, кто изучал ряды Фурье, знают общие правила вычисления коэффициентов a0, а1, b1, а2, b2, ... Если ограничиться лишь небольшим число гармоник, скажем, тремя, то коэффициенты можно приближенно найти интерполяцией. Взяв за узлы точки t = 0; 0,4 и 1 (выбор достаточно произволен) и решив систему трех линейных алгебраических уравнений,получим



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



Рис. 7.19. Периодический сигнал, подобный вырабатываемому генератором

строчной развертки в телевизоре

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



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



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



Преобразуем его к виду

(7.34)

где, как и выше,
, а
 (коэффициент 2 записан по традиции для К / 2т удобства). При малых колебаниях уравнение (7.34) превращается в

(7.35)

Его решение таково: затухающие колебания при к < ? и затухание без колебаний при к ? ?. Все это можно проверить в ходе численного моделирования, хотя уравнение (7.35) допускает аналитическое решение. Приведем его: при к < ?



где
, А -

амплитуда, ? - начальная фаза (А и ? легко выразить через начальные значения ?0 и v0). При k ? ?



где А и В также можно выразить через ?0 и v0.

Что же касается уравнения (7.34), то его аналитическое решение отсутствует, и при численном моделировании можно поставить ряд задач о том, насколько решения уравнений (7.34) и (7.35) различаются в зависимости от начальной амплитуды.

Вынужденные колебания. Если на маятник воздействует внешняя сила F(t),

меняющаяся со временем, то уравнения движения получаются из (7.34) добавлением F(t) к правой части. Рассмотрим лишь случай периодического внешнего воздействия: F(t) = F0 cos ?t, где ?

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

(7.36)

где
. При малой амплитуде результирующего движения уравнение (7.36) примет вид

(7.37)

Движение, описываемое уравнением (7.37), состоит из двух этапов. На первом оно складывается из двух колебательных движений: затухающих собственных колебаний с частотой
 (при к < ?) и вынужденных колебаний с частотой ?. На втором этапе, по истечении времени t >> 1/k, остаются лишь вынужденные периодические колебания, амплитуда которых зависит от соотношения частот ? и ?1 и резко возрастает при ? ? ?1 - явление резонанса, описанное в любом учебнике физики.


Численное интегрирование уравнения (7.37) необязательно, так как решение можно записать в виде формул, содержащих лишь элементарные функции:

(7.38)

А

и В — произвольные постоянные, находятся из начальных условий.

Исследования переходного процесса установления стационарных вынужденных колебаний, резонанса, биений, возникающих при k = 0 и ? ? ?1  (рис. 7.20-7.22), могут быть, конечно, проведены с использованием формул (7.38) простым табулированием с выводом результатов на экран компьютера в форме, удобной для восприятия; они же могут быть и объектами численного моделирования.



Рис. 7.20. Установление стационарных вынужденных колебаний маятника

при наличии трения при к = 0,5; ? = ?/2, ? = ?, f = 2?.



Рис. 7.21. Биения в системе с близкими частотами собственных колебаний

и с вынуждающей силой при k = 0; ? = 889?/9000, ? = ?/9, f = ?/70

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



Рис. 7.22. Возрастание амплитуды колебаний при прохождении через резонанс при k =

0; ? = 889?/9000, ? = ?/9, f = ?/70

Параметрические колебания.

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


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

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

(7.39)

где ?(t) — заданная функция, определяющая закон изменения частоты. Мы ограничимся простейшим случаем гармонического изменения ?2(t):



где ?

- частота изменения величины ?2(t).

При малых амплитудах колебаний и отсутствии трения уравнение (7.39) превращается в



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

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



Рис. 7.23 Фазовая диаграмма с зонами параметрического резонанса

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

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

Как можно численно установить границу зоны параметрического резонанса, например, первой? - Для этого надо задаться некоторыми значениями ?

(например, 0,1) и ?

(например, 0,3), не принадлежащими зоне неустойчивости, и проинтегрировать численно уравнение (7.40). Удобно предварительно обезразмерить время переменной ? = ?0t, после чего уравнение примет вид



(7.41)

Здесь
. Затем, медленно увеличивая ? (например, с шагом 0,01) и не меняя ?,

интегрировать уравнение (7.41), пока не попадешь в зону неустойчивости, и далее, пока не выйдешь из нее. Затем следует увеличить ? (например, взяв ? = 0,2) и снова повторить процедуру прохождения по значениям ? и т.д. - постепенно вырисуется картина границы зоны параметрического резонанса на фазовой плоскости.

Нарастание колебаний при параметрическом резонансе, описываемом уравнением (7.40), является неограниченным. Физически такого быть не может. Ограничение амплитуды колебаний наступает либо за счет учета трения, либо при возврате к sin? в уравнении (7.39), либо за счет обоих факторов. Следует учесть, что наличие трения не только ограничивает размах параметрических колебаний, но и «приподнимает» зоны параметрического резонанса над осью ?

на фазовой плоскости ?, ?, причем в разной мере. Моделирование этого и других явлений при параметрическом резонансе - интересная исследовательская работа.

Многогранность задачи об одномерных колебаниях. Колебания математического маятника одномерны в том смысле, что они описываются одной функцией ?(t) (хотя они и происходят в двумерном пространстве - плоскости, но жесткий стержень ликвидирует одну из степеней свободы, и в обычных декартовых координатах x(t), y(t) выражаются друг через друга).

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

(7.42)

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

• математический маятник:

• пружинный маятник, где сила, действующая на тело. определяется законом Гука;

• «физический» маятник-тело, свободно вращающееся около горизонтальной оси;

• крутильный маятник наручных часов - симметричное тело, совершающее колебания около вертикальной оси под действием спиральной пружины;

• ток в колебательном контуре;

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

• колебания магнитной стрелки компаса.

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

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


Содержание раздела