
книги из ГПНТБ / Кулесский, Р. А. электропривод постоянного тока с цифровым управлением
.pdfвлияния помех Q{t). На основании принципа суперпози ции запишем:
y(t) —t)w (t) +«/Q(0 .
где Ljw(t), г/Q ( /) — реакции системы рис. 2-23,6 соответ ственно на W(t) и Q(t).
Подобно рассмотренному в § 3-1 качество отработки полезного сигнала можно оценивать значением интеграла
Vi = ] [ y 9{ t ) - y v (t)}3dtt |
(4-1) |
6 |
|
которое определяет степень приближения |
переходного |
процесса отработки сигнала W(t) к желаемому процессу
y3(t) |
в эталонной системе. Влияние помех будем |
оцени |
|||
вать |
величиной |
дисперсии |
флюктуаций |
тока |
якоря |
.М[і2 (0] двигателя, |
вызванных |
действием |
Q{t). |
В соот |
ветствии с этим за критерий качества можно принять ми нимум интеграла (4-1) при ограничении среднего уровня
флюктуации тока якоря допустимым |
значением |
[Л. 37] |
V, = rain при М [г (0] = |
t'on. |
(4-2) |
Следует отметить, что ограничение М[і2(/)] не только улучшает энергетические характеристики электроприво да, но также снижает вероятность значительных по ве личине скачков ускорения, приводящих к механическим ударам в передаче, т. е. повышает надежность работы электропривода.
Структурная схема рис. 2-23,6 была разработана для следящего электропривода, хотя, как отмечено в § 2 -6 ,
она принципиально может быть распространена и на программный электропривод при некоторых типах по лезных сигналов. Однако такой привод детально не ис следовался, поэтому ниже все изложение будем строить применительно к следящему электроприводу. Примем, что период прерывания Т удовлетворяет неравенству
Т < 71 |
(4-3) |
С»С.т 1 |
|
где Ы .т — значение частоты, выше которой амплитудная частотная характеристика, определяемая как отношение амплитуды синусоидальных колебаний тока якоря к амп литуде задающих колебаний, приложенных к точке дей ствия помех, близка к нулю.
При выполнении (4-3) в соответствии с теоремой Ко тельникова (Л. 1] можно пренебречь импульсным харак-
130
те р е м сигналов в схеме рис. 2-23,6 и перейти для подав ляющего большинства промышленных электроприводов к структурной схеме рнс. 4-1. Действительно, при приня том в настоящее время подходе к проектированию вен тильных электроприводов верхнюю границу полосы про пускания контура регулирования тока якоря ограничи
вают значением |
іос.т = (1504-200) с-1. |
В то же время |
|
’ К ® |
h U) |
|
|
Ф ( р ) |
Ф(р) |
|
|
,R~(t) Ruit) |
|
cc(t) |
|
W(t) |
Wa(p) |
Wp(p) |
|
|
О б ъ е к т |
||
|
Wa(p) |
|
|
----------------------------------- |
G(p),g(v) ■ |
|
|
Рис. 4-1. Расчетная |
структурная |
схема. |
|
период прерывания основных аналого-цифровых преоб разователей позиционных электроприводов, как отмечено в § 1-5, не превышает 0,01 с. Вместе с тем после прибли женного синтеза может быть осуществлена проверка вы полнения (4-3) при окончательных значениях параметров регулятора. Для некоторых типовых решений допусти мые из удовлетворения (4-3) соотношения между пери одом прерывания и значениями параметров регулятора приводятся в § 4-6.
В соответствии с (4-2) искомый оптимальный регу лятор в системе рис. 4-1 обеспечивает наилучшее при ближение процесса отработки полезного сигнала к же
лаемому процессу |
в эталонной системе при допустимом |
|||||
в соответствии с |
величиной |
і"оп |
уровне |
флюктуации |
||
тока якоря, вызванных действием |
помех |
квантования. |
||||
Выбор передаточной W3(p) |
или |
весовой w3(т) |
функций |
|||
эталонной системы производится |
аналогично |
изложен |
||||
ному в гл. 3. |
|
|
|
|
|
|
Корреляционные функции Rx(x), |
Ry (т) и Q(t) в соот |
|||||
ветствии е § 2-6 и 2-7 определяются |
(2-92) |
и (2-96). По- |
9* |
131 |
лезпый сигнал W(t) согласно формулировке критерия качества (4-'2) используется как «пробный» для про верки степени приближения проектируемой системы к эталонной. В силу линейности обеих систем по условию целесообразно зафиксировать значение и, приняв его для наиболее тяжелого случая равным 1. Значениями ц и ѵ/2 в выражении для W(t) (см. § 2 -6 ) можно пренебречь.
Погрешность, вносимая при этом, не преззойдет допусти мой статической, а с точки зрения сравнения эталонной и проектируемой систем еще меньше. Тогда
W(t) = U2t. |
(4-4) |
Значения \ \ и М[і2(/)] в (4-2) при заданных харак теристиках объекта управления и известных характери стиках W(t), Rx{t) и Ry(t) зависят от весовой ^(т) или передаточной G(p) функций системы рис. 4-1. Определе ние g(x) и G(p), обеспечивающих выполнение (4-2), представляет собой задачу на условный экстремум, кото рая решается ниже методами классического вариацион ного исчисления. Искомое выражение весовой шр(т) или передаточной Wp(p) функций регулятора при известных g(i), G(p) определяется однозначно.
Перейдем к выводу уравнений, определяющих опти мальную в соответствии с (4-2) систему управления элек троприводом, характеризующуюся передаточной функци ей регулятора Wv {p) = W a(p) + W R(p)Wlb(p).
4-2. УРАВНЕНИЯ, ОПРЕДЕЛЯЮЩИЕ ОПТИМАЛЬНЫЙ СЛЕДЯЩИЙ ЭЛЕКТРОПРИВОД
Представим (4-2) |
в виде |
|
|
|
|
00 |
|
|
|
V= |
J 5, (() dt + гМ [г (0] = |
min, |
(4-5) |
|
где |
О |
|
|
|
|
|
|
|
|
(^)= Уа(0 |
2уэ (t) yw (t) -|- yw(t)\ |
(4-6) |
||
|
‘(0= t[Q *(0+ Q y(01; |
|
(4-7) |
|
г — неопределенный множитель Лагранжа. |
|
|||
С помощью |
(4-5) рассматриваемая |
задача |
сводится |
к определению условий существования безусловного экс тремума. Ход решения вначале включает определение принципиального решения для g( т) или G(p), достав-
132
Мющих минимум (4-5) при неопределенном значении При этом g(x) и G(p) являются функциями множителя г. Значение последнего определяется затем из дополни тельного условия M[i2 (Z)]= і2Д0П, после чего g( т) и G{p)
получают свой окончательный вид, однозначно опреде
ляющий передаточную Wv (p) или весовую &ур(т) |
функ |
|||
цию регулятора |
при известной передаточной |
функции |
||
W0(p) объекта. |
|
|
|
|
Пусть объект управления характеризуется-передаточ |
||||
ной функцией |
|
|
|
|
• |
W0( p ) = ^ W K.c (p), |
|
(4-8) |
|
где ИД.0 (р )— передаточная |
функция замкнутого |
конту |
||
ра регулирования скорости. |
весовые функции |
элементов |
||
Выразим Л4'[і2(^)] через |
системы рис. 4-1. Для этого сначала найдем изображе
ние і(/) |
по Лапласу. Как следует из § 1-3, |
|
|
|
|
i(p) = j±- p- a(p)=H(p)*(p) . |
|
(4-9) |
|
В свою очередь из схемы рис. 4-1 можно получить: |
||||
<х(р) |
|
( Р ) + Q y H l = Г Т Т О ) ф { Р ) Е |
|
|
т ^ ж Ь ) ф ( Р ) |
{ р ) ' |
|||
где |
|
|
|
(4-10) |
|
WAP) w*{P)w°iP) |
|
|
|
|
Я(р) = |
|
|
|
Преобразуем (4-10): |
|
|
|
|
- (п\ _ |
(?) W<bІР) Wo ІР) [^a ІР) + ИД ІР) |
(jo)] |
\ / |
|
КР) ~~ « Д (р )+ И Д (р ) wt ( р) (1 + И Д (р) [\Ѵа {p)+Wд ( » Й7ф ( я ) ) } * |
||||
|
X Ф (Р) Б (Р) = |
{Р) Ф (р) G(Р) Е (р). |
(4-11) |
Подставляя (4-11) в (4-9), получаем:
L(P) = G (р) ф (р) Н (р) X
X WAwJ p ) {P) ■Е (Р) = °(Р)р (Р) Е (Р)- (4-12)
Отметим, что при чисто цифровом регуляторе поло жения Wn(p) = 0 и Wv (p) =ѴРл(р)Ѵ?ф(р). При этом Л(р)=Я(р)Ф'(р).
133
ЕСЛИ №а ( р ) # '0 то определяя \Wp(p) из соотношения
G(P) _ |
Wo (p) |
(p) |
1 |
+ i'7o (p) |
(p) |
H подставляя его в выражение для F{p), получаем:
F (р) = \ѴЯ (р) Г ф (р) W0 (р) Ф(р)Н(р)[0-'(р)-\}. (4-13)
В соответствии с |
(4-12) |
|
||
|
со со |
|
|
|
IW = |
5 |
(У) / (4) [Q*V, У+ 4) + Qy (t. Y + 4)] dl dl,. |
||
|
ОО |
|
|
|
Используя (3-25), запишем: |
|
|||
/И [ г |
({.)] = |
М ( я (Y) % g (а) % / (7)) >1с / (в) % |
||
^ |
[Qx{t. У+ 4) -\~Qy У’У Ч~ 4)] [QxO*. а + |
s) -Ь |
||
|
|
|
+ Qy (t, а -)-£)]}. |
(4-14) |
Для дальнейших преобразовании изменим порядок операций интегрирования и определения математическо го ожидания в (4-14) и рассмотрим выражение
/ = M{[Q.V(/, у + л) +Q„(t, у+гі)]Х
X[Q.*(/, а + е) + Q y (t, а+е]}. |
(4-15) |
Осуществляя умножение под знаком математического ожидания и учитывая, что автокорреляционная и взаим ная корреляционная функции равны:
Кпп (т) =M[n(l )n(t +т)];
Кп т СО =M[n(t) т (/+т)],
представим (4-15) в виде
І = Кхх(у + г[, а+ е) +Кхѵ(-Y + Л. а + е) +
+Дд:у(а+8 , Y + T])+/Cra('Y + i], а + е).
Всилу независимости случайных сигналов Qx(t) и
Qy(t) корреляционные функции Кху (т) = Кух(х) = 0 . Если
учесть, что согласно изложенному в § 2-7 в |
расчетной |
схеме рис. 2-22,6 характеристики Rx(t) и Ry (t) |
совпада |
ют, то Л+Дт) = К ѴУ{т) = K QQ{X). При этом |
|
/ = 2/CQQ(Y + TI, а + е). |
(4-16) |
134
Подставляя |
(4-16) |
в (4-14), |
преобразуем |
последнее |
||||
к виду |
|
|
|
00 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
M [‘s(9] = |
2 jS 3 |
te.Y)rfr. |
(4-17) |
||||
где |
|
|
|
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
(g . Y) = |
g (Y) g (a) * |
f h) * |
f 0) * |
|
||||
|
* |
/,CQQ(Y + |
4.a + |
s)- |
|
(4-18) |
||
Согласно изложенному выше, |
в |
случае чисто цифро |
||||||
вого регулятора |
f (t) = |
h (Я) |
<р(х — Я). |
Если |
регулятор |
положения цифро-аналоговый, то в соответствии с (4-13) необходимо принять
fW = |
g-(Y )*M ß)*«>o (Ю*<Р(Я)* |
% Щ (а) щ)д (х — Y — ß - р. — 3 — я) — Шф (а) >(с |
|
Л (ß) % |
М % <р(Я) Шд (т: — р — И- — а — Я). (4-19) |
Здесь g~(у) весовая функция системы, обратной по отношению к системе с весовой функцией g(y).
Далее из (3-26) записываем входящее в (4-5) выра жение для Si, заменив x(t) на W(t)\
S , = W B (х) % Щ (з) % ^ ( t , *) Г (Л з) -
- 2 ® 3 M * £ ( Y ) > W . W * . Y) +
С учетом этого преобразуем первое слагаемое в (4-5),
С О |
СО |
|
S, ( t ) d t = |
^SJ2 (g,Y)dY> |
(4-20) |
6 |
oJ |
|
где |
|
|
Sn (ша, x) = w3(x) ou3 (o) >jc / 1F1F (x, o); |
|
|
(g. 4) = - 2 g (Y) |
(x) % / rft, (Y, *) + |
|
~Ь g (Y) g (°0 |
I\\7\\7(Y> ®)> |
|
|
e). |
|
135
На основе (4-20) представим (4-5) в виде
00
V = { d Y lS 13te.Y) + 2rSate.Y)] +
О
со |
|
+ I’d'cSn (шэ, х) = шіп. |
(4-21) |
Ь
Искомая весовая функция g(y), доставляющая мини мум (4-21), определяется из решения уравнения Эйлера
(4-22)
Здесь в соответствии с (4-20) и (4-21)
ß, = 2^ (а) % Iww (Y. а) - 2шэ (х) * Iv v (Y, х). (4-23)
Для В2 могут быть получены два выражения. В слу чае чисто цифрового регулятора, производя в (4-18) за мену /(п) =!і(І) >!<ф(т)—Я), получаем:
В., = 2g (а) ^ h (ß) >{< h (р) |
<? (ц - |
ß) % |
|
% Ч*(s — p) % /Cw (Y + |
-Ч.<* + |
e). |
(4-24) |
При цифро-аналоговом регуляторе положения f(r\) определяется (4-19). Подставляя последнее в (4-18), по сле преобразований находим:
В2 — 2g (а) ^ h (ß) j)< f (Я) % w0(|х) % Шф (а) % |
|
||||||
WA (Ъ ß + |
1 1 + 3 + |
*) % h (р) |
<р(0 ) % w0 (ѵ) % |
|
|||
^ щ {у.) % WA (s> р ~Ь v + z - Ь б) K QQ (Y ~Ь і)- а - | - £) - |
|
||||||
—2h (ß) % <р(Я) ^ ш0 |
(р.) % Шф (з) %tBA(T|,ß+ |
fi + |
|
||||
+ ° + |
А) |
(Р) * |
? (б) % |
( ѵ ) % |
а»* (х ) |
% |
|
* |
®д (е. Р + V+ X+ в) ^ QQ (Y+ |
-Ч. е )- |
|
(4-25) |
|||
Далее приведем (4-22) к виду |
|
|
|
|
|||
СО |
|
|
da — г (Y) = |
|
|
|
|
j & (а) /г (Y- а) |
0 при у > 0. |
(4-26) |
|||||
О |
|
|
|
|
|
|
|
із§
Здесь: а) для чисто цифрового регулятора |
|
||||||
/е (у, а) = І „ . (у, а) + 2 |
rh ф) % h (р) % ? (ц — р) |
||||||
|
% <Р(s — Р) % KQQ(Т + |
"Ч.а + |
е); |
(4-27) |
|||
|
r(4) = |
w3( , ) * I ww(b-)\ |
|
(4-28) |
|||
б) для цифро-аналогового регулятора |
|
|
|||||
k (у, а) = |
Ivw (у, а) + |
2г/г (р) % /г (р) >j< <р(Я) ^ |
|||||
% <Р(6) >t< |
М |
% |
( ѵ ) % Щ (а) >{< даф ('/) |
^ |
|||
* ® д (і. ß —{—JJ- — |
—j—Я) >j<оУд(е, р + ѵ+ |
х + |
0) * |
||||
|
|
|
(Т + 1]. а + е); |
|
(4-29) |
||
г (у) = щ (х) * |
Ітѵ (у, х) + 2 |
rh (р) * |
h (р) * |
||||
>к <РW |
>)< ¥ ( Ѳ) >1< а і0 (И-) % |
( ѵ ) |
3fc иуф (а) % % |
(х ) >)<] |
|||
* ®д fa. ß + |
Р + |
° + *) |
“ »д (е, р+Ѵ + Х + 6 ) >к TCQQ (у+т], е). |
||||
|
|
|
|
|
|
|
(4-30 |
Определение весовой функции замкнутой оптимальной системы привело к интегральному уравнению (4-'26) типа Винера—Хопфа. Решение этого уравнения в комплексной плоскости дается формулой (3-34). Оно будет осущест вляться достаточно просто, если Ф<(р), являющееся изоб ражением по Лапласу ср(т), может быть представлено дробно-рациональной функцией. Иначе процесс опреде ления [С]+ и [С]- согласно (3-34) может быть настолько сложным, что проще окажется найти приближенное ре шение (4-26) во временной области [Л. 23]. Таким обра зом, после определения Ф(д) необходимо осуществить
аппроксимацию ее дробно-рациональной |
функцией. |
В практически встречающихся случаях это |
возможно |
(Л. 12]. |
|
Для определения входящих в (3-34) компонентов осу ществим двустороннее преобразование Лапласа левой и
правой частей уравнений |
(4-27) — (4-30). Учитывая, что |
|
спектральная плотность |
СО |
|
|
|
|
S-'nn (р) = i |
j Knn (■=) е~р~' dt, |
(4-31 ) |
137
Мосле осуществления операций преобразования полу чаем:
К(р) =Iww{p) + 4 пгН(р)Н{—р)Ф(р) X
X <D (-P )S*M (/7); |
(4-32) |
R(p) = W3(p)Iw w (p); |
(4-33) |
К(р) —Iww(p) +4пгН(р)Н( —р)Ф(р) X |
|
Х'Ф (—р) Wo(p) Wo( - p) \Ѵф(р) W ^ - p ) |
X |
X WA(p)WR( - p ) S % Q(p)- |
(4-34) |
R(p) = W3(p)!ww(p) + 4 л г Н ( р ) Щ- р ) X
X Ф (p) Ф i- р ) Wo (p) W0 ( - P) W* (p) Wb (- p ) X
X\VK(p)WA( - p ) S * QQ(p). |
(4-35) |
|||
Выражение для Iww(p), |
входящее в (4-32) — (4-35), |
|||
может быть определено на основании |
(4-4) |
|
||
1WV ІР)= |
U2 р2 (_ ру • |
(4-36) |
||
Формулы (3-34), (4-32) |
и |
(4-33) |
определяют |
G(p) |
в случае чисто цифрового, |
а |
(3-34), |
(4-34) и (4-35) — |
в случае цифро-аналогового регуляторов положения как функцию неопределенного множителя г. Поэтому далее
необходимо задать г |
из условия ограничения среднего |
|
уровня флюктуаций тока якоря согласно |
(4-2). Подстав |
|
ляя |
|
|
K QQ (Т + 4 .а + |
Е)= = j S *QQ е ^ + П “ |
e)Pdp |
—00
в(4-17), после преобразований получаем:
|
СО |
М [Г (01 = |
j G ( p ) G ( - p ) F ( p ) F ( - p ) S*QQ(р) dp. (4-37) |
|
— СО |
В (4-37) |
G(p) является функцией г. Приравнивая |
М [і2 (/)] из (4-37) ідш , нетрудно определить значение г,
подстановка которого в принципиальное решение для G(p) завершает поиск оптимальной структуры замкну той системы. После определения G(p) передаточная функция регулятора Wp(p) при известной W0(p) нахо дится однозначно. Далее согласно последовательности
138
приближенного синтеза, изложенной в § 2 -1 , следует вы
бор передаточных функций Wa(p), \Ѵл(р) и W$(p). Вы бор передаточной функции дискретного фильтра W*^(z), определяющей программу .работы ЦВУ, производится из условия эквивалентности ее передаточной функции W)X(p) непрерывного-фильтра. Методика решения такой задачи приводится в [Л. 1]. Перейдем к рассмотрению примеров синтеза регуляторов.
4-3. СИНТЕЗ ЦИФРОВОГО РЕГУЛЯТОРА ПОЛОЖЕНИЯ
Синтез оптимальной передаточной функции цифро вого регулятора производится на основе формул (3-34), (4-32) и (4-33). Передаточная функция замкнутого кон тура регулирования скорости, определяющая W0(p) по (4-8), в соответствии с общепринятым подходом (Л. 11] выражается как
^к.с(Р) |
1 |
(4-38) |
|
А д .о* ( 4 7 - ^ p - l- I) |
|||
|
’ |
где Т — некомпенсированная малая постоянная времени
контура регулирования тока якоря.
В большинстве практических случаев полоса частот помех квантования шире полосы пропускания системы, входом которой является точка приложения помех Rx(t), Ry {t), а выходом — точка формирования сигнала тока якоря двигателя і (t) (рис. 4-1). Передаточная функция этой системы определяется согласно (4-12) соотноше нием
'■(р)/Е(р) = G(p)F(p), |
(4-39) |
где F(p) =Н(р)Ф(р) .
Как отмечалось в § 2-6, при приближенном синтезе
целесообразно |
пользоваться |
выражением (2-95) для |
|
S R R {Ü>). |
Тогда |
согласно (4-39) |
левая и правая границы |
полосы |
пропускания фильтра |
G(p)F(p) должны нахо |
диться в интервале частот от 2nU\/o до 2JT£/2/O (рис. 4-2).
Выберем So в |
(2-98) равным 0,5 5Л/?(сй). При этом со |
||
гласно |
(2-97) |
в полосе пропускания |
фильтра (4-39) |
Ф'(р) = 1. Такой подход к выбору Ф(р) |
существенно уп |
||
рощает |
процедуру синтеза регулятора |
и соответствует |
широко применяемому на практике в подобных случаях. Однако при этом после определения окончательных вы ражений для Wa(p), И7Ф(р) необходимо осуществить про верку справедливости принятого предположения о соог-
139