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

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

.pdf
Скачиваний:
8
Добавлен:
29.10.2023
Размер:
11.67 Mб
Скачать

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

В этом случае получим такие значения параметров:

 

 

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

 

Варианты

К 1

Ко

та

 

5

15,7

0 , 2 1 2

1,24

6

5,5

2,3

0,635

Значение постоянной времени сервомотора

Та =

—“

= 4,7 сек.; Тс2 — ~— = 0,43 сек.

1

0,212

2

2,3

Полученное значение постоянной времени 7с2=0,43 сек. близко к практически применяемому.

Коэффициенты уравнения в соответствии с (3. 13) и (3. 14) для 6 -го варианта будут такими:

« 0

=

0,318;

=

1,87;

« 2

«1

=

2,64;

ао =

4,15;

«2

=5,11; « 3 = 4,15;

=0,21; « 3 = 0,81.

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

ности. Нецелесообразно также при данном значении- Л/В= 1,5 выбирать изобра­ жающие точки с различными значениями А и В. так как во всех остальных слу­ чаях интегральная погрешность также будет увеличиваться. Полученные же зна­ чения коэффициента усиления A'i и время изодрома Т„ являются вполне прием­ лемыми.

Определим величины интегральных погрешностей для 4, 5 и 6 -го вариантов, воспользовавшись (3. 47) и полученными выше значениями для Т Ко и Ки

Соответствующие значения будут такими:

 

 

Варианты 4

5

6

 

 

 

 

 

 

 

Хо

0,25

1,10

0,55

 

 

 

Для проверки правильности

выбранных

параметров

построим переходные

процессы для 5 и 6 -го вариантов.

 

 

 

 

 

 

Для 5-го варианта имеем следующие коэффициенты уравнения:

 

а0 =

0,62;

« 1

= 1,87;

« 2 =

2,62;

«3’= 1 ,1 ;

 

«, =

1,35; «о=1.Ю ;

а2 =

0,41;

« 3 =

0,415.

 

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

 

(0,62/>з +

1,87/>2+2,62/7+1,1)*„ = (1,35р +

1,10) ХО.

 

Корни характеристического уравнения будут

 

 

 

X, = — 0,63;

Х2’= — 1,2 +

г-1,13;

Х3 =

— 1,2 — г-1,13.

 

Изображение искомой координаты при условии ХО =

1 [^] будет

 

1,355 +

1,10__________=

Ср

 

С,

Со

Сз

о ^,62S3 + 1,8752 +

2.62S + 1,1)

S

+ S — X, + S — Х2 +

S — Х3'

 

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

219

Огсюда находим значения постоянных Со, Сь С2, С3:

 

 

1,10

 

1,35*!+ 1,10

 

1,35*2+1,10

 

°

*1*2^3 '

1

*1 (*1 —*г) (^1 —*з) ’

2

*2 (*2 М) (^2

Ч)

 

 

 

 

1,35*3 + 1 ,1 0

 

 

 

 

 

 

3

4 i 4 — ^i) (^з— Ч)

 

 

Следовательно,

уравнение переходного процесса

будет таким:

 

 

 

 

Хп (О = Со + С / - ' + С2ех*' + С„еЧ

 

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

 

Х„ (/) =

1,0 — 0,31е“ °'63< — 1 .lie -1 '2* sin (1,13/ +0,65).

 

Соответственно для 6 -го варианта получим

 

 

 

 

X n (t) =

1,0 — О.ЗОе-1'27' — 1 ,le~2,3/sin (2,21/ + 0,65).

 

Соответствующие переходные

процессы даны яа

фиг.

3. 58.

 

Фиг.

3.58.

Переходные

процессы

с

изодромным регулятором.

Перенастройка

регулятора: /—вариант 5,

 

 

2—вариант 6.

5. 4—ва­

Изменение нагрузки: <?—вариант

 

 

риант 6.

 

Из сравнения кривых переходного процесса

видно,

что

лучшим

является

6 -й вариант, которому соответствует меньшая интегральная

погрешность.

Пример 2

З а д а н и е .

Проверить качество переходных

процессов системы, рассчитан­

ной в первом примере, но при возмущении на двигатель.

 

 

и 6 -го

И с х о д в ы е

д а н в ы е. Воспользоваться

результатами расчета 5

вариантов первого примера.

 

 

возмущении /° будет знакопо­

Р е ш е н и е .

Предполагая, что процесс при

стоянным, определим величины интегральных погрешностей согласно (3.48).

Для

5-го

варианта

получим

/у» =0,38; для

6 -го — /уо=0,20;

для

4-го —

Jjo =0,08.

 

 

 

 

 

 

 

 

 

Для проверки определим переходные процессы. Для 5-го варианта урав­

нение движения будет таким:

 

 

 

 

 

 

 

(0 ,62^з + 1 ,8 7 ^ 2 + 2,62/> + 1,1) Х п = р (0,41р +

0,41)/о.

 

 

Отсюда находим изображение искомой координаты при /о = 1 [/];

 

 

s

 

5(0,415 +

0,41)

_

С,

С2

С 3

 

"

S(0,62S3 +

1 ,8752 +

2 ,62S + 1.1)

 

5 — *i

S — *2 +

S — *3 '

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

Постоянные Съ С2, С3 определяем из таких выражений:

0 ,41Х]+ 0 ,4 1 0,4U 2 + 0,41 0.4П 3-{-0,41

Cl~ ( h ~ h ) ( h ~ h ) '' 2=

( h - h ) & 2 - h ) ''

