Автоматизированное формирование расчетных моделей турбогенераторов для программной среды FEMM
Представлены принципы автоматизированного формирования графической и физической моделей турбогенераторов
 на алгоритмическом языке Lua для расчета их магнитных полей и параметров в программной среде пакета FEMM.
 Возможности составленного скрипта Lua демонстрируются на примере реальн...
Saved in:
| Published in: | Електротехніка і електромеханіка |
|---|---|
| Date: | 2015 |
| Main Authors: | , |
| Format: | Article |
| Language: | Russian |
| Published: |
Інститут технічних проблем магнетизму НАН України
2015
|
| Subjects: | |
| Online Access: | https://nasplib.isofts.kiev.ua/handle/123456789/149284 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Journal Title: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Cite this: | Автоматизированное формирование расчетных моделей турбогенераторов для программной среды FEMM / В.И. Милых, Н.В. Полякова // Електротехніка і електромеханіка. — 2015. — № 4. — С. 7–14. — Бібліогр.: 4 назв. — рос. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| _version_ | 1860255907379675136 |
|---|---|
| author | Милых, В.И. Полякова, Н.В. |
| author_facet | Милых, В.И. Полякова, Н.В. |
| citation_txt | Автоматизированное формирование расчетных моделей турбогенераторов для программной среды FEMM / В.И. Милых, Н.В. Полякова // Електротехніка і електромеханіка. — 2015. — № 4. — С. 7–14. — Бібліогр.: 4 назв. — рос. |
| collection | DSpace DC |
| container_title | Електротехніка і електромеханіка |
| description | Представлены принципы автоматизированного формирования графической и физической моделей турбогенераторов
на алгоритмическом языке Lua для расчета их магнитных полей и параметров в программной среде пакета FEMM.
Возможности составленного скрипта Lua демонстрируются на примере реального турбогенератора.
Представлені принципи автоматизованого формування графічної і фізичної моделей турбогенераторів на
алгоритмичній мові Lua для розрахунку їх магнітних полів і параметрів у програмному середовищі пакету FEMM.
Можливості складеного скрипту Lua демонструються на прикладі реального турбогенератора.
Attention is paid to the popular FEMM (Finite Element Method
Magnetics) program which is effective in the numerical calculations of the magnetic fields of electrical machines. The main problem of its using - high costs in time on the formation of a graphical
model representing the design and on the formation of the physical
model representing the materials properties and the winding currents of machines – is solved. For this purpose, principles of the
automated formation of such models are developed and presented
on the turbogenerator example. The task is performed by a program written in an algorithmic language Lua integrated into the
package FEMM. The program is universal in terms of varying the
geometry and dimensions of the designed turbogenerators. It uses a
minimum of input information in a digital form representing the
design of the whole turbogenerator and its fragments. A general
structure of the Lua script is provided, significant parts of its text,
the graphic results of work's phases, as well as explanations of the
program and instructions for its use are given. Performance capabilities of the compiled Lua script are shown on the example of the
real 340 MW turbogenerator.
|
| first_indexed | 2025-12-07T18:49:12Z |
| format | Article |
| fulltext |
Електричні машини та апарати
ISSN 2074-272X. Електротехніка і Електромеханіка. 2015. №4 7
© В.И. Милых, Н.В. Полякова
УДК 621.313:681.3
В.И. Милых, Н.В. Полякова
АВТОМАТИЗИРОВАННОЕ ФОРМИРОВАНИЕ РАСЧЕТНЫХ МОДЕЛЕЙ
ТУРБОГЕНЕРАТОРОВ ДЛЯ ПРОГРАММНОЙ СРЕДЫ FEMM
Представлені принципи автоматизованого формування графічної і фізичної моделей турбогенераторів на
алгоритмичній мові Lua для розрахунку їх магнітних полів і параметрів у програмному середовищі пакету FEMM.
Можливості складеного скрипту Lua демонструються на прикладі реального турбогенератора. Бібл. 4, рис. 15.
Ключові слова: програма FEMM, метод скінченних елементів, турбогенератор, графічні та фізичні моделі, автомати-
зоване формування, Lua скрипт.
Представлены принципы автоматизированного формирования графической и физической моделей турбогенераторов
на алгоритмическом языке Lua для расчета их магнитных полей и параметров в программной среде пакета FEMM.
Возможности составленного скрипта Lua демонстрируются на примере реального турбогенератора. Библ. 4, рис. 15.
Ключевые слова: программа FEMM, метод конечных элементов, турбогенератор, графические и физические модели,
автоматизированное формирование, Lua скрипт.
Введение. Существуют программные продукты
для численных расчетов магнитных полей (МП) раз-
личных электротехнических устройств. Программы
позволяют учитывать достаточно точно реальную
геометрию устройств, а также нелинейные магнитные
свойства ферромагнитных сердечников. Среди таких
продуктов широкое распространение получила про-
грамма FEMM [1], основанная на методе конечных
элементов и ориентированная на расчет двухмерных
МП. Этого достаточно для электрических машин
(ЭМ) классической конструкции, имеющих плоскопа-
раллельное МП в пределах активной части.
Имеющиеся программные продукты уже не тре-
буют от пользователя разработки каких-либо матема-
тических моделей. Достаточно следовать инструкци-
ям при построении графических моделей, отобра-
жающих конструкцию ЭМ, а также физических моде-
лей, отображающих свойства материалов, токи обмо-
ток и параметры постоянных магнитов. Аналогично,
надо следовать инструкциям при извлечении резуль-
татов расчета МП и их последующем использовании.
При эксплуатации программы FEMM обуслов-
ленные действия выполнялись и еще выполняются
«вручную». Наибольшие затраты времени требуют
формирование графических и физических моделей и
преобразование результатов расчетов в форму, удоб-
ную для использования в процессе проектирования и
исследования ЭМ. Это обуславливает низкую произ-
водительность расчетчиков и часто становится пре-
пятствием ввиду необходимости освоения большого
объема инструкций и ошибок при следовании им.
Эти проблемы становятся еще серьезней при
многовариантных расчетах МП и электромагнитных
параметров ЭМ, что необходимо для варьирования
токами их обмоток, геометрией и размерами, взаим-
ным положением неподвижных и подвижных частей.
Проблемы можно преодолеть посредством авто-
матизации построения графических и физических
моделей ЭМ, управления процессом расчета и пере-
бора параметров, а также извлечения, обработки и
приведения в необходимую для дальнейшего исполь-
зования форму результатов. Такие возможности поль-
зователям программы FEMM предоставляет интегри-
рованный в нее алгоритмический язык Lua. И задача
теперь заключается в написании соответствующего
комплекса программ (скриптов). Эта задача весьма
сложна, так как требует комплексного владения зна-
ниями в областях ЭМ, программирования и теории
МП. Но эта задача вполне окупаема, так как програм-
ма пишется один раз, а используется бесчисленное
число раз и любым количеством пользователей.
Для обеспечения высокой эффективности расчетов
МП и электромагнитных параметров ЭМ в среде
FEMM, нами разработан комплекс скриптов Lua. Эти
программы позволяют многократно повысить произво-
дительность расчетчиков, а также обеспечивают воз-
можность получения таких параметров и характеристик
ЭМ, которые практически невозможно получить в
«ручном» режиме работы с программой FEMM.
Цель работы. Все программы показать в одной
статье не представляется возможным. Поэтому целью
данной работы является представление принципов
составления скрипта Lua для автоматизированного
формирования графической и физической моделей
турбогенератора (ТГ) – одного из наиболее крупных и
ответственных представителей семейства ЭМ [2].
Требованием к такой программе является ее уни-
версальность, с точки зрения варьирования геометрией
и размерами проектируемых ТГ, при минимуме вводи-
мой информации – и только в числовой форме.
Объект исследования. Для иллюстраций взяли
ТГ мощностью 340 МВт. Его электромагнитная сис-
тема представлена на рис. 1 поперечным сечением.
Здесь разной штриховкой стержней в пазах ста-
тора показаны фазные зоны его обмотки: А–А', В–В' и
С–С'. При построениях моделей ТГ используется
прямоугольная система координат x, y и принято:
МДС фазной обмотки А–А' и продольная ось ротора d
ориентированы по оси у (при дальнейших расчетах
ротор может вращаться); с осью х совпадают оси
симметрии двух пазов статора. Но, в зависимости от
относительного укорочения обмотки статора s, воз-
можно расположение по оси х и зубцов статора [3].
Это регулируется в программе посредством ширины
фазной зоны обмотки статора (в зубцовых делениях)
8 ISSN 2074-272X. Електротехніка і Електромеханіка. 2015. №4
)1( spss qz . (1)
Здесь используются число пазов на полюс и фазу
и полюсное деление статора (в зубцовых делениях)
s
s
s mp
Q
q
2
;
p
Qs
p
2
, (2)
где ms, p, Qs – числа фаз, пар полюсов и пазов статора.
Если zs – нечетное число, то по оси x на статоре
будет располагаться пазы, если четное – то зубцы.
rse
rre
de
lta
rs
irsn
rr
n
Рис. 1. Модель электромагнитной системы ТГ
Общая характеристика скрипта Lua. Тексты
файлов скриптов и файлов данных можно писать в
редакторе Блокнот. В строках этих файлов все после
двойного дефиса (--) программа воспринимает как
комментарии и для работы не использует.
В скриптах базовой является команда задания
параметров новой Задачи: ее тип (МП или другое),
единицы измерения размеров, глубина (активная дли-
на) по аксиальной оси, тип координат, точность и т.д.
Построение графических моделей посредством
Lua основано на командах расстановки узлов – коор-
динатных точек (x,y). На их основе своими командами
строятся прямые отрезки и дуги. В образующихся
замкнутых подобластях – блоках ставятся их метки,
через которые им приписываются физические свойст-
ва. Они задаются отдельными командами, в которых
указываются все необходимые параметры материа-
лов, токов, треугольников конечно-элементной струк-
туры и др. Есть также свои команды для выбора и
расстановки граничных условий. Перечисленные выше
объекты можно объединять в нумерованные группы и
оперировать со всеми объектами группы сразу.
Lua, как и другие языки программирования, позво-
ляет создавать функции и процедуры, оперировать с
циклами, условными операторами, делать вычисления
по формулам, стандартным функциям: ceil – целая часть
результата; sqrt – корень квадратный; sin, cos, tan, atan –
синус, косинус, тангенс, арктангенс и др.
Для тригонометрических функций углы надо за-
давать в радианах, а при построении геометрических
объектов Lua оперирует углами в градусах.
Скрипты вызываются (запускаются) из рабочего
окна в FEMM «кнопкой» Open Lua Script в меню File.
Конструкторская модель турбогенератора.
Представим исходные параметры ТГ, которые зада-
ются в файле данных для полного формирования гра-
фической и физической моделей ТГ. Далее исполь-
зуемые величины обозначены так, чтобы это соответ-
ствовало правилам идентификации в скрипте Lua.
К исходными схемно-конструктивным парамет-
рам ТГ (рис. 1) относятся: rre, rsi – радиусы окружно-
стей по поверхностям сердечников: внешней ротора и
расточки статора; has – высота спинки сердечника
статора; p, ms – числа пар полюсов и фаз; Qs, Ncs, as,
bss – числа пазов, эффективных проводников в
стержнях и параллельных ветвей, относительное уко-
рочение шага обмотки статора (s); Qr, Qru, Ncr –
числа пазов ротора обмотанных и условно-полное, а
также эффективных проводников в пазу; la – активная
длина статора; dsl – толщина листов его сердечника.
Кроме названного, детализации требуют зубцово-
пазовые структуры статора и ротора, общепринятые
формы пазов которых с обозначением необходимых
размеров и опорных точек даны на рис. 2 и рис. 3.
1as 2as
rsi
xg
xe
xd
xbxc
xf
xa
hks
ybya
hss
11
129
8
6
13 14
5
4
2
31
bk
s bs
rsn
hs
7
1b
1a
2a
x
y
10
Рис. 2. Конструкция паза статора
2
10
bk
r
br
3
2ar
57
9
11
1ar
hsr
812
hkr
hr
rre
rrn
xb
ybya
1b
2a
y
1
xa
4
14
xcxdxexg
1a
x
Рис. 3. Конструкция паза ротора
Паз статора имеет обозначения (выделены на
рис. 2): hs, bs– высоты и ширина; hss – высота шлица;
hks, bks – высота и ширина выемок под клин; as1, as2
– 1-ый и 2-ой углы их скоса. Для паза ротора приняты
обозначения (рис. 3): hr, br – высота и ширина; hsr –
высота шлица; hkr, bkr – высота и ширина выемок под
клин; ar1, ar2 – 1-ый и 2-ой углы их скоса.
На основе перечисленных исходных данных вы-
числяются еще параметры графической модели:
rse=rsi+hs+has – внешний радиус сердечниеа статора;
delta= rsi-rre – воздушный зазор;
ts =2*Pi*rsi/Qs – пазовый шаг по расточке статора;
tr=2*Pi*rre/Qru – пазовый шаг по поверхности ротора;
tsa =360/Qs – пазовый шаг статора (в градусах);
tra=360/Qru – пазовый шаг ротора (в градусах);
bts = ts-bs – ширина зубца статора по его поверхности;
btr = tr-br – ширина зубца ротора по его поверхности;
ISSN 2074-272X. Електротехніка і Електромеханіка. 2015. №4 9
rsn=sqrt(x*x+y*y) – радиус статора по углам дна паза,
где x=bs/2; y=hs+sqrt(rsi*rsi-x*x);
rrn=sqrt(x*x+y*y) – радиус ротора по углам дна паза,
где x=br/2; y=sqrt(rre*rre-x*x)-hr;
bzs=2*Pi*rsn/Qs-bs – ширина основания зубца статора;
bzr=2*Pi*rrn/Qru-br – ширина основания зубца ротора;
qsp=Qs/(2*ms*p) – число пазов ст-ра на полюс и фазу;
tp=Qs/(2*p) – полюсный шаг в зубцовых делениях;
zs=ceil(qsp+tp*(1-bss)) – ширина фазной зоны статора
по его крайним стержням в зубцовых делениях.
Координатные размеры элементов паза статора:
ya=bs/2; yb=bks/2; xa=sqrt(rsi*rsi-ya*ya); xb=xa+hss;
xe=xa+hks; xf=xa+(hks+hs)/2; xg=xa+hs; b1=yb-ya;
a1=b1*tan(as1); xc=xb+a1; a2=b1*tan(as2); xd=xe-a2.
Координатные размеры элементов паза ротора:
ya=br/2; yb=bkr/2; xa=sqrt(rre*rre-ya*ya); xb=xa-hsr;
xe=xa-hkr; xg=xa-hr; b1=yb-ya; a1=b1*tan(ar1);
xc=xb-a1; a2=b1*tan(ar2); xd=xe+a2.
Электромагнитные свойства и параметры
турбогенератора в его физической модели.
Магнитные свойства сердечников ТГ представ-
ляются их кривыми намагничивания – зависимостями
магнитной индукции от напряженности магнитного
поля, т.е. B(H). Эти кривые вводятся в программу
FEMM парами значений B и H (по сути – таблицами),
а затем используются программой автоматически.
В ТГ используются кривые намагничивания для
ярма сердечника статора, для его зубцов и для сер-
дечника ротора [2]. Условные имена кривых предва-
рительно вводятся из общего файла исходных дан-
ных, также как и имена соответствующих им файлов,
откуда берутся массивы кривых намагничивания.
Приняты (произвольно) такие обозначения: St_j_st,
St_z_st, St_rot – имена сталей для спинки и зубцов
сердечника статора и для сердечника ротора; bhj, bhz,
bhr– имена файлов с кривыми намагничивания для
этих сталей.
Для учета насыщения магнитопровода вводятся:
kFes – коэффициент заполнения сердечника статора с
учетом пакетирования; kFer – коэффициент приведе-
ния длины ротора к полной длине сердечника статора.
После формирования расчетной модели ТГ воз-
можен тестовый расчет его МП. Для этого вводятся: Ir
– ток обмотки ротора; Is – действующее значение фаз-
ного тока статора; beta – угол сдвига между направле-
ниями действия МДС обмоток ротора и статора в соот-
ветствии с принятой системой фазовых соотношений
ТГ [4]; fs – частота напряжения и тока статора.
Расчет МП проводится на конечно-элементной
структуре, которая генерируется программой FEMM
автоматически. Для варьирования ею вводятся мак-
симально допустимые размеры сторон треугольников:
fe1, fe2, fe3, fe4, fe5, fe6, fe7 – в зазоре; в ярме, зубцах
и пазах статора; в пазах, зубцах и бочке ротора.
Структура программы построения графиче-
ской и физической моделей турбогенератора.
Основные логические части программы таковы:
1. Задание имени файла исходных данных.
2. Открытие этого файла и чтение из него данных.
3. Задание имени расчетной модели ТГ.
4. Расчет дополнительных параметров на основе ис-
ходных данных.
5. Создание магнитной Задачи и ее параметров.
6. Создание внешней поверхности сердечника статора.
7. Создание зубцово-пазовой структуры статора.
8. Создание зубцово-пазовой структуры ротора.
9. Дорисовка поверхностей больших зубов ротора.
10. Создание свойства материалов с магнитной про-
ницаемостью воздуха.
11. Создание параметров сердечников – имен магнит-
ных свойств магнитопровода.
12. Ввод кривых намагничивания сталей.
13. Подготовка фазных токов обмотки статора.
14. Подготовка параметров геометрической структу-
ры обмотки статора.
15. Задание материала и токов в стержнях шести фаз-
ных зон обмотки статора.
16. Задание материала и токов в пазах ротора.
17. Задание магнитных свойств воздушного зазора.
18. Создание граничных условий Дирихле и назначе-
ние их внешней поверхности статора.
19. Задание материала для ярма сердечника статора.
20. Задание материала для зубцов сердечника статора.
21. В зависимости от укорочения обмотки статора и
ширины фазной зоны статора (1) поворачивается или
нет его зубцовая структура на ползубцового шага.
22. Задание материала для зубцов сердечника ротора.
23. Задание материала для бочки ротора.
24. Запоминание созданной модели турбогенератора.
25. Выполнение тестового расчета магнитного поля.
26. Переход к показу результатов расчета.
Файлы исходных данных для скрипта Lua.
Эти файлы могут начинаться со строк комментариев,
которые не участвуют в работе программы, но при
вводе информации главной программой должны быть
прочитаны – «сняты» как текстовые строки.
Программой вводятся числовые исходные дан-
ные, либо текстовые строки в двойных кавычках (" ").
Единица измерения геометрических размеров
указываются в данных для главного процессора про-
граммы FEMM – обычно принимаются миллиметры.
Для примера подачи информации в файле дан-
ных дадим только отдельные фрагменты исходных
параметров из перечня, представленного выше.
Примеры исходных данных в файле D_TG340.txt:
"TG340" -- name_TG - имя расчетной модели ТГ
"18.04.2015" -- data_dan - дата составления данных
…………………….
1 -- p - число пар полюсов
3 -- ms - число фаз обмотки статора
30 -- Qs - число пазов статора
36 -- Qr - число пазов ротора обмотанных
52 -- Qru - число пазов ротора общее (условное)
637.5 -- rsi - радиус расточки статора
…………………….
0.8 -- bss – относит. укорочение обмотки статора
2206 -- Ir - значение тока возбуждения, А
11547 -- Is - действующее значение тока статора, A
-160.47 -- beta - значение угла бета, градус
…………………….
"St3" -- St_rot - имя стали для сердечника ротора
"bhr" -- bhr -имя файла с кривой намагнич. для St_rot
…………………….
8 -- fe1 - макс. размер треугольника в зазоре
25 -- fe2 - макс. размер треугольника в ярме статора
…………………….
10 ISSN 2074-272X. Електротехніка і Електромеханіка. 2015. №4
Содержимое файла "bhr" с кривой намагничивания.
--Сталь Ст3 - для ротора: пары значений B-H (Тл,А/м)
0.1 80 0.2 160 0.3 240 0.4 320 0.5 400 0.6 488 0.7
584 0.8 682 0.9 792 1.0 924 1.1 1090 1.2 1290 1.3
1590 1.4 2090 1.5 2890 1.6 4100 1.7 5700 1.8 8000
1.9 13000 2.0 23500 2.1 42500 2.15 62000 2.2 95000
2.25 132000 2.3 170000 2.35 209000 2.4 249000 0 0
Нулевую точку кривых намагничивания про-
грамма FEMM добавляет сама, а нули в конце файлов
сигнализируют о конце вводимого массива кривой.
Скрипт Lua для формирования модели ТГ.
Далее приводятся фрагменты программы, которая
после ее вызова в среде FEMM формирует графиче-
скую и физическую модели ТГ и запускает программу
FEMM на расчет МП.
В тексте скрипта Lua команды, начинающиеся с
mi_, являются стандартными процедурами. Их описа-
ние можно найти через кнопку Help в рабочем окне
FEMM, войдя затем в раздел Lua Scripting.
Приведем часто встречающиеся команды:
1) добавление: mi_addnode(x,y) – узел;
mi_addsegment(x1,y1,x2,y2) – отрезок между узлами;
mi_addarc(x1,y1,x2,y2,angle,maxseg) – дуга между уз-
лами с опорным углом ‘angle’ из угловых сегментов
размером ‘maxseg’; mi_addblocklabel(x,y) – ярлык
блока, где (x,y), (x1,y1) и (x2,y2) – координаты узлов;
2) выбор ближайших к точке (x,y) объектов::
mi_selectnode(x,y) – узел; mi_selectsegment(x,y) – отре-
зок; mi_selectarcsegment(x,y) – дуга; mi_selectlabel(x,y)
– ярлык блока; mi_selectgroup(n) – n-ая группа;
3) копирование: mi_mirror(x1,y1,x2,y2,edit) – вы-
бранных объектов относительно линии, проходящей
через точки (x1,y1) и (x2,y2), edit – действие редакти-
рования: 0 – для узлов, 1 – отрезков, 2 – дуг, 3 – ме-
ток, 4 – групп; mi_copyrotate(bx, by, angle, copies, edit)
– copies – столько раз с поворотом на угол angle отно-
сительно базовой точки bx, by; edit – см. выше;
4) очистка выделения всех ранее выбранных
объектов: mi_clearselected();
5) установка свойства материала конкретной по-
добласти – блока: mi_setblockprop("blockname",
automesh, meshsize, "incircuit", magdir, group, turns),
где: "blockname" – имя ярлыка для магнитных свойств
блока, automesh – выбор способа треангуляции ко-
нечных элементов, meshsize – ограничивающий раз-
мер элемента, "incircuit" – символ тока, magdir – на-
правление намагничивания, group – номер группы,
turns – число эффективных проводников с током.
В программе используется нумерация точек па-
зов статора и ротора, приведенная на рис. 2 и 3.
По ходу программы вставлены рисунки, показы-
вающие результат ее выполнения на данный момент.
Текст программы.
name_prog="Lua_M_TG" --имя программы
F_dan="D_TG340.txt" --имя файла данных
f_d=openfile(F_dan,"r") ––его открытие
--считывание двух строк заголовка
for i=1,2,1 do str=read(f_d,"*l") end
--Ввод данных из File_dan по имени f_d
name_TG=read(f_d, "*n", "*l")
data_dan=read(f_d, "*n", "*l")
p= read(f_d, "*n", "*l") ms=read(f_d, "*n", "*l")
Qs=read(f_d,"*n","*l") Qr=read(f_d,"*n","*l")
Qru=read(f_d,"*n","*l") rsi=read(f_d,"*n","*l")
…………………….
bss=read(f_d, "*n", "*l") Ir=read(f_d, "*n", "*l")
Is=read(f_d, "*n", "*l") beta=read(f_d, "*n", "*l")
…………………….
-- имя стали для сердечника ротора
St_rot=read(f_d, "*n", "*l")
-- файл с кривой намагнич.для St_rot
bhr=read(f_d, "*n", "*l")
…………………….
fe1=read(f_d, "*n", "*l")
fe2=read(f_d, "*n", "*l")
…………………….
-- Досчет геометрических параметров
gr=Pi/180 --коэф.из градусов в радианы
rse=rsi+hs+has delta = rsi-rre
rd=rre+delta/2 ts =2*Pi*rsi/Qs
tr=2*Pi*rre/Qru tsa =360/Qs
tra=360/Qru bts = ts-bs btr = tr-br
x=bs/2 y=hs+sqrt(rsi*rsi-x*x)
rsn=sqrt(x*x+y*y)
x=br/2 y=-hr+sqrt(rre*rre-x*x)
rrn=sqrt(x*x+y*y) bzs=2*Pi*rsn/Qs-bs
bzr=2*Pi*rrn/Qru-br qsp=Qs/(2*ms*p)
tp=Qs/(2*p) zs=ceil(qsp+tp*(1-bss))
closefile(f_d)--закрытие файла данных
--Создаем функцию для ввода кривых
--намагничивания В(Н) для сердечников
function get_B_H(bh_name, Steel_name, dsl, kFe)
--материал участка магнитопровода
mi_addmaterial(Steel_name, 0, 0, 0, 0, 0, dsl, 0, kFe,0,0)
--открытие файла с кривой В(Н)
B_H = openfile(bh_name .. ".txt","r")
str=read(B_H,"*l") –снятие заголовка
repeat --повтор ввода по условию
B,Н=read(В_Н,"*n","*n") --чтение B и Н
--добавление точки кривой В(Н) в FEMM
mi_addbhpoint(Stel_name,B,H)
until (B==0) --условие конца кривой
closefile(B_H) end --function get_B_H
--Ввод кривых намагничивания сталей
--ротора,ярма и зубцов статора
mvi=get_B_H(bhr,St_rot,la,kFer)
mvi=get_B_H(bhj,St_j_st,dsl,kFes)
mvi=get_B_H(bhz,St_z_st,dsl,kFes)
--Составное имя модели ТГ - тип fem
geom_mod="Fe_Md_" .. name_TG .. ".fem"
-- НАЧИНАЕТСЯ РАБОТА ПРОГРАММЫ FEMM
Start_date= date()
create(0) --Создаем магнитную Задачу
-- Параметры магнитной Задачи
mi_probdef(0, "millimeters", "planar", 1E-8, la)
--СОЗДАНИЕ ГЕОМЕТРИИ ТУРБОГЕНЕРАТОРА
mi_addnode(0,0) --ставим точку центра
--точки верха и низа пов-ости статора
mi_addnode(0,rse) mi_addnode(0,-rse)
--размер сегмента окружности (град.)
max_seg = 5 --далее может меняться
--левая полуокружность пов-сти статора
mi_addarc(0,rse,0,-rse,180,max_seg)
--правая полуокружность
mi_addarc(0,-rse,0,rse,180,max_seg)
--Создание зубцовового слоя статора
ISSN 2074-272X. Електротехніка і Електромеханіка. 2015. №4 11
--координатные размеры паза статора
ya=bs/2 yb=bks/2 b1=yb-ya
xa=sqrt(rsi*rsi-ya*ya) xb=xa+hss
xe=xa+hks xf=xa+(hks+hs)/2 xg=xa+hs
a1=b1*tan(as1) xc=xb+a1 –-1 скос паза
a2=b1*tan(as2) xd=xe-a2 – 2 скос паза
--ставятся точки от 1 до 7-ой (рис.2)
mi_addnode(xa,ya) mi_addnode(xb,ya)
mi_addnode(xc,yb)
mi_addnode(xd,yb) mi_addnode(xe,ya)
mi_addnode(xf,ya) mi_addnode(xg,ya)
--линии между точками 1,2,3,4,5,6,7
mi_addsegment(xa,ya,xb,ya)
mi_addsegment(xb,ya,xc,yb)
mi_addsegment(xc,yb,xd,yb)
mi_addsegment(xd,yb,xe,ya)
mi_addsegment(xe,ya,xf,ya)
mi_addsegment(xf,ya,xg,ya)
--размечаются точки от 1 до 7-ой
mi_selectnode(xa,ya) mi_selectnode(xb,ya)
mi_selectnode(xc,yb) mi_selectnode(xd,yb)
mi_selectnode(xe,ya) mi_selectnode(xf,ya)
mi_selectnode(xg,ya)
--размечаем линии между тчк. от 1 до 7
mi_selectsegment(xa,ya) y=(ya+yb)/2
x=(xb+xc)/2 mi_selectsegment(x,y)
x=(xc+xd)/2 mi_selectsegment(x,yb)
x=(xd+xe)/2 mi_selectsegment(x,y)
x=(xe+xf)/2 mi_selectsegment(x,ya)
x=(xf+xg)/2 mi_selectsegment(x,ya)
--назначаем выделенным точкам и линиям
--группу 1 и снимаем выделение
mi_setnodeprop(" ",1)
mi_setsegmentprop(" "," ", 1,0,1)
mi_clearselected() --см. рис.4
Рис. 4. Центр и поверхность сердечника статора с опорными
точками дуг и стенка паза статора
mi_selectgroup(1) --выдделяем стенку паза
--и отображаем симметрично оси х вниз
mi_mirror(0,0,rse,0,4) mi_clearselected()
--проводятся и размечаются линии 5-12,
-- 6-13 и 7-14,им назначается группа 1
mi_addsegment(xe,ya,xe,-ya)
mi_addsegment(xf,ya,xf,-ya)
mi_addsegment(xg,ya,xg,-ya)
mi_selectsegment(xe, 0) mi_selectsegment(xf, 0)
mi_selectsegment(xg,0)
mi_setsegmentprop(" "," ", 1, 0, 1)
mi_clearselected() --процесс на рис.5
Рис. 5. Процесс построения паза статора:
а – верхняя стенка паза; б – ее отражение относительно оси х;
в – добавили линии дна паза и слоев обмотки
--Выделяем паз и поворачиваем его на
--ползубцового шага вверх над осью х
mi_selectgroup(1) mi_moverotate(0,0,tsa/2,4)
--готовим координаты точек 8 и 14 и
--копируем их симметрично под ось х
u1=atan(ya/xa) u2=0.5*tsa*gr-u1
x8=rsi*cos(u2) y8=rsi*sin(u2)
u3=atan(ya/xg) u4=0.5*tsa*gr-u3
x14=rsn*cos(u4) y14=rsn*sin(u4)
mi_selectnode(x8,y8) mi_selectnode(x14,y14)
mi_mirror(0,0,rse,0,0)
--дуги зубца от точек 8 и 14 вниз
max_seg=1 u2=2*u2/gr
mi_addarc(x8,-y8,x8,y8,u2,max_seg)
max_seg=2 u4=2*u4/gr
mi_addarc(x14,-y14,x14,y14,u4,max_seg)
--выделяем дуги и назначаем группу 1
mi_selectarcsegment(rsi,0) mi_selectarcsegment(rsn,0)
mi_setarcsegmentprop(1," ", 0, 1)
mi_clearselected() --итог на рис.6
Рис. 6. Процесс построения зубца статора
--выделяем и копируем рис.6 (Qs-1)-раз
--и получаем зубцово-пазовую структуру
mi_selectgroup(1) mi_copyrotate(0,0,tsa,(Qs-1),4)
Аналогично создается геометрия ротора и объе-
диняется все в группу 9. Результат формирования
графической модели ТГ представлен на рис. 7.
12 ISSN 2074-272X. Електротехніка і Електромеханіка. 2015. №4
Рис. 7. Сформированная графическая модель ТГ
--Создаем свойства материалов ТГ
--магнитная проницаемость воздуха
mi_addmaterial("mu0",1,1)
--ЗАДАНИЕ ТОКОВ В ОБМОТКЕ СТАТОРА
Im=Is*2^0.5/as --амплитуда токов
--Фазные токи в обмотке статора
IA=Im*cos(beta*gr)
IB=Im*cos(beta*gr -2*Pi/3)
IC=Im*cos(beta*gr+2*Pi/3)
--Параметры обмотки статора
--радиусы ее нижнего и верхнего слоев
rn=rsi+0.75*hs rv=rsi+hs/4+hks
--числа пазов нижнего и верхнего слоев под осью q
qsn=0.5*(qsp-tp*(1-bss)) qsv=qsp-qsn
--углы расположения первых стержней
--нижнего и верхнего слоев фазы А
an1=180-(qsn-0.5)*tsa av1=180-(qsv-0.5)*tsa
--Функция задания тока в фазной зоне
--обмотки статора и материала стержней
--cir_name-имя зоны, i_s -ее ток,
--j_beg, j_end --номера начального и
--конечного проводников
--N_group -задаваемый номер группы
function tok_faz(cir_name, i_s, j_beg, j_end, N_group)
--связка имени тока с его значением
mi_addcircprop(cir_name, i_s, 1)
--цикл перебора проводн.фазной зоны
for j=j_beg,j_end,1 do ai=(j-1)*tsa
--координаты нижних стержней
ani=(an1+ai)*gr
x=rn*cos(ani) y=rn*sin(ani)
--добавление и выделение метки свойств
mi_addblocklabel(x, y) mi_selectlabel(x, y)
--задание свойств стержней
mi_setblockprop("mu0",0,fe4,cir_name,0,N_group,Ncs)
mi_clearselected()
--координаты верхних стержней
avi=(av1+ai)*gr x=rv*cos(avi) y=rv*sin(avi)
--добавление и выделение метки свойств
mi_addblocklabel(x, y) mi_selectlabel(x, y)
--задание свойств стержней обмотки
mi_setblockprop("mu0",0,fe4,cir_name,0,N_group,Ncs)
mi_clearselected()
end --конец перебора стержней по j
end --конец function tok_faz
--Задание токов в шести фазных зонах обмотки
--статора и отведение им групп 2 - 7
zap=tok_faz("IA+",IA,1,qsp,2)
zam=tok_faz("IA-",-IA,3*qsp+1,4*qsp,5)
Рис. 8. Метки свойств в стержнях фазных зон фазы А
zap=tok_faz("IB+",IB, 2*qsp+1,3*qsp,3)
zam=tok_faz("IB-",-IB,5*qsp+1,6*qsp,6)
zap=tok_faz("IC+",IC, 4*qsp+1,5*qsp,4)
zam=tok_faz("IC-",-IC, qsp+1,2*qsp,7)
Рис. 9. Добавлены метки свойств в стержнях фазных зон
фаз В и С обмотки статора
--Материал и токи в пазах ротора
--Связка имен тока с его значениями Ir
mi_addcircprop("IF-",-Ir, 1) mi_addcircprop("IF+",Ir, 1)
--добавление и выделение метки в пазу
y=rre-0.75*hr x=0.4*tr
mi_addblocklabel(y,x) mi_selectlabel(y,x)
--задание свойств стержней и группы 10
mi_setblockprop("mu0", 0,fe5, "IF-", 0, 10, Ncr)
Рис. 10. Метка свойств стержня в пазу ротора
--задание токов "IF-" по пазам справа
--вверху, отражение вниз
mi_copyrotate(0,0,tra,(Qr/4-1),2)
mi_selectgroup(10) mi_mirror(0,0,rre,0,2)
--перевод пазов из 10-ой в 9-ую группу
ISSN 2074-272X. Електротехніка і Електромеханіка. 2015. №4 13
mi_selectgroup(10)
mi_setblockprop("mu0", 0, fe5, "IF-", 0, 9, Ncr)
mi_clearselected()
--добавл. и выделение метки тока "IF+"
mi_addblocklabel(-y,x) mi_selectlabel(-y,x)
--задание свойств стержней, группа 10
mi_setblockprop("mu0", 0, fe5, "IF+", 0, 10, Ncr)
--задание токов "IF+" по пазам слева
--вверху, и отражение вниз
mi_copyrotate(0,0,-tra,(Qr/4-1),2)
mi_selectgroup(10) mi_mirror(-rre, 0, 0, 0, 2)
--перевод пазов из 10-ой в 9-ую группу
mi_selectgroup(10)
mi_setblockprop("mu0",0,fe5,"IF+",0,9,Ncr)
mi_clearselected()
Рис. 11. Формирование свойств стержней в пазах ротора
--Ставится и выделяется метка, задание
--магнитных свойств зазора, группа 0
mi_addblocklabel(0,(rre+rsi)/2)
mi_selectlabel(0,(rre+rsi)/2)
mi_setblockprop("mu0",0,fe1," ",0,0,0)
mi_clearselected()
--Создаем граничные условия Дирихле
mi_addboundprop("Grn",0,0,0,0,0,0,0,0)
--выделяем дуги поверхности статора
mi_selectarcsegment(rse,0) mi_selectarcsegment(-rse,0)
max_sg=5
--Назначаем поверхности ст-ра условия
mi_setarcsegmentprop(max_sg,"Grn",0,0)
mi_clearselected()
--СОЗДАЕМ ПАРАМЕТРЫ СЕРДЕЧНИКОВ
--Ставится и выделяется метка, задание
--магнитных свойств ярма статора, гр.0
y=rsi+hs+has/2 mi_addblocklabel(0,y)
mi_selectlabel(0,y)
mi_setblockprop(St_j_st, 0, fe2, " ", 0, 0, 0)
mi_clearselected()
Рис. 12. Метки для зазора и ярма сердечника статора
--Задание метки и материала для зубца
--статора на оси х, группа 1
x=rsi+hs/2 mi_addblocklabel(x,0)
mi_selectlabel(x,0)
mi_setblockprop(St_z_st, 0, fe3, " ", 0, 1, 0)
--распространение по всем зубцам
mi_copyrotate(0,0,tsa,(Qs-1),2)
--деление ширины фазной зоны zs на 2
z1=zs/2 z2=ceil(z1) –-дробное и целое
--зависимо от этого поворачиваем или
--нет зубцовую структура на ползубца
if z1==z2 then str="zs-четное и нет"
else str="zs-нечетное и есть поворот")
mi_selectgroup(1) mi_moverotate(0,0,tsa/2,4) end
Рис. 13. Метки и свойств зубцов сердечника статора
--Метка св-тв, ее выделение и задание
--материала для зубцов ротора, гр. 9
x=rre-hr/2 mi_addblocklabel(x,0) mi_selectlabel(x,0)
mi_setblockprop(St_rot, 0, fe6, " ", 0, 9, 0)
Рис. 14. Метка свойств для зубца сердечника ротора,
она распространяется и на прилегающую кольцевую зону
и через нее на все остальные зубцы ротора
--Задание материала для ярма сердечника ротора
y=rre/2 mi_addblocklabel(0,y) mi_selectlabel(0,y)
mi_setblockprop(St_rot, 0, fe7, " ", 0, 9, 0)
mi_clearselected()
mi_saveas(geom_mod) --Запомнить модель ТГ
mi_analyze(1) --Выполнить расчет МП
--Переход к показу результатов расчета
mi_loadsolution()
writeto() --Закрытие файла результатов
mi_close() --Закрывается документ
--препроцессора магнетизма
14 ISSN 2074-272X. Електротехніка і Електромеханіка. 2015. №4
Рис. 15. Полная готовность модели турбогенератора
Иинструкция по использованиюк скрипта.
1. Файл скрипта с расширением *.lua и все файлы
с исходными данными должны находиться в одной
папке на диске компьютера.
2. Запускаем программу FEMM и в меню File вы-
бираем Open Lua script и в диалоговом окне запускаем
файл с именем программы в работу.
3. Дальнейшие построения программа выполняет
автоматически. В программу можно вставить коман-
ды Pause для ее останова в конкретном месте и про-
смотра текущего изображения на экране компьютера.
Клавишей Enter продолжить выполнение программы.
4. Команды mi_saveas(name) в программе позволя-
ют сохранить очередной файл с текущим результатом
построенной модели ТГ под именем name.
Вывод. Представленный скрипт Lua дает широ-
кие возможности пользователям программы FEMM в
создании расчетных моделей для обеспечения расче-
тов магнитных полей и электромагнитных параметров
ТГ. Этот скрипт универсален с точки зрения геомет-
рических размеров и структуры электромагнитной
системы ТГ в рамках их распространенной типовой
конструкции. На базе рассмотренных принципов
формирования расчетных моделей могут быть разра-
ботаны аналогичные скрипты для других типов ЭМ
(например, уже создан и используется скрипт для
формирования моделей асинхронных двигателей).
СПИСОК ЛИТЕРАТУРЫ
1. Meeker D. Finite Element Method Magnetics. FEMM 4.2
11 Oct 2010 Self-Installing Executable. Режим доступа:
www.femm.info/wiki/OldVersions.
2. Титов В.В., Хуторецкий Г.М., Загородная Г.А. и др.
Турбогенераторы. – Л.: Энергия. – 1967. – 895 с.
3. Милых В.И., Полякова Н.В. Численно-полевая оценка
эффективности укорочения обмотки статора турбогенератора
// Електротехніка і електромеханіка. – 2014. – №4. – С. 12-16.
4. Милых В.И., Полякова Н.В. Система направлений и фазо-
вых соотношений электромагнитных величин при численных
расчетах магнитных полей в турбогенераторе //
Електротехніка і електромеханіка. – 2011. – №5. – С. 33-38.
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. Titov V.V., Hutoreckij G.M., Zagorodnaja G.A., Vartan'jan
G.P., Zaslavskij D.I., Smotrov I.A. Turbogeneratory [Tur-
bogenerators]. Leningrad, Energiia Publ., 1967. 895 p. (Rus).
3. Milykh V.I., Polyakova N.V. Numerical field estimation of
turbogenerator stator winding shortening efficiency. Elektro-
tekhnіka і elektromekhanіka – Electrical engineering & elec-
tromechanics, 2014, no.4, pp. 12-16. (Rus).
4. 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. Elektro-
tekhnіka і elektromekhanіka – Electrical engineering & elec-
tromechanics, 2011, no.5, pp. 33-38. (Rus).
Поступила (received) 20.04.2015
Милых Владимир Иванович1, д.т.н., проф.,
Полякова Наталия Владимировна1, инженер,
1 Национальный технический университет
«Харьковский политехнический институт»,
61002, Харьков, ул. Фрунзе, 21,
тел/phone +38 057 7076514, e-mail: mvikemkpi@gmail.com
V.I. Milykh1, N.V. Polyakova1
1 National Technical University «Kharkiv Polytechnic Institute»,
21, Frunze Str., Kharkiv, 61002, Ukraine.
Automated formation of calculation models of turbogenerators
for software environment FEMM.
Attention is paid to the popular FEMM (Finite Element Method
Magnetics) program which is effective in the numerical calcula-
tions of the magnetic fields of electrical machines. The main prob-
lem of its using - high costs in time on the formation of a graphical
model representing the design and on the formation of the physical
model representing the materials properties and the winding cur-
rents of machines – is solved. For this purpose, principles of the
automated formation of such models are developed and presented
on the turbogenerator example. The task is performed by a pro-
gram written in an algorithmic language Lua integrated into the
package FEMM. The program is universal in terms of varying the
geometry and dimensions of the designed turbogenerators. It uses a
minimum of input information in a digital form representing the
design of the whole turbogenerator and its fragments. A general
structure of the Lua script is provided, significant parts of its text,
the graphic results of work's phases, as well as explanations of the
program and instructions for its use are given. Performance capa-
bilities of the compiled Lua script are shown on the example of the
real 340 MW turbogenerator. References 4, figures 15.
Key words: program FEMM, Finite Element Method, tur-
bogenerator, graphical and physical models, automated
formation, Lua script.
|
| id | nasplib_isofts_kiev_ua-123456789-149284 |
| institution | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| issn | 2074-272X |
| language | Russian |
| last_indexed | 2025-12-07T18:49:12Z |
| publishDate | 2015 |
| publisher | Інститут технічних проблем магнетизму НАН України |
| record_format | dspace |
| spelling | Милых, В.И. Полякова, Н.В. 2019-02-19T20:21:54Z 2019-02-19T20:21:54Z 2015 Автоматизированное формирование расчетных моделей турбогенераторов для программной среды FEMM / В.И. Милых, Н.В. Полякова // Електротехніка і електромеханіка. — 2015. — № 4. — С. 7–14. — Бібліогр.: 4 назв. — рос. 2074-272X DOI: https://doi.org/10.20998/2074-272X.2015.4.02 https://nasplib.isofts.kiev.ua/handle/123456789/149284 621.313:681.3 Представлены принципы автоматизированного формирования графической и физической моделей турбогенераторов
 на алгоритмическом языке Lua для расчета их магнитных полей и параметров в программной среде пакета FEMM.
 Возможности составленного скрипта Lua демонстрируются на примере реального турбогенератора. Представлені принципи автоматизованого формування графічної і фізичної моделей турбогенераторів на
 алгоритмичній мові Lua для розрахунку їх магнітних полів і параметрів у програмному середовищі пакету FEMM.
 Можливості складеного скрипту Lua демонструються на прикладі реального турбогенератора. Attention is paid to the popular FEMM (Finite Element Method
 Magnetics) program which is effective in the numerical calculations of the magnetic fields of electrical machines. The main problem of its using - high costs in time on the formation of a graphical
 model representing the design and on the formation of the physical
 model representing the materials properties and the winding currents of machines – is solved. For this purpose, principles of the
 automated formation of such models are developed and presented
 on the turbogenerator example. The task is performed by a program written in an algorithmic language Lua integrated into the
 package FEMM. The program is universal in terms of varying the
 geometry and dimensions of the designed turbogenerators. It uses a
 minimum of input information in a digital form representing the
 design of the whole turbogenerator and its fragments. A general
 structure of the Lua script is provided, significant parts of its text,
 the graphic results of work's phases, as well as explanations of the
 program and instructions for its use are given. Performance capabilities of the compiled Lua script are shown on the example of the
 real 340 MW turbogenerator. ru Інститут технічних проблем магнетизму НАН України Електротехніка і електромеханіка Електричні машини та апарати Автоматизированное формирование расчетных моделей турбогенераторов для программной среды FEMM Automated formation of calculation models of turbogenerators for software environment FEMM Article published earlier |
