Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
2023.pdf
Скачиваний:
85
Добавлен:
07.01.2021
Размер:
2.62 Mб
Скачать

При гармонической линеаризации интервал значений 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)Bsin ω′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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jt+ϕ )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]