Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
ТАУ лекции 2.pdf
Скачиваний:
573
Добавлен:
11.05.2015
Размер:
1.97 Mб
Скачать

6. ВРЕМЕННЫЕ ХАРАКТЕРИСТИКИ НЕЛИНЕЙНЫХ СИСТЕМ

Рассмотренные до сих пор переходные процессы в нелинейных системах часто носили качественный характер, поскольку время в явном виде не фигурировало. Для количественной оценки качества переходного процесса необходимо построить зависимость x = f (t) .

6.1. Построение переходного процесса методом фазовой плоскости

Пусть фазовая траектория системы имеет вид, приведенный на рис. 6.1, а. Возьмем начальную точку M 0 (x0 , x&0 ) . Разобьем ось x на участки (жела-

тельно равные) x1 - x0 , x2 - x1, x3 - x2 и т.д. Если интервалы xi - xi-1 малы, то можно считать, что на интервале скорость изменения координаты неизменна и определяется средней величиной

x&

=

x&i + x&i-1

=

Dxi

=

xi

- xi-1

.

(6.1)

 

 

 

 

 

 

 

 

ср.i

 

 

2

 

 

Dti

Dti

 

 

 

 

 

 

 

Отсюда

 

 

 

Dxi

 

 

 

xi - xi-1

 

 

 

 

Dt =

 

= 2

,

(6.2)

 

 

 

 

 

i

 

x&ср.i

 

 

x&i + x&i -1

 

 

 

 

 

 

 

 

 

 

 

ti +1 = ti + Dti.

 

 

 

(6.3)

Естественно, отрезки Dti в общем случае будут различными по величине. Точность построения зависит от величины интервалов разбиения оси на участки. Кривая переходного процесса x(t) в данном случае представлена на

рис. 6.1, б.

x&

x&8

x6

x7

x3 x2 x1 x0

x

 

x5

x4

 

 

x8

 

x&0

M 0

x&5

x&1 x&4

 

x&3

 

x&2

 

x

x0 g

x1 g

x2 g

x3 g

x4

0

 

 

g

 

 

 

 

x8g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x5

 

1

2

3

4

5

6

7 t

 

g

gx

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

x6 g

Рис. 6.1. Построение переходного процесса по фазовой траектории: а – фазовая траектория, б – кривая x = f (t)

117

6.2. Построение переходного процесса методом гармонической линеаризации

Рассмотрим колебательное звено второго порядка. В линейном звене свободные движения (при v(t ) = 0 ) соответствуют характеристическому уравне-

нию T 2l2 + 2zT l +1 = 0 , корни которого l = g ± jw . Эти движения определя-

ются выражением

 

 

 

 

1,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x = A egt sin (wt + j),

 

(6.4)

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

1 - x2

 

 

 

причем g =

 

= const, w =

 

 

 

= const , так как

показатель

колебательности

T

 

 

x = const .

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

В нелинейном звене g = g(t ) и w = w(t ), поэтому

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

òg(t )dt

æ t

ö

 

 

 

 

x = A0e0

sin ç òw(t)dt + j÷.

(6.5)

 

 

 

 

 

 

 

è 0

ø

 

В связи с этим коэффициенты гармонической линеаризации в общем случае будут функциями A, w, g. Выражение (4.15) с учетом этого примет вид:

ìR e [D (jw, A, g, П )] = 0,

(6.6)

í

îIm [D (jw, A, g, П )] = 0.

 

По этим выражениям строят диаграммы качества нелинейных колебаний. Это семейства линий g = const и w = const на плоскости (A, П), где П – параметр, например, коэффициент усиления системы.

Для этого построения из первого уравнения (6.6) выделяют w = f1(A, g, П) и подставляют во второе, получают А = f2 (П, g) . Придавая g постоянные значения, строят семейство, показанное на рис. 6.2, а. Затем, подставляя А = f2 (П, g) в w = f1 (A, g, П) , исключают g и получают А = f3 (П, w) - рис. 6.2, б.

Отметим, что для линейных САУ линии на рис. 6.2 вертикальны.

Из рис. 6.2, а следует, что при g = 0 можно найти амплитуду автоколебаний (для точки C – амплитуда AC ). Слева от заштрихованной линии– затухающие колебания, справа – расходящиеся.

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

Если параметр П соответствует точке L (затухающие колебания), то переходный процесс будет от начальной точки M0 с амплитудой A0 до A ® 0.

Если же П соответствует точке E, то характер процесса зависит от выбора начальных условий. Если это точка M02, то колебания расходятся до автоколе-

118

баний. Если это точка M03, то колебания затухают до автоколебаний. В обоих

