УДК 621. 313. 323
А. Б. Красовский
ИМИТАЦИОННЫЕ МОДЕЛИ В ТЕОРИИ
ВЕНТИЛЬНО-ИНДУКТОРНОГО
ЭЛЕКТРОПРИВОДА
Обоснована эффективность использования имитационного моделирования для исследования вентильно-индукторного электропривода. Рассмотрены общие принципы построения имитационных моделей и предложен вариант их реализации в среде МЛТЬЛБ-БШиЬШК. На этой основе выработаны рекомендации по методам управления приводом, обеспечивающим получение максимальной мощности, снижение пульсаций момента и исключение аномальных режимов. Достоверность результатов подтверждена экспериментальными исследованиями.
Вентильно-индукторный электропривод (ВИП) является одной из наиболее бурно развивающихся областей современной электромеханики. О сновные преимущества ВИП по сравнению с традиционными электроприводами обусловлены простотой, надежностью и высокой технологичностью электрической машины (ВИМ), на базе которой он выполняется [1]. Однако, несмотря на признанные практически всеми специалистами достоинства, массового распространения он еще не получил. Причина этого состоит в том, что потенциально присущие ВИП достоинства реализуются только при соответствующих алгоритмах управления. При этом специфика ВИП предполагает использование отличных от традиционных для классических систем электропривода инструментов их разработки.
Наиболее характерные особенности ВИП состоят в следующем [2]:
— дискретность в управлении;
— питание фазных обмоток двигателя осуществляется однополяр-ными импульсами тока, форма которых в общем случае претерпевает значительные изменения с изменением режима работы привода;
— изменение в достаточно широких пределах состояния магнитной системы двигателя в цикле коммутации каждой фазы; линейный режим работы (без насыщения) чередуется с режимом локального насыщения взаимодействующих зубцов, которое при определенных условиях может переходить к частичному общему насыщению магнитной системы двигателя;
— двухсторонняя зубчатость магнитной системы двигателя при ограниченном числе его фаз приводит практически к трапецеидальному закону изменения магнитной проводимости воздушного зазора меду статором и ротором.
Поэтому наличие методов и средств исследования, хорошо приспособленных для решения задач анализа и синтеза и адекватно воспроизводящих реальные процессы в приводе, для ВИП имеет особое значение.
Принципы построения моделей ВИП. В основу математического описания процесса преобразования энергии в ВИП могут быть положены, как обычно, уравнения электрического равновесия фаз двигателя
Uk — 1ф k Rk +
и уравнение электромагнитного момента
dt
M —
^ dWk
д©
k=l
(1)
(2)
1ф k=const
где ик, 1ф к, Як, Фк — напряжение, ток, активное сопротивление и по-токосцепление к-й фазной обмотки ВИМ, к = 1,... ,т, т — число фаз
-<ф к
ВИМ; Wk — Ф^;dlф k — коэнергия k-й обмотки; © — угловой сдвиг
между взаимодействующими зубцами статора и ротора.
Перечисленные выше особенности ВИМ обусловливают нелинейную зависимость потокосцепления Фк как от 1ф к, так и от в (см. рис. 1, а). Поэтому в развернутой форме уравнение (1) имеет вид
ТТ _ _ , дfrk(1фk, ©) dlфk , д(1фk, ©) d©
Uk — 1ф k Rk +--^--т:--+
dIt
ф k
dt
д©
dt
(3)
Слагаемые в правой части уравнения (3) — ЭДС самоиндукции
дФk (1фk, ©) d^
д1ф k dt от этих же переменных.
и ЭДС движения
дФk (1фk, ©) d©
д©
dt
нелинейно зависят
Ф ф ф <1
Рис. 1. Характерный вид магнитных характеристик ВИМ:
1 -1ф 1; 2-1ф 2 ; 3-1ф 3; 4-1ф 4
Аналогично, в выражении (2), преобразованном к виду
-ф k k м=£ / 1 = £ / ^Мц^/ф k, (4)
к=1 о к=1 о
фазная индуктивность Lk таюке зависит от тока и положения. На рис. 1, б показана трансформация кривой интегральной индуктивности L(©) одной из фаз ВИМ с ростом фазного тока (/ф 1 < /ф 2 < /ф з < /ф4).
Таким образом, даже без учета особенностей работы силовой и информационной частей системы управления ВИП его развернутое математическое описание получается нелинейным, громоздким и неудобным для анализа. Поэтому становятся неэффективными традиционные методы исследований электроприводов с использованием структурных схем, преобразований координат, векторных диаграмм и т. п.
В последние года появились специализированные компьютерные программы для расчета и исследования ВИП. Наиболее известные из них — Motor-CAD, SPEED, SRDaS и др. Стоимость их достаточно высока, набор решаемых с их помощью задач фиксирован, а внутренняя структура и реализованные в них алгоритмы явля тся собственность авторов и, как правило, не доступны для пользователя. Это исключает их развитие, модернизацию и адаптацию к новым задачам.
Использование численного интегрирования нелинейных дифференциальных уравнений дает возможность исследовать ВИП при минимальном количестве допущений. Однако специфичность ВИП не всегда позволяет записать в компактной аналитической форме функциональные связи между всеми его переменными. К тому же часто представляет интерес не столько количественная оценка эффективности работы привода, сколько его поведение в той или иной ситуации.
Поскольку ВИП находится в стадии активных исследований и в преддверии широкого практического применения в самых разных областях техники, проблемы исследования его будут неизбежно обновляться и углубляться, затрагивая все новые аспекты поведения В . Соответственно, модель должна обладать способностью гибкой адаптации к новым задачам. Это предполагает ее модульное построение в "открытой" для пользователя программной среде.
еречисленными свойствами обладает имитационная модель ( ) ВИП, воспроизводящая логику его функционирования во времени при различных сочетаниях параметров и сигналов управления [3-5]. Имитационное моделирование может быть реализовано с применением лю-бого языка программирования высокого уровня, например FORTRAN, C++ и т. п. Однако для этого необходимо развернутое математическое
описание работы всех функциональных блоков привода. Поэтому рассматриваются возможности использования для имитационного моделирования ВИП современных специализированных программ и приложений к ним, непосредственно ориентированных на моделирование сложных динамических систем: MATLAB-SIMULINK, MODEL VISION STUDIUM, MODELICA и т. п. Установлено, что наиболее полно поставленной цели отвечает широко известный математический пакет MATLAB с приложением SIMULINK.
Допущения. Наиболее важные допущения относятся к ВИМ. В соответствии с соотношением (4) процесс преобразования энергии в ВИМ связан с изменением фазных индуктивностей Lk (0). Кривые Lk (©) могут быть построены по экспериментальным зависимостям Фк (1ф k) при различных положениях ротора 0. В качестве примера на рис. 2, а сплошными линиями показано семейство экспериментальных зависимостей Ф(/ф) одной из фаз реальной ВИМ.
Из рис. 2, а видно, что насыщение магнитопровода ВИМ начинает появляться примерно при токе /нас = 25 А. Нижняя кривая соответствует совпадению осей зубца статора и паза ротора при 0 = 0мин = 0°, верхняя кривая при 0 = 0макс = 45° — положению, когда оси зубцов статора и ротора совпадают. В положении 0 = 0рас = 15° происходит переход от полностью рассогласованного положения зубцов к частичному их перекрытию. Соответственно, при 0 = 0согл = 41° частичное перекрытие зубцов переходит в их полное перекрытие.
Усредненная реальная кривая L(0) для 1ф < /нас показана на рис. 2, б ж ирной линией. Н а рис. 2, б выделено четыре зоны. Зоны I и II при —0рас < 0 < 0мин и 0мин < 0 < 0рас соответствуют рассогласованному положению зубцов. В зоне III при частичном перекрытии
Рис. 2. Аппроксимация магнитных характеристик ВИМ
зубцов Ьрас < L(0) < Lcorji. В начале полного перекрытия зубцов L(0) = Ьсогл. Максимальное значение индуктивности Ьмакс соответствует совпадению осей зубцов статора и ротора при угле 0макс. Из рис. 2, б видно, что в зоне III зависимость L(0) близка к линейной. В зонах I и II и на границе между зонами III и IV зависимость L(0) нелинейна.
Кривую L(0) можно построить по результатам полевых расчетов ВИМ. Однако это связано с достаточно большим объемом вычислений, усложнением модели и, как следствие, со снижением ее быстродействия при компьютерной реализации .Для упрощения модели часто оказывается достаточным использовать аналитическую аппроксимацию реальной кривой L(0), предварительно оценивая значения Ьмакс и Ьмин для исследуемого двигателя.
Аппроксимация реальной зависимости L(0) основывается на результатах предварительных исследований, согласно которым взаимное влияние фаз и падение МДС в стали ВИМ, как правило, можно не учитывать.
Учет локального насыщения зубцов в зоне их перекрытия производится представлением кривых Ф(1ф) в виде кусочно линейных аппроксимирующих зависимостей при различных значениях 0, как показано на рис. 2, а пунктирными линиями. Для L(0) при I < /нас имеем четыре характерные зоны (см. рис. 2, б).
В зоне I при — 0рас < 0 < 0мин имеем уравнение
L(0) = Km ln0 00ра0 + ßR, (5)
0 0рас
в зоне II при 0мин < 0 < 0рас — уравнение
L(0) = Km 1п0р0 — 00 ßR, (6)
0 рас - 0
где KM — масштабный коэффициент, ßR — угловой размер зубца ротора. Кривая, построенная по этим уравнениям, представлена на рис. 2, б тонкой линией. Как видно, аппроксимирующая и экспериментальная кривые достаточно близки.
В зоне III при 0рас < 0 < 0согл зависимость L(0) аппроксимируется прямой линией, уравнение которой имеет вид
L(0) = ЬРас + L™ — LT 0. (7)
0согл 0рас
В зоне IV нелинейность функции L(0) обычно несущественна, поэтому при 0согл < 0 < 0макс имеем
L(0) = Ьмакс = const. (8)
Фазные ток и момент определены, как в работе [4]. При 1ф < I имеем
*(в)
L(0)'
1ф(©) = ^, (9)
где Ф(0) = J и(0)dt, и(0) — фазное напряжение. Уравнение для фазного момента имеет вид
M (в) = 11ф ^. (10)
В режиме локального насыщения зубцов при 1ф > /нас получаем
/ф(в) = /нас+Ф(в) -/-с£(в) . (11)
^рас
Фазный момент, соответственно, находим по формуле
M(в)= (/Нас/ф(в) - I£с) ^LeT•
(12)
Базовая структура модели ВИП. На рис. 3 представлена имитационная модель ВИП с четырехфазной ВИМ, широко используемая автором при исследовании ВИП [3-5]. Поскольку модель включает большое количество элементов с весьма сложными взаимосвязями между ними, функционально законченные ее части выполнены в виде субсистем (Subsystem). Субсистемы Subsystem1-Subsystem4, формирующие законы изменения L(0) каждой фазы по формулам (5)-(8), а также закон изменения напряжения на каждой из них, выполнены идентично.
В блоках Subsystem5-Subsystem8 выполняется расчет фазных токов и моментов ВИМ по формулам (9)-(12). На выходе блока Suml формируется суммарный электромагнитный момент, который с учетом момента нагрузки (блоки Ms, Constant, Switch) дважды интегрируется блоками Integratorl и Integrator2 для получения сигналов скорости и положения. Имитация циклического изменения фазных индуктив-ностей обеспечивается с помощью специального блока коммутации (SubsystemlO).
Эта модель легко трансформируется в различные вида специализированных моделей в зависимости от принятого вида коммутации в реальном приводе (имеется датчик поло ения ротора или его нет) или способа аппроксимации функции L(0). Кроме того, поскольку структура и логика работы имитационных моделей должна быть наиболее близка к условиям работы реального привода, в них предусмотрена возможность имитации временной дискретизации сигналов управления и разброса электромагнитных параметров фаз.
Subsysteml О —►
Jl
G2 03
Constant
Subsysteml ~ -
Lpad G U1 L1
3ubsystem2
SubsysternS
SubsysterrU
Ms
SubsystemS
Lpad
U1
L.1
OJ
M1
EubsystemB
Lpac2 0 U 2 L2 Lpac2 U2 L2 М2 ш
г*
Г*
►
Subsystein7
ЬрасЗ Lpac3 из kj--L.3 М3 ш
G U3 i
L3 J —►
Subsystem8
Lpac4 L.pac4
G U 4 г» U4 ил
L:4 М4 ■■■i
L4 I
Switch
Surn2
Sumí
1 ts
lntegrator2 ►
Intégralo r1
Рис. 3. Структура базовой имитационной модели ВИП в среде MATLAB-SIMULINK
Таким образом, для все- МТТм стороннего исследования ВИП 0.25 разработан базовый вариант имитационной модели и набор взаимозаменяемых моде- 015 лей отдельных элементов при- 0.1 вода. Все эти модели в совокупности образуют дополнительную специализированную 0 библиотеку SIMULINK, ори- 2 з 4 5 б 7 Шс
ентированну на исследова-
Рис. 4. Иллюстрация адекватности имита-ние поведения ВИП в различ- ционной модели ВИП: ных режимах его работы. 1 — макет; 2 — модель
Иллюстрация адекватности модели. Испытаниям подвергался двигатель со следующими параметрами: Рном = 3 кВт, пном = = 3000 об/мин, количество полюсов на статоре NS = 12, количество полюсов на роторе N = 8, сопротивление фазы Дф = 0,89 Ом, число витков фазы w = 220, Ьмин = 0,11 Гн, Ьмакс = 1,1 Гн. Проведено сопоставление реальных кривых фазных токов и суммарного момента с кривыми, полученными при тех же условиях на модели. В качестве примера сплошной линией на рис. 4 показана восстановленная расчетами кривая суммарного момента для двух фаз по экспериментальным кривым фазных токов и напряжений, а пунктирной линией — аналогичная кривая, полученная на модели при напряжении Цф = 28 В и частоте вращения п = 1000 об/мин. Хорошее совпадение кривых свидетельствует о том, что модель достаточно точно воспроизводит характер изменения моментов коммутируемых фаз и может быть использована для исследования В .
Эффективность разработанных имитационных моделей ВИП подтверждена тем, что решены некоторые важные для мировой практики теоретические и практические задачи, вытекающие из продиктованных опытов применения этого привода и препятствовавшие его широкому распространени на современном этапе. алее кратко рассмотрены основные из этих задач.
Алгоритмы управления ВИП при максимальном использовании его энергетических и регулировочных возможностей. Установление связи между предельными энергетическими возможностями ВИП, параметрами ВИМ и управления является важнейшим этапом рационального проектирования привода, а неверно заданные управляющие воздействия могут существенно снизить эффективность его работы. При исследовании условий повышения энергетической эффективности ВИП проведен анализ механических характеристик М(ш),
и = Г/,™= const
Юбаз
Рис. 5. Механические характеристики ВИП
представленных на рис. 5 [6]. Весь скоростной диапазон работы ВИП может быть разделен на две характерные зоны. Граничное значение скорости между этими зонами принято за базовое значение скорости Шбаз (точка а на рис. 5).
Зона I соответствует относительно малым скоростям, где ВИМ обеспечивает постоянство выходного момента за счет работы токоогра-ничения. Алгоритмы управления и предельные возможности ВИП в этой зоне известны. В зоне II действие регулятора тока прекращается и к фазам ВИМ прикладывается неизменное напряжение питания ином = const. Здесь основными управляющими воздействиями являются углы включения ©вкл и коммутации ©ком (начала отключения работавшей фазы).
Рост ш относительно точки а (см. рис. 5) при неизменных параметрах коммутации ВИМ ведет к резкому снижению амплитуды 1ф (©). Это видно из сопоставления кривой 1 на рис. 6, соответствующей ш = шбаз на рис.5 и кривой 3 на рис.6 для ш > шбаз. При этом средний момент и средняя выходная мощность ВИП также снижаются. Вентильно-индукторная машина оказывается недоиспользованной
по выходной мощности (штриховая кривая на рис. 5). Для повышения выходной мощности ВИП, например, как на характеристике, показанной на рис. 5 сплошной линией без штриховки, необходимо увеличивать фазный ток двигателя за счет более раннего включения фазы (кривая 2 на рис. 6).
Теоретическая граница максимальной выходной мощности ВИП имеет место в режиме граничной коммутации ВИМ [6, 7], Рис. 6. Кривые фазного тока ВИМ когда
7вкл + 7раб + Тоткл = TR, (13)
где тя — полюсное деление ротора. При этом ток /ф(0) спадает до нуля к моменту очередного вкл чения этой е фазы, т. е. интервал протекания тока в фазе возрастает до 360°. Режим граничной коммутации является предельной теоретической границей нормальной работ ВИП, поскольку любое, даже очень малое увеличение либо 7вкл, либо 7раб по отношению к этому режиму приводит к одновременному увеличению 7откл. Это нарушает нулевые начальные условия для тока и потока включаемой фазы, т. е. начальное значение потока возрастает с каждым циклом коммутации. В результате среднее значение потока таюке возрастает, что ведет к общему насыщению магнитной системы В и сни ени входного момента.
Соотношение между 7вкл и 7раб в режиме граничной коммутации имеет вид
где 7вкл = Твкл^П 7раб = 7раб^/п.
Условие (14) является необходимым, но недостаточным для обеспечения максимального среднего момента фаз за цикл коммутации. Оно не дает однозначного соотношения между 7вкл и 7раб. Оптимальное соотношение между 7вкл и 7раб, обеспечивающее максимум среднего момента, зависит от многих факторов, в частности от 7вкл, характера изменения /ф(0) на интервале 7раб и от конфигурации и состояния магнитной системы . ля его определения использована имитационная модель ВИП, позволяющая имитировать различные пространственные сдвиги полюсов ВИМ. В результате установлено, что средний момент наиболее распространенной четырехфазной ВИМ при изменении Кь = Ьмакс/Ьрас в диапазоне от 5 до 10 достигает максимального
зНачеНИя, К°1Да 7Вкл = 7Вкл.опт = 1,0 И 7раб = 7р*аб.опт = 0,5 Эти ПаРа-
метры коммутации соответствуют предельному режиму максимальной выходной мощности ВИП.
Рассмотрена одна из важнейших для практики постановок задачи о максимуме выходной мощности ВИП: при полном использовании установленной мощности инвертора. Для этого необходимо, чтобы в любой точке рабочих характеристик привода при номинальном напря ении питания максимальное значение /ф было равно максимальному допустимому значению выходного тока инвертора.
Характер изменения выходной мощности ВИМ в относительных единицах Рср(ш*) при трех уровнях ограничения фазного тока показан на рис. 7. Здесь базовое значение момента определяется из соотношения
(14)
(15)
Рис. 7. Изменение выходной мощности ВИП при различных уровнях ограничения установленной мощности инвертора:
1_г * =3; 2_Т * =2; 3_I * =1
Ф.макс ' ф.макс ' ф.макс
где /нас — ток насыщения ВИМ, принятый за его базовое значение; w — число витков фазы.
Наиболее важной является зона при ш* > 1, где с увеличением скорости мощность Рс*р монотонно растет. Максимальное значение скорости для каждой кривой ограничивается режимом граничной коммутации (кривая со штриховкой). Как видно на рис. 7, наибольший скоростной диапазон работа ВИП в этом режиме соответствует /| макс = 1, т. е. работе ВИМ без насыщения, а увеличение /| макс приводит к резкому снижению этого диапазона.
Найдены условия обеспечения постоянства выходной мощности ВИП [8]. Установлено, что наибольшими возможностями алгоритм управления обладает при изменении только 7вкл. С применением имитационной модели ВИП на примере четырехфазной ВИМ с конфигурацией полюсов 8/6 получены наиболее важные для практики количественные оценки этого режима. Установлено, что диапазон изменения скорости при Рс*р = const и прочих равных условиях зависит от состояния магнитной системы ВИМ. Наилучшими возможностями в этом смысле обладает линейный режим, т. е. режим без насыщения. В частности, для наиболее распространенной ВИМ с типовыми параметрами этот диапазон может составлять 14 и более.
Алгоритмы компенсации пульсаций момента. Пульсации момента ВИП провоцируют такие нежелательные явления, как шум, вибрации, ограничения в точности воспроизведения заданных параметров движения (особенно для малоинерционных нагрузок), а таюке нарушения нормального пуска привода (при значительных провалах момента).
ричина пульсаций момента состоит в несогласованном изменении токов во включаемой /вкл и отключаемой /откл фазах ВИМ. Наибольшие
пульсации момента имеют место, когда одновременно оказываются под током три фазы, например в зоне высоких скоростей или при несимметричной коммутации фаз.
ценка физической реализуемости полного устранения пульсаций момента ВИМ проведена на имитационной модели при анализе двух частных случаев, имеющих практическое значение [9].
Формирование 1вкл при постоянном напряжении на отключаемой фазе иоткл = const. Ток 1вкл формируется так, чтобы сумма моментов Мвкл и Моткл включаемой и отключаемой фаз соответственно оставалась неизменной и равной М^ (рис. 8, а). Для этого, как видно из рис. 8, а, необходим нелинейный закон изменения напряжения ивкл на включаемой фазе с бесконечно большой его форсировкой в начале коммутации.
ля того чтоб пульсации момента полность отсутствовали, необходимо, чтобы момент Мвкл(0) включаемой фазы нарастал с такой же скоростью, с какой снижается момент Моткл(0) отключаемой фазы, т. е.
dMBI0I (0) ¿Моткл(0)
de de
(16)
В первое мгновение после начала коммутации ток /откл отключае-
„ , dIоткл(в)
мои фазы и производная-—— не равны нулю. 1 юэтому производ-
d0
dMоткл(в)
ная-—-такке отлична от нуля. Однако ток 1вкл при этом равен
d0
б
Рис. 8. Изменение параметров В^М при приложении к отключаемой (а) и включаемой (б) фазам прямоугольного импульса напряжения
нулю и условие (14) выполняется при бесконечно больших производ-
„ ¿/вкл(0) тт „
ной-—— и напряжении ивкл. 1юскольку реально такое условие не-
а©
выполнимо, сделан вывод о том, что если для отключаемой фазы в пер-
d/откл (©)
вое мгновение после начала коммутации производная-отлич-
а©
на от нуля, то никакой физически реализуемый закон изменения напряжения на включаемой фазе не обеспечит полного устранения пульсаций момента. В кривой суммарного момента М^ в начале коммутации фаз будет наблюдаться провал.
Формирование 1атя при постоянном напряжении на включаемой фазе UeKn = const. В этом случае, как видно из рис. 8, б, необходим нелинейный закон изменения напряжения иоткл при бесконечно большом возрастании его в конце коммутации. Поскольку это такке невыполнимо, в конце коммутации будет наблюдаться всплеск суммарного момента М^.
Сопоставление результатов рассмотренных случаев коммутации фаз позволило заключить, что алгоритмы компенсации пульсаций момента могут быть построены на их комбинации. Для расчета текущих значений напря ений на коммутируемых фазах из условия компенсации пульсаций момента в структуру управления приводом целесообразно включить его математическую модель. Основное требование к модели В в данном случае, кроме ее адекватности реальному приводу, состоит в том, что расчетные максимальные значения напряжений на фазах должны быть в пределах напряжения силового источника питания. Предложены два варианта построения таких алгоритмов.
Первый алгоритм управления с непосредственным формированием напряжений фаз. В начале коммутации на включаемую фазу подается напряжение силового источника питания ипит. Напряжение иоткл задается математической моделью привода. Это продолжается до тех пор, пока расчетное значение иоткл не достигнет уровня ипит. Как только это произойдет, для отключаемой фазы иоткл = ипит, а ивкл изменяется по закону, рассчитанному на модели с теми е условиями.
Чередование описанных режимов управления позволяет всякий раз подстраивать напряжение на одной из фаз под условие подавления пульсаций момента. Однако в окрестности точки начала перекрытия зубцов мо ет возникать ситуация с одновременным ограничением как напряжения ивкл, так и иоткл, в результате чего появляется отклонение в суммарном моменте М^. На рис. 9 пунктиром показаны такке кривые токов /вкл.расч и /откл.расч на включаемой и на отключаемой фазах, рассчитанные из условия полного подавления пульсаций момента, а сплошными линиями — кривые токов /вкл и /откл этих фаз, сформи-
Рис. 9. Компенсация пульсаций момента ВИМ
рованные по описанному выше алгоритму. Как видно, /вкл.расч и /вкл, а такке /отк расч и /откл не совпадет, что и является причиной отклонений в суммарном моменте М^.
Для минимизации ошибки в М^ предложено формировать управляющие воздействия непосредственно в функции отклонений фазных токов от их расчетных значений.
Второй алгоритм управления с непосредственным формированием токов фаз. На математической модели ВИП рассчитываются законы изменения токов /вкл и /откл, а такке напряжения иоткл. Как и в первом алгоритме, на начальном этапе коммутации на включаемую фазу подается полное напряжение ипит. Отличие состоит в том, что иоткл формируется в функции рассогласования между расчетным /откл.расч и реальным /откл токами этой фазы. Рассчитываемое на модели напряжение для отключаемой фазы является лишь критерием перехода на другой режим управления.
Наиболее просто этот алгоритм реализуется с применением релейных регуляторов. Результаты моделирования показали, что ошибка в токе любой из фаз, возникшая из-за неспособности источника питания обеспечить необходимый темп изменения тока, является кратковременной и компенсируется релейным регулятором.
Исключение аномальных режимов работы ВИП. Работа ВИП в нормальном режиме предполагает соответствие при допустимых отклонениях законов изменения его регулируемых координат (фазных токов, моментов, скорости и т.д.) исходным данным на разработку привода. Это основные расчетные режимы ВИП, по которым выбираются его элементы и определяются алгоритмы управления. В аномальных режимах отклонения от расчетных режимов бывают значительными, что ведет к нарушению хода обеспечиваемого ВИП технологического процесса, а иногда и к аварийной ситуации.
Рассмотрены две структуры управления ВИП: с коммутацией ВИМ по командам датчика положения и бездатчиковое управление [3]. Для них предложен удобный и информативный способ детектирования аномальной ситуации, основанный на информации о степени отклонения фазных токов 1ф ВИМ от расчетных значений. При датчиковом управлении наиболее существенное влияние на отклонения амплитуды 1ф оказывают неучет временного квантования сигналов управления и пренебрежение разбросом электромагнитных параметров фаз ВИМ [10].
Установлено, что при этом способе управления временная дискретизация сигналов мо ет привести только к запаздывани команд на коммутацию ВИМ, что приводит к снижению амплитуды тока 1ф относительно расчетного значения. При управлении ВИП от микропроцессорной систем период временной дискретизации сигналов управления Ткв определяется ее архитектурой, системой команд и тактовой частотой, а так е особенностями проектируемого алгоритма управления ВИП .Для большинства реальных ВИП значение Ткв оказывается равным периоду широтно-импульсной модуляции ТШИМ напряжения на фазах ВИМ.
Исследовано влияние квантования сигналов по времени на характер взаимных изменений амплитуд 1ф соседних фаз ВИМ при их последовательной коммутации, поскольку, кроме снижения развиваемого момента, наличие периодичности в этих изменениях мо ет провоцировать нежелательные вынужденные колебательные процессы в приводе. Установлено, что характер изменений амплитуд тока 1ф зависит от соотношения длительностей Ткв и Тком — периода коммутации фаз ВИМ. Если значение Тк*м = Тком/Ткв — целое число, то амплитуда токов 1ф соседних фаз не изменяется, а если Тк*м не является целым числом, то амплитуда тока 1ф в соседних циклах коммутации фаз изменяется, при-
чем характер этих изменений зависит от значения дробной части Тком.
Таким образом, при повышенных требованиях к точности формирования фазных токов В в ироком скоростном диапазоне работы в программной части ВИП рекомендован переход от значения Ткв, определяемого ТШИМ на низких скоростях, к формированию Ткв специальной подпрограммой обслуживания прерываний по фронтам сигналов ДП на высоких скоростях.
Разброс электромагнитных параметров фаз ВИМ проявляется, главным образом, в изменении максимального значения индуктивности фазы Ьмакс. Реальные отклонения Ьмакс от расчетного значения достигают ±(10... 15)% и более. Установлено, что к заметному изменению амплитуды 1ф приводит лишь уменьшение Ьмакс относительно расчетного значения только при условии, что в зоне перекрытия зубцов до момента отключения фазы ток 1ф близок к постоянному значению или нарастает. Если же в зоне перекрытия зубцов ток 1ф спадает (наиболее часто встречающийся случай), то уменьшение Ьмакс практически не сказывается на амплитуде тока 1ф. Изменяется незначительно лишь форма кривой 1ф в зоне перекрытия зубцов.
Исследован наиболее распространенный способ бездатчикового управления, основанный на сравнении текущего значения потокосцеп-ления фазы с заданным его значением. В результате моделирования установлено, что в нем, кроме временной дискретизации сигналов и разброса параметров фаз, к аномальным режимам могут приводить отклонения в значениях настроечных коэффициентов алгоритма. При определенных их сочетаниях из-за накопления ошибки по углу включения фаз В цикл ее коммутации сдвигается в сторону зон создания тормозного момента с резким увеличением фазного тока. Эта ситуация имеет место в зоне II для экспериментальных кривых фазных токов, представленных на рис. 10. В зонах I и III коммутация фаз происходит обычным образом. Исследования показали, что при отсутствии таких отклонений условия коммутации ВИМ в датчиковом и в бездатчиковом вариантах одинаков .
Установлено, что причиной смещения цикла коммутации ВИМ в зону тормозного момента в режимах, близких к холостому ходу, может быть наличие "мертвой зоны", вводимой для устойчивой работа алгоритма в области малых сигналов. С ростом нагрузки ВИП основной причиной перехода ВИМ в зону торможения является неблагоприятное сочетание масштабных коэффициентов при вычислении текущего значения потокосцепления фазы.
Выводы. 1. Предложен и экспериментально проверен метод исследования с использованием имитационных моделей, воспроизводящих логику работа во времени реального привода, ориентирован-
Рис. 10. Кривые фазного тока в аномальном режиме работы ВИП при бездатчи-ковом управлении
ный на поиск адекватных алгоритмов управления. Разработана базовая имитационная модель ВИП в среде МАТЬАВ-ЗГМиЬШК, являющаяся основой имитационных моделей основных режимов работы привода.
2. В режиме максимальной выходной мощности ВИП при ограничении амплитудного значения фазного тока с увеличением скорости момент спадает медленнее, чем растет скорость, и наблюдается рост максимальной выходной мощности ВИП по мере увеличения скорости. Это позволяет средствами управления обеспечить постоянство выходной мощности ВИП.
3. Определены условия полного устранения пульсаций момента средствами управления. Наилучшими показателями обладает алгоритм управления с формированием фазных напряжений в функции отклонений реальных фазных токов от расчетных значений, определенных в реальном времени из условия компенсации пульсаций момента.
4. При наиболее распространенном алгоритме бездатчикового управления ВИМ возможны по разным причинам как более ранняя, так и более поздняя коммутации фаз. айден наиболее неблагоприятные сочетания параметров ВИП, вызывающие резкое превышение фазным током ВИМ расчетного значения, и определены условия их исключения.
5. Возможность учета в имитационных моделях наиболее специфических особенностей ВИП подтверждает их широкие функциональные возможности. Это позволяет эффективно использовать их и при оптимизации параметров привода на стадии его проектирования, выборе
рациональной структуры и параметров замкнутого управления приводом, организации генераторного режима работы ВИМ в зависимости от конкретных условий, характера нагрузки и т. п. Значение предложенных методов и средств исследования ВИП при решении подобных задач будет неизбежно возрастать по мере распространения этого нового, перспективного типа электропривода.
СПИСОК ЛИТЕРАТУРЫ
1. Ильинский Н. Ф. Перспективы развития регулируемого электропривода // Электричество. - 2003. - № 3. - С. 2-7.
2. M i ll e r T. J. E. Switched Reluctance Motors and Their Control. - Oxford: Magna Physics Publishing and Clarendon Press, 1993. - 205 с.
3. Красовский А. Б. И митационные модели в теории и практике вентильно-индукторного электропривода: Автореф. дис. ... д-ра техн. наук. - М.: МЭИ, 2004.-40 c.
4. Красовский А. Б. Применение имитационного моделирования для исследования вентильно-индукторного электропривода // Электричество. - 2003. - № 3.
- С. 35-45.
5. К р а с о в с к и й А. Б. Базовая имитационная модель вентильно-индукторного электропривода // Системы управления и информационные технологии. - 2004.
- № 3.- С. 13-17.
6. Красовский А. Б. Получение максимальной выходной мощности вентильно-индукторного электропривода средствами управления // Электричество. - 2002.
- № 9.- С. 29-36.
7. Krasovskii A.B. Obtaining the maximum output power of a switched reluctance electric drive by control means // Electrical Technology Russia. - 2002. - № 3. -P. 107-120.
8. Красовский А. Б. Анализ условий формирования постоянства выходной мощности в вентильно-индукторном электроприводе // Электричество. - 2002. -№ 2. - С. 36-46.
9. Красовский А. Б., Б ы ч к о в М. Г. Исследование пульсаций момента в вентильно-индукторном электроприводе // Электричество. - 2001. - № 10. -С. 33^4.
10. Красовский А. Б. Аномальные режимы в вентильно-индукторном электроприводе при датчиковом варианте управления // Вестник МГТУ им. Н.Э. Баумана. Серия "Машиностроение". - 2003. - № 2. - С. 85-103.
Статья поступила в редакцию 25.04.2005
Александр Борисович Красовский родился в 1954 г., окончил в 1977 г. Московский энергетический институт. Д-р техн. наук, профессор кафедры "Электротехника и промышленная электроника" МГТУ им. Н.Э. Баумана. Автор около 70 научных работ в области дискретного электропривода.
A.B. Krasovsky (b. 1954) graduated from the Moscow Institute for Power Engineering in 1977. D. Sc. (Eng.), professor of "Electrical Technology and Industrial Electronics" department ofthe Bauman Moscow State Technical University. Author of more than 70 publications in the field of discrete electrical drives.