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

ПРОБЛЕМЫ ЭНЕРГЕТИЧЕСКОГО

МАШИНОСТРОЕНИЯ

УДК 621.313.322

ЭЛЕКТРОМАГНИТНАЯ МОДЕЛЬ И ОПТИМИЗАЦИЯ ПАРАМЕТРОВ ТОРЦЕВОГО ГЕНЕРАТОРА

© 2012 А.Л. Встовский, М.П. Головин, К.С. Федий, Н.А. Колбасина, Д.И. Морозов

Сибирский федеральный университет, г. Красноярск

Поступила в редакцию 16.03.2012

В статье изложены особенности расчетной модели активного объема низкоскоростного торцевого синхронного генератора с постоянными магнитами и программного комплекса синтеза и многокритериальной оптимизации геометрических параметров, как компонента интегрированной среды проектирования свободнопоточных микроГЭС.

Ключевые слова: торцевой генератор, микроГЭС, внешняя характеристика

Развитие малой гидроэнергетики, возобновляемых источников энергии потребовало разработки низкоскоростных герметичных генераторов, способных работать в погруженном состоянии. Эту проблему решает низкоскоростной торцевой синхронный генератор (НТСГ), имеющий высокотехнологичную конструкцию и низкую материалоемкость. В Политехническом институте Сибирского федерального университета по проекту Минобрнауки РФ (Госконтракт № 02.740.11.0056) разработан торцевой синхронный генератор, в котором в качестве обмоточной структуры использован активный распределенный слой статора, предложенный В.М. Казанским в 1970 г., а возбуждение осуществляется от поперечно намагниченных постоянных магнитов [1]. Особенности этого типа генератора позволяют его реализовать в низкоскоростном герметичном исполнении, необходимом для свобод-нопоточной микроГЭС. Результаты испытаний первых опытных образцов микроГЭС с НТСГ показали их работоспособность, соответствие условиям эксплуатации. В то же время некоторые из основных эксплуатационных параметров могут быть повышены за счет уточнения расчетных моделей при более эффективном выборе электромагнитных и геометрических параметров.

Встовский Алексей Львович, кандидат технических наук, профессор. E-mail: [email protected] Головин Михаил Петрович, кандидат технических наук, профессор. E-mail: [email protected] Федий Константин Сергеевич, кандидат технических наук, доцент. E-mail: [email protected] Колбасина Наталья Анатольевна, кандидат технических наук, доцент. E-mail: [email protected] Морозов Денис Иванович, старший преподаватель. Email: [email protected]

Примеры повышения эффективности НТСГ: увеличение сечения магнита при оптимальном его размещении в ограниченном пространстве; снижение пазового рассеяния.

Необходимость повышения единичной мощности микроГЭС (от 1-3 кВт до 20-30 кВт) и возможность использования на крутосклонных реках со скоростью течения до 3,5-5 м/с потребовали существенных изменений параметров:

- использование провода не только круглого, но и прямоугольного сечения вследствие повышения токов обмотки статора;

- увеличение высоты паза, приводящее к возрастанию пазового рассеяния, и поиску компромисса между радиальным габаритом генератора, массой магнитов и внешней характеристикой генератора;

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

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

модели, представлена на рис. 1.Выражение для составляющей вектора магнитной индукции в воздушном зазоре (зона 3) расчетной модели, полученная из решений уравнений магнитостатики [4], имеет вид

