книги / Системы управления летательными аппаратами и их силовыми установками
..pdfВывод: при статически неустойчивом ЛА СУС структурно неус тойчива, поэтому необходимо ввести в систему автомат стабилиза ции, который обеспечит устойчивость углового движения ЛА.
При статически устойчивой ЛА СУС находится на границе ус тойчивости. Если учесть влияние движения центра масс ЛА на угло вое движение, а также инерционность элементов системы, то угловое движение, как и в предыдущем случае, будет неустойчиво. Следова тельно, в данном случае также необходимо ввести автомат стабили зации, обеспечивающий устойчивость углового движения ЛА.
14.3. Обоснование закона управления СУС
Под законом управления будем понимать математическую за висимость между регулирующим параметром (в данном случае уг лом поворота рулевых органов) и регулируемым параметром (углом рыскания) без учета инерционности элементовавтомата ста билизации:
S = /(¥ )• |
(14.8) |
1. Пусть в СУС реализуется закон управления по отклонению регулируемого параметра (по координате):
Ь=КЧ1у. |
(14.9) |
Тогда будем считать, что передаточная функция АУС
frv (p) = Kv . |
(14.10) |
Передаточная функция разомкнутой скорректированной СУС, т.е. при учете АУС, запишется в виде
Щ р ) = к Л ь ~ —V • |
О4-11) |
P ^4'V |
|
Представим характеристическое уравнение замкнутой скоррек тированной системы:
\ + W2(p) =0 |
(14.12) |
P2+ КхА\‘Ь+ |
= 0 . |
( 14. 13) |
Как видно из уравнения (14.13), система будет находиться на |
||
границе устойчивости при любом знаке bvv, если К |
^ 6 >|èw | . Сис |
|
тема будет неустойчива при \ 1Ц1< 0 и КЧЬЧ8 < |ôw | |
|
Таким образом, при реализации в СУС закона управления по от клонению регулируемого параметра имеется качественная аналогия со случаем, когда управление движением отсутствует, что является
неудовлетворительным. |
|
|
2. |
Пусть в СУС реализуется закон управления по отклонению |
|
регулируемого параметра и скорости изменения данного отклонения |
||
(закон управления по скорости и координате): |
|
|
|
8,|(= К ,Д + ^М /- |
(14.14) |
В этом случае передаточная функция АУС примет вид |
|
|
|
W4(p) =Kv{TKp +\), |
(14.15) |
/С, где Тк = -4 -.
%
Тогда передаточная функция разомкнутой скорректированной
системы запишется следующим образом: |
|
|
Щ (Р) = |
■ |
(14.16) |
Р |
+04/1(1 |
|
Представим характеристические уравнения замкнутой системы: |
||
P2 + KvbvbTKp + КцЬу6 +0W =0. |
(14.17) |
Условие устойчивости системы в соответствии с критерием Гурвица примет вид
K4bvS+bv v >0. |
(14.18) |
Как видно из формулы (14.18), при ЬЩ1 <0 выполнение данного
условия зависит от выбора величины коэффициента ЛТ(|).
Итак, реализация в СУС закона управления по скорости и коор динате обеспечивает структурную устойчивость системы. Соответст
вующий выбор параметров АУС Kv , |
позволяет добиться тре |
буемого качества регулирования. |
|
14.4. Особенности стабилизации углового движения жесткого ЛА с БЦВМ в конторе управления
Осуществим анализ устойчивости канала рыскания дискретной СУС, для чего произведем выбор передаточной функции дискретно го вычислительного устройства исходя из обеспечения устойчивости и качества регулирования, оценим влияние основных параметров дискретного автомата стабилизации (периода квантования, коэффи циента передачи) на динамику системы, осуществим выбор данных параметров, учитывая требования к точности СУС.
1.Ст рукт урная схема дискретного канала ры скания СУС
Используя упрощенные уравнения движения ЛА, а также функцио нальную схему дискретной системы стабилизации (см. рис. 13.4), со ставим структурную схему канала рыскания дискретной СУС.
Для упрощения процедуры исследования устойчивости введены следующие допущения:
- будем считать ЛА статически нейтральным, т.е. положим ко эффициент = 0;
- будем полагать, что чувствительный элемент системы угловой стабилизации - трехосный гиростабилизатор - выдает информацию об углах отклонения ЛА от программных значений и имеет переда точную функцию, равную WT(p) =Kv;
рулевой привод будем считать безынерционным элементом с передаточной функцией fVn(p) = К„ .
Структурная схема канала рыскания дискретной СУС пред ставлена на рис. 14.5.
Рис. 14.5
Следует отметить, что в качестве экстраполятора в системе использовано запоминающее усройство нулевого порядка (фикса тор). D(z) - передаточная функция ДВУ.
2. Анализ устойчивости канала рыскания. Для анализа устой чивости определим z- и w-передаточные функции разомкнутой не скорректированной системы, затем, применив один из критериев ус тойчивости, оценим устойчивость системы.
Итак,
(14.19)
При анализе устойчивости будем полагать, что D(z) = 1. Используя таблицу z-преобразований (приложение 1), получим
W\(z) =КгКпЬуЪT0\ z +1) |
(14.20) |
|
|
2 (z -l)2 ' |
|
Выполнив подстановку |
1 + W |
w-передаточную |
Z — Z » определим |
||
функцию системы: |
1 - W |
|
|
|
|
Wl(w) =KrKab^TQ21- w |
(14.21) |
|
|
4w2 |
|
Для анализа устойчивости применим критерий Гурвица, для че го найдем характеристическое уравнение системы в виде
\ + W^(w)-Q\
4w2 - КгКпЬфТ2w+ КгКпЬ ^Т 2 = О
Осуществив анализ данного уравнения, можно отметить, что система структурно неустойчива.
3. Выбор передаточной функции ДВУ, построение областе устойчивости СУС. В связи с тем, что знак второго члена харак теристического уравнения отрицательный, в состав передаточной функции ДВУ должен входить член, содержащий w. Таким образом,
в общем случае дискретное вычислительное устройство должно представлять собой форсирующее звено.
Выберем передаточную функцию вида |
|
£>(w) = /:K(7 > + l). |
(14.22) |
Тогда w-передаточная функция разомкнутой скорректированной системы запишется следующим образом:
Ц М - W |
о - « O t e » * 0 |
(14 .23) |
|
4w |
|
или |
|
|
W2(w) = K0b STQ (I~ Wfe v ^+_0 |
( ,4 24) |
|
|
4w |
|
где KQ- коэффициент передачи АУС, |
|
|
|
К0 =КгКкК„. |
(14.25) |
Характеристическое уравнение системы в этом случае при |
||
мет вид |
|
|
(4 - K0bvST2TK)w2 + K0bwST2 (Тк - 1V + * А Л 2 = 0 • |
(14.26) |
|
Запишем условия устойчивости системы: |
|
4 - K obvST X > 0 , |
(14.27) |
||
TK>I- |
|||
|
|||
Из (14.27) получим зависимость для определения Тк |
|
||
1<Г„ < |
4 |
(14.28) |
|
|
кК о К Х '
Период квантования можно вычислить с помощью зависимости
Т0 < |
|
|
|
(14.29) |
Для |
анализа |
влияния |
периода |
|
квантования построим области устой |
||||
чивости системы в плоскости пара |
||||
метров |
Тк, |
К0, |
изображенные на |
|
рис. 14.6, для чего используем зави |
||||
симость |
(14.28). |
На |
рис. 14.6 |
|
TQ> TQ> T" ; |
K0m - |
максимально до |
пустимый коэффициент передачи оп ределяется из выражения (14.28) при
г,= 1 . |
|
к » , |
04.30) |
®Ч'5''0 Анализ областей устойчивости показывает, что увеличение пе
риода квантования, т.е. уменьшение частоты квантования системы отрицательно сказывается на качестве регулирования системы, так как происходит уменьшение областей устойчивости. Кроме того, увеличение коэффициента, характеризующего эффективность руле вых органов {ЬЦ1&), также отрицательно сказывается на качестве ре
гулирования системы.
Так как эффективность рулевых органов в процессе полета воз растает (см. рис. 13.5), то рабочую точку (значения Тк, Кк) необхо
димо выбирать для конца активного участка, что обеспечит выпол нение условий устойчивости на всей траектории полета ЛА.
14.5. Анализ точности дискретного канала рыскания СУС
Осуществим анализ точности канала рыскания СУС и оценим влияние на точность СУС коэффициента передачи и периода кванто вания. Будем считать, что на систему действует постоянное во вре мени возмущение
Âf¥0>) = — . |
(14.31) |
Р |
|
Используя теорему о конечном значении [2], определим ус |
|
тановившуюся ошибку системы: |
|
Z—1С М^(р)
„ ^ „2
|
(14.32) |
Здесь |
|
W1{z) = W[{z)D{z\ |
(14.33) |
D(z) определяется из зависимости (14.22) путем использования под-
становки z = 1+ w
1 —w
(14.34)
Итак,
(14.35)
Учитывая формулу (14.20), а также используя таблицу z-преоб- разований, в итоге получим
£ z i * & £ ± £ j |7 |
¥ |
|
z 2(z - |
1)3 |
|
z + 1 |
(rKH-l)z + l - r K |
|
2 (z -l)2 |
z + 1 |
_ —М„у . (14.36) *<А„5
Анализ зависимости (14.36) показывает, что уменьшение ус тановившейся ошибки, т.е. увеличение точности системы, происхо дит при увеличении коэффициента передачи К0. Установившаяся ошибка системы не зависит от периода квантования, однако учет влияния его на точность можно осуществить, рассмотрев зависи мость (14.30). Данная зависимость показывает, что увеличение пе риода квантования приводит к уменьшению максимально допусти мого значения коэффициента передачи, что отрицательно влияет на точность системы.
Таким образом, увеличение периода квантования отрицательно сказывается как на устойчивости, так и на точности системы.
Осуществим выбор коэффициента передачи, исходя из требо ваний к точности системы. Значение установившейся ошибки ог раничивается допустимым значением:
Vy ^4 V |
(14.37) |
Подставив (14.37) в (14.36), получим
М ш
(14.38)
14.6. Способ повышения точности стабилизации движения ЛА по каналу тангажа
1. Ст рукт урная схема дискретного канала т ангаж а СУС.
Запишем зависимость для отклонения по дальности, ограничиваясь линейными членами разложения в ряд Тейлора:
dL |
• |
( dL |
Ay + |
f dL \ |
Az + |
dL |
AL = |
Дх+ |
— |
|
Ax + |
||
dx j |
|
[dy |
p |
\d z j p |
dx |
|
|
dL\ |
|
( d £ |
|
^ dL^ |
(14.39) |
|
A ' |
Az + |
dt. |
|||
|
|
|
|
\d t j |
||
|
# Г |
+ |
<SzJD |
|
Анализируя данное выражение, можно отметить, что наиболь шее влияние на ошибку по дальности полета оказывают члены, учи тывающие влияние отклонений продольной и нормальной проекций вектора скорости от расчетных значений, т.е. члены
dL |
Лх + f dL^ Aÿ |
dx |
\ f y j |
На ряде ЛА уменьшение данных отклонений осуществляется, главным образом, с помощью системы регулирования кажущейся скорости, системы стабилизации движения центра масс в направле нии по нормали к траектории полета или сокращенно - системы нор мальной стабилизации. Однако на ЛА с двигателем, работающим на твердом топливе, система регулирования кажущейся скорости от сутствует.
Для того чтобы при этом точность попадания в заданную точку земного пространства не ухудшалась, усложняются алгоритмы наве дения ЛА (производится учет нелинейных членов разложения в ряд Тейлора), а также увеличивается точность стабилизации ЛА по углу тангажа, что обусловливает уменьшение отклонения угла наклона вектора скорости от расчетного значения.
Решение данной задачи обеспечивается путем введения интег рала в закон управления по углу тангажа:
6Э=К9AS + КйА&+Kf \АШ . |
( 14.40) |
Сравнив это выражение с зависимостью (14.14), видим, что в от личие от канала рыскания в закон управления по углу тангажа вво дится член ATj jAvk?/. Рассмотрим влияние данного члена на точность
и устойчивость канала тангажа дискретной СУС. Для решения зада чи представим структурную схему дискретного канала тангажа (рис. 14.7).
Как видно на рис. 14.7, дискретное вычислительное устройство включает в себя, кроме форсирующего звена с передаточной функ цией D(z), дискретный интегратор с передаточной функцией £>H(z), подключенный параллельно форсирующему звену.
Известно, что в дискретных вычислительных устройствах ин тегрирование осуществляется по правилам численного интегриро вания (см. приложение 2).
Наибольшее применение в системах управления ЛА нашел спо соб численного интегрирования по правилу трапеций, обес печивающий достаточно высокую точность интегрирования при сравнительно простой реализации. Передаточные функции интег ратора, интегрирующего по правилу трапеций, можно представить в виде
T„z + l |
тх и l + z~‘ |
А ,(*)=*« 2 Z - 1 |
=К« 2 l - z " ,; |
|
( 1 4 . 4 1 ) |
ПО |
|