ЧМТФ Лабораторная работа 6. Лабораторная работа Моделирование траекторий космических тел
 Скачать 0.97 Mb.
|
|
Лабораторная работа 6. Моделирование траекторий космических тел Цель работы:
Теоретические сведения Между любыми двумя материальными точками действуют силы взаимного притяжения, прямо пропорциональные произведению масс этих точек и обратно пропорциональные квадрату расстояния между ними (рис. 6.1)
где  - сила тяготения, действующая на точку с массой  ,  - радиус-вектор, проведенный из этой точки в точку, обладающую массой  ,  - расстояние между точками. Коэффициент  называется постоянной всемирного тяготения.Из рис. 6.1 видно, что вектор можно представить как разность радиус-векторов точек и , то есть  .Для компьютерного моделирования удобнее пользоваться следующей формулой закона всемирного тяготения  .Приведенные формулы справедливы не только для материальных точек, но и для двух твердых тел шарообразной формы. В этом случае - радиус-вектор, соединяющий центры первого и второго тел. Формула (6.1) выполняется также в том случае, когда одно из тел имеет произвольную форму, но его размеры во много раз меньше размеров второго тела.Потенциальная энергия взаимного притяжения двух материальных точек (или шаров) определяется по формуле  .Если  ,то вторая частица практически неподвижна и располагается в центре масс системы, а первая частица движется в центрально-симметричном поле, в котором выполняются законы сохранения энергии и момента импульса. Траектория частицы определяется величиной эксцентриситета  , где Eи L - полная энергия и момент импульса материальной точки соответственно.Компьютерное моделирование движения частиц, движущихся в плоскости, производится аналогично, как и рассмотренные ранее случаи одномерного движения, но при этом необходимо учесть следующие особенности:
Пример. Вычислить траектории движения двойной звезды, состоящей из двух звезд с массами  кг и  кг (рис. 6.2). В начальный момент времени расстояние между ними равно  м. Скорость первой звезды  м/с и направлена под углом  к отрезку, соединяющему звезды в начальный момент времени. Вычисления произвести в системе, связанной с центром масс двойной звезды (ц-системе), а также в системе, относительно которой центр масс движется горизонтально со скоростью  м/с (л-системе).
Решение. На рис. 6.3 приведен фрагмент листинга программы для ввода данных и начальных условий задачи.
Первая строка рис. 6.3 соответствует данным из условия задачи. В системе отсчета, связанной с центром масс (ц-системе), общий импульс равен нулю, отсюда можно найти начальную скорость второй звезды. Полученная формула приведена во второй строке. Во второй строке также введен малый промежуток времени  и опция, показывающая, что нумерация строк и столбцов матриц начинается с единицы. В третьей строке рис. 6.3 введены радиус-векторы и векторы скоростей рассматриваемых частиц в начальный момент времени в ц-системе.Для решения задачи в л-системе, согласно принципу относительности Галилея, нужно изменить начальные скорости (рис. 6.4).
Процедура решения задачи по методу Эйлера приведена на рис. 6.5.
Выходным параметром процедуры является матрица, состоящая из четырех столбцов, в которых располагаются декартовы координаты (компоненты радиус-векторов) звезд в различные моменты времени. Для построения графиков из матриц выделим соответствующие колонки
Траектории звезд, вычисленные в ц-системе, представлены на рисунке 6.7, в л-системе - на рис 6.8.
Порядок выполнения работы Задание 1. Приближение к Луне Ракета массой mподлетает к Луне (рис. 6.9 - 6.14). Когда расстояние до Луны становится равным а, а прицельное расстояние  , скорость объекта равна  . Рассчитать траекторию полёта ракеты вблизи Луны. Масса Луны равна  кг, радиус Луны равен  км. Изобразить траектории ракеты при  . При расчётах принять, что малым является промежуток времени  с.
Данные для расчета взять из табл. 6.1.
Для каждого значения начальной скорости вычислить значения энергии, момента импульса и эксцентриситета. Провести исследование, как зависит вид траектории от значений энергии и эксцентриситета. Результаты занести в табл. 6.2.
Задание 2. Посадка спутника в атмосфере. Ракета массой mприближается из космоса к планете массой Mи радиусом Rс начальной скоростью V0 (рис. 6.9 - 6.14). Начальное положение характеризуется параметрами а и p0. На планете существует атмосфера, эффективная толщина атмосферы которой H. Сила сопротивления в атмосфере характеризуется по модулю силой  и направлена в сторону, противоположную движению, где h- расстояние до поверхности. Найти координаты приземления ракеты без учёта сопротивления (x1, y1) и с учётом сопротивления (x2, y2). При расчётах принять, что малым является промежуток времени  с. Данные для расчета взять из табл. 6.3.
Задание 3. Построение траектории Луны. Построить траектории Земли и Луны в гелиоцентрической системе координат с использованием компьютерного моделирования. Начальную конфигурацию системы «Солнце (S) - Земля (Z) - Луна (L)» принять изображенной на рис. 6. 15. Время моделирования - 30 суток, малым является промежуток времени - 1 мин. При расчетах принять: Масса Солнца  кг, масса Земли  кг, масса Луны  кг, Среднее расстояние от Солнца до Земли  м, среднее расстояние от Земли до Луны  м, период обращения Земли вокруг Солнца  суток, период обращения Луны вокруг Земли  суток.
Контрольные вопросы
Как зависит траектория движения частицы в центрально симметричном поле от массы, прицельного расстояния? |
||||||||||||||||||||||||||||||||||||||||
,
