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

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

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

Рис. 5. 14. Определение коэффициента усиления в нелиней­ ной следящей системе

Рис. 5. 15. Определение величины входного сигнала k(y) нелинейного элемента

0,2 0,6 1

Z

у,град

Рис. 5. 16. Статическая характеристика нелинейного элемента

220

Для построения амплитудно-частотных логарифмических ха­ рактеристик (5. 15) достаточно (при определенных значениях 61= const и разных частотах) последовательно вычесть из амплитудно-частотных характеристик замкнутой системы

—(До; &(у))

амплитуды частотной

характеристики

6]

S^const

функции

k{y).

|иД(/ч>)| интегрирующего звена и амплитуды

Следовательно, значения k(y) и у стали известны. Тогда

стати­

ческая

характеристика нелинейного элемента определится из со­

отношения

(5.16)

 

y' = k(y)y-

На

рис. 5. 13—5. 16 показаны кривые графического

решения

указанной задачи.

Описанные выше способы определения частотных свойств нелинейных следящих систем дают возможность исследовать динамические процессы в нелинейных системах автоматического

управления.

 

5 . 3 . В О ЗМ О Ж Н Ы Е

СОСТОЯНИЯ САУ, ОБУСЛОВЛЕННЫ Е

НЕЛИНЕЙНЫ М И

Х А РА К Т Е Р И С Т И К А М И ПРИВОДОВ

Определим влияние параметров линейных и нелинейных эле­ ментов приводов, применяемых для построения внутренних вы­ сокочастотных контуров стабилизации. Известно, что от совер­ шенства этих контуров зависит возможность построения более точных и совершенных систем управления центром масс.

Рассмотрим систему стабилизации угловой скорости само­ лета— контур демпфирования. Его структурная схема показана на рис. 5. 1. Здесь учитываются два тона изгибных колебаний, и в качестве датчиков используются угловой акселерометр и дат­ чик угловой скорости.

Рассмотрим вначале случай, когда самолет жесткий, а для демпфирования угловых движений используется информация

только датчика угловой

скорости. При исследовании влияния

 

нелинейности привода воспользуемся модифицированным мето­

дом' корневого гбдографа.

 

 

 

Характеристическое

уравнение нелинейной гармонической

 

линеаризованной системы можно представить в виде

 

 

 

А (р) = А 2(р;

k(y); k ) + A l(p; k(y); k) =0,

(5.17)

 

или

Л2 {P\ k (у); k)

(5.18)

-

 

 

A (p; k (y); k)

Здесь варьируемыми параметрами являются коэффициент усиления привода k(y) и коэффициент усиления прямой цепи контура k. Для устойчивости этой нелинейной системы, как и для линейного случая, все нули лу полинома четной степени

221

A2(p\ k(y)- k) и полюсы Hi полинома нечетной степени Ai (p\ /г (у); k) должны лежать на мнимой оси плоскости р, а их

абсолютные

величины — удовлетворять соотношению

(4,50).

 

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

авто­

пилота, бустера и самолета будут

 

 

 

 

 

 

 

 

 

 

 

 

\

 

 

 

S1

