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

vostrikov

.pdf
Скачиваний:
182
Добавлен:
23.01.2019
Размер:
12.99 Mб
Скачать

10.2. Метод гармонического баланса

341

которое, если принять обозначения

 

q ( A,

)

 

b1

,

 

 

 

 

 

 

 

 

1

 

 

 

A

 

 

 

 

 

 

(10.10)

 

 

 

 

 

c1

 

 

 

q ( A,

)

 

 

,

 

 

 

 

 

 

 

 

2

 

 

 

A

 

 

 

 

 

 

 

можно записать в виде

 

 

 

 

 

 

 

 

 

u

q ( A, )

 

q2 ( A, )

.

(10.11)

 

 

 

1

 

 

 

 

 

 

 

 

Здесь q1(A, ) и q2 (A,

) – коэффициенты гармонической линеари-

зации.

 

 

 

 

 

 

 

 

 

Как видим, уравнение нелинейного звена (10.11) с точностью до высших гармоник является квазилинейным. При постоянных значениях амплитуды входного сигнала A коэффициенты гармонической линеаризации q1(A, ) и q2 (A, ) являются постоянными. Однако раз-

личным значениям амплитуды A соответствуют разные коэффициенты q1(A, ) и q2 (A, ) . В этом заключается отличие гармонической

линеаризации от обычной (см. разд. 8).

Таким образом, вместо нелинейного элемента с характеристикой (10.5) можно рассматривать эквивалентное линейное звено, поведение которого описывается уравнением (10.11). Оно может быть представлено в операторной форме

u q1( A, )

q2

( A, )

p .

(10.12)

 

 

Для гармонически линеаризованного нелинейного элемента можно записать передаточную функцию

WНЭ ( p, A, )

u

q1

( A, )

q2

( A, )

p

(10.13)

 

 

 

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

WНЭ (A, j ) q1(A, ) jq2 (A, ) .

(10.14)

342 Глава 10. АНАЛИЗ ПРОЦЕССОВ В НЕЛИНЕЙНЫХ СИСТЕМАХ

В случае статической нелинейной характеристики вместо (10.5) имеем

uf ( )

иуравнение (10.12) принимает вид [6]

u q1( A)

q2

( A)

p

,

(10.15)

 

 

где коэффициенты гармонической линеаризации

q1(A) и q2 (A) зави-

сят только от амплитуды. При этом получим передаточную функцию

WНЭ ( p, A, )

q1( A)

q2

( A)

p

(10.16)

 

 

 

и частотную характеристику

 

 

 

 

 

 

WНЭ (A, j )

q1(A)

jq2 (A)

 

(10.17)

статического нелинейного звена.

Для однозначной статической нелинейной характеристики коэффициент q2 (A) 0 , и вместо (10.14) получим

WНЭ (A, j ) q1(A) .

(10.18)

Коэффициенты гармонической линеаризации типовых статических нелинейных звеньев приводятся в литературе (например, в [6, 40]).

ПРИМЕР 10.2

Определить эквивалентную передаточную функцию нелинейного звена, которое представляет собой идеальное реле (рис. 10.5).

 

 

 

 

 

 

 

Поскольку идеальное реле имеет одно-

 

 

u

 

 

 

 

 

 

 

 

 

значную статическую

характеристику, вы-

 

 

c

 

 

 

 

 

 

 

ражение для его передаточной функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(10.16) имеет вид

 

 

 

 

 

–c

 

 

 

 

 

WНЭ ( A, j )

q1( A) ,

 

 

 

 

 

где коэффициент q1(A)

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

Рис. 10.5. Статическая харак-

 

теристика идеального реле

 

 

10.2. Метод гармонического баланса

 

 

343

q ( A)

1 2

f ( A, t) sin( t) d ( t)

4c

.

 

 

 

 

 

 

1

A 0

 

A

 

 

 

 

Далее, учитывая полученные выражения для передаточных функций гармонически линеаризованных нелинейных элементов (10.14), (10.17), рассмотрим соотношения метода гармонического баланса.

10.2.3. ОСНОВНОЕ УРАВНЕНИЕ МЕТОДА ГАРМОНИЧЕСКОГО БАЛАНСА

