Алгоритм оптимизационных расчётов высоковольтных асинхронных двигателей Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 621.313.333.001.5

АЛГОРИТМ ОПТИМИЗАЦИОННЫХ РАСЧЁТОВ ВЫСОКОВОЛЬТНЫХ АСИНХРОННЫХ ДВИГАТЕЛЕЙ

В.С. Баклин, И.И. Пушкарёв*

Томский политехнический университет *ОАО «ТЭМЗ», г. Томск E-mail: [email protected]

Предлагается алгоритм расчёта и поиска оптимального варианта расчёта высоковольтных асинхронных двигателей с коротко-замкнутым ротором с пазами статора и ротора прямоугольной формы. Алгоритм содержит минимальное количество варьируемых переменных: индукцию в воздушном зазоре, ширину стандартного прямоугольного провода, число эффективных проводников в пазу статора и ширину паза ротора, что позволяет существенно уменьшить количество вариантов расчёта и для поиска оптимального варианта, удовлетворяющего функции цели, использовать метод полного перебора.

Ключевые слова:

Алгоритм, оптимизационные расчёты, асинхронный двигатель, статор, ротор, варьируемые переменные, ограничения, функция цели, метод полного перебора. Key words:

Algorithm, optimized calculations, induction motor, stator, rotor, varying variables, limitations, purpose function, method of complete excess.

Создание новых конкурентоспособных асинхронных двигателей предполагает наличие современной вычислительной техники и программного обеспечения, позволяющего автоматизировать расчёты и поиск оптимального варианта расчёта, удовлетворяющего функции цели. Создание такого программного обеспечения является актуальной задачей для многих предприятий электротехнической промышленности, связанных с производством асинхронных двигателей.

В программах расчёта и поиска оптимального варианта электрической машины исходные данные подразделяются на постоянные - неизменные в течение всего расчёта-поиска, переменные - изменяемые в течение всего расчёта-поиска и лими-тёры (ограничения). Варианты, не удовлетворяющие хотя бы одному из лимитирующих показателей, исключаются из расчётов. Остальные варианты расчётов подвергаются оценке на степень оптимальности, т. е. соответствию функции цели.

Для задач проектирования асинхронных двигателей наиболее характерны [1-3] следующие критерии оптимальности:

1. Минимум приведенной стоимости машины, р/год.

2. Минимум массы О - критерий, часто используемый при разработке транспортных машин.

3. Максимум КПД или минимум потерь Ьр.

4. Совокупный критерий, например, О+3р, где 3 - коэффициент.

Лимитирующими показателями [1, 2] могут быть:

1. Массогабаритные по наружному диаметру корпуса, по объёму машины, по массе машины.

2. Эксплуатационно-технические по превышению температуры обмотки статора А$м1 над температурой окружающей среды, по кратно-стям начального пускового момента Кп и максимального момента Км, по кратности пускового тока К,, скорости нарастания температуры обмотки статора при заторможенном роторе.

3. Конструкционно-технологические - по ширине зубцов статора и ротора, минимальным пазовым открытиям, высоте пазов ротора и толщине спинки статора, размерам проводников обмотки статора и т. д.

Например, из конструкционно-технологических ограничений должны быть: минимальная ширина зубцов статора и ротора Ьгтк>1+0,01ВП, мм; высота спинки статора Н]1>0,055Вп, где Бп - наружный диаметр сердечника статора. Полная высота паза ротора Нп2 должна быть: при отсутствии аксиальных вентиляционных каналов Нп2>0,05Вп-0,5(Б1-йъ)-5; при наличии аксиальных вентиляционных каналов Ап2>0,03Д,-0,5(Д-йУ-5-42, где Д - диаметр расточки статора, мм, йъ - внутренний диаметр сердечника ротора под вал, мм, ёк2 - диаметр аксиального канала, мм, 5 - воздушный зазор, мм.

В проектировании электрических машин термин «синтез» означает процесс нахождения, по меньшей мере, одного варианта удовлетворяющего всем требованиям, предъявляемым к машине. Отличительные особенности основных методов проектирования представлены в [1, 3]. К ним относятся: синтез машины, подобной изготовленному образцу; синтез на основании уравнения Арнольда; синтез в направлении от электромагнитных нагрузок к размерам и т. д.

Выбор функции цели, лимитирующих показателей, метода проектирования определяется конкретной задачей проектирования, её спецификой. Задачей проектирования являются высоковольтные асинхронные двигатели с пазами статора и ротора прямоугольной формы. Обмотка статора в таких двигателях выполняется катушечной из провода прямоугольного сечения, что придает обмотке достаточно высокую механическую прочность, а двигателю - надежность в эксплуатации. Пазы на роторе закрытые, короткозамкнутая обмотка ротора сварная из меди. Двигатели специализированные, предназначены для привода вентиляторов

главного и местного проветривания, должны хорошо встраиваться в вентиляторную установку, удовлетворять ограничениям по наружному диаметру и длине сердечника статора.

При разработке алгоритма и создания программного обеспечения оптимизационных расчетов высоковольтных асинхронных двигателей нами использованы элементы метода проектирования «синтез в направлении от электромагнитных нагрузок к размерам». Данный метод [1] предполагает использование в качестве входных данных семи электромагнитных нагрузок. К ним относятся: магнитные индукции ВВ - в воздушном зазоре, БгХ и В2 - в узких сечениях зубцов статора и ротора, В]1 и В,2 - в спинках статора и ротора; отношение линейной нагрузки А1 к плотности тока /1, обозначенное символом П=А/11 и произведение линейной нагрузки на плотность тока А=А111. Считается заданным также наружный диаметр сердечника статора.

