
книги из ГПНТБ / Кулесский, Р. А. электропривод постоянного тока с цифровым управлением
.pdfношении полосы частот Rx(t), Ry(t) и полосы пропуска ния фильтра (4-39). Если соотношение, иллюстрируемое кривыми на рис. 4-2, не выполняется, то Ф(р) должна определяться согласно (2-98) и процедура определения
Рис. 4-2. Амплитудная частотная
характеристика |
[G(іы) F (іа>) | и |
спектральная |
плотность S nn(ш) |
при Ф(р) = 1. |
|
№д(р), ТС'ф(р) повторяется. В данном примере примем Ф(/?) = 1, что справедливо для большинства практиче ских случаев. Перейдем к определению принципиального решения для G(p). Подставляя IWw(p) и Н(р) из (4-36) и (4-9) при Ф(р) = 1 в (4-32), получаем:
К (Р ) |
4кгГІ |
6 s + / > 8 |
(4-40) |
k2f |
рЧ- PY |
а >
Здесь
(4-41)
Согласно (3-34) представляет К(р) ів виде произве дения К~{р) и К+(р), где
К- ( р ) = - f i - ^ - 2 .6 'Ьа3 + 3,416Ѵ-2,616Ѵ+&^ . ^ ^
А«Ра
Д-+ ^ ___ р4 |
2,616/?34- 3 , 4 1 4 - 2,616*p + b ' |
(4-43) |
|
|
Для определения R(p) необходимо выбрать эталон ную систему. Как указывалось в § 3-1, выбор эталонной системы должен осуществляться исходя из требуемого качества отработки полезного сигнала. Примем в каче стве эталонной идеальную по отношению к полезному сигналу систему с передаточной функцией W3(p) = L Подставляя ее в (4-33), находим:
= |
(4‘44) |
140
Далее записываем выражение для R(p)/K~(p) и, раз лагая его на простые дроби в соответствии с (3-34), по лучаем:
г |
R(p) |
] |
__2,61 b * p + |
Ь* |
|
|
К - {Р) |
\ + |
Р°- |
|
|
При этом искомое выражение для передаточной функ |
|||||
ции приобретает вид: |
|
2,61 |
|
|
|
|
|
|
|
|
|
_ |
|
|
ъ р + 1 |
|
|
^ (Р)= ~ 1і |
<Гё ] |
з " Т і |
2Тб1 |
* |
|
|
2,61 |
3,41 |
|
|
Ь* |
6 3 |
6 2 • Рг ■ |
Р + 1 |
Введем для удобства обозначение т4 |
= & - 4 и преобра |
||
зуем /+46): |
|
|
|
G(р) = ----------------- 2 ,6 |
ltp_-j-j____________ (4-47) |
||
x*J r‘ + |
2,61t3Jö3+ |
3,41tV2 + |
2 , 6 1 x p + 1 ‘ |
Выражение (4-47) определяет принципиальное реше ние для G(p). Здесь неопределенным остается значение постоянной времени т=Ь~1, являющейся функцией неоп ределенного множителя г (4-41), который должен выби раться согласно (4-2), исходя из ограничения М[і2(£)]. Выразим М[і2(/)] через найденное решение для G(p). Подставляя (4-47) и F(p)=H(p) в (4-37), получаем:
---Tl So |
1 |
Р — |
6 ,82х- (ко)6 + |
(іа,)* . |
' (4-48) |
fkl |
2 я |
J |
/і4 (ісо) К ( - |
i«) |
гд е
hi(ia) = т 4 (/со) 4 + 2,61т3(ісо)3 +
+ 3,41 т2(гсо)2+2,61т(ісо) + 1.
Пользуясь для вычисления этого интеграла формула ми табл. 3-1, находим:
X-5. (4-49)
Р«а Приравнивая (4-49) i^on, получаем:
(4-50)
141
Здесь So должно быть представлено в относительных единицах. Для этого вычисленное по формуле (2-95) зна
чение необходимо разделить на |
• |
Передаточная |
функция разомкнутой оптимальной системы рис. 4 - 1 бу
дет иметь вид:
J ( p \ — |
6 |
— _______ 2 , 6 Н р + |
1________ |
|
|
I — G (p) |
3 ,4 Ь 2р2 (0,294т2р2+ 0 ,7 7 т р + 1)' |
||
Пример 4-1. |
Определим |
минимальное из |
условия ограничения |
|
Af[i2 (0] значение |
г. Примем і~оп = 0,0!, |
Т ы = 0,06 с, кл = |
=0,01 м -рад-1, р = 0,075, £2И= 78,5 рад-с-1. Значение So |
по данным |
||||
примера 2-5 выберем равным 4- ІО-8 м2-с. Подставляя в |
(4-50) не |
||||
обходимые значения, |
получаем: |
|
|
|
|
т |
4,76-3,14-(0,Об)2-4 - ІО-8 |
|
|
||
(0,075)2-(0,01)2-0,01 |
(78,5)2 = 0,144 с. |
|
|||
Вычисленное значение т подставляем в (4-51) |
и находим пере |
||||
даточную функцию разомкнутой оптимальной системы: |
|
||||
J(p ) = |
|
0,776/7+ |
1 |
|
|
0,071p2 (0.00614/72 + |
0,1 11/7 + |
1) ' |
|
||
Логарифмическая |
амплитудная |
частотная |
характеристика |
||
(ЛАЧХ) L(cö) =20 lg |/(іш) I приведена |
на рис. 4-3 |
(кривая 1). |
Рис. 4-3. ЛАЧХ разомкнутых систем с оптимальными регу ляторами положения.
Передаточная функция оптимального регулятора по ложения может быть определена из соотношения
kp,.n*W-p.a (p )W 0 ( p ) = J ( p ) ,
142
Подставляя в него / (р) и tѴ0(р) из (4-51) и (4-3), находим:
2 .С ІѴ + 1 . |
^ p + l _______ |
W V* \ P) — 3)41т^ д |
0,294x^2+ 0,77т/> + 1 |
Выражение (4-52) сложно для реализации. С целью упрощения примем 47^ = 0,77 т. Тогда вторым сомножи телем в (4-5'2) можно пренебречь и передаточная функ ция регулятора положения принимает вид:
^р.п(р): |
1 3 ,6 7 ^ + 1 |
|
(4-53) |
||
где |
ТпР |
|
3,41т2&д>П;і.£а |
||
7 Н= |
||
|
Погрешность при переходе от (4-52) к (4-53), как следует из сравнения кривых 1 и 2 на рис. 4-4, незначи-
Рис. 4-4. Переходные процес сы отработки сигнала W(t)
в системе рис. 4-1.
l ~ W p n (p) |
определена по |
(4-52); |
|
2 . - Г р;п( р ) - п о |
(4-53); |
3 - |
|
\Гр п (р )— по |
(4-66); |
4 — сигнал |
ЩІ).
тельна. Согласно (4-53) оптимальный регулятор является интегрально-пропорциональным и по условию — цифро вым, т. е. Wa(p) = 0 и
Г р . п (р ) = И7д (р ) = 13, 6 T J T B+ 1 /7 > |
(4 -5 4 ) |
Чтобы избежать необходимости осуществления опе раций умножения в цифровой форме, принимаем соглас но (2-4):
W ^ ( p ) ^ l 3 ß T J T n; Г Д1(р ) = 1;
Ц7ф2( р ) = 1 / Г и; Г д (р ) = 1/р.
143
При этом очевидно, что и?*ді (г) = 1 . Определение
как отмечалось в § 2-1 и 2-7, должно произво диться из условия обеспечения эквивалентности переда
точных функций \Ѵяг{р) и \V*RZ(z). |
Здесь возможно не |
||
сколько подходов, |
позволяющих |
получить |
различную |
степень совпадения |
свойств WRZ(p) |
и W*RZ(z). Рассмот |
|
рим некоторые из них. |
|
|
|
Представим етР в виде |
|
|
|
|
етР ^ \ + Т р . |
|
(4-55) |
Тогда, учитывая, что z = e TP, получаем: |
|
||
Up~TJ{z—l ) = T z - 4 ( l —z- *) =W*v( z ) . |
(4-56) |
Выражение (4-56) для передаточной функции цифро вого интегратора пояснялось в примере 2-1. Реализация его требует запоминания предыдущих значений сигналов рассогласования и выходного цифрового интегратора. Более простой вариант построения цифрового интегра тора, требующий запоминания только предыдущего зна чения одного выходного сигнала, характеризуется пере даточной функцией (1-5) и получил широкое использова ние на практике. Для лучшего приближения свойств цифрового и аналогового интеграторов, используя ап проксимацию Пада [Л. 38], можно представить етР в виде
Тр _ |
1 |
+ |
|
Тр/2 |
|
(4-57) |
|
~ |
1 |
- |
|
Тр/2 |
• |
||
|
|
||||||
При этом |
|
|
|
|
|
|
|
1 + Z - * |
Г*дз(2 ). |
(4-58) |
|||||
1 |
— Z- |
1 |
|||||
|
|
|
|||||
Для реализации W*j$(z) |
|
по |
(4-58) |
в сравнении |
с (4-56) ЦВУ должно содержать дополнительный сумма тор. Возможны и более сложные цифровые интеграторы [Л. 1, 2], т. е. лучшее приближение свойств цифрового и аналогового интеграторов достигается путем усложнения ЦВУ. Выбор того или иного способа реализации цифро вого интегратора лучше производить на этапе уточнен ного анализа из компромисса между сложностью по строения ЦВУ и величиной эффекта, достигаемого этим усложнением. Особенно удобен при этом метод матема тического моделирования на цифровых и цифро-аналого вых вычислительных машинах.
144
Рассмотренный в данном параграфе синтез регулято ра положения справедлив при условии Ф (/?) = ] и \Ѵэ(р) = \. Для других выражении Ф(р) и ином выборе передаточной функции эталонной системы определение G(p) производится аналогично по (3-34), (4-32) и (4-33)
4-4. СИНТЕЗ ЦИФРО-АНАЛОГОВОГО РЕГУЛЯТОРА ПОЛОЖЕНИЯ
В соответствии с изложенным в § 4-2 синтез опти • мальной передаточной функции цифро-аналогового регу лятора должен производиться на основе (3-34), (4-34) и (4-35). Примем и в этом случае передаточную функцию эталонной системы Wa(p) — 1, выбрав в качестве эталон ной идеальную по отношению к полезному сигналу си
стему. Тогда из сравнения (4-34) и |
(4-35) следует, что |
|
R( p) =K{ p) и |
|
|
Г г ( п ) - |
[Я (/>)/*" (1)] + |
1 . |
w |
К+ (р) |
|
Отсюда можно заключить, что наличие помех кванто вания на входе цифрового регулятора не накладывает ограничений на структуру системы управления. Дейст вительно, передаточная функция оптимальной системы может быть выбрана равной передаточной функции эта лонной системы. При этом ограничение М{і2(^)] обеспе чивается соответствующим выбором Wa(p). Выбор струк туры в этом случае должен осуществляться на основе детерминистского подхода, изложенного в гл. 3. Ход ре шения при этом включает в себя: 1 ) определение прин
ципиального решения для G(p)\ 2) расчет параметров G(p), обеспечивающих допустимое в соответствии с (4-2) значение M[t2(/)].
Перейдем к решению этой задачи. Пусть передаточ ная функция замкнутого контура регулирования скоро сти, как и прежде, определяется (4-38). Зададим соглас
но § |
3-4 в (3-31) и (3-32) гі= |
1, 7 * 2 — г4~ 0, а / * 3 будем счи |
тать |
неопределенным. После |
того, как будет найдена |
G(p), 7 * 3 определим из условия Mi[i2(^)] = i2ÄOn. Ограниче
ние М[і2(^)] обеспечивается ограничением полосы пропу скания замкнутой системы. При этом энергия, расходуе мая в силовой части электропривода, при каждом значе нии полезного сигнала также будет определяться полосой пропускания. Таким образом, в (3-72) величина Л3
косвенно задается из условия ограничения флюктуаций
10—181 |
145 |
тока якоря, вызванных действием гіомек квантования. Значения А3, соответствующие различным параметрам полезного сигнала, образуют некоторую область значе ний, которую мы назвали допустимой. При этом опти мальноспроектированный электропривод при любом ко личестве расходуемой энергии из допустимой области значений обеспечивает наилучшее приближение процесса отработки к процессу в эталонной системе.
Найдем принципиальное решение для G(p). Подстав
ляем |
Ь а {р), Н ( р ) |
и Ixx(p) = h v w ( p ) |
из (3-1), |
(3-12) |
и |
|||||||||
(4-36) |
в (3-31). После |
|
преобразований |
получаем: |
|
|||||||||
|
|
К ( Р ) : |
|
|
|
4 |
+ |
р 1 |
|
|
|
(4-59) |
||
|
|
|
9% |
Р2 ( - Р ) 2 ' |
|
|
||||||||
|
|
|
|
|
|
|
|
|
||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4b* = |
U l Pk l f r 3Tl . |
|
|
|
(4-60) |
||||||
Раскладываем К (р) |
|
на сомножители: |
|
|
|
|||||||||
|
К' (Р): |
2гзт1 p 2 — 2bp + |
2b* |
|
(4-61) |
|||||||||
|
|
^Р |
|
(- РУ |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
К |
+ (р) = |
р- + |
2Ьр + |
2Ь2 |
|
|
(4-62) |
|||||
|
|
|
|
|
|
|
Р2 |
|
|
|
|
|
|
|
Подставляя |
в |
(3-32) |
Wa(р) = |
1 |
и |
|
г2 = |
гі = |
0, полу |
|||||
чаем: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Я ( р ) " |
|
2U\ |
|
|
|
|
(4-63) |
||||
|
|
|
’ р 2 (— р)2 |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
Записываем |
выражение для |
R(p)/K~(p) |
и, разлагая |
|||||||||||
его на простые дроби, определяем [R(p) /К~ (р)]+: |
|
|||||||||||||
|
|
Г R(P) ) |
— |
Zbp + 2b* |
|
|
|
/4б4ч |
||||||
|
|
L ^ - .(P )J+ |
|
Р2 |
|
‘ |
|
|
( |
’ |
||||
Подставляем (4-64) и (4-62) в (3-34). При этом иско |
||||||||||||||
мое принципиальное решение для G(p) |
приобретает вид: |
|||||||||||||
|
|
|
|
= |
|
|
|
|
|
|
|
(4‘ 65) |
где т = 1 /2 Ь.
В (4-65) т является функцией г3, как это следует из (4-60). Выбор г3, а следовательно, и т должен произво диться из условия ограничения флюктуаций тока якоря согласно (4-2). Выразим M[i2{t)] через найденное реше ние для G(p) с помощью (4-37). Для этого сначала най-
146
дем выражение для W4(p)— WR(p)W<b(p). Передаточная функция разомкнутой оптимальной системы при G(p) в виде (4-65) определяется (3-53). Тогда передаточная функция регулятора может быть найдена из (3-82) при подстановке Ф(р) и lFKC(p) из (3-53) и (4-38).
^р.п(р) |
(~*Р4- 1)(47’[J, Р + |
1) |
(4-66) |
.цДа |
|
||
|
|
|
Согласно (4-66) оптимальный регулятор положения является ПИД регулятором. В соответствии с изложен ным в § 2 - 1 в цифровой форме должна вычисляться ин
тегральная составляющая, в связи с чем принимаем:
Wn( p ) =HTup, |
(4-67) |
где
и—‘ k . С *
Выбирая, как и в § 4-3, Ф(р) = 1, подставляем Ф(р), Ц7ц(р), W0(p) и Н(р) из (4-67), (4-8), (4-9) и (4-38) в (4-13). После вычисления и подстановки F(p) в (4-37), находим:
|
С О |
(«■й2 |
|
|
|
Ж [г ( 0 1 |
И |
|
diu, |
(4-68) |
|
|
А, («о) h3(— іа) |
|
|
||
где |
|
|
|
|
|
К (Н = 8 хТ (т у + |
(8Г Т + 2 г ) ( ь у + |
(47; + |
2т) іш + |
І6 *д.п. W
Пользуясь для вычисления этого интеграла формула ми табл. 3-1, получаем:
А(2Т^ + х)
(4-69)
тз (в^ + ^Г^ + т*)
Искомое значение т, соответствующее допустимому
уровню флюктуаций тока якоря, определяется из уравне
ния M\P(t)\ = '? , которое в соответствии с (4-69) при
нимает вид: |
|
•‘в + 47> 4 + 8 Г Ѵ - г |
-ф— 27; = 0. (4-70) |
доп |
1доп |
10* |
147 |
Согласно правила Декарта уравнение (4-70) имеет один положительный корень, значение которого и опре деляет искомую величину т. Так как это уравнение не может быть разрешено в общем виде, для определения т необходимо использовать приближенные методы.
Пример 4-2. Определим значение -с, соответствующее ограниче
нию |
М [іг (/)] значением ідОП= 0,01 при |
7^ = 0,01 |
с, £д.п* = 0,9 и |
||||
исходных данных примера |
4-1. |
|
|
|
|
|
|
|
Определяем величину |
А; |
|
|
|
|
|
|
3,14 (0,06)=-4-10-» |
|
|
|
|
||
|
А — 16 - (0,9)а - (0.01)2 - (0,075)= |
0,01 - (78,5)2 |
^ |
І0~ |
|
||
|
Подставляем значения |
А, ідОП, 7" |
в (4-70) и |
приводим |
послед |
||
нее |
к виду: |
|
|
|
|
|
|
|
Т5+0,04т4+0,0008Тз— 10-4—2- 10-в= |
0. |
(4-71) |
||||
|
Решая (4-71) приближенно, находим |
искомое |
значение |
т, рав |
ное 0,095 с. Подставляя т=0,095 с в (3-53), определяем передаточ ную функцию разомкнутой оптимальной системы:
7 (р) = |
0,19/7+1 |
|
|
0,018/Р |
1 |
||
|
Логарифмическая амплитудная частотная характеристика (ЛАЧХ) L(co) =20 lg |7(im) | изображена на рис. 4-3 (кривая 2).
Сравнение частотных характеристик оптимальных систем, представ ленных на рис. 4-3, показывает, что использование цифро-аналогово го регулятора позволяет при той же заданной степени влияния по. мех квантования обеспечить существенно более широкую полосу пропускания. При этом быстродействие системы с цифро-аналоговым регулятором по сравнению с системой, имеющей цифровой регуля тор, возрастает примерно в 2 раза.
Согласно (4-66) регулятор положения является ПИД
регулятором. При этом в соответствии с |
изложенным |
в § 2 - 1 интегральная составляющая закона |
регулирова |
ния, определяющая статическую точность, должна вы числяться в цифровой форме, а пропорциональная и диф ференциальная— в аналоговой. Целесообразно принять
ѴРф(р) = \ / Т а, г д е Ги определяется (4-67), а WR(p) = \/p.
Определение W*R(z) в соответствии с отмеченным в §2-1 и 2-7 должно производиться из условия обеспечения эквивалентности передаточных функций WR(p) и Ц7*д(г). Некоторые способы выбора Т57*д(;г) рассмотрены в § 4-3. Детальное исследование вариантов построения цифровых интеграторов из условия эквивалентности их свойств не прерывным интеграторам дается в {Л. 1 , 2 ].
148
Так как (4-66) для Wixri(p) структурно не отличается от Wp,c (p) (4-56), то схема аналогового ПД регулятора положения в принципе не отличается от схемы регуля тора скорости. Разница будет заключаться лишь в их параметрах. Поэтому все изложенное по этому вопросу в § 3-3 будет справедливо и здесь. Схемная реализация аналоговых регуляторов детально изложена в [Л. 10, 11, 16, 17].
На рис. 4-4 представлен переходный процесс отра ботки полезного сигнала в системе по рис. 4-1 с опти мальным регулятором (4-66). Сравнение кривых пере ходных процессов, соответствующих оптимально спроекти рованным цифровому и цифро-аналоговому регуляторам (кривые 1 и 3), наглядно иллюстрирует преимущест ва последнего. Необходимо, однако, отметить, что ис пользованный при синтезе цифро-аналоговых регулято ров подход предполагает, что уровень помех аналоговых датчиков неизмеримо меньше, чем помех квантования. При большом уровне помех аналоговых датчиков необхо димо будет сужать полосу пропускания системы и тогда соотношение кривых переходных процессов на рис. 4-4 не будет столь существенным.
Значение -г промышленных позиционных приводов не удается получить менее (0,04-4-0,06) с. Поэтому если при расчете по (4-70), где х определяется с учетом влияния только помех квантования, значение его будет меньше, то это лишь будет говорить о практически несуществен ном влиянии помех квантования и необходимости учета помех аналоговых датчиков. Подход к проектированию при учете последних, изложенный в § 3-5, может быть использован и здесь.
4-5. ОЦЕНКА ВЛИЯНИЯ КВАНТОВАНИЯ ПО ВРЕМЕНИ
■При выборе типа дискретных элементов для построения узлов АЦП и ЦВУ, исходя из допустимой степени влияния квантования по времени, необходимо знать допустимое значение периода преры вания Т. Для определения его требуется осуществить проверку вы
полнения неравенства (4-3). Такая проверка может потребоваться и для подтверждения корректности расчета, выполненного согласно методике, изложенной в § 4-2. Проверку (4-3) на основе определе ния сос.т, данного в § 4-1, осуществить трудно. Это связано с тем, что выбор степени приближения нулю амплитудной частотной ха рактеристики, определяющей сос.т, достаточно произволен, так как лишь косвенно оценивает влияние квантования по времени на про цессы в системе. Целесообразно поэтому для ряда типовых структур оптимальных систем, в том числе характеризующихся (4-47) и
149