Формулировка задания
Смоделировать падение тела в поле тяжести Земли с учетом сопротивления воздуха.
Тело бросается с высоты H с начальной скоростью v0, направленной под углом 0 к горизонту. Помимо силы тяжести на тело действует сила сопротивления воздуха, имеющая вид Fc = kSv2, где k – некоторый численный коэффициент, зависящий, в частности, от формы тела и качества обработки его поверхности, - плотность воздуха, S - площадь сечения тела, v - скорость тела. Ускорение свободного падения g = 9.81 м/c2.
Для примера рассмотрите полет ядра, используя следующие значения параметров: m = 7 кг, S = 0.01 м2, H = 2 м, v0 = 30 м/c, = 1.18 кг/м3.
Выберите шаги h для метода Эйлера и для модифицированного метода Эйлера, решив задачу без учета сопротивления среды (k = 0).
Расстояние, на которое улетит тело по горизонтали без учета сопротивления среды определяется формулой
.
Максимум дальности полета достигается при угле бросания = 45. Исследуйте влияние сопротивления среды на дальность полета ядра и на оптимальный угол бросания. Для этого проведите серию расчетов без учета и с учетом сопротивления среды при различных углах бросания, задав коэффициент сопротивления k = 0.1.
Исследуйте влияние попутного и встречного ветра, направленного вдоль оси x, на дальность полета ядра и на оптимальный угол бросания.
Физическая модель
Тело имеет массу m и движется под действием двух сил: силы тяжести mg, направленной вертикально вниз, и силы сопротивления воздуха Fc, направленной в сторону противоположную скорости. Предполагаем, что начальная скорость тела достаточно мала, чтобы можно было пренебречь изменением величин g и при перемещении тела. Ускорение тела находится из второго закона Ньютона.
Математическая модель
Пусть ось x направлена горизонтально в направлении броска, а ось y - вертикально вверх. Поскольку обе силы, действующие на тело, лежат в плоскости (x,y), то траектория тела также будет лежать в этой плоскости. Используя второй закон Ньютона, получаем систему из четырех дифференциальных уравнений:
Тело бросают в момент времени t = 0 с начальными условиями
Численный алгоритм
Для нахождения решения полученной системы дифференциальных уравнений используем метод Эйлера с постоянным шагом по времени t. Получаем систему разностных уравнений:
Визуализация результатов расчетов
Разделим экран на четыре окна. В первом окне будем изображать траекторию движения тела y(x), в остальных окнах построим графики зависимостей от времени y(t), v(t) и полной механической энергии тела
В качестве теста используем вариант с отсутствием силы сопротивления воздуха. В этом случае механическая энергия в точной постановке задачи сохраняется.
Для большей наглядности при исследовании зависимости движения от параметров результаты расчетов, полученные для разных значений параметров, отображаются одновременно на экране разными цветами.
При создании проекта необходимо рассмотреть следующие пункты:
Формулировка задания
Физическая модель
Математическая модель
Численный алгоритм
Компьютерное моделирование
Визуализация результатов расчетов
Интерпретация результатов
В формулировке задания необходимо отразить содержание задачи в целом.
1,2,3,4,7 – должны содержаться в отчете о разработке проекта.
- Пример выполнения задания
- Формулировка задания
- Примеры тем для базовых проектов
- Изучение движения баскетбольного мяча
- Стрельба из сверхдальнобойной пушки
- Прыжок с парашютом
- Скатывание тела по сложному рельефу
- Скоростной спуск «Американские горки»
- Движение яхты
- Карусель смеха
- Моделирование взлета с Земли многоступенчатой ракеты
- Торможение космического аппарата при входе в атмосферу
- Посадка на Луну
- Движение частицы в центральном поле
- Распад спутника на орбите
- Система Земля, Луна, спутник Луны
- Полет на Марс
- Солнечная система
- Движение планеты в гравитационном поле двойной звезды
- Движение двух планет вокруг массивной звезды
- Рассеяние электрона на атоме водорода
- Гравитационная машина
- Формирование планетной системы из газового облака
- Формирование спиральных рукавов в галактиках
- Физический маятник
- Движение маятника в магнитном поле
- Исследование явления резонанса
- Шарик на пружинках
- Взаимодействие двух маятников
- Волны в упругой среде
- Движение биллиардных шаров
- Релаксация энергии в системе двухатомных молекул
- Силовые линии и эквипотенциали для системы зарядов
- Колебательный контур
- Циклотронный резонанс
- Демонстрация эффекта Доплера
- Остывание пластины
- Ядерный реактор