Материал: 6251

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

ω ≠ 1 , и эта даже небольшая разница с течением времени приводит к всё

бо́льшим отличиям получаемого решения от реального процесса. Ляпунов и независимо от него Рэйли предложили несколько видоиз-

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

пока значению

ω≠1. Поскольку при = 0

должно выполняться равен-

ство ω =1 или ω2 =1, частота ω должна зависеть от µ

 

 

 

 

 

 

 

 

 

 

 

ω2 = ω2 (µ)=1+ µh + µ2h + ...

,

 

 

 

 

(2.3.37)

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

 

 

где h1, h2 и т.д. – частотные поправки.

 

 

 

 

 

 

 

 

 

 

 

 

 

Перепишем уравнение (2.3.28) в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ɺxɺ+ ω2x = (ω2 1)x + µϕ(x, xɺ).

 

 

 

 

 

(2.3.38)

Подставив ряды (2.3.31) и (2.3.37) в уравнение (2.3.38), получим

ɺɺ

ɺɺ

 

2

ɺɺ

+ ...)+ ω

2

(x0

+ µx1 + µ

2

x2

+ ...)= (µh1

+ µ

2

h2 + ...)×

 

(x0 + µx1 + µ

 

x2

 

 

 

 

× (x

+ µx

+ µ2x

2

+ ...)+ +µϕ(x , xɺ

0

)+ µ2

∂ϕ

 

 

x

+ µ2

∂ϕ

 

xɺ

+ ... .

 

 

 

 

 

0

1

 

 

 

 

 

 

0

 

 

 

 

x

 

1

 

 

xɺ

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x=x0

 

 

 

x=x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x=x0

 

 

 

 

x=x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ɺ

ɺ

 

 

 

 

ɺ ɺ

 

Далее процедура известна – раскрываем скобки и приравниваем множители при µ в одинаковых степенях в правой и левой частях полученного выражения. В результате получаем бесконечную систему линейных дифференциальных уравнений

116

ɺxɺ

+ ω2 x

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ɺxɺ

+ ω2x = h x + ϕ(x

, xɺ

),

 

 

 

 

 

 

 

 

 

 

1

1

1

0

0

0

 

 

 

 

 

 

 

 

 

 

 

ɺxɺ

+ ω2 x

 

= h x

+ h x

+

∂ϕ

 

 

 

x + ∂ϕ

 

 

xɺ

,

(2.3.39)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

1 1

2 0

 

 

 

 

x= x

0

1

ɺ

 

x=x

1

 

 

 

 

 

 

 

 

 

x

 

ɺ ɺ

 

x

 

ɺ ɺ

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x= x0

 

 

 

x=x0

 

 

... .

Последовательно решая уравнения (2.3.39), получаем все больше и больше членов ряда (2.3.31). Частотные поправки hk вычисляются из условий равенства нулю вековых членов последовательно по всей цепочке уравнений. Попутно узнаём всё точнее неизвестную частоту ω в соответствии с формулой (2.3.37).

2 . 3 . 9 . Гар мо ническая линеар изация

Метод гармонической линеаризации предложили Н.М. Крылов и Н.Н. Боголюбов в 1934 г. В приложении к САУ этим методом занимались Л.С. Гольдфарб и Е.П. Попов. В учебной и технической литературе встречаются и другие названия этого метода – метод гармонического баланса, метод описывающих функций, метод эквивалентной линеаризации.

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

По-прежнему представляем систему состоящей из нелинейного звена

117

с характеристикой ϕ(x)и линейной части, заданной передаточной функ-

цией W(s) (рис. 2.27).

W(s)

y

φ(x) x

Рис. 2.27. Нелинейная система

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

x(t)= x0 + Asinωt ,

(2.3.40)

где x0 – постоянная составляющая, A – амплитуда, а ω – частота автоколебаний.

На самом деле автоколебания в нелинейной системе не являются моногармоническими вследствие искажения сигнала нелинейным звеном. Но предположение о моногармоническом сигнале на выходе линейной части (или на входе нелинейного звена, что одно и то же) основано на факте, что в большинстве случаев линейная часть является фильтром нижних частот (гипотеза фильтра). Это означает, что амплитуды высших гармоник много меньше, чем амплитуда основной частоты ω . Поскольку частота автоколебаний ω выясняется только в результате применения метода гармонической линеаризации и заранее неизвестна, проверку гипотезы фильтра возможно провести только в конце исследования. Если в

результате такой проверки гипотеза фильтра не подтверждается, к полу118

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

На выходе нелинейного звена сигнал уже не будет моногармоническим из-за нелинейных искажений (на рис. 2.28 показан пример прохождения синусоидального сигнала через усилительное звено с зоной нечувствительности).

Поскольку сигнал y на выходе нелинейного звена периодический, можно разложить его в ряд Фурье, а так как через линейную часть проходит только первая гармоника (гипотеза фильтра), то имеет смысл интересоваться именно этой гармоникой. В результате получим

y(t)= ϕ(x)= y0 + C1 sinωt + C2 cosωt ,

 

 

 

(2.3.41)

 

y

 

y

 

 

 

 

 

 

a

 

 

π/2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

π

 

 

 

 

α

x

2

 

π

ωt

 

 

 

 

 

 

 

 

π

А

 

 

 

 

 

 

2π

Рис. 2.28. Искажение сигнала нелинейным звеном

где постоянная составляющая y0 = 1 2πϕ(x)d(ωt), а коэффициенты 2π 0

Фурье С1 и С2 вычисляются по известным формулам

119

C1 = 1 2πϕ(x)sinωt d(ωt),

π 0

C2 = 1 2πϕ(x)cosωt d(ωt).

π 0

Выражение (2.3.41) удобнее представить несколько в другом виде. Найдём из формулы (2.3.40) sinωt и cosωt

sinωt =

x x0

, cosωt =

p(x x0 )

 

, ( p =

d

)

 

Aω

 

 

A

 

 

dt

и подставим в формулу (2.3.41). Получим

 

 

 

 

 

 

 

k

 

 

y = k

 

x

 

+ k

 

+

г

p

(x x )

 

 

 

 

 

г0

 

0

 

г

 

ω

 

0

или

 

 

k

 

 

y = y0

+ kг +

г

p

(x x0 ) .

ω

 

 

 

 

В приведённых соотношениях введены обозначения

 

 

 

 

 

y0

 

 

 

 

1

2π

 

kг0 =

=

 

 

 

ϕ(x)d(ωt),

 

 

 

 

 

x

 

2π

 

 

 

 

 

x

 

 

 

 

 

0

 

 

 

 

 

0

 

 

0

 

 

 

C

 

 

1

 

 

2π

 

 

kг

=

1

 

 

=

 

 

 

 

 

ϕ(x)sin(ωt)d(ωt),

 

 

 

 

 

 

 

 

 

A

 

 

Aπ 0

 

 

 

 

 

C2

 

 

 

1

 

 

2π

 

 

kг′ =

 

=

 

 

 

ϕ(x)cos(ωt)d(ωt).

 

 

 

 

 

 

 

 

A

 

 

 

Aπ 0

 

 

 

120

(2.3.42)

(2.3.43)

(2.3.44)