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

книги из ГПНТБ / Настенко Н.Н. Системы автоматического регулирования зерноуборочных комбайнов

.pdf
Скачиваний:
20
Добавлен:
24.10.2023
Размер:
10.1 Mб
Скачать

Поэтому исследование начнем с определения зависимости ча­ стоты Q автоколебаний от параметров системы.

Из уравнений (131) находим, что

 

 

 

 

( 1 3 2 )

Из этого выражения следует, что частота Q автоколебаний

системы зависит от постоянных времени объекта

Тмг

и регуля­

тора Тк1, а также от

постоянного

запаздывания

т,

т. е. Q =

= G ( Г к 1 , Тм1, т).

 

относительно

Q и разрешимо

Уравнение (132) трансцендентно

лишь приближенными

способами.

 

 

 

При значениях параметра т, имеющих место в САР загрузки, принимать для упрощения расчетов tg xQ ^ xQ нельзя. Так, при г = 0,75 с для САР 6-го порядка получено (см. § 13) Q 1 с - 1 , откуда tgr Q = 0,932 и tg xQ > T Q .

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

 

Преобразуя уравнение

(132),

запишем

 

 

 

 

 

 

 

 

т

 

I — TKLQ

tg тЯ

 

 

 

 

/1 оол

 

Задаваясь

различными

значениями

Q и

вычисляя

 

каждый

раз

Г м 1 , получаем

зависимость

Q = Q (Тн1)

при

Тк1

=

const.

Аналогично

можно

определить

Q = Q (Тк1)

при

Т м 3

=

const.

 

При расчетах по выражению

(133) задаваться следует

такими

значениями Q, которые

приводят

к положительным значениям

ТЫ1 и

Тк1.

 

Ты1

 

 

 

 

 

 

 

 

 

 

Для

обеспечения

>

0 необходимо, чтобы

 

 

 

 

 

1) 1 — 7 , K l Q t g x Q > 0 ,

откуда

t g r Q < ^ - ;

 

 

 

 

2)

T&Q -г tg TQ >

0, откуда tg TQ >

TKlQ.

 

 

f

Графическое определение интервалов

значений

частот

Q, вы­

текающее из указанных неравенств, показано на рис. 48, а. Иско­ мые интервалы частот Q соответствуют в каждом периоде изме­

нению

 

аргумента xQ функции

tg xQ, заключенной между пара­

болой

l/TKlQ

и линией — TKJQ

(на рис. 48, а эти участки выпол­

нены

сплошными

линиями).

 

 

Для

каждого периода изменения функции tg xQ с учетом най­

денных

интервалов частот Q,

обеспечивающих положительные

значения Тм1,

по

уравнению

(133) при Тк1

= 0,5 с и т =

1,25 с

определяем

зависимость

Q = Q (Тм1).

График этой зави­

симости

приведен

на рис. 48, б.

 

130

Из рассмотренного следует, что 1) при заданных параметрах системы имеется бесчисленное множество частот й автоколебаний, удовлетворяющих решению уравнения (133); 2) в пределах каж­ дого из интервалов частот Й, обеспечивающих условие положитель­ ности постоянных времени и удовлетворяющих решению уравне­ ния (133), с увеличением постоянных времени объекта ТЫ1 и регу­ лятора ТК1 частота й автоколебаний уменьшается.

Я,с-'

Рис. 48. Графо-аналитическое определение

зависимости

частоты

автоколебаний

от

постоянной

времени объекта и

постоянного

запаздывания

Если

ТМ1

~ 0,

что соответствует

снижению

порядка

системы

с 3-го до 2-го, то,

согласно уравнению

(132),

 

 

 

Из зависимости

(134)

следует,

что в системе 2-го порядка

без запаздывания автоколебания с частотой

й << оо

возникнуть

не могут.

Значит

при

исследовании

САР

загрузки

с т = О

учет постоянных времени объекта и регулятора имеет принци­ пиальное значение, как бы малы они не были. Пренебрежение одной из постоянных времени, пусть даже очень малой, приводит к уменьшению порядка системы, что вносит принципиальные ка­ чественные отличия (исключение автоколебаний при любых пара­ метрах системы) по сравнению с исследованием системы 3-го по­ рядка (в которой автоколебания возможны). При наличии запазды­ вания даже в системе 2-го порядка частота й ф 0.

При расчетах по уравнениям (133) и (134) находим, что для

системы 3-го

порядка

с параметрами ТМ1 0,063 с; ТК1

= 0,5 с

и т = 1,25 с

частота

й = 0,88 с - 1 (эти и последующие

значения

частоты взяты для первого интервала частот, нижняя линия на

рис.

48, б); для системы

