книги из ГПНТБ / Кашин Г.М. Автоматическое управление продольным движением упругого самолета
.pdfП ередаточн ая функция реш аю щ его усилителя определится ф орм улой
|
|
|
|
|
|
( Р ) |
1 |
|
1 |
|
|
|
|
(4. 40) |
||
|
|
|
|
|
|
В |
У\о.с |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где |
В = |
|
|
|
С \ В \ |
-- |
A \ D \ |
С |
---ЛтАт ^ 5 * 5 ----A5 D 5 |
|
||||||
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
В \ |
|
^ |
|
в ,з |
|
+ |
в 5 |
|
|
|
|
Если ветви входной цепи Z lt |
Z 3 |
и Z5 |
решающего |
усилителя реализуются |
|||||||||||
Т-образными цепочками Л?С, то |
|
|
1 |
|
|
|
|
|
|
|
|
|||||
В : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Z \ \ -Г Z i 3 - f Z n F ^ Z i.3 |
Z31 + Z 33 + 2 3 1 Г з2 ^ зз |
Z51 -г Z 5 3 _ + Z 5 1 Kg2 Z5 3 |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.41) |
|
При структуре входной цепи, показанной |
на |
рис. |
4. |
19, а, |
выражение |
для |
В |
|||||||||
запишется в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
В = |
|
1 |
|
|
|
*1*зЧ/* |
|
|
|
R\C^p |
(4. 42) |
||||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
2 * i ( * А /> |
+ |
1) |
2 * ! |
( * 3 С3/> + |
1) |
2 * x (R 5C 5p + 1 ) |
|
|
||||||
При условии R \ C X= * 3 C3 = R 5C 5, получим |
|
|
|
|
|
|
||||||||||
|
|
|
B |
= |
|
2*i (*iCijP + 1)__ |
|
|
(4. 43) |
|||||||
|
|
|
* 1 * з ф 2 + R \ C s p 4- 1 |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||
|
Если проводимость цепи обратной связи у 0.с |
(см. |
схему на рис. 4. |
15) |
вы |
|||||||||||
брать равной |
(4.26) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
У о.с = |
R\CчСа,р2 + *,С4/у 4- 1 |
|
|
(4. 44) |
|||||||||
|
|
|
|
|
2*2 (R2C2P + 1) |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
'(*) = |
|
* 2 |
(R2C2P + 1) (* 1*зС3/?2 + |
RlCeP-r 1) |
|
|
||||||||
|
евх |
* 1 |
( * А / > |
+ |
1) ( * 2 ^ 2 ^ 4 * 2 + |
R iC ^ p |
-f- 1) |
|
|
|||||||
|
|
|
|
|||||||||||||
Следовательно, при условии |
(см. рис. |
4. 21, а) |
|
|
|
|
|
|
|
|||||||
|
|
|
* А |
- R 2C 2 = * 3с 3 = R 5C 5, |
|
|
|
|
||||||||
|
|
|
(/>) = |
|
*2 |
*1*зС3Р2 + R1C5P + 1 |
|
(4. 45) |
||||||||
|
|
евх |
|
* 1 |
R y C y C ^ p ^ |
-f- R 2C 4 P “Е 1 |
|
|||||||||
|
|
|
|
|
|
|
||||||||||
|
Аналогичную передаточную функцию можно реализовать, если |
|
|
|||||||||||||
|
В = |
— |
1 |
1 |
|
|
|
|
|
2 * i ( R xC xp |
+ |
1) |
|
|
||
|
|
+ R \ C ? ,p |
* 2 C iC 3/> 2 + * iC 3p + 1 |
|
|
|||||||||||
|
|
|
|
|
|
|||||||||||
|
|
2* i ( R \ C xp + |
|
|
|
|||||||||||
|
|
1 ) |
2 * i |
|
|
|
|
|
|
|
|
|
||||
Это соответствует электрической схеме решающего усилителя, показанной |
на |
|||||||||||||||
рис. |
4. 21, б, и при той же цепи обратной связи и при R \ C \ = R 2 C 2 |
|
|
|||||||||||||
|
|
|
■(*) = |
* |
2 |
|
R \ C xC 3Pi |
+ R \ C ?Jp + |
1 |
|
(4.46) |
|||||
|
|
|
*i |
|
I&C2C4P2 + R2C4P |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
200
|
|
|
Т 1 |
т |
Эта схема имеет ограничения и может быть реализована при условии-^—= |
-у— |
|||
\ [см. формулу (4.34)]. Передаточную |
функцию |
(4.34) можно |
реализовать |
|
с меньшим количеством элементов. Возьмем в качестве базовой |
синтезирован |
|||
ную ранее схему на рис. 4. 17, в. Для |
получения |
необходимой |
передаточной |
функции требуется определить параметры дополнительной прямой цепи с поо-
1
водимостьюy-j = — (см. рис. 4. 21, б). На рис. 4. 21. г показан ее граф, состав ит
Рис. 4.21. Реализация' передаточных функций вида
евь,х |
* № + 2^ |
+ |
1) |
<?вх |
Т\р* + |
+ |
1 |
ленный из комбинации в прямой цепи ygyg и в цепи обратной связи hy гра фов. Передаточная функция определится по правилу Мезона (см. рис. 4. 21, г).
евых |
_ |
l . P kDk |
___________ 2567 + 167(1 + 232 + 2432) |
евх |
(Р) ~ |
D |
~ 1 — 7867 — 78432567 +232-7867 +2432-7867 |
(4.47)
(при й->-оо контуры и произведения контуров, не содержащие k, исключаются), тогда (см. рис. 4. 21, в)
^вых . . _ |
RiCip + R.2 [1 + R5 (Cl + |
С4 ) Р] |
-_______________________ ^7 |
||
евх Р |
C 1R2 C4 R5 P2 + Т?5 (Cj + |
С4) р + 1 |
Для получения передаточной функции (4. 34) необходимо, чтобы
1RjC7p + 1
=Ri '
201
тогда (см. рис. 4. 21, б) |
|
|
|
R2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
'(/>) = - |
vRix |
■C4)] /> + 1 |
X |
Rf, (С1 |
|
|
■C7) + Rf(C1 |
||
C 1R2 C4 |
RFP2 |
+ |
/?,5 (C1 + |
C4) p 4- 1 |
(4. 48) |
|
|
|
С помощью такого фильтра можно реализовать передаточную функцию, когда Ti>T2 и £i> £2,
или
%{ГХ- 2£2Г2 >
2С2^2
Необходимо отметить, что схема, показанная на рис. 4.21, а, является универсальной, ибо здесь каждая ветвь прямой цепи соответствует своему коэффициенту полинома числителя. Исключением соответствующих ветвей пря мой цепи схемы, показанной на рис. 4.21, а, можно реализовать следующие передаточные функции [см. формулу (4.42)]:
^вых |
^ ^ _ |
_ |
R% _____ R\C*Р |
\ 1____ |
|
||
евх |
|
|
R\ |
R^CyCqp'Z + R2C4P + 1 |
|
||
при RxCi = R5C$ — R2C2; |
|
|
(4. 49) |
||||
евых |
(p ) = |
|
R2 |
R\RsPIp2 + 1 |
|
||
---------------- |
------ |
------------- ----------------------------------------------------------------------------- |
|
||||
евх |
|
|
^1 |
R^Pj ^ |
4 P^ 4- R2 C4 P 4~ 1 |
|
|
при R XCX= R3 C3 = R 2 C2 - |
|
|
|
|
|||
4 . 3 . КРИТЕРИИ |
КАЧЕСТВ А |
И ОЦЕНКА |
ОПТИМАЛЬНОСТИ |
СИСТ |
|||
АВТОМ АТИЧЕСКОГО |
УПРАВЛЕНИЯ |
УПРУГОГО С А М О ЛЕТА |
|||||
Как отмечалось ранее, синтез структуры регулятора методом |
|||||||
инверсии графа предполагает |
получение передаточной функции |
||||||
(последовательно |
включенных основного и инверсного графов) |
||||||
равной единице. |
Такой |
регулятор |
реализовать трудно, да |
и |
в этом нет необходимости. При реализации приходится вводить некоторые операторы 0 (р). В ряде случаев целесообразно аппроксимировать свойства передаточных функций лишь в опре деленном диапазоне частот. Соблюдая условия реализуембсти, получим корректирующие устройства, отличающиеся от синте зированных по методу инверсии графа.
Критерием оптимальности спроектированной системы управ ления может явиться оценка качества регулирования. Такая оценка обычно осуществляется по качеству переходных процес сов или по точности стабилизации регулируемого параметра.
Для пилотируемых летательных аппаратов качество замкну той системы стабилизации и управления должно оцениваться с двух точек зрения. Во-первых, должна быть обеспечена необ ходимая точность стабилизации летательного аппарата при дей ствии на него возмущений турбулентной атмосферы и, во-вто-
202
рых, замкнутая система должна удовлетворять вполне опреде ленным требованиям управляемости.
В качестве примера на рис. 4.22 представлена схема стаби лизации и управления углом тангажа самолета. Здесь коррек тирующие устройства WI;2 и WI;3 должны быть выбраны из усло вия необходимого быстродействия и точности стабилизации угловых положений самолета при действии возмущений ветра W n ( t ) .
S, Коррек. Коррек. устройЧХ)* устрой ство ство
WK4 WK,
Коррек.
4
устрой Привод Самолет стВо
W,
-ДУС-
- ц г в ■
Рис. 4. 22. Структурная схема стабилизации и управления углом тангажа самолета
Корректирующее устройство WKi во входной цепи выбирается из условия удовлетворения требованиям управляемости само лета летчиком.
При разработке и исследовании САУ упругого самолета ока залось удобным производить оценку качества стабилизации регу лируемых параметров методом модифицированного корневого годографа, методами минимизации интегральных квадратичных оценок и интегральных оценок произведения весовых функций.
1. Оценка устойчивости и качества регулирования методом модифицированного корневого годографа
При исследовании влияния изменения параметров упругого самолета на его устойчивость, при определении критической ско рости флаттера, а также при определении параметров системы управления необходимо иметь метод, позволяющий сравнительно легко выявить динамические свойства системы высокого поряд ка. В этом случае может быть использован метод, вытекающий из анализа корневых годографов. При этом критерий устойчиво сти формулируется следующим образом: для устойчивости систе мы порядка п необходимо и достаточно, чтобы все корни tji не четной А\(р) и корни Xi четной А2(р) частей характеристическо го уравнения А(р) лежали на мнимой оси в плоскости р, а их абсолютные значения находились в соотношении
у0< х 1< у 1< х 2< у 2<хз< ■ |
. (4.50) |
где у0 находится в начале координат плоскости р.
20 3
Так, характеристическое уравнение вида |
|
|
|
||||
|
|
Л(р) = ^ a |
kpk--=0 |
|
|
(4.51) |
|
|
|
к=О |
|
|
|
|
|
можно представить следующим образом: |
|
|
|
||||
|
А (р) = А2(р) Л-А1(р) = 0, |
|
|
(4.52) |
|||
где в А2(р) и А\(р) |
входят соответственно |
четные и |
нечетные |
||||
члены уравнения (4.51). Разделив |
уравнение (4.52) |
на А\(р), |
|||||
получим |
|
|
|
|
|
|
|
|
|
А 2 ( р ) = |
_ |
! |
|
|
|
|
|
■Д (/>) |
|
|
|
|
|
или в развернутом виде |
|
|
|
|
|
||
А 2 (р ) |
С п ( р т |
+ а т _ 2Р т ^ + • • ■+ а'2Р 2 + Др) |
|
(4.53) |
|||
(/О |
р ( / |
+ a'I_ 2 pl ~ 2 + |
|
+ а'3 р 2 + |
а[) |
|
|
|
|
|
|||||
где Сл = ^=4; |
1= т — п — 1 — для |
нечетных |
значений га и |
||||
ч |
|
|
|
|
|
|
|
Сп = ———; 1= п — 2; гаг = га для |
четных значений |
га. |
|
||||
i |
|
|
|
|
|
|
|
Для исследования устойчивости необходимо определить нули и полюсы выражения (4.53). При нарушении указанной в (4.50) последовательности полюсов Уг и нулей лу система будет неустой
чивой. В случае, когда полюсы выражения —— (4.53) равны
М р )
нулям, получаются мнимые корни характеристического урав нения (4.51), и система оказывается неустойчивой. В случае близости полюса к нулю или больших (относительно) значений коэффициента Сп система оказывается слабодемпфированной.
При |
достаточном расхождении полюсов и нулей |
уравнения |
(4. 53) |
система получается хорошо демпфированной. |
При малых |
значениях коэффициента Сп корни характеристического уравне
ния (4.51) располагаются вблизи |
полюсов |
выражения |
(4.53), |
|||
а при больших значениях Сп корни |
(4.51) |
лежат вблизи |
нулей |
|||
(4.53). Для пояснения этого метода |
рассмотрим систему ше |
|||||
стого порядка, |
характеристическое |
уравнение которой |
будет |
|||
иметь вид |
|
|
|
|
|
|
А(р) =а6р6+ а5р5+ ... + а1р + ао = 0. |
|
|||||
Аналогично (4. 53) запишем |
|
|
|
|
||
(р) |
Ч (р6 + ^ - Р 4 + — р2 + — ) |
|
||||
V |
ч |
ч |
ч |
) |
|
М ( Р )
ЧР 1 Р 4 + — Р г +
ча5
2 0 4
Если, например, для СТС с учетом двух тонов упругих коле баний (3. 33)
.4 (/;) = р» + 92р* + 4325р4 + 70584/?3 + 762200р2 + 3994000р +
+ 30070000 = 0, |
то |
|
|
|
Л2 (р) _ |
рб + 4325/И + 762200^2 + 30070000 |
|
|
|
А г (р ) ~ |
92р5 + 70584/?з + |
3994000/? |
|
|
или |
|
|
|
|
х 3 + 4325x2 + 762200х + |
30070000 |
1, |
(4. 54) |
|
|
92i/o (У2 f 770у + 43500) |
|||
|
|
|
||
|
Q _ |
1 |
|
|
|
6~ а5 ~ |
92 |
|
|
Здесь в уравнении (4.54) произведена замена: в числителе р2= х и в знаменателе р2= у и р = Уо-
Нули и полюсы уравнения (4.54) имеют следующие зна чения:
Уо |
X, |
У\ |
•*2 |
2 |
х 3 |
|
|
У |
|
||
0 |
48 |
150 |
303 |
1190 |
2820 |
Такое расположение нулей х, и полюсов Уг свидетельствует об устойчивости системы.
Часто при анализе и синтезе систем автоматического управ ления требуются более детальные_представления о динамиче ских свойствах системы. Так, если система состоит из звеньев, каждое из которых характеризует свою степень свободы и опре деленные динамические свойства, то интересно знать, как на раз личных режимах полета изменяются характеристики каждого звена в системе. Такое, более четкое, представление о системе можно получить, выполнив графическое построение кривых взаимного расположения корней уравнения (4.54).
Для определения степени влияния некоторых параметров системы запишем характеристическое уравнение выше рассмат риваемой системы в виде
М Р ) = {Р 2+ 2^ % Р + |
“*) (Р2+ |
2!> з 0 + |
ш!) {Р2+ |
Р + “!) = °- |
||
Здесь |
|
|
|
|
|
|
+ |
= 61, |
ш2=184, |
и 0)2= 2687 в |
|
||
в- |
|
|
3 |
4 |
С 2 |
|
Рассмотрим |
поведение системы при |
изменении |
декрементов |
|||
затухания £*>; £3; |
£4. Для этого будем определять корни х* и у\ |
|||||
уравнения, (4. 54) |
^ |
для различных вариантов значений &■ |
||||
|
|
М Ш У) |
|
|
|
205
£з и £4. Затем выполним графическое построение характеристик устойчивости, как показано на рис. 4.23. Здесь по оси абсцисс на равных произвольных расстояниях произведена разметка кор ней i/j и I; в последовательно
У,:Л- |
сти, нужной для устойчивости |
||
|
системы, а по осп ординат от |
||
|
ложены их численные значения |
||
|
(последние целесообразно от |
||
|
кладывать в логарифмическом |
||
|
масштабе). Из анализа кривых |
||
|
можно |
установить |
некоторые |
|
свойства такой системы при из |
||
|
менении ее параметров. Для оп |
||
|
тимального варианта, когда де |
||
|
кременты затухания отдельных |
||
|
звеньев |
= ^3 = ^4=0,7 (кри |
|
|
вая 1, |
корни уи Xi |
находятся |
|
почти на прямой линии. |
Рис. 4. 23. Кривые расположения плю сов i/i нечетной Ai\y{) и нулей-х» чет ной A2 (Xi) частей характеристическо го уравнения Л(х,; уг) при различ ных значениях декрементов затухания
Неустойчивость системы наступает, когда соблюдается хотя бы одно из следующих соотношений:
|
— > 1 ; |
— |
1; — > 1 ; |
— !>1. |
(4.55) |
||
|
У1 |
*2 |
У2 |
|
*3 |
могут |
быть прибли |
Частоты слабодемпфированных колебаний |
|||||||
женно определены по формулам: |
|
|
|
|
|||
|
Х 1 |
■м2, |
У1+ -*2 |
. ...2_ |
Х2, + Уч |
(4. 56) |
|
ш —-----+ У\ |
2 |
|
4 |
2 |
|||
& |
2 |
|
’ |
|
|||
а декременты затуханий — по соотношениям: |
|
||||||
|
У1— *1 . |
*2 — Ух ; с4 |
Хз — У2 |
||||
Кроме того, |
Ml |
|
х 2 |
|
|
х3 |
определяется |
степень |
демпфирования |
системы |
еще величиной коэффициента Сп (см. рис. 4.23, кривая 5).
2. Оценка качества регулирования на основе интегральных квадратичных оценок и интегральных оценок произведения весовых функций
Для синтеза и анализа контуров САУ упругого самолета, когда динамика движения описывается системой дифферен
2 0 6
циальных уравнений высокого порядка, при выборе параметров элементов цепей самонастройки эффективным оказался метод минимальных квадратичных оценок и интегральных оценок ве совых функций:
|
|
о |
|
|
|
|
(4. 57) |
|
|
|
|
|
|
|
|
|
|
о |
|
|
|
|
|
где |
(t) и Шр.;- (t) |
— весовые функции, |
выражающие реакции, |
||||
г-й и /-й координаты |
на 6-импульс, поданный соответственно на |
||||||
v-й и p-й входы системы. |
|
|
|
|
|
||
|
Точками приложения 6-импульсов в реальных системах |
||||||
используются точки |
(источники) |
действий |
возмущений |
ветра, |
|||
управляющих сигналов и т. д. |
|
|
|
|
|
||
|
Точками измерения реакции системы являются точки (стоки) |
||||||
измерения параметров системы, т. е. места |
установки датчиков |
||||||
(ДЛУ, ДУС, ДУУит. д.). |
|
|
|
|
|
||
|
При синтезе отдельных блоков самонастраивающейся си |
||||||
стемы точками приложения 6-импульсов |
и измерения сигналов |
||||||
могут быть соответственно точки |
входа |
и |
выхода |
цепи |
само |
||
настройки. |
|
|
|
порядка |
могут |
быть |
|
|
Критерии (4.57) для систем высокого |
|
|||||
определены на аналоговой или цифровой |
вычислительных ма |
||||||
шинах. |
|
|
|
|
|
|
Глава V
АНАЛИЗ И СИНТЕЗ КОНТУРОВ СТАБИЛИЗАЦИИ САУ УПРУГОГО САМОЛЕТА С УЧЕТОМ НЕЛИНЕЙНОСТЕЙ ПРИВОДА
При выборе структуры САУ, корректирующих устройств и фильтров нужно учитывать влияние некоторых существенных нелинейностей. Такими нелинейностями являются нелинейности типа зоны нечувствительности и насыщения, которые имеют ме-
Рис. 5. 1. Схема управления угловой скоростью самолета
сто в характеристиках некоторых датчиков и исполнительных элементов. Наличие этих нелинейностей может приводить к уменьшению запаса устойчивости, возникновению автоколе баний или неустойчивости.
В данной главе рассматриваются возможные состояния не линейных систем и методы их исследования, предлагаются кор ректирующие устройства для устранения вредного влияния не линейностей.
При приближенных исследованиях нелинейных систем авто матического управления эффективными оказались два метода: метод логарифмических частотных характеристик [17], и моди фицированный метод корневого годографа. Они позволяют при
•208
исследовании нелинейных многоконтурных систем высокого |
по |
|||||||||
рядка получать обозримые результаты с минимальной затратой |
||||||||||
времени. |
5. 1 |
представлена схема |
управления угловой |
|
ско |
|||||
На рис. |
|
|||||||||
ростью самолета с учетом двух тонов |
упругих колебаний, |
где |
||||||||
показаны характерные нелинейности типа «зона нечувствитель |
||||||||||
ности» и «ограничение». |
них |
являются |
|
нелинейности, |
представ |
|||||
Существенными из |
|
|||||||||
ляющие собой статические скоростные |
характеристики y'=f(y) |
|||||||||
и e' = f(e) |
соответственно сервопривода |
автопилота и |
бустера. |
|||||||
Эти.нелинейности будут учитываться при анализе и синтезе |
||||||||||
контуров САУ. |
нелинейных |
систем |
будет |
использован |
метод |
|||||
При анализе |
||||||||||
аппроксимации |
нелинейных элементов |
эквивалентным |
коэффи |
|||||||
циентом усиления. |
|
|
|
|
|
|
|
|
||
5 . 1 . АППРОКСИМАЦИЯ |
Х А РА К ТЕ РИ С Т И К |
НЕЛИНЕЙНЫ Х |
ЭЛЕМ ЕНТ |