3 _ ( h ~ h ) ( h ~ h ) '

Уравнение же переходного процесса примет вид

 

х п(i) =

c ie1' 4 c 2tl!4 ^

‘'

После подстановки цифровых значений и преобразований для пятого вари­ анта получим

Х а (/) = О .Н е-0'63' — О.ббе-1'1^ cos (1,13< + 1,26).

Аналогично для шестого варианта получим

 

Х п (0 = 0 ,29е—1'29/ — 0,34е-2,32 cos (2,2It +

0,52).

 

На фиг. 3. 58 приведены и эти переходные процессы.

 

не может

Еще раз

подчеркнем, что возмущение

в виде практически

существовать,

так как изменение /° скачком

по существу

означает

изменение

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

Пример 3

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

И с х о д н ы е д а н н ы е . Уравнение движения двигателя принять из первого примера (стр. 83), а уравнение движения всей системы регулирования — по (3. 56).

Р е ш е н и е . Воспользуемся диаграммой Вышнеградского, приведенной на фиг. 3.56. Определим значения коэффициентов К,, Кс, Ки Ti, при которых переход­ ный процесс оставался бы монотонным, а время регулирования было бы наимень­ шим. Для этого выберем параметры Вышнеградского по кривой, соединяющей минимумы границы монотонности (кривая аЕ), и воспользуемся выражениями для коэффициентов Вышнеградского:

аг

Т + Т, р,

 

1

. A

f ' _

( Г Г ^ ^ А А ) 1'3

|

а2

*lAc (Аа + К\Т\) -)-Pl

 

1

 

 

(7Ti)1/3(*iAiAc)2/3

1

 

 

Г

 

 

 

 

1

.

° i4 '3

r P ^ A jA c ) 1'3

 

1

 

 

Из полученного видно, что параметры А», Кс, Ai, Т не определяются одно­ значно. Однако поставленную задачу можно решить, если задаться значением Т\, а из первых двух выражений (3. 68) определить общий коэффициент усиления

K=biKiKc и Аа, т. е.

+ 7~1Рт)3

Аа _

ТТ\А

Г

Д2ИГ7-, '

А Ч П \

Aj

T +

Tl9l

