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

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

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

где q {Alb) и q' {Alb) — безразмерные коэффициенты гармони­ ческой линеаризации нелинейности.

 

q(Ajb)

= я

А/Ь

\ V'

(A/bf

 

V

 

(A/mbf

 

 

 

(93)

 

 

 

 

 

(при

Alb s> 1);

 

 

 

 

 

 

 

 

 

g'(A/b)

=

- 2

l

(1 — m)

 

(при A/b-^

1).

 

 

 

 

 

 

 

 

л

(Л/6)2

 

 

 

 

 

 

 

 

 

 

 

Зависимости

q {Alb)

и

q' {Alb),

рассчитанные

по

формулам

(93)

при

I m I =

0,5,

приведены

на рис.

25 (кривые / ) .

 

 

 

 

 

 

 

 

 

 

 

 

Релейные характеристики,

 

 

 

 

 

 

 

 

 

представленные

на рис. 24, в

 

 

 

 

 

 

 

 

 

и г можно рассматривать как

 

 

 

 

 

 

 

 

 

частные случаи

релейной

ха­

 

 

 

 

 

 

 

 

 

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

 

общего

вида

 

 

 

 

 

 

 

 

 

при т — 1 (рис. 24, в)

и т =

 

 

 

 

 

 

 

 

 

=

1,

Ь =

0

(рис.

24,

г).

 

 

 

 

 

 

 

 

 

 

 

Релейная

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

 

 

 

 

 

 

 

 

 

с

зоной

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

 

 

 

 

 

 

 

 

 

(рис. 24,

в).

Значения

коэф­

 

 

 

 

 

 

 

 

 

фициентов гармонической ли­

 

 

 

 

 

 

 

 

 

неаризации для этой нелиней­

 

 

 

 

 

 

 

 

 

ности

получаем,

подставляя

 

 

 

 

 

 

 

 

 

в

выражения

(91)

и

(92)

Рис.

25.

Безразмерные

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

т = 1:

 

 

 

 

 

 

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

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

 

релейной

 

 

 

1-

 

Ь*_

 

 

характеристики общего вида (I), с зоной

 

 

 

 

 

 

 

 

пА

 

 

 

(94)

нечувствительности (II) и

идеальной (III)

 

 

 

 

Л 2

 

 

 

 

 

 

 

 

 

 

q

{А)

0

(при

А:

Ь).

 

 

 

Из выражения (94) находим коэффициент гармонической линеа­

ризации

в безразмерном

виде

 

 

 

 

 

 

 

 

 

 

 

 

q{A/b)

_4

1 _

 

 

 

 

 

(при A/b^

1).

 

(95)

 

 

я

А/Ь

 

 

(А/Ь)*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По формуле (95) находим, что

q{A/b)max

=

при A lb

=

У 2.

Зависимость

q {Alb)

приведена на рис. 25 (кривая

/ / ) .

 

 

 

Идеальная релейная характеристика (рис. 24, г). Для этой ре­

лейной характеристики имеем т =

1, b =

0, тогда в соответствии

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

(91)

и (92)

получаем

 

 

 

 

 

 

 

 

 

 

 

9 ( Л ) = - £ р ,

 

д'{А)

=

0

( п р и Л > 0 ) .

 

,

_(96)

80

Откуда коэффициент гармонической линеаризации в безраз­ мерном виде будет

q(A/b)=^~

(при А/Ь>0).

(97)

Зависимость q [Alb) показана на рис. 25 (кривая / / / ) .

Нелинейные характеристики с насыщением

Характеристика с насыщением и зоной нечувствительности

(рис. 26, а). Так как нелинейная характеристика однозначна, коэффициент q' (А) = 0. Прохождение сигнала через нелинейность возможно только при условии, что А >• Ь.

/(As/пф)

 

 

 

г

 

 

 

 

//'

 

//

 

4>,Ь *-Цг\\

[7

4> = 9t

 

b

х

 

 

6)

 

 

 

 

S)

 

F

 

Рис.

26.

Определение

<

 

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

гармо­

А

с

нической

линеариза­

0

 

ции нелинейной харак­

Ь,

х

теристики

с

насыще­

нием

и зоной

нечувт

 

 

ствительности (а), с зо­

г)

 

ной

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

 

 

сти без насыщения

(в);

 

 

с насыщением без зоны

 

 

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

(г)

