Статьи Соросовского Образовательного журнала в текстовом формате
С помощью маятников и их простых математических моделей демонстрируются динамический хаос, гистерезис, мультистабильность и другие феномены, присущие нелинейным колебательным системам.
НЕЛИНЕЙНЫЕ МАЯТНИКИ И ИХ МОДЕЛИБ. П. БЕЗРУЧКО
Саратовский государственный университет им. Н.Г. Чернышевского
ВВЕДЕНИЕ
Познавая природу, люди опираются на представления об объектах и явлениях, которые считаются понятными и приобрели статус эталонных. Для колебательных систем, движения в которых обладают той или иной степенью повторяемости во времени, такими исключительно удобными эталонами являются маятники. Их рассматривали основатели физики и создатели современной теории колебаний и волн. К ним обращаются и сейчас при изучении колебаний в курсах общей физики, но, как правило, ограничиваются узким классом гармонических и затухающих квазигармонических движений. Заглянем с помощью маятников в безграничный мир сложных нелинейных колебаний, не выходя при этом за уровень, доступный вооруженному компьютером школьнику.
ЧТО ПОНИМАЮТ ПОД СЛОВОМ "МАЯТНИК"
Обратимся к "Физическому энциклопедическому словарю": "Маятник - твердое тело, совершающее под действием приложенных сил колебания около неподвижной точки или вокруг оси. Обычно под маятником понимают тело, совершающее колебания под действием силы тяжести, при этом ось маятника не должна проходить через центр тяжести тела. Простейший маятник состоит из небольшого массивного груза, подвешенного на нити (или легком стержне). Если считать нить нерастяжимой и пренебречь размерами груза по сравнению с длиной нити, а массой нити по сравнению с массой груза, то груз на нити можно рассматривать как материальную точку, находящуюся на неизменном расстоянии от точки подвеса (рис. 1, а). Такой маятник называется математическим_" Приведенное в энциклопедии определение относится к объектам механики, а в настоящее время термин "маятник" широко используют в приложении к системам различной природы (см., например, [1]). Так, электрическим маятником называют цепь, состоящую из конденсатора и катушки индуктивности (колебательный контур, рис. 1, б ), химическим - смесь химикатов, вступающих в колебательную реакцию, экологическим маятником - две взаимодействующие популяции хищников и жертв. Этот же термин применяется к экономическим системам, в которых имеют место колебательные процессы.
Далее мы будем называть маятником систему, способную совершать колебания около положения равновесия. При этом ограничимся строго диссипативными системами, для описания состояния которых при отсутствии внешних воздействий достаточно указать две переменные величины: x, называемую координатой, и скорость ее изменения . Маятник, выведенный из состояния равновесия и предоставленный сам себе, совершает колебания, называемые собственными. Движения маятника при наличии внешнего воздействия именуют вынужденными.
СОБСТВЕННЫЕ КОЛЕБАНИЯ
Из куска пластилина, карандаша и иголки создадим конструкцию, напоминающую маятник, упомянутый в энциклопедии, и поэкспериментируем с ней. Что бы мы ни проделывали, но предоставленный самому себе наш маятник в конце концов оказывается в состоянии устойчивого равновесия, в котором стержень направлен вертикально, а груз располагается ниже оси (рис. 1, а). В зависимости от начальных положения и скорости, формы, массы, трения в оси этому предшествуют различные варианты движений.
Будем характеризовать состояние маятника величиной угла отклонения от вертикали x и рассматривать зависимости x и угловой скорости от времени t. Качественный вид графиков x(t) и (временных реализаций) приведен на рис. 1, а. На рисунке цветовым фоном выделены две области, в которых движения качественно различны. Синяя область соответствует большим отклонениям от устойчивого равновесия, груз может достигать верхнего положения или даже проскакивать его вращаясь. Временные реализации координаты и скорости здесь качественно отличаются друг от друга, и (обратите внимание) интервал времени между соседними максимумами кривых меняется с уменьшением размаха колебаний при их затухании. Эта особенность получила название "неизохронность". Отмеченные закономерности просматриваются еще более ярко в осциллограммах колебаний электрического маятника (рис. 1, б ), для него x(t) и имеют смысл напряжения и тока в цепи. Похоже ведет себя и груз на пружине (рис. 1, в).
Области, выделенные на рис. 1 желтым цветом, соответствуют малым колебаниям. Вид временных реализаций здесь близок к затухающей по экспоненте косинусоиде или синусоиде:
x(t) = x0 exp (- dt) cos (w1t + j),
где постоянная d определяет затухание, w1 = 2p / T1 - частота, T1 - период (вернее, квазипериод) малых колебаний, j - начальная фаза, которая определяет сдвиг кривой по оси t. Временные реализации координаты и скорости по форме подобны, но сдвинуты во времени на T1 /4 (по фазе на p /2). При малых амплитудах колебаний временные интервалы между соседними максимумами не меняются, то есть маятник обладает свойством изохронности, которое отсутствовало в случае больших отклонений. На эту особенность малых колебаний обратил внимание еще Галилей, наблюдая за люстрой в соборе.
ЛИНЕЙНЫЕ КОЛЕБАНИЯ
Рассмотренные малые изохронные колебания часто называют линейными. Название связано с математической природой уравнений, используемых для описания такого процесса. Например, монотонное затухание колебаний (как и безграничное возрастание) наглядно моделируется линейным отображением, которое можно получить из (1). Для этого следует перейти к рассмотрению колебательного процесса (1) в фиксированные моменты t = nT1 , где n = 0, 1, 2, _ - дискретное время. Обозначим значения x(t) в моменты nT1 и (n + 1)T1 как xn и xn + 1 и введем параметр a = exp (- dT1) < 1. Выразив с помощью (1) значение xn + 1 через xn , получим одномерное отображение
xn + 1 = axn .
График отображения представляет собой прямую линию, проходящую через начало координат. От графика линейного отображения, приведенного на рис. 3, б синей линией, он отличается лишь сдвигом. Если задать начальное значение x0 и по полученной формуле подсчитать x1 , затем подставить x1 , вычислить x2 и т.д., можно убедиться, что значение xn монотонно стремится к нулю. Не составляет труда проделать то же самое графически, проведя построения, аналогичные показанным стрелками на рис. 3, б, и убедиться, что изображающая точка приближается к нулю. Такое построение называют диаграммой Ламерея (см. [4]). При замене знака перед d (при этом a > 1) линия графика отображения расположится выше диагонали и изображающая точка убежит на бесконечность.
Кроме того, основанием назвать колебания в желтой области линейными является и то, что решение вида (1) имеет обыкновенное дифференциальное уравнение:
известное как уравнение линейного диссипативного осциллятора. Переменная x и ее производные входят в это уравнение в первой степени (линейно), поэтому уравнение удовлетворяет принципу суперпозиции. Это означает, что сумма двух любых решений уравнения или решение, умноженное на константу, тоже являются решениями. В виде (3) удается записать уравнения динамики различных колебательных систем. Широкий набор интересных примеров имеется в [1]. Проинтегрировав уравнение (3), можно получить (1) с
Поведение большинства маятников хорошо описывается линейными математическими моделями, если движение ограничено некоторыми пределами. Например, для груза на стержне значения х должны быть столь малыми, что sin x © x, а значения ограничены сверху условием пропорциональной зависимости сил вязкого трения от скорости.
НЕЛИНЕЙНЫЕ КОЛЕБАНИЯ
Для моделирования движений наших маятников в синей части рис. 1 необходимы нелинейные отображения и дифференциальные уравнения. Уравнение осциллятора (3) станет нелинейным, если в него войдут члены, содержащие любые нелинейные функции от x и его производных (квадраты, кубы и т.п.). Так, колебания математического маятника (рис. 1, а) хорошо моделируются уравнением (3), если x в последнем слагаемом заменить на sinx. В получившемся нелинейном уравнении в отличие от (3) суперпозиция решений уже не будет его решением. Нелинейности отображения xn + 1 = f (xn) соответствует отличный от линейного вид функции f.
По сравнению с единым стандартом (1) линейных колебаний многообразие типов нелинейного поведения очень велико. Так, в трех примерах на рис. 1 при полной идентичности линейной стадии характер нелинейных участков существенно различается. Это связано с особенностями каждого из маятников: характером их нелинейности, видом зависимости параметров системы от координат и скоростей. Например, в случае математического маятника (см. рис. 1, а) нелинейность определяется синусоидальной зависимостью от угла поворота величины момента силы тяжести относительно оси вращения. В электрическом маятнике с полупроводниковым конденсатором (варакторным диодом, см. рис. 1, б ) она вызвана зависимостью емкости и сопротивления диода от напряжения. В пружинном маятнике нелинейность обусловлена видом зависимости силы вязкого трения от скорости движения и силы упругости от величины деформации пружины. Например, при сильном сжатии витки пружины смыкаются, так что сила упругости резко возрастает по сравнению с предписываемой законом Гука (пружина становится более жесткой). Обратите внимание, что в этом случае период колебаний уменьшается при увеличении их размаха (см. рис. 1, в). По аналогии с грузом на пружине нелинейность маятников, приводящую к уменьшению периода колебаний с ростом их амплитуды, стали называть нелинейностью типа жесткой пружины, а в противоположном случае - типа мягкой пружины.
ВЫНУЖДЕННЫЕ КОЛЕБАНИЯ
Будем периодически подталкивать маятник, не давая ему вернуться в состояние устойчивого равновесия. Результат в значительной степени зависит от способа воздействия. Мы же выберем очень простой вариант - экспериментируя с грузом на стержне, будем периодически брать груз руками и поворачивать маятник на некоторый одинаковый угол А, а затем отпускать каждый раз с нулевой начальной скоростью. В колебательном контуре с диодом аналогичный способ воздействия реализуется с помощью источника импульсной ЭДС, обеспечивающего прямой ток через диод. При этом во время действия импульса собственные колебания быстро затухают, а по его окончании начинаются с близкой фазы (рис. 2, а, б ).
Если короткие импульсы воздействия следуют через равные интервалы времени Т (период воздействия), временная реализация вынужденных движений состоит из цугов затухающих собственных колебаний (последовательностей уменьшающихся отклонений от положения равновесия), как это показано на примере контура с диодом на рис. 2. В зависимости от амплитуды и периода воздействия, уровня затухания, вида нелинейности, а также начальных условий в системе могут устанавливаться различные повторяющиеся (периодические) движения - циклы. Времена повторения (периоды циклов) равны периоду воздействия или кратны ему - составляют величину kТ, где k принимает значения 1, 2, 3, _
Вариантов вынужденных периодических движений может быть очень много, поэтому их требуется как-то систематизировать. Например, для маятников с нелинейностью типа мягкой пружины циклы по признаку подобия можно условно разделить на две группы. В каждой из них сохраняются некоторые особенности формы временных реализаций, а также предельных циклов в фазовом пространстве. Первую группу составляют циклы периода воздействия 1Т, существующие в области низких частот f = 1/ T # f1 = 1/ T1 , - субгармонические. Так как период воздействия велик по сравнению с временным масштабом собственных движений, на временных реализациях этих циклов в общем случае имеется несколько максимумов (см. рис. 2, б ). Вторую группу составляют циклы периода kT, где k = 2, 3, _, которые наблюдаются при больших частотах воздействия f © (0,5-2) f1 . Временная реализация такого цикла (см. рис. 2, а) как бы собрана из k начальных кусков нелинейной стадии осциллограммы рис. 1, б. Так как при увеличении амплитуды воздействия смена этих режимов сопровождается последовательным увеличением k на единицу, их называют циклами последовательности добавления периода. Достаточно традиционно условное обозначение циклов - Gm / l . Здесь l соответствует отношению периода воздействия к квазипериоду собственных колебаний (его можно оценить по числу максимумов на интервале Т осциллограммы), а m - период цикла, измеренный в единицах T.
БИФУРКАЦИИ
При изменении параметров маятника и внешнего воздействия все циклы эволюционируют. При малых амплитудах колебаний они плавно переходят друг в друга, а при достаточно больших амплитудах, когда проявляется нелинейность, испытывают качественные перестройки - бифуркации [4]. Наиболее типична бифуркация удвоения периода (period doubling), в результате которой колебания начинают повторяться через раз, немного изменив свою форму. Для нелинейных маятников при изменении параметров характерны последовательности (каскады) бифуркаций удвоения периода: циклы последовательно увеличивают свой период до бесконечности. В итоге колебания становятся неповторяющимися во времени, беспорядочными, хаотическими [3]. При дальнейшем изменении параметров хаотические колебания, в свою очередь, видоизменяются, происходят переходы к новым циклам и новым видам хаотических движений и т.д.
Множество вариантов вынужденных колебаний, бифуркации, хаос и многие другие явления, свойственные маятникам, будем исследовать с помощью их математических моделей. Наиболее полную информацию дают дифференциальные уравнения. Мы же ограничимся рассмотрением лишь субгармонических колебаний, а для этого достаточно более простых и наглядных модельных отображений.
ДИСКРЕТНАЯ МОДЕЛЬ МАЯТНИКА
С МЯГКОЙ ПРУЖИНОЙ
Воспользуемся временной реализацией субгармонических колебаний маятника с мягкой пружиной (см. рис. 2, б ) и промоделируем ее последовательностью цугов собственных колебаний (рис. 2, в). Установим, что в цуге зависимость x(t) имеет вид (1), но квазипериод линейно зависит от начальной амплитуды цуга xn : теперь это не T1 , а T1(1 + bxn), где b = const > 0. Тогда в цуге, начинающемся в момент времени nT, зависимость переменной от времени можно представить как x(t) = = xn exp (- dt) cos [2pt / T1(1 + bxn)], если t отсчитывать от начала цуга. С учетом оговоренных особенностей способа внесения воздействия начальная амплитуда следующего цуга (см. рис. 2, в)
где А - амплитуда воздействия, N = T1 / T - нормированная частота воздействия, d = dT1 - коэффициент линейной диссипации, b - коэффициент нелинейности. График отображения (4) имеет вид изменяющейся по амплитуде и частоте синусоиды (рис. 3, а), ограниченной сверху и снизу прямыми, наклоненными под углами ? arctg (exp(- d / N )) к горизонтальной оси.
ВИДЫ КОЛЕБАНИЙ. МУЛЬТИСТАБИЛЬНОСТЬ
Изменение параметров A, N, d, b приводит к перемещению графика отображения (4) по вертикали и его деформации, к изменению числа, места и характера его пересечения с диагональю xn = xn + 1 . Например, при b = = 0 график превращается в прямую (рис. 3, б ) и имеет одно пересечение в точке С. Изображающая точка попадает в С из любого начального значения x0 (см. рис. 3, б ), то есть эта точка является единственным аттрактором [3, 4]. Она соответствует колебаниям маятника, повторяющимся через один период внешнего воздействия (циклу периода 1Т ). В общем случае точек пересечений может быть несколько. Например, на рис. 3, а их три, из которых одна (правая) является аттрактором, а две другие неустойчивы. Аттракторы, соответствующие устойчивым колебаниям, можно найти последовательно вычисляя по (4) значение xn + 1 . Начальное значение xn = x0 выбирается произвольно. Это можно сделать и графически с помощью диаграммы Ламерея (см. рис. 3, а). Через достаточно большое число шагов (например, n > 100 или 500) устанавливается колебательный режим с определенным повторяющимся набором или постоянными усредненными значениями переменных. Повторение значений xn через раз означает установление цикла периода 2. В циклах периода М значения xn повторяются через М шагов дискретного времени. Например, при значениях параметров, соответствующих рис. 3, а, в зависимости от x0 в модели устанавливаются колебания периода 3 или 1. Такая ситуация называется бистабильностью, а если возможных видов колебаний несколько, говорят о мультистабильности. Выбор одного из мультистабильных видов колебаний определяется заданием начальных условий. Неустойчивые циклы не реализуются в эксперименте, но оказывают существенное влияние на динамику маятника.
КАРТЫ РЕЖИМОВ. ПОДОБИЕ
Для описания множества типов колебаний маятника используем подходы, разработанные в нелинейной динамике. Наглядное представление о колебаниях, которые устанавливаются в системе после затухания переходного процесса, дают: 1) бифуркационные диаграммы - зависимости значений xn , получающихся после достаточно большого числа шагов счета, от величин параметров (рис. 3, в-д ); 2) карты режимов - на плоскостях параметров цветом или с помощью линий выделяются области значений, при которых устанавливается тот или иной колебательный режим (рис. 4).
Рассмотрим карту (см. рис. 4), построенную на плоскости параметров воздействия. Бросаются в глаза цветные языки, структуры которых качественно подобны. Это области сложных колебаний на базе различных субгармонических циклов. В крайней правой области эволюционирует цикл G: при этих значениях параметров существуют колебания периода 2, 4, 8, _ и хаос на его базе. Левее располагаются области режимов на базе циклов G1/2 , G1/3 и т.д. Эти области имеют качественно одинаковую структуру - они подобны друг другу и самоподобны. Самоподобие означает устройство по типу матрешки: основной конструктивный элемент, представленный на черно-белых фрагментах рисунка, дублируется во все более мелких масштабах; он просматривается во всех цветных языках и их фрагментах. Но в отличие от матрешки, границы составляющих элементов которой не пересекаются, области существования различных видов колебаний на плоскости параметров маятника при достаточно малых уровнях затухания накладываются друг на друга (см. рис. 4). При значениях параметров, соответствующих областям пересечения цветных языков, в системе могут установиться несколько различных видов колебаний - имеет место мультистабильность.
НЕЛИНЕЙНЫЙ РЕЗОНАНС. ГИСТЕРЕЗИС
Результат раскачивания маятника внешней силой зависит от временных масштабов воздействия. Это хорошо иллюстрируют зависимости установившихся значений xn от нормированной частоты N импульсов воздействия при неизменной величине А (см. рис. 3, в-д ). Оказывается, при некоторых частотах следования импульсов подталкивание маятника наиболее эффективно - наблюдается увеличение амплитуды колебаний (резонанс). Этот феномен более заметен, когда амплитуды колебаний незначительны | x | ! 1 или b = 0. В этом (линейном) случае существуют лишь регулярные колебания периода воздействия, установившееся значение xn единственно, и график xn(N ) имеет вид линии - резонансной кривой (см. рис. 3, г). Максимумы кривой соответствуют различным резонансам: основной, N = 1, наблюдается при равенстве периодов воздействия и собственных колебаний. Максимумы на частотах N < 1 (субрезонансы) соответствуют ситуациям, когда внешние толчки следуют точно через 2, 3, 4, _ периода собственных колебаний. При наличии нелинейности (b ? ? 0) явление резонанса сохраняется, но из-за бифуркаций и хаоса линии графиков xn(N ) раздваиваются или размываются (см. рис. 3, в). Качественное различие резонанса при малых и больших амплитудах колебаний (линейного и нелинейного) удобно проследить наблюдая за одним из максимумов резонансной кривой при увеличении А (см. рис. 3, д ). С ростом амплитуды кривая изгибается, резонансная частота уменьшается, а затем линия разрывается на ветви. Ветви соответствуют разным нелинейным видам колебаний. Например, на рис. 3, д при А = 2 верхняя ветвь принадлежит циклу G1 , а нижняя - G1/2 .
На частотах воздействия, соответствующих закрашенной области резонансной кривой на рис. 3, в, д и соответственно в выделенной цветом области черно-белого фрагмента на рис. 4, имеет место бистабильность. Изменение параметров при переходе через эту область сопровождается гистерезисом. Греческое слово "гистерезис" (отставание, запаздывание) используется физиками в ситуациях, когда величина, характеризующая состояние объекта наблюдения, неоднозначно зависит от внешних условий. Например, намагниченность образца ферромагнетика может быть разной при одной и той же индукции внешнего магнитного поля. У нас с этим словом связано то, что при попадании в область бистабильности на рис. 3, в, д слева устанавливается один цикл, а справа - другой. Другими словами, при увеличении N реализуется нижняя ветвь резонансной кривой, а при уменьшении - верхняя. В областях мультистабильности можно встретить более сложные варианты гистерезиса.
ХАОС
При достаточно больших b, когда амплитуда колебаний с ростом А становится значительной, в системе устанавливаются различные виды хаотических колебаний. Кроме беспорядочности такие движения обладают большой чувствительностью к изменению начальных условий. В этом можно убедиться, если, не меняя значений параметров, повторить расчеты xn(n) после сколь угодно малого изменения начальных условий x0 . В результате обнаруживается быстрое нарастание расхождения между новой и старой реализациями с течением времени, хотя вначале они почти совпадают. Этот факт очень важен для понимания феномена динамического хаоса [3]. Он был осознан научным сообществом сравнительно недавно, и открытие его признано одним их важнейших достижений науки в XX веке. Хаос занимает большую часть плоскости параметров (см. рис. 4). Двигаясь по карте режимов (рис. 4) из области периодических колебаний к хаосу, не пересекая границ листов (см. черно-белый фрагмент), можно убедиться, что хаосу предшествует последовательность бифуркаций удвоений периода некоторого цикла.
ЗАКЛЮЧЕНИЕ
Рассказанное про маятники на примере модели с нелинейностью типа мягкой пружины в той или иной степени присуще различным диссипативным системам, совершающим вынужденные колебания около положения равновесия. Об этом говорит сравнение с результатами численных исследований дифференциальных уравнений различных нелинейных осцилляторов и данными физических экспериментов. Так, универсальные элементы приведенной на рис. 4 структуры присутствуют на аналогичных плоскостях параметров маятников с жесткой пружиной и наблюдаются при других способах внесения воздействия.
В то же время приведенная картина не исчерпывает многообразия возможных вариантов нелинейного поведения маятников. Например, если смоделированный нами маятник подталкивать не в одну сторону от положения равновесия, а симметрично - в одну и другую, то перед удвоениями периода колебаний происходит так называемая бифуркация потери симметрии. При этом вместо одного симметричного движения в системе рождаются два зеркально-симметричных движения-близнеца. Естественно, что поведение более сложной колебательной системы может быть еще более изощренным. Поэтому, разглядывая переплетение цветов и линий на рис. 4, невольно задумываешься о сложности мироустройства, если примитивный груз на стержне ведет себя столь непросто.
Не так давно в журнале "В мире науки" вышла статья о хаотических колебаниях в контуре с диодом, где отмечалась доступность нелинейного электрического маятника и говорилось, что он "должен стоять в каждом доме на подоконнике". Каждому дому навряд ли, но физическим кабинетам и практикумам он действительно очень подходит (практические рекомендации можно найти в [5]). Что же касается простых модельных отображений, то их исследования вполне по силам каждому студенту или сильному ученику-старшекласснику. Запрограммируйте формулу (4) и с помощью карты режимов отправляйтесь в увлекательную прогулку по просторному зоопарку нелинейных колебаний маятника.
ЛИТЕРАТУРА
1. Трубецков Д.И. Колебания и волны для гуманитариев. Саратов: ГосУНЦ "Колледж", 1997. 392 с.
2. Данилов Ю.А. Нелинейность // Знание - сила. 1982. ╧ 11. С. 34.
3. Анищенко В.С. Детерминированный хаос // Соросовский Образовательный Журнал. 1997. ╧ 6. С. 70-76.
4. Белых В.Н. Элементарное введение в качественную теорию и теорию бифуркаций динамических систем // Там же. ╧ 1. С. 115-121.
5. Безручко Б.П., Жалнин А.Ю., Прохоров М.Д., Селезн╦в Е.П. Дискретные нелинейные модели периодически возбуждаемой RL-диод цепи // Изв. вузов. Прикл. нелинейная динамика. 1997. Т. 5, ╧ 2/3. С. 48-62.
Рецензент статьи В.П. Крайнов
* * *
Борис Петрович Безручко, доктор физико-математических наук, профессор кафедры электроники, колебаний и волн Саратовского государственного университета. Область научных интересов - теория колебаний и волн, нелинейная динамика, электроника сверхвысоких частот. Автор более 100 научных публикаций.