Спроектируйте демпфер отклонения от курса для струйного транспорта

Обзор этого Тематического исследования

Это тематическое исследование демонстрирует инструменты для классической системы управления путем продвижения через проект демпфера отклонения от курса для 747® реактивный транспортный самолет.

Создание струйной модели

Струйная модель во время рейса круиза в МАХЕ = 0.8 и H = 40 000 футов

A=[-.0558 -.9968 .0802 .0415;
     .598 -.115 -.0318 0;
    -3.05 .388 -.4650 0;
        0 0.0805 1 0];

B=[ .00729  0;
   -0.475   0.00775;
    0.153   0.143;
     0      0];

C=[0 1 0 0;
   0 0 0 1];

D=[0 0;
   0 0];

sys = ss(A,B,C,D);

Следующие команды задают эту модель в пространстве состояний как объект LTI и присоединяют имена к состояниям, входным параметрам и выходным параметрам.

states = {'beta' 'yaw' 'roll' 'phi'};
inputs = {'rudder' 'aileron'};
outputs = {'yaw' 'bank angle'};

sys = ss(A,B,C,D,'statename',states,...
                 'inputname',inputs,...
                 'outputname',outputs);

Можно отобразить модель LTI sys путем ввода sys. Эта команда приводит к следующему результату.

a = 
            beta      yaw     roll      phi
   beta  -0.0558  -0.9968   0.0802   0.0415
   yaw     0.598   -0.115  -0.0318        0
   roll    -3.05    0.388   -0.465        0
   phi         0   0.0805        1        0
 
b = 
          rudder  aileron
   beta  0.00729        0
   yaw    -0.475  0.00775
   roll    0.153    0.143
   phi         0        0
 
c = 
               beta   yaw  roll   phi
   yaw            0     1     0     0
   bank angle     0     0     0     1
 
d = 
                rudder  aileron
   yaw               0        0
   bank angle        0        0
 
Continuous-time model.

Модель имеет два входных параметров и два выходных параметров. Модули являются радианами для beta (угол заноса) и phi (угол банка) и радианы/секунда для yaw (уровень отклонения от курса) и roll (прокрутите уровень). Руководящий принцип и отклонения элерона исчисляются в радианах также.

Вычисление полюсов разомкнутого цикла

Вычислите полюса разомкнутого цикла и постройте их в s - плоскость.

>> damp(sys)
                                                                       
         Pole              Damping       Frequency      Time Constant  
                                       (rad/seconds)      (seconds)    
                                                                       
 -7.28e-03                 1.00e+00       7.28e-03         1.37e+02    
 -5.63e-01                 1.00e+00       5.63e-01         1.78e+00    
 -3.29e-02 + 9.47e-01i     3.48e-02       9.47e-01         3.04e+01    
 -3.29e-02 - 9.47e-01i     3.48e-02       9.47e-01         3.04e+01    

pzmap(sys)

Эта модель имеет одну пару слегка ослабленных полюсов. Они соответствуют тому, что называется "режимом Голландского шага".

Предположим, что вы хотите спроектировать компенсатор, который увеличивает затухание этих полюсов, так, чтобы получившиеся комплексные полюса имели коэффициент затухания ζ> 0.35 с собственной частотой ωn <1 рад/секунда. Можно сделать это использование аналитических инструментов Control System Toolbox™.

Анализ разомкнутого цикла

Во-первых, выполните некоторый анализ разомкнутого цикла, чтобы определить возможные стратегии управления. Запустите с ответа времени (вы могли использовать step или impulse здесь).

impulse(sys)

Импульсная характеристика подтверждает, что система слегка ослабляется. Но период времени слишком долог, потому что пассажиры и пилот более обеспокоены поведением в течение первых нескольких секунд, а не первых нескольких минут. Затем посмотрите на ответ за меньший период времени 20 секунд.

impulse(sys,20)

Посмотрите на график от элерона (введите 2) к углу банка (выводит 2). Чтобы показать только этот график, щелкните правой кнопкой и выберите I/O Selector, затем нажмите (2,2) запись. Селектор ввода-вывода должен выглядеть так.

Новую фигуру показывают ниже.

Самолет колеблется вокруг ненулевого угла банка. Таким образом самолет поворачивается в ответ на импульс элерона. Это поведение окажется важным позже в этом тематическом исследовании.

Как правило, демпферы отклонения от курса спроектированы с помощью уровня отклонения от курса, как обнаружено выведено и руководящего принципа как вход управления. Посмотрите на соответствующую частотную характеристику.

