ЧМТФ Лабораторная работа 6

Лабораторная работа 6.

Моделирование траекторий космических тел

Цель работы:

• повторение законов Ньютона и Кеплера;

• знакомство с методами компьютерного моделирования динамики двумерного движения, визуализация движения космических тел.

Теоретические сведения

Между любыми двумя материальными точками действуют силы взаимного притяжения, прямо пропорциональные произведению масс этих точек и обратно пропорциональные квадрату расстояния между ними (рис. 6.1)

,	(6.1)
Рисунок 6.1

где - сила тяготения, действующая на точку с массой , - радиус-вектор, проведенный из этой точки в точку, обладающую массой , - расстояние между точками. Коэффициент называется постоянной всемирного тяготения.

Из рис. 6.1 видно, что вектор можно представить как разность радиус-векторов точек и , то есть .Для компьютерного моделирования удобнее пользоваться следующей формулой закона всемирного тяготения .

Приведенные формулы справедливы не только для материальных точек, но и для двух твердых тел шарообразной формы. В этом случае - радиус-вектор, соединяющий центры первого и второго тел. Формула (6.1) выполняется также в том случае, когда одно из тел имеет произвольную форму, но его размеры во много раз меньше размеров второго тела.

Потенциальная энергия взаимного притяжения двух материальных точек (или шаров) определяется по формуле .

Если , то вторая частица практически неподвижна и располагается в центре масс системы, а первая частица движется в центрально-симметричном поле, в котором выполняются законы сохранения энергии и момента импульса. Траектория частицы определяется величиной эксцентриситета , где E и L - полная энергия и момент импульса материальной точки соответственно.

Компьютерное моделирование движения частиц, движущихся в плоскости, производится аналогично, как и рассмотренные ранее случаи одномерного движения, но при этом необходимо учесть следующие особенности:

1. Начальные условия должны быть представлены в виде двумерных матриц, компонентами которых являются проекции начального радиус-вектора и начальной скорости, рассматриваемой материальной точки на оси декартовой системы координат.

2. Сила, действующая на материальную точку, должна быть векторной величиной.

3. При соблюдении первого и второго условий вычисление ускорения, скорости и радиус-вектора материальной точки в ненулевые моменты времени методом Эйлера MathCAD автоматически производит в векторном виде.

Пример. Вычислить траектории движения двойной звезды, со­стоящей из двух звезд с массами кг и кг (рис. 6.2). В начальный момент времени расстояние между ними равно м. Скорость первой звезды м/с и направлена под углом к отрезку, соединяющему звезды в начальный момент времени. Вычисления произвести в системе, связанной с центром масс двойной звезды (ц-системе), а также в системе, относительно которой центр масс движется горизонтально со скоростью м/с (л-системе).

Рисунок 6.2

Решение. На рис. 6.3 приведен фрагмент листинга программы для ввода данных и начальных условий задачи.

Рисунок 6.3

Первая строка рис. 6.3 соответствует данным из условия задачи. В системе отсчета, связанной с центром масс (ц-системе), общий импульс равен нулю, отсюда можно найти начальную скорость второй звезды. Полученная формула приведена во второй строке. Во второй строке также введен малый промежуток времени и опция, показывающая, что нумерация строк и столбцов матриц начинается с единицы. В третьей строке рис. 6.3 введены радиус-векторы и векторы скоростей рассматриваемых частиц в начальный момент времени в ц-системе.

Для решения задачи в л-системе, согласно принципу относительности Галилея, нужно изменить начальные скорости (рис. 6.4).

Рисунок 6.4

Процедура решения задачи по методу Эйлера приведена на рис. 6.5.

Рисунок 6.5

Выходным параметром процедуры является матрица, состоящая из четырех столбцов, в которых располагаются декартовы координаты (компоненты радиус-векторов) звезд в различные моменты времени. Для построения графиков из матриц выделим соответствующие колонки

Рисунок 6.6

Траектории звезд, вычисленные в ц-системе, представлены на рисунке 6.7, в л-системе - на рис 6.8.

Рисунок 6.7	Рисунок 6.8

Порядок выполнения работы

Задание 1. Приближение к Луне

Ракета массой m подлетает к Луне (рис. 6.9 - 6.14). Когда расстояние до Луны становится равным а, а прицельное расстояние , скорость объекта равна . Рассчитать траекторию полёта ракеты вблизи Луны. Масса Луны равна кг, радиус Луны равен км. Изобразить траектории ракеты при . При расчётах принять, что малым является промежуток времени с.

Рисунок 6.9	Рисунок 6.10	Рисунок 6.11
Рисунок 6.12	Рисунок 6.13	Рисунок 6.14

Данные для расчета взять из табл. 6.1.

Таблица 6.1

Для каждого значения начальной скорости вычислить значения энергии, момента импульса и эксцентриситета. Провести исследование, как зависит вид траектории от значений энергии и эксцентриситета. Результаты занести в табл. 6.2.

Таблица 6.2

Задание 2. Посадка спутника в атмосфере.

Ракета массой m приближается из космоса к планете массой M и радиусом R с начальной скоростью V₀ (рис. 6.9 - 6.14). Начальное положение характеризуется параметрами а и p₀. На планете существует атмосфера, эффективная толщина атмосферы которой H. Сила сопротивления в атмосфере характеризуется по модулю силой и направлена в сторону, противоположную движению, где h - расстояние до поверхности. Найти координаты приземления ракеты без учёта сопротивления (x₁, y₁) и с учётом сопротивления (x₂, y₂). При расчётах принять, что малым является промежуток времени с. Данные для расчета взять из табл. 6.3.

Таблица 6.3

Задание 3. Построение траектории Луны.

Построить траектории Земли и Луны в гелиоцентрической системе координат с использованием компьютерного моделирования. Начальную конфигурацию системы «Солнце (S) - Земля (Z) - Луна (L)» принять изображенной на рис. 6. 15. Время моделирования - 30 суток, малым является промежуток времени - 1 мин. При расчетах принять: Масса Солнца кг, масса Земли кг, масса Луны кг, Среднее расстояние от Солнца до Земли м, среднее расстояние от Земли до Луны м, период обращения Земли вокруг Солнца суток, период обращения Луны вокруг Земли суток.

Рисунок 6.15

Контрольные вопросы

1. Сформулируйте закон всемирного тяготения.

2. Сформулируйте законы Кеплера.

3. Дайте определение первой и второй космических скоростей.

4. Как зависит вид траектории частицы в центрально-симметричном поле от значения полной энергии?

Как зависит траектория движения частицы в центрально симметричном поле от массы, прицельного расстояния?

8