Коэффициент гармонической линеаризации q (А) для рассма­ триваемой нелинейности находим при х = A sin Ш (А > bj. При колебаниях входной величины по гармоническому закону, выход­ ная величина нелинейности представляет периодическую функ­ цию F (A sin i|)), из графика которой (рис. 26, б) видно, что зна-

чение интеграла J F (A sin 1|з) sin -ф d\p в выражении коэффи-

о

циента q (А) будет одинаковым для каждой четверти периода. Учитывая, что в интервале 0 г|з =s£ F(A sin ij;) = 0 , по­

лучим

6 Н . Н. Н а с т е н к о

81

Разбиваем интервал

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

на два участка ^ 1 sg; - ф г ^

< Ф г И 1 | ) 2 < 1 ) ) < у И

определяем для

каждого из них значение

функции F (A sin г|)). Для первого интервала из графика рис. 26, а

находим, что F (х) = k (х Ь) или,

учитывая, что х — A sin tp,

имеем F (A sin tp)^= k (A sin tp — b).

Для второго интервала

F (A sin ty) = с.

 

Подставляя в подынтегральное выражение коэффициента q (А)

значения F (A sin tp),

с = k (bx

— b), ipx =

arcsin

и ap2 =

= arcsin - ~ , окончательно находим

 

 

q(A) =

I arcsin -j

arcsm-j

+

 

q' (A) = 0 (при A 5s bx).

Для приведения коэффициента q (А) к безразмерной форме воспользуемся допущением, что bx = nb, где п = 2, 3, 4, . . ., т. е., что ширина линейной зоны кратна ширине зоны нечувстви­ тельности, тогда

q {Alb,

n) = kq{A/b,

п),

 

где q {Alb, п) = ±

( arcsin ^

- arcsin

+

 

 

 

 

 

 

(99)

 

Зависимость

# (Л/6, «) показана на рис. 27 (кривая / ) .

 

Нелинейные

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

представленные

на

рис. 26, в

и

г, можно рассматривать как частные

случаи

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

с

насыщением и зоной нечувствительности.

 

 

 

Характеристика с зоной нечувствительности без насыщения

(рис. 26, б). Эту характеристику

можно

рассматривать

как част­

ный случай характеристики рис. 26, а при условии Л ^

Ьх. Тогда

 

Л

 

 

 

 

 

яр2 = -2"i и в соответствии с уравнениями (98)

 

 

 

q (А) = k - - f - ( arcsin \

+ ±

У 1 -

;

(100)

 

 

 

 

 

 

<7'(Л) = 0 (при Л^гЬ )

82

Из выражения

(100) находим

 

 

 

 

q {Ajb) =

kq(Alb),

 

где q {Alb) =

1

^- ( arcsin

A/b ^ Alb - V

{A/b)2

 

 

(при Alb =s 1).

(101)

Характеристика с насыщением без зоны нечувствительности

(рис. 26, г). Для звена с такой характеристикой при b = 0 и ipr = 0, пользуясь уравнениями (98), получим

9 (Л) =

\ arcsin' 1л

1 л |/

Л2

(102)

 

<?'(Л) = 0

(при

А^ЬХ).

 

При Л < £>х коэффициент с7 (Л) = /г, т. е. характеристика ли неиная.

Рис. 27. Безразмерные коэффициенты гармонической линеаризации нелинейной характеристики насыщения с зоной нечувствитель­ ности (1), зоны нечувствительности (II) и насыщения (III )

При амплитудах колебаний входной величины, захватывающих зону насыщения, данное звено можно считать линейным с тем меньшим коэффициентом усиления q (Л), чем больше ампли­ туда (в противоположность предыдущему случаю).

Из уравнений (102) имеем q{A/b1)=^kq{Alb1),

x

) = А

(a r c s l n

+ ^_

1

где q'{Alb

j /

 

 

(при

A/bi^l).

(103)

6*

83

Из выражений (102) и (103) следует, что

<?(Л)ш а х = К a ~q (A/bJmn = 1 (при Albx = 1).

Нелинейная характеристика типа люфта или зазора

Характеристика этого типа представлена на рис. 28, а. Приведем выражения коэффициентов q (А) и q (А) в оконча­

тельном виде [22]:

 

, № = i [ i + a r c s i „ ( l - i f )

+

 

 

 

 

+ Ф - т ) } / Т Т т г т ) ]

<"р«^^

<104)

Вводя

безразмерные

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

q

и

q ,

