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

книги из ГПНТБ / Кулесский, Р. А. электропривод постоянного тока с цифровым управлением

.pdf
Скачиваний:
13
Добавлен:
20.10.2023
Размер:
8.47 Mб
Скачать

влияния помех 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

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