Эти электромагнитные нагрузки [1] позволяют определить семь главных выходных величин: диаметр расточки статора Д, высоту спинки статора кп и ротора Ад, минимальную ширину зубца статора Ьй и ротора Ъй, расчетное число витков w1 фазы обмотки статора и расчетную длину /В1 сердечника статора. Предполагается, что остальные параметры, характеризующие структуру машины, либо зафиксированы, либо однозначно определяются главными выходными величинами.

При таком подходе к определению главных выходных величин, например, ширины зубца статора ЬгХ при заданной толщине пазовой изоляции в случае пазов статора прямоугольной формы не учитываются стандартные размеры прямоугольных проводов. Не учитывается, что при определении числа витков w1 число эффективных проводников ип1 в пазу статора должно быть целым, а в большинстве случаев еще и четным числом. Полученный оптимальный вариант без учёта этих факторов требует пересчёта и после пересчёта может оказаться далеко не оптимальным.

Суть, предлагаемого нами алгоритма программы расчёта и поиска оптимального варианта высоковольтного асинхронного двигателя с пазами статора и ротора прямоугольной формы заключается в следующем.

В качестве входных постоянных данных принимаются предельные значения электромагнитных нагрузок: магнитные индукции в узком сечении зубца статора ВЙ1ШХ и ротора Вйшш в спинке статора В,1тяу и ротора В,2тяу; линейная нагрузка А1тах; произведение линейной нагрузки на плотность тока А/гаах. Предельные значения электромагнитных нагрузок могут быть приняты на основании рекомендаций [2] из опыта проектирования машин подобного типа.

Линейная нагрузка А1тах и произведение линейной нагрузки на плотность тока А/гаах позволяют определить допустимую плотность тока в обмотке статора

Т = АТ / А

1тах тах 1тах

и допустимое отношение линейной нагрузки к плотности тока

^ тах А1тах / АТтах *

Допустимые значения плотности тока /11ШХ и отношения иш1Х используются в программе расчёта-поиска как лимитирующие показатели.

Алгоритмом программы предусмотрен ввод предварительных значений КПД т]х и коэффициента мощности ео8фх. В задачах проектирования, где критерием оптимальности является минимум потерь в двигателе, значения этих коэффициентов принимаются не ниже, чем у двигателей-аналогов. Значение КПД т]х используется как лимитирующий показатель. В программе предусмотрена также таблица стандартных размеров прямоугольных проводов, из которой проводится автоматизированный выбор толщины прямоугольного провода по заданной его ширине.

Расчётная длина сердечника статора /я принята фиксированной, а при необходимости может быть использована как независимая варьируемая переменная.

В качестве независимых переменных приняты: индукция в воздушном зазоре В5, стандартные размеры ширины прямоугольного провода Ьм обмотки статора, число эффективных проводников ип1 в пазу статора, ширина паза Ьп2 ротора.

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

Варьирование независимых переменных целесообразно начинать с их верхнего уровня. При принятых лимитирующих показателях это позволяет существенно уменьшить количество вариантов расчёта.

