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

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

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

В соответствии с выражением (152), учитывая, что т = О,

запишем:

 

 

Х(а,

о>) = клд(а)-(Тк1

+ Тм1)и*;

Y (а,

со) = со + К ^

со - Гк 1 Гм 1 со3 .

Вычисляя, согласно неравенству (83), соответствующие про­ изводные и учитывая знаки производных коэффициентов гармо­ нической линеаризации (см. рис. 25, кривую / ) , получаем

( ё ) ' = - 2 ( ^ + ^ 1 ) Я < 0 ;

[да )

k A da ) \ < 0 при Л/6 > 1,3;

( £ ) • = * . ( £ ) • > « - p .

Из второго уравнения системы (152) при т = 0 имеем

Учитывая это, найдем значение производной

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

 

(т* + т<о(% ) * > r K l r M l a ( f ) * .

 

(156)

При амплитуде Л lb

> 1 , 3 ,

как это

следует из

приведенного

 

дХ \*

 

 

 

 

 

 

 

( -g^-j , неравенство

(156)

удовлетворяется,

и перио­

дическое

решение устойчиво.

 

 

 

 

 

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

изменении амплитуды

1 =^ Л lb <

1,3

необходимо

определить

значения

производных

(-^-^ н

(~da)

'

Определим эти

значения

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

\da ) \ Да ) ' \da ) \ Да / '

150

Из графиков рис. 25 (см. кривую I) находим, что при 1 ^

< Alb < 1,3

Подставляя эти значения в выражение (156), имеем

тп

+ ты >

1.66Q.

 

Реализуемыми являются решения уравнений (155)

только

при х ^ 0. На основании

этого

из второго уравнения

системы

(155) находим максимальное значение частоты периодического решения

С учетом

полученного

значения Qmax

условие

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

периодического решения

при

изменении

1 «g Alb <

1,3

будет

 

7 ^ 1 + 7 \ ц > 1,66 ] / 7 ^ 7 V

 

 

(157)

Для основного варианта значений параметров

 

исследуемой

САР загрузки

Тк1 == 0,5 с, Тм1

= 0,063 с критерий

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

выполняется.

Следовательно,

вся кривая

AAhjb (kn,

с,

b) ^ 1

соответствует устойчивому периодическому решению (автоколе­ баниям).

Для сравнения на рис. 56 приведены зависимости A^Jb X

X

(k„, с, b)

и Q (kn,

с, b) для релейной характеристики с зоной

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

и релейной

характеристики

общего вида.

Из

рисунка

видно,

что наличие

зоны возврата

m < 1 (неодно­

значность характеристики) существенно сужает область устой­ чивости.

Так при изменении m от 1 до 0,5 критическое значение пара­ метра х к р уменьшается с 27,9 до 6,35, т. е. в 4,4 раза. При этом амплитуда автоколебаний увеличивается, а частота уменьшается. При m <С 1 частота автоколебаний зависит от изменения пара­ метров kn, с, Ь, т. е. Q — Q (k„, с, Ь), а при m = 1 Q = const. Очевидно, что при т >> 0 область устойчивости еще более сужается.

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

151

Исследование САР с учетом люфта или зазора в регуляторе

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

Выясним влияние нелинейности типа люфта или зазора в ре­ гуляторе на устойчивость и автоколебания САР, описываемой системой уравнений (127). Так как влияние запаздывания т на

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

выяснено,

то

для упрощения расчетов

примем

т = 0,

с

учетом

сделанных

раньше

замечаний относительно

нового

смысла

постоянной времени

Тм1.

 

 

Нелинейность типа люфта или зазора без ограничения пред­

ставлена на рис. 28, а; это неоднозначная

нелинейная

характе­

ристика с

гистерезисной

петлей.

 

 

 

Гармоническая линеаризация нелинейной характеристики типа

люфта

или зазора

дает

 

 

 

 

 

 

 

 

F(Ah3)=\q(A)+^-

Q

р

