Иван Шунейко - Пилотируемые полеты на Луну

10. Роwers W. F., Mc Dannell J. P. Switching conditions and a synthesis technique for the singular Saturn guidance problem. AIAA Paper № 70—965, ЭИ АиР, 1971, № 15; РЖ, 1971, 3.41.59

11. Luh J. Y. S., Maguiraga M. Minimum trajectory sensitivity of a large launch booster control system. IEEE Trans. Aerospace and Electron. System., 1969, 5, № 2, ЭИ AиР, 1969, № 39; РЖ, 1969, 11.41.200

12. Horn Н. J., Chandler D. C., Buckelew V. L. Iterative guidance applied to generalized missions. J. Spacecraft and Rockets, 1969, 6, № 1, ЭИ АнР, 1969, № 39: РЖ, 1969, 10.41.264

13. Мartin F. Н., Battin R. H. Computer—controlled steering of the Apollo' spacecraft. J. Spacecraft and Rockets, 1968, 5, № 4, (ЭН АиР, 1968, № 30); РЖ, 1968, 10.41.252

14. Mc Ruer D. Т., Weir D. Н., Klein R. Н. A pilot-vehicle systems approach to longitudinal flight director design. AIAA Paper № 70—1001, (ЭИ АиР, 1971, № 13)

15. Сhen P. P. Real—time Kalman filtering of Apollo LM/AGS rendezvous radar data. AIAA Paper, № 70—957, ЭИ АиР, 1971, № 10; РЖ 1971, 2.41.273

16. Satin A. L., Рixley P. T. Statistics of state—vector corrections for Apollo onboard computers. AIAA Paper, № 70—162, ЭИ АиР 1970 № 33; РЖ, 1970, 10.41.255

17. Bielkowicz P., Horrigan R. C., Walsh R. C., Manual onboard methods of orbit determination. AIAA Paper, № 70—159, (ЭИ АиР, 1970, № 33)

18. Salinger S. N., Brandstaller J. J. Application of recursive estimation and Kalman filtering to Doppler tracking. IEEE Trans. Aerospасе and Electron. Syst, 1970, 6, № 4, ЭИ АиР, 1970. № 45; РЖ, 1970, 12.41.220

19. Cox K. J. Apollo reaction control systems. IEEE Transection automatic control,IEEE №9C41-AC, Aug.4, 1969

20. Stubbs G.S., Penchuk A., Schlundt R.W Digital autopilot for thrust vector control of the Apollo CSM and CSM/LM vehicles. AIAA Paper № 69—847

21. Miller J. E., Laats Ain. Apollo guidance and control system flight experience. AIAA Paper № 69—891

22. Stengel R. F. Manual attitude control of the Lunar Module. AIAA Paper № 69—892

23. Mason W. L., Wedekind D. E. Prediction and measurement of strap-down inertial measurement unit performence on lunar missions. AIAA Paper № 70—1028

Задача прицеливания на траектории выведения к Луне состоит в определении параметров старта с Земли и участка разгона с околоземной орбиты (независимые переменные) для заданного набора параметров прицеливания (зависимые переменные). Параметрами прицеливания являются радиус периселения окололунной траектории Rm, ширина периселения в лунной системе координат Lm и высота условного перицентра траектории возвращения RE. В качестве трех независимых переменных рассматриваются время старта Tl, продолжительность движения на промежуточной околоземной орбите tc и удельная энергия на траектории к Луне С3. Эти переменные, будучи определенными с помощью итеративного процесса, устанавливают 3 важных зависимых параметра задачи: время старта для заданного азимута, время до второго включения ступени S-IVB при разгоне с околоземной орбиты (на втором или третьем обороте) и удвоенную удельную энергию эллиптической траектории полета к Луне.

При вычислении независимых переменных используется метод Ньютона-Рафсона для системы нелинейных уравнений. Линеаризованные уравнения, записанные в матричной форме, имеют следующий вид:

где [X] —вектор-столбец поправок ?Хj к независимым переменным; [В] – вектор-столбец ошибок зависимых переменных (Yi—Yi); .[А]—якобиан (матрица частных производных ошибок зависимых переменных по поправкам, dYi/dXi).

Для заданного азимута запуска траектория выведения на орбиту ИСЗ оптимизируется независимо от расположения Земли и Луны. Однако участок разгона с орбиты зависит от расположения Земли и Луны, которое определяет требования к изменению плоскости движения при втором запуске ступени S-IVB. Поэтому участок выведения на траекторию полета к Луне должен оптимизироваться совместно с определением независимых переменных. Схема, выбранная для вычислительной программы прицеливания ракеты-носителя на участке выведения к Луне, основана на аппроксимации по методу наименьших квадратов оптимальных параметров активного участка полета ступени S-IVB, выражаемых через параметры гиперповерхности. Это позволяет независимо оптимизировать выведение на траекторию полета к Луне в процессе итерационного вычисления зависимых переменных. Гиперповерхность, показанная на рис. 31.1, образована путем состыковки конических сечений для двух притягивающих центров.