2-го

порядка

с ТМ1

= 0;

ТК1

=

0,5 с

их—

1,25 с частота й ж 1,14

с - 1 ; для безынерционной

системы

с ТМ1

= 0; ТК1 = 0 и т =

1,25

с частота

Й

1,258

с - 1 .

Из

этих

9*

131

данных следует, что с уменьшением инерционности системы увеличивается частота автоколебаний.

Для основного варианта значений параметров системы 6-го по­ рядка было получено, что Q = 0,957 с - 1 . Эта частота превышает частоту автоколебаний системы 3-го порядка примерно на 8%. При оценке этого расхождения следует учитывать, что на изме­ нение частоты Q действуют два фактора: 1) снижение инерцион­ ности системы, в результате которого частота автоколебаний по­ вышается; 2) пренебрежение коэффициентом (7, вследствие чего частота автоколебаний снижается. Результирующее действие этих

факторов привело к сниже­ нию частоты, что свидетель­ ствует о более эффективном влиянии на нее изменений коэффициента U, чем суммар­ ного действия постоянных времени Т\, Т2и, Тли

Учитывая, что tg xQ яв­ ляется периодической функ­ цией аргумента тй, на осно­ вании зависимости (132) за­ пишем

Рис. 49. Зависимость

частоты автоколе­

tg (TQ -}- /от)

=

1

7'K1TM1Q2

 

 

 

 

баний от постоянного

запаздывания при

( Я = 0,

1,

2,

...). (135)

различной постоянной

времени и объекта

Из этого выражения следует, что одному и тому же значению запаздывания т (при Тк1 const и Ти1 = const) соответствует бесчисленное множество значений частоты у . Практическое зна­ чение, как было показано, имеет первый интервал (п = 0) измене­ ния Q при г = var. Для этого случая зависимость Q (т) на основа­ нии уравнения (132) представим в следующем виде:

x = „arctg

1 — 7'к 1 УМ 1 й2

(136)

 

Q (TKI + тМ1) •

 

На рис. 49 приведена построенная по уравнению (136) зави­ симость Q = Q м1, т). Из рисунка следует, что при заданных параметрах системы, с увеличением постоянного запаздывания частота автоколебаний уменьшается. При этом с уменьшением инерционности зависимость частоты Q от изменения -запаздыва­ ния т резко усиливается. Поэтому для малоинерционных САР загрузки влияние постоянного запаздывания на частоту й зна­ чительно сильнее, чем для существенно инерционных.

 

Для

основного варианта

значений

параметров системы Г м 1 =

=

0,063

с и Тк1 = 0,5 с изменение

запаздывания

до значений

т ^

1 с

вызывает резкое

изменение частоты

автоколебаний

132 "

(рис. 49); при дальнейшем увеличении т зависимость Q (х) стано­

вится

плавной.

 

 

 

 

 

 

 

При заданных параметрах системы максимальную

частоту Q m a x

определяем

при т = 0. Из

уравнения

(136)

находим,

что

 

 

 

" m ax

w -

f

 

 

 

 

 

 

 

V 1

K I '

M I

 

 

 

При т =

0 зависимость

Q ( Г к 1 ,

Г м 1 ) не

является периодиче­

ской. Следовательно, периодическое изменение частоты Q проис­

ходит из-за наличия в системе запаздывания т и появления

вслед­

ствие

этого

в выражении

(132)

периодической функции

tg xQ.

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

метров линейной части

Q к1, Тм1,

х).

 

 

 

Определим влияние основных параметров системы на ампли­

туду автоколебаний, т. е. установим

зависимость А (/гл, с,

Ь,

х,

Тк1, Ты1).

Общее выражение взаимозависимости указанных

пара­

метров,

найденное из

уравнений (131),

имеет вид .

 

 

 

Г ( T K I + Т ш ) Ц2

Ъ

< 1 3 7 .

 

4

cos тй

 

kjfi'

'

'

Для решения уравнения (137) могут быть использованы сле­ дующие способы. Если выражение q(Alb) разрешимо относительно амплитуды А, то из уравнения (137) можно определить искомую зависимость в общем виде. Если неразрешимо, то из уравнения (137) определяем q в функции любого параметра системы при постоянстве остальных, например, q = q (с), откуда с учетом выражения q {Alb) находим амплитуду А автоколебаний. Для взятого примера получим А = А (с). Аналогично может быть опре­ делена зависимость амплитуды А от изменения каждого параметра системы.

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

