CC BY
50
УДК 729.78
DOI: 10.18698/0236-3941-2017-2-4-14
ДИНАМИКА ОТДЕЛЕНИЯ ГОЛОВНОГО БЛОКА ОТ РАКЕТЫ-НОСИТЕЛЯ В СЛУЧАЕ АВАРИИ НА АТМОСФЕРНОМ УЧАСТКЕ ВЫВЕДЕНИЯ
В.В. Кокушкин1 М.К. Хомяков2 Н.Ю. Овсянникова2 Н.В. Богомолов2
[email protected] [email protected]
1 МГТУ им. Н.Э. Баумана, Москва, Российская Федерация
2 Ракетно-космическая корпорация «Энергия» им. С.П. Королёва, г. Королёв, Московская обл., Российская Федерация
Аннотация
Исследована динамика отделения головного блока вместе со спускаемым аппаратом космического корабля «Союз» в случае аварии ракеты-носителя на атмосферном участке выведения. Выполнен сравнительный анализ различных математических моделей тела переменной массы — отделяемого головного блока, выбрана оптимальная модель, позволяющая проводить сквозное моделирование движения от момента отделения головного блока от ракеты-носителя до участка автономного движения спускаемого аппарата. Проведена оценка безударности отделения головного блока от ракеты-носителя, представлены кинограммы относительного движения, исследован участок его автономного движения. Установлены основные закономерности переходных процессов движения с оценкой устойчивости отделяемого головного блока в потоке. Выполнен статистический анализ кинематических параметров на момент отделения спускаемого аппарата при неблагоприятных условиях отделения от ракеты-носителя
Ключевые слова
Отделяемый головной блок, аварийное отделение, динамика отделения, относительное движение, статическая устойчивость
Поступила в редакцию 04.05.2016 © МГТУ им. Н.Э. Баумана, 2017
Введение. В случае возникновения аварийной ситуации на атмосферном участке выведения пилотируемого космического корабля (КК) «Союз» предусматривается спасение экипажа. Отделяемый головной блок (ОГБ), в состав которого входит спускаемый аппарат (СА), после выключения маршевых двигателей носителя уводится от ракеты-носителя (РН) с помощью твердотопливной двигательной установки (ДУ).
Проектирование и разработка системы аварийного спасения (САС) экипажа (1960-е гг.) проводились в консервативной постановке для узкого диапазона углов атаки при движении ОГБ. В настоящее время выполняется верификация доработок системы, а также оптимизация условий движения, что потребовало создания быстродействующих моделей расчета в статистической постановке.
Аварийный стык
III
В статье рассматривается динамика отделения и автономного полета ОГБ в случае аварии РН на атмосферном участке выведения [1-3]. На этом участке при возникновении сигнала «Авария» КК фиксируется относительно ОГБ в поперечном направлении тремя опорными ложементами САС и подается команда на отключение маршевых двигателей РН. Затем силовые связи по аварийному стыку разрываются (рис. 1), подаются команды на раскрытие решетчатых стабилизаторов (РС) и включение ДУ САС. Под действием тяги центрального ракетного двигателя (ЦРД) ОГБ отделяется от РН, а управляющие ракетные двигатели (УРД) отклоняют головной блок из плоскости стрельбы (I—III плоскость стабилизации). Затем с некоторой задержкой (порядка 20.30 с), определяемой временем аварии по траектории РН, подается команда на отделение СА и задействуется циклограмма его посадки.
Для обеспечения статической устойчивости головного блока в потоке на всем атмосферном участке в диапазоне углов атаки а < 20° в верхней части ДУ САС устанавливается балансировочный груз (БГ), масса которого определяется
из условия устойчивости ОГБ на момент полной выработки топлива ДУ САС.
Анализ участка автономного движения ОГБ необходим для подтверждения безударности отделения ОГБ от РН, получения начальных условий отделения СА и обоснования массы БГ.
При расчете движения учитывали следующие факторы: разбросы массовых, центровочных и инерционных характеристик (МЦИХ) ОГБ, аэродинамических характеристик, различные типы аварийных ситуаций, а также разбросы параметров траектории и атмосферы.
Методы исследования. Основной метод исследования динамики отделения ОГБ от РН — это численное интегрирование уравнений пространственного движения абсолютно твердых тел (АТТ) и анализ их относительного движения [4—6].
В ходе исследования сравнивались результаты расчетов, выполненных для трех математических моделей ОГБ.
Модель 1: ОГБ—АТТ («тяжелый» ОГБ — МЦИХ соответствуют характеристикам ОГБ на момент отделения от РН (модель 1а); «легкий» ОГБ — МЦИХ соответствуют характеристикам ОГБ на момент выгорания топлива ЦРД и УРД (модель 1б)).
Z' 1 Головной обтекатель II _ РН
-J-
Рис. 1. Головной блок
Использование такой модели не дает точного представления о движении реального ОГБ, масса которого изменяется по мере выгорания топлива, но позволяет оценить диапазон значений кинематических параметров ОГБ на момент отделения СА.
Модель 2: ОГБ-тело переменной массы (МЦИХ изменяются по известному линейному закону).
В ходе расчета с использованием данной модели на каждом шаге интегрирования проводится пересчет МЦИХ ОГБ. При этом перемещение собственно центра масс головного блока в инерциальном пространстве не учитывается, что вносит искажение в картину относительного движения РН и ОГБ. Однако такой подход позволяет оценить вращательную составляющую движения ОГБ.
Модель 3: ОГБ представлено двумя телами (АТТ и телом переменной массы), которые соединены связью типа «заделка».
Данная модель лишена недостатков приведенных вариантов представления ОГБ, к тому же в перспективе может быть использована для реализации сквозного расчета отделения СА: отделение ОГБ от РН, автономный полет ОГБ и затем отделение СА от ОГБ.
Рассмотрим третий вариант моделирования ОГБ подробнее. Имеем два тела (рис. 2): тело А — АТТ и тело В — тело, МЦИХ которого изменяются по известному закону в зависимости от времени. Поскольку по условию эти тела представляют собой один рассматриваемый объект, то объединим тела А и В связью типа «заделка» в месте их механического соединения, предотвращая смещение и поворот тел относительно друг друга. В этой точке прикладываем к телам силы и моменты реакции связи.
R ^
z/c X
Рис. 2. Схема для определения реакций в точке связи
Движение тела А определяется численным интегрированием уравнений движения, а кинематические параметры тела В переменной массы на каждом шаге интегрирования — следующими зависимостями [7]:
®2 =®1; У2 = Vl +Й1 • С1С2,
где ®1, ®2 — угловые скорости вращения тел А и В относительно их центров масс; С1С2 — вектор, проведенный из центра масс тела А к центру масс тела В; ^ь — скорости центров масс тел А и В.
Для определения реакций в точке связи для каждого тела запишем уравнения движения в виде [8]
ч -г
m -—- = Fi; dt2
т dOi - / - ч -Ji—- + -(/i-Oi ) = LQ , dt
(1) (2)
где Г), Ьа — главные векторы внешних сил и моментов, I = 1, 2.
Внешние силы и моменты, действующие на ОГБ, приведем к центру масс тела А. Силу и момент реакции в точке связи обозначим как Я и М, получим:
Fi = X Fk + R;
Lei =X Lk + M;
(3)
(4)
3 = -Я; 1с2 = - м. Уравнения связи типа «заделка»:
% + Г1 = г02 + г2; 6 = 62,
(5)
(6)
(7)
(8)
где г01, г02 — радиус-векторы положения центров масс тел А и В относительно инерциальной системы отсчета СХХХ (см. рис. 2); г1, г2 — радиус-векторы положения точки связи соответственно в связанных с телами А и В системах координат С1 и С2Х2У"2Z2; 61, б2 — угловые ускорения вращения тел А и В относительно их центров масс.
Продифференцируем дважды уравнение (7), учитывая изменение положения центра масс тела В, и, используя выражения (1)-(6), выделим из уравнений связи (7) и (8) силу и момент реакции связи:
(т-1 + т-1) Я + [/-1М + /-1 (п х Я)] х п + [/-1М + /2-1 (г2 х Я)] х г2 =
= 62 х(Й2 X Г2 ) + 202 X Г2 - Ji 1
XLk -Oi x(Jic51)
X ri -mi iXFk-Oi x(Oi Xri);
J-1 ( x R)+/—i ( X R)+(/Г1+/—i) M = -Ji-i
XLk -Oi X(JiOi )
k
Решая указанную систему уравнений, получаем силу и момент реакции в точке связи, которые используем для интегрирования уравнений движения тела.
Корректность такого подхода к моделированию движения тела переменной массы может оцениваться сравнением результатов расчета по первой и третьей моделям ОГБ с фиксированными МЦИХ.
Диапазоны кинематических параметров ОГБ на момент отделения СА определяли методом статистических испытаний [9, 10].
Результаты моделирования. Выявлено, что при возникновении аварии на большей части атмосферного участка выведения (от старта до высот ~30 км и М < 4) автономный полет ОГБ характеризуется устойчивым положением в потоке с колебаниями по углам атаки не более 10°. Основные источники угловых возмущений — аварийный завал РН, приводящий к выключению ее ДУ и началу аварийной циклограммы, и программный разворот ОГБ — создают стартовые возмущения, которые в процессе движения до момента отделения СА демпфируются аэродинамическим потоком. На рис. 3 приведены графики изменения компонента угловой скорости (юх, юг, юг) вращения ОГБ и угла атаки
<ах, %
со у, %
Рис. 3. Графики компонента угловой скорости (а-в) вращения и пространственного угла атаки (г) при отделении ОГБ в случае аварии РН на 70-й с полета
а (г) во времени при отделении головного блока от РН на участке максимального скоростного напора. На каждом графике показаны результаты моделирования с использованием трех приведенных моделей ОГБ.
При этом на всем рассматриваемом участке траектории в случае аварии РН обеспечивается безударность отделения на начальном участке относительного движения ОГБ и РН (рис. 4), а также интенсивный увод ОГБ от носителя (рис. 5). На рисунках приведены результаты моделирования отделения ОГБ в случае аварии РН на участке максимального скоростного напора (70 с полета) для крайних случаев начальных возмущений РН (шу0 = ±10°/с, шz0 = ±10°/с ).
Рис. 4. Кинограммы относительного движения ОГБ в системе координат,
связанной с РН
Исследование отделения ОГБ в случае аварии РН на высотах более 30 км показало, что в связи с интенсивным падением плотности атмосферы для определенных сочетаний возмущающих факторов и разбросов параметров начальное возмущение, полученное ОГБ, не может парироваться аэродинамическим демпфированием (рис. 6). Таким образом, выбранная из условия статической устойчивости ОГБ в диапазоне углов атаки а < 20° масса БГ не гарантирует устойчивость движения головного блока при его отделении от РН на всем атмосферном участке выведения.
X, км X, км
Рис. 5. Положение центра масс ОГБ в связанной с центром масс РН системе координат (X0Z — а, Y0Z — б), оси которой параллельны осям инерциальной системы координат
Рис. 6. Графики компонента угловой скорости (а-в) вращения и пространственного угла атаки а (г) при отделении ОГБ в случае аварии на 110-й с полета: ----модель 1а; —о--модель 1б; .......— модель 2;--модель 3
Анализ моделирования позволяет сделать вывод о существенных отличиях результатов, полученных при исследовании движения ОГБ с использованием первой модели (ОГБ-АТТ). В то время как результаты, полученные с помощью моделей ОГБ 2 и 3, очень близки.
-100 -50 0 50 шх,°/с -50 -25 0 25 а у, % -50 -25 0 25 со2,7с 0 45 90 135 а,'
J
11 ■ ,_-гтгПтГ
4» =110 с
тИтут,
-ГП-ггЛ
-гтну
I ■ Ц 1|1
-100 -50 0 50 а>х, 7с -50 -25 0 25 а>¥, 7с -50 -25 0 25 со2,7с 0 45 90 135 а,'
= 114 с
ш
и I
-100 -50 0 50 <сх, 7о -50 -25 0 25 7с -50 -25 0 25 а2,7с 0 45 90 135 а, °
Рис. 7. Гистограммы распределения компонента угловой скорости и пространственного угла атаки ОГБ при его отделении
от аварийной РН
На основе модели 3 была разработана программа расчета и проведена серия статистических испытаний отделения ОГБ в случае аварии РН в конце рассматриваемого атмосферного участка. Анализ движения ОГБ на этом участке позволяет сделать вывод о том, что отделение СА от ОГБ возможно при произвольном положении изделий в потоке и их интенсивном вращении. На рис. 7 приведены гистограммы распределения компонента угловой скорости для наиболее неблагоприятного случая отделения от РН при потере управляемости.
Заключение. Подтверждена безударность отделения ОГБ от РН на всем рассматриваемом участке выведения.
Исследование автономного движения головного блока при аварии РН на большей части атмосферного участка выведения (на высотах не более ~30 км) подтвердило сохранение устойчивости движении ОГБ в потоке благодаря аэродинамическому демпфированию начальных возмущений, что обеспечивает благоприятные условия отделения СА. Выбор БГ из условия статической устойчивости ОГБ не гарантирует устойчивости в случае его отделения на участке с интенсивным падением скоростного напора. Так, при аварии РН на высотах более 30 км возможна потеря устойчивости головного блока в процессе его автономного движения, что определяет отделение СА при произвольной ориентации ОГБ в потоке.
Разработанная модель движения позволяет имитировать движение ОГБ с переменными МЦИХ и проводить сквозное моделирование движения, включая последующее отделение СА.
ЛИТЕРАТУРА
1. Колесников К.С., Кокушкин В.В., Борзых С.В., Панкова Н.В. Расчет и проектирование систем разделения ступеней ракет. М.: Изд-во МГТУ им. Н.Э. Баумана, 2006. 376 с.
2. Спасение экипажа при авариях в плотных слоях атмосферы // Космические аппараты: веб-сайт. URL: http://kmapp.narod.ru/st064.htm (дата обращения: 01.04.2016).
3. Средства аварийного спасения КК «Союз» // Космические аппараты: веб-сайт. URL: http://kmapp.narod.ru/st066.htm (дата обращения: 01.04.2016).
4. Борзых С.В., Кокушкин В.В., Петров Н.К., Яськов В.В. Разработка и моделирование процессов отделения крупногабаритных ракетно-космических блоков // Космическая техника и технологии. 2013. № 1. С. 44-55.
5. Виттенбург Й. Динамика систем твердых тел. М.: Мир, 1980. 292 с.
6. Кокушкин В.В., Хомяков М.К., Овсянникова Н.Ю. Динамика движения элементов головного блока пилотируемого корабля в аварийной ситуации // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2016. № 1. С. 4-16. DOI: 10.18698/0236-3941-2016-1-4-16
7. Лойцянский Л.Г., Лурье А.И. Курс теоретической механики: В 2 т. Т. I. Статика и кинематика. М.: Наука, 1982. 352 с.
8. Лойцянский Л.Г., Лурье А.И. Курс теоретической механики: В 2 т. Т. 2. Динамика. М.: Наука, 1983. 640 с.
9. Вентцель Е.С. Теория вероятностей. М.: Физматлит, 1962. 564 с.
10. Соболь И.М. Численные методы Монте-Карло. М.: Физматлит, 1973. 312 с.
Кокушкин Вячеслав Вячеславович — д-р техн. наук, профессор кафедры «Теоретическая механика» МГТУ им. Н.Э. Баумана (Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5).
Хомяков Михаил Кимович — канд. техн. наук, специалист РКК «Энергия» им. С.П. Королёва (Российская Федерация, 141070, Московская обл., г. Королёв, ул. Ленина, д. 4а).
Овсянникова Наталья Юрьевна — специалист РКК «Энергия» им. С.П. Королёва (Российская Федерация, 141070, Московская обл., г. Королёв, ул. Ленина, д. 4а).
Богомолов Николай Вячеславович — специалист РКК «Энергия» им. С.П. Королёва (Российская Федерация, 141070, Московская обл., г. Королёв, ул. Ленина, д. 4а).
Просьба ссылаться на эту статью следующим образом:
Кокушкин В.В., Хомяков М.К., Овсянникова Н.Ю., Богомолов Н.В. Динамика отделения головного блока от ракеты-носителя в случае аварии на атмосферном участке выведения // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2017. № 2. С. 4-14. БО!: 10.18698/0236-3941-2017-2-4-14
DYNAMIC OF EMERGENCY FAIRING SEPARATION FROM LAUNCH VEHICLE IN ATMOSPHERIC ASCENT
V.V. Kokushkin1 М.Ю Khomyakov2 N.Yu. Ovsyannikova2 N.V. Bogomolov2
[email protected] [email protected]
1 Bauman Moscow State Technical University, Moscow, Russian Federation
2 S.P. Korolev Rocket and Space Corporation Energia, Korolev, Moscow Region, Russian Federation
Abstract
The purpose of the article was to examine dynamics of the emergency separation of the fairing along with the spacecraft "Soyuz" reentry module during the rocket atmospheric ascent. The study suggests different mathematical models of variable mass body — separate fairing. First, we chose the optimal model allowing us to perform end-to-end simulation from fairing separation to reentry module self-contained motion. Then, we assessed non-collision of fairing separation from the launch vehicle, recorded the separation process and examined the fairing self-contained motion. Next, we found the regularities of the transient process and determined the fairing flow stability rating. Finally, we made a statistic analysis of the fairing self-contained motion
Keywords
Separate fairing, emergency separation, separation dynamics, relative motion, stability rating
REFERENCES
[1] Kolesnikov K.S., Kokushkin V.V., Borzykh S.V., Pankova N.V. Raschet i proektirovanie sistem razdeleniya stupeney raket [Calculation and designing of submissile compartmentation system]. Moscow, Bauman MSTU Publ., 2006. 376 p.
[2] Spasenie ekipazha pri avariyakh v plotnykh sloyakh atmosfery [Crew rescue in case of crack-up in dense atmosphere]. Kosmicheskie apparaty: website. Available at: http://kmapp.narod.ru/ st064.htm (accessed 01.04.2016) (in Russ.).
[3] Sredstva avariynogo spaseniya KK "Soyuz" [Emergency recovery equipment of "Soyuz" spacecraft]. Kosmicheskie apparaty: website. Available at: http://kmapp.narod.ru/st066.htm (accessed 01.04.2016) (in Russ.).
[4] Borzykh S.V., Kokushkin V.V., Petrov N.K., Yas'kov V.V. Development and simulation of large space rocket stages separation processes. Kosmicheskaya tekhnika i tekhnologii [Space Engineering and Technology], 2013, no. 1, pp. 44-55 (in Russ.).
[5] Vittenburg Y. Dinamika sistem tverdykh tel [Rigid body system dynamics]. Moscow, Mir Publ., 1980. 292 p.
[6] Kokushkin V.V., Khomyakov M.K., Ovsyannikova N.Yu. Dynamics of manned spaceship fairing elements motion in emergency. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Mashinostr. [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng.], 2016, no. 1, pp. 4-16 (in Russ.). DOI: 10.18698/0236-3941-2016-1-4-16
[7] Loytsyanskiy L.G., Lur'ye A.I. Kurs teoreticheskoy mekhaniki. T. I. Statika i kinematika [Engineering mechanics course. Vol. 1. Statics and kinematics]. Moscow, Nauka Publ., 1982. 352 p.
[8] Loytsyanskiy L.G., Lur'ye A.I. Kurs teoreticheskoy mekhaniki. T. 2. Dinamika [Engineering mechanics course. Vol. 1. Dynamics]. Moscow, Nauka Publ., 1983. 640 p.
[9] Venttsel' E.S. Teoriya veroyatnostey [Probability theory]. Moscow, Fizmatlit Publ., 1962. 564 p.
[10] Sobol' I.M. Chislennye metody Monte-Karlo [Monte Carlo numerical methods]. Moscow, Fizmatlit Publ., 1973. 312 p.
Kokushkin V.V. — Dr. Sc. (Eng.), Professor of Theoretical Mechanics Department, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation).
Khomyakov M.K. — Cand. Sc. (Eng.), specialist of S.P. Korolev Rocket and Space Corporation Energia (ul. Lenina 4a, Korolev, Moscow Region, 141070 Russian Federation).
Ovsyannikova N.Yu. — engineer of S.P. Korolev Rocket and Space Corporation Energia (ul. Lenina 4a, Korolev, Moscow Region, 141070 Russian Federation).
Bogomolov N.V. — engineer of S.P. Korolev Rocket and Space Corporation Energia (ul. Le-nina 4a, Korolev, Moscow Region, 141070 Russian Federation).
Please cite this article in English as:
Kokushkin V.V., Khomyakov M.K., Ovsyannikova N.Yu., Bogomolov N.V. Dynamic of Emergency Fairing Separation from Launch Vehicle in Atmospheric Ascent. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Mashinostr. [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng.], 2017, no. 2, pp. 4-14. DOI: 10.18698/0236-3941-2017-2-4-14
