CC BY
16ДИСКРЕТНАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ЦИФРОВОГО СЛЕДЯЩЕГО ЭЛЕКТРОПРИВОДА С СИНХРОННЫМ ИСПОЛНИТЕЛЬНЫМ ДВИГАТЕЛЕМ
С.Л. Лисин, А.В. Стариков
Самарский государственный технический университет 443100, г. Самара, ул. Молодогвардейская, 244
Разработана математическая модель следящего электропривода с синхронным исполнительным двигателем, учитывающая процесс квантования по времени. Найдены дискретные передаточные функции замкнутого электропривода для двух значений периода дискретизации. Показано, что уменьшение периода дискретизации значительно увеличивает быстродействие электропривода.
Ключевые слова: следящий электропривод, синхронный электродвигатель, дискретная передаточная функция.
Все современные следящие электроприводы имеют микропроцессорную техническую реализацию. Очевидно, цифровой способ обработки информации приводит к тому, что на работу систем управления таких электроприводов оказывает влияние процесс квантования сигналов по времени.
Целью данной статьи является исследование особенностей цифрового следящего электропривода с синхронным исполнительным двигателем, структурная схема непрерывного прототипа которого имеет вид, приведенный на рис. 1.
Рис. 1. Структурная схема непрерывного прототипа следящего электропривода с синхронным исполнительным двигателем
Электропривод имеет три контура: внутренний контур скорости и два контура положения. Для организации обратной связи по скорости сигнал безынерционного датчика положения с коэффициентом передачи kдn дифференцируется звеном с передаточной функцией
Ксс (p) = ,
где kосс - коэффициент передачи обратной связи по скорости.
Регулятор внутреннего контура (контура скорости) выбран пропорционально-
Сергей Леонидович Лисин, начальник отдела кадров.
Александр Владимирович Стариков (к.т.н., доц.), доцент кафедры «Электропривод и промышленная автоматика».
дифференциальным с передаточной функцией
W* (p) = knd (Tndp+1),
где knd - коэффициент передачи; Tnd - постоянная времени регулятора.
Регулятор второго контура представляет собой пропорциональное звено с коэффициентом передачи kn, а регулятор третьего (внешнего) контура принят интегральным с передаточной функцией
Wu (p) =^,
Tup
где Tu - постоянная времени интегрирования.
Синхронная машина, работающая в режиме бесколлекторного двигателя постоянного тока, на структурной схеме представлена колебательным звеном [1]
Woy1( p) = -
k„
(тк2 p2 + 2^1^ +1) р’
где koy - коэффициент передачи объекта; Tк - постоянная времени колебательного
звена; - его коэффициент демпфирования.
Силовой преобразователь принят апериодическим звеном с передаточной функцией
W ( p) = —^— ,
сп \г / гг 7 ?
TспP + 1
где kсп и ^ - коэффициент передачи и постоянная времени силового преобразователя соответственно.
Для разработки дискретной математической модели рассматриваемого цифрового следящего электропривода перейдем к z-преобразованиям. Тогда его структурная схема примет вид, приведенный на рис. 2.
<pXz)
(->
С-)
(->
w*(z) +
w, (z)
<Pde(Z)
к-
'дп
Рис. 2. Структурная схема следящего электропривода с синхронным исполнительным двигателем при переходе к z-преобразованиям
Непрерывный объект управления с учетом экстраполятора нулевого порядка представлен дискретной передаточной функцией W0(z). Поскольку функцию экстраполятора в рассматриваемом электроприводе выполняет силовой преобразователь, то передаточная функция непрерывной части системы имеет вид
Кпку
(TKp + 2$Xp+1) p
Дискретная передаточная функция такого объекта с учетом экстраполятора ну-
левого порядка находится по известным правилам [2] и в конечном итоге запишется следующим образом:
Ко( г) =
Рдв (г)
= к к
аг2 + Ьг + с
^пд(г) сп °У (г -1)(г2 - 2zdcos[Т + й2)
(1)
1 - 2<52
где а = Т - 2%кТк (1 - й С08 [Т)-------------------------— й 8ІП [Т ;
Ь = 2
[
1 - 2Е1
%кТк (1 - й2) +------^ й 8Іп [Т - Тй со 8 [Т
[
с = Тй2 + 2^кТк (й2 - й со8 [Т) -
1 - 2£ [
й 8Іп [Т;
Т - период дискретизации (квантования по времени).
В состав цифрового управляющего устройства рассматриваемого следящего электропривода с синхронным исполнительным двигателем входят пропорционально-дифференциальный, пропорциональный и интегральный регуляторы. Дискретная передаточная функция пропорционально-дифференциального регулятора при вычислении производной как первой обратной разности будет иметь вид
( Т \
кпд (Тпд + Т)' г -
Кд (г) =-
Тпд + Т ,
(2)
Тг
Передаточная функция регулятора второго контура представляет собой пропорциональное звено с коэффициентом передачи
К (г) = кп. (3)
Дискретная передаточная функция интегрального регулятора внешнего контура, осуществляющего вычисление интеграла как полной суммы, будет выглядеть следующим образом:
К (г) = Тг . (4)
Ти (г -1) ' '
Так как дифференцирующее звено также использует в своем алгоритме работы первую обратную разность, его передаточная функция будет иметь вид
к°сс (г - 1)
К°сс (г) =
Тг
(5)
Поскольку правила преобразования структурных схем цифровых систем управления аналогичны правилам преобразования структурных схем непрерывных систем, формулы (1) - (5) позволяют найти дискретную передаточную функцию замкнутого следящего электропривода:
Кз(г)=
Ь05 г5 + Ь15 г4 + Ь25 г3 + Ь35 г2
г6 + а15 г5 + а25 г4 + а35 г3 + а45 г2 + а55 г + а6
(6)
, к2 а(Тпд + Т)
где Ьо5 = ------; Ь,5 =
15
-дпТи
-дпТи
КЬ - а)Тпд +ЬТ ];
ь25 =-[(Ь - с)Тпд - сТ]; Ь35 =- ;
кдпТи кдпТи
2 + 21008рт - №+ к.Т)а<Тпд + Т) - к2а(Т,д + Т)
а25 = 1 + ё2 + 41 008 РТ +
та Ти
к1 [ф - 3а)Тпд +(Ь - 2а)Т ] +
Т2
+ к2 [(Ь - 2а)Тпд + (Ь - а)Т] + к2 [(Ь - а)Тпд + ЬТ] ’
Т Ти
Г 2 к [ (3Ь - 3а - с)Тпд + (2Ь - а - с)Т]
а35 = -Г 212 + 21 008 рТ + -^--------------------------------—---------------------— ■
, к2 [(2Ь - а - с)Тпд +(Ь - с)Т] , к2 [(Ь - с)Тпд - сТ]Г
а 12 , к1 [(3Ь - а - 3с)Тпд + (Ь - 2с)Т], к2 [(Ь - 2с)Тпд - сТ] ^сТпд .
а45 1 + 2 + .
^ т 2 т т
_[ к1 [(Ь - 3с)Тпд - сТ ] - к2 сТпд 1 . а =-к1сТпд а55 — I о 1 ; аб5
Т 2 т [ ? 65 Т 2
Дискретная передаточная функция (6) позволяет исследовать устойчивость и реакцию рассматриваемого цифрового следящего электропривода с синхронным исполнительным двигателем на изменение управляющего воздействия. Выражение (6) имеет характеристический полином шестого порядка, поэтому применять алгебраические критерии устойчивости нецелесообразно. Граничные значения параметров будем искать численными методами путем оценки модуля корней уравнения
2 + аи 2 + а^52 + а35 z + а45 z + а55 z + а^ _ 0.
С позиций параметрического синтеза регуляторов интересно определение граничного, с точки зрения устойчивости, значения Тр периода дискретизации. Применительно к следящему электроприводу с синхронным электродвигателем 5FK70605AF71 определена величина Тр для двух наборов параметров регуляторов (см. таблицу).
Зависимость граничного значения Тгр периода дискретизации от величин Т , кпд, Тпд, к , Т и к
сп 5 пд ^ пд^ п1 и осс
и
Т , с сп кпд Тпд, с кп Ти, с косс , с Тр , с
0,0002 256 0,0045 32 0,001 0,0256 0,00023
0,0016 2 0,1175 3 0,0128 0,0128 0,0022
Данные таблицы показывают, что рассматриваемый цифровой электропривод всегда устойчив, если при выборе параметров регуляторов за постоянную времени Тсп силового преобразователя принимать величину периода дискретизации Т .
Промоделируем разработанный цифровой следящий электропривод в про-
граммной среде Ма1кЬ Simulink для двух значений периодов дискретизации Т = 0,0016 с и Т = 0,0002 с и соответствующих наборов параметров регуляторов.
При Т = 0,0016 с дискретная передаточная функция замкнутого электропривода принимает следующие численные значения:
„г/. 0,003894z5 + 0,0111234 - 0,011163г3 - 0,003552г2
Щ ( z) = —------------5-----------4----------3-----------2----------------------. (7)
z6 - 3,712493г5 + 5,656237z4 - 4,35954z3 + 1,5500462 - 0,058714z - 0,075774
Передаточная функция (7) позволяет исследовать реакцию разрабатываемого электропривода на изменение управляющего воздействия. График переходного процесса (рис. 3) показывает, что при Т = 0,0016 с время входа в 2%-ю зону составляет 1пп = 0,0368 с, а перерегулирование - а = 0,946 %.
From: Input Pointl То: Output Pointl
Time (sec)
Рис. 3. Переходный процесс в цифровом следящем электроприводе по управляющему воздействию при периоде дискретизации T = 0,0016 с
Для периода дискретизации T = 0,0002 с передаточная функция следящего электропривода становится равной
„г/. 0,002716z5 + 0,008211z4 - 0,007662z3 - 0,002575z2
( z) = 6----------5------------4-----------3------------2-----------------------. (8)
z6 - 3,841663z5 + 6,115285z4 - 4,9625z3 + 1,867273z2 - 0,074699z - 0,103012
На основании формулы (8) построен график переходного процесса (рис. 4). Время переходного процесса в этом случае составляет tnn = 0,0048 с, а перерегулирование - а = 0,514%.
Сравнение полученных результатов показывает, что уменьшение величины периода дискретизации при выбранном структурном построении цифрового следящего электропривода с синхронным исполнительным двигателем значительно увеличивает его быстродействие. Таким образом, разработанная дискретная математическая модель подтверждает эффективность применения способа построения следящего электропривода с минимальными вычислительными затратами.
0 1 2 3 5 6r с
Time (sec) ~.
A *T
Рис. 4. Переходный процесс в цифровом следящем электроприводе по управляющему воздействию при периоде дискретизации T = 0,0002 с
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Лысов М.С., Стариков А.В., Стариков В.А. Линеаризованная математическая модель синхронного электродвигателя при различных способах управления его скоростью // Вестник Самарского государственного технического университета. Серия: Физико-математические науки, 2008. - № 1 (16). -Самара: СамГТУ, 2008. - С. 102-107.
2. Микропроцессорные системы автоматического управления / В.А. Бесекерский, Н.Б. Ефимов, С.И. Зиатдинов и др.; под общ. ред. В.А. Бесекерского. - Л.: Машиностроение, 1988. - 365 с.
Статья поступила в редакцию 29 ноября 2012 г.
DISCRETE MATHEMATICAL MODEL OF THE DIGITAL OF THE SERVO DRIVE WITH THE SYNCHRONOUS MOTOR
S.L. Lisin, A. V. Starikov
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100
The mathematical model of the servo drive with the synchronous motor, considering quantization process on time is developed. Discrete transfer functions of the closed electric drive for two values of the sampling time are found. It is shown that reduction of the sampling time considerably decreases response time of the electric drive.
Keywords: the servo drive, the synchronous motor, the discrete transfer function.
Sergey L. Lisin, Head of HR.
Alexander V. Starikov (Ph.D. (Techn.)), Associate Professor.
