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

?=(2/T)tg^-1u (22.10)

В результате получаются графики D*(j?) G*r(j?), которые модифицируются с помощью функций M(j?) и Н(j?) для определения разомкнутой частотной характеристики аппарата как твердого тела

G0r(j?)=D*(j?) M(j?) G*r(j?)Н (j?) (22.11)

Контур компенсации эксцентриситета вектора тяги

Влияние контура компенсации эксцентриситета вектора тяги на характеристики разомкнутой цепи ЦАП по конструктивным соображениям ограничивается областью ниже 2 рад/сек. Это значительно ниже частоты квантования интегратора этого контура, равной 2 гц (12,56 рад/сек), и частоты квантования ЦАП (25 гц для основного блока и 12,5 гц для корабля Apollo). Можно показать, что влияние указанных частот квантования пренебрежимо мало в диапазоне частот до 2 рад/сек, и контур компенсации эксцентриситета вектора тяги может быть аппроксимирован передаточной функцией непрерывного сигнала вида

где Км – коэффициент усиления контура компенсации эксцентриситета вектора тяги;

Тм – постоянная времени низкочастотного фильтра.

Контур управления траекторией полета

Динамика контура управления траекторией полета зависит от времени до окончания работы ЖРД, tgo. Для больших значений tgo этой зависимостью можно пренебречь при выводе передаточных функций контура управления траекторией полета. Эти функции даже в приближенном виде весьма полезны для понимания влияния контура управления траекторией полета на работу ЦАП в целом.

Контур управления траекторией полета выполняет следующие операции.

1. Приращения скорости, измеряемые акселерометром, накапливаются и вычисляется текущая скорость V;

2. Каждые 2 сек вычисляется разность Vg между требуемой скоростью Vr и текущей скоростью V.

3. Каждые 2 сек находится векторное произведение Vg и ?V, где ?V – изменение скорости в течение последних 2 сек.

4. Результат векторного произведения нормируется по отношению Vg и ?V и затем умножается на коэффициент усиления Ksteer для получения вектора команды скорости ориентации.

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

6. Управляющие команды по угловой скорости аппарата умножаются на период квантования ЦАП Т, чтобы получить-приращения, которые подаются в ЦАП через каждые Т сек.

Приближенные аналитические соотношения, описывающие-динамику контура управления траекторией полета выводятся при следующих допущениях:

1) все тригонометрические функции заменяются их приближенными значениями для малых углов;

2) аппарат рассматривается как твердое тело с инерци-альной измерительной платформой, установленной в центре тяжести;

3) эффекты квантования пренебрежимо малы;

4) ось X аппарата и ось ЖРД первоначально выставлены параллельно вектору скорости Vr; в этом случае угол между вектором тяги и вектором Vr определяется как ? – ? (? – ориентация аппарата относительно вектора Vr, ? – отклонение ЖРД от начального направления; считается положительным, если создает положительное ускорение ?);

5) время tgo постоянно;

6) вектор Vr постоянен по величине и направлению;

7) запаздывание при вычислении в контуре управления траекторий полета пренебрежимо мало.

При сделанных допущениях угол между вектором тяги и вектором требуемой скорости определяется как ? – ?, а соответствующие углы в плоскости тангажа и в плоскости рыскания для векторов Vg и ?V аппроксимируются выражениями

(Ts – период квантования контура управления траекторией полета 2 сек).

Нормированное векторное произведение векторов Vg и ?V дает разность ?vg – ??V, которая после умножения на коэффициент усиления Kst дает сигнал по угловой скорости ориентации .

Эта команда по угловой скорости поступает в ЦАП в виде приращения ?cT, находится разность между командным и измеренным приращением для определения ошибки ориентации.

Таким образом, ?c образуется путем квантования с периодом 2 сек непрерывного интеграла от ?—? и последовательного преобразования в приращения, которые суммируются на интервале в Т сек. Такая комбинация двух периодов квантования и аналогичного интегрирования усложняет задачу учета влияния контура управления траекторией полета на устойчивость ЦАП.

Для упрощения анализа влияния управления траекторией полета автопилот можно рассматривать как элемент с непрерывным сигналом в контуре управления траекторией полета. Тогда остается только одна частота квантования 0,5 гц, связанная с формированием управляющего сигнала ?c. Частотную характеристику разомкнутого контура управления траекторией полета можно представить ?-преобразованием

где F1(?)-преобразование

F2(?) представляет процесс получения угла векторного произведения ?vg – ??V из интеграла (?—?); F1(?) описывает умножение этого угла на коэффициент усиления Kst для получения ?c, последующее интегрирование ?c для получения ?c', образования разности (?—?) автопилотом и ее интегрирование.

Взаимное влияние автопилота и процесса управления траекторией полета проявляется более наглядно с помощью приближенного аналитического метода, заключающегося в добавлении к частотной характеристике разомкнутого контура ЦАП влияния управления траекторией полета.

При этом предполагается, что эффекты транспонирования частот отсутствуют и частотная характеристика F1(?) может быть заменена частотной характеристикой в области реальных частот

а также имеет место

где Кr – коэффициент эффективности управления системы ЖРД-аппарат.

В предположении, что

F2(?) примет вид

Частотная характеристика разомкнутого контура управления траекторией полета может быть представлена в виде

где

В области низких частот, где применима функция Gst(j?), влияние управления траекторией на характеристики разомкнутого контура ЦАП можно аппроксимировать путем прибавления Gst(j?) к единичному коэффициенту обратной связи по углу ориентации

Следует заметить, что выражение H(j?) является хорошей аппроксимацией только для низких значений ?, меньше 0,2 рад/сек. Однако именно этот диапазон частот представляет наибольший интерес при анализе эффектов управления траекторией полета.