Ah,

 

 

 

 

 

 

 

Г\

"'"3.

 

где q (А) и q' (А) — коэффициенты, описываемые уравнениями (104).

Характеристическое уравнение гармонически линеаризован­ ной системы имеет вид

TKlTulffi

+ (TKl + 7-м 1 2 + [1 + кл ЗЖ]

р + клЯ(А) = 0. (158)

Подставив в это уравнение р = jQ, получим уравнения для определения амплитуды и частоты периодического решения:

kaQ(A)-(Ta

+ TuJQ*

=

0;

)

 

а + клд,(А)-Тк1Тла*

=

0.

J

( [ Щ

Уравнение (159) не может быть удовлетворено при значениях амплитуды Alb < 1; значит, периодическое решение возможно только при Alb 5г 1, так как при уменьшении амплитуды ко­ лебаний ведущего звена до значений, меньших половины люфта, прекратится передача движения на ведомое звено.

Выделим параметр нелинейности k (см. рис. 28, а), для чего

воспользуемся

преобразованиями q (Л) = kq

[Alb) и q'

(Л) =

= kq' (Alb).

kQ = k„k, где

k0 — коэффициент

 

 

Обозначим

усиления

линей­

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

системы.

 

 

Уравнения (159) позволяют определить влияние основных параметров системы на границу области устойчивости и пара­ метры автоколебаний. Эти уравнения трансцендентны относи­ тельно Л и Я, но разрешимы относительно k0, Тк1 и Ти1.

152

Поэтому, найдя параметр k0

из

первого

и

второго

уравне­

ний (159) и подставив выражения

q {Alb)

и

q' {Alb),

получим

я ( Т к 1

+ Г м 1 )

 

 

 

 

я ,

.

(.

2

 

 

 

 

Alb

J

У

A/b\

A/b

)

(160)

я

 

(Q-TK1TU1QS)

 

 

 

 

J -

 

i

L _ \

 

 

 

 

Alb

\

 

Alb )

 

 

 

 

Методика определения зависимости амплитуды и частоты периодического решения от параметров системы, а также нахож­ дения критического значения коэффициента k0 к р графо-анали- тическим путем аналогична подробно изложенной при исследова­ нии релейной характеристики общего вида.

Задаваясь одинаковыми значениями Q и изменяя значения Alb

в уравнениях

(160), получим два семейства

кривых

k0 = kQ {A/b)

при Тк1

=

const

и Тм1

= const. Точки пересечения этих

кривых

с одинаковыми значениями Q будут графическим решением урав­

нений

Alb

= f {k0)

и

Q =

Q (k0).

 

 

 

 

 

 

Результаты графо-аналитического решения уравнений (160)

при Тк1

=

0,5 с и Тм1

=

0,063 с следующие:

 

 

 

й , с " 1

. . .

1,6

1,8

2,3

2,35

2,4

2,5

3

4

 

Alb

 

 

1,2

1,35

1,95

2,0

2,1

2,2

2,9

5,25

 

k0

 

 

9,3

7,0

6,1

6,05

6,2

6,45

7,65

10,4

 

По этим данным

построен график рис. 57, из которого видно,

что для данной системы в зависимости от значения

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

усиления k0 имеется три области состояния. При малом

коэффи­

циенте

усиления

0 <

k0

<

k0l к р

имеется

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

равно­

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

Аналитически

коэффициент

k0l

к р

можно

вычислить,

если из уравнений (160),

исключив

Q,

получить

выражение

k0

= k0 {А) и

приравнять

производную - j | нулю. Для рассматриваемой САР это сложно,

и коэффициент k0l. к р находим

приближенно в результате графо­

аналитического

решения.

 

 

 

kob

к р ^

Средним значениям общего коэффициента усиления

kQ

s-S ko2

к р

соответствует область периодических

решений,

в которой

имеются две ветви для амплитуды и частоты.

 

 

За

областью периодических

решений

при k0

> ko2, к р

нахо­

дится

область

неустойчивости.

Значение

ko2_ к р

равняется

зна-

153

чению коэффициента усиления krp на границе устойчивости ли­ нейной модели рассматриваемой системы без учета люфта в регу­ ляторе:

k

— h

Tig -f- Тш

 

"•02. кр

"-гр

 

Действительно, для линейной системы q (А) — k, и из пер­

вого

уравнения системы

(159)

находим гр

(TKi + Т м 1 )

со

 

 

 

 

 

г р ,

а из

второго уравнения

согр

1

 

 

 

 

 

VTK1TU

 

 

лг/fcmj

3,268

2,451

1,634

0,817

Рис. 57. Границы областей устойчивости и неустойчивости и изменение параметров автоколебаний астатической САР с нелинейностью типа люфта или зазора в плоскости общего коэффициента усиления системы

Равенство ko2 к р = %тр можно вывести исходя из следующих соображений. Числитель второго уравнения системы (160) об­

ращается

в нуль

при Q2 = -=—~—, а

амплитуда

Alb

при

конеч-

ном kn

неограниченно

 

'

К1 У Ml

 

 

 

А

=

оо и

Q2

=

возрастает. Подстановка

1

 

в первое уравнение системы (160) и дает

значение

k,02. кр-

 

 

Следовательно,

 

при

k,02.

кр А

=

О О ,

^max ~

 

y-f

 

 

По-

этому

при kQ > ko2,

к р

имеются

расходящиеся

процессы,

т.

е.

находится

область

неустойчивости.

 

 

 

 

 

 

 

Из

равенства

ko2

к р

=

krp

следует,

что граница

неустойчи­

вости

нелинейной

системы

совпадает

с

границей

колебательной

устойчивости линейной системы. Значит, область периодических решений занимает часть области устойчивости линейной системы без учета люфта.

Из рассмотренного вытекает, что исследование данной САР загрузки как линейной без учета люфта или зазора в регуляторе далеко не отражает действительной картины процессов и дает 154

ошибочное представление об области устойчивости, так как внутри области устойчивости при наличии люфта или зазора возможны установившиеся периодические режимы (автоколебания) с прак­

тически неприемлемыми

амплитудами.

 

 

 

Из рис. 57 видно, что периодическое решение AAhJb

(k0) имеет

две

ветви.

Исследуем

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

периодического

решения

с помощью

аналитического

критерия

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

 

В

соответствии с уравнением

(158)

получаем

 

 

 

X (а, со) = knq

(а) -

к1

+

Тя1) со2;

 

 

 

Y (а, со) =

со +

К

со -

Гк 1 Тм 1 со3 .

 

Вычисляя соответствующие производные и учитывая знаки производных коэффициентов гармонической линеаризации, оче­ видные из графиков рис. 28, б, запишем:

 

( - £ ) * =

- 2

^ 1

+ 7 м 1 ) й < 0 ;