Алгоритм варьирования переменными следующий. Задаются максимальное В8ша и минимальное В^д значения индукции в воздушном зазоре. Изменение индукции в зазоре Вад от своего максимального до минимального значения проводится в цикле £=1,...,« с шагом АВЗ=(В5пш-В5пЛпк.

Рассчитываются вспомогательные коэффициенты

С = 1 -

^1(4) 1

В

В( к)

В^тах кс1

И С 1(4 ) = 1 +

В

В (к)

В 1тах Ркс 1

где кс1 - коэффициент заполнения сердечника статора железом; р - число пар полюсов двигателя.

Из таблицы стандартных размеров прямоугольных проводов, представленных в программе, в цикле /=1,...,«; задается ширина провода Ьм(1) от максимального до минимального значения и рассчитывается соответствующая ему ширина паза статора по формуле

кт _ (Ъш{1)+А4и + Арг + \ )• 1(Г3, (1)

где АЬп - двухсторонняя толщина пазовой изоляции с учетом разбухания и припуска на укладку обмотки; Арг - двухсторонняя толщина изоляции провода; А5( - припуск на сборку сердечников статора и ротора. Минимальная ширина проводника Ьш{)) должна быть такой, чтобы ширина паза Ьп10, рассчитанная по формуле (1), из конструкционно-технологических ограничений была больше 7...8мм [3].

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

_ паК кес*т

п1( к,;)

Ва( к )Ь,

(к)" „1(1)

и.

нф

В

'&(к)

*ао = 1 - сг10,к) К'

диаметр расточки статора

Л"„1(0

пС.

*1( 1 )

шаг по пазам

*1( и)

и минимальная ширина зубца статора

*1( 1)

_ г

1( 1г)

п 1(1)

Если В^В^ или Ьт)<0,001+0,01Вп (конструкционно-технологические ограничения), то вариант расчета с числом проводников ип1(1.0 исключается, и проводится расчёт с уменьшенным по формуле (3) числом проводников в пазу.

При выполнении условий Вти,))<Вт1тях и Ьг1(,)>0,001+0,01^„ рассчитываются предварительно полная высота прямоугольного паза статора

К

1( 1

_ 0,5 Д. -

0,521 С

11( к)

пС

1( 1,0

„1(0'

где а1 - число параллельных ветвей; д1 - число пазов на полюс и фазу обмотки статора; 2Х - число пазов статора; ке - коэффициент отношения ЭДС к номинальному фазному напряжению; ки - коэффициент,

и высота паза занятого обмоткой, мм

КиШ) _ К„1(1,0103 -АКп ,

где АН„ - высота паза статора не занятая обмоткой, включающая в себя толщину изоляции по высоте паза, толщину клина, высоту шлица, припуски на сборку сердечника и укладку обмотки статора, мм.

Предварительно определяется толщина эффективного проводника, мм

к _-

и

где инф - фазное напряжение; км - обмоточный коэффициент.

В зависимости от требований, предъявляемых к обмотке статора, число эффективных проводников и'„1(к)), рассчитанное по формуле (2), программно округляются в меньшую сторону до ип1(к)) целого числа или до целого четного числа.

Числа эффективных проводников в пазу статора изменяется в цикле у=1,...,«,-. При номере шага у=1 проводится присвоение и^^и^. При номерах шага у>1 число эффективных проводников уменьшается на величину Аип по формуле

ип1Ш) _ ипЮ-1,0 -Аип , (3)

где Аип - шаг уменьшения числа проводников, причем принимается Аип=1 при нечетном и Аип=2 при четном числе эффективных проводников в пазу.

Для каждой ширины паза и соответствующего ей числу эффективных проводников в пазу уточняется вспомогательный коэффициент

с _ и Вч к Ако

1,0 "п1( ],() г г ■

па1ки ке

Рассчитываются индукция в зубцах статора В,

Кпи( 1, г)

•*и1(/,о

-А„

(4)

п1( 1, г)

где Апр - толщина двухсторонней изоляции провода.

Минимальная толщина стандартного прямоугольного провода 0,8 мм [4]. Двухсторонняя толщина изоляции прямоугольного провода обмотки статора в высоковольтных асинхронных двигателях составляет порядка 0,45 мм. С уменьшением толщины провода возрастает влияние толщины изоляции провода на коэффициент заполнения паза медью, на плотность тока. Чем меньше коэффициент заполнения паза медью, тем меньше КПД двигателя, поэтому для обмоток статора высоковольтных асинхронных двигателей целесообразно ограничить минимальную толщину прямоугольного провода. С учётом стандартных размеров прямоугольных проводов в программе минимальная толщина провода принята 1,18 мм.

При толщине провода а'„у<1,18 мм вариант расчёта с шириной паза ЬпЩ) исключается и осуществляется переход к расчётам с последующей шириной паза Ьп1т.

Если а'мЫ>1Д8 мм, то соответственно ширине провода Ьщ из таблицы стандартных размеров прямоугольных проводов программно выбирается стандартная толщина провода амМ, ближайшая (в меньшую сторону) предварительно определенной по формуле (4) толщине а „(,,) провода.

Уточняются высота паза статора, м,

■п К ю г) +АК ) -10-3,

КшШ) _ (К,г)