Т{Гг,

 

Т 1 -[- То

 

 

 

 

 

 

 

к (у)'

к(1) 2 + k (у)

 

 

 

 

 

Лр У

 

 

 

 

 

(5.

19)

 

 

 

 

 

 

 

 

 

 

к (е)

 

 

 

 

 

 

 

 

 

 

02Л12

1

к^ (ТvP +

1)

 

 

 

 

 

H2

Р

 

v + i

 

 

 

 

 

 

 

7’i ^ + 24

)

 

 

 

 

 

 

 

 

 

 

 

Здесь

соблюдается

условие 6 Ша х * ^ ф т а х , при котором

бустер

работает

в

режиме

максимальных

коэффициентов

усиления

fe(e) и минимальных фазовых сдвигов. И в этом случае постоян­ ная времени бустера может быть аппроксимирована постоянной

величиной

—^—

=Тз. Тогда характеристическое уравнение зам-

 

к{г)

 

 

 

 

 

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

 

Аъ{р'-Ь(уУ>к) =

аврв +

аьр5+ . . .

а^р-\-а0=

 

= Т{ГгТ3Т2. р*+ \т \ (7\7’а+ 7’1Г, +

7’,7’,) +

 

+ 2 ^ 7 ’17’я7’3] p8 +

[7’J ( r i -f-7’a+ 7 's) +

 

+ 2^

\ Т{Г2+ ТХТз+ Т3Т3) ; О

Д

р* +

+ [Т2.(1 + T3k{у)) + 2

(7\ + Г2 + Т3) + Т{Г2+ Т{Г3 + Т2Т33 +

+

[ ( ^ + 2 ^ 7 ’з7’3| А(у) +

2^7’з.+ 7'1-(-7'2-|-7’з| р2+

 

+

+

 

+

1 ]

р-\-

 

 

+ (1 + V )* (Y ) = 0.

 

(5.20)

Для исследований этой нелинейной системы определим гра­ ницы устойчивости в области параметров k и k(y) на двух ха­ рактерных режимах полета: № 1 (режим больших скоростных напоров) и № 2 (режим малых скоростных напоров). [См. при­ ложение А, табл. А. 3.]

222

Функции для определения границ устойчивости будут для режима № 1

Л2 (х; &(у); £)^=х! + 9150л:4-г [885000 + 130000£(у)] х +

+ [7700000 4- 5040000^] k (у) = 0;

T+y+lY K к)

//„ /г [805

13,9* (y)] ^ +

 

 

+ 43500 + (2830 + 30800Ш (у) = 0

(5.21)

для режима № 2

 

 

 

 

А2[х\ k (y); k) = х3 -|- 8300л:4 + [9030012600* (у)] л: +

 

 

+ [745000+ 104000*] £(y) = 0;

 

 

Ai [у; k (у); к)

у{) {у2[736 +

14,5* (у)] у +

 

 

+ 4250 + 300* {y)-t-4250*-* (у)} = 0 .

(5.221

Определим границы устойчивости для k(y) =4,

5, 9 и 16 и для

fe = var. Результаты решения уравнений

(5.21) и

(5.22) приве­

дены на рис. 5.17. Здесь граница устойчивости изображена

кри­

вой k(y)=f(k),

которая проведена по точкам пересечения

кри­

вых движений

корней у \(k) и х2(/г) для различных k (у) = const.

Кривые движения других корней здесь не показаны. Левее кри­ вой k (у) =f(k) находится область устойчивости.

Из кривых рис. 5. 17, а видно, что меньшим значениям k(y) соответствуют большие значения k. Для этого случая

dfe(V)

дк ^

При достижении граничных значений k> ka в системе нач­ нется режим автоколебаний. Частота автоколебаний будет зави­ сеть от величины коэффициента усиления в прямой цепи k кон­

тура

стабилизации. Точки

аи а2, а3, а4 пересечения

кривых

y\(k\

k(y)) и x2{k\ k(y)) определят частоту автоколебаний. Так,

для режима № 1

 

 

<оа =

]/"+ — \/~ х 2 = У 320 ~

18 рад/с при &= 0,28 с и * (у) =

18 с-1;

соа — У 130 ~ 11 рад/с при £ = 0,35 с и k (у) = 4,5 с- ! .

Кривая k(y)=f(k) крутопадающая. Здесь частота автоколе­ баний весьма чувствительна к изменениям коэффициента усиле­ ния k. На режиме № 2 (см. рис. 5. 17, б) кривая k(y) = /(£) также крутопадающая, но для случая

дк

* 0

дк

223

при k~>k&в системе возникают колебания возрастающей ампли­ туды и уменьшающейся частоты. Система становится неустойчи­ вой в большом. Это явление недопустимо с точки зрения безо­ пасности полетов.

Рис. 5. 17. Кривые границ устойчивости замкнутой системы стабилизации угловой скорости без кор­ ректирующего устройства

Как видно, модифицированный медод

корневого годографа

и способ построения кривых определения

границ устойчивости

k(y) =f(k) по точкам пересечения годографов основных корней Уг и Xi в области двух параметров k(y) и k позволяют при малом объеме вычислений исследовать динамику и возможные состоя­ ния нелинейной системы стабилизации и управления самолета. Данный метод дает обозримые -результаты при исследовании

22 4

систем автоматического управления высокого порядка с не­ сколькими нелинейностями.

При наличии в системе нелинейных элементов с неоднознач­ ными характеристиками типа «гистерезис», «люфт» или элемен­ тов чистого или распределенного запаздывания эти нелинейные элементы аппроксимируются эквивалентным коэффициентом усиления вида k ( y ) = g ( у)—jb(у) и исследование состояний не­ линейной системы осуществляется аналогично.

Рассмотрим состояние замкнутой системы при наличии нелинейных элементов в последовательно включенных контурах.

В качестве существенной нелинейности, кроме скоростной

характеристики привода автопилота у'(у),

будем

учитывать

нелинейную

статическую характеристику

бустера

е'(е) (см.

рис. 5.1, 5.5,о и 5.9).

в точке

А системы

Частотную

характеристику разомкнутой

жесткого самолета (см. рис. 5.1) при использовании информа­ ции от ДУСа можно записать в виде выражения

W{j^) = k тгт2

 

\ + Т 2

1

 

X

 

 

к (у)

 

k(y)

Р k (у) Р \

4-1

 

 

1

 

 

к^ (

T v p +

1 )

 

(5. 23)

 

X к(в)р + 1

 

 

 

 

 

 

7 xI

2+

^

V

+ 1

'

 

 

 

Важными параметрами статических характеристик нелиней­

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

являются

 

максимально

возможные скорости

аВТОПИЛОТа ^ m a x ^ Y m a x И

буСТвра

' - Р т а х = =

£ т а х -

ЕСЛИ бщ ах^фтах»

то п р и н а л и ч и и

р е з к и х в о з м у щ е н и й ,

п о д а в а е м ы х н а

си стем у , и л и

п р и б о л ь ш и х п о м е х а х б у с те р б у д е т в ы х о д и т ь н а р е ж и м н а с ы щ е ­

н и я по ск о р о сти , т. е. н а

р е ж и м

м а л ы х к о э ф ф и ц и е н т о в у с и л е н и я

к(г), п р и

к о т о р ы х

и м ею т м ест о

б о л ь ш и е

ф а з о в ы е з а п а з д ы в а н и я

(см . рис.

5.10). В т а к о й

си с т е м е м о ж е т

б ы ть

н е с к о л ь к о а в т о ­

к о л е б а т е л ь н ы х р е ж и м о в .

 

 

 

 

 

 

На рис. 5. 18 показаны амплитудно-фазовые

характеристики

системы

(5. 23) для режима полета при больших скоростных на­

порах, когда

 

 

 

 

 

 

 

 

»

I л

1,5(0,43р + 1)

 

 

 

{рУ-

1

 

 

1

 

 

 

V

 

 

 

р + 1

 

 

 

-------------/ > 2

+

2 - 0 , 4 —

Ю

 

 

 

 

102-

г

 

 

 

Кривая 1 фазовой характеристики соответствует случаю, когда бтах^фтах и бустер не выходит в режим насыщения по скоро­

сти. Кривые 2 и 3 получены при условии, когда бтах>фтахИз кривой 2 следует, что в данном случае возможны автоколеба­ ния с частотами сон и coB, а из кривой 3 — что в системе воз­

8

3819

225

можны автоколебания только низкой частоты сон и большой

амплитуды.

Поэтому в системах с несколькими последовательно вклю­ ченными контурами, содержащими нелинейные элементы, пара­ метры нелинейных (или линейных) элементов должны быть вы-

/

г

3

5

Ю

20 30 со,раО/с

Рис. 5, 18. Частотные характеристики системы при нали­

 

чии нелинейного внутреннего контура

браны так, чтобы ни один из контуров не выходил на режим

больших фазовых запаздываний.

ранее,

должно

В рассматриваемом' случае, как указывалось

быть соблю дено уСЛОВИе ф ш ах ^ бщах-

 

 

5 . 4 . СИНТЕЗ КОРРЕКТИРУЮ Щ И Х УСТРОЙСТВ

ДЛЯ

УСТРАНЕНИЯ

НЕУСТОЙЧИВЫ Х СОСТОЯНИЙ НЕЛИНЕЙНОЙ СИСТЕМЫ СТАБИЛИЗАЦИИ

Наличие инерционных звеньев и нелинейностей приводов приводит к тому, что на некоторых режимах полета в замкну­ той системе угловой скорости могут возникнуть расходящиеся (до определенных пределов) колебания с уменьшающейся ча­ стотой, т. е. система становится неустойчивой в большом. Упру­ гие колебания, измеряемые датчиками, усугубляют этот процесс. Кроме того, при неудачном размещении датчиков область устой­ чивости системы стабилизации самолета с учетом нежесткости

226

его конструкции значительно сужается, а при больших возму­ щающих воздействиях система может стать даже неустойчивой вследствие явлений «захвата» (образования положительной обратной связи).

Следовательно, система стабилизации угловой скорости дол­ жна проектироваться так, чтобы было исключено вредное влия­ ние нелинейностей и инерционности приводов, а также влияние помех от упругих колебаний при широком изменении динамиче­ ских свойств самолета.

Привод, состоящий из последовательного соединения серво­ привода и бустера, имеет полосу пропускания 2—2,5 Гц с фазо­ вым запаздыванием выходного сигнала не менее 90° (см.

рис. 5. 12).

Нелинейная статическая скоростная характеристика приво­ дов типа «зона нечувствительности» и «насыщение» сущест­ венно искажает фазовые характеристики замкнутой следящей системы и делает эти характеристики неоднозначными. Поэтому необходимую точность стабилизации и качество регулирования на некоторых режимах полета получить не удается. Расширение же полосы пропускания сервопривода приводит к уменьшению запаса устойчивости' замкнутого контура из-за влияния упругих колебаний. Поэтому целесообразен другой путь создания кон­ тура стабилизации угловой скорости самолета как жесткого тела.

При наличии привода с минимально необходимой полосой пропускания требуется сформировать управляющий вектор с не­ обходимым фазовым опережением, корректируемый по частоте. Для этого должен быть выбран определенный состав датчиков и разработаны корректирующие устройства, которые должны:

1)обеспечивать фазовое опережение сигналов обратной связи

вдиапазоне частот короткопериодического движения самолета. Это даст возможность скомпенсировать фазовые запаздывания привода;

2)существенно уменьшать среднюю квадратическую инте­ гральную оценку параметра ■в’ или if);

3)обеспечивать расширение границ устойчивости по отноше­ нию к коэффициенту усиления k\

4)исключать возможность появления неустойчивости си­ стемы в большом (собственная частота системы на всех режи­ мах полета должна быть в пределах 3—3,5 Гц);

5)не усиливать помехи от упругих колебаний (должна быть обеспечена амплитудная стабилизация тонов упругих колебаний с запасом 8—10 д Б );

6)быть простыми в реализации.

