CC BY
122
621.313:681.3
doi: 10.20998/2074-272X.2016.1.05
В.И. Милых, Н.В. Полякова
ОПРЕДЕЛЕНИЕ ЭЛЕКТРОМАГНИТНЫХ ПАРАМЕТРОВ И ФАЗОВЫХ СООТНОШЕНИЙ В ТУРБОГЕНЕРАТОРАХ АВТОМАТИЗИРОВАННЫМ РАСЧЕТОМ МАГНИТНОГО ПОЛЯ В ПРОГРАММНОЙ СРЕДЕ FEMM
Надані теоретичні основи розрахунку електромагнітних величин і їх фазових співвідношень для турбогенераторів. Це реалізується шляхом числових розрахунків магнітного поля у програмному середовищі пакету FEMM. Розрахунки автоматизовані за допомогою створеної програми на алгоритмічній мові Lua Вона управляє розрахунками і організовує видачу результатів в текстовий файл. Представлений скрипт Lua є доволі універсальним і може послужити основою щодо створення аналогічного програмного забезпечення для інших типів електричних машин. Бібл. 9, рис. 6.
Ключові слова: програма FEMM, турбогенератор, електромагнітні параметри, фазові співвідношення, автоматизовані розрахунки, Lua скрипт.
Представлены теоретические основы расчета электромагнитных величин и их фазовых соотношений для турбогенераторов. Это реализуется путем численных расчетов магнитного поля в программном пакете FEMM. Расчеты автоматизированы посредством созданной программы на алгоритмическом языке Lua. Она управляет расчетами и организует выдачу результатов в текстовый файл. Представленный скрипт Lua универсален и может послужить прототипом аналогичного программного обеспечения для других типов электрических машин. Библ. 9, рис. 6. Ключевые слова: программа FEMM, турбогенератор, электромагнитные параметры, фазовые соотношения, автоматизированные расчеты, Lua скрипт.
Введение. Широко известная программа FEMM [1] эффективна для расчетов двухмерных магнитных полей (МП) электрических машин (ЭМ) классической конструкции при практически плоскопараллельном МП в пределах их активной части. Например, это относится к турбогенераторам (ТГ) [2]. Эта программа освободила расчетчиков от разработки реализуемых математических моделей и программирования и обладает достаточно удобным интерфейсом.
Дополнительным средством облегчения труда пользователей программы FEMM являются скрипты на алгоритмическом языке Lua, интегрированном в нее. Так, например, в [3] описывается наша программа, автоматизировавшая построение расчетных моделей ТГ, в [4] дана программа автоматизированных расчетов динамики их электромагнитных процессов.
Эффективность скриптов Lua заключается в том, что программы пишутся один раз, а используются бесчисленное число раз и любым числом пользователей.
Данная статья продолжает серию работ, начатых в [3, 4]. И ее целью является представление скрипта Lua, автоматизирующего получение ряда электромагнитных параметров и фазовых соотношений в ЭМ посредством численных расчетов МП в программной среде FEMM. Это делается на примере ТГ - одного из крупнейших и ответственных их представителей [5].
Объект исследования. Для иллюстраций здесь, как и в [3, 4], взяли ТГ мощностью 340 МВт. Его параметры есть далее в файле исходных данных.
Электромагнитная система ТГ представлена на рис. 1 его поперечным сечением. Выделены фазные зоны двухслойной укороченной трехфазной обмотки статора А—А, В—В и С-С'. Показаны используемые прямоугольная (x, у) и полярная (r, а) системы координат, продольная d и поперечная q оси ротора.
Основой разработанного скрипта Lua является заранее подготовленная для среды FEMM расчетная модель электромагнитной системы ТГ. У нас это де-
лается посредством скрипта Lua, представленного в [3], а при отсутствии такого - в «ручном» режиме в соответствии с инструкциями к программе FEMM.
Рис. 1. Модель электромагнитной системы ТГ в его поперечном сечении
Вид вызванной на экран компьютера в среде FEMM (тип файла - Femme Document) расчетной модели ТГ представлен в целом на рис.2, а на рис. 3 выделен ее фрагмент. Видны метки физических свойств блоков - ограниченных участков области расчета. Метки Fe1, Fe2 и Fe3 задают кривые намагничивания ферромагнетиков, m0 - магнитную проницаемость |д0=4-л-10-7 Гн/м, метки iA, iB, iC - фазные токи, Ir -ток ротора. Добавление + или - к метке тока говорит о его направлении: если ток в данный момент времени положительный, то в данном месте он берется с соответствующим знаком. Цифра после двоеточия обозначает число эффективных проводников в данном блоке. Окружности вокруг меток показывают
© В.И. Милых, Н.В. Полякова
26
ISSN 2074-272X. Електротехніка і Електромеханіка. 2016. №1
максимально допустимый размер треугольников в структуре конечных элементов, формируемой программой FEMM.
Фрагмент на рис. 3 взят с правой половины полной модели по рис. 2, в остальной части все сформировано аналогично в соответствии и моделью на рис.1. В частности, на левой половине знаки токов противоположны тому, что видим на рис. 3.
Рис. 2. Исходная физико-геометрическая модель ТГ
Рис. 3. Фрагмент физико-геометрической модели ТГ
Упомянутые блоки в FEMM-модели объединены в группы для оперирования с ними, как с единицами области расчета. Из назначенных номеров групп [3] здесь используются: 2, 3, 4, 5, 6, 7 - фазные зоны А, В, С, А, В, С обмотки статора, 9 - весь ротор (рис.1).
Теоретические основы. Расчет стационарного МП программой FEMM проводится на основе численного решения дифференциального уравнения:
=kJz, (1)
где Az, Jz - аксиальные составляющие векторного магнитного потенциала (ВМП) и плотности тока; р -абсолютная магнитная проницаемость; к - орт.
rot
-rot
(k Az)
1
M-
На внешней границе области расчета для ВМП принято известное условие Дирихле Az = 0.
В модели ТГ (рис. 2) для расчета мгновенного распределения МП конкретного установившегося режима надо иметь соответствующие токи. В обобщающем режиме нагрузки (РН) [9] задается постоянный ток ротора Ir и мгновенные значения симметричной системы фазных токов статора: iA=Imcos(P);
iB=ImCOs(P-2rc/3); z'c=ImCOs(P+2rc/3), где Im=42 Is - их амплитуда и Is - действующее значение; р - их начальная фаза [6-9], задающая смещение МДС обмотки статора Fs от МДС обмотки ротора Ff (рис. 1). Здесь же дана условная результирующая МДС для МП в РН
Fi = Ff+Fs. (2)
В режиме холостого хода (ХХ) задается только ток ротора, в режиме установившегося короткого замыкания - только токи статора. Направления токов (рис. 1) регулируются по системе, принятой в [7].
Для ТГ результатом является распределение МП в его поперечном сечении в виде функции Az(x,y). Картина МП в области расчета представляется структурой силовых линий - линий равного ВМП Az=const. Пример этого для ТГ в режиме его номинальной нагрузки (НН) есть на рис. 1.
Одной из ключевых величин при анализе электро -магнитных параметров ТГ является магнитное потокосцепление (МПС) [2] и оно находится по распределению ВМП. Так, например, для любой из шести фазных зон (рис. 1) МПС определяется на основании
4= jAzdS * Nk£ Az,av, jASj , (3)
йф S S9 j=1
°ф J
где S<p - площадь сечения по несущим ток элементам фазной зоны; Кф - число элементов ее дискретизации; Az,avj - среднее значение ВМП в j-том элементе.
В формуле (3) учитывается активная длина ТГ la и число витков фазной обмотки статора Ns, а для определения S<p и интеграла в скрипте Lua есть соответствующие процедуры [1].
Исходной основой для выявления фазовых соотношений в любом из расчетных режимов является временная фаза или направление вектора МПС фазной обмотки А-А. МДС этой базовой обмотки на рис.1 ориентирована по продольной оси ротора d.
Направление же вектора МПС фазной обмотки А-А определяется структурой МП конкретного расчетного режима. Так на рис.1 это направление проявляется ориентацией силовых линий МП и в целом соответствует направлению вектора МДС Fl (2). Но пока этот вектор построен априори с учетом тех результатов, которые еще будут получены далее на основе методики, которая опубликована в [4-6]. Здесь ее изложим кратко для понимания смысла разработанного и представленного далее скрипта Lua.
Упомянутую методику изложим применительно к режиму ХХ, а затем она будет иметь силу для любого другого режима. Конкретно, будем считать, что МП в режиме ХХ рассчитано программой FEMM. Структура этого поля показана на рис. 4.
ISSN 2074-272X. Електротехніка і Електромеханіка. 2016. №1
27
на фоне рассчитанного МП для сбора значений МПС
В исходной угловой позиции aj фазной зоны А на рис.4,а на основании формулы (3) определяется МПС Теперь переместим условно эту зону на зубцовое деление в позицию a2 (рис. 4,6) и найдем соответствующее значение МПС ¥2. То же повторим в позиции a3 (рис. 4,е) и т.д. - вплоть до позиции a30 (рис. 4,г). В общем же случае число позиций равно числу пазов статора Qs. А каждое смещение фазной зоны делается на угол Aa=360/Qs (в градусах).
В итоге получится дискретная угловая функция МПС одной фазной зоны на ее периоде Т, т.е.
(ak); ak = (k-1) -Aa; k = J, 2, 3,... Qs, (4)
где k - номера угловых позиций этой зоны.
Каждой позиции фазной зоны в пределах от J до QJ2 есть диаметрально противоположная позиция в пределах от QJ2+J до Qs. Они попарно создают условно перемещающуюся фазную обмотку А-А, для которой угловая функция МПС на ее полупериоде:
^A,k = ¥ -Vk+Qs/2; к = J,2,3,..Q/2. (5)
По полученным результатам, с учетом периодичности величин в ТГ типа ¥(a + T /2) =-¥(a), дополним полупериод МПС (5) до полного периода:
^A,k+Qs/2 =-4A,k; k = 1,2,3,...Qs/2. (6)
Числовой ряд МПС (5) проиллюстрирован на рис. 5 совокупностью выделенных точек на кривой ХХ в диапазоне k от J до J5. На основании (6) кривая дополнена до полного периода (точки от j6 до 30), а точка 3j повторила первую (точка j) через период Т.
Аналогичные действия были проведены и после расчета МП в режиме НН, а результаты представлены кривой НН на том же рис. 5. Для проявления начальных фаз функций МПС кривые продолжены влево.
Рис. 5. Угловые зависимости МПС фазной обмотки ТГ
Полученную численную полную угловую функцию ¥Ak(a^), k = J,2,3,...Qs, можно разложить в гармонический ряд [6, 7]. Ввиду характерного для ТГ условия периодичности ¥(a + T /2) = -¥(a) в ряду нет постоянной составляющей и четных гармоник. Принятый в [6] косинусный ряд для МПС имеет вид
¥A = Yj¥m,vcos(7)
v=1,3,5,...
где амплитуды и аргументы (начальные фазы)
¥m,v = д/SV + c2; Yv = -arctg ^ (8)
cv
находятся по синусно-косинусным коэффициентам:
2 Qs 2 Qs
sv = Q Z¥A,k sin(vak ) ; cv = Q Z¥A,k cos(vak ) .
Qs k=1 Qs k=1
В (8) функция арктангенса arctg(sv/cv) должна быть расширенной с получением значения в диапазоне от -180° до +180° в зависимости от знаков cv и sv. В Lua есть соответствующая функция atan2(svcv).
В ЭМ, в том числе и ТГ, анализ фазовых соотношений проводится для первых гармоник. Для первой косинусной гармоники МПС в (7) получены при ХХ амплитуда и начальная фаза ¥m,/=51,97 Вб и уу =0, при НН - ¥mj =53,89 Вб и уі = -35.75° (в обозначениях номер гармоники не показываем). Остальные гармоники, в принципе, оказались пренебрежимо малы.
Иллюстрация углов уу и у1 есть на рис. 5, где видим, что амплитуда, а значит и вектор МПС в режиме НН сдвинуты в сторону отрицательного значения угла a по отношению к МПС в режиме ХХ.
Угол сдвига © = уу - у1 является углом нагрузки ТГ (в нашем случае он имеет положительное значение). Аналогичный сдвиг имеют и векторы F1 и Fy на рис. J.
С учетом того, что в ТГ угол поворота связан с временем, т.е. a=Q-t, где Q - угловая скорость, можно перейти к общепринятым временным векторным диаграммам (ВД) для иллюстрации и других фазовых соотношений, а также определения прочих параметров. Такая ВД представлена на рис. 6, где за начало отсчета углов принята вертикально расположенная продольная ось ротора d, как и углов на рис. J.
28
ISSN 2074-272X. Електротехніка і Електромеханіка. 2016. №1
Рис. 6. Векторная диаграмма, соответствующая режиму нагрузки ТГ
Вектор фазного тока h проведен относительно оси d под обусловленным ранее углом р (в данном случае он равен -160,43°). Так же направлен вектор МПС Tа фазной обмотки, вызванного МП реакции якоря, соответствующий вектору Fs. Векторы МПС Tf от МП обмотки ротора и результирующего МПС T проведены под уже определенными углами Jf и Jh
Вообще то, выделение Ta и Tf из МПС T является в РН условным, причем здесь значение Tf совершенно другое по сравнению с режимом ХХ. Доли МПС определим благодаря принятой здесь методике.
А именно, принимаем за базу уже определенное МПС T и из конца его вектора проводим линию параллельно вектору Ta до пересечения с линией, ориентированной, как и Tf под углом jf. Это пересечение в точке а и выявляет конкретные длины векторов Ta и Tf и их значения в масштабе вектора Т/.
На основе известной теории [2] по амплитуде МПС при частоте fs имеем действующее значение первой гармоники фазной ЭДС, в том числе и в РН:
EI = 42 л fs Tmj. (10)
Вектор этой ЭДС, как известно, отстает от вектора своего МПС Т на -90°. Векторы ЭДС Еа и Ef от МПС Ta и Tf определяются и строятся аналогично.
Получается, что треугольник ЭДС Е;, Еа и Ef подобен треугольнику МПС T Ta и Tf. Поэтому между векторами Ei и Ef будет тот же угол нагрузки ©.
Для того, чтобы определить фазное напряжение Us, необходимо учесть падение напряжения UR=Rsls на активном сопротивлении Rs и ЭДС Ev=-jXv I, от потока лобового рассеяния фазной обмотки, где Xv -соответствующее индуктивное сопротивление.
Очевидно, что вектор UR параллелен вектору I,., а вектор Ev - перпендикулярен (в сторону отставания). Так они и пристроены к концу вектора Et и дают в итоге вектор фазного напряжения Us = Et+Ev - UR.
Известно, что в ТГ величина UR пренебрежимо мала по сравнению с Et и Us, но здесь фигурирует для демонстрации полнофакторного подхода.
Из геометрических соотношений на ВД, фрагмент которой в повернутом виде показан подробнее внизу справа (без соблюдения пропорций), активная и реактивная составляющие напряжения Us, его действующее значение и фазовый сдвиг от вектора тока 1s:
Us,а = Et COS ф/ - UR; Us,r = E/ sinФ/ - Ev; (11)
Us =yj Ula + Ulr ; Ф, = arctg(Us,r / U s,а ), (12)
где по ВД угол между I, и Ei составляет ф; = -р -90 +jt.
В итоге, по полевому расчету определяется целый ряд уже рассмотренных электромагнитных параметров, а также активная электрическая мощность ТГ:
Ра = m,UsIsCOSфs. (13)
Кроме этого, можно определить электромагнитный момент ТГ, что делается через тензор магнитного натяжения [2] в соответствии с формулой
M
em
/
а
Н 0(rs - rr )
2л rs
j j r ■ Br ■ Ba
0 Гг
г ■ dr ■ da.
(14)
где rr и rs - радиусы окружностей, ограничивающие зазор со стороны ротора и статора; Br и Ba - радиальная и угловая составляющие магнитной индукции.
Магнитный поток через некоторую поверхность, ограниченную контуром /, имеет выражение [2]:
Ф = 4 -42)• /а , (15)
где Az1, 4z2 - значения ВМП в точках 1 и 2 в плоскости расчета, через которые проходят в аксиальном направлении боковые стороны упомянутого контура /.
Все обусловленные расчеты и призван автоматизировать представляемый далее скрипт Lua с соответствии с целью данной работы.
Файл исходных данных ТГ для скрипта Lua.
Требованиями к созданной программе, как и в [3, 4], являлись ее универсальность с точки зрения геометрии и размеров ТГ, а также установившихся режимов их работы - при минимуме вводимых данных.
Значения различных параметров ТГ и обслуживающих локальных констант и переменных можно задавать непосредственно в скрипте Lua, либо вводить из заранее подготовленных текстовых файлов. В этой работе все параметры, характеризующие данную задачу и возможно меняющиеся для разных вариантов ТГ, вводятся из отдельного файла.
Текстовые файлы Lua-скриптов и файлов данных пишутся в редакторе Блокнот.
Приведем файл данных, который может иметь произвольное имя с расширением txt. Данные имеют форму чисел или текстовых строк в кавычках " " и сопровождаются комментариями. Для обозначения начала комментариев выбрали двоеточие.
Текст файла исходных данных.
"FemC_TG340" : name_fem - имя модели в FEMM "RezC_TG340" : name_rez - имя файла результатов 3 : ms - число фаз обмотки статора 1 : p - число пар полюсов 50 : fs - частота, Гц
ISSN 2074-272X. Електротехніка і Електромеханіка. 2016. №1
29
-160.47 : beta - начальная фаза токов статора, град 3151.4 : Ir - ток возбуждения ротора, А 11547 : Is - фазный ток статора, А
0.8 : bs - относительное укорочение обмотки статора
0.00266 : Rs - активное сопрот. обмотки статора, Ом 0.063 : Xv - индукт. сопрот. лоб. части обм. стат., Ом 0.0 : gf - начальная фаза МПС в режиме ХХ, градус 30 : Qs - число пазов статора
10 : Ns - число последовательных витков в фазе стат.
1 : as - число параллельных ветвей обмотки статора 5.308 : la - активная длина турбогенератора, м 637.5 : rsi - радиус расточки сердечника статора, мм 699 : rv - средний радиус верхн. слоя обм. статора, мм 778 : rn -средний радиус нижнего слоя обм.статора,мм Общая характеристика языка Lua. Язык Lua, как и другие языки программирования, позволяет создавать циклы, условные операторы, процедуры, выполнять вычисления по математическим формулам, по стандартным и созданным самостоятельно функциям. Есть функции: sqrt - корень квадратный; sin, cos, atan2 - синус, косинус, арктангенс и другие. Есть встроенное число л - Pi. Тригонометрические функции оперируют с радианами, а при работе с геометрическими объектами Lua использует градусы, прямоугольные координаты измеряются в миллиметрах, что предопределяется при формулировке задачи FEMM в меню Problem, что сделано в [3].
В скрипте Lua команды, начинающиеся с mo_ или mi_, являются его процедурами. Их описание можно найти через кнопку Help в рабочем окне FEMM, войдя затем в раздел Lua Scripting [1].
Строки программы или их правая часть, начинающиеся с двойного дефиса --, являются комментариями и не влияют на ее работу.
Внешнее имя файла данных скрипт Lua запрашивает после его запуска оператором с «подсказкой»: File dan=prompt("Введите имя файла")
Это имя набирается с клавиатуры и вводится нажатием Enter. Скрипт открывает файл для чтения ("r") и присваивает ему внутреннее имя f_d, т.е. f_d = openfile(File_dan .. ".txt","r")
Оператор чтения read(f_d, "*n", "*l")
из очередной строки считает значения стольких числовых или строковых данных, сколько раз в нем фигурирует опция "*n", а опцией "*l" будет сделан перевод на новую строку, т.е. оставшийся в строке любой текст игнорируется и используется расчетчиком как комментарий (у нас отделяется символом :).
В операторах write для записи в файл результатов опция "\r\n" дает переход на новую строку, а использование символов типа %6.3f является форматом (format) печати, т.е. в данном случае под число выделяется 6 позиций, в том числе 3 - под дробную часть.
Далее приводится программа, которая после ее вызова в среде FEMM организует ввод исходной информации и расчет МП, формирует угловую функцию МПС фазной обмотки, определяет ее амплитуду и начальную фазу, а также ряд электромагнитных параметров и фазовых соотношений разных величин в ТГ. Скрипт Lua.
--Задание файла исходных данных File dan=prompt("Введите имя файла")
--Открытие файла данных
f_d = openfile(File_dan .. ".txt","r")
--ИСХОДНЫЕ ДАННЫЕ из текстового файла name fem=read(f d,"*n","*l") --имяFEMM
--Имя файла для записи результатов name_rez=read(f_d,"*n","*l") ms=read(f d,"*n","*l") --число фаз
p=read(f d,"*n","*l") --пары полюсов
fs=read(f d, "*n", "*l") --частота
Таким же способом считываются все данные из приведенного выше файла данных вплоть до rn, т.е: rn=read(f d,"*n","*l") -- средний
--радиус нижнего слоя обмотки статора closefile(f d) --закрытие файла данных --Организация файла результатов writeto(name rez) --name rez имя файла --Имя геометрической модели ТГ-тип fem name fem=name fem .. ".fem" mi saveas(name fem) --запомнить модель “-ДОПОЛНИТЕЛЬНЫЕ РАСЧЕТНЫЕ ДАННЫЕ: gr=Pi/180 --множитель градус=>радиан beta=beta*gr --значение угла beta, рад Ur=Rs*Is --падение напряжения на Rs Ev=Xv*Is --падение напряжения на Xv ts=360/Qs --угол зубц. шага статора tp=Qs/(2*p) --полюсный шаг в зуб. дел.
qsp=tp/ms --число пазов на полюс и фазу qsn=0.5*(qsp-tp*(1-bs)) --число пазов --нижнего слоя обм. ст-ра над осью q qsv=qsp-qsn --число пазов верхнего -- слоя обмотки статора над осью q --Угол первого нижнего стержня фазы А an1=180-(qsn-0.5)*ts
--Угол первого верхнего стержня фазы А av1=180-(qsv-0.5)*ts Im=Is*2A0.5/as --амплитуда фазного -- тока в параллельной ветви --Запись некоторых параметров ТГ write("Данные возбуждения ТГ","\г\п") write(format(" Is= %5.0f",Is), format(" Ir= %5.0f",Ir), format(" beta=%7.2f",beta/gr),"\r\n")
F = {} --создание массива для хранения -- значений МПС фазной обмотки -- Фазные токи обмотки статора с=2*Рі/3 Ia=Im*cos(beta) Ib=Im*cos(beta-^ Ic^m^cos^eta+с)
-- Задание токов фаз и ротора в FEMM mi modifycircprop("iA+",1,Ia) mi modifycircprop("iA-",1,-Ia) mi modifycircprop("iB+",1,Ib) mi modifycircprop("iB-",1,-Ib) mi modifycircprop("iC+",1,Ic) mi modifycircprop("iC-",1,-Ic) mi modifycircprop("Ir+",1, Ir) mi modifycircprop("Ir-",1,- Ir)
—В FEMM-модели есть имена меток токов -- Расчет МП в FEMM и переход к показу --и извлечению результатов расчета mi analyze(l) mi loadsolution() --Цикл условного перемещения фазной -- зоны A с пазовым шагом ts for k=1,Qs,1 do ak=(k-1)*ts --угол --Цикл расч. углов проводников фазы А for i=1,qsp,1 do ai=ak+(i-1)*ts --Расчет координат проводников --нижнего слоя и их выделение
30
ISSN 2074-272X. Електротехніка і Електромеханіка. 2016. №1
ani=(an1+ai)*gr x=rn*cos(ani) y=rn*sin(ani) mo selectblock(x,y) --Расчет координат проводников верхнего слоя и их выделение avi=(av1+ai)*gr x=rv*cos(avi) y=rv*sin(avi) mo selectblock(x,y) end --завершение цикла по i
--Считывание площади выделенных блоков S = mo blockintegral(5)
--Значение интеграла ВМП по площади S A = mo blockintegral(1) mo clearblock() --очистка выделения --МПС фазной зоны обмотки А F[k] = Ns*A/S
end --завершение пазовых шагов по k --Получение МПС фазной обмотки по МПС --фазной зоны и продление на период q2=Qs/2 --половина позиций фазной зоны for i=1,q2,1 do F[i]=F[i]-F[i+q2] F[i+q2]=-F[i] end
--Расчет амплитуды и фазового сдвига --первой гармоники МПС по ряду Фурье Fsin=0 Fcos=0 --обнуление накопителей for k=1,Qs,1 do —цикл по массиву МПС ak=ts*(k-1)*gr --угловая позиция --Амплитуды sin и cos составляющих Fsin=Fsin+2*F[k]*sin(ak)/Qs Fcos=Fcos+2*F[k]*cos(ak)/Qs end --перебора k элементов массива Fm1=sqrt(FsinA2+FcosA2) --амплитуда gl=atan2(Fsin,Fcos) --начальная фаза --Электромагнитные параметры ТГ --Действующее значение фазной ЭДС El=Pi*2A0.5*fs*Fm1
--Фазовый сдвиг ЭДС El от тока Is fil=-(beta+Pi/2-gl)
--Активная и реактивная составляющие --фазного напряжения
Usa=El*cos(fil)-Ur Usr=El*sin(fil)-Ev --Действующее знач. фазного напряжения Us=sqrt(UsaA2+UsrA2)
--Фазовый сдвиг между Us и Is --и коэффициент мощности ТГ fis=atan(Usr/Usa) cosfi=cos(fis) teta=gf*gr-gl --угол нагрузки ТГ Pa=ms*Us*Is*cosfi --активная мощность --Запись результатов в файл write(" МПС, фазовые сдвиги, ЭДС, напряжение и акт. мощность","\r\n")
-- Пересчет углов ( радианы в градусы) gl=gl/gr fil=fil/gr fis=fis/gr teta=teta/gr beta=beta/gr write(format(" Fm1=%5.2f",Fm1), format(" gl=%7.2f",gl),"\r\n") write(format(" El=%5.0f",El), format(" fil=%7.2f",fil),"\r\n") write(format(" gf=%4.1f",gf), format(" teta=%6.2f",teta), format(" Us=%5.0f",Us),"\r\n") write(format(" fis=%7.2f",fis), format(" cosfi=%4.2f",cosfi), format(" Pa=%6.1f",Pa*1e-6),"\r\n") write("Электромагные момент и мощность","\r\n")
--Выделение группы блоков ротора mo groupselectblock(9)
--Считывание электромагнитного момента Mem=mo blockintegral(22) Mem=Mem*1e-3 mo clearblock() --очистка выделения --Электромагнитная мощность Pem=Mem*2*Pi*fs*1e-3
write(format("Mem=%5.0f",Mem)," кН*м", format("Pem=%6.1f",Pem)," МВт","\г\п") --Магнитный поток на полюс --в зазоре на расточке статора --Цикл угловых шагов в da градусов for i=1,180,1 do da=1 --Координаты двух точек рядом ai=da*(i-0.5)*gr ai0=da*(i+0.5)*gr
x=-rsi*sin(ai) y=rsi*cos(ai) x0=-rs^sin(ai0) y0=rsi*cos(ai0)
--Считывание ВМП в двух точках A=mo getpointvalues(x,y)
A0=mo getpointvalues(x0,y0)
--Поиск максимума ВМП на окружности if A0<A then Pot=2*A*la i=181 end end --перебора i точек по окружности write("Магнитный поток на полюс, Вб", format(" Pot=%5.3f",Pot),"\r\n") writeto() --Закрытие файла результатов mi close() --Закрытие документа --препроцессора магнетизма
Указания по использованию программы.
Для использования скрипта полевого расчета параметров и фазовых соотношений ТГ необходимо:
1. Файл скрипта с расширением *.lua и файл с исходными данными с расширением *.txt должны находиться в одной папке на диске компьютера.
2. Запускаем программу FEMM - модель ТГ. Редактор Lua запускается из рабочего окна в FEMM «кнопкой» Open Lua Script в меню File, а затем в открывшемся окне запускается необходимый файл lua.
3. По запросу программы в открывающееся окно вводим имя файла данных и нажимаем клавишу Enter. Дальнейшие действия программа выполняет сама.
4. Результаты расчета находятся в текстовом файле с заданным в файле данных именем RezC_TG340 и читаются редактором Блокнот.
Содержимое файла результатов.
Приводимые данные соответствуют режиму НН и уже иллюстрировались на рис. 5 и рис. 6. Для удобства представления в формате данной статьи текст файла результатов несколько переформатирован. Данные возбуждения ТГ Is=11547 beta=-160.43 Ir=3153
МПС, фазовые сдвиги, ЭДС, напряжение и активная мощность Fm1= 53.89 gl= -35.75 El=11971 fil= 34.68 gf= 0.0 teta= 35.75 Us=11547
fis= 31.79 cosfi=0.85 Pa=340.0 Электромагнитные момент и мощность Mems=-1086 кН*м; Pems=-341.2 МВт
Магнитный поток на полюс статора, Вб Pot=6.608
В полной версии скрипта Lua вычисляются и выводятся в текстовый файл и другие результаты.
Вывод. Представленные теоретические основы и скрипт Lua дают широкие возможности пользователям программы FEMM для автоматизированного получе-
ISSN 2074-272X. Електротехніка і Електромеханіка. 2016. №1
31
ния электромагнитных параметров и фазовых соотношений в ТГ на основе расчета магнитного поля. Этот скрипт является универсальным с точки зрения структуры электромагнитной системы ТГ в рамках их распространенной конструкции, а также установившихся режимов их работы. Разработанный скрипт Lua может явиться прототипом для аналогичного программного обеспечения расчетов ЭМ других типов.
СПИСОК ЛИТЕРАТУРЫ
1. Meeker D. Finite Element Method Magnetics. FEMM 4.2 11 Oct 2010 Self-Installing Executable. Режим доступа: www. femm. info/wiki/OldVersions.
2. Милых В.И., Полякова Н.В. Численно-полевые расчеты электромагнитных параметров турбогенераторов // Вісник НТУ «ХПІ». - 2014. - №38(1081). - С. 3-18.
3. Милых В.И., Полякова Н.В. Автоматизированное формирование расчетных моделей турбогенераторов для программной среды FEMM // Електротехніка і електромеханіка. - 2015. - №4. - С. 7-14.
4. Милых В.И., Полякова Н.В. Автоматизированные расчеты в программной среде FEMM динамики электромагнитных процессов турбогенераторов // Електротехніка і електромеханіка. - 2015.- №6. - С. 24-30.
5. Титов В.В., Хуторецкий Г.М., Загородная Г.А. и др. Турбогенераторы. - Л.: Энергия, 1967. - 895 с.
6. Милых В.И., Полякова Н.В. Анализ фазовых соотношений электромагнитных величин в турбогенераторе на основе численных расчетов магнитных полей // Електротехніка і електромеханіка. - 2003. - №4. - С. 59-64.
7. Милых В.И., Полякова Н.В. Система направлений и фазовых соотношений электромагнитных величин при численных расчетах магнитных полей в турбогенераторе // Електротехніка і електромеханіка. - 2011. - №5. - С. 33-38.
8. Милых В.И., Полякова Н.В. Основы численного анализа фазовых соотношений электромагнитных величин в турбогенераторе // Электрика. - 2012. - №3. - С. 31-33.
9. Милых В.И., Полякова Н.В. Организация численного расчета магнитного поля турбогенератора в режиме нагрузки с обеспечением заданных его выходных параметров // Електротехніка і електромеханіка. - 2012. - №1. - С. 36-41.
REFERENCES
1. Meeker D. Finite Element Method Magnetics. FEMM 4.2 32 bit 11 Oct 2010 Self-Installing Executable. Available at: www.femm.info/wiki/OldVersions (accessed 10 March 2014).
2. Milykh V.I., Polyakova N.V. Numerically-field calculations of the electromagnetic parameters of turbogenerators. Visnyk NTU «KhPI» - Bulletin of NTU «KhPI», 2014, no.38(1081), pp. 3-18. (Rus).
3. Milykh V.I., Polyakova N.V. Automated formation of calculation models of turbogenerators for software environment FEMM. Elektrotekhnika і elektromekhanika - Electrical engineering & electromechanics, 2015, no.4, pp. 7-14. (Rus).
4. Milykh V.I., Polyakova N.V. Automated calculations of the dynamics of turbogenerator electromagnetic processes in software environment FEMM. Elektrotekhnika і elektromekhanika - Electrical engineering & electromechanics, 2015, no.6, pp. 24-30. (Rus).
5. Titov V.V., Hutoreckij G.M., Zagorodnaja G.A., Vartan'jan G.P., Zaslavsky D.I., Smotrov I.A. Turbogeneratory [Turbogenerators]. Leningrad, Energiia Publ., 1967. 895 p. (Rus).
6. Milykh V.I., Polyakova N.V. Analysis of phase relationships of electromagnetic parameters in a turbogenerator on the basis of numerical calculation of magnetic fields. Elektro-tekhnika і elektromekhanika - Electrical engineering & electromechanics, 2003, no.4, pp. 59-64. (Rus).
7. Milykh V.I., Polyakova N.V. A system of directions and phase relationships for electromagnetic parameters at numerical calculations of magnetic fields in a turbogenerator. Elektrotekhnika i elektromekhanika - Electrical engineering & electromechanics, 2011, no.5, pp. 33-38. (Rus).
8. Milykh V.I., Polyakova N.V. Bases of numerical analysis of phase correlations of electromagnetic sizes are in a turbogenerator. Elektrika - Electrician, 2012, no.3, pp. 31-33. (Rus).
9. Milykh V.I., Polyakova N.V. Organization of numerical calculation of turbogenerator magnetic field under load with specified output parameters control. Elektrotekhnika i elek-tromekhanika - Electrical engineering & electromechanics, 2012, no.1, pp. 36-41. (Rus).
Поступила (received) 22.07.2015
Милых Владимир Иванович1, д.т.н., проф.,
Полякова Наталия Владимировна1, инженер,
1 Национальный технический университет «Харьковский политехнический институт»,
61002, Харьков, ул. Фрунзе, 21,
тел/phone +38 057 7076514, e-mail: [email protected] V.I. Milykh1, N.V. Polyakova1
1 National Technical University «Kharkiv Polytechnic Institute», 21, Frunze Str., Kharkiv, 61002, Ukraine.
Determination of electromagnetic parameters and phase relations in turbo-generators by the automated calculation of the magnetic field in the software environment FEMM.
The theoretical bases of calculation of electromagnetic quantities and time-phase relationship are presented for the turbogenerators. This is done by numerical calculations of the magnetic field in the software environment package FEMM (Finite Element Method Magnetics). A program which controls calculations and organizes the issuance of the results to a text file is created on the algorithmic language Lua. The program is universal in terms of a turbo-generator models, as well as steady-state modes of their work with a minimum of input data. The exciting current of the rotor and the phase currents of three-phase stator winding in accordance with their initial phase are given for the calculation of the magnetic field. The key function for the analysis of electromagnetic parameters is the calculated angular function of the magnetic flux phase stator winding. The expansion in the harmonic series is carried out and amplitude and initial phase are received for this function. Next, the phase EMF and voltage, phase shifts between all values, active power, electromagnetic torque, the magnetic flux in the gap and other parameters are determined. The presented Lua script is a prototype for a similar calculation software of electric machines of other types. References 9, figures 6.
Key words: Finite Element Method Magnetics, program FEMM, turbo-generator, electromagnetic parameters, phase relationships, automated calculations, Lua script.
32
ISSN 2074-272X. Електротехніка і Електромеханіка. 2016. №1
