CC BY
37
ЭНЕРГЕТИКА, ЭНЕРГЕТИЧЕСКОЕ МАШИНОСТРОЕНИЕ И ЭЛЕКТРОТЕХНИКА
УДК 621.313.3
Ю. В. Домрачева, С. Ю. Логинов
ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ СИНХРОННОЙ ИНДУКТОРНОЙ БЕСПОДШИПНИКОВОЙ МАШИНЫ
Рассматриваются принципы построения модели синхронной индукторной беподши-никовой машины, позволяющей проводить исследование влияния допущений, принятых при построении математического описания, на работу системы управления. Приведены основные соотношения для создания программной модели машины. Проведена сравнительная оценка возможностей моделей различных типов.
Ключевые слова: бесподшипниковая электрическая машина, программная модель, имитационное моделирование.
Преимущества электрических машин с полным бесконтактным магнитным подвесом ротора позволяют использовать их в тех областях, в которых применение машин с традиционными опорами не возможно из-за наличия смазки и механического контакта между вращающимися и неподвижными частями механизма.
В настоящее время в конструкции таких машин широко используются активные магнитные подшипники (АМП). Для осуществления полного магнитного подвеса ротора необходимо два радиальных и один осевой АМП или два конических АМП (см. рис. 1.а). Установка АМП требует увеличения длины ротора, что затрудняет использование такой машины в высокоскоростных установках из-за уменьшения значений критических скоростей [1].
а) Осевой
АМП ЭМ АМП д]ущ
11' II |~М1 Г
АМП ЭМ АМП
б) Осевой
БЭМ БЭМ
АМП
БЭМ БЭМ
Рис. 1. Варианты компоновки электрической машины с полным магнитным подвесом ротора
Следующим шагом в развитии электрических машин с полным магнитным подвесом ротора являются бесподшипниковые электрические машины (БЭМ). Магнитное поле в зазоре модуля БЭМ создаёт как вращающий момент, так и управляемые радиальные силы, удерживающие ротор в центральном положении. Для создания такого поля в пазы статора укладывается две обмотки: обмотка привода и обмотка подвеса. Благодаря тому, что функции электродвигателя и АМП объединены в одном модуле, длина ротора, при той же мощности привода и подвеса, может быть сокращена [2]. Уменьшение габаритов связано с тем, что магнитное поле привода используется как составная часть поля подвеса. Это также позволяет снизить энергопотребление. Варианты компоновки БЭМ показаны на рис. 1. б).
Существенным недостатком БЭМ является взаимное влияние магнитных сил и вращающего момента. Управлять этими величинами раздельно позволяет специальная конструкция модулей БЭМ. В 1974 году был предложен принцип построения конструкции электрической машины, совмещающей функции электродвигателя и АМП [3]. Согласно ему одна из обмоток статора выполняется с числом пар полюсов равным р, а вторая — р+1. Ротор должен иметь то же число полюсов, что и обмотка привода. При таком сочетании количества полюсов суммарная электродвижущая сила (ЭДС), наведённая одной из обмоток во второй, равна нулю, а магнитодвижущая сила (МДС) обмотки подвеса не создаёт вращающего момента. Однако это условие будет выполняться только в том случае, если ротор находится в центральном положении. При смещении ротора распределения МДС обмоток статора становятся несимметричными, вследствие чего обмотка подвеса создаёт вращающий момент, а обмотка привода — радиальную силу. На рис. 2 показан пример расположения обмоток в пазах статора БЭМ.
А — фаза обмотки привода, а — фаза обмотки подвеса
Взаимное влияние привода и подвеса можно уменьшить, используя специальные алгоритмы управления токами обмоток. Для синхронной индукторной БЭМ (рис. 2) предлагалось два варианта построения таких алгоритмов. В одном из них использовался специальный компенсатор, рассчитывающий сигналы задания для токов привода с учётом величины токов подвеса и наоборот [4]. В 2010 году была предложена структура системы управления, построенная на базе инверсной системы [5].
Во всех предложенных системах управления конструктивные параметры машины представлены в виде индуктивностей Ld и Ь(1, которые определяются
магнитными проводимостями зазора над полюсами ротора и между ними соответственно. В [6] зависимости Ld и Ь(1 от перемещений ротора х и у предложено определять экспериментально. Такой подход позволяет получить достаточно точную математическую модель машины для построения системы управления. Но отсутствие аналитических зависимостей индуктивностей от конструктивных параметров БЭМ не даёт возможности составить методику проектирования.
При составлении аналитического описания индукторных БЭМ, кроме обычных для электрических машин допущений, так же принимаются следующие:
1. Величина воздушного зазора описывается приближённым выражением.
2. Магнитная проводимость воздушного зазора между полюсами ротора равна нулю.
Необходимая точность построения математической модели машины определяется работоспособностью системы управления, построенной на её основе.
Проверка работоспособности системы управления не может производиться на модели машины, построенной с теми же допущениями. Для этой цели удобно использовать такую модель, в которой любой из неучтённых признаков можно исключить или добавить в рассмотрение. В данной статье рассматривается один из вариантов построения такой модели для исследования системы в среде МЛ^ЛВ 81ши11пк. Все вычисления выходных переменных блоков модели будут производиться с помощью языка программирования С++, в связи с чем в дальнейшем будем называть её программной моделью.
Основные соотношения, используемые в модели синхронной индукторной БЭМ.
Модель синхронной индукторной бесподшипниковой машины (СИБМ) можно разделить на две подсистемы: электромагнитную и механическую (см. рис. 3). Электромагнитная подсистема в зависимости от положения ротора и токов обмоток статора формирует вращающий момент и радиальные силы, которые в свою очередь являются входными сигналами для механической подсистемы. Выходными координатами модели являются угол поворота р и радиальные перемещения ротора х и у.
М1
12а
12Ъ
^4а
иъ
магнитная
подсистема
М
Л
Я,
V
Механическая
подсистема
х
Л
Рис. 3. Структура модели БЭМ
В блоке электромагнитной подсистемы момент и радиальные силы должны вычисляться по заданным токам с учётом конфигурации обмоток и сердечников статора и ротора. Эти величины вычисляются, как частные производные магнитной энергии по соответствующей координате. Чтобы получить аналитическую зависимость магнитной энергии от угла поворота и смещений ротора необходимо иметь выражение для величины воздушного зазора, учитывающее её изменение вдоль этих координат. Как уже отмечалось, для построения системы управления пользуются приближенным выражением. Однако в модели должна быть возможность использования любого способа вычисления величины зазора. Поэтому электромагнитный момент и радиальные силы будем определять как равнодействующие моментов и сил, действующих на некоторые элементарные части поверхности ротора.
Рис. 4. Фрагмент воздушного зазора СИБМ
Разобьём поверхность ротора на элементы, длина которых равна длине активной части сердечника, а ширина равна длине дуги, опирающейся на угол Аф (см. рис. 4). Угол Аф должен быть настолько мал, чтобы с определённой степенью точности можно было считать, что величина зазора над элементом одинакова во всех его точках. Величина зазора, в зависимости от угловой коор-
176
динаты ф и радиальных перемещений x и у, описывается следующим выражением [1].
5 = 50 - x • cos ф- у • sin ф = 50 - 5', (1)
где 80 — величина зазора при центральном положении ротора,
5' = x • cos ф + у • sin ф — величина постоянная в пределах Дфг. Электромагнит-
ный момент и радиальные силы, действующие на элемент Дф., находятся численно по следующим выражениям.
Mг =^, (2)
г Дф
г ДЩ г ДЩ
Fxi =---L, Fvi =-----, (3)
xi Дх уг Ду
где Mt — электромагнитный момент, Fxi и Fyi — радиальные силы, ДЩ —
магнитная энергия, запасенная на участке воздушного зазора машины, ограниченном углом Дф .
1 N N
ДЩ = 2tt^ • ij ■ ik, (4)
2 j=1 k=1
где N — число обмоток машины, i — ток обмотки, ALt — собственная или взаимная индуктивность обмоток, определяемая элементом зазора Дф-.
. ii0 • Дф. • r • i
* 8 -8--------Wj 'Wi, (5)
где i0 = 4п • 10-7 — магнитная проницаемость вакуума, S = r • Дф • l — площадь элемента Дф, r — радиус ротора, l — длина активной части сердечника,
w — число витков соответствующей обмотки.
Результирующее значение электромагнитного момента и радиальных сил
M = £м1, (6)
г=1
Fx = itFx., Fy =±Fy,, (7)
i=1 i=1
где n — общее число элементов Дфг.
По уравнениям (1)-(7) производятся расчеты выходных величин электромагнитной подсистемы на каждом шаге расчета модели.
Угол поворота ф и радиальные перемещения ротора в блоке механической подсистемы можно найти с помощью дифференциальных уравнений движения.
m • x = Fx - Fxf ,
m • у = Fy - Ff , (8)
J • р = M - Mf,
где т — масса ротора, J — момент инерции ротора, х и у — ускорения ротора в направлении осей х и у соответственно, ф — ускорение вращательного движения, FXJ■, Fyf и М ^ — возмущающие силы и момент. Из уравнений (8)
получим выражения для расчета перемещений и угла поворота на 7-ом шаге.
ЛТ7
Х = — (Fx - Fxf )+ 2Х j-1 - Х
Vі x ± xf ^ j-1 лі-2'
m -k
= —t(Fv - Fvf)+ 2Vj-1 - Vj-2. (9)
у'j=m - '+2 уj-1- у*-2,
ф j = J M - M f )+ 2 ф j-1 - ф j-2 ’ где Tk — шаг расчета, с, xj-1, уj-1 и ф j-1 — значения перемещений и угла поворота на (j-1)-ом шаге расчёта, xj _2, уj _ 2 и ф j _ 2 — значения перемещений и
угла поворота на (j-2)-ом шаге расчёта.
Проверка алгоритма системы управления СИБМ на программной модели, построенной с помощью уравнений (1)-(9), даёт результаты более приближённые к результатам испытаний на реальной машине, чем при использовании упрощенных аналитических выражений. Для повышения точности расчёта выходных параметров модели необходимо, по возможности, уменьшать шаг расчёта и величину угла Дф .
Кроме того, программная модель даёт возможность учитывать или исключать из модели те признаки, влияние которых на работу системы необходимо исследовать. Рассмотрим способы учёта некоторых допущений, которые обычно принимаются для индукторных бесподшипниковых машин.
Приближённое выражение для обратной величины зазора.
Для определения аналитических зависимостей индуктивностей обмоток статора от радиальных перемещений ротора приходится использовать приближённое выражение для обратной величины зазора машины. Связано это с тем, что результат интегрирования при нахождении этих выражений по точной формуле неудобен для их использования в системе управления.
Точное выражение для обратной величины зазора в зависимости от смещений и угловой координаты выглядит следующим образом.
1 =----------1-----------= —1—. (10)
5 50 - x • cos ф- у • sin ф 50 - 5'
Приближённые выражения для обратной величины зазора получают, представляя её как сумму геометрической прогрессии.
1 1 11 1 U F с“2
-= = — (1 + 5 + 5 +.
= — .(1 + S + S2 +...), (11)
0 0 1 - S S ’
где S = S' /S0 < 1.
S So - S' So 1 - S So
Обычно для описания БЭМ в выражении (11) используются только слагаемые первого порядка [6]. В этом случае величина зазора будет вычисляться достаточно точно лишь при малых перемещениях ротора («20% от д0). В остальных случаях возникает существенная погрешность. При учёте слагаемых второго порядка, можно рассчитывать зазор при перемещениях « 30% от д0 и т. д.
Чтобы при моделировании исключить влияние погрешности расчёта обратной величины зазора, необходимо определять её по тем же выражениям, что и система управления. В предлагаемой программной модели, это можно сделать с помощью пересчёта значения д'. Приравнивая точное значение обратной величины зазора приближённому, получаем следующее соотношение.
д' = дод. (12)
1 + д
Аналогично при учёте в (11) слагаемых второго порядка.
д' = до •дд11. (13)
0 1 + д + д2
Магнитная проводимость воздушного зазора между полюсами ротора.
В математической модели СИБМ, используемой для построения системы управления, магнитная проводимость воздушного зазора между полюсами ротора не учитывается. В предлагаемой модели это означает, что суммирование (6) и (7) производится только по элементам, находящимся на поверхности полюса: от угла ф1 до ф2 (см. рис. 5).
Рис. 5. Фрагмент воздушного зазора СИБМ
Для учёта проводимости межполюсного пространства, в суммы (6) и (7) необходимо включить слагаемые, учитывающие соответствующие воздействия. Например, при учёте боковых проводимостей полюсов, суммирование производится так же от угла ф2 до ф3. В этом случае величина зазора (1) увеличится на , которая в свою очередь может быть рассчитана исходя из геометрии возможного пути протекания потока [7]. В пределах одного элемента Аф., величи-
на И. является константой, поэтому выражения для расчета модели не изменятся. Следует только учитывать, что
д = д0 - д = д0 + \ - д'. (14)
Высшие гармоники магнитодвижущей силы.
Вычисление выходных переменных электромагнитной подсистемы происходит численно, поэтому изменение формы МДС не требует доработки алгоритма расчета модели. Необходимые гармоники добавляются к основной в выражении для магнитной энергии, как дополнительные слагаемые.
Зубцовые гармоники могут быть учтены добавлением к МДС гармоники порядка, равного числу пазов статора. В случае неравномерности распределения пазов статора вдоль зазора или существенной разницы между шириной паза и зубца эту гармонику следует учитывать изменением проводимости зазора, т. е. добавлением к д0 в (1) соответствующего слагаемого.
Насыщение стали сердечников.
Магнитное сопротивление стали сердечников должно быть учтено включением последовательно сопротивлению зазора. В предлагаемой модели удобно привести его к величине зазора.
Rд + Кс = + Rc = д0 - д'+ 00 •5 • К , (15)
°0-5 °0-5
где Rд — магнитное сопротивление воздушного зазора, Кс — магнитное сопротивление сердечников машины. Тогда величина зазора
д = д0 - д' + М0 • 5 • Rc, (16)
Насыщение стали можно учесть с помощью изменения её магнитного со-
противления в зависимости от прикладываемой МДС.
Рс/ = р-------°0 •5 • Пс'-1--, (17)
с/ 2 д0 - д + 00 • 5 • Rcj-1
где р — падение МДС в сердечниках статора и ротора на /-ом шаге расчета, р — суммарная МДС обмоток, Я--1 — магнитное сопротивление стали сердечников на (/-1)-ом шаге расчёта.
Зависимость Rc = /(Рс) составляется по кривой намагничивания стали. При расчете модели программно, её удобно задавать как численно, так и с помощью аппроксимирующей функции. После расчёта р определяется новое
значение Яс/-.
К/ = / (р//), (18)
Вычисление р по (17) вносит погрешность в определение сопротивления стали из-за использования Кс, определённого на предыдущем шаге расчёта. Поэтому при резком изменении , для повышения точности вычисления необходимо повторить операции (17) и (18) несколько раз, используя новые значения К*..
Ч
Сравнение возможностей методов моделирования.
Исследование электромеханических систем можно производить с помощью различных типов моделей: аналитической модели, модели, построенной методом конечных элементов, программной модели. Каждая из них имеет свои преимущества и недостатки. Сравним их с точки зрения исследования влияния допущений на работу системы.
Таблица 1
Сравнительная оценка типов моделей
^"^^^^^Модель Допущения"^^^^ Аналитическая модель Модель, построенная МКЭ Программная модель
Упрощённые выражения для величины зазора сложно невозможно просто
Проводимость зазора между полюсами сложно возможно возможно
Высшие гармоники МДС сложно сложно просто
Насыщение стали сердечников сложно просто возможно
В табл. 1 дана оценка возможности учёта и исключения из моделей различных признаков исследуемой системы.
Составление аналитической модели подразумевает получение системы уравнений, описывающей все необходимые свойства объекта. Добавление или исключение одного из них означает получение новой системы уравнений. В случае СИБМ этот процесс сложен даже при всех сделанных допущениях. Получение аналитических описаний при различных сочетаниях свойств машины — нерационально, т. к. есть более простые способы моделирования.
МКЭ является самым распространённым и простым способом моделирования электромагнитных полей. Модель, построенная этим методом, позволяет достаточно точно рассчитывать силы электромагнитного взаимодействия, не упрощая при этом никаких конструктивных признаков машины, свойств материалов и формы МДС в зазоре. Но исключение из модели каких-либо признаков не всегда возможно или представляет значительные затруднения.
С помощью программной модели, благодаря гибкости алгоритма расчёта, можно достаточно просто организовать изменение состава признаков иссле-
дуемого объекта. При этом, конечно же, не обойтись и без других методов моделирования. Например, для построения зависимости Rc = f (Fc) удобно пользоваться МКЭ, который позволит учесть не только общее, но и локальное насыщение стали сердечников.
Использование программной модели объекта управления расширяет возможности исследования влияния его свойств на работу системы.
Литература
1. Журавлёв Ю. Н. Активные магнитные подшипники: теория, расчет, применение / Ю. Н. Журавлёв. СПб.: Политехника, 2003. 206 с.
2. Логинов С. Ю. Сравнение габаритных показателей бесподшипниковой индукторной машины и двигателя с активными магнитными подшипниками / С. Ю. Логинов, Ю. В. Дом-рачева // Электротехнические комплексы и системы управления. 2012. № 2. С. 20-24.
3. Hermann P. K. «A Radial Active Magnetic Bearing Having a Rotating Drive», London Patent No. 1 500 809, February 9, 1974.
4. Design and Simulation of Control System for Bearingless Synchronous Reluctance Motor / Hannian Zhang [и др.] // Electrical Machines and Systems, 2005. ICEMS 2005: proceedings of the Eighth International Conference. 2005. Vol. 1. P. 554-558.
5. Decoupling control of bearingless synchronous reluctance motor based on inverse system method / Huangqiu Zhu [и др.] // Robotics & Control Systems: Control and Decision Conference (CCDC). Chinese, 2010. P. 2120-2125.
6. Chiba A. Magnetic Bearings and Bearingsless Drives / A. Chiba, T. Fykao, O. Ichikawa, M. Oshima, M. Takemoto, D.G. Dorrell; ELSEVIER. 2005. 381 р.
7. Сливинская А. Г. Электромагнитны и постоянные магниты / А. Г. Сливинская. М.: Энергия, 1972. 248 с.
Об авторах
Домрачева Юлия Вячеславовна — старший преподаватель кафедры электропривода и систем автоматизации ФГБОУ ВПО ПсковГУ.
E-mail: juli-politeh@yandex.ru
Логинов Сергей Юрьевич — старший преподаватель кафедры электропривода и систем автоматизации ФГБОУ ВПО ПсковГУ, канд. техн. наук.
Y. V. Domracheva, S. Y. Loginov
SIMULATION TECHNIQUE OF SYNCHRONOUS RELUCTANCE
BEARINGLESS MACHINE
Construction principles for synchronous reluctance bearingless machine model creation are considered. This model allows investigating influence of mathematical modeling assumptions to control system performance. Basic relations for program model creation are presented. Comparative assessment of various machine model types is made.
Keywords: computer simulation, bearingless electric machine, simulation technique, program model.