[

(3.69)

(Т- + Г1Р1)2.

Связь с параметрами Вышнеградского

получим,

если перемножим первое

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

2 21

и третье выражения (3.68), т. е.

АС = Г + Г1;>1- .

(3.69а)

Т

 

Далее на кривой, соединяющей минимумы гравии монотонности диаграммы Вышнеградского, выбираем точку, соответствующую найденному значению произ­

ведения ЛС_П|ри заданном значении Т\, которая и определяет значения парамет­ ров А, В и С.

Фиг.

3.59. Зависимость вели­

Фиг.

3.60.

Переходные процессы.

чины

интегральной погрешно­

Т,

АГа=1,4;

ЛГ,=4,5;

Г,-0,2;

 

сти от значения 'Ti.

2: J—ЛГа=2,9;

Л",=12,9;

Г,=0,1;

 

 

 

 

3; 6 - А"а=5,8;

/<-,=41,5;

Г,=0,05.

Параметры, подсчитанные для различных Т\ при Кс=2,3, имеют следующие значения:

 

 

 

 

 

Параметры

 

 

 

Тх

АС

А

В

С

К

Кс

к а

*1

К ,

 

К 1

 

 

 

 

 

 

 

 

 

0,2

1,4

2,15

2,60

0,65

3,4

2,3

1,4

4,5

0,31

0,1

1,2

2,05

2,50

0,58

9,8

2,3

2,9

12,9

0,23

0,05

1,1

2,02

2,45

0,54

31,6

2,3

5,8

41,5

0,14

Значения

интегральной

погрешности,

вычисленные по выражению

 

 

 

 

 

Х°

l +

ftlKaftc

 

 

 

 

 

 

 

biKiKc

'

 

 

приведены на графике фиг. 3. 59.

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

Значения коэффициентов уравнений (3.56) в соответствии с полученными выше значениями параметров регулятора будут такими:

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

т1

“ 1

а 0

а2

а3

*0

« 1

CL^

«3

 

0 , 2

0 , 6 8

3,4

0,066

0,33

0 , 1

0,7

2,75

3,4

0 , 1

0,98

9,8

0,033

0,33

0,05

0 , 6

4,18

9,8

0,05

1,58

31,6

0,016

0,33

0,025

0,55

6,96

31,6

Переходные же процессы, построенные аналогично тому, как это уже дела­ лось в предыдущих примерах, приведены на фиг. 3. 60. Уравнения движения соот­ ветственно для Ti = 0,2, /i = 0,l и /i=0,05 будут:

Х п it) = 1,0 - 0 ,9 4 e - 1'8W + 0,54e-2'56',sin (3,36/ + 0,09);

Х п (/) = 1,0 — 0 ,93е~3-69' + 2,6 8 ? ~4’16' sin (5,97/ + 0,025);

Хп (О = 1.0 —2,41е-7'4< + 2 ,13а-7,3< sin (10,84/ + 0,72).

Из приведенных переходных процессов видно, что они монотонные.

Пример 4

З а д а н и е . Проверить качество переходных процессов

системы,

рассчитан­

ной в 3-м примере, но при возмущении /° на двигатель.

 

расчета

предыду­

И с х о д н ы е

д а н н ы е .

Воспользоваться

результатами

щего примера.

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

что процесс при возмущении

/° остается

знако­

Р е ш е н и е .

постоянным, определим величины интегральных погрешностей по (3.63).

 

 

 

7-1 =

0,2;

0,1; 0,05;

 

 

 

 

 

 

//„ =

0,098;

0,034;

0,011.

 

 

 

 

Кривые переходных процессов, построенные аналогично тому, как это уже

делалось в предыдущих примерах, для

случая

/ i = 0 ,2 ; / 1= 0, 1

и /1=0,05

приве­

дены на фиг. 3. 60, уравнения которых соответственно выразятся:

 

 

Х п (/) =

0, 17й_ 1,88/ — 0,24е—2’56f‘ cos (3,36/ +

0,83);

 

 

Х п (/) =

0 ,12g—3,697— 0 ,15й—4,167 cos (5,97/ +

0,64);

 

 

Х п (/) =

0, Обе- 7 , 47 — 0,098а-7,37cos (10,84/ + 0,87).

 

 

Из кривых видно, что они остаются знакопостоянными, причем лучшему про­

цесс) соответствуют параметры