/

дУ \*

 

 

 

 

 

 

\

да> )

 

 

 

 

 

 

(dX_\*-k

(JLY\<0

при

1 <

Alb < 2;

Уда)

к

Л

а а

) { >

0 п р и А

/ Ь >

2 ;

Из второго уравнения (159) имеем

1 + * л ^ - =

ТаТж1№.

С учетом этого найдем значение производной

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

 

(TKl

+ TMl)(^y>TKlTMlQ(-^y.

 

 

(161)

Из

графиков

q (А)

и <?' (А) (см. рис. 28,

б) видно,

что для

ветвей

этих коэффициентов с большими амплитудами

при Alb —*

оо q(A) и q' (А) одинаково стремятся к постоянным

значениям

и, следовательно,

их

производные по а для

больших

амплитуд

155

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

Т Т:

- > " ш а х -

Ml

М.'

 

М

 

Так

как

Q m a x ==

1

, то условие устойчивости

периоди-

ческого

решения

 

У

ТмТм

 

 

 

 

 

будет

 

 

 

 

 

 

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

A\hjb

(k0)

при Alb

> 2

критерий

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

удовлетворяется.

Для

нижней

ветви

найденного

периодического

решения

при

Alb <С 2 критерий

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

не выполняется,

значит нижняя

ветвь

 

принадлежит

неустойчивому периодическому

решению,

т. е. является границей устойчивости системы «в малом».

 

 

Перенося результаты, полученные на границе области перио­

дических

режимов

на

области

О <С k <t k0i.

кр

и

& >

Къ. К р >

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

§

16. УСТОЙЧИВОСТЬ И АВТОКОЛЕБАНИЯ

С

УЧЕТОМ ОГРАНИЧЕНИЯ СКОРОСТИ КОМБАЙНА

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

(TKlp+ \)Av = F(AH).

Если в качестве начальных условий при исследовании принять значения параметров системы и постоянного внешнего возмуще­ ния в установившемся режиме с расчетной подачей, обеспечивае­ мой при среднем значении скорости на данной передаче путем настройки регулятора, то можно считать, что нелинейность F (АН) имеет симметричную, однозначную характеристику насы­ щения (см. рис. 26, г).

Для рассматриваемой системы нелинейные звенья 6 и 1 регу­ лятора (см. рис. 16, а) оказались соединенными последовательно. ~ Вместе с тем по условиям метода гармонической линеаризации выходной сигнал нелинейного, звена должен фильтроваться, проходя инерционную линейную часть, чтобы выполнялись для системы условия обобщенного свойства фильтра.

156

Допустимость для рассматриваемой САР (с последовательно

соединенными релейным усилителем и нелинейным

регулирующим

органом) метода гармонической линеаризации проще

показать

при записи уравнения звена 6 в такой форме: ТврАН

= F (Ah3),

где Т 6 — постоянная времени звена, отнесенная

к единичному

перемещению поршня

гидроцилиндра исполнительного меха­

низма;

Тв - i - . При

этом для релейной характеристики коэф­

фициент

гармонической

линеаризации

Выделяя в нелинейных звеньях 6 и / линейные части и нели­ нейности (безынерционные элементы нелинейных звеньев), полу­ чаем:

А6 = F (Ah3); ТврАН = А9; Дг|> = F (АН);

к1р + 1) Av = Дф.

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

применение

метода

гармонической линеаризации.

Гармоническая линеаризация звеньев 6 и 1 соответственно

дает:

 

 

 

 

 

 

F(Ah3)

=

q6(Ae)Ah3,

 

 

F(AH)

=

<7i г) АН,

