CC BY
100
УДК 623.54:623.451.08
МЕТОДИКА ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ ПРОДОЛЬНОГО ДВИЖЕНИЯ КОСМИЧЕСКОГО АППАРАТА ПРИ ПОСАДКЕ НА ПОВЕРХНОСТЬ МАЛОГО НЕБЕСНОГО ТЕЛА
В.П. Казаковцев, Чень Даньхэ, В.В. Корянов
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация
e-mail: [email protected]; [email protected]; [email protected]
Разработана методика определения параметров движения космического аппарата при его посадке на поверхность малого небесного тела при учете влияния реакции от грунта поверхности на подвижные опоры и работы прижимных двигателей. Рассмотрено применение методики для исследования динамики движения космического аппарата на примере его посадки на поверхность Луны. Для случая вертикальной посадки выбрано оптимальное значение силы в амортизаторах, позволяющей при умеренных силах реакции грунта на опоры получать рациональные значения углов подкосов после посадки. Проведена оценка влияния начального углового положения космического аппарата к поверхности Луны на изменение положения опор, нормальных сил реакции грунта на опоры, деформации амортизаторов при посадке. Рассчитаны графики продольной и вертикальной скоростей космического аппарата с учетом влияния подвижности опор при их неодновременном касании грунта.
Ключевые слова: посадка на малые небесные тела, метод динамики движения космического аппарата, оценка влияния возмущающих факторов, условия для безопасной посадки.
METHOD FOR DETERMINING THE LONGITUDINAL MOTION PARAMETERS OF THE SPACECRAFT DURING LANDING ON THE SURFACE OF A SMALL CELESTIAL BODY
V.P. Kazakovtsev, Chen Danhe, V.V. Koryanov
Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: [email protected]; [email protected]; [email protected]
The method is developed for determining the motion parameters of a spacecraft during its landing on the surface of a small celestial body, taking into account the influence of the response from the soil on the movable supports and the effect of clamping engines. The method application for studying dynamics of the spacecraft motion is considered by the example of its landing on the Moon surface. For the case of vertical landing, the optimal value of the force in shock absorbers is selected, which makes it possible to obtain rational values of angles of struts after the landing with moderate forces of soil response on supports. The influence of the initial angular position of the spacecraft with respect on the Moon surface on a change in the support position, normal forces of the soil response on supports, deformation of shock absorbers during the landing is estimated. Plots of longitudinal and vertical spacecraft speeds are calculated taking into account the influence of movability of supports when they touch the soil non-simultaneously.
Keywords: landing on small celestial bodies, the method for study of dynamics of spacecraft movement, evaluation of influence of disturbing factors, conditions for safe landing.
Программа исследования различных небесных тел солнечной системы предполагает в ближайшие десятилетия осуществление поле-
26 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №1
тов космических аппаратов (КА) к ним. В ряде случаев для получения информации о небесном теле предполагается посадка КА на его поверхность.
Особый интерес представляют малые небесные тела, у которых отсутствует атмосфера и имеется достаточное поле притяжения. К таким телам относятся астероиды типа Церера, Паллада и другие спутники планет (Луна у Земли, Ио и Европа у Юпитера и др.).
На этапе эскизного проектирования КА, предназначенного для посадки на поверхность малого небесного тела, требуется методика для определения параметров движения КА, оценки силовой нагрузки на посадочные опоры. Это позволяет определить основные характеристики посадочных опор, выбрать параметры амортизаторов, оценить условия безопасной посадки на поверхность малого небесного тела.
В опубликованных ранее работах [1-6] показано, что при посадке на поверхность Луны наиболее опасным случаем является продольное движение КА с боковой составляющей скорости при посадке на одну опору на поднимающийся склон поверхности. Поэтому на ранних этапах проектирования возможно использование методики расчета параметров посадки при продольном движении КА. Это позволяет оценить наибольшие нагрузки на посадочные опоры КА.
Сформулируем математическую модель движения КА на заключительном этапе посадки на поверхность малого небесного тела.
Основные допущения:
а) рассматриваем продольное движение КА с тремя опорами;
б) аэродинамические силы отсутствуют;
в) ускорение притяжения является постоянным;
г) трение в шарнирах опор не учитываем;
д) считаем жестким участок поверхности небесного тела, на который производится посадка КА.
При формировании математической модели движения КА используются следующие системы координат:
1) поверхностная система координат (ПСК) X, Y, Z связана с поверхностью небесного тела (оси X и Z лежат в плоскости поверхности посадки, при этом ось X лежит в плоскости первой опоры);
2) связанная с КА система координат (ССК) XcYcZc направлена вдоль продольной оси КА, ось Xc — перпендикулярна оси Yc и находится в плоскости XY ПСК, ось Zc — дополняет СК до правой);
3) гравитационная система координат (ГСК) XgYgZg (ось Yg направлена вдоль линии действия силы притяжения, ось Xg перпендикулярна оси Yg и лежит в плоскости XY ПСК).
На рис. 1 показан вид сверху на КА. Цифрами 1, 2, 3 обозначены точки касания опор поверхности небесного тела; точками Pi,P2,P3
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 1 27
обозначены места расположения прижимных двигательных установок (ДУ). Углы 72,73 определяют разворот второй и третьей опор относительно оси OcXc; точки B и D являются местами крепления подкосов.
Размеры всех опор одинаковые. На рис. 2 приведены основные размеры первой опоры: l0(AE) — проекция подкосов AB, AD на плос-
Рис.1. Вид Сверху на КА с тремя кость XcYc; laM(АС) — длина стерж-опорами
ня с амортизатором; lE — расстояние от центра масс (ц.м.) КА до точки Е по направлению оси OcXc; h0 ,hi — расстояния от ц.м. КА до точек крепления подкоса и амортизатора соответственно; H — расстояние от ц.м. КА до точки А по направлению оси OcXc; а0 — угол наклона подкоса к линии параллельной оси OcXc [7].
Основные силы, действующие на КА при посадке показаны на рис. 3.
Pi ,P2 ,P3 — значение тяги ДУ; mg — сила тяготения; FNi, FN2, FN3 — проекции реакции грунта по направлению оси OY ПСК; FTi, FT2X, FT3X — силы трения трех опор по направлению оси OX ПСК.
Также обозначено: $ — угол наклона КА к оси OX; 0g — угол между гравитационной вертикалью и осью OY ПСК.
Дифференциальные уравнения, описывающие продольное движение КА при посадке на поверхность небесного тела, в проекциях на оси поверхностной СК сформулируем в следующем виде [8]:
Рис. 2. Основные размеры крепления стержня и амортизатора к КА
28 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №1
Рис. 3. Схема действия сил на КА при посадке
dVx
dt
dVy
dt
dx
dt
1
m
1
m
[Fti + Ft2x + Ftsx — (Pi + P2 + P3) sin F] —g sin вд; [FN1 + FN2 + FN3 + (P1 + P2 + p3) cos F] — g cos 6g;
Vx;
dy
dt
dwz
dt
Vy;
1 г
— (Ft 1 + FT2X + FT3X )y + Fn 1 x1—
z
— (Fn2 + Fn3)x2 + Pixp — P2xp cos 72 — P3xp cos 73 ;
dF
dt
w
z
(1)
Размеры опор h0, h1, H, lE, l0 — выбираются заранее. Начальные значения следующих величин определяются по формулам
а0 — arcsin
H — ho lo
1ам0 — у/(l0 cos «0)2 + (H + h1)2;
в0 — arcsin I
l0 cos a0
l
ам0
(2)
Рассмотрим алгоритм решения задачи посадки КА на поверхность небесного тела для первой (второй) опоры, направленной по оси OX в ПСК.
Запишем условие касания первой (второй) опорой поверхности в виде
8у1 = У — У1 — о, (3)
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 1 29
где
yi(2) = ho cos $ - Ie sin $ + lo(ao - $1(2)); (4)
$2 = arcsin(cos(n — y2) sin $) (5)
— угол наклона КА для второй опоры).
Для Sy1(2) = 0 определим текущие значения для у и $, полученных при интегрировании уравнений (1):
У + lE sin $ — h0 cos $
a1(2)
= arcsin
lo
+ $1(2);
Hi(2) = lo sin a1(2) + Ie ;
L1(2) = lo cos a1(2) + ^e ;
^ам1(2) = \J(H + h1(2))2 + (L1(2) — lE)2;
01(2) = 1ам1(2) — 1ам1(2) — изменение длины амортизатора;
x1 = lo cos(a1 — $) + lE cos $ + ho sin $;
X2 = cos 72(lo cos(«2 — $2) + Ie cos $2 + ho sin $2);
e • (L1 — Ie \
01(2) = arcsrnl —------- .
\ 1ям1(2) /
В соответствии с рис. 4 по 51 определяем значение силы в амортизаторе:
0 01(2) ^ 0jk;
(7)
(6)
17 /с\ ) т-t 01(2) — 0&1 г г е
FL1(2) (0) = \ Fo-------------J-------, j < 01(2) < °fc1(2);
0o
Fo, 01(2) > 0k1(2),
где 0o — начальное значение упругой зоны деформации амортизатора; 0k1j — значение упругой зоны деформации амортизатора; j — начальное значение деформации при j-м касании опоры поверхности, Fo — сила при разрушении амортизатора.
Рис. 4. График изменения силы в амортизаторе
30 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. № 1
Запишем уравнения равновесия сил в точке А первой (второй) опоры как по оси OX, так и по оси OY. После преобразований получим формулу для определения реакции грунта на опору:
F1(2) = FL1(2) COs(ai(2) + ^1(2))/COs(a1(2) + ^1(2) — A), (8)
где ^1(2) — коэффициент трения на первой (второй) опоре.
Реакция грунта на первую (вторую) опору по оси OY
FN 1(2) = F1(2) c°S ^1(2). (9)
Модуль силы трения на второй опоре
Ft2 = F2 sin ^2.
Для определения силы трения на первой (второй) опоре вычислим скорость ее перемещения по направлению оси OX:
Кг1(2) = Vx + ^Vx1(2), (10)
где VX — скорость движения ц.м. КА по направлению оси OX; ^VX1(2) = — Vy tg(a1(2) — А) — скорость перемещения опоры при изменении длины амортизатора; SVx2 = —SV2 cos y2 и 5Vz2 = —SV2 sin y2 — проекция и скорости перемещения второй опоры на оси OX и OZ.
Сила трения на первой и второй опорах определяется как
Ft1 = Ft1x = —F1 sin(^) ж1 . (11)
abs(Vx 1)
Третья опора расположена симметрично второй опоре, поэтому в продольном движении при посадке КА на поверхность небесного тела силы от реакции грунта на вторую и третью опоры одинаковые:
FN 2 = FN 3, FTX2 = FTX3. (12)
Рассмотрим пример применения предложенной методики по определению параметров продольного движения КА при посадке на поверхность Луны [9, 10].
Конструктивные параметры КА и опор: m = 900 кг, Iz = = 800 кг • м2, h0 = 0,38 м, h1 = 0,20 м, H = 1,06 м, l0 = 0,8 м, lE = 1,0 м, 80 = 0,001 м , F0 = 3000 Н.
Оценим влияние различных сил в амортизаторах при вертикальной посадке КА при следующих условиях (y0 = 1,4 м, Vy0 = 0,1, Vx0 = 0, ^1 = ^2 = 0,2, А0 = 0°, 90 = 0°). На рис. 5. представлены графики изменения нормальных сил реакции грунта на опоры КА при F0 = 4000 Н (1-й вариант), F0 = 3000 Н (2-й вариант), F0 = 2000 Н (3-й вариант).
Выполнив анализ этих графиков, выявили, что при первом варианте сила реакции грунта на опоры имеет очень большое значение и
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 1 31
Рис. 5. Графики изменения нормальных сил реакции грунта на три опоры КА в точках (1, 2, 3) касания поверхности небесного тела при различных силах в амортизаторе
достаточно большой угол наклона подкоса (ак = 47,8°) после посадки КА. В третьем варианте сила реакции грунта на опоры небольшая, но маленький угол наклона подкоса (ак = 32,3°) после посадки КА. Это неприемлемо, так как предельный угол наклона подкоса после посадки КА не должен превышать акпред = 30,0°. При учете влияния возмущающих факторов для третьего варианта силы в амортизаторах угол наклона подкоса после посадки КА может быть меньше допустимого значения. Поэтому для опор рассматриваемых размеров и конструктивных величин КА более приемлемым является второй вариант, при котором угол наклона подкоса после вертикальной посадки КА ак = 43,4°.
Оценим различные параметры посадки КА при следующих начальных условиях: уо = 1,4 м, Vy0 = 0,1, Vx0 = 0, ^ = у,2 = 0,2, = 10°,
вд = 0°. На рис. 6 представлены графики нормальных сил реакции грунта на опоры КА.
Видно, что первоначально грунта касаются одновременно второй и третьей опорами. Затем КА разворачивается на этих опорах до момента касания грунта первой опорой. После касания грунта первой опорой нормальная сила реакции грунта на вторую и третью опоры дважды уменьшается до нулевой.
Рис. 6. Графики изменения нормальных сил FN1, FN2 реакции грунта на первую (•) и вторую (■) опоры
32 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. № 1
Рис. 7. Графики изменения деформации D\ и D2 амортизаторов первой и второй опор при посадке КА
Рис. 8. Графики скоростей Vx и Vy КА при посадке
На рис. 7 приведены графики деформации амортизаторов при посадке КА, из которых следует, что деформация первой опоры наибольшая. Угол наклона подкоса на первой опоре после посадки КА
«к = 33,4°.
На рис. 8 представлены графики скоростей Vx и Vy при посадке КА. Видно, что при касании грунта второй и третьей опорами начинается рост Vx, обусловленный изменением длины их амортизаторов, и уменьшение по модулю Vy. По окончании посадки КА эти скорости равны нулю.
Отметим, что после окончания посадки угол наклона КА к поверхности Луны $к = -2,87°.
Заключение. 1. Разработана методика определения параметров движения космического аппарата при его посадке на поверхность малого небесного тела при учете влияния реакции от грунта поверхности на подвижные опоры и работы прижимных двигателей.
2. На примере посадки космического аппарата на поверхность Луны проведена оценка влияния некоторых возмущающих факторов на изменение положения опор, сил в амортизаторах; рассмотрены условия безопасной посадки: для случая вертикальной посадки КА выбрано значение силы в амортизаторах, позволяющее при умеренных силах реакции грунта на опоры получать рациональные углы подкосов после посадки; проанализировано влияние начального угла наклона КА к
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 1 33
поверхности Луны на динамику его движения при посадке; выяснены
особенности работы опор, изменения горизонтальной и вертикальной
скоростей КА, угловое положение КА после его посадки.
ЛИТЕРАТУРА
1. Проектирование спускаемых автоматических космических аппаратов / под ред.
B. М. Ковтуненко. М.: Машиностроение, 1985. 264 с.
2. Баженов В.И., Осин М.И.Посадка космических аппаратов на планеты. М.: Машиностроение, 1978. 158 с.
3. Егоров В.А. Пространственная задача достижения Луны. М.: Наука, 1965. 224 с.
4. Виноградов А.П. Передвижная лаборатория на Луне ЛУНОХОД-1. Том 1. М.: Наука, 1971. 128 с.
5. Гродзовский Г.Л., Иванов Ю.Н., Токарев В.В. Механика космического полета. Методы оптимизации. М.: Наука, 1975. 702 c.
6. Базилевский А.Т., Григорьев Е.И., Ермаков С.Н.Проектирование спускаемых автоматических космических аппаратов. М.: Машиностроение, 1985. 264 с.
7. Борзых С.В., Воронин В.В. Динамика взаимодействия трансформируемого механического посадочного устройства возвращаемого космического аппарата с посадочной поверхностью // Вестник МГТУ им. Н.Э. Баумана. Машиностроение. Спец. вып. “Крупногабаритные трансформируемые космические конструкции и материалы для перспективных ракетно-космических систем”. 2012. C. 69-79.
8. Корянов В.В. Исследование динамики движения спускаемого аппарата при жесткой посадке на поверхность планеты // Полет. Общеросс. науч.-технич. ж. 2010. № 1. C. 42-49.
9. Пичхадзе К.М., Воронцов В.А., Защиринский А.М., Пономарев В.А. Система спуска с орбиты и аварийного спасения на основе надувного тормозного устройства // Полет. Общеросс. науч.-технич. ж. 2003. № 8. C. 9-13.
10. Казаковцев В.П., Корянов В.В. Метод исследования динамики углового движения космического спускаемого аппарата с надувным тормозным устройством // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2012. № 3 (88).
C. 39-46.
REFERENCES
[1] Kovtunenko V.M. Proektirovanie spuskaemykh avtomaticheskikh kosmicheskikh apparatov [Constructional design of landing unmanned space vehicles. Moscow, Mashinostronie Publ., 1985. 264 p.
[2] Bazhenov V.I., Osin M.S. Posadka kosmicheskikh apparatov na planet [Spacecraft landing on the planety]. Moscow, Mashinostronie Publ., 1978, 158 p.
[3] Egorov V.A. Prostranstvennaya zadacha dostizheniya Luny. [Spatial problem of reaching the Moon]. Moscow, Nauka Publ., 1965. 224 p.
[4] Vinogradov A.P. Peredvizhnaya laboratoriya na Lune LUNOKHOD-1. Tom 1 [The mobile laboratory lunar vehicle LUNOKHOD-1. Vol. 1]. Moscow, Nauka Publ., 1971. 128 p.
[5] Grodzovskiy G.L., Ivanov Yu.N., Tokarev V.V. Mekhanika kosmicheskogo poleta. Metody optimizatsii [The space-flight mechanics. Optimization methods]. Moscow, Nauka Publ., 1975. 702 p.
[6] Bazilevskiy A.T., Grigor’ev E.I., Ermakov S.N. Proektirovanie spuskaemykh avtomaticheskikh kosmicheskikh apparatov [Constructional design of landing unmanned space vehicles].Moscow, Mashinostronie Publ., 1985. 264 p.
34 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2014. №1
[7] Interaction dynamics of reshaped mechanical landing gear of re-entering body with the landing surface. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Mashinostr., Spetsvyp. “Krupnogabaritnye transformiruemye kosmicheskie konstruktsii i materialy dlya perspektivnykh raketno-kosmicheskikh sistem” [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng., Spec. Issue “Large-sized reshaped space structures and materials for advanced rocket-and-space systems”], 2011, pp. 69-79 (in Russ.).
[8] Koryanov V.V. Dynamic study of lander moving with an impact landing on the planet surface. Obshcherossiyskiy nauchno-tekhnicheskiy zhurnal “Polet” [All-Russ. Sc.-Techn. Journ. “Flight”)], 2010, no. 1, pp. 42-49 (in Russ.).
[9] Pichkhadze K.M., Vorontsov V.A., Zashchirinskiy A.M., Ponomarev V.A. Deorbit system and emergency rescue based on the inflatable braking device. Obshcherossiyskiy nauchno-tekhnicheskiy zhurnal “Polet” [All-Russ. Sc.-Techn. Journ. “Flight”)], 2003, no. 8, pp. 9-13 (in Russ.).
[10] Kazakovtsev V.P., Koryanov V.V. Research technique of the attitude dynamics of a landing unmanned space vehicle with an inflatable braking device. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Mashinostr. [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng.], 2012, no.3 (88), pp. 39-46.
Статья поступила в редакцию 22.05.2013
Виктор Поликарпович Казаковцев — д-р техн. наук, профессор кафедры “Динамика и управление полетом ракет и космических аппаратов” МГТУ им. Н.Э. Баумана. Автор более 120 научных работ в области баллистики и динамики полета космических и спускаемых аппаратов.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
V.P. Kazakovtsev — Dr. Sci. (Eng.), professor of “Dynamics and Flight Control of Rockets and Spacecrafts” department of the Bauman Moscow State Technical University. Author of more than 120 publications in the field of ballistics and flight dynamics of space and descent vehicles.
Bauman Moscow State Technical University, Vtoraya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Чень Даньхэ — аспирантка кафедры “Динамика и управление полетом ракет и космических аппаратов” МГТУ им. Н.Э. Баумана. Специализируется в области моделирования динамики движения спускаемых аппаратов.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Chen Danhe — post-graduate of “Dynamics and Flight Control of Rockets and Spacecrafts” department of the Bauman Moscow State Technical University. Specializes in the field of simulation of dynamics of descent vehicle movement.
Bauman Moscow State Technical University, Vtoraya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Всеволод Владимирович Корянов — канд. техн. наук, доцент кафедры “Динамика и управление полетом ракет и космических аппаратов” МГТУ им. Н.Э. Баумана. Автор более 10 научных работ в области моделирования баллистики и динамики движения космических и спускаемых аппаратов.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
V.V. Koryanov — Cand. Sci. (Eng.), assoc. professor of “Dynamics and Flight Control of Rockets and Spacecrafts” department of the Bauman Moscow State Technical University. Author of more than 10 publications in the field of simulation of ballistics and dynamics of motion of space and descent vehicles.
Bauman Moscow State Technical University, Vtoraya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2014. № 1 35