Рассмотрим следящую систему автопилота (или демпфера). На рис. 5. 19 показан граф прохождения сигнала в следящей си­

8*

227

ст ем е п р и в о д а а в т о п и л о т а и и н в ер сн ы й г р а ф . К о э ф ф и ц и е н т п е р е ­ д ач и и н в ер сн о го г р а ф а б у д ет

Тоо'= 1 4~ 1Т2р1-{~(7\ХТ2) р-f 1] —■1 p = -L(p)'

(5.24)

д у с ____________________________________________ 8

Д У У

Из^ соотношения (5.24) следует, что сигнал корректирующего устройства Т00, является суммой сигналов двух путей. Для ком­ пенсации запаздывания привода в контуре стабилизации угло-

 

/

/

 

]_

 

O'

T,P+1

Т2р+1

 

Р

и

*1

 

 

г

f

с?

а)

 

 

 

II

 

 

 

\ &

 

 

 

Т1Р+1

7

W Р

 

 

\ ггР+1

0

4

"

l Т

1 м7' 71

г/

 

____________ 1______L

 

.

6)

' 6 (т1р2+г1тзр+1)(Цр+1)

Рис. 5. 19. Основной и инверсный графы сигналов привода автопилота

вой скорости целесообразно использовать два датчика: датчик угловой скорости (ДУС) и датчик углового ускорения (ДУУ), причем для компенсации влияния нелинейной статической ха­ рактеристики привода коэффициент усиления ДУУ должен быть

