книги из ГПНТБ / Шевяков, Алексей Андреевич. Автоматика авиационных силовых установок учебник для авиационных вузов
.pdf28 |
Глава I. Поршневые двигатели легкого топлива |
|
|
|
|
||||
|
|
|
Сравнивая |
(1-43) |
|
с |
|||
|
|
(1.33), можно видеть пол |
|||||||
|
|
ную их аналогию, поэтому |
|||||||
|
|
приведенные |
ранее |
выводы |
|||||
|
|
полностью |
относятся |
и |
к |
||||
|
|
этой системе регулирования. |
|||||||
|
|
Заметим, что и в этом слу |
|||||||
|
|
чае |
стабилизация |
процесса |
|||||
|
|
регулирования |
происходит |
||||||
|
|
лишь за счет положительно |
|||||||
|
|
го самовыравнивания в объ |
|||||||
|
|
екте регулирования. |
|
|
|||||
|
|
|
Дальнейшее |
улучшение |
|||||
|
|
процесса |
регулирования |
рк |
|||||
|
|
можно осуществить введени |
|||||||
|
|
ем жесткой |
обратной |
связи |
|||||
|
|
между регулирующим |
орга |
||||||
|
|
ном и золотником, как это |
|||||||
Фиг. |
1.11. Принципиальная схема астати |
показано на схеме фиг. |
1. |
12. |
|||||
|
ческого регулятора иаддува. |
В |
этом |
случае |
уравнение |
||||
нее |
соответствовавшее1 астатическому |
движения сервомотора, ра |
|||||||
звену, превращается в ста |
тическое звено. Если коэффициент жесткой обратной связи взять равным р, то вместо (1.32) получим
|
|
|
(1.44) |
Разрешая (1.19), (1.42) и (1.44) |
относительно Хр , |
получим |
|
f,,+ ( r*. + Вш/С3) г |
К 2Ка f |
1 |
|
■ х°+ |
’ + l) |
F„ |
(1.45) |
Из полученного видно, что при возмущении, обусловленном как перенастройкой регулятора Х°, так и нагрузкой /> к, система стала статической, ввиду того что
X <> (оо) = Нт
р- °
Х х ‘(оо) = Нш Р-0
(
7V
—ТРК^1.2+
^— |
Р + 1уFPK |
|
СО QQ- |
4-7” \ |
к2к, |
/ 1 |
||
\№3^ Рк Г |
р +1 |
|
|
к2к. |
|
|
----- х° |
|
( 1 |
р |
к2ка |
\ |
||
W 3 +ТРк) р + |
р |
к .
К2К,Л-Р
_ к2к„
Р+ к *к .
4. Динамика системы регулирования |
29 |
Изменением значения р можно получить заданную |
высотную |
характеристику двигателя, т. е. заданный (линейный) закон изме-
нения рк от высоты полета. |
|
|
|
|
|||||
Поэтому |
жесткая обрат |
|
|
|
|
||||
ная связь в рассматрива |
|
|
|
|
|||||
емой системе регулирования |
|
|
|
|
|||||
рк обычно применяется лишь |
|
|
|
|
|||||
для получения заданной вы |
|
|
|
|
|||||
сотной характеристики |
дви |
|
|
|
|
||||
гателя, а не для улучшения |
|
|
|
|
|||||
переходного |
процесса, |
так |
|
|
|
|
|||
как за счет большого |
само- |
|
|
|
|
||||
выравнивания объекта регу |
|
|
|
|
|||||
лирования |
переходные |
про |
|
|
|
|
|||
цессы получаются |
приемле |
|
|
|
|
||||
мыми и без |
дополнительной |
|
|
|
|
||||
стабилизации. |
|
|
|
|
|
|
|||
Теперь рассмотрим нели |
|
|
|
|
|||||
нейную систему регулирова |
|
|
|
|
|||||
ния р*, учитывая, что |
нели |
|
|
|
|
||||
нейность образуется за счет |
|
|
|
|
|||||
наличия |
перекрытия |
окон |
|
|
|
|
|||
золотника для системы, ра |
Фмг. |
1. 12. Принципиальная схема |
регуля |
||||||
ботающей с жесткой обрат |
|||||||||
ной |
связью, |
согласно |
тора |
наддува с жесткой |
обратной |
связью. |
|||
фиг. |
1. 12. |
|
|
такого |
вида |
нелинейности |
приведена на |
||
Характеристика |
|||||||||
фиг. |
1. 13, |
из которой следует, |
что движение сервомотора начнется |
Фиг. 1. 13. Характеристика нелинейного элемента.
лишь в том случае, если величина управляющего сигнала |Хупр| где b — величина перекрытия окон золотника.
Для решения такой задачи воспользуемся методом фазовой плоскости. Исходными уравнениями будут
30 |
Глава 1. Поршневые двигатели легкого |
топлива |
|
|
|||||||||
|
(7ркр + \ ) Х Рк= - К ,Х , - , |
|
|
|
|
|
|
|
|
|
|||
|
|
К о Л= к 2х Рк. |
|
|
|
|
|
|
|
|
|
||
|
|
^уПр= ^"зол |
|
|
|
|
|
|
|
|
(1.46) |
||
|
|
|
|
|
|
|
|
|
|
|
|||
|
р Х а= 0, |
когда I A'yjjpl < |
b\ |
|
|
|
|
|
|
|
|
||
|
р Х , = К3 (X |
УПР |
— b) sign X |
, когда |
I |
\Х |
|
\>Ь, |
|
||||
|
|
|
} W1&xi * 1ynpl |
|
|
упр I |
|
|
|
||||
где |
sign Хупр обозначает знак перед координатой .Уупр. |
|
|||||||||||
|
Будем рассматривать движение фигуративной точки на фазо |
||||||||||||
вой плоскости, |
образуемой координатами рХр |
= у |
и ХРк. |
||||||||||
|
|
|
|
|
Из фиг. 1.13 видно, что от |
||||||||
|
|
|
|
|
ключение регулятора происходит |
||||||||
|
|
|
|
|
в тот момент, |
когда Х упр= ± Ь , |
|||||||
|
|
|
|
|
поэтому |
|
вся |
|
фазовая |
плоскость |
|||
|
|
|
|
|
должна |
разделиться |
на три об |
||||||
|
|
|
|
|
ласти. Уравнение линий, разде |
||||||||
|
|
|
|
|
ляющих области фазовой плос |
||||||||
|
|
|
|
|
кости, получим при совместном |
||||||||
|
|
|
|
|
решении |
|
(1.46) с |
учетом, что |
|
|
|
|
РХр„ =У> т- е- |
|
|
|
|
|
|
К . Ь |
к0к2 X. . |
|
|
|
|
|
У = : утр |
||
|
|
|
|
Рк 1Тр |
Рк* |
|
Фиг. |
1. 14. |
Кривые на |
фазовой |
пло |
|
0-47) |
|
|
|||||
|
|
скости. |
|
Это уравнение двух парал |
||
ным |
|
|
|
лельных прямых с отрицатель |
||
коэффициентом наклона, которые отсекают на оси абсцисс |
||||||
отрезок, |
равный ± |
° |
, как это показано |
на фиг. |
1.14. |
|
|
|
|
г ”Т“ К 2*^ ст |
|
|
Фазовые траектории внутри полосы А — В найдем, если про дифференцируем первое уравнение (1.46), принимая рХр^ = у , и
второе уравнение (1.46). В результате получим такие два урав
нения: |
_ |
у _ . d x |
|
d y |
= У - |
||
d t |
|
d t |
|
Исключая время из (1.48) и интегрируя, получим
X
у = - - ^ + С .
(1.48)
(1.49)
Это также уравнение прямых, но с меньшим отрицательным коэффициентом, чем в (1.47).
4. Динамика системы регулирования |
31 |
Значение постоянной интегрирования С зависит от начальных условий.
Направление движения вдоль траекторий определяется тем, что выше оси Хр находятся положительные значения у , как это показано на той же фиг. 1. 14.
Для отыскания фазовых траекторий правее полосы А — В, где
p X a = K 3(Xynf — b), |
решив |
это уравнение совместно с (1.46), по |
|||||||
лучим |
|
|
|
|
|
|
|
|
|
ТРк df t |
+ (ТРк К $ + |
1)у + К3(Р- |
К2Кс) ХРк = К3КаЬ. |
(1.50) |
|||||
Разделив |
(1.50) |
на р Х Рк = у , |
получим |
|
|
|
|||
|
|
dy |
а \ — |
а 2 Х р к |
а3, |
|
|
(1.51) |
|
|
|
dX Рк |
|
У |
|
|
|
||
где |
|
|
|
|
|
|
|
||
„ |
К 3К аЬ _ _ |
к 3 / 0 , |
^ |
ч |
Ч К з Р + 1 |
|
|||
а\ ---- Z |
> а2—Z. |
(P + ^V^°)' а3— |
~ |
|
|||||
|
1 Рк |
|
1 р к |
|
|
|
|
‘ Р* |
|
Интегрируя это уравнение (переменные разделяются с по |
|||||||||
мощью подстановки ХРк= г + Н и выбора |
такого h, когда аг— |
||||||||
— аф — 0, а |
далее с помощью |
подстановки |
y /z = u ), |
получим |
|||||
такие уравнения фазовых траекторий: |
- а3+ |
Яд —4а3 |
|
||||||
|
|
|
|
|
|
|
|||
|
«з — |
а\ — 4а2 |
|
|
|
1У “3 ~ 4 |
|
||
|
' |
а2/ |
|
|
X |
|
|||
у + - |
2 |
' |
|
|
|
|
|||
|
|
|
|
|
|
о3+ У а |
|
|
|
|
3 + jX a3—4я2 |
|
Ui \ |
гУ 4 - |
4а 2 |
|
|||
|
|
|
|
|
(1.52) |
||||
X у + - |
О |
' |
V |
я2/. |
|
|
|
||
для случая, когда корни уравнения |
|
|
|
|
|||||
|
ТР Х * + (TP K S + |
1) |
К3(Р - |
К2К,) = 0 |
(1 .5 3 ) |
комплексные, и
для случая, когда корни уравнения (1.53) действительные.
32 |
Глава I. Поршневые двигатели легкого топлива |
На фиг. |
1. 14 нанесены фазовые траектории, соответствующие |
(1.52), из рассмотрения которых следует, что все они стягиваются к отрезку на оси ХРк, ограниченному полосой А—В, который
и является местоположением равновесных точек системы. Отсут ствие замкнутых фазовых траекторий (предельных циклов) указы вает на отсутствие автоколебаний в системе. На фиг. 1. 15 нанесены фазовые траектории, соответству
ющие (1.53а).
Следовательно, рассматрива емая система регулирования, с точки зрения ее устойчивости, яв ляется приемлемой.
Двигатель с ВИШ и ТК
Рассмотрим систему регулиро вания, когда включены регулято ры числа оборотов вала винта и наддува. Раздельное рассмот рение контуров регулирования хо тя и проще, но получаемые ре зультаты менее точны.
Напомним, что из выявленных выше свойств рассматриваемого двигателя как объекта регулирования следует, что для регулиро вания величины наддува необходимо применять регуляторы со ста билизирующими устройствами.
С и с т е м а р е г у л и р о в а н и я , д в и ж е н и е к о т о р о й о п и с ы в а е т с я л и н е й н ы м д и ф ф е р е н ц и а л ь н ы м
ур а в н е н и е м
Вкачестве регулятора числа оборотов вала винта примем преж нюю схему, приведенную на фиг. 1.8; схема регулятора наддува показана на фиг. 1. 16, из которой ясен способ ее работы. Стаби-
Фиг. 1. 16. Принципиальная схема регулятора над дува с сигналом от п\.
4. Динамика системы регулирования |
33 |
I
лизирующим устройством является рассмотренное выше устрой ство, дающее сигнал, пропорциональный величине ускорения вра щения турбокомпрессора.
Фиг. 1. 17. Структурная схема.
Общая структурная схема всей системы регулирования двига теля представлена на фиг. 1. 17, где указаны обобщенные коорди
наты. Там же показаны возмущения Х° и Х Рк, обусловленные из
менением заданного значения числа оборотов вала винта, и соот ветственно величины наддува. В качестве уравнений движения объекта регулирования и звеньев регуляторов используем ранее полученные уравнения (1.24), (1.31), (1.32), (1.26), (1.29), (1.30), (1.42), (1.39), (1.40):
(.Т'гР+ |
P i) Х п— К 1РкХРк— K ivX ? - f K lPrX Pr\ |
|
|||
|
|
X . » = к 2 (Ха- Г пУ, |
|
||
|
|
рХу = Х 3Х.П1, |
|
||
{Т2р + р2) ХпХ= КРгХРг — К2РкХр^ -J- KaW\ |
|
||||
K3niXnl + КЗРкХРк—К3пХп= 0; |
(1.54) |
||||
^4ркХРк+ |
XtnX„ — K.i°Xz— K4PrXPr —0 ; |
||||
|
|||||
* . з о л = К 3(ГРк- Х Рк); |
|
||||
р Х а = К3Х\ |
|
||||
у __ у |
___у |
|
|||
'*■ |
'Мзол |
ул буксы» |
|
||
•^бу к сы |
^ 4 Р^ n1' |
|
3 207
34 |
Глава I. Поршневые двигатели легкого топлива |
Рассмотрим поведение системы, когда возмущение наносится мгновенным заданием нового значения числа оборотов Х°п и но
вого значения наддува Хрк. |
|
|
Решив (1. 54) относительно координат Хп и Х.Рк, |
получим |
|
(а0р* + а У + а2р2+ а3р - f а4) ХП= |
Л |
|
= (b0P2+ b lP + а4) х : + {Ъ2р + Ьг) рХ0?', |
(1.55) |
|
( а у + а У + а У + а3р + а4) х Рк = |
||
|
||
= ( Ь У + Ь5р + а4)Х®к + (Ь6р-\-Ь7)р'Х°„. |
] |
Здесь коэффициенты левых и правых частей уравнений образованы различной комбинацией коэффициентов системы (1.54).
Из (1.55) видно, что относительно регулируемых параметров система остается астатической, так как
*„(оо ) = Х0; Л ^ ( о о ) = Х ° , к .
Рассмотрим, как изменяется запас устойчивости системы при /изменении зна
чений коэффициентов усиления К4, Кз/Кз и коэффициента самовыравнивания р 2- Свободные колебания системы определяются левыми частями уравнений, и они определяют устойчивость системы /регулирования. Если, все коэффициенты характеристического уравнения положительны, то условия устойчивости можно
записать в таком виде:
«з (ага2 — а0а3) — а\а4 > 0 . |
( 1 .56) |
Так как в эти коэффициенты входят все коэффициенты исходной системы уравнений, при этом в очень сложной зависимости, для исследования влияния интересующих нас коэффициентов на устойчивость будем рассматривать раздель но влияние каждого коэффициента, принимая остальные за постоянные.
Если таким образом преобразовать (1.56) относительно коэффициента Рг, определяющего величину самовыравнивания турбокомпрессора, то полечим нера венство вида
( 6 ip| 4" б'грг + 6 3 ) Р2 -f- > (С5Р2 + С6) р2+ б 7,
из которого следует, что с увеличением р 2 запас устойчивости системы увеличи вается. С увеличением высоты полета значение р 2 уменьшается, поэтому на боль ших высотах полета запас устойчивости системы релулирования уменьшается.
Если аналогичные действия проделать относительно коэффициента Kt, опре деляющего величину сигнала по ускорению вращения турбокомпрессора, то по лучим неравенство вида
6V64 -f- С7К\-\- СЛК\ -j- 69/64 —|—C’io> C i 1AT4 -f- С12К 4 + С 13,
из которого также следует, что с увеличением этого коэффициента запас устой чивости увеличивается.
Если положить, что Kt —■0, т. е. отсутствует сигнал, пропорциональный уско рению вращения турбокомпрессора, то для обеспечения устойчивости требуется, чтобы Cjo> C is, что практически_на больших высотах полета не соблюдается.
Введем обозначение 1 =Кз/Кз, что соответствует отношению коэффициентов усиления в уравнениях движения сервомоторов ВИШ и регулятора наддува. За
менив Ка=ЧКз и считая в дальнейшем все коэффициенты системы уравнений за данными, за исключением у , получим условия устойчивости в виде
«з (аха2 — айа3)> a\Ci.
4. Динамика системы регулирования |
35 |
Из этого неравенства следует, что чем меньше значение Y, тем запас устойчиво сти больше. В действительности f <1, поэтому уменьшение его значения приводит к тому, что оба регулятора начинают работать на существенно различных часто тах или с существенно различными скоростями, и тем в меньшей степени один регулятор влияет на другой. В пределе, когда '! -* 0, запас устойчивости будет наибольшим. Это соответствует случаю, когда один из регуляторов не работает.
Для еще большей стабилизации системы регулирования целесообразно в рассматриваемую схему регулятора наддува ввести жесткую обратную связь. В этом случае регулятор величины наддува будет работать с остаточной неравно мерностью, что иногда практически и делается.
С и с т е м а р е г у л и р о в а н и я , д в и ж е н и е к о т о р о й о п и с ы в а е т с я н е л и н е й н ы м д и ф ф е р е н ц и а л ь н ы м у р а в н е н и е м
Схемы регуляторов числа оборотов и величины наддува оста вим прежними, но учтем нелинейность, образующуюся за счет зоны нечувствительности и ограничения скорости движения серво мотора ВИШ. Форма нелинейной зави симости приведена на фиг. 1. 18.
Определим условия возникновения автоколебаний в системе регулирования, для чего воспользуемся методом гармо нической линеаризации уравнения движе ния нелинейного звена.
Для проведения такого анализа предполо жим, что колебания сигнала на входе в нелиней ный элемент в замкнутой системе будут почти синусоидальными и что высшие гармоники сигна ла на выходе из нелинейного звена не учитыва ются. Практически это означает, что сигналы, вхо
дящие из нелинейного звена и соответствующие высшим гармоникам, проходя через линейную часть системы, настолько уменьшаются по амплитуде, что ими можно пренебречь.
Уравнения движения, входящие в (1.54), останутся прежними, за исключением третьего уравнения, которое теперь запишется в виде
pX9 = F (X J . |
(1.57) |
Так как отыскиваются условия возникновения автоколебаний, то принимаем Х° = Хрк = 0.
Когда система находится в режиме автоколебаний, все ее об общенные координаты находятся в этом же режиме; поэтому не имеет значения, относительно какой обобщенной координаты оты скивать автоколебания, если не учитывать величины амплитуды. В данном случае используем координату Хзал, входящую в линей ную и нелинейную части уравнения движения.
Решая |
систему уравнений |
движения относительно координа |
ты Х30Л, получим |
|
|
(аоД3 + аУ |
+ а2р + а 3)рХ30Л+ |
(а4р2 + а5р + а6) F (XaJ = 0. (1.58) |
3* |
|
|
36Глава /. Поршневые двигатели легкого топлива
Вэтом уравнении F(Xзал) является нелинейной функцией. В со ответствии со сказанным выше о том, что колебания в системе будут почти синусоидальными, принимаем
Хзм = A sin at,
(1.59)
F (X 30!I) = F (A sin at).
Разлагая второе выражение (1. 59) в ряд Фурье и ограничиваясь
первой гармоникой, получим
2с
F(A sin о>t)= F(A sin и) s z — f F(Asinn) du-f- 2n J
о
4K
— j F (A sin u) sin и du sin и +
\_
2*
+ 71 ^ ( A sin u)cos и du cos u, 0
где u = mt.
Здесь первый интеграл представляет некоторую постоянную ве-
личину. Для упрощения примем J F(A sin и) du — О, что принципи-
о
ально не нарушает дальнейших общих рассмотрений. Из вида нелинейности, приведенной на фиг. 1. 18, видно, что функция являет ся отрицательной, поэтому члены ряда Фурье, содержащие cos, бу дут равны нулю. Окончательно получим такое выражение:
F (A sin at) * |
F (Л) sin at = ^ A л ;зл, |
(1.60) |
||
где |
2 г. |
|
||
|
|
|||
Д(А) = — |
^ Д (A sin ш£) sin atdt; |
(1.61) |
||
или |
о |
|
||
2 г. |
|
|
||
F(A) |
|
|
||
|
sin at) sin atdt. |
(1.61а) |
||
А |
т.А |
|||
|
|
Уравнение (1.61) является выражением для коэффициента ряда Фурье. Коэффициенты, входящие в выражение ряда Фурье, при нимаются постоянными потому, что постоянными приняты ампли туда А и частота со автоколебаний. Для переходного же процесса, когда A =var и co=var, значения этих коэффициентов также изме няются.
4. Динамика системы регулирования |
37 |
Таким образом, второе выражение (1.59) можно записать так:
F (Хзм) =* F(A sin w t ) ^ F ( A ) sin wt. |
(1. 62) |
Подставив первое выражение (1.59) и (1.62) в (1.58), полу чим
(а0р3 + а^р1+ а2р + а3) рА sin ast +
+ (а4р2 + а5р + а6) F (A) sin cat = 0. |
(1.63) |
Продифференцировав соответствующее число раз A sin а>(, как этого требует (1.63), и разделив его на два уравнения, содержа щие только sin и только cos, получим
а0ш4 —а2си2 + ( а 6 —а4ш2) ^ ^ = 0 ; |
(1.64) |
а3ш —а 1ш3 + а6ю ^ -^ = 0. |
(1-65) |
Из этих двух уравнений можно определить амплитуду А и часто ту (о автоколебаний, т. е.
|
F{A) |
Дз |
_____ а\ |
^3^4 — |
|
||
|
|
|
|
2a5(a0ab— a1ai) [^5^2 """ ^1^6 |
|
||
|
: Y |
al al + 4 К а 5 - |
+ а3аА— а2а5)] ==р,; |
(1.66) |
|||
|
|
О |
|
1 |
\a2as— axaQ— a3aA± |
|
|
|
|
|
2 (ДоЛ5 — |
|
|||
|
|
|
|
|
|
||
|
Y |
al al + |
4 К «5 - |
а 1а4)(а1а6+ а3а4- а2а5)] = р2. |
( 1. 67) |
||
Из |
(1. 64) |
и (1. 65) |
видно, что при определенных условиях в си |
||||
стеме |
автоматического |
регулирования могут быть |
автоколебания |
с четырьмя значениями частот и амплитуд.
Реально существующим автоколебаниям должно соответство вать действительное значение Л и со, определяемые (1. 66) и (1. 67). Отсюда следует, что для исключения автоколебаний параметры регулятора, входящие в коэффициенты правых частей (1.66) и (1.67), должны так выбираться, чтобы не получалось действи тельного значения Л и со.
Пользуясь (1.66) и (1.67), |
эти условия |
можно записать так: |
Р ,< 0 ; |
р2< 0 . |
(1.68) |
Совершенно ясно, что можно применять какое-либо одно усло вие (1.68).
Граничным значением условий возникновения автоколебаний
будет |
(1.69) |
Pi “ О; Р2= 0. |