
книги из ГПНТБ / Настенко Н.Н. Системы автоматического регулирования зерноуборочных комбайнов
.pdfгде 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: |
|
|
|
|
|
|
|||||||
гармонической |
линеаризации |
|
релейной |
|
|
4с |
|
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, б) видно, что зна-
2л
чение интеграла 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, г) при несимметричных колебаниях дает следующие значения постоянной составляющей q° и коэффициента гармони ческой линеаризации:
- 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