Если в системе, изображенной на рис. 10.3, гармонически линеаризовать нелинейный элемент, заменив его эквивалентной передаточной функцией WНЭ ( p, A, ) , то она становится линейной (рис. 10.6). Сле-

довательно, в этом случае для анализа свойств системы можно применять методы линейной теории управления.

v = 0

 

 

 

 

 

 

u

 

y

 

 

W

( A,

, p

)

 

 

Wл(p))

 

 

 

 

 

 

 

 

 

 

 

 

НЭ( A,

 

 

 

 

 

 

НЭ

 

 

 

 

л

 

Рис. 10.6. Структурная схема гармонически линеаризованной системы

Как известно, в линейной системе (при отсутствии синусоидального сигнала на входе) незатухающие колебания будут возникать лишь в том случае, когда она находится на границе устойчивости. Таким образом, для определения автоколебаний в исходной системе (см. рис. 10.3) необходимо рассмотреть условие границы устойчивости линеаризованной системы. В соответствии с критерием Найквиста в этой ситуации амплитудно-фазовая характеристика разомкнутой системы должна

проходить через точку 1,

j 0 , т. е.

 

 

 

 

Wр ( A, j )

1 .

 

 

Учитывая, что Wр ( A, j )

WНЭ ( A, j

)Wл ( j

) ,

запишем условие гра-

ницы устойчивости в виде

 

 

 

 

WНЭ (A, j )Wл ( j )

1.

(10.19)

344

Глава 10. АНАЛИЗ ПРОЦЕССОВ В НЕЛИНЕЙНЫХ СИСТЕМАХ

Это уравнение и представляет собой основное уравнение метода гармонического баланса, из которого можно определить параметры автоколебаний. Если (10.19) не имеет положительных вещественных решений относительно A и , то автоколебательный режим в нелинейной системе не возникает.

Для решения основного уравнения метода гармонического баланса были предложены различные способы, которые мы далее последовательно и рассмотрим.

10.2.4. АНАЛИТИЧЕСКИЙ СПОСОБ ОПРЕДЕЛЕНИЯ АВТОКОЛЕБАНИЙ

В этом случае для определения границы устойчивости линеаризованной системы удобнее воспользоваться критерием Михайлова. С этой целью решим характеристическое уравнение системы (см. рис. 10.6)

F( p, A, ) 1 WНЭ ( p, A, )Wл ( p) 0,

(10.20)

заменив в котором p на j , получим условие границы устойчивости согласно критерию Михайлова

F(A, j ) 1 WНЭ (A, j )Wл ( j ) 0 .

(10.21)

Заметим, что (10.21) есть преобразованная форма записи основного уравнения метода гармонического баланса (10.19).

Выделяя вещественную и мнимую части F( A, j ) , получим соотношения

Re F ( A, j ) 0,

(10.22)

Im F ( A, j ) 0,

которые позволяют аналитически вычислить параметры автоколебаний

A и .

ПРИМЕР 10.3

Определить параметры автоколебаний в системе (см. рис. 10.3), если нелинейный элемент представляет собой идеальное реле (см. рис. 10.5) с уровнем ограничения c , а линейную часть описывает передаточная

2

функция Wл ( p) p( p2 p 1) .

10.2. Метод гармонического баланса

 

 

 

 

 

 

 

 

 

345

 

Запишем передаточную функцию гармонически линеаризованного иде-

 

 

ального реле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W

( A, j

)

 

 

4c

4

,

 

 

 

 

 

 

 

 

 

НЭ

 

 

 

 

A

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а затем характеристическое уравнение (10.18)

 

 

1

 

 

4

 

 

 

2

 

 

 

 

 

0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A p( p2

 

p

1)

 

 

 

 

 

 

 

 

 

 

 

которое преобразуем к виду

 

 

 

 

 

 

 

 

 

 

 

A( p3

p2

 

p) 8 0 .

 

Заменив здесь p на j , получим

 

 

 

 

 

 

 

 

 

 

A( j

3

2

 

j

)

 

 

8

0 .

 

В результате расчетные соотношения имеют вид

 

 

 

 

 

 

 

A 2

 

8

0,

 

 

 

 

 

 

 

 

 

 

3

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, в нелинейной системе будут возникать периодические

 

движения со следующими параметрами:

 

1 и A 8 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10.2.5. ВЛИЯНИЕ ПАРАМЕТРОВ СИСТЕМЫ НА ПЕРИОДИЧЕСКИЕ ПРОЦЕССЫ

В некоторых случаях возникает необходимость оценить влияние одного из параметров системы (обозначим его ) на автоколебания. При этом уравнение (10.21) принимает вид

F(A, j , ) 1 WНЭ (A, j )Wл ( j , ) 0 .

Следовательно, соотношения (10.22), кроме параметров A и , содержат :

Re F ( A, , ) 0,

(10.23)

Im F ( A, , ) 0.

346

Глава 10. АНАЛИЗ ПРОЦЕССОВ В НЕЛИНЕЙНЫХ СИСТЕМАХ

Решив уравнения (10.23) относительно A и , получим параметри-

ческие зависимости

 

 

 

 

A A(

),

(10.24)

 

(

)

 

 

и построим соответствующие графики (рис. 10.7).

A

A0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

гр

0

 

 

 

 

 

 

 

 

 

 

 

гр

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 10.7. Пример влияния параметра α на периодические процессы

По этим графикам можно выбрать значение

 

 

0 , при котором в сис-

теме будут возникать периодические процессы с определенной амплитудой и частотой ( A0 и 0 ).

10.2.6. СПОСОБ ГОЛЬДФАРБА

Решение основного уравнения метода гармонического баланса (10.19) относительно амплитуды и частоты автоколебаний можно получить графически.

В способе Гольдфарба [6, 19] предлагается решить основное уравнение относительно частотной характеристики линейной части системы следующим образом:

Wл ( j )

1

 

.

(10.25)

 

 

WНЭ ( A, j

)

 

 

 

Затем на комплексной плоскости строятся амплитудно-фазовая характеристика Wл ( j ) и характеристика, соответствующая нелинейному элементу, т. е.

10.2. Метод гармонического баланса

347

1

.

(10.26)

 

 

WНЭ ( A, j )

Если эти две характеристики не пересекаются, то периодических процессов в нелинейной системе не возникает.

При наличии пересечений частота автоколебаний определяется по частотной характеристике линейной части системы Wл ( j ) , а амплитуда – по характеристике нелинейного элемента в точке пересечения.

Поскольку в общем случае точек пересечения Wл ( j ) и характери-

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

Устойчивость найденного колебательного режима позволяет оценить следующее правило (оно не является строго обоснованным, но зачастую оказывается полезным). Если при движении по обратной частотной характеристике нелинейного элемента в сторону увеличения амплитуды происходит пересечение амплитудно-фазовой характеристики линейной части «изнутри наружу», то этой точке пересечения соответствуют устойчивые колебания (автоколебания). В противном случае колебания будут неустойчивыми.

На рис. 10.8 характеристики W ( j ) и

W (A, j ) 1

пересекают-

л

НЭ

 

ся в двух точках. Это означает, что в системе могут возникать два вида колебаний.

 

ImIm

11

1

 

НЭ( A)

1

WНЭ ( А)

 

 

eRe

 

Wл ( jj ))

 

