- •А.И. Ходанович математическое моделирование на компьютере
- •История математического моделирования на компьютере
- •Задачи аналитического и численного моделирования
- •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
- •Литература
- •Ходанович
2.4 Баллистическая кривая
Тело брошено под углом к горизонту с начальной скоростью. Данные ввести самостоятельно. Построить траекторию, если сила сопротивления F=0.001v2. Сравнить ее с параболической траекторией, получающейся без учета сопротивления воздуха. Показать на экране траектории при разных параметрах задачи.
Рис.5. Моделирование баллистической кривой
Указание. Уравнения движения по координатным осям имеют вид [8]:
2.5 Модели космической динамики
Пусть скорость тела v перпендикулярна к направлению на центр Луны (а=0). Прицельное расстояние р=5000 км. Как будет двигаться тело при скоростях 500, 1000, 1500 м/с ? Радиус Луны равен 1700 км.
Рис. 6. Динамика в гравитационном поле
Указание. Необходимые формулы численного моделирования:
Задавая различные значения величин р,а,v, исследуйте вопрос о движении тела вблизи Луны или Земли. Постройте траектории движения и определите космические скорости.
Исследуйте влияние сопротивления воздуха на движение спутника, вращающегося вокруг Земли. Предположите, что сила сопротивления пропорциональна квадрату скорости. Как меняются со временем полная энергия и скорость спутника ?
2.6 Посадка «Гюйгенса» на Титан
Уникальный эксперимент с посадкой зонда «Гюйгенс» на поверхность загадочного спутника Сатурна стал кульминационным моментом всей беспрецедентной для современной космонавтики программы исследований Сатурна и его системы исследовательской станцией «Кассини-Гюйгенс» (2005 г.)
Зонд вошел в атмосферу планеты. В это время он находился на расстоянии около 1270 км от поверхности Титана. Скорость летательного аппарата в начальный момент времени составляла примерно 6000 м/с. После четырехминутного снижения на высоте ~180 км и на скорости 400 м/с раскрылся первый из трех парашютов необходимых для снижения скорости зонда до расчетной, а также для приведения в дейст-вие механизмов отстрела защитных поверхностей, предохранявших зонд и научную аппаратуру от перегрева.
Рис.7. Схема запуска зонда «Гюйгенс» на Титан
На протяжении всего спуска в атмосфере, который, по оценкам ученых, должен был продолжаться примерно два часа, «Гюйгенс» передавал научную информацию на борт «Кассини», который затем должен был ретранслировать ее на Землю.
По данным источников информации рассчитать траекторию зонда за первые 4 мин. с момента вхождения летательного аппарата в атмосферу Титана до раскрытия первого парашюта.
Рис.8. Иллюстративный материал к
учебному вычислительному эксперименту
Указание. Рассмотреть физическую модель падения тела в поле тяготеющей массы с учетом силы сопротивления среды:
.
m – масса «Гюйгенса»;
M – масса Титана;
х – расстояние от входа в атмосферу до поверхности Титана;
h – расстояние от раскрытия первого парашюта до поверхности Титана;
Rт – радиус Титана;
F – сила тяготения;
– сила сопротивления среды.
Данные, найденные в источниках информации о посадке Гюйгенса на Титан необходимые при решении данной задачи:
G= 6,67∙10-11(ед.СИ)- гравитационная постоянная;
Rт= 2575 км – радиус Титана;
М= 1.8∙Млуны=1.8∙7.33∙1022 кг =1.35∙1023 – масса Титана;
m= 320 кг – масса Гюйгенса;
х=1270 км – расстояние от входа в атмосферу до поверхности Титана;
h= 180 км – расстояние от раскрытия первого парашюта до поверхности Титана;
= 6000 м/с – скорость, при входе в атмосферу; = 400 м/с – скорость при раскрытии первого парашюта; = 4 мин = 240 с- время до раскрытия первого парашюта (рис. 9).
Рис.9. Расчет динамики летательного аппарата на основе
данных астрофизических наблюдений