книги из ГПНТБ / Основы автоматизации для металлургов
..pdfТак как отображение суммы функций равно сумме отображе ний, то
L |
(52<р" + 5 І <р' - И 0 < Р ] = L [р.]. |
В соответствии с табл. 18 находим: |
|
L |
[s2cp"J = s2p2<? (p) -s2p2<? (0) -s2p<?' (0); |
L |
[ s l T ' l =siP<? (P) — ^і/^Т (0); |
Z. [s0cp] =S0cp (/>); |
|
|
|
|
|
||||
Z [ц] = |
|i(/>). |
|
|
|
|
|
|||
Подстановкой этих выражений в первоначальное дифференци |
|||||||||
альное уравнение получим алгебраическое |
уравнение |
||||||||
S2P2<? |
(P) - s 2 |
P \ (0) —SiPy' |
( 0 ) ( p ) - s l p f |
(0)+s0 cp (/?) = (x (p). |
|||||
Преобразованием уравнения получим |
|
|
|||||||
~ t n \ |
= |
|
|
1 |
/ „\ |
I |
*2І2 ? (0) + |
s 2 P r ( 0 ) + ^ІРУ (0) |
|
т1 ^ |
|
S 2 / 7 2 + |
slP + Sq |
Р ^ 1- |
S2p2 + |
S l P + |
S q |
||
Член |
|
|
|
|
|
|
|
|
|
F ( / , ) = - SlP2 |
+ 1Sip + |
So |
|
|
|
|
|||
является отображением передаточной функции, а член |
|||||||||
SiP2f |
(0) |
+ |
s2/>-f ' (0) + slP-f |
(0) |
|
|
|
||
|
|
SlP2 |
+ SlP + |
s0 |
|
|
|
|
|
зависит |
от начальных условий. Подстановкой заданных начальных |
||||||||
условий |
получим |
|
|
|
|
|
(?(P)=-^pT+SlF7^0 •
При помощи табл. 18 можно получить оригинал решения:
S2P1P2 |
Р\ — Р2 |
|
где pi, р2 |
— корни характеристического уравнения. |
|
|
|
Алгебра передаточных функций |
Если |
сложную |
систему регулирования разделить на отдельные |
звенья, то свойства целого можно определить на основании иссле дования свойств отдельных звеньев и их связи. Изучаемую систему можно представить в виде блок-схемы, на которой в виде отдельных блоков изображены отдельные звенья, а стрелки показывают на
правление воздействия их выходных величин. |
В прямоугольни |
||
ках |
показаны передаточные характеристики |
соответствующего |
|
звена в графическом |
виде или в виде математического выражения. |
||
На |
рис. 77 показано |
последовательное соединение двух звеньев Я |
|
с запаздыванием. |
|
|
146
Задача заключается теперь в определении результирующей пе
редаточной |
функции контура, состоящего из нескольких звеньев |
с различной |
схемой их соединения. В качестве первого случая рас |
смотрим последовательное соединение двух звеньев, показанное на
рис. 77. Выходная величина |
первого звена с передаточной |
функ |
|||||||||||
цией |
Fi(p) |
|
является |
входной |
величиной |
второго |
звена с / ^ ( р ) , |
||||||
тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
Передаточной функцией является отношение выходного си |
|||||||||||||
гнала к входному, т. е. |
|
|
|
|
|
|
|
||||||
Л (Р) |
912 |
(Р) |
|
912 |
(р) . |
|
|
|
|
(112) |
|||
911 |
|
(р) |
|
91 (Р) |
' |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|||||
|
|
922 (р) |
_ |
92 (Р) |
|
|
|
|
|
(113) |
|||
|
|
912 |
(Р) |
|
921 (Р) |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||||
а |
|
|
|
|
|
|
_ |
|
|
|
|
|
|
|
|
- С |
|
|
|
|
|
|
|
^//1 |
|
|
|
|
|
|
|
|
|
(С |
|
|
|
П |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
û |
|
|
|
|
|
|
|
|
|
|
i L |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
<Р,=Ѣ |
. |
А |
/ |
|
|
|
|
|
|
ь, |
|
|
|
|
|
/*Т,р |
|
|
|
|
|
|
|
|
|
||
Рис. |
77. |
Блок - схема |
последовательного |
соедине |
Рис. 78. |
Параллельное |
сое |
||||||
|
|
|
|
|
ния |
звеньев: |
|
|
|
динение звеньев |
|
||
а — графическое |
и з о б р а ж е н и е |
переходных |
харак |
|
|
|
|||||||
теристик |
звеньев; |
б — аналитическое |
выражение |
|
|
|
|||||||
|
переходных |
характеристик |
звеньев |
|
|
|
Общая передаточная функция системы имеет вид
Для последовательно включенных звеньев результирующая пе редаточная функция определяется произведением передаточных функций отдельных звеньев. Для графического построения произ ведения двух передаточных функций в комплексной плоскости ко ординат используем выражение
FnU*) = \4n\-\9v\ei{*l+**\ |
(115) |
т. е. произведение двух передаточных функций определяется про изведением амплитуд и суммой фазовых углов.
Параллельное включение показано на рис. 78. Аналогично пре дыдущему случаю
Fi+2{p)=FAP)+F2{p). |
: 116) |
Результирующая передаточная функция параллельного соеди нения равняется сумме отдельных передаточных функций. В комп лексной плоскости координат результирующая передаточная
10* |
147 |
функция строится как сумма векторов передаточных функций от
дельных |
звеньев. |
|
|
|
|
Для встречного включения (рис. 79) |
|
||||
ЛР) |
F\ |
(Р) |
|
(117) |
|
1 + Fi |
ІР) F2 |
ІР) " |
|||
|
|
||||
Полупараллельная |
схема соединения показана на |
рис. 80. |
|||
В контур |
поступают две входные величины, при этом |
выходные |
величины суммируются. Результирующая передаточная функция
имеет |
вид |
|
|
|
|
|
|
|
F(P) |
= Fi(P) |
+ kF2(p), |
|
|
|
|
(И8) |
|
где |
|
|
|
|
|
|
|
|
k |
= Т21 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
S<2 |
|
|
|
|
|
|
|
С |
|
|
|
|
|
<Ргі |
hip) |
<Ргг |
/*1 |
|
|
|
|
|
|
F R |
||||
|
|
|
|
|
|
|
||
Рас. |
79. |
Встречное |
соедине |
Рис. 80. Полупараллельное |
Рис. 81. |
Разомкнутый контур |
||
|
|
ние звеньев |
|
соединение |
звеньев |
|
регулирования |
|
|
|
|
|
Система регулирования |
|
|
||
|
Д ля упрощения исследования разомкнем |
контур |
регулирования |
|||||
в |
одном определенном |
месте. Как |
показано |
на рис. |
81, контур ре |
гулирования превращается в последовательное соединение двух звеньев •— системы и регулятора. Так как воздействие регулятора должно быть направлено против направления отклонения регули руемой величины, передаточную функцию разомкнутого контура
можно выразить следующим |
образом: |
|
|
||||
F = |
-FSFR. |
|
|
|
|
|
(119) |
Если в эту передаточную функцию внести выражения для об |
|||||||
щей |
передаточной |
функции |
системы |
и регулятора, то |
получится: |
||
F(p) |
= (1 + zlP + |
+ г0 |
+ |
г 1 р |
|
+ г0 + rip |
|
. . . ) |
( S 0 |
4- S l P |
+ . . . ) |
S 0 4- (TjSo + Sl)p+ |
... • (12°) |
Дифференциальное уравнение разомкнутого контура регули рования при символической записи получится тогда в следующей форме:
2 2 х л - * 5 * { р 2 я ) = = — |
m = |
2 Гт<?І |
(121) |
п—0k=0 |
— 1 |
|
|
|
|
148
В замкнутом контуре выполняются соотношения:
?1=С Р2 = |
СР И Л |
И |
гЧ = |
г12 = |
г1, |
|
(122) |
||
откуда |
вытекает, |
что F(p) |
= — 1 . Это означает выполнение |
условия |
|||||
-ЪгтРт |
|
|
1 |
|
|
|
(123) |
||
2 IIP1 |
S |
SkP" |
|
|
|
||||
|
|
|
|
|
|
||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
+ г0 + |
Гір |
|
1. |
(124) |
||
(1 + ЧР |
+ . . . ) |
( s 0 |
+ s l |
P + |
...) |
||||
|
|
Рассмотрим теперь явления, происходящие в замкнутом кон туре регулирования. С точки зрения практического применения регулятора особый интерес представляют передаточные функции возмущения регулирующего воздействия.
б г |
•к |
9 |
ß |
|
|
ß |
|
|
|
|
|
|
|
|
|
|
<р |
|
Fr |
|
|
А, |
} |
Рис. |
82. П еред аточн ая функция |
возму |
Рис. |
83. Исходна я блок-схема |
|
|
щения |
|
|
|
|
Передаточная функция |
возмущения |
имеет |
важное значение |
с точки зрения определения изменений регулируемой величины под влиянием возмущений. При определении передаточной функции исходят из упрощенной схемы, представленной на рис. 82. Затем, исподьзуя выражение для передаточной функции полупараллель ной схемы соединения, получим
Fc |
|
|
|
(125) |
|
Р. (Р) = - |
|
|
|
||
тогда |
|
1 + |
ьр |
+ |
|
FAP)- |
|
||||
Sip + |
. . .) (1 + |
ххр |
+ ...) + пр + г0 + |
||
(S0 + |
|||||
|
2 |
^ір1 |
|
(126) |
|
Ii skpk |
2 т-іР1 + S |
rmp" |
|||
|
Аналогично определяем передаточную функцию регулирующего воздействия. Исходная блок-схема показана на рис. 83. Подставив выражения передаточных функций регулируемой системы и регу лятора в следующее равенство:
F A |
P ) = |
,1*р% |
F |
|
1 |
= |
— |
( |
т |
|
|
F |
|
|
|
|
|
||
|
|
' |
S |
R |
F0 |
|
|
|
|
где |
F 0 = |
—FrFs |
|
|
|
|
|
||
|
|
|
|
|
|
149
получим
На основании передаточных функций и регулирующего воздей ствия можно определить передаточную функцию регулируемой си стемы и регулятора. Очевидно, что:
Fs |
= |
- |
^ Y |
, |
(129) |
F R |
= |
- |
^ . |
|
(130) |
а
Величину коэффициента усиления можно установить или не посредственно из дифференциального уравнения, рассмотрев уста новившееся значение (все производные примем равными нулю при t^oo), или из передаточной функции, приняв р = 0. Для статиче ской системы (so^O)
т. е. коэффициент усиления равен обратному значению коэффици
ента So- |
|
|
|
Аналогично для |
регулятора |
|
|
l i m J i i £ ) = 1 |
і т і Ц 4 |
= Г о = /г |
(132) |
т. е. коэффициент усиления равен константе регулятора. Из этого для передаточной функции разомкнутого контура регулирования вытекает
F0=Fs,FR=-&. |
(133) |
Устойчивость системы регулирования
Устойчивость системы регулирования является одним из основ ных требований, предъявляемых к системе регулирования. Устой чивостью системы называется явление перехода в новое установив шееся состояние после появления отклонений (под влиянием воз мущения или управляющего воздействия).
Устойчивость системы можно исследовать, проанализировав дифференциальные уравнения, описывающие систему.
Решение однородного уравнения
«»«P^ + |
e » - ! ^ " " ' ^ • • • +а,т >'-Мо<Р = 0, |
(134) |
где ап, |
. . -, а0 — реальные числа, выражающие физические |
свой |
|
ства контура регулирования, имеет вид |
|
п |
|
|
150
где pu — корни характеристического уравнения, a Ck — постоянная интегрирования.
Из выражения (135) видно, что если контур регулирования должен быть устойчивым, то отклонение ср должно уменьшаться, т. е. слагаемые выражения (135) должны с увеличением времени снижаться до нуля. Это возможно лишь в том случае, если корни характеристического уравнения будут представлять собой отрица тельные числа или действительная часть комплексных корней бу дет отрицательной (рис. 84).
|
Отсюда вытекает и условие для значений коэффициентов ag,.. ., |
|||
..., |
ап, которые должны |
быть положительными, |
чтобы было выпол |
|
нено условие существования отрица |
|
л |
||
тельных действительных |
частей кор |
/ |
||
ней |
уравнения. |
|
|
|
Выражения для корней уравне ния могут определять и другие ха рактерные свойства переходных про цессов. Если комплексные корни получатся с отрицательной действи тельной частью, то процесс регули рования протекает в форме зату хающих колебаний. Положительные действительные части означают не устойчивость и постоянный рост ам плитуды, чисто мнимые корни озна чают гармоническое колебание с по стоянной амплитудой.
• |
Рг |
|
|
Р, |
YYs |
+ |
|
- |
|||
|
-J
Рис. 84. Передаточн ая функция регу л и р у ю щ е й величины:
При анализе устойчивости конту |
/ — устойчивая область; / / — неустой- |
чивая область |
|
ра регулирования нет необходи |
|
мости всегда определять действительные значения корней харак теристического уравнения. Чтобы решить, является контур устой чивым или нет, достаточно только установить, является ли отрица тельной действительная часть комплексных корней.
Знак действительной части корней проверяют при помощи кри терия Гурвица, который формулирует условия устойчивости сле дующим образом: все определители Гурвица, составленные из коэффициентов характеристического уравнения при а п > 0 , дол жны быть больше нуля.
Всего составляется п определителей Гурвица: Д0 , А\, .. ., Д п - ь Определитель До имеет п столбцов и п строк; по главной диаго нали выписываются коэффициенты характеристического уравнения по возрастающим индексам, далее, правее главной диагонали, все строки заполняются коэффициентами по убывающим индексам, левее — по возрастающим. Оставшиеся пустые места заполняются нулями. Из До последовательным вычеркиванием крайних левых столбцов и верхних строк получаются определители Д ь Дг, ... ,
..., Дп-ь При этом
Д в _ ! = <*„_!. |
(136) |
151
Например, для п — 4 |
|
|
|
|||||
|
а0 |
О |
О О |
|
а; а 0 |
О |
|
|
|
а2 |
ах |
а0 |
О |
Ai |
а3 а 2 |
ûi |
(137) |
|
а 4 |
а 3 |
а2 |
а{ |
|
О а 4 |
а 3 |
|
|
|
|
||||||
|
О О а4 |
а3 |
|
|
|
|
||
Д , = |
а2 |
ах |
|
А, |
|
|
|
|
а 4 |
аа |
|
|
|
|
|
||
|
|
|
|
|
|
|
||
Положительность |
коэффициентов ( При Эп >0) есть необходи |
|||||||
мое условие |
устойчивости. |
|
|
|
||||
Условия |
устойчивости |
Нейквиста—Михайлова оценивают поло |
жение частотной характеристики разомкнутого контура регулиро вания по отношению к критической точке (—1, /0). Контур регу-
n
Pue. 85. |
Критерий |
Нейквиста: |
Рис. |
86. |
Логарифмическая амплитудно - фа |
||
/ — устойчивая область; |
/ / — неустой |
зовая |
характеристика |
неустойчивого |
кон |
||
|
|
тура |
|
||||
|
чивая область |
|
|
|
|
|
|
лирования |
будет устойчивым, если точка |
(—1, /0) |
лежит вне |
замк |
нутой области, ограниченной частотной характеристикой, постро енной для диапазона частот ( X ^ C Ü ^ OO [критическая точка (—1, /0) должна лежать слева от частотной характеристики, если движение по ней осуществляется в направлении повышения ча стоты]. Пример устойчивого и неустойчивого контура приведен на рис. 85.
При анализе логарифмических амплитудно-фазовых характери стик условие устойчивости можно сформулировать следующим об разом: контур регулирования является устойчивым, если логариф мическая амплитудная характеристика пересекает ось 0 слева от частоты, при которой, согласно фазовой характеристике, происхо
дит фазовый сдвиг — я |
(рис. |
86). |
|
|
|
Для |
лучшей оценки |
вводится понятие «избыток фазы Ф„» (Phase |
|||
margin), определяемый |
как |
Ф„ = Ф + я; если |
избыток |
фазы для |
|
точки |
log Л (со) = 0 является |
положительным, |
то контур |
устойчи |
|
вый, в противном случае он является неустойчивым. |
|
152
Повышение устойчивости системы регулирования
Устойчивость процесса регулирования можно улучшить не сколькими способами, а именно: уменьшением коэффициента уси ления регулятора; использованием дифференцирующего звена; уменьшением чистого запаздывания контура регулирования; ста билизацией отдельных параметров регулируемой величины; ис пользованием следящего регулирования.
Рассмотрим кратко основные способы улучшения хода про
цесса регулирования. |
|
У л у ч ш е н и е у с т о й ч и в о с т и с и с т е м ы |
р е г у л и р о в а |
н и я у м е н ь ш е н и е м к о э ф ф и ц и е н т а у с и л е н и я р е г у л я т о р а . Для неустойчивого контура частотная характеристика
Рис. 87. Влияние усиления на |
Рис. 88. |
Схема теплообменника: |
||
устойчивость |
/ — термометр; |
2 — с о с у д |
д л я кон |
|
|
||||
|
денсата; |
3— управляющий |
вентиль; |
|
|
4 — змеевик |
теплообменника |
находится слева от критической точки. При уменьшении коэффи циента усиления уменьшается модуль вектора частотной характе
ристики, в результате |
чего можно достичь вывода критической |
точки Ръ. = — 1 за |
пределы этой частотной характеристики |
(рис. 87). Однако усиление регулятора нельзя уменьшать до про извольно малой величины, так как малое усиление приводит к за медлению процесса регулирования и большому отклонению регу
лируемой величины. |
|
|
У л у ч ш е н и е у с т о й ч и в о с т и |
с и с т е м ы |
р е г у л и р о |
в а н и я и с п о л ь з о в а н и е м д и ф ф е р е н ц и р у ю щ е г о з в е - н а. При использовании дифференцирующего звена фазовый сдвиг увеличивается, в результате чего улучшается значение избытка фазы. При этом фазовый сдвиг, равный я, происходит при боль
шей частоте, и, следовательно, система из |
неустойчивой может |
|||
стать устойчивой. |
|
|
|
|
У л у ч ш е н и е у с т о й ч и в о с т и п р о ц е с с а |
р е г у л и р о |
|||
в а н и я |
у м е н ь ш е н и е м |
з а п а з д ы в а н и я |
в с и с т е м е . |
|
На рис. |
88 показан теплообменник. Воздух |
нагревается паровым |
змеевиком, пар в змеевике конденсируется и отводится через ре зервуар 2. Постоянная температура воздуха поддерживается ре гулятором, который воздействует на подачу пара в змеевик.
153
Температурный датчик / расположен в потоке нагреваемого воз духа. При таком регулировании имеется значительное запаздыва ние между моментами изменения температуры воздуха и воздей
ствия на |
температуру |
воздуха |
через змеевик. При этих условиях |
|||||||||
|
|
|
|
і |
] |
температура |
будет |
сильно |
ко- |
|||
|
|
|
|
|
лебаться, |
так что |
амплитуда |
|||||
|
|
|
|
V |
1 |
колебаний будет большой, а ре |
||||||
|
|
|
|
|
|
гулирование — длительным. |
||||||
|
|
|
|
|
|
На рис. 89 показана си |
||||||
|
|
|
|
|
|
стема, в которой |
запаздывание |
|||||
|
|
|
|
|
|
полностью |
устранено. |
Подвод |
||||
|
|
|
|
|
|
воздуха разделен на две ветви, |
||||||
|
|
|
|
|
|
по одной из которых проходит |
||||||
|
|
|
|
|
|
холодный |
воздух, |
а |
по |
дру |
||
н и е . |
89. Подогреватель: |
|
гой •— воздух, |
нагретый змееви |
||||||||
/ — термометр; |
2 — с о с у д |
д л я |
конденсата; |
3 — |
ком. В этом |
случае |
регулятор |
|||||
|
р е г у л и р у ю щ и й |
клапан |
|
воздействует |
на |
|
клапан |
3, |
||||
|
|
|
|
|
|
|
||||||
при помощи которого устанавливается требуемое |
соотношение |
|||||||||||
теплого |
и |
холодного |
воздуха, |
т. |
е. требуемая |
температура |
воз |
|||||
духа. |
|
|
|
|
|
|
|
|
|
|
|
|
Оптимальный процесс регулирования
Кроме устойчивости контура регулирования, при расчете ре гулятора необходимо достичь высокого качества процесса регули рования.
Основными показателями качества являются (рис. 90) длитель ность переходного процесса (тп ), интегральная оценка отклоне
F-J |
(<p-(poa)dt |
о |
|
Рис. |
90. Характерные параметры про |
|
Рис. 91. Интегральная |
оценка |
отклонения |
||
|
цесса регулирования |
|
|
|
|
|
|
ния |
регулируемой |
величины |
от |
установившегося |
значения |
( ф ш ) |
|
и максимальное |
значение |
ее переменной составляющей |
в пе |
||||
реходном режиме |
(Афтах) - |
ЭТИ показатели должны быть как |
|||||
можно ниже. Продолжительность |
переходного процесса |
измеря |
ется от момента возникновения возмущения до момента, когда от
клонение регулируемой величины станет постоянным |
или будет |
||
ниже |
установленного |
предела (обычно выбирают ± 5 % |
регулируе |
мого |
значения ф) или |
ниже предела чувствительности |
регулятора. |
154
Т А Б Л И Ц А 18
Оптимальные константы регулятора в соответствии с квадратичным критерием
Регулятор
Иs0sl
2so
п и
ПИД
И |
S0Si |
s 0 5 3 |
52 |
|
|
|
S2 |
|
ПИ |
J _ |
_ L |
|
9 |
S3 |
ПИД |
|
|
И |
|
|
|
|
_ 1 _ |
|
|
«3 |
ПИ |
|
|
ПИД |
|
sfs3 |
16 |
|
|
|
|
|
ПИ |
|
|
|
1 |
s 2 |
|
s |
4 |
ПИД |
• S2S5) |
u = r0+s0
г > > 0
•?1^2 sa
So/"
S2f > 53 И
— $ 2 М
3 |
s3 |
X |
X U 2 - |
S1S4 \ , |
|
, 1 |
Sl$4 |
|
9 |
' |
c? |
|
|
S2S3 |
16 |
s4 |
s4 |
_3_ |
|
_ 3 . . - £ | . ( 4 s 3 s 4 - |
|
16 |
s4 |
||
A 4 |
|||
|
|
||
|
|
• 3s2 s5 ) |
155