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

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

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

П ередаточн ая функция реш аю щ его усилителя определится ф орм улой

 

 

 

 

 

 

( Р )

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

ЭК В И ВА Л Е Н Т Н Ы М КОЭФФИЦИЕНТО М УСИЛЕНИЯ
При исследовании САУ, содержащей нелинейные элементытипа зоны нечувствительности и насыщения, необходимо иметь простую и достаточно точ­ ную аппроксимацию нели­ нейных элементов.
Аппроксимация нели­ нейных элементов с по­ мощью эквивалентного ко­ эффициента усиления поз­ воляет так упростить мате­ матические уравнения не­ линейной системы, что легко удается определить частот­ ные характеристики многих замкнутых нелинейных си­ стем, произвести исследова­ ние их устойчивости частот­ ным методом, модифициро­ ванным методом корневого годографа и т. д.
Физический смысл экви­ валентного коэффициента усиления состоит в следую­ щем. Из статической харак­ теристики y'( y) нелинейного
элемента, показанной на
рис. 5.2, а, эквивалентный Рис. 5 .2. Статические характеристики нелинейного элемента привода yAy) 11 коэффициент усиления &(у) кривые изменения эквивалентного коэф­
определится как отношение фициента усиления k (у)
2С9

исследовании нелинейных многоконтурных систем высокого

по­

рядка получать обозримые результаты с минимальной затратой

времени.

5. 1

представлена схема

управления угловой

 

ско­

На рис.

 

ростью самолета с учетом двух тонов

упругих колебаний,

где

показаны характерные нелинейности типа «зона нечувствитель­

ности» и «ограничение».

них

являются

 

нелинейности,

представ­

Существенными из

 

ляющие собой статические скоростные

характеристики y'=f(y)

и e' = f(e)

соответственно сервопривода

автопилота и

бустера.

Эти.нелинейности будут учитываться при анализе и синтезе

контуров САУ.

нелинейных

систем

будет

использован

метод

При анализе

аппроксимации

нелинейных элементов

эквивалентным

коэффи­

циентом усиления.

 

 

 

 

 

 

 

 

5 . 1 . АППРОКСИМАЦИЯ

Х А РА К ТЕ РИ С Т И К

НЕЛИНЕЙНЫ Х

ЭЛЕМ ЕНТ

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