
книги из ГПНТБ / Кулесский, Р. А. электропривод постоянного тока с цифровым управлением
.pdf(рис. 2-5,а). При этом выходная координата регулятора Для дис кретных моментов времени определяется выражением
s {пТ)]= c,s„ («7') + |
c2sn (пТ) + c3Tsa (пТ) = |
|
|
||||||
|
Т )1 + |
с |
2 |
|
|
л-1 |
|
(кТ). |
|
|
|
|
2 |
|
|||||
= с, [в* (пТ) - е* (пТ - |
|
|
е * |
(пТ) + |
с, Г |
А = 1 |
в* |
|
Умножение сигналов sR(nT), sn(nT) и sB(nT) на коэффициенты Сі, с2, сз осуществляется в аналоговой форме блоками ЦАП, кото
рые также производят и суммирование корректирующих сигналов
Рис. 2-5. Блок-схемы ЦВУ при параллельном (а) и прямом (б) про граммировании вычислений.
NKX(t) и после их преобразования непрерывными фильтрами \Ѵщ(р), где і= 1, 2, 3. Если же привести (2-16) к виду (2-3), то
принимая, например \^ф(д) = 1, из (2-15) получаем:
W*Hz) = ö0z2 + a,z + а2 |
(2-18) |
где ао = сі+ с2; ах=СзТ—с2—2сі; а 2 = С і .
Передаточной функции (2-18) соответствует разностное урав
нение
s(nT) —anе* (пТ) +аіе* (пТ—Т ) +
+ а2е* (пТ—2Т) + s (пТ—Т).
Согласно последнему ЦВУ можно построить из цифровых бло ков умножения, задержки на один период прерывания БЗ и сумми
рования (рис. 2-5,6).
После решения задачи синтеза, когда становятся из вестными передаточная функция оптимальной системы и закон изменения е(0 в линейной зоне, оценка влия-
60
иия помехи Ѵц (t) на работу электропривода может быть осуществлена на основе структурной схемы, приведенной на рис. 2-6. Выбирая в качестве индикатора степени
У ѵ ( п Т )
Рис. 2-6. Структурная схема для оценки влияния помехи от кванто вания по времени.
влияния помех уровень флюктуации некоторой 'коорди наты объекта управления yv(t), определяем передаточ ную функцию Я (р) линейного фильтра из условия
Н (р) = Y v (p)/Y (р). |
(2-19) |
Если yv( t ) = i ( t ) , то согласно изложенному |
в § 1-2 |
для электроприводов, предназначенных для регулирова
ния скорости и положения, выражения для Я(р) |
будут |
|
иметь вид: |
|
|
HJp) = |
Ta pl р; |
(2-20) |
Ha{p) = |
Tu fflpka . |
(2-21) |
С учетом (2-19) можно записать:
Yv {z,q) = ^ * { z ) \ ^ \ { z ) ^ Z
z — 1 7 гг д (р) |
(р) а (р) н (р) |
z |
PW^P) |
Г^ ф ІР) G (Р) Н |
CP) I |
|
1 |
Pwv (Р) |
J |
]}= E **(z) К* (*,?).
(2-22)
Определение yv(tiT+qT) по (2-22) может быть осу ществлено известными методами [Л. 1, 3, 6]. Так как в* (■£)— квантованная по уровню произвольная функция времени, удобнее определять yv(nT, qT) путем решения разностного уравнения, соответствующего (2-22). Если передаточная функция -К* (z, q) может быть представле на отношением двух полиномов от z
К* (z, q) |
(г, д) |
g0zm + ахг-т-х-\- |
•• + a m |
, (2-23) |
|
Е** (z) |
z ‘ + ö,21_1 + .. |
+ bi |
|
61
то соответствующее ей разностное уравнение имеет вид:
Уѵ(пТ, qT)= а0гЩп—1 + пі) Т] +
+ аіе*[ (п—1+ пг—1) Т] + ... + ат е*{(/г—/) Т]—
—biyA(n—l)T + q T ] - ... —b,y£(n—l)T + qT\. (2-24)
Выражения (2-23) п (2-24) могут быть использованы при вычислении уѵ(пТ, qT) при l^-ni. Если 'W*R(z) вы брана из условия ѵц(пТ)=0, то в установившемся режи ме при I е(оо) I <>а в* (оо) = 0 и Гц(оо) =0.
Пример 2-2. Определим характер изменения vn(t) при следую
щих исходных данных:
ß(t)=e~'\ \Ѵя (р) = 1/р\ \Ѵф(р)=.1; |
сг=0,1; |
Г=0,1; |
алгоритме пре |
||
образования |
e(/) в цифровую форму |
(I-3) |
н двух |
типах №*д (г): |
|
|
^ д . ( г ) = T / ( z - 1); |
\F V _ ( z) = |
Tz/(z- |
!). |
|
Принимая G(p)H(p)jWp{p) = l |
и |
используя табл. 2-1, из (2-22) |
|||
определяем |
выражение для K * i ( z , |
q) |
при \Ѵ*я (г) = \Ѵ*л і (г): |
K*(z,q) = K*l |
---- |
||
Z —■1 |
Tz |
= — qT. |
|
z |
[ ( 2 - I)2 + 4 |
||
|
Согласно (2-24) получаем:
V m ( i i T + q T ) = — q T s * ( n T ) .
При \V*R(z)=W *v{z) имеем:
T ( l - q ) z - T ( l - q )
K*(z,q)=K* 2 (z,q) = |
z — 1 |
|
|
|||
|
|
|
|
|
|
|
Тогда с |
учетом (2-24) можно |
записать выражение |
|
|||
|
|
Оц2 (пТ+ qT) = Т (1—<7 ){в* (пТ) — |
|
|
||
|
|
- е * |
(пТ—Т)]+ у9(nT—T + q T ). |
|
|
|
Графики |
Цці(t) и |
Ѵцг(і) представлены на рис. |
2-7. Как видно |
|||
из рисунка, |
|
при передаточной функции \V*n(z) = lK*Äi(z), |
удовле |
|||
творяющей |
(2-15), установившееся |
значение помехи |
vni (t) |
от кван |
тования по времени равно нулю. В то же время установившееся значение ццг(0 отлично от нуля.
Количественная оценка влияния компенсации помех квантования на динамические, свойства системы, как указывалось, является задачей уточненного анализа, осуществляемого методами моделирования. Однако
62
в некоторых простейших случаях для приближенной оценки динамических свойств электроприводов с цифро вым управлением может быть использован метод фазо-
s(t)
1,0
0,8
0,6
о,If
0,2
Рис. 2-7. Расчетные кривые к примеру 2-2.
вой плоскости {Л. 8]. Воспользуемся им для анализа ра боты позиционного электропривода. Если принять W,l(p) = l, lFa(p )= 0 и 0, что характерно для пози ционного электропривода в режиме отработки заданных
перемещений, |
а |
переда |
|
|
||
точную |
функцию |
после |
|
|
||
довательно |
соединенных |
|
|
|||
фильтра |
й/ф(р), запоми |
|
|
|||
нающего элемента |
W3(p) |
|
|
|||
и |
объекта |
управления |
|
|
||
принять равной |
Рис. 2-8. Структурная схема при |
|||||
|
|
|
|
|||
|
|
|
|
приближенной |
оценке динамиче |
|
|
|
|
|
ских |
свойств |
электропривода. |
то |
структурная |
схема |
к виду рис. 2-8, допу |
|||
рис. 2-1,а может быть приведена |
||||||
скающему относительно простое |
использование метода |
|||||
фазовой плоскости. |
|
(рис. 2-8) имеет вид: |
||||
|
Уравнение движения системы |
|||||
|
|
|
Л/ + У= k (х* — У*)■ |
(2-26) |
Заменяя рассмотрение движения системы под дейст вием скачка х рассмотрением свободного движения при
63
ненулевых начальных условиях, принимаем л; = 0 и пере ходим к уравнению
*У-\-У = — ky*> |
(2-27) |
где у * согласно (1-3) может принимать значения 0, ± о ,
± ’2сг ...
Приведя (2-27) к виду |
|
|
|
|
|
тz- \ - z = |
— ky*; |
|
|
(2-28) |
|
ij = |
z |
|
|
|
|
|
|
|
|
||
и разделив уравнения (2-28) |
одно на |
другое, получим |
|||
дифференциальное уравнение |
|
|
|
|
|
dz/dy = — (ky* + z)lxz. |
|
|
(2-29) |
||
Решением (2-29) является уравнение фазовых траек |
|||||
торий |
|
|
|
|
|
У^Уо + i Z o - t z - k y l n |
+ |
. |
(2-30) |
||
характеризующее движение |
изображающей |
точки |
в плоскости параметров z, у в интервалах между сосед ними значениями у*, отличающимися на дискретную единицу сг. Поэтому фазовая плоскость, отображающая движение, будет состоять из ряда листов, каждый из которых заполняется кривыми (2-30). При этом у каж дого листа (см. рис. 2-9,а) у * принимает одно из сле дующих значений 0, ± 0 , ±2ст, ±Зо ... В свою очередь траектории движения состоят из отрезков 'кривых (2-30) при соответствующих значениях у*. На рис. 2-9,а пока заны фазовые траектории, а на рис. 2-9,6 — кривые пере ходных процессов, построенные по соотношению (2-30) для т=0,5 с, о=1 при двух значениях k. Фазовая траек тория 1 и соответствующие ей переходные характеристи
ки yi(t), Zi(t) |
построены при k —\, а аналогичные кри |
вые 3 и уз(і), |
z3(t) — при k — A. |
При компенсации помех квантования структурная |
схема рис. 2-8 преобразуется в схему рис. 2-4, в которой передатоиная функция последовательно соединенных ре гулятора и объекта — R{p). Рассматривая свободное движение при ненулевых начальных условиях, записыва ем уравнения движения в виде
У + ~ ! г У + ~ У = У + ^ (О'у + ^У = 0, |
(2-31) |
64
Приводим дифференциальное уравнение второго по рядка к системе двух уравнении первого порядка
z -ф- 2уш.г -ф- игу = 0;
(2-32)
y = z.
Дифференциальное уравнение фазовых траектории
dz __ |
2уи ? + &-у |
(2-33) |
|
7/г/ |
г |
||
|
полученное делением уравнений (2-32) одного на другое, является однородным и может быть решено методом разделения 'переменных после замены u= z/y.
Рис. |
2-9. Фазовые траектории (а) |
и кривые переходных |
процес |
||
сов |
(б) для системы |
рис. 2-8. |
|
|
|
Решение (2-33) при представляющих наибольший |
|||||
интерес значениях у<1 |
дает |
уравнение фазовой |
траек |
||
тории |
|
|
|
|
|
|
(2 + W )2+ |
• y2 = |
Cexp{^=drctg -^ Ц Ш ,(2 -34) |
||
где |
Д= 4со2(1—у2). |
|
0,5ijVh |
f |
|
|
|
|
|||
|
Постоянная интегрирования С находится при подста |
новке начальных значений координат изображающей точки z = z 0, у = уо- На рис. 2-9,а показаны фазовые тра-
5—181 |
65 |
екторпи, а на рис. 2-9,6 — кривые переходных процессов, построенные по (2-34) при тех же значениях т и к, что и в цифровой системе без компенсации помех квантова ния. Фазовая траектория 2 и соответствующие ей пере
ходные |
характеристики ih(t), |
z2(t) построены |
при /е = 1, |
а аналогичные кривые 4 и |
zi(t) — при |
/г = 4. Так |
|
как для |
системы рис. 2-8 U7*Ä'(z) = W:{(p) = 1, то соглас |
но (2-13) Оц(/)=0. В связи с этим фазовые траектории (2-34) характеризуют движение при цифровом управле нии с компенсацией помех квантования и при отличном от нуля значении периода прерывания Т.
Сопоставление |
кривых |
переходных |
процессов |
(рис. 2-9) показывает, что для |
системы рис. 2-8 компен |
сация помех квантования позволяет устранить отрица тельное влияние квантования сигналов по уровню и приблизить динамические свойства систем с цифровым управлением к динамическим свойствам систем при не прерывном управлении. Необходимо отметить, что пред ставление структурной схемы позиционного электропри вода при его работе в линейной зоне изменения коорди нат в виде рис. 2-8 возможно лишь при целом ряде упрощающих допущений относительно структурной схе мы объекта (рис. 1-10). В связи с этим метод фазовой плоскости применим здесь лишь как метод качественно го исследования динамических свойств.
2-3. ОСОБЕННОСТИ КОМПЕНСАЦИИ ПОМЕХ КВАНТОВАНИЯ ПРИ НЕЛИНЕЙНОЙ ХАРАКТЕРИСТИКЕ ИЗМЕРИТЕЛЯ ПРОИЗВОДНОЙ РЕГУЛИРУЕМОЙ КООРДИНАТЫ ОБЪЕКТА УПРАВЛЕНИЯ
При рассмотрении в § 2-2 вопросов компенсации помех от квантования сигнала в обратной связи пред полагалось, что датчики производной ДП в схемах фор мирования корректирующих сигналов NK (рис. 2-2) обладали линейной характеристикой. Вместе с тем в ря де случаев это не выполняется. Например, в электропри воде, предназначенном для регулирования скорости (рис. 1-9), динамический ток
с, = |
ш |
(2-35) |
|
д |
Р<Ро |
ѵ |
’ |
связан линейно с производной скорости лишь при по стоянном значении магнитного потока возбуждения сро- В соответствии с этим при изменяющемся потоке воз буждения датчик динамического тока при использовании
66
его в качестве ДП будет обладать нелинейной характе ристикой. Для позиционного электропривода скорость двигателя линейно связана с производной утла поворота выходного вала редуктора лишь при отсутствии люфта, сухого трения, проскальзывания и т. п. Только при этом условии датчик э. д. с. (скорости) двигателя, используе мый в качестве ДП, будет обладать линейной характе ристикой. Нелинейность характеристики ДП может при водить к нарушению условий компенсации помех кванто вания и должна быть учтена соответствующими измене ниями в схемах рис. 2-2. Могут иметь место два вариан та преобразования последних.
Первый вариант предусматривает при нелинейной за висимости tj = F{yn.„) включение на выход ДП в схемах рис. 2-2,а, 6 нелинейного преобразователя с характерис тикой Г= F. Так, при использовании датчика ід в качест ве датчика производной скорости двигателя нелинейный преобразователь согласно (2-35) будет представлять собой блок умножения г/д.1Тна значения сигнала ср0р/Гм.
Если нелинейная зависимость связывает значения интегралов ^ у dl — ѵр j удЛ1 dt, что имеет место, напри-
/t
мер, в позиционном электроприводе при нелинейности в механической передаче, то функциональный преобра зователь с характеристикой Q = VF включается между выходом БИ и входом БНУ (рис. 2-2,6). Использование схемы рис. 2-2,а в этом случае нецелесообразно. Изме нение в соответствии с изложенным функциональных схем рис. 2-2 позволяет обеспечить компенсацию помех квантования с той же точностью, как и в случае линей ной характеристики ДП, так что эквивалентная схема рис. 2-4 сохраняет свою силу. Однако для некоторых ти пов нелинейных зависимостей Y можно достичь улуч шения динамических и статических характеристик без изменения схемы рис. 2-2,6, т. е. не вводя преобразо ватель с характеристикой Q. К числу таких зависимос тей относятся зависимости типа зоны нечувствительно сти, люфта и т. п.
Рассмотрим позиционный электропривод при наличии люфта в механической передаче. Структурная схема электропривода представлена на рис. 2-10. Наличие ком пенсации помех от квантования сигнала обратной связи отражено введением на вход сигнала Выражение
5’ |
67 |
для R(p) соответствует (2-25). Для анализа восполь зуемся методом фазовой плоскости при следующих до пущениях:
а) маховые массы, присоединенные к выходному в лу редуктора, достаточно малы, так что после выбора
Рис. 2-10. Структурная схема позиционного элек тропривода при наличии люфта в механической передаче.
люфта выходной вал мгновенно приобретает скорость
входного; |
компенсация помех квантования в цепи обратно |
б) |
связи осуществляется устройством рис. 2-2,6 и в качест
|
|
ве |
ДП |
используется |
датчик |
|||
|
|
э. д. с. двигателя; |
|
|
||||
|
|
в) |
|
|
|
|
||
|
|
люфт выбрани г/го = 2Уо = ?/о —до |
||||||
|
|
определим условия возник |
||||||
|
|
новения |
периодических |
режи |
||||
|
|
мов. |
Свободное движение в |
|||||
|
|
плоскости у, z при ненулевых |
||||||
|
|
начальных условиях у~уо, z = |
||||||
|
|
= z 0 |
и |
компенсации |
|
помех |
||
|
|
квантования |
характеризуется |
|||||
|
|
(2-34). Принимая для опреде |
||||||
|
|
ленности |
положение |
изобра |
||||
Рпс. 2-11. Фазовые траекто |
жающей точки 0 с координата- . |
|||||||
ми го, уо в четвертом квадран--. |
||||||||
рии |
движения в системе |
|||||||
рпс. |
2-10. |
те, можно |
построить |
участок |
||||
|
|
0/1 |
фазовой |
траектории |
z(y) |
|||
(рис. 2-11). В точке А -координата zT= yr изменяет |
знак |
и начинается вы-бор люфта. При этом t/*=const и, как следует из функциональной схемы рис. 2-2,6, импульсы Са не вырабатываются, а ДП с блоком интегрирования
68
БИ обеспечивает жесткую обратную связь по координате у,- Поэтому движение изображающей точки в плоскости zr, ут будет по-прежнему определяться (2-34). Если со вместить оси координат z, у и zr, у,-, то фазовая траекто рия zr, у,- будет изображаться отрезком АЕ при том же значении постоянной С=СПв (2-34), что и 0/1.
В дочке D, имеющей координаты г,-о и у,-2= г/і+ 2б,
заканчивается выбор люфта. При этом выходная коор
дината у |
мгновенно |
приобретает |
скорость у г, т. е. |
z = z r |
||||||
п изображающая точка (г, |
у) попадает в точку F. |
Одна |
||||||||
ко очередной |
импульс |
Со |
моиіет появиться лишь в мо |
|||||||
мент равенства у целому |
числу |
единиц |
дискретности |
|||||||
а. |
Если |
уі = Ігв+ ао, |
где |
а< 1, |
/г = 0, —1, —2 |
..., то |
||||
импульс |
Сг |
появится |
в момент |
времени, |
когда |
изобра |
||||
жающая точка (г,-, уг) |
совпадает с точкой Е, имеющей |
|||||||||
координаты гг3, y r3= |
(/i+ 2ft-fа(1— и). От |
момента |
окон |
|||||||
чания выбора люфта до появления импульса С0, |
фазо |
|||||||||
вой |
траектории z ( y ) |
соответствует отрезок F G , |
удовле |
творяющий (2-34). В этом случае постоянная интегриро
вания С= С[ и определяется из |
(2-34) подстановкой |
в него координат точки F. |
|
При появлении импульса Са с |
помощью БИУ на вы |
ходе БИ устанавливается значение сигнала, равное ну лю, и для системы рис. 2-10 становится справедливой эквивалентная схема рис. 2-4. Фазовая траектория z(y) представлена для этого случая отрезком GH и опреде ляется уравнением (2-34) при том же значении постоян ной С=С 1 , что и участок FG. В точке Н zT изменяет
знак, начинается следующий цикл выбора люфта и даль нейшая процедура построения фазовой траектории ана логична описанной выше. Таким образом, фазовая траектория z (у) будет представлять собой кривую
OAFHIK ...
При равных значениях координат точек пересечения фазовой траектории с осью ординат, т. е. при \zi\=Z2,
возможно возникновение периодического цикла — авто-» колебаний. Определим условия их существования. Под ставляя в (2-34) г/= г/і и г = 0, находим:
С0 = |
d•exp j —- arclg |
j , |
(2-36) |
где d = ш; b — yur, |
d i = d 2— b 2. |
|
|
69