CC BY
12
УДК 621.313
Численная модель коммутации коллекторных электрических машин в динамических режимах работы
А.И. Тихонов, Ю.Б. Казаков, доктора техн. наук, И.М. Лашманов, ст. преп.
Разработана численная модель процессов коммутации коллекторных электрических машин в динамических режимах работы с использованием результатов расчетов магнитных полей методом конечных элементов.
Ключевые слова: коммутация, коллекторные машины, динамические режимы работы.
Numerical Model of Switching Commutator Machines in Dynamic Models
A.I. Tikhonov, Doctor of Engineering, Yu. B. Kazakov, Doctor of Engineering, I. M. Lashmanov, Senior Teacher
The authors developed the numerical model of switching processes of commutator machines in dynamic modes based on calculation results of magnetic fields with finite-element method.
Key words: switching, commutator machines, dynamic models.
Коммутационная надежность остается основной проблемой коллекторных машин (КМ), работающих в динамических режимах, с перегрузками, с частыми включениями и реверсированиями, особенно для неявнополюсных машин постоянного тока (НПМПТ) с распределенными обмотками индуктора [1] и при работе на переменном токе.
В НПМПТ смешанного возбуждения параллельная и последовательная обмотки возбуждения (ОВ), имеющие свои токи, числа витков и провода, могут иметь разновитковые катушки, разные их числа и шаги, быть распределены по разным пазам. Формы распределения намагничивающих сил ОВ для разных их размещений в зависимости от тока якоря меняются.
Конструктивные отличия НПМПТ - распределенность ОВ главных полюсов (ГП) на несколько катушек; распределенность магнитного потока ГП; малый и одинаковый зазор под ГП и добавочными полюсами (ДП) оказывают влияние на величину потока ГП и характер коммутации. Из-за разного числа пазов индуктора и якоря, неравномерного распределения компенсационной обмотки (КО) по пазам индуктора полной компенсации реакции якоря со стороны КО получить не удается. Возникающие местные некомпенсированные поля обмотки якоря (ОЯ) ухудшают коммутацию.
В НПМПТ зачастую наблюдается затрудненная коммутация даже в режиме холостого хода. Это указывает на влияние поля ГП на коммутацию. Для оценки влияния особенностей НПМПТ на коммутацию целесообразно применение численных методов расчета магнитного поля, например метода конечных элементов (МКЭ) (рис. 1). Так, для НПМПТ [4] мощностью 40 кВт выполнены расчеты магнитных полей и индукции в зоне коммутации (рис. 2).
Рис. 1. Конечно-элементная модель КМ в окне библиотеки конечно-элементного моделирования магнитного поля
1 ь 1 1 1 1 изк . ,
1
й /
1 о о ^— _ J
г
J
-10 -7.5 -5 -2.5 0.0 2.5 5 У, град
Рис. 2. Распределение индукции потока ГП НПМПТ в зоне коммутации
Искрение на коллекторе определяется ЭДС, наводимыми в коммутирующих секциях. В коммутируемых секциях ОЯ от поля ГП наводится ЭДС, усложняющая коммутацию. Добавка к реактивной ЭДС при номинальной нагрузке может достигать 10-15 %. При этом эта ЭДС за период коммутации Т может ме-
нять знак четыре раза. Это затрудняет настройку коммутации.
Для генераторного режима явнополюсных машин постоянного тока (ЯМПТ) поле ГП во второй половине периода коммутации во время Т/2 < t < Т наводит в коммутируемой секции ЭДС такого направления, которое способствует ускорению коммутации, а для двигательного режима -замедлению коммутации. Наоборот, у НПМПТ в двигательном режиме в период Т/2 < t < Т при ЬдП = ЬЗК в коммутируемых секциях наводится ЭДС, способствующая ускорению коммутации, а для генераторного режима - замедлению коммутации. В период 0 < t < Т/2 все наоборот. То есть, у неявнополюсных двигателей постоянного тока (НПДПТ) из-за влияния поля ГП происходит снижение плотностей тока и под сбегающим, и под набегающим краями щетки. У неявнополюсных генераторов постоянного тока (НПГПТ) из-за влияния поля ГП происходит увеличение плотностей тока и под сбегающим, и под набегающим краями щетки.
Актуальной является задача расчета ЭДС в коммутируемых секциях с учетом особенностей конструкции КМ. Перспективным методом расчета является комбинирование численных моделей электрической цепи и магнитного поля машины, где модель поля используется для расчета параметров электрической цепи на каждом шаге интегрирования ее по времени [2]. При этом основная сложность, возникающая при моделировании коммутации, состоит в постоянном изменении схемы электрической цепи в ходе расчетного процесса, вызванном переключениями секций ОЯ на коллекторе при его перемещении.
Для расчета параметров электрической цепи КМ использована динамически подключаемая библиотека конечно-элементного моделирования магнитного поля [3]. С помощью этих функций в среде МаШЬаЬ организована численная модель процессов коммутации.
Электрическая схема замещения КМ с учетом коммутации при числе пар полюсов р = 2 в общем случае представлена на рис. 3.
Рк+11 <> 1---------=*-----------1
\Е2 ет-2
іГ , —-—----- V--> I
----—----------1
Рис. 3. Схема замещения КМ (Яь И2 - сопротивления параллельных ветвей обмотки якоря; Еь Е2 - ЭДС, наводимые в параллельных ветвях ОЯ; рд - сопротивление контакта между і-й коллекторной пластиной и щеткой из д-й группы; гэ, еэ - сопротивление и ЭДС э-й секции; И., Ег- сопротивление и ЭДС обмотки возбуждения; и - напряжение сети)
Так как по мере поворота якоря электрическая цепь КМ перестраивается, система диф-
ференциальных уравнений, описывающих динамику машины, должна формироваться на каждом шаге интегрирования.
Алгоритм формирования системы дифференциальных уравнений включает следующие шаги:
1. Формируются списки Пд пластин коллектора, перекрываемых щетками (в простой волновой обмотке все щетки, чередуясь, делятся на две группы по р щеток в каждой).
2. Формируются списки Сд секций, замкнутых щетками д-й группы накоротко.
3. Формируются списки Ву номеров секций, входящих в параллельные ветви ОЯ (простая волновая ОЯ содержит две ветви).
4. Первые два уравнения описывают механические процессы:
м - м0 - мп
(1)
где М - электромагнитный момент; О - частота вращения; J - момент инерции якоря.
5. Два первых контура (по методу контурных токов) проходят через контакты с наименьшим сопротивлением рк1 и рт>2, определяемым по спискам Пд, источник напряжения и и параллельные ветви ОЯ (рис. 4).
6. Остальные контуры проводятся в пределах каждой щеточной группы сопротивления рк1 или рт>2 и остальных контактов данной группы. Направление каждого контура выбирается так, чтобы в контактах рк1 и рт>2 оно совпадало с направлением контуров 1 и 2. Количество таких контуров пд в пределах каждой д-й щеточной группы равно количеству короткозамкнутых секций в д-й группе, номера которых сохранены в списках Сд.
7. Короткозамкнутые секции в каждой группе делятся на две подгруппы Сду с количеством элементов пду, хранящих номера секций, расположенных по отношению к сопротивлениям рк1 и рт>2 со стороны v-х ветвей ОЯ (рис. 4).
8. Формирование списков номеров контуров Кду, образованных секциями из списка С
ду-
9. Формирование подмножеств К(
дуэ
списков Кду, из которых исключаются первые э'-1 элементов, где 5' - порядковый номер, с которым контурообразующая секция э присутствует в списках Сду и Кду.
10. Формирование подмножеств Сдуэ, состоящих из первых э' элементов списка
дуэ-
С
Для расчета ЭДС секций используется библиотека конечно-элементного моделирования магнитного поля, функции которой используются непосредственно в алгоритме интегрирования системы уравнений динамики КМ. Значения токов в секциях и обмотке воз-
буждения, а также угол поворота якоря берутся из расчета системы уравнений. После расчета магнитного поля библиотека возвращает потокосце-пления секций и обмотки возбуждения.
Так как при вращении якоря потокосцепле-ние 5-й секции % = '¥3(І1, і2, іп, а) зависит от токов во всех п контурах и от угла поворота ротора а, можно воспользоваться преобразованием
= --
ГТ
Гі
=-!
р=1 дір — да —
п —о
=-! -ср,
р=і -
где
АТ АТ
I - э Г' .. э —эр - * ; , Сэ
Ар
Аа
Аналогично для ОВ
Е,-Т. = -ї ^ - с, а,
> Г- /—! ,р гИ- ’
- р=1
Гі
(2)
(3)
(4)
где Т,= Т(іьі2,.., іп,а) - потокосцепление ОВ. ЭДС параллельных ветвей ОЯ
nsv пэ^
е=-Т '
э=1 э=1
п (nsv а>т( '\ ,
s,v _
Гі
дТ
=-!
р=Л э=1 и1р
гр
-і
п
I
э=1
дТ 17Т э,у
да
Г а
ГГ
(5)
■Д пр
=-1 —гр-гр - С□, р=1 -
где
пэ\ АТэ
^р *I-fL, С,-I
э=1 Аір
э=1
А^э,у
Аа
(6)
Система уравнений, описывающих динамику КМ с учетом коммутации, имеет вид
-X = иг*.
-
(7)
Матрица коэффициентов М имеет вид
М:
" 1 0 0 0 0 . . 0 "
0 1 0 0 0 . . 0
С1 0 -11 -12 —13 . . —1п
С2 0 —21 2 2 ^1 —23 . . —2п
С3 0 —31 —32 —33 . . —3п
.Сп 0 -п1 -п 2 -п3 . . —пп _
(8)
Вектор неизвестных X имеет вид
X =[а □ і 1 і 2 ... іп-1 іп] ,
(9)
где
—1
Я
АТ, - I АТэ,1 _ I АТэ _ I
и А і А і
АТ
АіЯ АіЯ э^с11 АіЯ эеС21 АіЯ
С ^ АТ^ - I АТэ - I АТэ
1 Аа “ Аа “ Аа “ Аа
э=1 эеСц эеС21
Я=1..п
і АТ, -^-2 АТэ 2 Х-' ат^ ^ і э
—2я -------------------------------------------------------------------------L + I-— + I -* + I -э
2Ч *; /—і *; /—і *; /—і
АТ
АТ
АіЯ ІТІ АіЯ /г? АіЯ /г1 АіЯ
Я э=1 ц эеС12 Я эеС22 Я
Я=1. .п
С2
пэ2
,+I ^—1
АТ пэ2 АТэ2
Аа э=1 Аа
^ АТ_ ^ АТэ
- I э э,
„г Аа — Аа эе^2 эеС22
-I
АТС
эеСш АіЯ
Сп
Я=1..п рєКц
I
АТ
эеС11э
Аа
рєК11
I
АТ
рЯ А
эеС^ Ш Я
Я=1..п рєК12
Ср - I
АТ
эеС12э
Аа
рєК12
рЯ
-I
АТС
эеС21э АіЯ
Я=1..п рєК21
Ср--I
АТэ
эеС21з
Аа
рєК2і
-ря - I
АТ
рЯ А
эеС22э Я
Я=1..п
рєК22
Ср - I
АТ
эеСо'
Аа
рєК22
(10)
Вектор правых частей Р имеет вид
□
М - Мо - мп
^ (1 + і-2)+ + I Гэ I і-1 + I ^|+ I Гэ + I
эеСц ^ У^Кцэ J эеС21 ^ У^К21э
+ рЛ,11 і-1 + і2 + I і] + I ^|+рщ2 I і1 + ^ + I і] + Iі
УєКц УЄК12 J ^ УЄК21 УЄК22
~и + Я, (1 + /2)+ Я22 + I гэ I і_2 + I ^|+ I Гэ М2 + I іу
эеС12 I МК12э ) эеС22 I УеК22э
УєКц уєК12
+ рт,2 1 + ^ + I уєК21 іУ + I 0 Ує К22
+ I іу УєК12 ]* э!.'- [ '1 + I іі УєК11э
I іу УєК12 )* I '• 1 ^ •єС12s ^ 2 + I < уеК12э
+ I УєК22 ']+ I 'э ^ •єС21s і1 + I убК21э
+ I іу ує К22 ]+ I 'э ^ •єС22s і2 + I ч ує К22э
рєК12
рєК21
рєК22
(11)
Расчет матриц Р и М осуществляется на каждом ?-м шаге цикла интегрирования. Алгоритм интегрирования (7) состоит из следующих шагов:
1. Определенные на ?-м шаге интегрирования контурные токи 1д пересчитываются в токи в цепи ОЯ /а, параллельных ветвей ОЯ 1а1 и а2 и коммутирующих секций /5:
Іа1 = І1, Іа2 = І2, Іа = І1 + І2,
is-C11s i*1 + Іj ,
S-C1
's-Cn
'S-Co
i^iis
= -I2 - у /j
j-K12s
=/1 + у /j,
j-K2lS
;=-l2 -Z
j-K22s
(12)
j
2. Рассчитываются элементы вектора Р.
3. Значения токов пересылаются в полевую модель, рассчитывается магнитное поле, определяются потокосцепления ОВ Та и всех секций ОЯ Тэ на ?-м шаге интегрирования и основной поток ф(.
4. Каждому ц-му контурному току дается приращение АіЯ. Рассчитывается магнитное поле, и в секциях ОЯ и ОВ определяются приращения потокосцеплений АТэ- и АТ,(. По (10) рассчитываются элементы —рЯ матрицы М.
5. Якорь поворачивается на угол Аа при неизменных контурных токах, вновь рассчитывается магнитное поле и определяются приращения АТэ и АТ,( секций и ОВ. Рассчитываются элементы Ср матрицы М.
6. По (7) рассчитывается вектор производных по времени от интегрируемых величин. Определяется вектор интегрируемых величин
„ -X ..
+-Г, . <13)
t+At
7. Элементы вектора пересчитыва-
ются в реальные величины, характеризующие процессы на ?+Д^-м шаге интегрирования:
dX1
t+At
а
t+At '
а
dt
dX
•At,
dt
•At, i„
dX
■' I„ +-
2+q
dt
•At.
(14)
К полученным контурным токам применяется преобразование (12).
8. Рассчитывается электромагнитный момент:
Mt = Ст ф1 4(. (15)
9. Процесс повторяется, начиная с п. 1.
Приведенный математический аппарат
позволяет осуществить расчет динамического режима КМ с учетом процессов в коммутирующих секциях, особенностей конструкции магнитной системы и формы питающего напряжения.
В среде МаНаЬ разработана программа, позволяющая формировать и решать систему дифференциальных уравнений, описывающих динамику КМ с учетом процессов коммутации. Индуктивные параметры секций находятся из расчета магнитного поля. Пото-косцепление секции определяется углом поворота ротора, токами якоря в параллельных ветвях, токами возбуждения и коммутируемых секций. Расчет поля и переформирование системы осуществляется на каждом шаге интегрирования. Разработана интерактивная визуальная подсистема для проверки правильности формирования всех списков.
Список литературы
1. Казаков Ю.Б., Лашманов И.М. Особенности и моделирование коммутации в динамических режимах неявнополюсных коллекторных электрических машин: мат-лы междунар. науч.-техн. конф. «Состояние и перспективы развития электротехнологии». Т. 2 / Иван. гос. энерг. ун-т. - Иваново, 2009. - С. 57.
2. Технология численного исследования электрических машин с использованием библиотеки конечно-элементного моделирования магнитного поля / А.И. Тихонов, С.Ю. Кучеров, И.М. Лашманов, Д.В. Рубцов // Вестник ИГЭУ. - 2006. - Вып. 3. - С. 5-8.
3.Тихонов А.И., Казаков Ю.Б. Программирование численного эксперимента с использованием конечно-элементной модели магнитного поля в объектах электромеханики / ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина». - Иваново, 2008.
Лашманов Иван Михайлович
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
старший преподаватель кафедры электромеханики,
адрес: г. Иваново, ул. Рабфаковская, 34,
телефон (4932) 26-97-06,
e-mail: elmash@em.ispu.ru
Тихонов Андрей Ильич,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
доктор технических наук, доцент, зав. кафедрой физики,
адрес: г. Иваново, ул. Рабфаковская, 34,
телефон (4932) 26-97-06,
e-mail: elmash@em.ispu.ru
Казаков Юрий Борисович,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
доктор технических наук, профессор, зав. кафедрой электромеханики,
адрес: г. Иваново, ул. Рабфаковская, 34,
телефон (4932) 26-97-06,
e-mail: elmash@em.ispu.ru
t