§ 15. ВЛИЯНИЕ ТИПОВЫХ НЕЛИНЕЙНОСТЕЙ РЕГУЛЯТОРА НА УСТОЙЧИВОСТЬ И АВТОКОЛЕБАНИЯ

Регулятор с нелинейностью релейного типа и зоной нечувствительности

Для САР загрузки с такой нелинейностью сравнительно просто обеспечивается изменение параметров kn и с. При исследовании подобных систем в ряде случаев удобно рассматривать указанные

133

параметры как один совокупный настроечный параметр х = кл~^-,

из значений которого определяются оптимальные значения каж­ дого из настроечных параметров.

Влияние параметров системы на границу области устойчи­ вости и параметры автоколебаний определяем по следующему уравнению, полученному из выражений (131):

1

к1т)1

(138)

х = п

cos тй

 

Согласно выражению (138) окончательно зависимость ампли­ туды А от параметров системы для рассматриваемого случая будет

 

 

^ _

26 ] / к (2х =

 

(139)

 

 

 

я

 

 

 

где

 

 

 

 

 

 

 

 

 

В =

cos TQ

 

 

 

 

 

 

При

этом

частоту

Q находим

из выражения (132).

Как

было

установлено,

Q =

Q к1, Тм1,

т), и имеется мно­

жество интервалов изменений частоты Q. Очевидно, что значения В

будут разными для различных интервалов частот Q, при этом для

обеспечения х >> О следует

брать cos xQ >

0.

Aig/Ь.

AS/(с мм)

6.532

4,902

3,268

' 1.634

2

 

141

Область 1—устойчивости равновесия

-||

|

1

1

1

|

1

1

1

|

1

Ь.58в\

i

//

'м!

i

Я,с-

 

 

 

 

Я

T„<

 

 

 

 

 

\

 

 

 

0,8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V'"

 

Рис. 50. Граница об­

 

 

'MI

0,6

ласти

устойчивости и

fatr

 

изменение

параметров

0,063с

 

 

'

1

 

 

 

автоколебаний релей­

 

Пйлпгтн

 

0,4

ной астатической САР

 

авпиколебаний

 

в плоскости

совокуп­

 

 

 

ного

настроечного па-

 

 

 

 

 

 

Li" 0,5c

 

 

раметра

kn

с

I

 

0,2

при

z

= 1,25c

 

различных

 

значениях

 

 

 

 

 

 

 

 

 

 

постоянной

времени

 

 

 

 

 

 

объекта

8

Ь

С учетом характера

изменения

коэффициента

q {Alb) (см.

рис. 25, кривая

/ / ) при условии

В =

const (т. е. при заданных

параметрах Тк1,

Тм1

и т) из выражения (139) следует, что с

изме­

нением параметра

х в пределах x m l n

-s^ х

оо амплитуда

Alb

изменяется в пределах оо ^

Alb ^

1. При этом экстремум

функ­

ции Alb (к) будет в точке

с

координатами

х = и г а Ш ,

Alb — ] / 2 .

134

 

 

 

 

 

 

 

 

 

Зависимость Alb (и)

при

В = const имеет

форму

кривых

(рис. 50), построенных

для

нижнего интервала

частот

Q (см.

рис. 48, б). Для остальных интервалов частот вследствие увеличе­

ния

частоты Q и, значит, В кривые

Alb (У) будут

располагаться

на

рис.

50 значительно

правее

и не

будут

иметь

практического

смысла.

Поэтому все

расчеты

проводим

только

для нижнего

интервала

частот. Как видно из рис. 50, зависимость Alb (%)

имеет две

ветви периодического решения. Установим, какая из

них соответствует устойчивому периодическому решению (авто­

колебаниям). Исследование проведем с помощью аналитического

критерия устойчивости периодического решения

(см. гл. I I I , § 10).

Для

этого, подставляя в характеристическое

уравнение (129)

р = /со,

получаем

 

 

 

X (а,

со) = knq cos тсо — (Г к 1 + Тм1)

со2;

 

Y (а,

со) = со — kaq sin тсо — TKlTma>3.

Найдем частные производные от X и Y по а и со:

(ш ) ' = ~ s i n т й ~ 2 + Т ^ Q ;

j£T = -^sinTQ(S)*;

u

/дХ\*/дУ\*

/дХ\*/дУ\*

п

Левая

часть выражения ( ж )

- ( ^ )

( _ j > 0 , o n p e -

деляющего устойчивость найденного периодического решения,

имеет

вид

 

 

 

 

 

 

k »

- 3 w s )

c o s t Q - K»QX

- 2

+ T * J Q S I N

T Q ] -

Согласно

уравнениям (130),

