- •«Сибирский государственный аэрокосмический университет им. Академика м.Ф. Решетнева»
- •1. Общие понятия.
- •1.1. Гипотеза непрерывности среды.
- •1.2. Основные физические свойства жидкости.
- •2. Предмет курса механики полёта беспилотных летательных аппаратов
- •2.1. Этапы исследования полёта летательных аппаратов.
- •2.2. Общие сведения об управлении полётом.
- •2.3. Силы, действующие на летательный аппарат.
- •Фиг.1 Силы, действующие на летательный аппарат в полёте
- •Силы, действующие на летательные аппараты при прямолинейном полёте.
- •Аэродинамическая нормальная сила.
- •Силы, действующие на летательный аппарат при горизонтальном манёвре без скольжения (вид сзади),
- •Аэродинамически осесимметричные схемы летательных аппаратов.
- •Нормальная сила, создаваемая двигателем.
- •Пример установки стартового двигателя под углом к оси летательного аппарата.
- •3. Основные соотношения теории реактивного движения.
- •3.1. Движение точки переменной массы. Тяга реактивного двигателя.
- •Тяга реактивного двигателя
- •Фиг.11 Силы, действующие на закреплённую ракету.
- •Фиг.12. Распределение сил давления по поверхности ракеты и по внутренней поверхности камеры реактивного двигателя.
- •Удельная тяга или удельный импульс
- •4. Органы управления летательными аппаратами.
- •Управление с помощью аэродинамических сил.
- •Фиг.14. Воздушные рули на задних кромках крыльев.
- •Фиг.15. Основные типы интерцепторов.
- •Фиг.16. Управляемый снаряд с поворотными крыльями
- •Фиг.17. Воздушные рули в схеме «утка»
- •Фиг.18. Пример схемы «бесхвостка»
- •Управление с помощью реактивных сил
- •Фиг.19. Схема управления летательным аппаратом посредством поворота ракетного двигателя.
- •Управление креном
- •Фиг.19. Силы, возникающие при отклонении элеронов
- •5. Система управления летательным аппаратом.
- •6. Уравнения движения жидкости и газа. Законы истечения.
- •6.1. Уравнение постоянства расхода (уравнение неразрывности).
- •6.2. Дифференциальные уравнения движения идеальной жидкости (уравнение Бернулли).
- •6.3. Геометрический и энергетический смысл уравнения Бернулли.
- •6.3.1. Геометрический смысл уравнения Бернулли.
- •6.3.2. Физический смысл уравнения Бернулли.
- •6.4. Уравнение Бернулли для потока реальной (вязкой) жидкости.
- •6.5. Применение уравнения Бернулли на практике.
- •Трубка полного напора (трубка Пито)
- •7. Аэродинамика и газодинамика. Характеристики потока.
- •7.1. Стандартная атмосфера (са)
- •7.2. Сжимаемость газов. Скорость распространения звука в газе.
- •Фиг. 20 Распространение слабых возмущений в неподвижной среде.
- •Фиг. 21 Распространение слабых возмущений в дозвуковом потоке.
- •Фиг. 22. Распространение слабых возмущений в сверхзвуковом потоке.
- •Фиг. 24. Огибающая семейства ударных волн в сверхзвуковом потоке. Приемник воздушного давления (пвд).
- •7.3 Теория ламинарного течения в круглых трубах.
- •7.3.1 Расход при ламинарном течении жидкости
- •7.3.2. Средняя скорость
- •7.3.3. Потери напора на трение
- •7.4. Турбулентный режим движения жидкости.
- •7.4.1. Особенности турбулентного течения. Пограничный слой.
- •8. Ламинарный и турбулентный пограничные слои.
- •8.1. Основные понятия пограничного слоя.
- •8.2. Выводы по разделу
- •8.3. Ламинарный и турбулентный режимы течения в пс.
- •Фиг. 27. Структура пс при переходе ламинарного течения в турбулентное.
- •8.3. Интегральное соотношение для установившегося течения в пограничном слое несжимаемой жидкости.
- •8.4. Пограничный слой и сопротивление трению плоской пластины в несжимаемой среде для ламинарного режима течения.
- •Фиг.31. Изменение напряжения трения и толщины пограничного слоя по длине пластины.
- •8.5. Отрыв течения в пограничном слое и образование вихрей. Аэродинамически удобообтекаемые и неудобообтекаемые тела.
- •Фиг.34. Развитие области отрывного течения за цилиндром.
- •Фиг.35. Развитие области отрывного течения за сферой.
- •Фиг.37. Зависимость коэффициента лобового сопротивления Сх от числа Rе для удобообтекаемого тела.
- •9. Элементы теории подобия потоков.
- •9.1. Геометрическое, кинематическое и динамическое подобие. Коэффициенты подобия.
- •9.2. Полное и частичное динамическое подобие. Критерии динамического подобия.
- •10. Динамика полёта.
- •10.1 Предмет и задачи курса динамики полёта.
- •10.2. Системы координат.
- •Фиг. 41 Геоцентрическая система координат
- •Фиг. 42 Геодезическая система координат
- •Фиг. 43 Земная система координат
- •Фиг. 44 Стартовая система координат
- •Фиг. 45 Ориентация связанных осей на старте летательного аппарата
- •Начальная стартовая система координат.
- •Фиг. 46 Местная географическая система координат
- •Геоцентрическая система координат.
- •Скоростная и полускоростная системы координат.
- •10.3. Косинусы углов между осями систем координат
- •10.3.1. Косинусы углов между осями связанной и начальной стартовой систем координат.
- •Фиг.47 Переход от начальной стартовой системы координат к связанной.
- •Фиг.48 Последовательные повороты на углы:
- •10.3.2. Направляющие косинусы между осями земной и связанной систем координат.
- •Фиг.49 Переход от земной системы координат к связанной системе координат.
- •10.3.3 Косинусы углов между осями полускоростной системы координат и местной географической системой координат.
- •Фиг.50 Переход от местной географической системы координат к полускоростной системе координат.
- •10.3.4. Косинусы углов между связанной и скоростной (поточной) системами координат.
- •Фиг.51 Связанная система осей координат.
- •Фиг.52 Скоростная система координат .
- •Фиг.53. Переход от скоростных осей к связанным осям.
- •11. Уравнения движения ла.
- •11.1. Принцип составления уравнений полёта реактивного ла.
- •11.1.1. Теорема о количестве движения.
- •11.1.2. Теорема о моменте количеств движения.
- •11.1.3. Принцип затвердевания.
- •11.1.4. Сила тяги реактивного двигателя (рд).
- •11.2. Уравнения движения ла в векторной форме
- •11.2.1. Уравнения движения центра масс.
- •11.2.2. Уравнение вращательного движения относительно центра масс.
- •11.2.3. Векторные уравнения движения ла относительно Земли.
- •11.3. Уравнения движения летательного аппарата в скалярной форме.
- •11.4. Общая система уравнений движения летательного аппарата.
- •Кинематические уравнения движения центра масс летательного аппарата.
- •11.1.5.Связи, накладываемые на движение летательного аппарата системой управления.
- •11.1.6. Уравнения системы стабилизации.
- •11.1.7. Уравнения системы наведения.
- •12. Траектории движения летательного аппарата.
- •12.1. Траектории полета баллистических ракет(бр) и ракет-носителей(рн).
- •12.2. Участки траектории полета баллистической ракеты и рн.
- •12.2.1 Участок выведения. Номинальные параметры и возмущенное движение
- •Фиг.65 Функции номинального и возмущенного движений.
- •Уравнение баллистики.
- •Фиг.66 Отсчёт координат от теоретической вершины ракеты.
- •Программа выведения.
- •12.2. Полёт летательного аппарата в центральном поле тяготения за пределами атмосферы.
- •Фиг.68 к выводу уравнений движения в полярной системе координат.
- •Траектория движения в общем виде
- •Фиг.69 Траектории свободного полета при различных скоростях выведения.
- •12.3. Атмосферный участок (входа в атмосферу).
- •12.4. Уравнения возмущенного движения.
- •12.6. Передаточные функции и их свойства.
- •12.7. Частотные характеристики и частотный критерий устойчивости.
- •12.8.Структура автомата стабилизации.
- •12.9 Эффективность органов управления.
Уравнение баллистики.
Посмотрим, как изменяются уравнения (144)- (149) при отсутствии возмущений и какие соотношения устанавливаются между номинальными параметрами, определение которых и представляет собой предмет баллистики.
Фиг.66 Отсчёт координат от теоретической вершины ракеты.
Наибольшие изменения
претерпевает уравнение (146). Т.к.
номинальный угол тангажа φн
изменяется весьма медленно, то в О
обращается левая часть уравнения (
).
Так же могут считаться равными О
кориолисов момент Mкор
и демпфирующий момент Мд,
величины которых пропорциональны малой
величине
,
поэтому
.
Согласно обозначениям на фиг. 66 получим:
.
Переходя от скоростной системы координат к связанной, в уравнениях (144) и (145) силы X и Y можно заменить осевой X1 и нормальной Y1 силами:
;
.
Исключая Yp , получим:
Согласно выражениям (148) и (149) имеем:
и
.
Исключая
,
приводим уравнения к окончательному
виду:
Эти уравнения дополняются еще двумя (148):
;
.
Уравнения (149) дают:
,
где закон изменения угла тангажа φ(t) в простейших случаях может считаться заданной функцией φпр(t), называется программой угла тангажа.
В дополнение к
этим уравнениям для окончательного
суждения о форме траектории есть
выражения ортодромной дальности
и местной высоты H=r-R
.
Интегрирование полученной системы уравнений дает закон движения ракеты как материальной точки. В результате получаем номинальные параметры траектории центра масс ракеты, определение которых и представляет собой основную задачу баллистических расчетов.
Программа выведения.
Программа выведения
– совокупность управляющих функций,
которые регламентируют поведение
вектора тяги в номинальном движении,
задают закон изменения секундного
расхода
и угловую ориентацию ракеты в пространстве
и времени. Программа выведения автономно
выполняется системой управления в
соответствии с функциями, выбранными
на стадии проектирования ракеты.
Эти функции суть: закон изменения расхода и законы изменения углов тангажа φпр(t) , рысканья ψпр(t) и крена γпр(t) на участке выведения.
Рассмотрим пример построения простейшей программы выведения для одноступенчатой ракеты.
Чтобы не устанавливать ракету на старте по линии прицеливания и освободить стартовую систему от тяжелых поворотных механизмов, вводится программа разворота по крену. Ракета на стартовой позиции устанавливается в раз и навсегда принятом положении (ориентация плоскости 1-3 стабилизатор), независимо от направления предстоящего пуска. Азимут вычисляется заранее и передается бортовому запоминающему устройству (БЗУ). После старта ракеты (на вертикальном участке траектории) автомат стабилизации разворачивает ее относительно вертикальной оси на линию прицеливания. На это уходит 10 – 15 сек., после чего начинается разворот ракеты по тангажу.
Угол ψпр(t) при старте ракеты с Земли, как правило, равен 0, но чтобы обеспечить, например, падение отделяемых элементов конструкции летательного аппарата в заданный район, может быть введена программа изменения угла рысканья.
Выбор программы выведения представляет собой поиск некоторого оптима. Но на этот поиск наложены ограничения, связанные с конструктивными и эксплуатационными особенностями ракеты.
Баллистические ракеты стартуют вертикально. Из этого вытекает начальное условие для изменения угла тангажа. Необходимо, чтобы при t=0 φпр=90о ,и продолжительность вертикального подъема должно быть не меньше того времени, которое потребуется двигательной установке (ДУ) для выхода на режим полной тяги. В случае программного разворота по крену надо выждать, пока и эта операция не будет закончена полностью.
Но и с разворотом
по углу тангажа нельзя медлить. Увеличение
вертикального подъема приводит к отходу
от оптимальных условий выведения, хотя
бы по дальности. Но главное- с увеличением
времени вертикального участка поворот
ракеты придется проводить при больших
скоростях полета, т.е. при большом
скоростном напоре
,
при опасных углах атаки. Каждому
скоростному напору для данной ракеты
соответствует свое максимально допустимое
значение угла атаки, величина которого
определяется прочностью корпуса на
изгиб под действием поперечного
нагружения. Чем больше скоростной напор,
тем больше допустимый угол атаки. Значит,
разворот ракеты надо начинать, пока
скоростной напор мал.
При скорости, приближающейся к скорости звука, происходит быстрое изменение аэродинамических сил и моментов, смещается центр давления. Это создаёт трудности для автомата стабилизации, поэтому этот участок траектории желательно проходить при нулевых углах атаки. Этого требуют и условия прочности, естественно.
На участке 0-t1 угол φпр остаётся неизменным (φпр=900). С момента t1 начинается разворот ракеты (уменьшается φпр) и возникает угол атаки. Надо сделать так, чтобы максимальное значение (по абсолютной величине) αmax достигалось в начале участка разворота , а при достижении скорости близкой к скорости звука, движение продолжалось с нулевым углом атаки. Ограничение по углу α может быть снято лишь после того, как ракета выйдет за пределы атмосферы, где скоростной напор ничтожно мал.
Для уравнения (153) α=φ=υ=χ, безразлично, какой из двух углов задан – φпр или α, а какой определяется в результате интегрирования. Поэтому, на стадии выбора программы вывода удобно аппроксимировать угол атаки на участке разворота функцией α(t), удовлетворяющей наложенным ограничениям, а программный угол получить в результате интегрирования уравнений движения.
Примером такой аппроксимирующей функции может служить быстрозатухающая функция вида:
.
Параметр определяет скорость затухания и выбирается таким, чтобы при числах М = 0,7…0,8 угол атаки был примерно равным нулю.
Неопределённым остаётся максимальный угол атаки αmax. Это и есть основной варьируемый параметр программы от которого зависит закон последующего движения.
Приняв указанную аппроксимирующую функцию для α, можно при некотором фиксированном значении угла αmax произвести интегрирование уравнений движения и получить закон изменения программного угла φпр (фиг.67).
Когда ракета выйдет за пределы атмосферы, ограничения на угол атаки снимаются, а с момента t=t2 можно принять φпр= φпр.к. неизменным На программе устанавливается “полка”.
При неизменном φпр угол атаки возрастает, а угол наклона траектории υ=φпр+ χ – α уменьшается. К моменту выключения двигательной установки t=tк угол υк должен соответствовать заданному краевому условию баллистической задачи.
В итоге получена довольно простая программа угла тангажа с двумя варьируемыми параметрами – αmax и φпр.к. (или t2, что одно и то же). Наличие свободного параметра φпр.к позволяет, помимо заданного краевого условия, выполнить ещё одно из экстремальных условий – обеспечить максимальную дальность или минимальное рассеивание.
В заключение отметим, что программный угол φпр и фактически реализуемый угол тангажа φ, различаются. Это различие зависит от статической характеристики автомата стабилизации и учитывается в уточнённых баллистических расчётах.