случаях амплитуда автоколебаний будет A4.

 

 

 

 

 

 

 

 

При этом частота изменяется в соответствии с графикомA = f3 (П,w)

по

вертикали. На рис. 6.3, б линии w = 0 соответствуют монотонные переходные

процессы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

A

 

g < 0 g = 0

б

A

w = 0

 

w0

 

 

 

 

 

 

 

A

M

g > 0

 

 

 

 

 

 

 

 

 

01 g0 M 03

 

 

 

 

w

 

 

 

 

 

 

A1

0

 

г

 

 

 

 

A

 

 

 

1w2

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

w3

 

 

 

 

 

 

 

 

 

 

 

C

A2

 

 

 

 

 

 

 

 

AC A

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

A3

 

 

 

 

 

 

 

 

 

 

A3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

g4

 

 

 

 

 

A4

 

 

 

w4

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

E

M

02

П

 

П = const

П

 

 

 

 

 

 

П = const

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 6.2

 

 

 

 

 

 

 

 

На основании этих рассуждений можно приближенно оценить переход-

ной процесс, определив при П = const огибающую и частоту колебаний.

 

 

 

Из

графика

на

рис. 6.2, а определяют T (A ) = -

1

,

а из рис.6.2, б –

 

 

 

 

 

 

 

 

 

 

 

 

i

i

gi

(Ai )

 

 

 

 

 

2p .

 

 

 

 

 

 

 

 

 

 

 

t

(A ) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

i

gi (Ai

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Огибающую можно построить, используя свойство экспоненты, заклю-

чающееся

в

том, что

 

длина

подкасательной

равна

постоянной

времени

(рис. 6.3).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A2

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

142443

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T0 (A0 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 42 43

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T1 (A1 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 6.3

 

 

 

 

 

 

T (A ) = - 1

; T (A ) = - 1

; T (A ) = - 1 ; и т. д.

 

 

 

 

 

 

0

0

 

g0

1

1

 

g1

 

2

2

g2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

119

 

 

 

Аналогично оценивается и период колебаний:

t

(A ) =

2p

;

t (A ) =

2p

;

 

 

 

 

 

 

 

 

 

 

 

0

0

w0

 

1

1

w1

 

 

 

 

 

2p

 

 

 

 

 

 

 

t

 

(A ) =

и т.д.

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

w2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончательно зависимость A = f (t) имеет вид рис. 6.4.

A

A @ A4 , w @ w3

0

2p/w0

4p/w1

6p/w2

t

8p/w3 10p/w4 12p/w5

T0 (A0 )

T1(A1)

T2 (A2 )

T3 (A3 )

Рис. 6.4

6.3. Компьютерное моделирование переходного процесса

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

Рассмотрим систему, структура которой приведена в примере5.1 на рис. 5.3, а ( M C учитывать не будем). Выполним исследование системы при подаче на вход гармонического воздействия. Будем считать, что усилитель K1 имеет характеристику типа «насыщение».

Схема моделирования системы в SIMULINK изображена на рис. 6.5. Нелинейный элемент представлен звеном, имеющим коэффициент

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

При подаче на вход системы сигнала v(t ) = 0,18sin 0, 25t осциллограммы выходного сигнала изображены на рис. 6.6.

120

 

1

 

 

 

v

 

 

 

НЭ

+

x

 

a

 

 

+

 

-a

 

 

 

a

 

 

K1

= 23000

-a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.052s+1

 

1

 

.0047

 

 

 

 

s

 

 

 

.12s+1

 

 

 

.003s+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 6.5. Схема моделирования нелинейной системы в SIMULINK

Рис. 6.6, а соответствует значению "a" =10 , при котором НЭ не входит в

насыщение, при этом ошибка x

= 4 ×10-4 . Рис. 6.6, б соответствует значению

max

 

 

"a" = 8 , при котором начинает сказываться

насыщение. Выходной сигнал от-

личен от синусоиды, ошибка x

= 7,5 ×10-3 ; рис. 6.6, в соответствует "a" = 6 –

max

 

 

НЭ входит в глубокое насыщение.

 

 

а

б

в

Рис. 6.6

При подаче на вход системы сигнала1(t) при значении "a" = 25000 , нелинейный элемент не входит в насыщение, переходной процесс и ошибка соответствуют линейной системе рис. 6.7, а. При a =10000 , нелинейный элемент входит в насыщение; переходной процесс затягивается рис. 6.7, б. При дальнейшем уменьшении "a" НЭ входит в глубокое насыщение, система становится «вялой», t p возрастает, перерегулирование уменьшается рис. 6.7, в.

а

б

в

Рис. 6.7

121