имеем

 

 

 

 

 

 

C O S T Q =

( ! k ± W ;

 

 

 

 

 

sinxQ=: '

" У "

^

Q.

 

С учетом

этого предыдущее выражение будет

 

 

*-(!)'

-(Tk1+J;i)Q2 о

+ г

к ^ 2 ) ]

. ( н о )

135

Выражение в квадратных скобках здесь отрицательное, а из рис. 25 (см. кривую / / ) видно, что

(агУ>0 п р и

KA/b<V%

( • g - ) * < 0 при

A/b>Т/2.

Отсюда вытекает, что критерий устойчивости периодического

решения удовлетворяется только при Alb >> ]/ 2 и, следовательно,

только верхние ветви кривых

(и) (см. рис. 50) соответствуют

автоколебаниям (устойчивому периодическому решению), чтс показано сходящимися стрелками. Нижние ветви кривых Alb (х)

соответствуют неустойчивому

периодическому

решению. Значит,

в рассматриваемой системе

автоколебания возможны при х >

>> x m l n и начальном отклонении амплитуды

А

1,416, т . е .

здесь имеется «жесткое» возбуждение автоколебаний. При началь­ ной амплитуде А •< 1,416 и х > x m l n в системе возникает неустой­ чивый периодический режим с переходом к устойчивости равно­ весия.

В связи с рассмотренным отметим, что при расчетах по урав­

нению (138) большая из двух амплитуд будет амплитудой

автоко­

лебаний.

 

При х < x m l n в системе невозможны периодические

решения

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

Влияние изменения параметров системы на границу области устойчивости определяем по следующим зависимостям:

при

В = В (ТкЪ Тм1, т) = const

 

 

 

и„п =

 

В:

 

 

"кр

О max

 

при

х = х

- | - j == const

 

 

При расчетах по уравнениям (138) и (139) получаем амплитуду автоколебаний входной координаты А1г3 нелинейного звена. Зная передаточные функции звеньев, можно определить амплитуду автоколебаний любой координаты системы. Перейдем к определе­ нию амплитуды автоколебаний регулируемой величины Ag.

Связь амплитуд автоколебаний этих координат находим через передаточные функции звеньев 4 и 5 Сем. рис. 16, а):

AAha = \Wi(s)\\W6(s)

| , = / а A&g.

(141)

136

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

стемы уравнений (127),

получаем

АД/ц

АAg'

и окончательно имеем

 

 

(142)

Из рис. 50 видно влияние инерционности объекта на границу области устойчивости равновесия в плоскости параметров &л , с и b и на параметры автоколебаний (на которые аналогично

 

 

 

 

 

кг/(с-ин)

 

 

1

 

 

1

/ г

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9,с-

 

 

 

 

 

 

 

 

 

1

 

Область /автоколебаний

 

 

 

 

 

 

 

6,532

i

 

 

 

 

 

1.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р и с

51.

Граница

об­

 

 

 

 

 

г-1.25 с

 

 

 

 

ласти устойчивости

и

4,902

 

1

1

/6

 

 

0,9

 

изменение

параметров

 

 

 

 

 

автоколебаний

релей­

 

 

• i

i n

 

 

 

 

 

 

ной

астатической САР

3,268

 

 

 

! /

Ч:=0,75с

 

0.6

 

в плоскости

совокуп­

 

l!

 

 

 

TKI-0,5c

 

 

 

и

 

с

па-

 

 

 

 

 

 

 

 

ного настроечного

 

 

 

 

 

 

 

 

 

 

 

 

раметра

кл

 

при

1,634 Ь 2

 

1 1 /

 

 

\

 

 

0.3

 

различных

значениях

141

 

l

У

 

 

Тм,-0.063с

 

 

 

 

 

 

 

 

 

 

постоянного

запазды­

 

 

 

 

 

 

\

 

 

 

 

 

 

 

вания

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

н

 

 

влияет

инерционность

регулятора

Тк1

при

Г м 1

= 0,5

с

и

т

=

=

1,25

с). Из рисунка следует также, что с повышением инерцион­

ности

сужается

область устойчивости

по параметрам

k„,

с,

b

и

увеличивается амплитуда А автоколебаний. Поэтому выбор пара­ метров по идеализированной математической модели САР, без учета некоторых постоянных времени, может привести в реальной системе к автоколебаниям. На рис. 50 амплитуда автоколебаний регулируемой величины рассчитана по выражению (142) для основ­ ного варианта значений параметров системы.

На рис. 51 показано влияние постоянного запаздывания на гра­ ницу области устойчивости в плоскости параметров kn, с, b и на параметры автоколебаний. Из рисунка следует, что с увеличением запаздывания т существенно сужается область устойчивости и повышается амплитуда автоколебаний.

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

