- •А.И. Ходанович математическое моделирование на компьютере
- •История математического моделирования на компьютере
- •Задачи аналитического и численного моделирования
- •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.6 Хаотическое поведение маятника и отображение Пуанкаре
Фазовая плоскость определяется как множество точек (x,v). (Иногда вместо скорости v используется импульс p=mv). Когда движение периодично его орбита описывает в фазовой плоскости замкнутую кривую. Вынужденные колебания нелинейной системы с кубичным упругим элементом могут иметь орбиту с самопересечениями, но, тем не менее, замкнутую. В этом случае следует ожидать присутствие субгармоник.
Напротив, орбиты хаотических движений никогда не бывают замкнутыми, не повторяются. Такие орбиты стремятся заполнить некоторую область фазового пространства. Хотя подобное блуждание орбит указывает на хаос, непрерывные графики на фазовой плоскости мало информативны, и следует использовать более совершенный метод исследования фазовых портретов, называемый отображением Пуанкаре.
При математическом исследовании динамических систем отображением называют временную выборку данных , для которой вводят обозначение Понятие отображения переносится и на боль-шее число переменных. В простом детерминированном случае отображение представляет систему разностных уравнений. Если и , то эта последовательность точек фазового пространства представляет собой двухмерное отображение:
.
Слайд 18. Отображение Пуанкаре в нелинейной динамике
Если моменты выборки подчиняются определенному правилу, например, кратны характерному периоду системы, это отображение называется отображением Пуанкаре [45].
Для многих механических систем характер движения можно определить приблизительно, комбинируя аналитические и численные методы. При определенных условиях детерминированные системы проявляют хаотическое поведение, чувствительное к начальным данным и к приближениям, сделанным в ходе вычисления траектории. В данном случае мы не в состоянии предсказать долговременное поведение динамической системы.
Указание. Исследовать поведение нелинейного маятника на фазовой плоскости с вертикально осциллирующей точкой подвеса. Уравнение движения имеет вид (Маклафлин, Персиваль, Ричардс (1981 г.):
.
3.7 Циклоидальный маятник
|
Обратимся к задаче, решение которой позволило нидерландскому ученому Х. Гюйгенсу сконструировать в 1673 г. «точные» часы. Эта задача заключается в построении в вертикальной плоскости такой кривой, чтобы время, необходимое для спуска по ней до фиксированного горизон |
та материальной точки, находящейся в начальный момент времени в состоянии покоя, не зависело от исходного положения точки на этой кривой. Как показал Гюйгенс, такой изохронной кривой оказалась циклоида: [3].
Циклоида обладает многими интересными свойствами (например, длина арки циклоиды выражается через диаметр целым числом 4), имеющими исключительное значение для физики и техники. В частности, профили зубьев шестерен, очертания многих типов эксцентриков, кулачков и иных деталей машин имеют форму именно циклоиды.
Р ис.24. Изохронные колебания циклоидального маятника
Практически решение этой задачи можно провести, используя парамет-рические уравнения циклоиды: .
Пусть координаты исходного положения шарика, т.е. точки М, а соответствующее ей значение параметра (рис.24).
Когда шарик скатится из положения М в некоторое положение , он снизится по вертикали на расстояние h, которое определится следующим образом: и, соответственно, скорость тела . С другой стороны, и для циклоиды , и поэтому дифференциальное уравнение принимает вид: .
Слайд 19. Колебания на циклоиде
Интегрируя это уравнение в соответствующих пределах, получаем, что . Если не учитывать трение, то из закона сохранения энергии шарик будет совершать колебания с периодом . Отличительным свойством циклоидального маятника по сравнению с обычным (круговым) является то, что период его колебаний не зависит от амплитуды (изохронные колебания).
При малых колебаниях контур циклоиды можно считать круговым и, положив L= 4r, получаем период малых колебаний математического маятника .