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

книги / Теория автоматического управления

..pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
16.96 Mб
Скачать

(9.47)

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

D = 1— 4Ab2.

(9.48)

Уравнение (9.47) имеет действительные корни х2та (а в системе соот­

ветственно возможны автоколебания с амплитудой хта) только в том случае, если дискриминант D > 0, т. е. если

2kj\cTQTд л (T о Т д) b 0.

(9.49)

При заданных в условии примера параметрах неравенство (9.49) не вы­ полняется, и следовательно, автоколебания в системе не могут возникнуть.

Из условия (9.49) можно определить предельное значение коэффициента линейной части

kji. пр — Tib {То _1“ Тд)/2сТ0Тд.

(9.50)

При большем значении коэффициента линейной части автоколебания возможны. Например, при kn = 0,4 уравнение (9.47) имеет два действитель­ ных решения: хта = 2,1 °С и хта = 12 °С.

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

ОР/дьу |(о=а)а = - 2 (Т0 -|- 7д)

соа;

(9.51)

 

dQldxm = О; до/дс* |ш=(0а = 1

- ЗГ0Гдсо2.

 

Подставляя значения частных производных (9.51) в условие (9.38) и учи­ тывая одновременно выражение (9.45) для частоты, получим

 

4клс (2Ь2

,) (1 -3 )/ю 4 ,а д Д т а - 62 > °

(9.52)

или

 

 

 

 

 

 

Х т а > Ь

У

2 .

 

 

(9.53)

 

Этому

условию удовлетворяет только второе решение. Следовательно,

в системе будут возникать устойчивые

автоколебания

с амплитудой хт& =

=

12 °С и

частотой соа =

\1л/Т0Тп =

0,45 с-1. Период колебаний Га =

=

2я/0,45 =

14 с.

 

 

 

 

Покажем, что при отсутствии зоны нечувствительности реле в системе

всегда будут

возникать

устойчивые

автоколебания.

Действительно, при

b = 0 уравнение (9.47) всегда имеет одно действительное решение

 

Хта — 4&JJC7*оТ д/я Q

7 д).

 

(9.54)

Нетрудно видеть, что при 7д 0 амплитуда автоколебаний будет стре­ миться к нулю, а их частота согласно выражению (9.45) — к бесконечности.

Пример 2. Определим амплитуду и частоту автоколебаний системы, состоящей из трехпозиционного реле (см. № 2 в табл. 9.1) и линейной части с запаздыванием

(9.55)

321

Параметры

системы:

/гл = 1,

Т — 10 с, т — 5 с, с — 25,

Ъ — 1.

Поскольку

система

обладает

запаздыванием, применим

графический

способ решения.

На рис. 9.10, б построены а. ф. х. линейной части и обратная характе­

ристика реле

 

 

 

(9.56)

Характеристика (9.56) дважды проходит вдоль

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

Р (©). При

xm = Ъ — 1 характеристика начинается

от — с», при хт =

= b У2 «

1,4 достигает максимального значения я6/2с, а затем при хт ->-оо

опять стремится к — оо.

Характеристики линейной и нелинейной частей пересекаются в двух

точках. Согласно приведенному выше правилу точка

(на верхней

ветви

обратной характеристики) соответствует неустойчивым

колебаниям, а

точка

