Модель одноканальной системы ориентации космического аппарата.

Цель работы

Построить модель одноканальной реактивной системы ориентации космического аппарата (КА) в SIMULINK.

Задачи

  1. Построить модель одноканальной системы ориентации космического аппарата.
  2. Задать следующие начальные условия \(\varphi_0 = 20^\circ\), \(\dot{\varphi}_0 = 15^\circ/c\).
  3. Построить фазовый портрет (\(\varphi, \dot{\varphi}\)).
  4. Построить график изменения угла \(\varphi\) по времени.
  5. Построить график изменения угловой скорости \(\dot{\varphi}\) по времени.
  6. Построить график изменения момента, создаваемого двигателями, по времени.
  7. Построить фазовый портрет системы при \(k_\omega = 0\).
  8. Построить график изменения угла \(\varphi\) при \(k_\omega = 0\).
  9. Построить график изменения угловой скорости при \(k_\omega = 0\).

Реактивная система ориентации

В состав системы ориентации КА входят датчик угла, датчик угловой скорости и исполнительные органы. Датчик угла определяет угол между осью \(x_c\), связанной КА, и некоторой осью \(x_0\) фиксированной в пространстве (неподвижная система координат). Датчик угловой скорости определяет угловую скорость КА по отношению к неподвижной системе координат.

Исполнительными органами системы ориентации являются два реактивных двигателя, создающие постоянный момент относительно центра масс КА. Сила \(P_-\) создает отрицательный момент, направленный в сторону уменьшения угла \(\varphi\). Сила \(P_+\) создает положительный момент, направленный в сторону увеличения угла \(\varphi\). Далее будем считать, что \(P_+ = P_- = P\) и плечи этих сил относительно центра масс КА равны (рисунок 1). Система ориентации должна поддерживать угол и угловую скорость КА близкими к нулю.

Рисунок 1 - Работа реактивной системы ориентации

Сигналы, поступающие с датчика угла и угловой скорости, суммируются с заданными коэффициентами: \[s = k_\varphi \phi + k_\omega \dot{\varphi}\]

\(k_\varphi\) и \(k_\omega\) – коэффициенты усиления угла и угловой скорости.

Итоговый сигнал \(s\) используется для управления включением и выключением двигателей ориентации. Модель функции \(s\) в SIMULINK, выполненная в виде отдельной подсистемы показана на рисунке 2.

Рисунок 2 - Модель подсистемы “Функция s”

Когда значение функции (сигнала) \(s\) по модулю превышает некоторое пороговое значение \(s_{max}\) происходит включение двигателя. В процессе работы двигатель создает постоянный момент относительно центра масс КА. Знак этого момента зависит от знака величины \(u\): \[u(s) = \begin{cases} +1, & s \geq s_{max} \\ -1, & s \leq -s_{max} \\ 0, & |s| \leq s_{max} \\ \end{cases}\]

Если значение сигнала \(s\) меньше некоторого порогового значения, то функция \(u\) равна нулю, таким образом у рассматриваемой системы ориентации есть зона нечувствительности (\(\mid s\mid \leq s_{max}\)). График функции переключения \(u(s)\) показан на рисунке 3.

Рисунок 3 - Функция управления u(s)

Величина управляющего момента будет определяться выражением \[M = - u(s) \cdot P \cdot h.\]

Модель функции переключения в SIMULINK в виде подсистемы показана на рисунке 4.

Рисунок 4 - Подсистема “Функция u(s)”

Основой модели является блок “реле”. Настройки блока “Реле” приведены на рисунке 5. Выходной сигнал блока “Реле” равен значению, указанному в поле ввода “Output when on”, если входной сигнал, увеличиваясь, превышает значение “Switch on point”. При уменьшении входного сигнала до уровня ниже “Switch off point” выходной сигнал равен “Output when off”. На вход блока “реле” подается модуль функции u(s). Для учета знака функции u(s) определяется ее знак при помощи блока “Sign” (Знак): +1 или -1, далее это значение умножается на результат работы блока “реле”.

Рисунок 5 - Настройки блока “Реле”

Значение сигнала с выхода подсистемы “Функция u(s)”, который может принимать значения минус 1, 0 или +1 умножается на момент силы тяги. Результат этой операции (момент силы тяги) подается в подсистему, моделирующую динамику КА.

Схема подсистемы, описывающей модель КА, приведена на рисунке 6. На вход подается величина момента, создаваемого двигателями. Момент делится на момент инерции КА относительно центра масс, для этого используется блок “усилитель” (Gain) с коэффициентом усиления 1/J, где J – момент инерции КА. Далее сигнал дважды интегрируется: таким образом интегрируется дифференциальное уравнение движения (вращения) КА в плоскости: \[J \ddot{\varphi} = M(s).\]

Рисунок 6 - Подсистема “Модель КА”

Общий вид модели показан на рисунке 7.

Рисунок 7 - Модель реактивной системы ориентации КА

На рисунке 8 показан пример фазового портрета системы.

Рисунок 8 - Фазовый портрет

Васильев В.Н. Системы ориентации космических аппаратов / В.Н. Васильев. - М.: ФГУП «НПП ВНИИЭМ», 2009. - 310 с.


© 2020. All rights reserved.

Powered by Hydejack v9.1.4