где с76

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

гармонической линеаризации релей­

ной характеристики

с зоной нечувствительности; qx г) — коэф­

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

Характеристическое уравнение гармонически

линеаризован­

ной системы имеет

вид

 

 

 

 

TKiTMlPs

+

к1

+

Тк1) р 2 +

р +

W e e ^ ' W i

= 0. • (162)

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

метров системы

вводим

преобразования:

 

 

9* =

-

f <?е(Ав/Ь);

gi =

fhffi(Aifbi).

 

157

Подставляя в уравнение (162)

р = /Q и е-'та

= cos xQ —

— / sin TQ, С учетом сделанных

преобразований

получаем два

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

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

и

параметры

 

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

 

 

Х =

й л С ^ 1 с о 8

т О - 6 ( 7 ' к 1 + 7'и 1 )£ 2 я =0;

1

Y

=

bQ — kacqeqi

sin xQ — bTKlTMlQ3

= 0.

J

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

TO-TQ =

1 ~ T K I T M I Q 2

1 б Т " -

0 ( Г к 1 + Г м 1 ) '

что аналогично ранее полученному выражению (132). Это понятно, так как было установлено, что параметры нелинейности не влияют на частоту автоколебаний. Поэтому все выводы, сделанные при

анализе

САР

с

релейным

регулятором

в отношении зависимости

Q (Тк1,

Ты1,

х),

остаются

в силе и для

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