нелинейны м -^-. Передаточную функцию (5.24) можно реали­

зовать в виде корректирующего устройства в цепи ДУСа, когда датчик углового ускорения отсутствует.

В ы б е р е м о п е р а т о р 0 с п е р е д а т о ч н о й ф у н к ц и ей

{ Т \ р 2 + + l ) ( Т 4р + 1 ) а 3рЗ + д 2 р 2 + a ip + j

Т о гд а п е р е д а т о ч н а я ф у н к ц и я к о р р е к т и р у ю щ е г о у с т р о й с т в а з а п и ш е т с я в в и д е

7°0' = 7W6 = f 1 +

[Т1г 2^ +

(7\ + Т 2) р +

1] - Р - )

х

1

 

 

k(y) 1

 

1

 

 

 

( Т 3 Р2

+ 2£зТ3р +

l) (ТАр + 1)

(5.25

где

^ 4 < 7 'з < 1 / г 17'2;

 

228

N= -----= —------нелинейный элемент, необходимый в основном

My) y '

для компенсации эффекта зоны нечувствительности.

Корректирующее устройство 7оо' можно представить пере­ даточной функцией вида:

N T { T 2p i 4- N (Ту + T 2) p 2+ N p + 1

 

b ^ p 3 + b ^ p 2 + b ^ p + {

(5. 26)

а-зр3 + а 2р 2 + а хр - f - 1

 

р 3+ \ Ч 2 р 2 +

\Ьу р + 1

 

 

Коэффициенты

а3; а2\

fli

выбираются

из

соотношений

а3='к3Ьз\ а2

ау = \Ь?,гце

с точки зрения реализации

X за­

дается в пределах 0,5—0,6.

 

 

у

которого

пара­

Для

привода контура угловой скорости,

метры

Г] = 0,01 с,

72= 0,04

с и My) (в режиме до насыщения)