толщина спинки статора

_ 0,5(Дп - Дюг) -2К 1,1 ^

коэффициент

С

11( 1 ,г)

_ 1 +

пС^1( 1 ,0 1 1.0 0,521Ьп1(г)

Рассчитывается индукция в спинке статора

B

B

5 (к)

ркЛСм,л -1)

Варианты с индукцией БущУБ^ и толщиной спинки статора кт<0,055Бп. Исключаются из расчётов.

Рассчитываются сечение эффективного проводника, м2,

Чи ш) = К , До ))10-6, сечение меди в пазу статора

Чип (] ,;) = Чи (] ^) ип1(_] ^)

и коэффициент заполнения паза статора медью

Чип (] ,1)

k

3 м ( j ,i)

1 А

n(i/'ni( j ,i )

Рассчитывается отношение линейной нагрузки к плотности тока

UAJ(j,i) 0,5k3M(ji) Cz i(j i)

D -

ZXC

ji( j ,•)

nC,

ili)

i( j ,i)

Лимитирующим показателем допустимых значений отношения линейной нагрузки к плотности тока является ипш=А21шш/А]Пш. Вариант Ц/адУ/шх из дальнейших расчётов исключается.

Исключается также вариант, в котором плотность тока в обмотке статора превышает допустимую плотность тока /1тях. С этой целью рассчитывается линейная нагрузка

Ai(j,o ="

npPe

(

C.

\2

z1( j i )

Zxt

"1(i)

(5)

2,22 fikwihiBs(k) В формуле (5) расчетная мощность двигателя

= кр„io3 ,

e Г1Х cos px '

где Рн - номинальная мощность двигателя, кВт.

Определяется ожидаемая плотность тока в обмотке статора

Ji( j,i)

A t

j,i)1 zi( j,i)

h

ЛШ) = 0,03Dn -0,5(Di(ji) - d2)-(5 + d2)i0-

где й2 - внутренний диаметр сердечника ротора под вал, м.

Для определения пазовой геометрии ротора в качестве независимой переменной в цикле /2=1,.. ,,и2 задается ширина прямоугольного паза ротора Ьод, м.

В первом приближении индукция в узком сечении зубца ротора принимается Бде)=Бг2шях. Рассчитывается коэффициент

C

B

5( k)

5

z 2(1- 2)

Bz 2(i2)kc 2 l2

и полная высота паза ротора

(

hn 2(i 2) 0'5 D2( j,i)

2

2n

C

2(i2)

n2(i2)

C

ni(i)

zi( j,i )

(6)

(7)

где 12 - длина сердечника ротора.

Если А„2(й)УА„2(д ) (предельной высоты паза ротора), то индукция Бт) в зубце ротора уменьшается на величину АБЙ=0,01 Тл, и проводится пересчет по формулам (6), (7) до выполнения условия

^п2(й)—^п2а,1).

Определяются высота стержня, м,

Кю 2) = Кп 2(/ 2) - + А 2)10-3>

его ширина

и сечение

"c 2(i 2)

qc2(i 2)

n 2(i 2) ^st 2i0

= h

c 2(i 2) <c 2(i 2) '

Чип (] ,1)

Вариант, в котором плотность тока в обмотке статора ^1(,-,-)У^1шах или линейная нагрузка АщУА11ШХ исключается из расчетов.

Независимые переменные - индукция в воздушном зазоре Б, стандартные размеры ширины Ъм прямоугольного провода обмотки статора, число эффективных проводников ип1 в пазу статора полностью определяют внутренний диаметр расточки, пазовую геометрию и обмоточные данные статора.

Для определения размеров и пазовой геометрии ротора рассчитывается наружный диаметр сердечника ротора