с двумя

нелинейностями.

 

 

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

баний

находим из системы

уравнений

(163):

 

 

 

 

 

 

 

 

 

"

 

( Г

- - +

Г м 1 )

Q 2

 

 

 

(164)

 

 

 

 

 

<76?IC0STQ

 

 

 

 

 

 

Критические

значения

настроечных

параметров

kn,

с и

Ъ

получим,

подставляя в

это

выражение

(<76<7i)max- Коэффициенты

(<7e<7i)max

можно

найти

обычным

путем

отыскания

экстремума

функции. Для этого необходимо иметь

зависимость

амплитуд

Ав

и А1

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

входе звеньев

6 и 1.

 

звена

6

 

Из

уравнения

гармонически

линеаризованного

 

 

 

 

AH

=

 

-^qe(A6)Ah3

 

 

 

 

 

 

находим

выражение передаточной функции этого звена

 

 

 

 

 

We(s)

=

 

-j-q,{At).

 

 

 

 

 

 

Учитывая, что

Ah3 =

Л 6

sin

Ш,

а

АН

=

A^in

(Qt

- f Р), по­

лучаем искомую зависимость амплитуд Ав

и

Аг:

 

 

 

 

 

 

 

A1

=

 

\Wa(s)\s=iaA,

 

 

 

 

 

 

158

или, подставляя сюда

значение модуля передаточной

функции

W6

(s), окончательно

имеем

 

 

 

 

 

А = -q- <7в (Л) Л-

(165)

' q6qx

Учитывая выражение (165)

и обозначая произведение

функций

через

Я, запишем

 

 

 

 

 

4 A )

=

<7e(A)fc(A).

 

 

Теперь,

задаваясь

из

условия допустимой установившейся

ошибки регулирования различным значением параметра Ь, на­

ходим Ав

из равенства

= 0, при котором

lqe в) qx

6 ) ] т а х .

Однако

ЭТОТ путь

является трудоемким.

Найдем

(<76<h)max

иным путем. Очевидно, что максимальное значение произведения коэффициентов qe(AJb) и q1(A1/b1), являющихся для рассма­ триваемой системы функциями одного аргумента, будет равно произведению максимальных значений этих коэффициентов (см.

рис. 25, кривую /7 и рис.

27, кривую III):

[q,{Ajb)qx

(AjbJ]

 

 

max —

Я 6 max<?l max-

Из рис. 27 (см. кривую / / / ) видно,

что <7l m a x = 1 при А11

sg; 1, что соответствует линейной зоне нелинейности с насыщением.

Следовательно, при Axlbx^s~

1, т.

е. по существу для линейного

звена 1 (см. рис. 16, a)

(q9qi)max

= Яцтах и Д л я определения

критических значений настроечных параметров получается урав­ нение, совпадающее с уравнением для системы, имеющей только релейный усилитель. Это естественно, поскольку в обоих слу­ чаях звено 1 рассматривается как линейное с qx х) = kv

Определим, при каких значениях параметров Ьх и сг нелиней­ ности насыщения действителен полученный вывод. Предвари­ тельно уточним некоторые условия. Для исследуемого реального объекта — комбайна статический коэффициент усиления регули­ рующего органа ^ п р и работе на определенной передаче остается постоянным. Поэтому при анализе нелинейной характеристики регулирующего органа (звено /) можно принять kx = const и исследовать влияние на устойчивость и автоколебания других параметров нелинейности: ширины линейной зоны Ьх и величины насыщения сг (см. рис. 26, г), соответствующего максимальной поступательной скорости комбайна. При kx = const эти пара­ метры связаны однозначно cxlbx = const.

Для определения влияния ограничения регулирующего воз­ действия — поступательной скорости комбайна на устойчивость и автоколебания при расчетах по выражению (164) с учетом выра­ жений для коэффициентов qe и qx удобнее пользоваться пара­ метром blt от которого легко перейти к параметру сх = kxbx.

159

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