Далее...
Алгоритм расчета ориентации
Структура алгоритма расчета ориентации была выбрана на основе калмановской теории фильтрации. Принцип работы алгоритма состоит в сравнении экстраполированных и измеренных значений угла.
Для расчета ориентации аппарата необходимыми ЦАП исходными измерениями являются углы кардана инерциальной стабилизированной платформы, которые выдаются каждые 0,1 сек. Чтобы выделить любые смещения углового ускорения, являющиеся следствием тяги посадочного ЖРД от углового ускорения, создаваемого ЖРД РСУ, в блок расчета ориентации необходимо ввести дополнительную информацию. Вводя информацию о работе ЖРД РСУ и информацию об управлении вектором тяги посадочного ЖРД, можно осуществить фильтрацию основных сигналов управления ориентацией, не прибегая к сложным цепям расчета угловой скорости и смещения углового ускорения.
Изменения угловых скоростей аппарата от действия ЖРД РСУ на интервале последнего управления рассчитываются по формулам
(24.9)
Смещение углового ускорения, происходящее вследствие команд на отклонение вектора тяги на интервале последнего управления, подсчитывается по уравнению
(24.10)
где T=0,1 сек период квантования; uq, ur—(+ 1,0—1) командные сигналы на поворот вектора тяги относительно осей Q и R.
Углы кардана инерциальной стабилизированной платформы, полученные в процессе предыдущего выполнения расчетов, хранятся в блоке памяти бортовой ЭЦВМ. Измерив углы кардана в данный момент, можно рассчитать изменение углов и перевести в изменение углов аппарата.
Разность между измеренным изменением ориентации и предсказанным изменением называется «необъяснимым» изменением ориентации. «Необъяснимое» изменение вычисляется и прибавляется к суммарному предыдущему «необъяснимому» изменению, образуя «необъяснимое» изменение ориентации в данный момент
СМ, В ССЫЛКЕ
где qp, qq, qr—компоненты «необъяснимого» изменения ориентации по осям Р, Q, R.
Выражения в скобках в правой части уравнений (24.11) можно определить как компоненты предсказанного изменения ориентации. Заметим, что предсказанное изменение ориентации исключается вследствие команд на привод кардана ЖРД и точное выражение для изменения ориентации в результате работы ЖРД РСУ аппроксимируется более простым выражением.
«Необъяснимое» изменение ориентации используется для уточнения расчетов угловой скорости и углового ускорения аппарата. Но вначале логика используется, чтобы исключить измеренный квантованный шум. Вследствие того, что вероятное распределение этого шума не гауссовское, а прямоугольное, в каждой из осей (Р, Q, R) блока расчета ориентации шум может быть исключен нелинейным логическим фильтром.
Для каждой из осей Р, Q, R, если компонент «необъяснимого» изменения ориентации меньше, чем пороговая величина Qmax=0,14 град, в этом случае коррекции смещения sw и sa для расчетов угловой скорости и углового ускорения считаются нулевыми.
«Необъяснимое» изменение ориентации не обнуляется и когда «необъяснимое» изменение ориентации превышает пороговую величину, вводится коррекция. Нулевые коррекции sw и sa рассчитываются по уравнениям
(24.12)
После введения коррекции «необъяснимое» изменение ориентации для этой оси вновь устанавливается равным нулю.
Расчет угловой скорости и смещения углового ускорения теперь производятся по уравнениям
Заметим, что смещение углового ускорения вокруг оси Р принимается равным нулю и на последнем интервале управления не учитывается изменение угловой скорости вследствие отклонения кардана ЖРД.
На пассивных участках траектории полета аппарата считается, что компоненты смещения углового ускорения по осям Q и R равны нулю.
Динамические характеристики блока расчета ориентации сильно зависят от выбора коэффициентов усиления цепей фильтра Кw и Кa. Выбор этих коэффициентов основывается на компромиссе между быстротой расчета и исключением колебаний из-за плескания топлива.
Учитываемые алгоритмом расчета ориентации динамические эффекты: плескание топлива, изгибные колебания, силовое взаимодействие выхлопных струй ЖРД РСУ с конструкцией корабля, запаздывание тяги ЖРД РСУ по отношению к командам, эксцентриситет тяги ЖРД РСУ, ускорение выхлопных струй ЖРД РСУ при перемещениях по У и Z, не обнаруженные неисправности ЖРД РСУ, запаздывание сервомотора кардана ЖРД, разброс моментов инерции, упругие деформации сервопривода, точность модели измерения количества топлива в баках.
Законы управления ЖРД реактивной системы управления
Быстрейшая скорость повторения расчетов для ЦАП, возможная в пределах расчетной нагрузки бортовой ЭЦВМ лунного корабля, составляет 10 цикл/сек. Однако для легкой взлетной ступени 2 ЖРД РСУ могут создать угловое ускорение 50 град/сек2.
Простейшим законом управления мог бы быть циклический закон, ЖРД РСУ включено — выключено. Но в этом случае точность управления угловой скоростью взлетной ступени была бы только 5 град/сек.
Для осуществления управления с необходимой точностью с помощью фактора эффективности управления точно определяется продолжительность включения ЖРД РСУ, потребная для необходимого изменения угловой скорости.
Порядок расчета времени работы ЖРД РСУ для случая автоматического управления одним лунным кораблем, запрограммирован для одной оси как функция:
1) ошибки ориентации Qе и ошибки угловой скорости wе
2) параметров, рассчитываемых в цепи настройки, которые определяют кривизну и положение парабол фазовой плоскости;
3) оценки выгоды получения с одного или двух ЖРД требуемого импульса момента.
Обращаясь к этой программе и вводя на входе соответствующую информацию, определяют время работы ЖРД РСУ для каждой из осей.
Рис. 24.5. Закон управления ЖРД реактивной системы
лунного корабля на пассивных участках траектории полета
Логика управления в пассивном полете иллюстрируется диаграммой в фазовой плоскости рис. 24.5. Фазовая плоскость выше оси Qе делится на 5 зон, ограниченных параболами.
В логике управления в фазовой плоскости Qе—wе используется также угловое ускорение, складывающееся из углового ускорения от ЖРД РСУ и возмущающего углового ускорения. Крутизна парабол, ограничивающих зону 2 на фазовой плоскости, соответствует условию, когда на траектории работают ЖРД РСУ. Значения ускорений от ЖРД РСУ, которые определяют эту параболу, подсчитываются в контуре настройки. Парабола, разделяющая зоны 4 и 5, не траекторная. Это пологая кривая, поддерживающая управление с малой угловой скоростью, однако не настолько, чтобы существовали большие Qе.
Форма кривой сохранена параболической, чтобы логика границ была общей и удобной для компактности кодирования. Выбранная пологость соответствует угловому ускорению 1,4 град/сек2.
Пересечение парабол с осью Qе является функцией зоны
нечувствительности Qdb, выбираемой астронавтом или программой, выполняющей автоматический маневр. Допустимые ошибки ориентации области зоны нечувствительности могут быть 0,3, 1 и 5°.
Фазовая плоскость сконструирована таким образом, что управляющее действие из любых начальных условии с помощью двух главных импульсов приводит к заданной ориентации. На рис. 24.8 это показано сегментом AВС.
После перехода в зону нечувствительности состояние аппарата определяется минимальным предельным импульсным циклом DEFG с единичным зажиганием одного ЖРД РСУ каждый раз при пересечении зоны 3.
Рис. 24.6. Закон управления ЖРД реактивной системы
лунного корабля на активных участках траектории полета
На активных участках траектории полета используется логика управления, показанная на рис. 24.6. Точки пересечения парабол с осью Qe передвигаются в. зависимости от величины расчетного смещающего углового ускорения. Точки пересечения (1°, —2°, 0,75°), показанные на рисунке, типичны для активного участка взлета, когда вектор тяги ЖРД взлетной ступени смещен. Крутизна четырех парабол устанавливается по четырем разностям угловых ускорений. Верхняя левая парабола определяется минимальным ускорением amin= 1,4 град/сек2, как и в случае пассивного полета. Верхняя правая парабола определяется располагаемой разностью ускорений между ускорением от смещающего момента и противоположным по знаку ускорением от ЖРД РСУ.
Нижняя правая парабола определяется только смещающим угловым ускорением (ЖРД РСУ не работают).
Фазовая плоскость спроектирована так, что аппарат управляется низкочастотным предельным циклом, в котором один раз за цикл дается команда на зажигание ЖРД РСУ и отрабатывается единичный управляющий импульс, противоположный по знаку смещающему моменту.
Хотя ЦАП должен выполнять одни и те же функции управления аппаратом, когда лунный корабль состыкован с основным блоком и если он летит один, однако для состыкованной компоновки режим ЦАП разрабатывался отдельно.
Особенность проблемы заключалась в том, что в состыкованном виде корабль Apollo имеет 3 формы низкочастотных изгибных колебаний и слабый по прочности переходник, стыкующий основной блок и лунный корабль. Поэтому была разработана специальная логика торможения, исключающая, возможность зажигания ЖРД РСУ с частотой собственных изгибных колебаний корабля Apollo.
Кроме законов автоматического управления ориентацией аппарата с помощью ЖРД РСУ были разработаны специальные законы, обеспечивающие ручное управление. По законам ручного управления в зависимости от отклонения рукояток управления ориентацией вырабатываются сигналы на вход. в ЦАП.
При доводке характеристик ЦАП с использованием ортогональных осей координат Р, U, V применительно к законам управления ЖРД РСУ в фазовой плоскости относительно каждой отдельной оси было обнаружено перекрестное влияние каналов управления.
Распределение массы лунного корабля таково, что главные оси моментов инерции проходят вблизи осей Q и R, а не U и V. В результате, момент от ЖРД оси V порождает угловое ускорение не только относительно оси V, но и относительно оси U. Величина одновременно возникающего перекрестного ускорения такова, что в худшем случае вектор углового ускорения отклоняется на 15° от действующего вектора момента.
Чтобы устранить перекрестное влияние каналов управления ЖРД РСУ, была введена неортогональная система осей координат U' и V' (рис. 24.7).
Направление осей U' и V' определялось единственным требованием, чтобы направление U' было ортогонально угловому ускорению, возникающему от вектора момента оси V или вектора момента оси Р, направление V было ортогонально угловому ускорению, возникающему от вектора момента оси U или момента оси Р. Управления, определяющие угол s имеют вид
Находя вектор ошибки ориентации и вектор ошибки угловой скорости на осях U' и V' и используя компоненты U' и V' для определения требуемых векторов моментов ЖРД РСУ по осям U и V, исключается перекрестное влияние каналов управления.
Рис. 24.7. Система неортогональных осей координат лунного корабля
Закон управления направлением вектора тяги
При проектировании управления карданом ЖРД для изменения направления вектора тяги посадочной ступени лунного корабля предусматривалось использование этого управления для совмещения вектора тяги с центром масс аппарата и уменьшения таким образом расхода топлива на ЖРД РСУ. Так как управление ориентацией должно обеспечиваться ЖРД РСУ, при проектировании не требовалось задавать большую угловую скорость изменения направления вектора тяги, и был выбран маломощный и легкий привод кардана, обеспечивающий изменение угла отклонения ЖРД со скоростью 0,2 град/сек. Привод связан с ЦАП простым принципом включено-выключено. Для обеих осей Q и R, вокруг которых можно поворачивать вектор тяги, ЦАП может давать команды на угловую скорость 0,2 град/сек; —0,2 град/сек или ноль.
Однако минимизация расхода топлива на ЖРД РСУ была основной проблемой, и так как в процессе торможения и посадки ЖРД посадочной ступени работает все время, искали закон управления ориентацией с использованием посадочного ЖРД и без включения ЖРД РСУ по каналам U и V. Выбранному закону соответствует минимальное время управления.
Дифференциальное уравнение, связывающее сигнал, управляющий карданом ЖРД, с отклонением лунного корабля от требуемой ориентации относительно осей Q и R, имеет вид
(24.17)
Первая и вторая производные от ошибки ориентации по времени есть ошибка угловой скорости и ошибка углового ускорения. Предполагая, что все переменные состояния, используемые законом управления, измеряются без шума и без ошибок, оптимальное управление можно определить как функцию состояния системы в данный момент следующим образом
Параметр С имеет размерность — время и обращает переменные состояния (ae, we, Qe) в безразмерные переменные (x1, X2, X3). Оптимальный управляющий сигнал uoptдается в безразмерных величинах состояния системы.
Работа цифрового автопилота при первой посадке на Луну
В процессе первой посадки на Луну ЦАП в начале управлял лунным кораблем в автоматическом режиме и в конце по командам от ручки управления ориентацией; при этом ориентация вектора тяги ЖРД и лунного корабля изменилась от горизонтальной в начале торможения до вертикальной при посадке.
Первые 4 мин активного участка торможения после начального неустановившегося режима закон управления ориентацией вектора тяги работал успешно, медленно изменяя ориентацию без помощи ЖРД РСУ по каналам U и V. Затем из-за плескания топлива в баках возникли колебания большой амплитуды, выходящей за пределы зоны нечувствительности закона управления ЖРД РСУ. ЦАП вырабатывал команды управления ЖРД РСУ, ограничивавшие амплитуду колебаний лунного корабля.
Колебания угловой скорости тангажа с частотой 0,5 гц из-за плескания топлива видны на рис. 24.8.
Рис. 24.8. Угловая скорость тангажа на активном
участке траектории посадки лунного корабля Apollo-11
СМ. В ССЫЛКЕ
Автоматическое управление вело лунный корабль на посадку в кратер размерами с футбольное поле с большим количеством огромных камней. Н. Армстронг переключил ЦАП на ручное управление, изменил курс корабля, перелетел кратер, выбрал ровное место, и посадил корабль с помощью ручного управления.
После взлета с Луны управление с помощью ЦАП тангажом взлетной ступени показано на рис. 24.9. После старта и вертикального подъема в течение 10 сек была подана команда на автоматическое управление тантажом со скоростью 10 град/сек и выход на угол тангажа 52°.
После окончания маневра по тангажу наблюдался типичный низкочастотный предельный цикл изменения ориентации.
Когда было израсходовано все топливо взлетной ступени, ее центр масс переместился из заднего в переднее положение. При взлете центр масс был сзади вектора тяги, после выгорания топлива он переместился в положение впереди вектора тяги. Это явилось причиной изменения компонентов смещения углового ускорения вдоль обеих осей U' и V'. Компонент V' смещения углового ускорения, подсчитанный за последние 150 сек активного участка взлета, показан на рис. 24.10.
ЦАП настраивал свои параметры в пределах закона управления ЖРД РСУ и в соответствии с изменяющимся ускорением смещения.
Рис. 24.9 Команды управления по тангажу и
фактическое изменение тангажа при взлете с Луны
Apollo-11.
Рис. 24.10. Подсчитанный компонент V' углового
ускорения в конце активного участка траекто-
рии взлета с Луны Apollo-11.
Рис. 24.11. Число включений ЖРД оси V реак-
тивной системы управления на активном участке
траектории взлета Apollo-11
Общее число зажиганий ЖРД РСУ оси V в течение всего активного участка взлета показано на рис. 24. 11.
После взлета во время маневра по тангажу управление ориентацией вокруг оси V поддерживалось около 200 сек зажиганием только + V ЖРД.
Позднее амплитуда предельного цикла увеличилась и для поддержания ориентации стали включаться +V ЖРД и —V ЖРД. Хотя —V момент создает угловое ускорение в том же направлении, в котором действует ускорение смещения от ЖРД взлетной ступени, однако +V момент и —V момент создаются ЖРД РСУ, дающими тягу вверх (+Х) и поэтому не происходит потери AV. Возникновение моментов + V и —V объясняется колебанием взлетной ступени вследствие плескания топлива в баках.
Перед концом активного участка взлета после перемены знака смещения углового ускорения не подавалось команд на зажигание + V ЖРД и управление ориентацией поддерживалось вокруг оси V только включением —V ЖРД. [19.]
Ник