получим

 

 

 

 

q (А)

=

kq {Alb),

 

q'(A)

=

kq~' {Alb).

 

 

 

Выражения

для

определения

q (A/b)

и

q' (Alb)

получаем

из

формулы

(104).

Зависимости

q(Alb)

и

q'

(Alb) приведены

на

рис. 28,

б.

 

 

 

 

 

 

 

 

 

 

 

 

Нелинейная характеристика с

насыщением

 

 

 

 

 

при несимметричных

колебаниях

 

 

 

 

 

 

 

 

Гармоническая линеаризация

нелинейности

с

насыщением

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

- T / l - < V ^ £ ) +

(», +

х», arcsin A + i L

_

 

— (&х — х°) arcsin h ~ A

X ° ]

(при A bx +

| хй

|).

(105)

9 (Д

x°) =

arcsin

 

(- arcsin - i - j

 

^

 

 

 

 

 

*B)

 

 

 

 

 

 

 

Л2

 

 

 

A

+ *

] / 1 _ i ^ + £ ! ) i

(при Л ^ 6 1 +

-0|).

(106)

84

Рис. 28. Нелинейная харак­ теристика типа люфта или зазора (а) и безразмерные коэффициенты гармонической линеаризации этой характе­

ристики (б)

Рис. 29. Коэффициенты гармонической линеаризации нели­ нейной характеристики насыщения при несимметричных колебаниях

На рис. 29, а, б приведены соответственно зависимости П Р И А/Ьх — const и q(Alb1) при = const [22]. Из

рис. 29, а видно, что при наличии колебаний входной величины нелинейного звена его статическая характеристика для медленно меняющегося воздействия (функция смещения) сглаживается. Из рис. 29, б видно, что увеличение смещения центра колебаний при­ водит к уменьшению коэффициента усиления q для колебательной составляющей.

§ 12. МЕТОД МОДЕЛИРОВАНИЯ^ НА АНАЛОГОВЫХ ВЫЧИСЛИТЕЛЬНЫХ МАШИНАХ

Метод математического моделирования на электронных цифро­ вых (ЦВМ) и аналоговых (АВМ) вычислительных машинах широко применяется при исследованиях динамики систем автоматического регулирования и управления.

Простота программирования и набора задачи, высокое быстро­ действие обусловили использование аналоговых вычислительных машин для решения задач как анализа, так и синтеза автомати­ ческих систем [14]. Большую роль при этом играет возможность визуального наблюдения (на экране электроннолучевого инди­ катора) за протеканием исследуемых процессов и удобство варьи­ рования значениями параметров системы. Точность решения на аналоговых машинах удовлетворяет в основном требованиям к точ­ ности инженерных расчетов САР мобильных сельскохозяйствен­ ных агрегатов. При более жестких требованиях к точности рас­ четов исследования на АВМ целесообразны для предварительной оценки влияния параметров системы на качество процесса регу­ лирования.

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

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

86

Методика подготовки и решения задач анализа динамических качеств систем автоматического регулирования и управления на аналоговых машинах кратко состоит в следующем.

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

Критерием достоверности модели и точности решения может служить совпадение результатов моделирования на аналоговой машине и реального процесса при натурных испытаниях.

На стадии предварительных расчетов точность результатов решения моделированием оценивают сравнением их с результа­ тами решения этой же задачи аналитическим методом. При совпа­ дении результатов математическая модель считается достоверной.

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

Постановка задачи

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

Исследование динамических качеств САР загрузки методом моделирования включает в себя решение следующих основных вопросов:

1) устойчивости, автоколебаний и качества переходных про­ цессов в одноконтурных астатической и статической (нелинейный

илинейный варианты) САР загрузки с запаздыванием;

2)влияния на устойчивость, автоколебания и качество пере­ ходных процессов одноконтурной астатической САР загрузки кор­ ректирующих устройств, обеспечивающих различные законы ре­ гулирования;

87

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

4)точности следующих структурных вариантов САР загрузки при случайных внешних возмущениях:

а) одноконтурной системы; б) одноконтурной системы с корректирующими устройствами,

обеспечивающими реализацию различных законов регулирования; в) системы с регулированием по возмущению; г) комбинированной системы регулирования; 5) точности решения на модели.

Математические модели исследуемых систем

