книги из ГПНТБ / Шевяков, Алексей Андреевич. Автоматика авиационных силовых установок учебник для авиационных вузов
.pdf228 Глава 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.
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+;