2

 

2

Рис. 10.8. Иллюстрация способа Гольдфарба

348 Глава 10. АНАЛИЗ ПРОЦЕССОВ В НЕЛИНЕЙНЫХ СИСТЕМАХ

Причем первой точке пересечения соответствуют устойчивые колебания (автоколебания) с амплитудой A1 и частотой 1 , а второй точке – неустойчивые.

ПРИМЕР 10.4

Определить параметры колебаний и проверить их устойчивость для системы, изображенной на рис. 10.6. Здесь нелинейный элемент представляет собой идеальное реле (см. рис. 10.5) с уровнем ограничения c , а передаточная функция линейной части следующая:

 

 

Wл ( p)

 

 

2

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p( p2

 

 

 

 

 

 

 

 

 

 

 

 

p 1)

 

 

 

 

Получим выражение для

амплитудно-частотной характеристики

(рис. 10.9) в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wл ( j

)

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j (

2

j

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

Wл ( j )

 

 

2

 

2

 

j2

(1

2 )

.

4

2 (1

2 )2

 

4

 

2 (1

2 )2

 

 

 

 

Запишем выражение для частотной характеристики нелинейного элемента

ImIm

-

 

1 1

 

W ( A)

 

WНЭ

( A)

 

 

 

 

 

 

 

 

 

-2

 

0 0

 

Re

 

–2

 

 

 

Re

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W ( j ))

