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

книги из ГПНТБ / Клюев, А. С. Автоматическое регулирование

.pdf
Скачиваний:
137
Добавлен:
22.10.2023
Размер:
17.2 Mб
Скачать

замыкании столбиком ртути контактов 2—3 - ртутного

 

электрокон-

тактного термометра Т-РК срабатывает реле /V

При

замыкании

кон­

тактов

1—2 срабатывает

реле

Р и прекращая

 

подачу

энергии

в объект

вследствие размыкания размыкающего контакта Р\. Одио-

 

 

 

 

временно реле

 

 

получает

пи­

 

 

 

 

тание

через

свой

замыкающий

 

 

 

 

контакт

и

ранее

замкнувшийся

 

 

 

 

замыкающий

контакт

реле

Яг.

 

 

 

 

Когда

 

температура

 

шкафа

 

 

 

 

уменьшается,

 

при

 

размыкании

 

 

 

 

контактов

2—3 отпадает

якорь

 

 

 

 

реле Я2, в результате чего реле

 

 

 

 

рі теряет питание и подача

 

 

 

 

энергии

в

объект

через

за ­

 

 

 

 

мкнувшийся размыкающий кон­

 

 

 

 

такт

Я 1

возобновляется.

З а ­

 

 

 

 

данное

значение

 

 

температуры

 

 

 

 

Ѳ3

находится

в

середине м еж ­

 

 

 

 

ду контактами 1 и 2. Темпера­

 

 

 

 

туру,

соответствующую делени­

 

 

 

 

ям шкалы термометра, на ко­

 

 

 

 

торыхустановлены электрокон­

 

 

6)

 

такты

1

и

2,

обозначим

соот­

 

 

 

ветственно

Ѳі

и

 

Ѳ2.

Следова­

 

 

 

 

тельно,

Ѳ з =

(Ѳі + Ѳг)/2.

 

 

Рнс. 7-14. Нелинейная АСР тем­

 

Пусть

передаточная

функ­

пературы

сушильного

шкафа

ция

сушильного

шкафа

имеет

с двухпозиционным релейным эле­

вид:

 

 

 

 

 

 

 

 

 

 

 

ментом

без зоны нечувствительно­

 

 

 

 

 

 

 

 

 

 

 

 

 

сти и ее структурная схема.

 

 

 

^ М

= Т Г Г ѵ

(?-п )

а нелинейная функция релейного элемента

м ож ет

быть

представле­

на как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пи= Я п(ДѲ).

 

 

 

 

 

 

 

 

 

 

 

 

Дифференциальное уравнение линейной части системы запишет­

ся как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сіАѲ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т dt

+ ДѲ = fcojÄU n ,

 

 

 

 

 

 

 

 

(7-12)

где Д Ѳ = Ѳ * —Ѳа — отклонение текущего значения температуры

Ѳ* от

заданного

Ѳ3; AUB=iUaU* — отклонение

напряжения

UB, фактиче­

ски подводимого к нагревательному элементу, от напряжения Ч3,

которое необходимо

для поддержания

заданной температуры

Ѳ3;

&об — коэффициент

передачи объекта

в абсолютных единицах;

Т

постоянная времени объекта.

П редположим, что напряжение, подводимое к нагревательному элементу, в 2 раза больше напряжения, которое требуется для обес­ печения заданного значения температуры, т. е. lh=2Ua.

Н иж е даны значения AUU в зависимости от величины отклоне­ ния температуры ДѲ

310

 

U„ =

и,

Ш ѵ = и з при

ѲК < Ѳ 2, т.

е. при

)

 

 

 

ДѲ:

Ѳ, — Ѳ.

 

Ѳ2 — н 3;

 

 

 

 

 

 

■=

 

 

 

Uи =

О,

ДUn = — U3 при

В* >

0 і ,

т. е. при

 

 

 

 

 

 

Ѳ2 — ѳ ,

 

 

 

 

 

 

 

Д 0 ^