М 2 (на нижней ветви) — устойчивым. Амплитуда автоколебаний прибли­ зительно равна хта « 8, а частота (оа « 0,375 с"1.

9.5. Оценка абсолютной устойчивости с помощью критерия Попова

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

Пусть в контуре нелинейной системы (см. рис. 9.2) содержится нелинейный элемент с характеристикой f (хн), имеющей любую конфигурацию, но не выходящей за пределы определенного сектора [0, kH] (рис. 9.11, а). Состояние равновесия нелинейной системы называется абсолютно устойчивым, если оно асимптотически устой­ чиво (см. гл. 5) при любой нелинейности, относящейся к определен­ ному классу. Нелинейности считаются одного класса, если их ха­ рактеристики f (хн) располагаются в секторе между осью абсцисс и прямой с угловым коэффициентом kn (см. рис. 9.11, а).

На первый взгляд, может показаться, что для оценки абсолют­ ной устойчивости нелинейной системы, состоящей из линейной части с Wn {р) и нелинейности с f (хн), достаточно оценить устойчи­ вость линейной системы с передаточной функцией kHWn (/?). Од­ нако это предположение выполняется только для некоторых част­ ных видов Wn (р)у а в общем случае оно несправедливо, и требуются специальные критерии.

Удобный критерий для суждения об абсолютной устойчивости нелинейных систем предложил в 1959 г. румынский уче­ ный В.—М. Попов. Его критерий основан, как и критерий Най­ квиста, на использовании а. ф. х. и имеет простую геометрическую трактовку.

Рис. 9.11. Критерий абсолютной устойчивости Попова

В формулировке критерия используется понятие модифициро­ ванной а. ф. х. Пусть линейная часть устойчива и имеет следующую а. ф. х.

Г л (/со) = Р(со) + /(Э(со).

(9.57)

Образуем из этой обычной а. ф. х. следующую видоизмененную

а. ф. х.

 

 

Гл ( /с о ) = Р (с о )

+ /co7MQ ( со) ,

(9.58)

у которой мнимая

часть получена

умножением Q (со) на соГм, где

Ты = 1

с — нормирующий множитель. Характеристика (9.58) и на­

зывается модифицированной.

 

к р и т е ­

Теперь можно

следующим образом сформулировать

р и й

а б с о л ю

т н о й у с т о й ч и в о с т и

равновесия

нели­

нейной

системы,

которая состоит из линейной

части

с а.

ф. х.

1РЛ(/со) и нелинейного элемента с характеристикой / (хн),

распо­

ложенной в секторе [О, Ан]:

 

чтобы

для абсолютной устойчивости равновесия достаточно,

модифицированная характеристика

(/со) не охватывала

точку (— Uku, 0) и через эту точку можно было провести прямую, не пересекающую характеристику W'n (/со) (последняя

лежит справа от прямой).

На рис. 9.11, б показан случай, когда критерий Попова выпол­ няется, а на рис. 9.11,6, г — случаи, когда не выполняется.

С помощью критерия Попова решают и обратную задачу: строят заданную характеристику W* (/со), затем проводят как можно

ближе к этой характеристике прямую так, чтобы получить наимень­ ший отрезок [0; — 1/&н], и таким образом находят допустимое зна­ чение углового коэффициента /гн.

По наклону прямой Попова, «прижатой» к кривой W*n (/о),

можно судить о допустимом классе нелинейности: если прямая вер­ тикальна, то нелинейность может быть только однозначной, а если она наклонена, то нелинейность может быть и однозначной и не­ однозначной (с гистерезисом).

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

UM/ ц ) - — -/й))

________________

(9.59)

( 7 \/ с о + 1) (7\,/со + 1)

D„ (/ш)

( Г л/со + 1)

модуль которой с ростом со монотонно убывает. Исследование абсо­ лютной устойчивости таких систем можно свести к анализу устой­ чивости соответствующей линейной системы: допустимое значение углового коэффициента kH находить через предельное значение общего передаточного коэффициента (см. 5.6).

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

элемента (рис. 9.12, а) с параметрами

b =

1, с = 8 и линейной части

(/со) =

*л е—/£0Т/(Т/со +

1)

 

(9.60)

с параметрами

&л = 0,25, Т =

10 с,

т =

10 с.

а

Представим а. ф. х. (9.60) в виде суммы действительной и мнимой частей:

 

(/<»)=

