Пользовательская переменная масса 6DOF ветер (кватернион)

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

Библиотека

Уравнения Motion/6DOF

Описание

Блок Custom Variable Mass 6DOF Wind (Quaternion) рассматривает вращение зафиксированного ветром координатного кадра (_Xw, _Yw, _Zw) о плоском Наземном кадре ссылки (_Xe, _Ye, _Ze). Источник зафиксированного ветром координатного кадра является центром тяжести тела, и тело принято, чтобы быть твердым, предположение, которое избавляет от необходимости рассматривать силы, действующие между отдельными элементами массы. Плоский Наземный кадр ссылки рассматривается инерционным, превосходное приближение, которое позволяет силам из-за движения Земли относительно “фиксированных звезд” быть пропущенными.

Переводное движение зафиксированного ветром координатного кадра приведено ниже, где приложенные силы [_Fx, _Fy, _Fz] ^Tнаходятся в зафиксированном ветром кадре. Vre _w является относительной скоростью в осях ветра в который массовый поток ( $\dot{m}$ ) извлекается или добавляется к телу.

$\begin{array}{l} {\bar{F}}_{w} = [\begin{array}{l} F_{x} \\ F_{y} \\ F_{z} \end{array}] = m ({\dot{\bar{V}}}_{w} + {\bar{ω}}_{w} \times {\bar{V}}_{w}) + \dot{m} \bar{V} r e_{w} \\ A_{b e} = D C M_{w b} \frac{[{\bar{F}}_{w} - \dot{m} V_{r e}]}{m} \\ {\bar{V}}_{w} = [\begin{array}{l} V \\ 00 \end{array}], {\bar{ω}}_{w} = [\begin{array}{l} p_{w} \\ q_{w} \\ r_{w} \end{array}] = D M C_{w b} [\begin{matrix} p_{b} - \dot{β} \sin α \\ q_{b} - \dot{α} \\ r_{b} + \dot{β} потому что α \end{matrix}], {\bar{w}}_{b} = [\begin{array}{l} p_{b} \\ q_{b} \\ r_{b} \end{array}] \\ A_{b b} = D C M_{w b} [\frac{\bar{F} w - \dot{m} V_{r e}}{m} - {\bar{ω}}_{w} \times {\bar{V}}_{w}] \end{array}$

Вращательные движущие силы зафиксированного телом кадра приведены ниже, где прикладные моменты [L M N] ^T и тензор инерции, I относительно источника O. Тензор инерции I намного легче задать в зафиксированном телом кадре.

$\begin{array}{l} {\bar{M}}_{b} = [\begin{array}{l} L \\ M \\ N \end{array}] = I {\dot{\bar{ω}}}_{b} + {\bar{ω}}_{b} \times (I {\bar{ω}}_{b}) + \dot{I} {\bar{ω}}_{b} \\ A_{b b} = [\begin{array}{l} {\dot{U}}_{b} \\ {\dot{V}}_{b} \\ {\dot{W}}_{b} \end{array}] = D C M_{w b} [\frac{\bar{F} w - \dot{m} V_{r e}}{m} - {\bar{ω}}_{w} \times {\bar{V}}_{w}] \\ I = [\begin{matrix} I_{x x} & - I_{x y} & - I_{x z} \\ - I_{y x} & I_{y y} & - I_{y z} \\ - I_{z x} & - I_{z y} & I_{z z} \end{matrix}] \end{array}$

Интегрирование скорости изменения вектора кватерниона приведено ниже.

$[\begin{array}{l} {\dot{q}}_{0} \\ {\dot{q}}_{1} \\ {\dot{q}}_{2} \\ {\dot{q}}_{3} \end{array}] = - \frac{}{12} [\begin{matrix} 0 & p & q & r \\ - p & 0 & - r & q \\ - q & r & 0 & - p \\ - r & - q & p & 0 \end{matrix}] [\begin{array}{l} q_{0} \\ q_{1} \\ q_{2} \\ q_{3} \end{array}]$

Параметры

Основной

Units

Задает модули ввода и вывода:

Модули	Силы	Момент	Ускорение	Скорость	Положение	Масса	Инерция
`Metric (MKS)`	Ньютон	Ньютон-метр	Метры в секунду придали квадратную форму	Метры в секунду	Метры	Килограмм	Килограммометр придал квадратную форму
`English (Velocity in ft/s)`	Фунт	Фунт ноги	Ноги в секунду придали квадратную форму	Ноги в секунду	Футы	Краткий заголовок	Отложите нога придала квадратную форму
`English (Velocity in kts)`	Фунт	Фунт ноги	Ноги в секунду придали квадратную форму	Узлы	Футы	Краткий заголовок	Отложите нога придала квадратную форму

Mass Type

Выберите тип массы, чтобы использовать:

`Fixed`	Масса является постоянной в течение симуляции.
`Simple Variable`	Масса и инерция отличаются линейно как функция массового уровня.
`Custom Variable`	Масса и изменения инерции настраиваемы.

Выбор Custom Variable соответствует ранее описанным уравнениям движения.

Representation

Выберите представление использованию:

`Wind Angles`	Используйте углы ветра в рамках уравнений движения.
`Quaternion`	Используйте кватернионы в рамках уравнений движения.

Выбор Quaternion соответствует ранее описанным уравнениям движения.

Initial position in inertial axes

Трехэлементный вектор для начального местоположения тела в плоском Наземном кадре ссылки.

Initial airspeed, sideslip angle, and angle of attack

Трехэлементный вектор, содержащий начальную скорость полета, начальный угол заноса и начальный угол нападения.

Initial wind orientation

Трехэлементный вектор, содержащий начальные углы ветра [банк, курс полета, и направляющийся], в радианах.

Initial body rotation rates

Трехэлементный вектор для начальной буквы зафиксированные телом угловые уровни, в радианах в секунду.

Include mass flow relative velocity

Установите этот флажок, чтобы добавить массовый скоростной порт родственника потока. Это - относительная скорость, в которой масса аккумулируется или удаляется.

Include inertial acceleration

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

Атрибуты состояния

Присвойте уникальное имя каждому состоянию. Используйте имена состояния вместо путей к блоку в течение процесса линеаризации.

Чтобы присвоить имя к одному состоянию, введите уникальное имя между кавычками, например, 'velocity'.
Чтобы присвоить имена к нескольким состояниям, введите разграниченный запятой список, окруженный фигурными скобками, например, {'a', 'b', 'c'}. Каждое имя должно быть уникальным.
Если параметр пуст (' '), никакое присвоение имени не происходит.
Имена состояния применяются только к выбранному блоку с параметром имени.
Количество состояний должно разделиться равномерно среди количества имен состояния.
Можно задать меньше имен, чем состояния, но вы не можете задать больше имен, чем состояния.
Например, можно задать два имени в системе с четырьмя состояниями. Имя применяется к первым двум состояниям и второму имени к последним двум состояниям.
Чтобы присвоить имена состояния с переменной в рабочей области MATLAB^®, введите переменную без кавычек. Переменная может быть вектором символов, массивом ячеек или структурой.

Position: e.g., {'Xe', 'Ye', 'Ze'}

Задайте имя положения для состояния.