close

Вход

Забыли?

вход по аккаунту

?

Параметрический синтез цифровых систем управления с широтно-импульсными преобразователями и эталонными моделями.

код для вставкиСкачать
ПАРАМЕТРИЧЕСКИЙ СИНТЕЗ ЦИФРОВОЙ СИСТЕМЫ
УПРАВЛЕНИЯ СКОРОСТЬЮ ДВИГАТЕЛЯ ПОСТОЯННОГО
ТОКА С ЭТАЛОННОЙ МОДЕЛЬЮ
А.В. Гурьянов, В. В. Кротенко
Сформулирована методика синтеза цифровой системы управления скоростью двигателя постоянного тока
с эталонной моделью при заданных периоде дискретности вычисления управления и частоте коммутации
силовых ключей широтно-импульсного преобразователя.
Наиболее перспективными и практически применимыми на современном этапе
являются двухэтапные методики синтеза цифровых систем управления с широтноимпульс-ными преобразователями (ШИП), ориентированные на ПЭВМ, при которых
предварительный синтез осуществляется на основе приближенных эквивалентных моделей, а окончательный анализ с уточнением результатов предварительного синтеза
производится на ПЭВМ с использованием имитационных моделей, достаточно точно
отражающих как электромагнитные процессы в энергетической подсистеме с ШИП,
так и дискретный характер сигналов в информационной подсистеме. Именно такая методика используется при синтезе рассматриваемой системы.
Прецизионные электроприводы строятся по широко распространенной схеме
подчиненного регулирования координат на базе моментных исполнительных двигателей постоянного тока, питаемых от транзисторных ШИП. Как правило, позиционный
электропривод содержит контуры регулирования тока, скорости и положения [1].
Быстродействие позиционного электропривода по управляющему воздействию в
случае настройки его контура регулирования скорости на «оптимум по модулю» практически в два раза выше быстродействия электропривода, контур регулирования скорости которого построен по классической схеме с настройкой на «симметричный оптимум» [2, 3], однако первый при компенсации возмущений обладает, как правило, неудовлетворительными по качеству переходными процессами. Сочетать показатели настроенных на «симметричный оптимум» систем при отработке ими возмущений с высоким быстродействием по управлению систем, настроенных на «оптимум по модулю»,
позволяют системы, построенные по схеме с эталонными моделями [4].
Задачей синтеза будем считать определение при заданных частоте коммутации
силовых ключей ШИП Тк и периоде дискретности вычисления управления То параметров цифровых регуляторов, обеспечивающих в цифровом контуре регулирования скорости с эталонной моделью процессы, по характеру близкие к стандартным процессам
в ее эквивалентной линейной модели.
Структурная схема системы изображена на рис. 1. Система содержит исполнительный ДПТ, силовой каскад СК, инерционный датчик скорости с коэффициентом передачи Кдс и постоянной времени Тдс, безынерционный датчик тока с коэффициентом
передачи Кдт, реализованные на микроЭВМ эталонную модель ЦЭМ, цифровые регуляторы скорости ЦР1, ЦР2 и тока ЦР3, а также широтно-импульсный модулятор ШИМ.
В представленной системе можно выделить основной и дополнительный контуры
регулирования скорости и контур регулирования тока.
Основной контур регулирования системы включает ДПТ, усилитель мощности и
цифровой регулятор ЦР1, с помощью которого формируются заданные динамические
показатели системы по управляющему воздействию.
Дополнительный контур регулирования содержит эталонную модель ЦЭМ и цифровой регулятор ЦР2. В случае применения эталонной модели с динамическими характеристиками, идентичными характеристикам осинового контура регулирования, в системе сохраняются заданные динамические качества по управляющему воздействию,
139
при этом заданные динамические качества по возмущающему воздействию формируются независимо с помощью цифрового регулятора ЦР2.
ωзад(n, zс)
ЦЭМ
ε2
ωэм
ε3
ЦР2
ε1
Uцр1
ЦР1
ωос(n, zс)
ЦР3
ШИМ
γ(N,zцап)
iос(n,zт)
Кдс
СК
Кдт
-Мс/См
Тдс p+1
ω
Uцр2
Rя
СеТм p
i
1/Rя
Тя p+1
Се
Рис. 2. Структурная схема цифровой системы регулирования скорости ДПТ
с эталонной моделью
Синтез цифрового контура регулирования тока рассмотрен в [5], цифровой ПИрегулятор тока, реализуемый на основании выражений
ε1(n) = Uцр1 (n) + Uцр2 (n) - iос(n),
(1а)
Uтп (n) = ε1(n) Ктп,
(1б)
Uти1 (n) = Uти1 (n-1) + Кти ε1 (n-1),
(1в)
U цр3 (n) = Uтп (n) + Uти1 (n),
(1г)
где Uцр1(n), Uцр2(n) – сигналы задания на входе контура регулирования тока; Uцр3(n) –
выходной сигнал регулятора тока; Ктп и Кти – коэффициенты передачи пропорционального и интегрального каналов регулятора тока, определяемые соответственно на основании выражений
(2а)
Ктп = Rя (1- exp(-То / Тт))/ Кдт Кст (1- exp(-То / Тя)),
Кти = Rя (1- exp(-То / Тт))/ Кдт Кст,
(2б)
обеспечивает при отработке скачков задания переходные процессы, по характеру близкие к экспоненциальному процессу с постоянной времени Тт ≥ То. На рис. 2 представлены предельные по быстродействию переходные процессы в контуре при
Тк = Тт = 0.001 с.
Оптимизация основного контура системы осуществляется с использованием эквивалентной модели со структурной схемой, изображенной на рис. 3. Настройку контура на «оптимум по модулю» обеспечивает регулятор П-типа [6], при этом в модели
контур регулирования тока представлен звеном в виде апериодического звена первого
порядка с постоянной времени Тт и коэффициентом передачи 1/Кдт.
140
Рис. 2. Переходные процессы в контуре регулирования тока
ωзад
Кпо
ωосэ
1/Кдт
iэ
Тт p + 1
Rя
ωэ
СеТм p
Кдс
(Тзпс p + 1)(Тдс p + 1)
Рис. 3. Структурная схема ЭЛМ, используемая при оптимизации системы
по управляющему воздействию
Малая некомпенсированная постоянная основного контура по управляющему воздействию Тμu определяется выражением
Тμu = Тт + Тдс + Тзпс,
(3)
где Тзпс =То/2 – постоянная времени, учитывающая запаздывание, вносимое ЦР1 П-типа.
Выбор коэффициента передачи Кпо на основании соотношения
Кпо = Тм Се Кдт / (2 Тμu Rя Кдс)
(4)
обеспечивает указанную настройку, передаточная функция основного контура эквивалентной модели характеризуется выражением
ω ( p)
1
W ( p ) = ОСЭ
=
.
(5)
2 2
ω ЗАД ( p) 2Tμu p + 2Tμu p + 1
Коэффициент передачи Кп1 цифрового П-регулятора (ЦР1) определяется на основании выражения
Кп1 = Кпо = Тм Се Кдт / (2 Тμu Rя Кдс).
(6)
На рис. 4 представлены переходные процессы в основном контуре цифровой системы регулирования скорости с ШИП и цифровой эталонной модели (соответственно
кривые ω, ωэ) при отработке скачков задающего и возмущающего воздействий. При
отработке задающего воздействия переходные процессы в основном контуре системы и
ЦЭМ совпадают, при отработке возмущающего воздействия основной контур обладает
статической ошибкой, в ЦЭМ такая ошибка отсутствует.
141
Рис. 4. Переходные процессы в основном контуре системы
Оптимизацию системы по возмущающему воздействию следует осуществлять с
использованием полной эквивалентной линейной модели [6], структурная схема которой изображена на рис. 5,а. Полная ЭЛМ содержит объект управления, регулятор в виде последовательно включенных ПИ-регулятора с коэффициентами передачи Кп и Ки и
двух ПД-регуляторов с временами упреждения Тд1 и Тд2, три апериодических звена
первого порядка с постоянными времени Тμ1 , Тμ2 , Тзпс, учитывающими неполную компенсацию ЦР2 постоянных Тт , Тдс и вносимое им запаздывание.
Мс/См
ωэ
Rя
iэ
Кдс
СеТм p
(Тзпс p + 1)(Тдс p + 1)
1/Кдт
(Кп p + Ки)(Тд1 p + 1)(Тд2 p + 1)
Тт p + 1
ωосэ
p(Тμ1 p + 1)(Тμ2 p+ 1)
а)
Мс/См
ωэ
Rя
Тμf p + 1
СеТм p
iэ
Кдс/Кдт
Кп p + Ки
p
б)
Рис. 5. Структурные схемы ЭЛМ, используемые при оптимизации системы
по возмущающему воздействию
Расчетная эквивалентная линейная модель представлена структурной схемой, изображенной на рис. 5,б, ее настройка на «симметричный оптимум» осуществляется в
случае выбора коэффициентов передачи Кп и Ки ПИ-регулятора на основании соотношений
142
(7)
Кп = Тм Се Кдт / (2 Тμf Rя Кдс) ,
Ки = Кп /Ти = Тм Се Кдт / (8 Тμf 2Rя Кдс) ,
(8)
передаточная функция эквивалентной модели по возмущающему воздействию характеризуется выражением
(4Tμ . f p + 1) / С М
i( p)
W ( p) =
=
,
3 3
М ЗАД ( p ) 8Tμ . f p + 8Tμ . f 2 p 2 + 4Tμ . f p + 1
где Тμf – малая некомпенсированная постоянная времени системы по возмущающему
воздействию.
В табл.1 представлены типы регуляторов, применяемых при оптимизации ЭЛМ
по возмущающему воздействию, и соответствующие им постоянные времени Тμf.
Тип АР
Компенсируемая
постоянная
И
-
ПИ-ПД
Тт
ПИ-ПД-ПД
Тт, Тдс
Параметры АР
Кп=0, Ки ≠ 0,
Тд1=0, Тμ1 = 0,
Тд2=0, Тμ2 = 0.
Кп ≠ 0, Ки ≠ 0,
Тд1=Тт, Тμ1= Т0/2,
Тд2=0, Тμ2 = 0.
Кп ≠ 0, Ки ≠ 0,
Тд1=Тт, Тμ1=Т0/2,
Тд2=Тдс, Тμ2=Т0/2.
Таблица 1
Малая некомпенсированная
постоянная Тμf ЭЛМ
Тμf = Тт+Тдс+Тзпс =
= Тт+Тдс+ Т0/2. (9)
Тμf = Тμ1+Тдс+Тзпс =
= Тдс+Т0. (10)
Тμf = Тμ1+Тμ2+Тзпс =
= 1.5Т0.. (11)
Если в режиме парирования возмущений не осуществляется компенсация постоянных Тт и Тдс, то в этом случае достаточно применения в качестве АР2 И-регулятора с
передаточной функцией W(p)=Ки/p и коэффициентом передачи Ки, определяемым на
основании (8). Коэффициент передачи Ки2 ЦР2 И-типа с учетом выражения (8) определяется на основании выражения
Ки2 = То Ки = То Тм Се Кдт / (8(Тт+Тдс+То/2)2Rя Кдс),
(12)
при этом малая некомпенсированная постоянная времени Тμf определится на основании
выражения (9).
На рис. 6 представлены переходные процессы в цифровой системе регулирования
скорости с ЦР2 И-типа. Ее быстродействие по управляющему воздействию в два раза
выше быстродействия рассмотренной в [2, 3] классической системы подчиненного регулирования скорости ДПТ, настроенной на «симметричный оптимум», по возмущающему воздействию быстродействие обеих систем одинаково.
Увеличить быстродействие системы в режиме парирования возмущений можно
при использовании в качестве ЦР2 регуляторов, осуществляющих компенсацию либо
одной постоянной Тт, либо двух постоянных Тт и Тдс. В случае компенсации одной постоянной Тт достаточно применения ПИ-ПД-регулятора, в случае компенсации двух постоянных Тт и Тдс, необходимо применение ПИ-ПД-ПД-регулятора.
Аналоговый ПИ-ПД-ПД-регулятор полной ЭЛМ может быть преобразован в эквивалентный регулятор с разделенными каналами на основании выражения
W(p) = (Кп+Ки/p)(Тд1p+1)( Тд2p+1) = Кэд2p2+Кэд1p+Кэп+Кэи/p,
(13)
где
Кэд2 = Кп Тд1 Тд2,
(14)
Кэд1 = Кп (Тд1 + Тд2) + Ки Тд1 Тд2,
(15)
Кэп = Кп +Ки (Тд1 + Тд2),
(16)
Кэи = Ки.
(17)
143
Рис. 6. Переходные процессы в системе с эталонной моделью и ЦР2 И-типа
Реализация ЦР1 и ЦР2 , в случае выполнения последнего по схеме с разделенными
каналами со структурой, соответствующей структуре АР (13), осуществляется на основании выражений
ε 2(n) = ωзад(n) – ωос(n),
(18а)
Uцр1 (n) = Kп1 ε 2(n),
(18б)
ε 3(n) = ωэм(n) – ωос(n),
(18в)
(18г)
Uп2 (n) = Kп2 ε 3(n),
Uи2 (n) = Uи2 (n-1) + Ки2 ε 3(n-1),
(18д)
(18е)
Uи2т(n) = Uи2 (n) +0.5 Ки2 ε 3(n),
(18ж)
Uд21 (n) = Kд21 ( ε 3(n) - ε 3(n-1)),
Uд22 (n) = Kд22 ( ε 3(n) - 2 ε 3(n-1)+ ε 3(n-1)),
(18з)
(18и)
U цр2(n) = Uп2(n) + Uи2т (n)+ Uд21 (n)+ Uд22 (n),
где Uцр1(n), Uцр2(n) – выходные сигналы ЦР1 и ЦР2; Uп2 (n) , Uи2т(n) , Uд21(n), Uд22 (n), Kп2,
Ки2, Kд21, Kд22 – сигналы П-, И- и Д-каналов ЦР2 и их коэффициенты передачи, соответственно; ωэм(n) – выходной сигнал ЦЭМ.
Расчет коэффициентов передачи ЦР2 Кп2, Ки2, Кд21 и Кд22 следует осуществлять с
учетом соотношений (14)–(17) на основании выражений
(19)
Кп2 = Кэп - Кп0 = Кп +Ки (Тд1 + Тд2) - Кп0,
Ки2 = Т0 Кэи = Т0 Ки,
(20)
(21)
Кд21 = Кэд1 / Т0 = (Кп (Тд1 + Тд2) + Ки Тд1 Тд2,)/Т0 ,
2
2
Кд22 = Кэд2 /Т0 = Кп Тд1 Тд2 /Т0 .
(22)
На рис. 7,а представлены переходные процессы в цифровой системе регулирования скорости с ЦР2 ПИ-ПД-типа, на рис.7,б – с ЦР2 ПИ-ПД-ПД-типа. Быстродействие
систем в режиме парирования возмущения характеризуется малыми некомпенсированными постоянными времени, определяемыми, соответственно, выражениями (10) и
(11). В последнем случае быстродействие максимально и может быть увеличено лишь
за счет уменьшения периода дискретности То. Рассмотренная система относится к классу беспоисковых самонастраивающихся систем с эталонной моделью, в ней при изменении параметров объекта управления обеспечиваются процессы, близкие по характеру
к процессам в эталонной модели, при этом ее адаптационные свойства улучшаются с
ростом коэффициента передача цифрового регулятора ЦР2 [4].
На рис. 7,в представлены переходные процессы в цифровой системе регулирования скорости с ЦР2 ПИ-ПД-ПД-типа при суммарном моменте инерции ДПТ с нагруз-
144
кой, трехкратно превышающем расчетный, период дискретности То в рассматриваемом
случае составлял величину То=0.0001 с.
а)
б)
в)
Рис. 7. Переходные процессы в системе регулирования скорости
с эталонной моделью с ЦР2 ПИ-ПД-типа (а) и ПИ-ПД-ПД-типа (б,в).
145
ЦЭМ настроенной на «оптимум по модулю» системы может быть построена на
основании структурной схемы, представленной на рис.8,а, ее расчет осуществляется на
основании выражений
ε 1(n) = xзад (n) - Uэм(n) ,
(23а)
ε 2(n) = ε 1(n) - Uос (n-1) ,
(23б)
Uи1(n) = Uи1(n-1) + Ки1 (n) ε 2(n),
(23в)
(23г)
Uос (n) = Uи1(n),
Uc (n) = (Uи1 (n-1) +Uи1 (n))/2,
(23д)
Uи2(n) = Uи2 (n-1) + Ки2Uc (n-1),
(23е)
(23ж)
Uи2т (n) = Uи2 (n) +0.5 Ки2Uc (n),
Uэм (n) = Uи2т (n),
(23з)
где xзад(n) – сигнал задания на входе ЦЭМ, ε 1(n) и ε 2(n) – сигналы рассогласования,
Uи1(n) и Uи2т(n) – выходные сигналы интеграторов, Uэм(n) – выходной сигнал ЦЭМ.
МикроЭВМ
xзад(n)
ε1(n)
ε2(n)
Uи1(n) Uи1(n)+Uи1(n-1)
Uc(n)
Uэм(n)
2
Uос(n-1)
а)
xзад
Киа
xэм
(Tνp + 1)p
xосэ
1
2Tзап p +1
б)
Рис. 8. Структурные схемы эталонной цифровой модели системы, настроенной на «оптимум по модулю» (а), и ее эквивалентной линейной модели (б)
ЦЭМ может быть поставлена в соответствие расчетная эквивалентная линейная
модель со структурной схемой, представленной на рис.8,б. Модель содержит интегрирующее звено с коэффициентом передачи Киа, и два апериодических звена первого порядка с единичными коэффициентами передачи и постоянными времени Тν и 2Тзап, последнее учитывает вносимое двумя интеграторами цифровой модели запаздывание.Оптимизация линейной эквивалентной модели осуществляется с использованием
расчетной эквивалентной модели, получаемой путем замены двух апериодических
звеньев с постоянными Тν и 2Тзап одним звеном с малой некомпенсированной постоянной Тμ, равной
Тμ = Тν + 2Тзап.
(24)
Выбор коэффициента передачи Киа интегрирующего звена на основании соотношения
Киа = 1 / 2 Тμ
(25)
обеспечивает в расчетной эквивалентной модели настройку на «оптимум по модулю».
В случае реализации цифровых интеграторов на основании выражений (23в) и (23,е,ж)
постоянная Тзап определяется на основании соотношения Тзап= То /2.
146
Выбор коэффициента передачи Ки1 осуществляется на основании выражения
Ки1 = 1- exp(-To / Tν),
(26)
в этом случае интегратор, охваченный отрицательной обратной связью, обладает динамическими характеристиками, близкими к характеристикам апериодического звена
первого порядка с постоянной времени Tν. Расчет коэффициента передачи Ки2 может
быть осуществлен с учетом (25) на основании выражений
Ки2 = То / 2 (Тν + То,).
(27)
На рис. 9 представлены переходные процессы в эталонной цифровой модели
системы xэм(n) и ее эквивалентной линейной модели xэм(t) при отработке скачка задающего воздействия.
Рис. 9. Графики переходных процессов в эталонной цифровой модели системы
с настройкой на «оптимум по модулю»
Параметры ДПТ и исходные данные для расчета системы представлены в
табл.2.
Таблица 2
Двигатель постоянного тока типа ДПМ-160-7-Д09
60
Максимальное напряжение питания Еп, В
6
Сопротивление обмотки якоря Rя, Ом
Номинальный ток Iном, А
4
Номинальный момент Мном, Нм
7
Электрическая постоянная цепи якоря Тя, мс
5
Механическая постоянная суммарная Тм, мс
50
0.56
Скорость холостого хода ωхх, рад/с
1.75
Постоянная по моменту, См, Нм/А
107.14
Постоянная по противоЭДС Се, Вс/рад
Транзисторный ШИП с управлением, осуществляющим изменение длительности
импульса дважды на периоде коммутации Тк
Период коммутации цепи якоря Тк, мс
0.4
Период дискретности То, мс
0.2
Стат. коэффициент передачи ШИП Кст,
6
Эквивалентная постоянная контура тока Тт, мс
1
Коэф. передачи датчика тока Кдт, В/А
1
Коэф. передачи датчика скорости Кдс, Вс/рад
17.857
Постоянная времени датчика скорости Тдс, мс
0.5
Приведенные результаты имитационного моделирования системы подтверждают
корректность предлагаемой методики синтеза цифровых систем управления с эталонными моделями и широтно-импульсными преобразователями.
147
Литература
1. Шрейнер Р.Т. Системы подчиненного регулирования электроприводов. Часть 1.
Электроприводы постоянного тока с подчиненным регулированием координат.
Екатеринбург, 1997.
2. Кротенко В.В. Параметрический синтез и анализ цифровой системы управления
прецизионного электропривода лидарной станции. Научно-технический вестник
СПбИТМО. Выпуск 11, С.86-96, СПб , 2003.
3. Кротенко В.В., Толмачев В.А., Томасов В.С., Синицин В.А. Синтез микропроцессорной системы управления электропривода опорно-поворотного устройства. // Изв.
вузов. Приборостроение. 2004. Т.47. №11. С.23–30.
4. Солодовников В.В., Шрамко Л.С. Расчет и проектирование аналитических самонастраивающихся систем с эталонными моделями. М.: Машиностроение, 1972.
5. Гурьянов В.А., Кротенко В.В. Цифровая система управления источника тока, построенного на основе транзисторного ШИП. / Научно-технический вестник
СПбИТМО. Выпуск 3. СПб , 2001. С.120–125.
6. Кротенко В.В. Параметрический синтез цифровых систем управления с транзисторными широтно-импульсными преобразователями. // Изв. вузов. Приборостроение.
2003. Т.46. №6. С.25–31.
148
1/--страниц
Пожаловаться на содержимое документа