-

- - 9 -

-

=

Ѳ, -

®3;

 

 

и и =

0;

Дия = — из при

Ѳ2 — Ѳ3 <

ДѲ <

Ѳ, — Ѳ3 и

(7-13)

 

 

 

 

 

 

 

dbSd

 

 

 

 

 

 

 

 

 

 

dt

< 0;

 

 

 

Uu=

и;

Ш п— иэ п р и

Ѳ2

Ѳ3 <;ДѲ <

© ! — Ѳ3 И

 

 

 

 

 

 

сгдѳ

° ‘

 

 

 

 

 

 

 

 

 

~ З Г >

 

 

 

 

 

Таким

образом

если

за входную величину

релейного

элемента

на рис. 7-14 примять отклонение температуры от заданной ДѲ, а за выходную — абсолютное значение <Ua напряжения, подводимого к нагревательному элементу, то характеристика релейного элемента является несимметричной и имеет вид, изображенный в правом ниж­ нем углу на рис. 7-10,г. Если же за выходную величину принять AU „, то характеристика релейного элемента будет симметричной без зоны нечувствительности (рис. 7-7,а) и с зоной неоднозначности

2 (-©г—Ѳ31= 2]Ѳ)—Ѳз I.

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

ния

на нагревательный

элемент. Разделив

уравнение (7-12) на

Ѳ уст, получим:

 

 

 

 

 

 

 

 

 

ДѲ

 

 

 

 

 

 

 

Т

Ѳ

 

ДѲ

 

U

A U n

 

 

dt

 

уст

^об а

' г/

 

 

 

 

 

в у і

и

 

ПЛИ

 

 

 

 

 

 

 

 

 

 

 

Т

dx.

 

 

 

 

 

 

 

dt + х

 

 

 

 

где

х= А Ѳ /Ѳ уСТ — текущее

относительное отклонение

регулируемого

параметра (температуры)

от заданного значения Ѳ3; k = k 0^Ujhyct

передаточный коэффициент

системы

в относительных

единицах.

21*

311

Подставив в это уравнение значения АUu/U из формул (7-13), получим систему дифференциальных уравнении:

dx

Т

+ х =

ки при л < — Ь\

 

 

dx

 

—- ки при

X ^ 6;

 

 

Т -JJ- -{- X =

 

 

_ dx

 

 

 

 

dx

(7-14)

+ х ~

ки при ~

ь < х < Ь

 

Т ~dt

и - j f > J J ;

 

dx

-f- X =

ки

при

 

dx

 

Т

b < X < b и -jjj < О,

 

где u—U3jU — относительное

отклонение

величин напряжения, по­

даваемого на нагревательный элемент, от напряжения, необходимого

для

поддержания заданного -значения

регулируемого

параметра;

2 6 = (Ѳі— Ѳ2) /Qyc- — зона неоднозначности релейного

элемента

в относительных единицах.

 

 

 

Переходный процесс в системе мож ет быть найден путем реше­

ния

дифференциальных уравнений (7-14)

по участкам, границы м еж ­

ду которыми определяют относительные переключающие значения выходной величины х (относительное отклонение температуры от заданного значения): Хі = — Ь, х2= Ь.

При решении дифференциальных уравнений на каж­ дом участке за начальные условия принимаются конеч­ ные значения переменных величин на предыдущем участке. Такой метод нахождения переходного процесса в системе путем решения дифференциальных уравнений на отдельных участках изменения исследуемого парамет­ ра X при условии согласования между собой значений са­ мого параметра х и его производных dx/dt в конце пре­ дыдущего и начале последующего участков называется

методом припасовывания или сшивания.

Характеристическое

уравнение системы,

изображенной

на

рис. 7-14, 771-1-1= 0 имеет

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

корень р = — 1/7.

Со­

гласно формуле (2-37) общее решение однородного дифференциаль­

ного уравнения (7-14) в таком случае имеет вид: __ і_

x(t) — Cte Т.

Частное решение неоднородного дифференциального уравнения при подаче напряжения на нагревательный элемент, т. е. при £/„=.£/,

можно найти как установившееся значение параметра, т. е. прини­ мая, что dxjdt—0. Тогда х п.Уст = ки.

Общее решение неоднородного дифференциального уравнения

складывается из суммы общего решения его

однородного уравнения

и частного решения неоднородного уравнения:

 

t

 

—У

+ ки.

X = X (/) + x B.yoi= C Iß

312

Обозначив х = — л'оіі при 1=0 и при Un = U, получим:

С 1= ----Von— k lL .

Таким образом, переходный -процесс в системе при \Ua=iU опре­ деляется выражением

 

t

___<_

 

х =

/е и ^ 1 — е

хопе

т,

(7-15)

которое справедливо

на участках,

где x=g:— b и — 6 < х < 6

при

dx/dt>0. Частное решение неоднородного дифференциального урав­ нения при U,1 = 0

Хоуст = —

Обозначив через Хоо начальное значение параметра х в диффе­

ренциальном уравнении при

Н н= 0 ,

аналогично

предыдущему

най­

дем выражение для переходного процесса:

 

 

 

__<_

__t_

 

x = ku^e

т

l j - f x 00e

г .

(7-16)

Этим выражением процесс регулирования определяется на уча­ стках, где х ^ Ь и — Ь<х<Ь при dxjdt<0.

 

Рис.

7-15. Переходный

процесс системы

 

регулирования, приведенной на рис. 7-14.

И з

уравнений (7-15) и (7-16) видно, что переходный процесс на

каж дом

участке

представляется как

сумма двух экспоненциальных

зависимостей.

Полная кривая процесса регулирования, построенная по урав­

нениям

(7-15) и .(7-16), представлена на

рис. 7

-15. Предположим, что

объект

включен в работу при

х=Х оп=Хоуст;

так

как при

этом

х < — Ь,

то на первом участке

процесс

регулирования определяется

выражением (7-15). Поскольку на первом

участке

кривой регули­

руемый

параметр будет возрастать, имеем

dx/dt>0,

поэтому

пере-

313

кліочемIіо рслеішого элемента произойдет только при .ѵ = 6. С этого момента начнется второй участок переходного процесса, который определяется выражением (7-16). Начальное значение параметра для второго участка ,ѵ0о=Ь. Так как в конце первого участка реле Рі срабатывает (см. рис. 7-14) и подача энергии в объект прекращается, то температура сушильного шкафа начнет уменьшаться и dxjdt<Si. Следовательно, процесс регулирования будет определяться выраже­

нием (7-16) до достижения регулируемым

параметром значения

х —

==— Ь. В этот момент согласно уравнениям

(7-14) переключается

ре­

лейный элемент, отключается реле Р\ и возобновляется подача энер­ гии в объект. С этого момента начнется третий участок, опреде­ ляемый, как II первый, выражением (7-15), но при начальном зн а­ чении параметра х Ь.

Третий участок закончится при ,ѵ=6.

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

Таким образом, начиная со второго участка в нелинейной си­

стеме по рис.

7-14

установятся незатухающие колебания,

которые

при

U=2U3 II

|Ѳ 2— ё>з! =

IѲі— Ѳ3|

определяются

выражениями

 

,

 

_J_V

 

t

 

 

 

 

'

 

 

 

 

 

Лх

 

 

 

 

 

 

 

 

 

 

X — ku 1 1- е

 

Л

be

т при — b <

X < b н

> 0 ;

 

 

 

 

 

 

t

 

 

~Tt

 

 

 

,

t

 

 

b > X

> — Ь и

dx

 

 

II

 

 

 

) + Ь е

т при

A О

 

 

 

 

 

 

 

 

 

dt

 

 

 

М етодом

 

припасовывания можно определить переходный про­

цесс

любой

нелинейной

системы

по участкам

м еж ду

соседними

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

к

