- •ПРЕДИСЛОВИЕ
- •ВВЕДЕНИЕ
- •1. НЕЛИНЕЙНЫЕ АВТОМАТИЧЕСКИЕ СИСТЕМЫ УПРАВЛЕНИЯ
- •1.1. Основные сведения о нелинейных системах
- •1.2. Структура нелинейной системы
- •1.4. Типовые нелинейные звенья
- •1.6. Примеры нелинейных автоматических систем
- •Контрольные вопросы и задания
- •2. АНАЛИЗ НЕЛИНЕЙНЫХ СИСТЕМ УПРАВЛЕНИЯ
- •2.1. Основы анализа нелинейных систем
- •2.2. Фазовое пространство
- •2.3. Фазовые траектории систем
- •2.4. Метод изоклин
- •2.8. Метод точечных преобразований
- •2.9. Гармоническая линеаризация нелинейностей
- •2.11. Метод гармонического баланса
- •2.12. Запаздывание гистерезисного реле
- •2.13. Частотный метод Попова
- •2.14. Прямой метод Ляпунова
- •2.15. Определение границ абсолютной устойчивости нелинейной системы
- •2.16. Анализ одночастотных вынужденных колебаний
- •Контрольные вопросы и задания
- •3. КАЧЕСТВО РЕГУЛИРОВАНИЯ НЕЛИНЕЙНЫХ СИСТЕМ УПРАВЛЕНИЯ
- •3.3. Показатели качества регулирования нелинейных систем
- •3.4. Анализ качества регулирования одноконтурной нелинейной системы
- •Контрольные вопросы и задания
- •4. КОРРЕКЦИЯ НЕЛИНЕЙНЫХ СИСТЕМ УПРАВЛЕНИЯ
- •4.1. Линейная коррекция нелинейных систем
- •4.2. Компенсация нелинейной характеристики
- •4.4. Вибрационная линеаризация релейного элемента
- •4.5. Псевдолинейная коррекция
- •4.7. Нелинейные системы с переменной структурой
- •Контрольные вопросы и задания
- •ЗАКЛЮЧЕНИЕ
- •ПРИЛОЖЕНИЯ
При гармонической линеаризации интервал значений 0 ≤ q(A) ≤ ∞ в действительности имеет место только для идеального реле. Для нелинейной характеристики, лежащей в интервале угла, значение q(A) заключено в пределах
q(A)min ≤ q(A)≤ q(A)max . |
(2.202) |
В таком случае дополнительную область устойчивости определяют в зависимости от максимального и минимального значений коэффициента гармонической линеаризации для каждого вида нелинейной характеристики.
Вначале границу области абсолютной устойчивости определяют при 0 ≤ q(A) ≤ ∞ согласно (2.191). Это решение дает достаточные усло-
ческой линеаризации. Затем получают дополнительную область ус-
вия устойчивости при любой форме однозначной нелинейной характеристики и при любом интервале измененияИкоэффициента гармони-
тойчивости за счёт ограничения интервалаДфактически возможного
изменения коэффициента гармонической линеаризации (учет границ
нелинейности). Для этого в уравнения (2.191) надо подставить значе-
и достаточные условия устойчивостиАпри заданной форме однозначной нелинейной характеристикии . Кроме этого, по условию (2.189) при различных варьируемых параметрах линейной части можно определить область устойчСвого автоколебательного режима (см. рис. 2.16).
ния q(A)min и q(A)max. Чем более узкий сектор, где разместилась вся не-
линейная характеристика, тем больше будет область дополнительной границы абсолютной устойчивостиб. Такое решение дает необходимые
2.16. Анализ одночастотных вынужденных колебаний
Исследование вынужденных колебаний нелинейных систем представляет собой в общем сложную задачу в связи с отсутствием суперпозиции отдельных решений, а также существенным изменением поведения системы в зависимости от размера амплитуды колебаний, с наличием не единственного установившегося режима и возможностью перескоков с одного режима на другой, с особенностями высших гармоник, субгармоник, комбинационных частот и с многими другими факторами.
Рассмотрим случай одночастотных вынужденных колебаний, т.е. колебаний нелинейной системы с частотой внешнего периодического воздействия, и найдем условия их существования.
92
Форма колебаний, как и прежде, на основании свойства фильтра линейной части системы будет считаться близкой к синусоидальной для переменной х, стоящей под знаком нелинейной функции. При рассмотрении вынужденных колебаний во многих случаях возникают ограничения, накладываемые на амплитуду и частоту внешнего периодического воздействия (зависящие также и от параметров системы) и обусловливающие существование одночастотных вынужденных колебаний в нелинейной системе. Будем их кратко называть условиями захватывания (в указанном широком смысле). Особое значение эти условия приобретают для автоколебательных систем при частотах, близких к частоте автоколебаний и выше.
Пусть нелинейная АСУ (см. рис. 1.1) с одной нелинейностью
F(x) имеет передаточную |
функцию линейной части, обладающую |
||||||
свойством фильтра низких частот. |
|
|
И |
|
|||
|
|
|
|
||||
Рассмотрим влияние на систему внешнего возмущающего гар- |
|||||||
|
|
|
|
Д |
|
||
монического воздействия f с известной амплитудой B' и частотой ω': |
|||||||
|
|
|
f (t) = B |
sin ω t . |
|
(2.203) |
|
|
|
|
′ |
′ |
|
|
|
|
|
А |
|
|
|
||
Уравнения линейной части системы и нелинейного элемента в |
|||||||
операторной форме: |
|
б |
|
|
|
|
|
|
Q (p)y = R(p)z + S (p) f ; |
(2.204) |
|||||
|
z = F(x) = −F( y). |
|
|
||||
Уравнение д нам ки замкнутой системы в операторной форме |
|||||||
при u(t) = 0 имеет в д |
|
|
|
|
|
|
|
Qи(p)y = −R(p)F(y) + S (p) f . |
(2.205) |
||||||
Возможные вынужденные колебания выходной величины с час- |
|||||||
тотой ω', амплитудойСА' и фазой φ' запишем как |
|
||||||
|
|
y = А′sin (ω′t +ϕ′). |
(2.206) |
||||
После гармонической линеаризации нелинейного элемента его |
|||||||
уравнение в операторной форме имеет вид |
|
|
|||||
|
|
|
|
′ |
p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′ y , |
(2.207) |
||
|
|
z = − q(A) +q (A) |
|||||
|
|
|
|
|
ω |
|
где коэффициенты гармонической линеаризации q(A) и q'(A) определяются согласно формулам (2.75) и (2.78).
93
Тогда уравнение замкнутой системы (2.205) с учетом формул
(2.203), (2.206) и (2.207) примет вид:
|
|
p |
|
|
|
Q(p)+ R(p) q(A) +q′(A) |
|
|
А′sin (ω′t +ϕ′)= S(p)B′sin ω′t. |
(2.208) |
|
′ |
|||||
|
|
ω |
|
|
Для отыскания периодического решения вынужденных колебаний выходной величины используем символический метод определения, перейдя в частотную область и заменив p = jω', а вместо sinω't подставив выражение ejω't. Тогда получим
Q ( jω )+ R( jω )(q(A) + jq (A)) |
А e |
j(ω t+ϕ ) |
= S( jω )B e |
j(ω t) |
; |
(2.209) |
|||||||||||||||||
|
′ |
|
|
′ |
′ |
|
′ |
|
|
|
|
′ |
|
|
′ |
|
|
||||||
|
|
|
|
|
|
|
|
′ |
|
′ |
|
|
|
|
|
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Q( jω )+ R( jω )(q(A) + jq (A)) |
|
|
|
|
|
|
j(ω t) |
|
|
|
|||||||||||
|
|
|
|
′ |
′ |
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
′ |
′ |
e |
|
′ |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
′ ; |
|
(2.210) |
||||||
|
|
|
|
|
S ( jω ) |
|
|
|
|
|
А |
= B |
|
e |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
j(ω t+ϕ ) |
|
|
|
||||||
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
Q ( jω )+ R( jω )(q(A) + jq |
|
И |
|
|
|
|
|
|
|
|||||||||||
|
|
|
(A)) |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
′ |
′ |
|
|
′ |
|
|
|
′ |
|
′ |
|
− jϕ′ |
. |
|
|
(2.211) |
|||
|
|
|
|
|
′ |
|
|
|
|
|
|
|
А = |
B e |
|
|
|
|
|||||
|
|
|
|
|
S ( jω ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Левая |
часть (2.211) обозначается |
|
Z(A) и представляет собой |
|||||||||||||||||||
|
|
|
|
|
А |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
комплексный годограф замкнутой системы с переменным параметром |
амплитуды А, определяемый при постоянной частоте внешнего гар- |
|||||||||||
монического воздействия ω'. |
|
Д |
|
|
|
||||||
Выражение (2.211) |
|
остается справедливым и для нелинейных |
|||||||||
|
|
|
и |
|
|
|
|
|
|||
систем с временным запаздыванием τ и принимает вид |
|
||||||||||
|
Q |
( jω )+e |
− jω τ |
R( jω )(q(A) + jq (A)) |
′ ′ − jϕ′ |
|
|||||
|
|
′ |
|
|
′ |
′ |
|
|
|
||
|
|
|
|
′ |
б |
|
|
|
|
||
|
|
С |
|
|
А = B e . |
(2.212) |
|||||
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
S ( jω ) |
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
|
|
|
Уравнение (2.211) с двумя неизвестными А' и φ' можно решить графически, как показано на рис. 2.38. Правая часть (2.211) изображается на комплексной плоскости [+1; j] в виде окружности радиуса В, а левая часть Z(A) строится как кривая по точкам при постоянной частоте ω' и различных значениях амплитуды А. Точки пересечения окружности с кривой Z(A) дают периодическое решение, причем величина амплитуды вынужденных колебаний определяется в точке пересечения по кривой Z(A), а фаза – по величине угла φ' (см. рис. 2.38).
На рис. 2.38 окружность пересекает кривую Z(A) только при радиусе, большем некоторого порогового значения В > Вкр. Следовательно, в этом случае одночастотные вынужденные колебания (2.206) возможны только при достаточно большой амплитуде В, а при меньшей амплитуде В внешнего воздействия будет иметь место сложное движение, включающее в себя и собственную частоту системы.
94
jV(ω)
ВВкр
0 |
U(ω) |
A |
φ' |
Z(A)
Рис. 2.38. Графическое определение параметров вынужденных колебаний
Построив серию кривых Z(A) по формуле (2.211) для разных значений частоты внешнего гармонического воздействия ω' (рис. 2.39, а),
можно определить зависимость порогового значения Вкр от ω' (рис. |
||||||
2.39, б), где ωa – частота автоколебаний данной системы. Тогда мы |
||||||
получим область значений В и ω', в которой существуют одночастот- |
||||||
ные вынужденные колебания, т.е. область захватыванияИ |
. |
|||||
jV(ω) |
|
а |
|
б |
Область |
|
|
|
|
Д |
|||
Вкр3 |
|
|
АВ |
захватывания |
||
|
|
|
|
|||
|
|
|
|
|
||
|
Вкр1 |
бU(ω) |
|
|
||
|
|
|
|
|||
0 |
и |
|
|
|
|
|
|
|
|
|
|
||
Вкр2 |
С |
|
|
|
|
|
|
|
|
|
|
|
|
|
Z3(A) |
|
Вкр(ω') |
|
|
|
Z1(A) |
Z2(A) |
0 |
|
ωа |
ω' |
|
|
|
Рис. 2.39. Определение области захватывания
Явление захватывания состоит в том, что при В > Вкр собственные колебания (автоколебания) срываются и система переходит целиком на одночастотные вынужденные колебания с частотой внешнего воздействия ω'. Данные одночастотные вынужденные колебания в соответствии со свойством фильтра низких частот линейной части системы для переменной х будут только близки к синусоидальным.
95
Рассмотрим анализ одночастотных вынужденных колебаний на |
|||||||||||||||
примере нелинейной АСУ (см. рис. 1.1) с нелинейной частью, обла- |
|||||||||||||||
дающей типовой нелинейной характеристикой F(x) «Двухпозицион- |
|||||||||||||||
ное реле с гистерезисом» (см. табл. 1.1, п. 8) с параметрами Bр и bр. |
|||||||||||||||
Уравнения линейной части системы и нелинейного элемента в |
|||||||||||||||
операторной форме: |
|
|
|
|
|
|
|
|
|
|
|
|
|||
(T p +1)(T p +1)py = Kz +(T p +1)(T p +1)pf ; |
(2.213) |
||||||||||||||
1 |
|
|
2 |
|
|
|
|
|
1 |
|
2 |
|
|
||
z = F(x) = −F( y). |
|
|
|
|
|
|
|
|
|
|
|||||
Уравнение динамики замкнутой системы в операторной форме |
|||||||||||||||
при u(t) = 0 имеет вид |
|
|
|
|
|
|
|
|
|
|
|
||||
(T1 p +1)(T2 p +1)py = −KF( y) +(T1 p +1)(T2 p +1)pf . |
(2.214) |
||||||||||||||
Пусть внешнее возмущающее гармоническое воздействие f на |
|||||||||||||||
систему является гармоническим с амплитудой B' и частотой ω': |
|||||||||||||||
|
|
|
|
|
|
|
|
′ |
|
′ |
|
|
|
|
(2.215) |
|
|
|
|
|
f (t) = B |
sin ω t . |
|
|
|
|
|||||
Тогда, согласно формуле (2.211), получимИ |
|
|
|
||||||||||||
|
|
|
|
K (q(A) |
+ jq |
′ |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
′ |
′ |
|
|
|
|||||
|
1 |
− |
|
|
(A)) |
|
− jϕ′ |
. |
(2.216) |
||||||
|
(T1 |
+T2 )ω′− jω′(1 |
|
ДА |
= B e |
||||||||||
|
|
|
−T1T2 |
ω′2 ) |
|
|
|
|
|||||||
|
|
|
|
|
|
А |
|
|
|
|
|
|
|||
Для заданной частоты ω' |
|
|
|
|
|
|
|||||||||
= 10 рад/с и заданных параметров ли- |
|||||||||||||||
нейной части системы K = 10; T1 = 0,01; T2 = 0,02 и реле Bр = 10; bр = 4 |
|||||||||||||||
кривая Z(А) |
|
|
|
|
б |
|
|
|
|
|
|
|
|
|
|
изображена на р с. 2.40, а. Проведя окружности разных |
|||||||||||||||
радиусов В, |
по точкам пересечения определим зависимости А'(B') и |
||||||||||||||
|
|
|
и |
|
|
|
|
|
|
|
|
|
|
||
φ'(B') (рис. 2.40, б) для вынужденных колебаний при частоте ω'. |
|||||||||||||||
jV(ω) |
|
|
а |
|
б |
25 |
|
|
|
А'(B') |
150° |
||||
|
С |
|
|
|
|
|
φ'(B') |
|
|||||||
|
|
|
|
|
|
А' |
|
|
|
|
|
φ' |
|||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
15 |
|
25 |
35 |
|
15 |
|
|
|
|
|
|
90° |
|||
|
φ' |
|
|
U(ω) |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
10 |
|
|
|
|
|
|
60° |
||||
A |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
Z(A) |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
Акр5 |
|
|
|
|
|
|
30° |
||||
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
0 Bкр 5 |
10 |
15 |
20 |
25 |
B' |
|||
|
Рис. 2.40. Определение области захватывания |
|
|
96