Моделирование процессов управления ЛА. Управление самолетом «Боинг - 747» при посадке.

Общие теоретические сведения.

Математическая модель САУ обычно записывается в форме векторного дифференциального уравнения, где свойства системы неявно выражены коэффициентами некоторых матриц. Эта форма записи является традиционной и предназначена для применения математического аппарата ТАУ. Получение матриц, характеризующих объекты реального мира, является задачей по сложности не уступающей, а по материальным затратам иногда превосходящей исследование уже готовой математической модели. Поэтому, перед тем как начать исследование математической модели реального летательного аппарата (ЛА) приведем общее описание пути построения ее самой.

Движение ЛА происходит под действием аэродинамических сил и моментов, сил и моментов двигателей и гравитационных сил. При определении движения ЛА необходимо в общем случае решить следующие задачи:

найти угловые и линейные скорости ЛА, обусловленные действием на него перечисленных сил и моментов;

определить углы ориентации ЛА относительно набегающего на него потока и осей координат , связанных с Землей;

определить перемещение самолета относительно Земли.

Обычно считается, что ЛА представляет собой абсолютно жесткое тело, масса и моменты инерции которого на рассматриваемых интервалах времени неизменны. Конфигурация ЛА имеет плоскость симметрии и массы распределены симметрично по отношению к этой плоскости. Используются следующие прямоугольные системы координат(рис.)

Нормальная земная система координат OXgYgZg. Эта система осей координат имеет неизменную ориентацию относительно Земли.

Начало координат совпадает с центром масс ЛА. Оси OXg и ОZg лежат в горизонтальной плоскости, ось OYg направлена вверх по местной вертикали.

Связанная система координат OXYZ. Начало координат расположено в центре масс ЛА. Ось ОХ лежит в плоскости симметрии и направлена вдоль линии хорд крыла к носовой части ЛА. Ось ОY лежит в плоскости симметрии самолета и направлена вверх (при горизонтальном полете), ось OZ дополняет систему до правой.

Положение связанной системы осей координат OXYZ относительно нормальной земной системы координат OXgYgZg может быть полностью определено так называемыми углами Эйлера (рис.2):

· угол рыскания y,

· угол тангажа J,

· угол крена g.

Последовательно поворачивая связанную систему на каждый из углов Эйлера, можно прийти к любому угловому положению связанной системы относительно осей нормальной системы. Положение связанной системы относительно вектора скорости центра масс ЛА определяется углами :

-атаки a,

-скольжения b.

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

1) Уравнения движения центра масс (уравнения поступательного, линейного движения; уравнения сил):

clip_image002 clip_image004

clip_image006 clip_image008

где clip_image010 - вектор скорости вращения связанных с ЛА осей;

clip_image012 - вектор гравитационных сил;

clip_image014 - сила тяги двигателей;

clip_image016 - главный вектор аэродинамических сил.

Символом clip_image018 обозначена полная производная (производная в земной системе), символом clip_image020 - местная производная (в связанной системе), clip_image022- любой вектор.

Тогда уравнения перемещения в проекциях на связанные оси запишутся так:

clip_image024

clip_image026

clip_image028

2) Уравнения вращения вокруг центра масс (уравнения вращательного, углового движения; уравнения моментов):

clip_image030 clip_image032

clip_image034- для плоской симметрии,

где clip_image036- момент количества движения;

clip_image038 - главный момент внешних сил;

Т - тензор инерции ЛА относительно связанных осей;

Jx,y,z - моменты инерции.

Уравнение вращения в проекциях на связанные оси будут иметь вид:

clip_image040

clip_image042

clip_image044

clip_image046

3) Закон сложения угловых скоростей:

clip_image048

Отсюда:

clip_image050

и

clip_image052

4) Закон сложения линейных скоростей:

clip_image054

clip_image056

clip_image058 clip_image060 clip_image062- плотность воздуха при Н=0

5) Проекции скорости на связанные оси:

clip_image064

6) Проекции сил на связанные оси.

Проекции силы тяжести:

clip_image066

Проекции силы тяги двигателя:

clip_image068

clip_image070

clip_image072

Проекции аэродинамической силы:

clip_image074

clip_image076

clip_image078

7) Проекция момента сил на связанные оси:

clip_image080

clip_image082

clip_image084

8) Безразмерные коэффициенты аэродинамических сил и моментов:

сx,y,z или mx,y,z = f(clip_image086),

clip_image088- число Маха,

clip_image090 - число Рейнольдса.

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

Уравнения продольного движения:

clip_image092

Уравнения бокового движения:

clip_image094

clip_image096

Требуется рассчитать переходный процесс для модели бокового движения тяжелого самолета «Боинг - 747» в посадочной конфигурации. Его линеаризованная модель движения в конфигурации , используемой при посадке , без применения системы автоматического повышения устойчивости имеет вид :

clip_image098clip_image098[1]clip_image100 ( 2.1)

Здесь приняты следующие обозначения для фазовых переменных :

u - скорость бокового скольжения ;

wY - угловая скорость рыскания ;

wX - угловая скорость крена ;

g - угол крена ;

y - угол рыскания ;

Z - боковое отклонение .

Управлениями в системе являются :

dЭ - угол отклонения элеронов ;

dН - угол отклонения руля направления .

Возмущением в системе является ветер :

v - скорость ветра .

Уравнения записаны в безразмерных величинах. В качестве единицы для углов принято 0.01 рад, а для скоростей - 0.305 м/с.

Для удобства построения характеристик, которые будут рассматриваться далее, приведем уравнения , описывающие данную модель , в матричной форме :

clip_image102 (2.2)

где

clip_image104 (2.3)

clip_image106

Вектор фазовых переменных обозначим через x ;

clip_image108.

Управление организовано в виде обратной связи

clip_image110 , (2.4)

Где clip_image112

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

Предлагаю ознакомиться с аналогичными статьями: