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

книги из ГПНТБ / Шевяков, Алексей Андреевич. Автоматика авиационных силовых установок учебник для авиационных вузов

.pdf
Скачиваний:
9
Добавлен:
29.10.2023
Размер:
11.67 Mб
Скачать
Фиг. 3. 68. Структурная схема регулято­ ра температуры газа.

228 Глава III. Системы автоматического управления ГТД

В качестве регулятора чисел оборотов возьмем регулятор с идеаль­ ным акселерометром, а в качестве регулятора температуры газов — регулятор температуры, схема которого приведена на фиг. 3. 27. Предполагаем, что имеется объединенное управление обоими регу­ ляторами, как это показано на той же фиг. 3. 27.

Рассмотрим указанную систему регулирования, когда она воз­ мущена перенастройкой регуляторов с помощью объединенного управления.

Выведем уравнение движения регулятора температуры газов. Структурная схема регулятора температуры показана на фиг. 3. 68, где указаны обозначения всех обобщенных координат.

Уравнение движения термо­ пары, как это рассматривается в курсе теории автоматическо­ го регулирования, возьмем в

виде апериодического

звена,

г. е.

 

 

(7 > + 1 ) ^ = 7 ^ .

(3.79)

Здесь Ti— постоянная

времени

термопары;

Кi— коэффициент

усиления;

Xi— ТЭДС

термо­

пары.

Из (3.79) видно, что чем меньше величина 7 л, тем мень­ ше динамические погрешности термопары, и, следовательно, в любой момент времени сигнал на выходе меньше отличается

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

менее 5—6 сек., а для апериодической—■1—2 сек.

магнитного

Усиление сигнала в первом каскаде,

состоящем из

и электронного усилителей, происходит

с некоторым

запаздыва­

нием за счет магнитного усилителя.

 

 

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

(Тгр + 1) ^ВЫХ1 = К Х ВЛ.

Уравнение же движения электронного усилителя принимаем как безынерционное звено в виде

Хвш2= КХк 2.

Фиг. 3. 69. Схема дифференцирующей ячейки.

2. Исследование систем автоматического управления

229

Поскольку ^ вьрц= ^ ВХ2. то общее уравнение усилителя можно запи­ сать так:

(Т2р + \ ) Х 2 = К2Х и

(3.80)

где К2 = КК.

Дифференцирование сигнала происходит с помощью R C — контура, схема которого показана на фиг. 3. 69.

Уравнение движения этого контура, полученное из уравнения, связывающего напряжения на входе и на выходе

R\Ri

CL + 12

R\Ri c i x +

R\ + R^

‘ ‘

#i + #2

преобразуется к такому виду:

(Tnp + \ ) X 3= (T up + R ) X 2,

где

(3.81)

 

п __ /?2

 

RiR2C t

 

Tu- R\~\~R2

R=

t'iS

R\-\-R2

Ri

R\ + #2 lи

V D

Подбирая значения сопротив­ лений и емкости, можно получить различные значения параметров Т\\ и R. Следует иметь в виду, что с уменьшением значения R точ­ ность дифференцирования увели­ чивается.

Следовательно, рассматриваемый дифференцирующий контур дает сигнал на выходе, зависящий от входного сигнала и первой производной от него.

Усиление сигнала во втором усилителе происходит по аналогич­ ному с (3. 80) уравнению, т. е.

(Т3Р + \ ) Х , = К , Х 2.

(3.82)

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

 

- P m- P » p - Q + Pp=0>

(3.83)

где РПр —усилие,

развиваемое пружиной;

 

Q — вес якоря и золотничка;

силами;

Рр—усилие,

развиваемое электромагнитными

Рт— инерционное усилие, обусловленное массой и ускорением подвижных частей.

230

Глава III. Системы автоматического управления ГТД

При малых отклонениях якоря (отклонения действительно бе­ рутся малыми) можно выбрать такое его положение, при котором величина электромагнитных сил Рр практически не зависит от по­ ложения якоря. В этом случае сама величина электромагнитных сил определяется таким выражением:

Яр = const г2,

где I — ток в катушке соленоида.

Положим, что г = г0 + А(. Тогда при условии (Дг)2яа0 получим

Яр = const (г2 + 2/0Дг).

Выражение для А/ определим, применяя уравнение электриче­ ского равновесия цепи катушки соленоида с учетом того, что потокосцепление <|) является функцией как тока I, так и положения яко­ ря соленоида /, т. е.

и= « -f- ф = гг+ — / + — — .

1 т

d l

d i dt

Здесь г — омическое сопротивление катушки соленоида; и — напря­ жение, снимаемое с усилителя.

Учитывая, что — = ai, и переходя к малым приращениям, полу-

d i

чим такое выражение:

Au = r\i-\- aiaAl + aAil-\- L — , dt

где L — d^/di — индуктивность катушки соленоида, a = const.

Обычно скорость движения якоря относительно мала; поэтому

членом a A i l можно пренебречь как величиной второго порядка ма­ лости. Вводя оператор дифференцирования, из полученного выра­ жения найдем значение Аг, т. е.

^ _ Аи — a i 0p \ i r + Lp

Подставляя полученное в выражение для Рр, получим

Я = const г2 + const 2 г' 0 Д“ ~

а1°р .

r +

L p

Подставляя это выражение для Яр в (3.83) и считая, что Япр=

=

Я0 4-Спрд/;

Pm = mp1Kl, после преобразования и сокращения на

члены, учитывающие статическое равновесие, получим

/

mL з

m

2

Cnp L — const-f2a/p

Л A t_const 4-гг02г0 Дгг

\ б"прг

Cnp

t7npr

/ /q

I0Cnpr

Uq

2. Исследование систем автоматического управления

231

Вводя обозначения

Г

L

п .г

 

тп

'

т _ const + 2 a i\

JK—

r :

/м ~

 

r

 

r

 

Г

 

 

C n p

 

 

^ n p '

Д, = const_+2«oip .

у

_

Д/

Дм

4

lnCmrО'-’пр'

 

5

/0

«о

получим окончательно

г2 „з I

-г2

(3.84)

+

rip + (7"K+ 74)p + 1 ] X5= K 4X 4.

Обычно 7), очень мало, поэтому полагая 7'Kss0, получим

(Tip2+ Т4р + 1) Х 5= К 4Х 4.

Для более грубых расчетов можно принимать и Т^~0, тогда получим

(Т4р + \ ) Х 5= К 4Х4.

(3.85)

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

тивного сопла.

Следовательно, уравнение движения гидромотора примет вид

p X F = KclX 5.

(3.86)

Из фиг. 2. 10 видно, что в некоторой области работы двигателя для получения экономичных режимов целесообразно изменять ре­ жим его работы за счет температуры газов, сохраняя неизменным значение степени повышения давления в компрессоре *-*. Будем

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

оставаться неизменным, что приблизительно

(при пологих харак­

теристиках компрессора) и будет отвечать

условию тг* =const.

Тогда вместо (3. 79) примем

 

= 7Т/)+ 1 (Х° - ^ 4 ) -

(3-87)

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

__________ К-!К2К3К4Ки(ТиР + Я)__________

(Ттр + 1 ) СТпР + 1 ) (Т2Р + 1 ) (ТзР + 1 ) (Т4р + 1

_________K2K3K4Kci(TuP + R) ______х о

Р (ТпР + 1) (Т2Р + 1) (ТзР + 1) (Т4р + 1)

232 Глава III. Системы автоматического управления ГТД

или, обозначив

ф

КтКцКяК^Ку (TUP + R)___________ .

 

1

( U p + 1) (Тир + 1) ( U p + 1) (Г3/>+1) ( Г 4р + \)р'

 

ф

K jK sK tK aiT u p + R)

 

 

2 Р (ТиР +1) (Up +1 )(Up +1) (Up + 1) ’

 

получим

 

 

 

- Ф 1Аг7-4 + Хр = - Ф 2Х°.

(3.88)

Передаточная функция регулятора числа оборотов с идеальным акселерометром, как это следует из рассмотренного выше, будет такой:

Ф

= +

— •

(3-89)

х п

 

Р

 

Следовательно, для решения поставленной задачи необходимо решить совместно систему (2.30), в которой четвертое уравнение заменить (2.50), с учетом (3.88) и (3.89).

Для дальнейшего упрощения систему уравнений, описывающую движение объекта регулирования [система (2.30) без четвертого уравнения и (2.50)], сведем к двум уравнениям, включающим координаты Х„ и XTi, т. е.

(Тхр + mQ) Х п+ mxXTi m2X0т = 0;

ЩХп — Щ Х ц + ЩХатХ Р= 0.

Здесь коэффициенты пц; mx ...ms выражаются через коэффициенты основной системы.

Итак, исходную систему уравнений принимаем в следующем виде:

(Т\р + щ ) х п + щ Х та— щ х от= 0; m3X n- m 4XTi-\-m5XGT- X F= 0;

- Ф 1^ Г4 + Л > = - Ф #Х°;

 

<ЬаХ я + Х а = 0 .

. I

Знаки перед Фх, Фг и Ф„ взяты с учетом характера регулирования, т. е. необходимости уменьшения расхода топлива при увеличении числа оборотов и необходимости увеличения проходного сечения реактивного сопла при увеличении температуры газов. В соответ­ ствии с этой системой уравнений укрупненная структурная схема

2. Исследование систем автоматического управления

233-

Фиг. 3. 70. Укрупненная структурная схема.

системы будет такой, как показано на фиг. 3. 70. Решение системы представим в таком виде:

Д*„=АЛХ° и ЬХТ4= Дт-д0,

Т\Р + т0

 

т х

0 /и2

Щ

 

- т 4

- -1

т 5

0

 

-Ф ,

1

0

®„

 

0

0

1

0

 

т ,

0

т 2

0

-

т 4

1

Щ

-Ф а

-

Ф.

1

0

0

 

0

0

1

Т,р + т0

 

0

0 — т 2

Щ

 

0

- 1

т5

0

 

-®а

1

0

Ф

 

0

0

1

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

234 Глава III. Системы автоматического управления ГТД

л магнитные усилители главным образом выполняют функции пре­ образователей тока.

Таким образом, вместо (3.80) и (3.82) получим

 

Х2 = К2Хр,

Х 4= К 3Х 3.

(3.92)

После раскрытия определителей и замены Ф^ Фг и Ф„ их значе­

ниями с учетом (3. 92) получим

 

 

{(7> + 1 ) (Тпр + 1 ) (Т4р + \)р [Кс(АТ, + КЛр) (mi/n5-

m2m4) -

— (?\р + mQ)pm4- t n ^ p ] - K TK2K3K4Kcl (Tup + R) [Kc {K{-f Kap)m2+

+ (Tlp + m0)p]} X n = -

туК2К3К4КЛ (Tnp + R)pW.

{(Trp + l ) ( T np + \ ) ( T 4p + \ ) p [ K c(Ki + Kj>) (OTjOTg — m2rn4) -

— (ТхР + т0)рт4 — т3т1р]—КТК2К3К4Ксх (Tnp Jr R )X

X [KAKi + KP) Щ + ^ р + т^р)} X T4~

= - {K2K3K4KZl {TuP + R) [(TlP + m0)p + m2Kc{Kx + Kap)]} X».

После преобразования получим следующие уравнения:

(а0рь -f агръ+ а2рк+ а3р3 + а4р2-f а5р + а6) Х п = (Ь1р + Ь2)рХ°. (3. 93)

{а0р6 + а ^ 5 + а2р* + а3р3 + а4р2 + аБр + а6) Х Т4 =

 

= (b3p3 + b4p2 + bsp + b6)W ,

(3.94)

где

QqаТгт4,

а1—ag ЬТ,/я4;

a2= b g + aKJipn — с'Грп4; a3 = cg + b K iKcm — KT1Tn;

а4 —Kg1+ сКхКрп\ as= Kg2,

a6= K K cK1m2R; b4= т^КТхх,

b2 — mxKR\

b3 = KTxTxx\

b^K fJxR+ Txxma + m^K^Jxxy, b5= K [m2Kc{KtR+ TxxKx) + Rm0]; be= K m 2KcKxR;

2. Исследование систем автоматического управления

235

а — ^т^4^11>

Ь = ТтТ4+ Т пТ4+ ТпТ4\

с ~ ^4+ Тп + ^-п

g^=KcKdm — >п0пг4 - т3тх-,

Si = Я Д .Т пот2 ~ ТпЩ +

£2 = К^КДт 2 — КСК4Тпт2 + Rm0;

т= т4тъ — т2т4,

/С=/<т/<2/<3ЛУ<с1;

К = К 2К3К4КС].

Из полученных уравнений видно, что при рассматриваемом воз­ мущении Х°, обусловленном перенастройкой регулятора темпера­ туры газов, для устойчивой системы после затухания переходных составляющих число оборотов остается прежним, а температура газа устанавливается заданной.

Действительно, предполагая, что система устойчивая, при еди­ ничном возмущении найдем конечное значение параметров Хп и XTi.

X (oo) = lim

•S2(6,.S + 62)

 

= 0 ;

.5 (a0Se + a 1S5+ . .

. +

 

s-о

a6)

(oo) =

lim I ■

(6-yS3 + 6.},S2 • ;- 6;;,S +

6,l)^

*6.

 

—о I

-S (йо-S6+«1-S5+

. .

й6)- +

“s

или после подстановки значений коэффициентов получим

 

 

ХТ4(оо) = Х°.

 

 

Предполагая,

что переходные процессы по Хп знакопостоянные,

а по X Ti процессы монотонные, определим величины интегральных погрешностей.

Предварительно определим начальные условия. Для уравнения

(3. 93) получим

 

 

 

(6,S + 62)S

 

 

 

 

X (0) = lim

5

 

 

 

= 0 ,

 

 

ао^е+

+

. . , + a 6)S

 

 

S~+ОО

 

 

и далее А'1(0) = АГ1Г(0) =

Л"111 (0) = 0;

 

 

 

 

AIV(0) = lim

5s —

 

(6!S+ft2)S

 

*L;

 

 

s —~

 

 

 

+ a i^5+ • • - + а б) ^

йо

 

A'v (0) = lim

S 6 ■

(6jS +

62) 5

.

^0)

62 _

M i

aiS5 +

• - + ae)s

 

«0

a0

S —►*»

(aos6 +

 

236

Глава III. Системы автоматического управления ГТД

Учитывая, что общее выражение для интегральной погреш­

ности

имеет вид

J = ^ [ a 0X ^ ( 0 ) + aiX ^ ( 0 ) + . . . + а „ _ 1.У(0)1, (3.94а)

то в соответствии с этим выражением для J получим

=_ Л

 

аб

й0 \ а 6

 

/

 

Заменяя коэффициенты

их выражениями, получим

 

J n= m 1 — ^

( К

1

\

1 1

(3. 95)

 

 

K-rKcKimi

 

 

 

 

 

Начальные условия для уравнения (3.94), с учетом переноса начала координат будут такими:

X (0) =

lim

s

 

b3S3+ b4S*+bsS + b 6

 

 

s^o

(a0S6-1-aiS5+ . . . + ae)S

 

X I (0) =

Нт

S2-

b4S2~b bsS

bg

■0 =

 

S-*ao

 

( a „S 6 + a i S 5 + . . .

+ a 6 ) S ~

 

и далее

 

 

 

X й (0) = 0;

 

 

 

 

 

 

 

 

(0) =

limLS4

 

b,S3-[- b4S%-j- b§S -(- bg

- 0

 

|

 

(a0S6 -|- a,S5-j- .

+ a6) s

 

X lv (0) = lim ■S5

 

-f- b4S2-j- b5S -f- bs

-5ЛГШ(0)

S->°о

(uqS^ -)- «iSS-)- . . .-|-ag)S

 

 

be

«6 ’

0

_ .

«о

b$ ®\ b;, '

О 9 a0 ag

(0) - lim ^6

+

b /jS 2 -f- b5S -f- b e

_ -S*Xlll(0) — S X lv(0)

 

( aoS6 +

f l l S 5 + .

. , +

e 6) s ‘

 

 

 

 

b3a2+

aifr4

aib3

 

 

«0

 

a0

 

2

 

 

 

 

clq

 

В соответствии с найденными начальными условиями напишем

выражение для интегральной погрешности

 

 

 

J T4— '

2a]b3

а5Ь6

 

 

 

ао

2

 

 

 

 

ав

«6

 

Заменяя коэффициенты их выражениями, получим

 

 

 

 

 

а\К

 

Тц

(3.96)

 

J TA=

2 T,Tn T\T\m\

KTR

Из

полученных выражений для интегральных ошибок

для Л,

и ] Ti

следует, что для уменьшения погрешностей необходимо уве­

2. Исследование систем автоматического управления

237

личивать коэффициент усиления регулятора числа оборотов (KiKc)

и уменьшать коэффициент усиления К регулятора температуры газов.

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

Чтобы составить представление о характере переходных про­ цессов при изменении некоторых параметров системы, разберем конкретный пример.

Пример

З а д а н и е . Рассчитать систему регулирования одновалыного ТРД с регу­ лируемым реактивным соплом (определить основные параметры) при работе на максимальном режиме так, чтобы при перенастройке регулятора температуры газа за турбиной характер изменения температуры был бы близким к монотон­ ному, а изменение чисел оборотов не превышало бы 5"/о за минимальное время регулирования.

И с х о д н ы е д а н н ы е . Регулятор числа оборотов принять с идеальным акселерометром, рассчитанным в прошлом примере. Данные по двигателю при­ нять из первого примера. Регулятор температуры принять по схеме фиг. 3. 27.

Инерционность магнитных усилителей не учитывать.

 

исполь­

Р е ш е н и е . Определим коэффициенты уравнений (3.90), для чего

зуем (2. 30), в котором четвертое уравнение заменим

(2. 50), и числовые значения

коэффициентов, полученных в примере 1 на стр.

85 и в примере 3 на стр.

88. По­

лучим следующие коэффициенты:

 

 

 

7-4 = 0,83; тир = — 4,32; ml = — 4,0;

m2 = — 1,78;

 

m3 = 5,66; m4 = — 3,84;

ms = — 2,25,

 

Подставив найденные коэффициенты, получим такие уравнения движения объекта регулирования:

(0,83р — 4,32) Х п — 4,0Хт4+ 1 ,78Х0т = 0;

5,667+ + 3,8 IXи —Х Р — 2,25 Х„, = 0.

В соответствии с (3.93), .(3.94) и 7(с=2,3; 7(i=8,25; 7(а=3,4 из прошлого примера, получим такие значения коэффициентов уравнений (3.93) и (3.94):

а0 = 3, WTrTu T4; а, = 3,19 (7ТГ4 + Тп Г4+ 7+7+ + 23ГТ7 П74;

*2 = 3,19 (7-4 + 7+ + 7+ + 23 (ТтТ4+ Тп Т4 + 7Т7+) + 41,37++7+

*3 = 3,19 + 23 (7-4 + 7+ + 7+ + 41,3 (ТтТ4+ Тп Т4+ 7-т7п ) -0,837(7+ ;

а4 = 23,0 + 41,3 (74 + 7+ + 7+ + 18,327(7+ - 0,83/С/?;

*5 = 41,3 + 18,327(7?+ 3 4 /0 + ; а 6 = 347(7?; 61 = 4/(7+;

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