Одноконтурная астатическая САР. Система описывается урав­ нениями (56); структурная схема ее приведена на рис. 16, а. В линейном варианте этой системы, при исследовании моделиро­ ванием, нелинейная функция F (A/i? ) заменена линейным выра­ жением k& A/i 3 . При моделировании принято Т\ = 0 и ТА\ = О, так как влияние этих постоянных времени на динамические ка­ чества системы, как показали расчеты, незначительное.

Для расчета исходного варианта блок-схемы математической модели исследуемой САР записываем значения коэффициентов уравнений, относительно которых будут варьироваться параметры системы при моделировании:

г 2 = о , 1 с 2 ; г к 1 = о,5с; т1 = 8,9- Ю" 3 с2 ; тм1 =

6 , 3 - к г 2 с; .

£/ =

0,286 с; kt

= 0,0148м/(с-мм); kt = 3,6 кгс-м-с/кг;

k2

= 2 кг/м;

kH = 0,075 м2 /с; k5 = 0,34 мм/(кгс-м).

 

Коэффициенты k2

и kK определены при J5p = 5 м, Q =

40 ц/га

и vcp = 1,5 м/с. Принятый диапазон изменений

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

усиления

линейной

части нелинейной системы kn = 0,01-^-0,06

(где kn =

k1k2kik5)

охватывает возможные изменения

каждого

из коэффициентов усиления звеньев системы. Диапазон изменения транспортного запаздывания т = 0-М,25 с. Постоянные времени изменяются на ± 2 5 % относительно приведенных значений. Ско­ рость исполнительного механизма ограничивается наименьшим значением его постоянной времени, которую, исходя из эксплуа­ тационных требований, принимают равной примерно 2 с. В реаль­ ном объекте максимальное изменение скорости комбайна огра­ ничено и равно при работе на второй передаче примерно 1 м/с (см. рис. 5).

Одноконтурная статическая САР. Система описывается урав­ нениями (57); структурная схема ее приведена на рис. 16, б. В линейном варианте этой системы нелинейная функция F (Ahy)

88

заменена линейным выражением ke>

A / i y , т. е. не учитывается

ограничение изменения координаты

Д # .

Значения коэффициентов уравнений первого, второго и третьего звеньев объекта аналогичны значениям, приведенным для аста­

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

 

Коэффициенты

уравнений

остальных

звеньев

имеют следующие значения:

 

 

 

 

 

 

1,9. ю - 2 с2 ;

Т д 1 =

7,1

( Г 2 с;

 

 

 

Т\ = 0,148 с2 ;

T„i =

0,661 с;

 

 

 

£4 = 7 мм - с/кгс.

 

 

Коэффициент

усиления

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

нелинейной

системы

(/гл = klk2kik5)

варьируется

в пределах 0,07—0,42 изменением на

модели коэффициента усиления &4 . Диапазон изменения транс­ портного запаздывания т = Он-0,8 с. Постоянные времени из­

меняются на

± 2 5 % относительно приведенных

значений.

 

—,йЛлд

&9т

AM

 

AQ

 

 

 

 

 

 

 

 

 

 

АН

М

 

 

 

A9z

AM

 

 

 

 

 

 

 

\Ах

 

 

 

 

 

 

 

а)

 

Al

 

{вариант

 

 

 

H Вдриант_

 

AQ

 

 

I вариант

/[вариант

AQ

 

 

&v\Ag

&9г

AM

 

 

 

 

 

Ail

 

 

4fc

AM

 

Ah,

 

 

 

 

 

 

Ah

 

 

 

 

Ак

 

 

 

 

 

 

 

 

В)

 

 

 

 

 

 

 

Рис. 30. Функциональные схемы CAP:

 

 

 

о д н о к о н т у р н о й с

к о р р е к т и р у ю щ е й

о б р а т н о й с в я з ь ю (а);

о д н о к о н т у р н о й с

введением

п р о ­

и з в о д н о й р е г у л и р у е м о й величины

(б); к о м б и н и р о в а н н о й (в); р а з о м к н у т о й

с р е г у л и р о в а ­

 

нием

по в о з м у щ е н и ю

(г)

 

 

 

При исследовании на модели рассматриваем симметричную характеристику нелинейности с насыщением, что соответствует настройке регулятора, при которой оптимальная подача обеспе­ чивается при среднем значении управляющего воздействия равном 40 мм (половине хода поршня цилиндра вариатора ходовой части).

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

Различные законы регулирования осуществляются путем вы­ бора уравнения обратной связи, охватывающей интегрирующий

89

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