sys11=sys('yaw','rudder') % Select I/O pair.
bode(sys11)

От этой диаграммы Боде вы видите, что руководящий принцип оказывает значительное влияние вокруг слегка ослабленного режима Голландского шага (то есть, рядом ω = 1 рад/секунда).

Проект корневого годографа

Разумная цель проекта состоит в том, чтобы обеспечить порцию затухания ζ> 0.35 с собственной частотой ωn <1,0 рад/секунда. Поскольку самый простой компенсатор является статическим усилением, сначала попытайтесь определить соответствующие значения усиления с помощью метода корневого годографа.

% Plot the root locus for the rudder to yaw channel
rlocus(sys11)

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

rlocus(-sys11)
sgrid

Это выглядит лучше. При помощи простой обратной связи можно достигнуть коэффициента затухания ζ> 0.45. Нажмите на синюю кривую и переместите маркер данных, чтобы отследить значения затухания и усиление. Чтобы достигнуть 0,45 коэффициентов затухания, усиление должно быть приблизительно 2,85. Этот рисунок показывает маркер данных с подобными значениями.

Затем закройте обратную связь SISO.

K = 2.85;
cl11 = feedback(sys11,-K);  % Note: feedback assumes negative 
                            % feedback by default

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

impulse(sys11,'b--',cl11,'r',20)

Ответ с обратной связью обосновывается быстро и не колеблется очень, особенно когда по сравнению с ответом разомкнутого цикла.

Теперь замкните круг на полной модели MIMO и смотрите, как ответ от элерона смотрит. Обратная связь включает вход 1 и вывела 1 из объекта (используйте feedback с векторами индекса, выбирающими эту пару ввода/вывода). В подсказке MATLAB® ввести

cloop = feedback(sys,-K,1,1);
damp(cloop)   % closed-loop poles

         Pole              Damping       Frequency      Time Constant  
                                       (rad/seconds)      (seconds)    
                                                                       
 -3.42e-01                 1.00e+00       3.42e-01         2.92e+00    
 -2.97e-01 + 6.06e-01i     4.40e-01       6.75e-01         3.36e+00    
 -2.97e-01 - 6.06e-01i     4.40e-01       6.75e-01         3.36e+00    
 -1.05e+00                 1.00e+00       1.05e+00         9.50e-01    

Постройте импульсную характеристику MIMO.

impulse(sys,'b--',cloop,'r',20)

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

Создание фильтра размытия

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

H(s)=ss+α

Фильтр размытия помещает нуль в начале координат, которое ограничивает спиральный полюс режима оставаться около источника. Мы выбираем α = 0.2 какое-то время постоянный из пяти секунд и используем метод корневого годографа, чтобы выбрать усиление фильтра H. Сначала задайте фиксированную часть s / (s + α) размытия

H = zpk(0,-0.2,1);

Соедините размытие последовательно с моделью проекта sys11 (отношение между входом 1 и выходом 1), чтобы получить модель разомкнутого цикла

oloop = H * sys11;

и чертите другой корневой годограф для этой модели разомкнутого цикла.

rlocus(-oloop)
sgrid

Создайте и перетащите маркер данных вокруг верхней кривой, чтобы определить местоположение максимального затухания, которое является о ζ = 0.3.

Этот рисунок показывает маркер данных в максимальном коэффициенте затухания; усиление - приблизительно 2,07.

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

K = 2.07;
cl11 = feedback(oloop,-K);
impulse(cl11,20)

Ответ обосновывается приятно, но имеет меньше затухания, чем ваш предыдущий проект. Наконец, можно проверить, что фильтр размытия решил спиральную проблему режима. Сначала сформируйтесь, полное размытие фильтруют kH (s) (размытие + усиление).

WOF = -K * H;

Затем замкните круг вокруг первой пары ввода-вывода модели MIMO sys и симулируйте импульсную характеристику.

cloop = feedback(sys,WOF,1,1);

% Final closed-loop impulse response
impulse(sys,'b--',cloop,'r',20)

Угловой ответ банка (выводит 2) из-за импульса элерона (вводит 2) теперь имеет желаемое почти постоянное поведение по этой кратковременной системе координат. Чтобы смотреть ответ более тесно, используйте Селектор ввода-вывода в контекстном меню, чтобы выбрать (2,2) пара ввода-вывода.

Несмотря на то, что вы действительно не совсем соответствовали спецификации затухания, ваш проект увеличил затухание системы существенно и теперь позволяет пилоту обычно управлять самолетом.

Для просмотра документации необходимо авторизоваться на сайте