участку изменяется структурная схема системы, а

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

и

динамические характеристики линейной системы,

эквивалентной

на данном участке исходной нелинейной.

 

6) Метод гармонической линеаризации

При подаче на вход линейной системы устойчивых гармонических колебаний *ВХ=А sin af на выходе си­ стемы также возникают гармонические колебания, но с другой амплитудой и фазой лгВых= ^і sin (co^+ ept). Отно­ шение выходной гармонической величины к входной, вы­ раженное в комплексной форме, определяет АФХ си­ стемы.

Если на вход нелинейного элемента подавать такие же синусоидальные колебания, то на его выходе также получим периодические колебания, но по форме сущест-

314

венно отличающиеся от синусоидальных. В качестве при­ мера на рис. 7-16 показан характер изменения выходной переменной хвых нелинейного элемента с релейной ха­

рактеристикой (рис. 7-6,а) при поступлении

на. его вход

синусоидальных колебаний

 

xBX—As\ nat .

(7-18)

Как известно, любая периодически изменяющаяся пе­ ременная величина путем разложения ее в ряд Фурье может быть представлена как сумма некоторой средней

постоянной и бесконечного множества гармонических со­ ставляющих

 

HnSinmrf + BnCOSflmf),

(7-19)

 

Я = 1

 

где Л0,

Лі, ßi, Л2, В2, . . . , АПі Вп — постоянные

коэффи­

циенты

ряда Фурье; <и = 2я/Г — частота колебаний пер­

вой гармоники (основная частота), равная частоте вход­

ных'синусоидальных колебаний; Т — период

колебания

первой гармоники, разный периоду входных

синусои-

315

далы-іых колебаний; 2ш= со2 , Зсо= саз,. . /гсо= со„ — часто­ ты, кратные основной частоте to.

Автоматические системы регулирования, как правило, обладают существенной инерционностью. Так если в си­ стеме на рис. 1-5 при установившемся значении темпе­ ратуры медленно перемещать движок трансформатора по синусоидальному закону, то температура объекта также будет изменяться по синусоидальному закону с той же частотой. Но если движок перемещать по синусоидаль­ ному закону с большой частотой, то в силу инерцион­ ности объекта его температура не сможет та’к быстро изменяться, а будет иметь постоянную величину, соот­ ветствующую среднему положению движка. Это эквива­ лентно тому, как если бы движок вообще не перемещал­ ся. Иными словами, входная величина с большой часто­ той практически не проходит через инерционную систему и не оказывает на нее существенного влияния. Таким образом, линейная часть системы является фильтром по отношению к высокочастотным гармоническим состав­ ляющим. В связи с этим, а также учитывая, что ампли­ туда выходных колебаний уменьшается с ростом часто­ ты, для приближенной оценки выходной величины, рас­ сматриваемой как спектр гармоник, в большом числе случаев достаточно учесть только первую ее гармониче­ скую составляющую.

С учетом этого при отсутствии постоянной составляю­ щей в выходных колебаниях выражение (7-19) прибли­ женно можно записать в виде

Хвы-а=Аі sin at + Bi cos at.

Так как производная входной величины dxTafdt= A a cos соі, то, подставляя значения cos at в (7-20), получаем:

__Ai

I В1

dx„x

A -^вх-Г Aa

dt ■

(7-20)

(7-18) sin соt и

(7-21)

Таким образом, нелинейная зависимость выходной величины от входной в нелинейном элементе приближен­ но заменяется линейной зависимостью, определяемой вы­ ражением (7-21).

С учетом этого отношение изображения выходной ве­ личины к изображению входной величины

31 6

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

Передаточная функция 1Уц(р, со, А) нелинейного эле­ мента имеет существенное отличие от передаточной функ­ ции линейной системы W(p), заключающееся в том, что Wn(p, со, А) зависит от амплитуды и частоты входного сигнала, так как выходная величина нелинейного эле­ мента (7-21) зависит от амплитуды и частоты входного сигнала.