| spellingShingle | Автоматизированное формирование расчетных моделей турбогенераторов для программной среды FEMM Милых, В.И. Полякова, Н.В. Електричні машини та апарати |
| title | Автоматизированное формирование расчетных моделей турбогенераторов для программной среды FEMM |
| title_alt | Automated formation of calculation models of turbogenerators for software environment FEMM |
| title_full | Автоматизированное формирование расчетных моделей турбогенераторов для программной среды FEMM |
| title_fullStr | Автоматизированное формирование расчетных моделей турбогенераторов для программной среды FEMM |
| title_full_unstemmed | Автоматизированное формирование расчетных моделей турбогенераторов для программной среды FEMM |
| title_short | Автоматизированное формирование расчетных моделей турбогенераторов для программной среды FEMM |
| title_sort | автоматизированное формирование расчетных моделей турбогенераторов для программной среды femm |
| topic | Електричні машини та апарати |
| topic_facet | Електричні машини та апарати |
| url | https://nasplib.isofts.kiev.ua/handle/123456789/149284 |
| work_keys_str_mv | AT milyhvi avtomatizirovannoeformirovanierasčetnyhmodeleiturbogeneratorovdlâprogrammnoisredyfemm AT polâkovanv avtomatizirovannoeformirovanierasčetnyhmodeleiturbogeneratorovdlâprogrammnoisredyfemm AT milyhvi automatedformationofcalculationmodelsofturbogeneratorsforsoftwareenvironmentfemm AT polâkovanv automatedformationofcalculationmodelsofturbogeneratorsforsoftwareenvironmentfemm |