Рис. 31.1. Гиперповерхность траекторий полета к Луне

Она представляет собой семейство конических сечений, которые начинаются у Земли и заканчиваются вблизи сферы действия Луны. Гиперповерхность определяют следующие параметры: вектор цели Т°, удвоенная удельная энергия С3, угол между вектором цели и радиусом-вектором перигея ?.

Рис. 31.2. Параметры участка выведения на траекторию полета к Луне: а – вид в плоскости выведения; б – проекция на плоскость, перпендикулярную плоскости выведения

Параметры гиперповерхности используются в качестве независимых переменных полиномов, описывающих активный участок ступени S-IVB. С помощью этих полиномов определяются параметры участка выведения к Луне и вектора состояния. На рис. 31.2 показаны участок выведения и геометрические соотношения для определения гиперповерхности.

При использовании полиномов необходимо знать удельную энергию на траектории к Луне С3, угол между вектором цели и радиусом-вектором перигея ?, а также склонение вектора цели относительно плоскости промежуточной орбиты ?. С помощью полиномов вычисляются параметры: ? – угол между радиусом-вектором точки начала выведения и проекцией вектора цели на плоскость промежуточной орбиты (Т'); ? – угол между радиусом-вектором точки начала выведения и узлом орбиты (?); ? —истинная аномалия радиуса-вектора точки конца выведения; Rp – радиус перигея участка выведения к Луне. Параметры ?, ?, ? и Rp задают вектор состояния участка выведения.

В системе уравнений указанные величины используются для определения параметров активного участка и вычисления переменных, соответствующих моменту выключения двигателя. По полиномам также вычисляется приращение характеристикой скорости ?V при повторном включении ступени S-IVB. Величина этого приращения необходима для определения веса аппарата и времени работы двигателя. Из-за неточности учета влияния сжатия Земли и модели изменения тяги расчеты по полиномам не совпадают с результатами оптимизации активного участка методом вариационного исчисления на основе ожидаемых параметров отлета от Земли. Поэтому полиномы тарируются, чтобы обеспечить данные, точно совпадающие с результатами расчета активного участка методом вариационного исчисления. Постоянные поправочные члены, необходимые для тарировки, вычисляются как разница между результатами расчета методом вариационного исчисления и величинами, полученными путем оценки полиномов в первом приближении при С3=С3g (где C3g – приближенное значение), ?=0 (компланарный случай) и ?=6°:

где индекс «вар» относится к результатам расчета методом вариационного исчисления.

В компланарном случае ?=?+? и ??=??. Указанные поправочные члены получены для обеих возможностей запуска и должны использоваться всякий раз, когда параметры активного участка вычисляются по аппроксимирующим полиномам. Полиномы для участка выведения к Луне тарируются путем добавления вычисленных поправочных членов к приближенным величинам, полученным при подстановке в полиномы текущих величин С3, ? и ?.

Полином, определяющий величину ?V при втором включении ступени S-IVB, не тарируется. Однако при каждом расчете по полиному вычисляется разница между компланарным значением ?V, основанным на параметрах первого приближения (C3g, ?=0, ?=6°), и значением ?V, определяемым текущими величинами указанных параметров:

При планировании задачи полета на Луну определенное преимущество достигается в случае двух возможностей отлета с околоземной орбиты. Вторая возможность появляется приблизительно через 90 мин после первой (т. е. через один оборот на промежуточной орбите ИСЗ) и оказывается полезной в тех случаях, когда не все системы ракеты-носителя и космического корабля проверены и готовы к повторному включению двигателя для выведения на траекторию полета к Луне. В процессе подготовки полета принимается решение о том, сохранять ли время перелета к Луне для второй возможности таким же, какое требовалось для первой (класс 1) или уменьшить время полета для второй возможности на 90 мин (класс 2). Уменьшение времени перелета на 90 мин при использовании второй возможности позволяет сохранить время прибытия к Луне приблизительно таким же, как для первой возможности.

Рис. 31.3. Выбор компромиссного времени старта: 1 – первая возможность; 2 – вторая возможность; 3 – компромиссное время старта; 4 – моменты времени, соответствующие компланарному выведению

В процессе прицеливания ракеты-носителя и выбора времени запуска исследуются оба класса перелетов с целью получения максимального веса на траектории к Луне для обеих возможностей запуска. Рис. 31.3 иллюстрирует логику выбора времени запуска для двух рассматриваемых классов.