/<л (cos (от — Т ш sin ют)

.

кл (sin сот + Ты cos сот)

(9-61)

 

 

 

Т 2со2 Н- 1

1

Т 2ш2 +

1

'

Соответствующая (9.61) модифицированная характеристика (рис. 9.12, б)

 

/ .,л

1<п (cos сот — Т со sin сот) .

кл (sin сот +

Гео cos сот)

(!Ш>

v -

 

---------- F w T i----------'

-------- №

f1 ------

Через

точку

с абсциссой, равной \J k „ =

0,125, проведена

прямая

/,

ко.

торая

не

пересекается с кривой W n (/со). Следовательно,

при

заданных

па-

раметрах

равновесие системы абсолютно

устойчиво.

 

 

 

 

Решим теперь обратную задачу: определим допустимое по условию устой­ чивости равновесия значение зоны нечувствительности Ъ. Для этого прове­

дем прямую 2, «прижатую» к характеристике

(/со). Она пересекает дейст­

вительную ось в точке

с абсциссой «0,11. Отсюда допустимое значение для

углового

коэффициента kn ^

1/0,11 « 9, а для

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

Ь > elkн =

8/9 « 0,9.

При

b < 0,9 состояние

равновесия системы будет

неустойчивым. -

9.6. Основные сведения о вибрационной и статистической линеаризации

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

Изменение характеристики линейной части достигается спосо­ бами, изложенными в гл. 7.

Изменение характеристики нелинейного элемента может быть получено несколькими способами. Некоторые статические нели­ нейности можно компенсировать при помощи соответствующих обратных нелинейностей. Для этого параллельно или последова­ тельно с основной нелинейностью f (х„) включают компенсирующую нелинейность, имеющую обратную характеристику f " 1 (хн). Экви­ валентное соединение при этом будет линейным.

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

Рассмотрим этот способ подробнее.

Вибрационная линеаризация. Вибрационная линеаризация

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

Рассмотрим математическую сущность вибрационной линеари­ зации и проиллюстрируем получение линеаризующего эффекта на примере вибрационной линеаризации идеального реле.

Пусть на входе нелинейного элемента с однозначной характе­ ристикой f (хн), включенного в контур системы (см. рис. 9.6, б), действует медленно меняющийся основной сигнал х0 (/) и допол­

нительный высокочастотный сигнал хд (t) с амплитудой хт

и ча­

стотой сод (рис. 9.13, а), т. е.

 

хн (0 = х0(t) + хд (0 = х0 (t) + хдт sin (o„t.

(9.63)

При этом частота сод должна быть достаточной для того, чтобы основной сигнал x0{t) в течение нескольких периодов дополни­ тельного сигнала хд (t) оставался практически постоянным (см. рис. 9.13, а), а амплитуда хДт — достаточной для срабатывания

Рис. 9.13. Вибрационная линеаризация идеального реле

326

релейного элемента (рис. 9.13, б), т. е. хЛт> \х0 (/)|. Кроме того, необходимым условием работоспособности такой системы является подавление линейной частью составляющей с частотой сод и более высоких гармоник:

I Wn (/©д) |< | W„ (/©„) |,

(9.64)

где (Оп — частота

пропускания линейной части.

 

Сигнал у„ (t)

на выходе нелинейного элемента

также будет

представлять собой сумму медленно меняющейся составляющей y0(t) (рис. 9.13, в), равной среднему значению сигнала уя (/), и быстро меняющейея уд (t):

Ун (0 = Уо (О + Уя (0.

(9.65)

где Уо (0 = /о [*о (t) 1, yA (t) = /д

[хд(0 ], причем функции /0 Ф

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

Функции /о и /д можно приближенно определить, если разло­ жить суммарный выходной сигнал yH(t) в ряд Фурье (2.18). При этом полагают, что на интервале разложения, равном периоду Т сигнала хд (t), составляющая y0(t) практически не меняется, т. е. y0(t) « const. Тогда по аналогии с основной процедурой гармони­ ческой линеаризации [см. формулу (9.27) 1 можно вместо нелиней­ ной зависимости

yn{f) = f[Xo, *я (01

(9-66)

записать

гармонически линеаризованное

уравнение:

Ун (t)

« Уо + Ц (ХЛ т) Хя (0.

(9.67)

Коэффициент g = fя определяется аналогично — по формуле (9.26). Составляющая у0 соответствует нулевому члену разложе­ ния периодического сигнала ун (t) в ряд Фурье, т. е.

О

7/2

(9.68)

|/о= = - | -

s f(x0, xA(t)] d t

т-772

Очевидно, что этот член будет зависеть только от вида функции f, медленной составляющей х0 и от амплитуды Хят высокочастотного сигнала, т. е.

£{/о= F (х0, Хят) ИЛИ y0 = F o {Хо).

(9-69)

Так, например, для рассматриваемого двухпозиционного реле без зоны нечувствительности (см. рис. 9.13,6) зависимость (9.69) принимает вид

я arcsin

*0

(9.70)

*дт

 

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

Для большинства однозначных нелинейностей график функции F0 (х0) имеет достаточно плавную форму (даже для релейных эле­ ментов). На рис. 9.13, г показана характеристика F0 (х0) идеаль­ ного реле при различных амплитудах хлт. Эффект сглаживания нелинейностей для медленной составляющей, достигаемый наложе­ нием высокочастотных колебаний, называется вибрационной линеа­

ризацией.

Плавность функции F0 (лг0) позволяет осуществить также обыч­ ную линеаризацию

у0 ж кпх0,

(9-71)

где kH= (dFoldxo)Xo=0 — передаточный коэффициент

нелиней­

ного элемента для медленной составляющей. Например, для двух­

позиционного реле передаточный

коэффициент

ft,, = 2с/лхЛт.

(9.72)

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

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

ний. При соблюдении определенных условий выбора

частоты сод

и амплитуды х ^ система под воздействием сигнала

хд (t) может

перейти из режима автоколебаний с частотой <оа в режим вынуж­ денных колебаний с частотой <ад. Это явление называется захваты­ ванием частоты. Так как вынужденные колебания имеют обычно меньшую амплитуду и большую частоту, чем автоколебания, то такое устранение автоколебаний также улучшает качество нели­ нейных систем.

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

Одним из вариантов осуществления вибрационной линеариза­ ции за счет внутренних свойств системы является создание так на­ зываемого скользящего режима.

Вибрационная линеаризация в скользящем режиме. Проиллю­

стрируем сущность вибрационной линеаризации в скользящем режиме на примере релейной системы стабилизации температуры (см. рис. 9.1).

Рис. 9.14. Фазовый портрет скользящего режима

Для простоты будем полагать, что зона нечувствительности реле b = 0 и постоянная времени двигателя Тл = 0. Согласно схеме (см. рис. 9.6, б) входной сигнал нелинейного элемента (с учетом внутренней обратной связи)

 

x » = — [l + ko.c(ToP+\)/kJio]x

(9.73)

или

после

преобразований

 

 

хн =

— ki (1 + Тгр) х,

(9.74)

где

к ^ =

1

( к0. с/k3kQ\ Т j —- k0. cT0/{k0. с

“I- k3k0).

Учитывая выражение (9.74), можно записать дифференциальное уравнение замкнутой системы:

T°Jtr+i r - - k"l(lh*+k'T'Jfr)-

( 9 - 7 5 )

В системе без внутренней обратной связи линия переключения на фазовой плоскости описывается выражением хп = 0 и совпадает с осью ординат (рис. 9.14, а). Если на вход реле кроме основного сигнала х подается дополнительный сигнал, пропорциональный производной от х у то уравнение линии переключения для рассмат­ риваемого идеального реле без зоны нечувствительности имеет вид

T l -±*- + x= Q .

(9.76)

d t

 

В системе координат «хх» уравнению (9.76) соответствует пря­ мая M N (рис. 9.14, б), имеющая наклон 1/7\.

Фазовые траектории по обе стороны линии переключения опре­ деляются такими же уравнениями (9.17), как и в системе без внут­ ренней обратной связи. Но на наклонной линии переключения всегда образуется особый отрезок АВ , который в своих крайних точках касается двух фазовых траекторий (показаны штриховыми линиями). Вне этого отрезка фазовые траектории подходят к линии переключения с одной стороны, отходят с другой, а на самом от­ резке АВ — с обеих сторон. Это означает, что изображающая точка, попав на отрезок АВ, не сможет «уйти» с него, а так как производ­

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

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

x(t) = x0e - ^ T\

(9.77)

Убывание сигнала х по экспоненте означает, что в скользящем режиме нелинейная система подчиняется линейным законам, и ее можно приближенно рассматривать как линейную систему первого порядка. К этому выводу можно прийти и путем следующих рассуждений. Средняя составляющая сигнала хн в скользящем режиме приблизительно равна нулю, а среднее значение сигнала х конечно. Поэтому передаточный коэффициент прямой цепи на рис. 9.6, б можно считать равным бесконечности. Соответственно передаточная функция замкнутой системы в скользящем режиме согласно основ­ ному свойству предельной системы [см. формулу (2.145) ] будет равна обратной передаточной функции звена обратной связи:

ф (Р) = X (р)1хз(р) « (l/k1)/(T1p + 1).

(9.78)

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

Статистическая линеаризация. Линеаризующий эффект в не­

линейных системах получается также при воздействии на них высокочастотных случайных сигналов.

330