изменяется в пределах Му) = 4—25 оператор 0 может быть реа­ лизован передаточной функцией вида

________________1_________________

( “ 7 2 + 2-0,5 + Р + 1)

Передаточная

функция

(5.26)

может

быть

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

вана функцией более низкого порядка

 

 

 

Ь?р* + Ь%р2 + Ь?р+. 1 _

В ^ р 2+ В^р + \

 

/?3 +

Н“

Р +

1

A<ip^ +

р + 1

 

Пользуясь методом приближения функций [6], запишем

 

 

 

-f ^) (Л 7 а + А ,Р - И ) ~

( В 2 P 2 J r B i

Р -\~

О (а з 7 3 + а 2 7 2_)_ а 1 /7 +

!)•

Приравнивая последовательно коэффициенты при рп (п=\, 2,

3 и 4), получим уравнения для определения

В2 ;

Ву\ А2 и А у-.

В ? = А у — ау + №

 

 

 

 

В2 =

■jN

N .

 

 

,JV

 

dyВ\ -f- b\ Ay

-(- A2d2-\- b

 

Ay:

a2 £ f +

 

hN

 

ДЗ

hN

,N

hN

 

hN

 

uN

 

 

hN

 

«3 D ' V

I

1*2 D ‘

b"

л

 

 

 

3

 

 

 

 

 

 

bN

Л1*

 

 

 

 

 

 

°2

 

 

 

Решая эти уравнения, находим £г, coi.h g2; сог для различных значений k (у) =4, 18 и 25 с-1.

8

3819

229

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