л

Рис. 10.9. Расположение характеристик (10.23) для примера 10.4

4с 4

WНЭ ( A, ) q1( A) A A ,

а затем построим годограф (рис. 10.9)

1

 

A

.

WНЭ ( A)

4

 

Как видим, эти характеристики пересекаются в одной точке, которая соответствует автоколебаниям. Для определения их параметров найдем координаты точки пересечения, для чего приравняем нулю мнимую часть Wл ( j ) :

10.2. Метод гармонического баланса

349

 

ImWл ( j )

 

2

(1

2 )

 

0 .

 

4

 

2 (1

2 )2

 

 

 

 

 

Отсюда следует, что

a

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При найденном значении частоты получим

 

 

 

 

 

ReWл ( j

 

a )

2 .

 

 

Из условия

 

 

 

 

 

 

 

 

 

 

ReWл ( j a )

2

 

A

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

определим амплитуду автоколебаний:

Aa

8 .

 

 

10.2.7. СПОСОБ КОЧЕНБУРГЕРА

Способ Коченбургера представляет собой второй вариант графического решения основного уравнения метода гармонического баланса (10.19). В этом случае его предлагается решить относительно характеристики нелинейного элемента системы следующим образом:

WНЭ ( A, j )

1

 

.

(10.27)

 

 

 

Wл ( j

)

 

 

 

Как и в способе Гольдфарба, точки пересечения двух характеристик согласно (10.27) свидетельствуют о наличии в системе колебательного режима. Причем частота колебаний определяется по обратной частотной характеристике линейной части

системы Wл 1( A, j ) , а амплитуда – по

характеристике нелинейного элемента WНЭ ( j ) в точке пересечения (рис. 10.10).

Процедура определения автоколебаний аналогична способу Гольдфарба, однако правило формулируется следующим образом. Если при движении по характерис-

ImIm

-

 

1 1

 

W ( jw )

 

Wл л ( jw)

WWНЭ((A, jwjw))

НЭ

ReRe

Рис. 10.10. Определение автоколебаний способом Коченбургера для примера 10.5

350

Глава 10. АНАЛИЗ ПРОЦЕССОВ В НЕЛИНЕЙНЫХ СИСТЕМАХ

тике нелинейного элемента в сторону увеличения амплитуды происходит пересечение обратной частотной характеристики линейной части «снаружи внутрь», то этой точке пересечения соответствуют автоколебания. В противном случае колебания будут неустойчивыми.

ПРИМЕР 10.5

Определить параметры колебаний методом Коченбургера и проверить их устойчивость для системы, показанной на рис. 10.6. Здесь нелинейный элемент представляет собой идеальное реле (см. рис. 10.5) с уровнем ограничения c , а передаточная функция линейной части следующая:

Wл ( p)

12

 

.

(5 p 1)(3 p2

 

 

2 p 1)

Определим необходимую характеристику реле (см. пример 10.4)

4

WНЭ ( A, ) WНЭ ( A) A

и построим ее на плоскости (рис. 10.11). Запишем выражение для обратной частотной характеристики линейной части

1

13

2

1

j

15

3

7

 

 

 

 

 

 

 

 

 

 

Wл ( j )

 

 

12

 

 

 

 

12

 

 

 

 

 

 

 

 

 

и найдем несколько ее значений (табл. 10.2).

Im

1

( j )

 

–1

 

 

 

W

 

-

 

л- Wл ( j,ω)

 

W ( A, j )

WНЭн э (Aj)

Re

Рис. 10.11. Иллюстрация способа Коченбургера

Соседние файлы в предмете Основы Теории Управления