книги из ГПНТБ / Кашин Г.М. Автоматическое управление продольным движением упругого самолета
.pdfРис. 3. 14. Определение передаточной функции системы по методу одинарных некасающихся контуров
150
Необходимо построить общий граф G системы и подграфы G0 и Gh. Для
пояснения принципа |
построения |
общего |
графа |
систему |
|
(3 |
14) |
представим |
||||||||||
а матричной форме |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
1 |
0 |
0 |
|
0 |
0 |
0 |
0 |
0 |
0 |
0 |
|
0 |
|
0 |
|
0 |
|
1~' |
— в |
A\ |
0 |
|
0 |
0 |
0 |
0 |
0 0 0 |
|
0 |
|
0 |
|
0 |
0 |
|||
0 |
-- 1 |
1 |
|
0 |
0 |
0 |
0 |
1 0 0 |
|
0 |
|
1 |
|
0 |
|
0 |
||
0 |
0 |
~ B 4 a 4 0 |
0 |
0 |
0 0 0 0 |
|
0 |
|
0 |
0 |
||||||||
0 |
0 |
0 |
|
0 |
An |
0 |
0 |
0 0 0 |
|
0 |
|
0 |
|
0 |
|
0 |
||
0 |
0 |
0 |
— 1 — 1 |
1 |
0 |
0 0 0 |
|
0 |
|
0 |
|
0 |
|
0 |
||||
0 |
0 |
0 |
|
0 |
о --- B s |
-4s |
0 0 0 0 0 |
|
0 |
0 |
||||||||
0 |
0 |
0 |
|
0 |
0 |
0 |
— в |
-4.3 |
0 |
0 |
|
0 |
|
0 |
|
0 |
|
0 |
0 |
0 |
0 |
|
0 |
0 |
0 ~ B ( |
0 |
-46 |
0 |
|
0 |
|
0 0 |
|
0 |
|||
0 |
0 |
0 |
|
0 |
0 |
0 |
0 |
0 |
~ B 7 At |
|
0 0 0 0 |
|||||||
0 |
0 |
0 |
|
0 |
0 |
0 |
0 |
0 |
- B a 0 |
-48 |
|
0 0 |
0 |
|||||
0 |
0 |
0 |
|
0 |
0 |
0 |
0 |
0 |
0 |
0 |
- B g Л5 |
0 |
|
0 |
||||
0 |
0 |
0 |
|
0 |
0 |
0 |
0 |
0 . 0 |
0 |
— Вio |
0 |
A io |
0 |
|||||
0 |
0 |
0 |
|
0 |
0 |
0 |
0 |
0 0 — 1 |
0 |
|
0 — 1 1 |
|||||||
|
|
|
|
|
|
X \ |
|
x 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x 2 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
*3 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
*4 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x 5 |
B 2x о |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
*6 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x 7 |
|
0 |
|
|
|
|
|
|
|
|
(3. 30) |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x 9 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x 10 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X 11 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-«12 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x 13 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
_ x 14_ |
|
0 |
|
|
|
|
|
|
|
|
|
На рис. 3. 14, б построен общий граф G. Элементы диагонали матрицы (3.30), являющиеся знаменателями передаточных функций отдельных звеньев схемы, представляют собой коэффициенты передачи собственных контуров па раметров х, системы.
|
На рис. 3. |
14, в и г построены подграфы |
Gо и Gk для вычисления |
знаме- |
|
|
х |
х ю |
|
нателя и числителя передаточной функции (3. 30) — . |
|
|||
|
Из графа |
|
*0 |
|
. |
G0 знаменатель (характеристический полином) передаточной |
|||
хю |
|
|
|
|
функции — определится из следующего выражения: |
|
|||
|
*0 |
|
|
|
2 ( |
- 1)*' П = |
( Ь Д г Ь Д з - А г 1 - А 5- А 6. А 7 - А 3- А д - А 10- D ( - I ) 13 + |
|
|
|
+ |
(Bn,B4B$-A§A7A%AgAiQ- 1 -А \-1) ( — I)9 + |
|
|
|
-h (BiB5B6BsBQ-1-АхАъ'АтAw-1) ( — l)7 + |
|
||
|
+ ( — В ^ - В 5В6-В ф 10-A3-AjAq) ( — l)4 + |
|
||
|
+ |
( В & В ц В ф г - А з А в А э А ю ) ( - |
D 5-. |
(3- 31) |
151
Из |
графа |
Gi, |
(рис. 3. 14, а) |
определится числитель передаточной |
функции |
||
*10 |
|
|
|
|
|
|
|
— в виде |
|
|
|
|
|
|
|
Х ° |
N |
п* = |
( 1 М ! • 1 • |
;;At • 1• Л5Л6 • Л3Л9А 1;) ■1 ) ( _ |
1) 1-4 ~ |
|
|
V ( - |
1) ^ |
|
|||||
|
|
-г (В4В5В3А 8А8А9АЮ- 1• 1 ■Ai • Л-2) ( - |
1)9 |
|
|
||
|
|
-Г (В4В ф ф %В9■.4ю• 1 • 1 • /11 ■Л2■ т4з) ( |
1)' |
-г |
|
||
|
|
+ |
( - В1В4В5В6В8В10-А2А3Лд) ( - 1)4 -1- |
|
|
||
|
|
-1- ( — В2В5В1:В7-1-Ау1-А 3-А4-А8-Ад-Аиу1 )( — 1)1" + |
|||||
|
|
+ ( — В ф ф ф ф т А 2Л3А8АдАю-1) t — I)7- |
(3. 32) |
Выражения (3.31) и (3.32) показывают, из каких контуров графов обра зуются составляющие П и П*. При наличии определенных навыков вычисле ние передаточных функций этим способом можно осуществить непосредственно по структурной схеме рис. 3. 14, а.
О п р е д е л е н и е п е р е д а т о ч н ы х ф у н к ц и й м ет од о м п р е о б р а зо в а н и я г
В ряде случаев при анализе и синтезе систем автоматического управления необходимо бывает изучить влияние отдельного участка цепи на поведение всей системы. При этом передаточ ную функцию части схемы удобно определить методом преобра зования графов. Определение передаточной функции части графа сводится к исключению некоторых узлов исходного графа.
На рис. 3. 15 приведены некоторые варианты преобразования графов [36].
И нверсия
Инверсия — это изменение направления пути направленного графа, при этом источник и сток этого пути меняются местами. Посредством изменения направления пути и тем самым измене ния структуры исходного графа представляется возможным раз решить систему относительно переменной, соответствующей но вому стоку на конце инвертируемого пути.
Поскольку в дальнейшем для синтеза регулятора будут использоваться методы инверсии графов, целесообразно пока-' зать два способа инверсии [24, 36].
На рис. 3.16, б, в выполнена инверсия графа, показанного на рис. 3. 16, а, соответственно способом, сохраняющим узлы, и спо собом, сохраняющим ветви. Инверсия, сохраняющая узлы, сохраняет все узловые сигналы, тогда как инверсия, сохраняю щая ветви, сохраняет расположение всех ветвей, но изменяет направление некоторых из них. Для некоторых графов нет необ ходимости удлинять узлы, чтобы произвести инверсию, сохра няющую ветви. Так, если узел имеет не более одной входящей и одной выходящей ветви, то граф уже будет обладать конфигу рацией, пригодной для инверсии, сохраняющей ветви.
152
1 \ J 9 j |
JL |
1-9 |
Рис. 3. 15. Упрощение графов
153
Рис. 3. 16. Инверсия графов
154
Блок-схемы автоматического регулирования обычно состав ляются в такой форме, когда имеются «точки суммирования» (где сходится некоторое число ветвей, по входит только одна ветвь) и «точки разветвления» (где одна входящая ветвь разде ляется на две или большее число выходящих путей).
Сигнал, связанный с «точкой суммирования», не сохраняется при инверсии, выполненной по способу сохранения ветвей. На рис. 3. 16, г, д, е показаны некоторые варианты инверсии графов.
2. Передаточные функции параметров движения самолета
- Представим уравнения (3.5) в виде направленного графа. На рис. 3. 17 показан общий граф прохождения сигналов корот копериодического движения самолета с учетом четырех топов упругих колебаний.
Рис. 3. 17. Граф прохождения сигналов продольного движения самолета с уче том четырех тонов упругих колебаний
На рис. 3. 18, а, б показан общий и существенный графы само лета с учетом двух тонов упругих колебаний.
На рис. 3. 18, в показан структурный граф четвертого порядка системы, характеризующий динамику движения самолета с уче том двух тонов упругих колебаний.
Для системы, предстайленной на рис. 3.18, характеристиче
ский полином запишется в виде |
|
А12Э4= (явР6+ аъръ+ а4 о4 + айръ+ а2р2 + ахр + а0) |
(3. 33) |
Значения коэффициентов а6; as; ...; ао приведены в приложе нии В. Числители передаточных функций ~ZPhAh определяются следующими выражениями.
155
Для пути прохождения сигнала в точку 1
0-»-1 или 6->-а
V /?u^?234= \пЬа(С\рь L С\р* -f С\рл-\- С\р2-|-
{в \р ^ -\- в \ р з Д- B lp 2Jr в \р ~ {-В о )-f- |
|
( O jp‘ + D lp ° + D l p t + D l p + Dl) + |
|
+ Щ ,ЛЕ\,*+Е1р>+ Е У ± E\l>+ E \,)]jr ■ |
(3.34) |
Рис. 3. 18. Общий, существенный и структурный графы параметров движения самолета с учетом двух тонов упругих колебаний
Для пути 0 —>3 или 8 —>q.A
V Л И ш 4 - [л«« ( c U 3 Д Cl»* д с \ р - \- с2)Д
Д я5*(Дтд д в1 о2д в 1 » - г в 1 ) +
|
|
3 ( d \ t f + Dip* Д |
L )\ o2 Д D \ p - f u l ) |
- r |
|
|
|
Д я ^ (ТтД3 Д E 32 p * + |
E'ln-EE'l)} — . |
|
(3. 35) |
|
|
|
pi |
|
|
Для |
пути |
0 — -7 или 8 —1>q± |
|
|
|
|
|
V Р ЫА \ \ М ^ [ п и { C t p 3 Д С2У Д C t p + C t ) + |
|
||
|
|
Д ( B t p 3 + в \ р 2 Д Д 5о4) д |
|
|
|
|
|
Д гао?3 ( А Д Д Е>\р2Д |
/Д рД /Д ) Д |
|
|
|
|
Д ^ Д ^ Д Д ^ з Д Д ^ Д Д ^ Д Д о ) ] |
— • |
(3.36) |
|
|
|
|
|
pi |
|
Для |
пути |
О —>2 или 8 —»)> |
|
|
|
^ P mAfm = [«5«(С24Д + С23Д Д Clp3+ С?рДС§) Д.
Д /г69 (В\р5Д В\р* Д В\р3Д В\рг-\- В\р-\-В20) Д |
|
Д K-bqi (7Д Д Д Е ^ з Р 3Д 7ДД Д D \ p - \ -D q) Д |
|
Дя-6qt {Eipi Д £ 3Д Д Ечр2Д Е\р Д /Го)] —г • |
(3. 37) |
р5 |
|
Коэффициенты С, В, D, Е определяются согласно |
алгорит |
мам, указанным в примере 1. |
|
Каждое добавление в систему нового тона упругих колебаний |
|
увеличивает порядок характеристического полинома и |
всех со |
ставляющих полиномов числителей передаточных функций на два.
На рис. 3. 19 показан структурный граф параметров |
движе |
ния самолета с учетом трех тонов упругих колебаний. |
|
Здесь |
|
Апт = (а8Р8ДачрчДа9рвДаърьДа4Д Д |
|
Д а з Д Д а з Д Д а ^ Д а ,,) - ^ . |
(3.38) |
р8 |
|
157
Числители передаточной функции (3.9) имеют следующую структуру.
О —* 1 или 8 —>а
V Я01Л?Д45= [пъЛсЬ1-!- ClP*+ .. . -г С\р + Со1) +
~ « а (Blp« + |
+ • • • + |
В\ о + Bl) + |
|
|
~ф Пц3(Мз »6Л- 0\г>ЪАг ■• - |
D\P~r Dо) —(— |
|
||
+ nb(u{Elps |
E |
lP5Jr- ■- + е 1р Ё £ о) + |
|
|
Н-я8?>( ^ в |
+ ^ |
6+ . - - + |
АГЬ+'^о)] у • |
(3.39) |
Для пути 0 —* 5 или 8 —* q5
/°15Л]2345= [«8а (Cfp5 4*' С4.04Д~ . . . -f- Ci о4Со) +
~г 14a {в\р° ~Ь В\р^“г • • • ~Ь в \р -)- Яо) -f
т" й»?з (-^5/ Л А Д ц4 Д-. ■■~\-D\р-\~ Do) Д-
-у Пъдл {Е \о° 4- Е\р*-\-. . |
. -j- Е\р-\- ^о) Д~ |
|
Д- fl'oqs {Кьр6~\- К1рЪ-\~ ■■.+ К Ь + К!)\ — • |
(3. 40) |
|
Аналогично определяются |
^^02^12345; ^ Я13Л12345; ^ |
12345- |
Графы и структура числителей передаточных функций будут иметь такой же вид и в случае, если источник возмущения — ветровые порывы WB(t). Изменяться будут только коэффи циенты эффективностей nWq. при соответствующих полиномах
числителей [см. (3.39) и (3.40)].
В табл. 3. 6 представлены полиномы числителя и знаменателя передаточных функций параметров СТС с учетом трех тонов упругих колебаний для скорости полета v = 220 м/с (Я = 1000 м).
158
1 > Р оА |
5 |
0 - 1
Вэ—а
2^02^12345 0 - 2
Вэ—D
2^*03^12345 0 - 3
Таблица 3.6
Полиномы передаточных функций
{0,37 (/7+25,6/6+3,8-103/5+7-104/4+5-106/3+35-106/2+15-108/+42-108) } — -
Р 1
= {0,37 (р + 2 ,85) (/>2+ 2 -0,035-22,6/1+22,62) (р 2+ 2■0,04-45,4/1+45,42) (/+ 2 -0 ,2 9 -3 9 ,1 /+ 3 8 ,12)}—
{5,6 (/i7+19,7/16+4-103/5+5.104/4+4,8-106/13+1,5-107/12—1,5-108/?+7,5- 108)} — |
— |
Р ‘ |
|
(5,6 (р +0,47) (/12+ 2-0,03-22,5 /Н + 2 ,52) (/i2+2-0,32-39,4р+-39,42) (/2 —2-0,03-45,1 р - |
45 , +) } — |
Р 1
%
{—35 (/6+ 60/5 + 4-103/4+17-104/3+4,3.106/2+10,3-106/+42-106) } — =, /6
= { -3 5 (/2 + 2 -0 ,3 1 -4 ,2 8 /+ 4 ,282) ( /2 + 2 -0.28-34/+342) (/2+2-0,008-44,9 /+ 4 4 ,92) ) — /6
*