
- •А.И. Ходанович математическое моделирование на компьютере
- •История математического моделирования на компьютере
- •Задачи аналитического и численного моделирования
- •1. Симметрия и модели с числовыми последовательностями
- •1.1 Формула пути равнопеременного движения
- •1.3 «Золотое сечение» и числа Фибоначчи
- •1.4 Сопротивление бесконечной электрической цепи
- •1.5 Сопротивление конечной электрической цепи
- •1.6 «Золотое сечение» в поле тяготения
- •Равенство тяжелой и инертной массы в свободном падении
- •2.2 Реальная динамика тела в среде
- •2.3 Цепная линия и метод кратных скобок Пуассона
- •2.4 Баллистическая кривая
- •2.5 Модели космической динамики
- •2.6 Посадка «Гюйгенса» на Титан
- •Задача Кеплера
- •Слайд 10. Компьютерно-графическое моделирование в небесной механике
- •2.8 Орбита в пространстве скоростей
- •2.9 Конические сечения в поле тяготения
- •2.10 Прецессия орбиты при малом возмущении
- •2.11 Задача трех тел в небесной механике
- •2.12 Задача о брахистохроне
- •Слайд 11. Иллюстрация экстремальной траектории в поле тяготения
- •2.13 Управление таймером в режиме реального времени
- •2.14 Физический опыт и модель измерений
- •2.15 Водяные часы
- •2.16 Фракталы в комплексной плоскости
- •Слайд 14. Комплексные отображения в нелинейной динамике
- •3. Колебания и адиабатические инварианты
- •3.1 Гармонический осциллятор
- •3.2 Задача «Бездонный колодец»
- •3.3 Фигуры Лиссажу
- •3.4 Параметрические колебания
- •3.5 Маятник Капицы
- •3.6 Хаотическое поведение маятника и отображение Пуанкаре
- •3.8 Задача с соударениями
- •3.9 Хаотические колебания при соударениях
- •3.10 Представление колебаний рядом Фурье
- •3.11 Два тела на пружине
- •3.12 Цепочка Ферми-Паста-Улама
- •3.13 Линейная цепочка связанных осцилляторов
- •3.14 Колебания мембраны в интерактивной графике
- •3.15 Вращение частицы в трехмерной графике
- •4. Вероятностно-статистическая линия математического моделирования
- •4.1 Нормальное распределение и статистический критерий
- •4.2 Распределение расстояний между молекулами газа
- •4.3 Радиоактивный распад изотопа
- •4.4 Флуктуации и биномиальный случай
- •4.5 Статистическое распределение с бесконечной дисперсией
- •4.6 Закон больших чисел и центральная предельная теорема
- •4.7 Статистические усреднения в эксперименте
- •4.8 Диффузия в модели случайных блужданий
- •4.9 Иерархия временных масштабов в броуновской динамике
- •4.10 Метод Монте-Карло с генератором случайных чисел
- •4.11 Алгоритмы числа
- •4.12 Модели плоских фигур
- •4.13 Площадь круга с малым возмущением
- •4.14 «Игла Бюффона»
- •4.15 Выделение сигнала на фоне «шумов»
- •4.16 Регрессионная модель эксперимента и метод наименьших квадратов
- •4.17 Модели линейной регрессии в планировании эксперимента
- •5. Дидактические игры математического моделирования
- •5.1 «Телепат»
- •5.2 Игрушка «Раскидай»
- •5.3 Игрушка «Ванька-встанька»
- •5.4 Игра Баше
- •5.5 «Перевернутая» игра Баше
- •5.6 Игра в 15
- •Литература
- •Ходанович
3.4 Параметрические колебания
Если один из параметров колебательной системы (например, приведенную длину маятника) изменять с частотой вдвое большей, чем собственная частота колебаний этой системы, то энергия от источника двойной частоты перекачивается в энергию собственных колебаний.
При модуляции длины подвеса, с частотой, близкой к удвоенной собственной частоте, в отсутствие трения маятник раскачивается, наблюдается параметрический резонанс.
Без ограничения
общности достаточно рассмотреть
уравнение вида:
.
Математическое моделирование
параметрического резонанса приводит
к уравнению Матьё:
при
.
Указание.
В качестве
примера можно рассмотреть «качели» с
собственной частотой математического
маятника
.
Покажите, что если приседать в крайних
точках два раза за период, так что длина
будет в этих точках увеличиваться по
сравнению с начальной длиной, то амплитуда
колебаний (качаний) будет возрастать.
Напишите уравнение вращательного
движения качелей
,
где J-
момент инерции,
-
угол,
-
момент силы. Положите
,
где а<<1.
Рис.22. Эксперименты с параметрическими колебаниями
3.5 Маятник Капицы
Маятником Капицы называется математический маятник, прикрепленный к точке подвеса, которая совершает быстрые вертикальные и горизонтальные осцилляции. При определенных условиях, наложенных на амплитуды и частоты осцилляций, возникают дополнительные положения устойчивого равновесия маятника. Маятник Капицы весьма удобен для изучения явления параметрического резонанса, он достаточно легко воспроизводится в лабораторных условиях и обстоятельно описан теоретически. В 1951 году П.Л. Капица предложил метод усреднения по быстро осциллирующим переменным, суть которого продемонстрируем на примере движения маятника с быстро осциллирующей точкой подвеса.
Случай 1.
Т
очка
подвеса маятника совершает вертикальные
колебания с большой частотой
.
Потенциальная энергия груза задается выражением
,
где m – масса груза, g – ускорение свободного падения, l – длина нити, угловая переменная осцилляций (отклонение от вертикали).
Запишем уравнения движения маятника по осям X и Y:
Определим
лагранжиан:
,
.
Выполним несложные преобразования:
Кинетическая
энергия:
.
Опускаем
полную производную по времени
.
Отсюда
следует, что
Подставляем в выражение для кинетической энергии и опускаем все полные производные по времени:
Таким
образом, лагранжиан равен:
Учитывая,
что
,
а
,
или
Подставляя f в выражение для компоненты, отвечающей за быстрые осцилляции, получаем:
и
находим среднеквадратичную скорость
этих колебаний:
. Таким образом:
.
Каждому устойчивому положению равновесия отвечает минимум эффективной потенциальной энергии.
Условия
глубокой высокочастотной модуляции:
,
приводят к появлению нового устойчивого
равновесного положения
=
(вертикально вверх).
Случай 2
Т
очка
подвеса совершает горизонтальные
колебания с большой частотой.
Запишем
уравнения движения маятника относительно
осей:
,
.
Определим
лагранжиан:
,
.
Выполним
несложные преобразования для выводы
формулы кинетической энергии:
Итак, кинетическая энергия:
Опускаем
полную производную по времени
;
Отсюда
следует, что
.
Подставляем в выражение для кинетической энергии и опускаем все полные производные по времени:
.
Лагранжиан
равен:
.
Учитывая,
что
,
а
,
где
.
Подставляя f в выражение для компоненты, отвечающей за быстрые осцилляции, получаем:
и находим среднеквадратичную скорость этих колебаний.
.
. Таким образом:
.
Условия устойчивости положений равновесия при разных соотношениях параметров модуляции можно сформулировать так:
1.
При
устойчивое
положение
;
2.
При
устойчивое
положение
,
т.е. во втором случае устойчивому
равновесию отвечает значение ненулевого
угла .
Случай 3
Т
очка
подвеса совершает колебания по
вертикальной окружности с постоянной
частотой.
Запишем уравнения движения маятника относительно осей:
Определим лагранжиан: ; .
Выполним несложные преобразования для вывода формулы кинетической энергии:
Кинетическая
энергия:
.
Опускаем
полную производную по времени
;
Отсюда
следует, что
.
Подставляем в выражение для кинетической энергии и опускаем все полные производные по времени:
.
Лагранжиан
равен:
.
Учитывая, что
,
а
,
где
.
Подставляя
f
в выражение для компоненты, отвечающей
за быстрые осцилляции, получаем:
и
находим среднеквадратичную скорость
этих колебаний:
.
Это верно в случае, если вертикальная и горизонтальная компоненты частоты колебания подвеса одинаковы, и амплитуды колебания тоже равны. . Таким образом, для этого случая:
.
Автоматизация вычислений на ПК в задачах теоретической физики и построение поверхностей эффективной потенциальной энергии при движении в быстро осциллирующих полях оказывается весьма полезным инструментом для изучения условий возникновения параметрических колебаний. Трехмерные графики, реализованные в среде компьютерной математики, позволяют исследовать перестройку положений устойчивого равновесия и представить результат в наглядном виде.
Рис.23. Поверхность эффективной энергии