/(«=5,8; Ki=41,5; 7i=0,05.

 

 

 

 

Пример 5

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

ходных процессов,

обусловленных перенастройкой регулятора с

акселерометром

и изменением нагрузки при идеальном акселерометре (7i = 0 ).

 

И с х о д н ы е

д а н н ы е . Уравнение движения двигателя принять по преж­

нему примеру, а всей системы регулирования — в соответствии с

(3. 60).

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

223

Р е ш е н и е . Условием монотонности для рассматриваемой системы пли на­ чальных условиях Х„(0) = —1; ХП (0)=0 будет

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

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

\2

АЬ\К\КС

(3.70)

 

— Ч г—^ > 0 .

Примем по-прежнему, что /(== 2,3;

тогда определению подлежат

парамет­

ры К, и А».

 

 

Фиг. 3.61. Зависимость ве­ личины интегральной по­ грешности от значения Аь

Фиг. 3. 62. Переходные процессы.

I ; 3 - K a- 2,1; А',=4,1;

Z 4 - K - Z A Л'.-в.гб.

Будем выбирать величины этих параметров из условия минимума интеграль­ ной погрешности, определяемой по (3.64), т. е.

 

,

1 +

b\KcKa

 

х"

i-iATj/Cc

На границе монотонности получим

 

 

/ ^ i * c * a -f- ^ \ 2

4 Ь \ К \ К с

\

Т

) ~

Т

откуда, считая по-прежнему заданными bt

и Т,

получим

 

Аа = — 1,33 ±

1,63/ А л

Задаваясь значениями Ai и определяя значения А», найдем величину инте­

гральной погрешности. На

фиг. 3.61

приведены

кривые изменения значений Jjr„

и Л’« для различных значений Ai, из которых

виден монотонный характер умень­

шения

о с увеличением Ал Практически

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

Ai>30-:-40, так как процесс от этого существенно не улучшится, а коэффициенты усиления тахометра Ai и акселерометра А» получаются очень большими, практи­ чески трудно выполнимыми. Переходные процессы для случаев А«=2,1 и /(„=3,4 приведены на фиг. 3.62, уравнения которых соответственно напишутся

Х п (/) =

1,0 + 0 ,17е-4'54'

-

1,17е °.66/;

* „ ( 0 1 =

1 ,0 -)- 0,19е—6,18/

1 ,19е—0,S8<-

224

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

 

 

 

Выбирая большие значения Кп, получим лучшие переходные процессы. При

возмущении /° на

объект регулирования величины

интегральных

погрешностей

по

(3. 6 6 ) при тех

же значениях Ki определяются

кривой, приведенной на

тон

же

фиг. 3. 61, а переходные процессы — на фиг. 3. 62, уравнения

которых

со­

ответственно будут

Х п (t) = 0,26 (/~ 0'Ш — /-4.54/);

Х„ ( 0 = 0 . 1 9 ( / - ° ’98/ — / -6 .1W ).

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

Из рассмотренных примеров видно, что посредством регулято­ ров чисел оборотов исследованных схем можно получить заданные технические условия в части динамики процесса для максималь­ ных режимов работы одновального Т Р Д с нерегулируемым соплом.

Пример 6

З а д а н и е . Определить в общем виде условия возникновения автоколеба­ ний в системе регулирования числа оборотов турбореактивного двигателя с не­ регулируемым реактивным соплом и изодромным регулятором, приведенным на фиг. 3 .1 1 .

 

 

Фиг. 3. 63. Характер нелинейности.

 

И с х о д н ы е

д а н н ы е . Уравнения движения двигателя принять по

(2 . 34)

и (2.41),

а тахометра по (3.7).

золот­

Р е ш

е н и е .

Рассматриваемая система является нелинейной, так как

ник 4 поршня обратной связи (изодрома) 9 при самых малых его перемещениях открывает значительно большие проходные сечения, нежели сечение жиклера 2. Следовательно, истечение рабочей жидкости в полость изодрома и обратно прак­ тически происходит с постоянной скоростью, независимо от величины открытия окон золотника обратной связи. Предполагая, что окна золотника изодрома имеют некоторое перекрытие, характеристика нелинейного элемента будет такой, как показано на фиг. 3. 63, где Ъ— величина перекрытия окон золотника, а Кг— ско­ рость передвижения изодрома.

Составим уравнения движения регулятора, для чего воспользуемся схемой, приведенной на фиг. 3. 64, и структурной схемой на фиг. 3. 65. Для линейной части

регулятора

входом является координата рХ3— скорость изодрома, а выходом —

координата

— положение буксы обратной связи.

примем в виде

Уравнение движения сервомотора как астатического звена

 

Г-К4 = К С(Х 1— Х^),

(3.71

где К<— коэффициент усиления сервомотора.

Уравнение движения изодрома получим в предположении, что скорость дви­ жения поршенька изодрома (координата Х2) отличается от скорости движения

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

22”)

поршня сервомотора (координата X,) на постоянную величину скорости изодро­ ма рЛ'з, т. е.

Рх 4 р Х 3 = р Х 2.

h

 

Замечая, что Х 2 = — Х 2 — К Х 2, получим окончательно

 

‘1

(3,72)

р Х ^ — р Х з = КрХ2.

Таким образом, уравнения движения линейной и нелинейной частей системы будут такими:

(7> +

1)Хл = - * Л ;

 

/ А =

Кс (Xj — Х2);

(3.73)

/>*4 - , Р * з = * №

Рх з = ^ №)■

Передаточная функция линейной части системы, полученная из первых че­ тырех уравнений системы (3.73), выразится так:

* 2 _________Tpl + p + K.btKc___________

(3.74)

р Х з р[КТр* + (К + КсТ ) р + К с (/Ci*iAT + 1)]

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

Частотная характеристика линейной части системы W(ia>) согласно (3.74) будет такой, как это показано сплошной кривой на фиг. 3. 6 6 . Что касается ха­ рактеристики нелинейной части, то для заданного вида нелинейности определим ее аналогично нелинейной системе на стр. 33, т. е. предполагая, что автоколе­ бания носят синусоидальный характер. Производя гармоническую линеаризацию pX3—F(A sin at) разложением в ряд Фурье и ограничиваясь первой гармоникой,

15 207

226

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

Фиг. 3. 65. Структурная схе­

Фиг. 3. 66. Частотные харак­

ма к фиг. 3. 64.

теристики.

получим приближенную передаточную функцию нелинейной части вида WH=F(A)/A. По фиг. 3.63 определим интеграл (1.61а), т. е.

1

2г.

2

М|

и2—К—К,

ТГ

.(

 

 

 

т.А

 

 

 

 

0

 

 

 

 

Первый и третий интегралы равны нулю, так

как в этом случае F(Xt)= 0.

Кроме того, видно, что sin ui = Ь/А.

Следовательно, вместо (1.61а) можно на­

писать

 

 

 

—I

 

 

 

Г (А)

2

 

WH=

Кг sin и du,

А

пА

 

 

 

 

где и — iet.

 

 

 

 

 

После интегрирования

получим

 

 

 

 

 

F(A)

4АГ2 1 У j

h -

Wн — А

TlA

У

А* '

Это и есть приближенная передаточная функция нелинейной части, которая за­ висит лишь от амплитуды А. Следовательно, частотная характеристика нелиней­ ного элемента будет располагаться на действительной оси, как это показано на той же фиг. 3.66, не доходя до начала координат на величину тсЬ/2Кг, где Ь— половина зоны перекрытия окон золотника, а Кг— постоянная скорость изодрома.

Как видно из фиг. 3. 66, в случае неучета запаздывания сгорания топлива в двигателе автоколебания в системе регулирования не возникнут, так как ча­ стотные характеристики линейной ш нелинейной частей системы не пересе­ каются.

В случае отсутствия зоны перекрытия окон золотника, когда нелинейная за-

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

управления

227

влсимость имеет характер, приведенный на фиг. 3. 67,

Р*з

 

частотная характеристика такого элемента будет рас­

 

 

полагаться также ва

действительной оси, но дойдет

 

 

до начала координат,

так как при Ь 0, величина от­

 

 

резка жЬ/2Кг=0.

 

 

 

 

Теоретически в этом случае также не могут воз­

 

 

никнуть автоколебания, так как частотные характе­

 

 

ристики не пересекаются, но учитывая возможные

 

 

ошибки приближенного метода решения, гарантиро­

 

 

вать отсутствия автоколебаний нельзя.

 

 

 

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

 

 

ре сгорания линейная часть системы регулирования

 

 

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

 

 

 

(7> + 1)«^л-л = -&,*<;

 

Фиг. 3. 67. Нелинейная

X, = к х п-

(3.75)

p X ^ K c ^ - X J ;

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

золотни­

 

ка.

 

Рх 4 — Рх ъ = КрХ2.

Передаточная же функция линейной части системы будет такой:

х_2_ =

________ ТрЪ + р +

_

р х з

р [KTpi-\-(K + K J ) р + КСК ib.Ke

(3.76)

1р - К е]

Частотная характеристика линейной части системы будет другой, чем для случая т=0, так как для каждого значения о) вектор повертывается на угол то. В сбшем случае частотная характеристика линейной части будет такой, как это показано пунктиром на фиг. 3. 66. В этом случае частотная характеристика может пересекаться с характеристикой нелинейного элемента и тогда в системе будут автоколебания.

Следовательно, условием отсутствия автоколебаний будет

nb

 

(3. 77)

с <

 

Ы<2 '

 

 

а граничным условием будет

 

 

 

 

(3.78)

где с — отрезок на действительной оси, отсекаемый

частотной

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

линейной части.

 

 

Величину отрезка с можно определить из (3.76), где необходимо подставить

р = (ю, частотную характеристику выразить так:

 

 

W (/ш) = Q (ш) -f- iR (и).

 

 

Полагая R(со)=0, определим и и подставим его значение в

Q(со).

Очевидно, что граничному условию (3.78) будет

соответствовать граничное

значение величины запаздывания т.

 

 

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

Системы регулирования для одновального ТРД с регулируемым реактивным соплом

Регулируемыми параметрами для такого двигателя являются число оборотов и температура газов перед или за турбиной. Струк­ турная схема подключения регуляторов показана на фиг. 3. 1, б.

15*

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