B2, =*>§ = №(Mk = С

^ch k^S-sin1^-x | =

0 sin1^-x

2y t

(1)

где Мк - постоянные интегрирования, определяемые из условия непрерывности потенциала; т - полюсное деление по среднему диаметру; 5 -величина воздушного зазора.

^_;_;_;_>л_;_;_;__

■ у ц = ............'...........|.................................Остов ротора .)...

Конструктивный зазор

~1 , г

Зубцовый с^ой статора

+j Воздушный зазор

Е '

Рис. 1. Расчетная модель активной части генератора

Поле возбуждения В2у, созданное магнитами ротора, и поле реакции якоря В3у относительно неподвижны. Пространственно эти поля сдвинуты друг относительно друга на внутренний угол нагрузки ©. Тогда поле реакции якоря этой же поверхности будет описываться выражением

К (х,у) = !Сзу + , (2)

где С3у - постоянная интегрирования являющаяся функцией настила тока:

J = ■

ипсл!ф

t

1

1 ф.

j - ю - Шф - B 2 y - /5 - т

+ j - хл + j - a- Шф -15

■т-«псл -(y - X + Z )-

гН

(3)

где Y,X,Z - функции внутреннего угла нагрузки © и геометрических параметров генератора:

воздушного зазора, полюсного деления, размеров магнита. Так как числитель представляет собой ЭДС магнитного поля ротора, то синхронное индуктивное сопротивление равно:

= j-ф-W - ls-r-uncn-(Y - X + Z)+j -

(4)

то есть представляет собой сумму индуктивного сопротивления реакции якоря и пазового рассеяния.

Основной характеристикой синхронного генератора является внешней характеристики -зависимость выходного напряжения генератора от тока якоря при постоянной угловой скорости и характере нагрузки U=f(I) (при Q=const, cos 9=const). Тогда используя выражение для фазного тока (3) можно записать:

U = -j-a-W^- Взу- S

1ф-rj + 1ф- j- xл + 1ф-j-a-Wф-lg-х-ипсл -(Y-X + Z)j (5)

Электромеханическая мощность, развиваемая ротором синхронной машины с беспазовым статором, равна приращению действительной части потока вектора Пойнтинга при переходе через границу ротор - воздушный зазор, с учетом (1) и (2) выражение электромеханической мощности примет вид:

b * P = 2p/suRe ]>омпсл-(y -X + Z)sin(©) -Hdx; (6) о ф

где uпсл - число проводников в пазу статора, Ь -зубцовый шаг статора.

Полученное из уравнения, составленного по второму закону Кирхгофа, значение фазного тока статора с учетом (1) и (2) после преобразований имеет вид

Используя это выражение, можно аналитичес-ки рассчитать угловую характеристику мощности PH при различных значениях угла © проанализировать влияние геометрических параметров на перегрузочную способность генератора. Расчетная модель генератора позволяет сформировать поток данных, необходимых для его проектирования и конструирования:

- энергетические характеристики (мощность, КПД, напряжение, ток статорной обмотки);

- массогабаритные показатели, необходимые для производителя и потребителя машины (массы электротехнической стали, обмоточной меди, магнитов, общая масса машины, ее габаритные размеры);

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

- разработка целевой функции качества системы;

3

x

Остов ротора

- формулирование критериальных ограничений;

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

В программном модуле предоставлена возможность комбинировать в любых сочетаниях методы оптимизации, предварительно определив:

- непрерывные диапазоны их изменений

min max

a| < а| < a| ,

(11)

где а1 - варьируемый параметр; а;т1П, а;тах -начало и конец диапазона; - необходимое количество функциональных ограничений в виде:

&(xk) ^ 0

(12)

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

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

где g(x0,...,x) - функция нескольких переменных.

- критериальные ограничения системы (критерии качества), их минимальные и/или максимальные значения в виде:

k < f (Хо .-Х) ^ MIN, k < f (x0 ...xn) ^MAX , (13)

где f (x . •••Хи) - функция заданная аналитическим выражением, k - максимальное (минимальное) значения критерия,.

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

Рис. 1. Основные элементы интерфейса программы КТ80еп

Количество точек ЛП-последовательности определяется по формуле:

c = ■

L

(14)

где Lmax - длина наибольшего диапазона значения параметра; p - шаг разбиения наибольшего диапазона, определяющий общее количество точек, а значит точность решения задачи.

Общее количество точек:

k = ф div 1)n

(15)

где п - количество параметров; с - количество точек (14).

Количество отрезков в сетке для одного измерения определяется так:

d =

k

c div 1

(16)

где k - общее количество точек (15); п - количество параметров; с - количество точек ЛП-последовательности (14). Значение параметра для одного решения:

, Ъ — a

k = a л---z,

d

(17)

где а - начало диапазона, Ь - конец диапазона, d - число определяется по (16), г - случайное число, обеспечивающее уникальных значений параметров.

Рис. 3. Ввод исходных данных и граничных условий

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

При определении коэффициента оптимальности точки по критерию, который должен стремится к максимуму или минимуму, используются формулы:

R = D

MAX, — Ct ' MAX, — MIN,

R = D •-

C — MIN,. MAX. — MIN. (18)

В (19) и (20) МАХ1, МШ1 - максимальное или минимальное значение /-го критерия, С, -значение /-го критерия точки, О, - значение

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

Таблица 1. Выходные параметры генератора

Выходные параметры Много-крите-риаль-ная оптимиза- Одно-нофак факторный анализ % откло кло-нения

ция

индукция в воздушном зазоре, Тл 0,68 0,7 -2,9

активный диаметр, м 0,632 0,67 -6,01

активная длина, м 0,02 0,021 -5

напряжение, В 238,32 218,21 8,44

максимальная мощность, Вт 6451 5864 9,1

масса генератора, кг 36,235 40,26 -11,1

масса магнитов, кг 1,226 1,326 -8,16

индуктивное сопротивление, Ом 4,245 5,224 -23,06

В качестве примера, иллюстрирующего возможности предлагаемых моделей и методов, использован генератор мощностью 3 кВт с числом полюсов 24. Здесь метод многокритериальной оптимизации позволил значительно улучшить выходные параметры генератора по сравнению с использованной ранее методикой одно-факторного анализа (табл. 1).

Применение метода многокритериальной оптимизации позволило снизить (знак минус) массу дорогостоящих магнитов на 7,6% и суммарную массу генератора на 7,4%, а также увеличить (знак плюс) энергетические показатели машины на 8,4% за счет снижения влияния сопротивления пазового рассеяния. Алгоритмы, разработанные на основе методик электромагнитного расчета и многокритериальной оптимизации, реализованы в виде программного комплекса, который позволяет автоматизировать твердотельное моделирование в САВ-среде, сокращая сроки разработки изделия от технического задания до комплекта рабочей документации. В состав программного компонента интегрированной среды проектирования свободнопоточ-ной микроГЭС, кроме описанных в статье алгоритмов, будут на следующем этапе включены модули анализа отказоустойчивости наиболее нагруженных элементов генератора: жесткость диска остова статора и ротора, сопряженного с валом, долговечность подшипников ротора и др.

Выводы: созданный коллективом авторов инструментарий автоматизированной разработки микроГЭС, опыт их проектирования и испытаний позволили перейти к сертификации изделий и технологической подготовки промышленного производства.

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

1. Патент 2246168 РФ, МПК7 Н 02 К 21/24. Торцевая электрическая машина / А.Л. Встовский, М.П. Головин и др.; опубл. 10.02.05, Бюл. № 2.

2. Головин, М.П. Автоматизация проектирования свободнопоточной микроГЭС / М.П. Головин, А.Л. Встовский и др. // Вестн. Краснояр. гос. техн. ун-т. Машиностроение. - Красноярск, 2005. Вып.40. С. 89-101.

3. Инкин, А.И. Электромагнитные поля и параметры электрических машин: учеб. пособие. - Новосибирск: ЮКЭА, 2002. 464 с.

4. Соболь, И.М. Выбор оптимальных параметров в задачах со многими критериями / И.М. Соболь, Р.Б. Статников. - М.: Наука, 1981. 328 с.

ELECTROMAGNETIC MODEL AND OPTIMIZATION PARAMETERS OF THE FACE GENERATOR

© 2012 A.L. Vstovsky, L.S. Golovin, K.S. Fediy, N.A. Kolbasina, D.I. Morozov

Siberian Federal University, Krasnoyarsk

In article features of calculation model of the active volume of the low-speed face synchronous generator with permanent magnets and a program complex of synthesis and a multicriteria optimization of geometrical parameters, as a component of the integrated medium of projection the free-stream microhydroelectric power stations are explained.

Key words: face generator, microhydroelectric power station, external characteristic

Aleksey Vstovskiy, Candidate of Technical Sciences, Professor. E-mail: val [email protected]

Mikhail Golovin, Candidate of technical Sciences, Professor. E-mail: [email protected]

Konstantin Fediy, Candidate of Technical Sciences, Associate Professor. E-mail: [email protected]

Nataliya Kosbasina, Candidate of technical Sciences, Associate Professor. E-mail: [email protected]

Denis Morozov, Senior Teacher. E-mail: [email protected]