Обозначив Аі/А*=Ь(А) и Ві/В = с(А) , запишем:

W* (р, ю, А) = b (А)'+ с (А).

(7-22)

Заменяя в (7-22) р на /со, по аналогии с линейными системами получим выражение для АФХ нелинейного элемента:

WB(A,jv>)=b(A)+jc{A). (7-23)

В выражении (7-23) величины b (А) и с (А) определя­ ются как отношение коэффициентов ряда Фурье для пер­ вой гармоники выходных колебаний к амплитуде вход­ ных колебаний.

Их значения для наиболее распространенных типовых нелинейностей приведены в Приложении 3.

7-3. ИЗОБРАЖЕНИЕ ПРОЦЕССА РЕГУЛИРОВАНИЯ НА ФАЗОВОЙ ПЛОСКОСТИ

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

Фазовой плоскостью называется плоскость, на кото­ рой изображается изменение какой-либо переменной ве­ личины X в функции скорости ее изменения y —dxjdt.

Рассмотрим изображение незатухающих гармониче­ ских колебаний x=y4sin<a£ на фазовой плоскости.

Скорость изменения переменной х выразится как

у = ‘^ г = Аю cos аit,

317

Исключав переменную t, найдем:

у = Аш cos arcsin -А- — Дш cos arccos

откуда

х2

Дг

1.

(7-24)

Д2

Д2со2

 

 

Следовательно, незатухающие гармонические колеба­ ния изображаются на фазовой плоскости в виде эллипса (рис. 7-17) с полуосями А и /4©.

Рис.

7-17. Н езатухаю щ ее гармоническое колебание

и его

изображение на фазовой плоскости.

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

При заданной амплитуде А для частот, изменяющих­ ся от 0 до ± о о , получим семейство эллипсов с общей осью по оси координат х и изменяющейся длиной оси 2Лсо; при сй = 0 эллипс вырождается в прямую, совпадаю­ щую с осью 2 Д ; при © = ± о о он распадается на две пря­ мые, параллельные оси координат у и отстоящие от на­ чала координат па расстояние ±А.

При изменении величины х соответствующая ей точка М на фазовой плоскости, называемая изобраоісающей точкой, будет перемещаться тздоль эллипса по часовой стрелке с периодом колебания Тк.

Кривая движения изображающей точки на фазовой плоскости называется фазовой траекторией.

Если колебательный процесс будет расходящимся (рис. 7-18), то амплитуда колебаний А будет увеличи­ ваться и изображение процесса на фазовой плоскости бу­ дет иметь вид расходящейся логарифмической спирали.

Наоборот, затухающий колебательный процесс

318

(рис. 7-19) на фазовой плоскости изображается в виде логарифмической спирали, сходящейся к началу коорди­ нат. Таким образом, по виду фазовой картины (фазового портрета) можно наглядно судить об устойчивости си­ стемы.

Рис.

7-19. Затухаю щ ее

колебание

и его изображ ение на

фазовой

плоскости.

 

 

Так

как

в верхней

половине

фазовой плоскости у =

= dxjdt> 0, то изображающая точка вдоль фазовой тра­ ектории движется в сторону увеличения х. В нижней по­ ловине у < 0, следовательно, изображающая точка дви­ жется вдоль фазовой траектории в сторону уменьшения

X . Так как на

фазовой плоскости в пересечениях фазовых

траекторий с

осью х производная

y=dxfdt = 0, то фазо­

вые траектории пересекают ось х

под прямым углом.

В качестве примера построим фазовую картину нелинейной АСР температуры сушильного шкафа, изображенной на рис. 7-20. Система содерж ит регулятор с нелинейным элементом, имеющим однозначную релейную характеристику с зоной нечувствительности (см. рис. 7-5,о ). Входной величиной релейного элемента является отклонение х температуры в объекте от заданной, а выходной — на­

пряжение 'І7д, подаваемое на электродвигатель Д.

319

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