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

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

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

(рис. 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

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