Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Кашин Г.М. Автоматическое управление продольным движением упругого самолета

.pdf
Скачиваний:
9
Добавлен:
24.10.2023
Размер:
10.78 Mб
Скачать

Рис. 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

*

Соседние файлы в папке книги из ГПНТБ