137

1) определено влияние изменения основных параметров на границу области устойчивости равновесия, т. е. установлена

зависимость

 

(kn-^-^J^

= f(ТыЪ

%) при

Тк1 =

const;

 

 

 

 

 

 

 

 

=

0,5с

1

1 ь,- 0,5с

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

Область

 

 

 

 

 

 

 

 

 

^^автоколебаний

 

 

 

 

 

 

 

0класть

 

 

 

 

 

 

 

 

 

• устойчивосни

 

 

 

 

 

 

 

 

 

 

рабновеа/я

 

 

 

 

 

 

 

 

 

 

0,4

0,8

0,8

1,0

1.2

г,с

 

 

 

 

 

 

 

 

В)

 

 

 

Рис.

52. Зависимость критических значений

совокупного

настроечного

пара­

метра

( кл

)

релейной астатической

САР

от

инерционности

 

объекта

(а)

\о /к р

ипостоянного запаздывания (б)

2)определена зависимость амплитуды Л автоколебаний от

параметров системы: А = А (кл, с, Ь, т, Тм1) при Тк1 = const. Влияния изменения значений параметров системы на частоту Q

и амплитуду Л автоколебаний для исследуемой САР загрузки видны из графиков рис. 49—51.

Рис. 53. Зависимость критических значений отношения параметров

нелинейности релейного типа (Ь/с)Кр отпостоянного запаздывания при различ­ ных значениях коэффи­ циента усиления линейной

части системы

 

0,2

0,4

0,6

0,8

10

1,2

IX

 

 

Зависимости

критических

значений настроечных

параметров

от

инерционности объекта

Тм1, постоянного запаздывания т и

коэффициента

усиления линейной

части системы kn

приведены

на

рис. 52, а,

б и 53.

 

 

 

 

 

138

Полученные результаты позволяют определять для заданной САР значения настроечных параметров, обеспечивающие устой­ чивость системы или требуемые параметры автоколебаний.

Регулятор с нелинейностью насыщения и зоной нечувствительности

При работе этого регулятора возможны следующие три режима. Первый, когда амплитуда колебаний на входе нелинейности А ^ by, что соответствует работе регулятора с нелинейностью насы­ щения и зоной нечувствительности (см. рис. 26, а). При этом коэффициент гармонической линеаризации q (А, Ь, &,) описывается уравнением (98). Второй, когда Ь •< A < by, что соответствует работе с нелинейностью типа зоны нечувствительности (см. рис.

26, в).

В этом случае выражение для определения

коэффициента

q (А,

Ь)

имеет вид

уравнения

(100). И,

наконец,

третий, когда

А «=с b,

то сигнал

вообще не

проходит

через нелинейность.

Исследуем устойчивость и автоколебания САР, регулятор ко­ торой имеет нелинейность с насыщением и зоной нечувствитель­ ности.

Из выражения (138) можно записать общее уравнение для определения влияния параметров системы на границу области

устойчивости и параметры

автоколебаний

в виде

 

knk=

J ^

\ + T f Q 2

(143)

 

q (A,

b, by) cos rQ

 

 

Определить амплитуду

автоколебаний

из зависимости

(143)

в явном виде невозможно. Поэтому задачу нахождения зависи­ мости A (kn, с, b, by, т, Тк1, Тму) следует решать путем последова­ тельного определения зависимости амплитуды от изменения каж­ дого (или нескольких) параметров системы в отдельности, т. е. вторым из рассмотренных ранее способов.

Найдем влияние изменения настроечных параметров kji на амплитуду автоколебаний для основного варианта значений пара­

метров САР:

Тк1

= 0,5 с; Тм1

=

0,063

с;

т =

1,25 с, при

ширине

линейной зоны,

равной

ширине

зоны

нечувствительности, т. е.

при by — 2b.

В этом случае

выражение

(143)

запишем

так:

 

 

£ Л _

 

+ ^ M I )

Q2

 

(144)

 

 

л

b

q(A,

6)COSTQ

 

 

Полученная зависимость отличается от уравнения (138) для релейной характеристики с зоной нечувствительности только выражением для коэффициента q.

Зависимость Alb (kn, с, b), вычисленная по формуле (144), приведена на рис. 54. Эта зависимость имеет две ветви периоди­

ческого

решения [в соответствии с

двумя

ветвями изменения

q (Alb,

by) на рис. 27]. Определим,

какая

из них соответствует

139

Соседние файлы в папке книги из ГПНТБ