Аао = »ю;) - 25'10-3-

С учётом конструкционно-технологических ограничений рассчитывается допустимая высота паза ротора, м,

где НА - высота шлица в случае полузакрытых или высота ферромагнитного мостика в случае закрытых пазов ротора; А5Й - припуск на сборку сердечника и сварной короткозамкнутой обмотки ротора.

Дальнейшие расчёты параметров, магнитной цепи, рабочих и пусковых характеристик, тепловой расчет двигателя проводятся по известным методикам [2, 4]. В ходе расчётов проверяются эксплуатационно-технические лимитирующие показатели: кратности пускового и максимального моментов, кратность пускового тока, превышение температуры нагрева обмотки статора, КПД двигателя, значение которого должно быть не меньше чем т]х. Варианты, не удовлетворяющие лимитирующим показателям, из расчётов исключаются.

КПД является общей оценкой двигателя наиболее близкой к экономическому критерию и в частности к энергоэффективности двигателя. Повышение КПД асинхронных двигателях для привода вентилятора при ограниченных размерах наружного диаметра и длины сердечника статора может достигаться снижением переменных потерь за счёт более напряженных магнитных нагрузок в железе магнитопро-вода. Но это может сопровождаться не обоснованным увеличением массы меди обмотки статора и снижением коэффициента мощности двигателя.

Критерием оптимальности выбран обобщенный параметр, включающий в себя КПД двигателя т]х, значение которого должно быть не ниже КПД

двигателя-аналога при минимальной массе меди обмотки статора проектируемого двигателя. Так как КПД двигателя т]х используется в качестве лимитирующего показателя, то оптимальным считается вариант расчёта с минимальной массой меди обмотки статора.

В программе вариант расчёта, полностью удовлетворяющий лимитирующим показателям, запоминается. Если в полученном варианте КПД двигателя больше значения ц„, а плотность тока в обмотке статора меньше допустимой плотности тока /1гаах, то из таблицы стандартных размеров провода берётся толщина провода а^ на ступень ниже, чем в предыдущем варианте. Проводится расчёт с уменьшенной толщиной провода. Если вариант с уменьшенной толщиной провода удовлетворяет лимитирующим показателям, то из двух вариантов расчёта предыдущего и последующего запоминается вариант с минимальной массой обмотки статора. Уменьшение толщины провода проводится до тех пор, пока КПД двигателя больше значения ц„, а плотность тока в обмотке статора меньше допустимой плотности /11ШХ.

В поисках оптимального варианта расчёта, удовлетворяющего функции цели - обеспечение КПД проектируемого двигателя не ниже КПД двигателя-

СПИСОК ЛИТЕРАТУРЫ

1. Домбровский В.В., Зайчик В.М. Асинхронные машины: Теория, расчёт, элементы проектирования. - Л.: Энергоатомиздат, 1990. - 368 с.

2. Копылов И.П., Клоков Б.К., Морозкин В.П. и др. Проектирование электрических машин. - М.: Высшая школа, 2005. -767 с.

аналога при минимальной массе меди обмотки статора, используется метод полного перебора [1, 3].

Выводы

1. Разработан алгоритм расчёта и поиска оптимального варианта расчёта высоковольтного асинхронного двигателя для привода вентилятора.

2. Выбраны варьируемые переменные, которые позволяют найти оптимальный вариант расчёта двигателя с учётом целого числа эффективных проводников в пазу статора и стандартных размеров обмоточного провода.

3. Выбранный критерий оптимальности - КПД двигателя - не меньше, чем КПД двигателя-аналога при минимальной массе меди обмотки статора, позволяет спроектировать двигатель более энергоэффективным по сравнению с двигателем-аналогом, что повышает конкурентоспособность вентиляторной установки.

4. Принятые эксплуатационно-технические и конструкционно-технологические лимитирующие показатели, малое число независимых переменных и их варьирование, начиная с верхнего уровня, позволяют существенно уменьшить количество вариантов расчёта и использовать для поиска оптимального варианта метод полного перебора.

3. Бочаров В.И., Захаров В.И., Коломейцев Л.Ф. и др. Тяговые двигатели электровозов. - Новочеркасск: Агентство Наутилус, 1998. - 672 с.

4. Гурин Я.С., Кузнецов Б.И. Проектирование серий электрических машин. - М.: Энергия, 1978. - 480 с.

Поступила 08.11.2010 г.