Отчет по вычислительному практикуму
«Исследование траектории полета ЛА»
Вариант 2
Столбец 2
Выполнил студент гр. 030301 Путин В.В.
Машиностроительного факультета
Проверила доцент каф. ракетного Фомичёва О.А.
вооружения
Тула 2012
Содержание
Введение 4
1. Физическая модель полета ЛА 5
2. Математическая постановка задачи 6
3. Метод решения задачи 7
4. Результаты расчета 9
5. Анализ результатов 10
6. Библиографический список 11
Введение
Важнейшей особенностью современного этапа научно технического прогресса является стремительно возрастающая роль ЭВМ во всех областях инженерной деятельности - от систем автоматизированного проектирования и управления до контроля технологическими процессами. В настоящее время ЭВМ стали необходимым оборудованием НИИ, конструкторских бюро, заводов. Это позволило от простейших расчетов и оценок различных конструкций или процессов перейти к новому уровню инженерной деятельности - детальному математическому моделированию (вычислительному эксперименту), которое существенно сокращает потребность в натурных экспериментах, а в ряде случаев может их заменить. В связи с этим современный специалист с высшим образованием должен обладать не только высоким уровнем подготовки по профилю своей специальности, но и хорошо знать математические методы решения инженерных задач, ориентироваться на использование вычислительной техники, практически освоить принципы работы на ЭВМ.
Физическая модель полета ЛА
Траектория полета ракеты представляет собой сложную кривую. Для упрощения расчетов не будем учитывать влияние кривизны и вращения Земли на движение ЛА, а также примем, что движение происходит без колебаний вокруг продольной оси ракеты. Траектория полета делится на два участка: активный (с работающим двигателем) и пассивный (с неработающим двигателем).
С учетом принятых допущений центр масс ракеты будет совершать поступательное движение под действием сил, действующих в одной плоскости: Р – силы тяги, Х – силы лобового сопротивления, G – силы тяжести. На рис. 1 показана система сил, приложенных к центру масс ракеты на траектории в произвольный момент времени.
Рис. 1. Схема сил, действующих на центр масс ЛА
ЛА запускается с земли (x 0 ,y 0)=(0,0) под углом θ=θ 0 к горизонту со скоростью V=V 0. Изменением ускорения свободного падения по высоте пренебрегаем (g=const), тягу принимаем постоянной в течение всего времени активного участка (ta), расход массы определяется соотношением mt /ta, где mt – масса топлива.
Методика расчета траектории базируется на основном законе механики (втором законе Ньютона), связывающем ускорение, испытываемое точкой массой m под действием силы F:
.
Если на тело действует несколько сил, то 
– равнодействующая сила:
.
Математическая постановка задачи
В соответствие с принятой физической моделью в основе математической модели, описывающей траекторию полета неуправляемого ЛА, лежит система уравнений, описывающих движение тела переменной массы, брошенного под углом к горизонту.
(1)
где 
– масса ЛА;
– координаты центра масс ЛА на траектории;
– коэффициент аэродинамического сопротивления;
– плотность воздуха;
– площадь миделя (площадь поперечного сечения ЛА);
– стартовая масса ЛА.
Дифференциальные уравнения второго порядка в рассмотренной системе (1) могут быть заменены на дифференциальные уравнения первого порядка, заменив 
и 
.
(2)
Метод решения задачи
Полученную систему дифференциальных уравнений первого порядка c заданными начальными условиями можно решить различными способами.
Одним из наиболее распространенных методов решения обыкновенных дифференциальных уравнений является метод конечных разностей (МКР). Рассмотрим применение МКР для численного решения на ЭВМ простейшего дифференциального уравнения первого порядка:
с начальными условиями X 0, Y (X 0) = Y 0.
Решение 
будем искать в интервале [ X 0 , b ] и будем полагать, что функция на данном интервале удовлетворяет условиям гладкости.
Разобьем область аргумента Х на множество отрезков длиной Δ X и разложим функцию Y в ряд Тейлора в окрестности произвольной точки Xi из области существования функции:
.
Отбрасывая члены ряда, содержащие производные второго и высшего порядков, получаем конечно-разностное выражение первой производной
.
Отсюда 
Вычисляя последовательно от начального значения Y 0 значения Y 1, Y 2, Y 3,... по данной формуле, находим искомое решение.
На рис. 2 показана форма численного решения, получаемого с помощью таких вычислений
Рис.2. Схема приближенного решения методом Эйлера
Данный метод решения обыкновенного дифференциального уравнения носит название метода Эйлера. При достаточно малых величинах шага 
.Чем меньше шаг 
, тем точнее получается приближённое решение. Метод Эйлера дает решение с большой точностью, так как погрешность близка к 0(
) на каждом шаге процесса.
Для решения системы уравнений (2) методом Эйлера запишем ее в конечно-разностном виде:
,
где i – номер текущего шага по времени.
Выразим в этих уравнениях в явном виде 
. (3)
Решая систему уравнений (3) на каждом шаге по времени, можно последовательно вычислить все точки траектории полета ЛА.
Результаты расчета
5.Анализ результатов
Пишите, что разработали численную модель полёта ЛА, потом, что получили его траекторию, не забудьте упомянуть дальность, время, максимальную и конечную скорость полёта.
6.Библиографический список
1. Калиткин Н.Н. Численные методы М.: Наука, 1978.-512с.
2. Вычислительная техника в инженерных и экономических расчетах: учеб. для вузов. - 2-е изд., перераб. и доп./ Петров А.В., Алексеев В.Е., Титов М.А. и др.; Под ред. А.В.Петрова. - М.: Высшая школа, 1984. - 320 с.
3. Архангельский А.Я. Программирование в С++ Builder 6. – ЗАО "Издательство БИНОМ", 2002 г. – 1152 с.