УДК 621.313
Г.В. Пуйло, Е.П. Насыпаная, Д.М. Левин
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ЭЛЕКТРОМАГНИТНЫХ ХАРАКТЕРИСТИК ЭЛЕКТРОТЕХНИЧЕСКИХ СТАЛЕЙ
Обґрунтована необхідність математичного моделювання та автоматизованого формування аналітичних моделей характеристик електротехнічних сталей для дослідницького проектування, синтезу, аналізу та оптимізації електротехнічних пристроїв. Сформульовані вимоги до цих математичних моделей та запропоновано алгоритм їх автоматизованого формування на основі пологих апроксимуючих функцій з експоненціальнимядром, що дозволяє отримувати математичні моделі у вигляді багатократно диференціюємих функцій, які задовольняють заданим вимогам по точності апроксимації.
Обоснована необходимость математического моделирования и автоматизированного формирования аналитических моделей электромагнитных характеристик электротехнических сталей для исследовательского проектирования, синтеза, анализа и оптимизации электротехнических устройств. Сформулированы требования к этим математическим моделям и предложен алгоритм их автоматизированного формирования на основе пологих аппроксимирующих функций с экспоненциальным ядром, что позволяет получать математические модели в виде многократно дифференцируемых функций, которые удовлетворяют заданным требованиям по точности аппроксимации.
Современная электротехническая промышленность производит множество различных видов электротехнических устройств (ЭУ) (электромагнитных и электромеханотронных преобразователей энергии, электроизмерительных систем, электросварочных агрегатов, электрических аппаратов и т.д.), в которых используются самые разнообразные виды магнитных система (МС) из электротехнической стали (ЭТС). Анализ тенденций совершенствования электротехнических устройств, как одних из важнейших составляющих развития мировой энергетики и энергетики Украины [6, 10] показал, что основными их направлениями являются снижение потерь и расхода электротехнических материалов.
Наиболее дорогим и нежелательным видом потерь в ЭУ являются потери в МС, обусловленные, прежде всего, уровнем удельных потерь руд в используемой ЭТС. Жесткая необходимость снижения уровня потерь в МС потребовала от мировой электротехнической промышленности создания новых марок ЭТС с существенно улучшенными магнитными свойствами и со сниженными удельными потерями [11]. Разработаны методы создания ЭТС с заданными свойствами [5, 7]. Сейчас в мире производится несколько десятков марок ЭТС, которые применяются в производстве ЭУ.
Существенное снижение потерь энергии в МС и расхода ЭТС в ЭУ возможно только при эффективном использовании электротехнических и энергетических свойств используемой ЭТС.
Однако производители ЭТС приводят в каталожных данных только минимальную информацию об их характеристиках. Это, как правило, представленные в табличной форме для промышленной частоты / зависимости:
• характеристика намагничивания (зависимость индукции В от напряженности магнитного поля Н), т.е. В(Н);
• зависимость руд = у1(5);
• зависимость удельной намагничивающей мощности от индукции дуд = у2(8);
• толщина и состав ЭТС;
• вид электроизоляционного покрытия.
Следует отметить, что эти зависимости приводятся по экспериментальным данным заводов-производителей ЭТС [8, 11] ив них могут содержаться неизбежные неточности, обусловленные как технологическими отклонениями параметров ЭТС, так и погрешностями эксперимента, что и отражается соответствующими допусками на значения каталожных данных параметров ЭТС.
Однако при исследовательском проектировании, проектном синтезе, анализе и оптимизации ЭУ необходим достаточно большой объем дополнительной информации об электромагнитных свойствах применяемой марки ЭТС, а во многих случаях и о ферромагнитных конструкционных материалах, применяемых в конструктивных компонентах ЭУ.
Кроме того, эта информация должна быть представлена в виде математического описания (математической модели). Это необходимо для возможности моделирования распределения и расчета магнитного поля в МС, определения нагрева МС, удельных потерь и намагничивающей мощности при любых (не только табличных) значениях магнитной индукции и частоты магнитного поля и т.д. В ряде эффективных методов оптимизации и расчетов магнитного поля, переходных процессов в ЭУ, требуется использование непрерывных и дифференцируемых функций магнитной проницаемости ц в зависимости от магнитной индукции, а также ее производных ц' = Л(В), ц" = /2(В), зависимости Н(В) и.т.д.
В связи с ростом количества новых видов ЭТС, перспективой их дальнейшего совершенствования и необходимостью эффективного использования их электромагнитных свойств проблема совершенствования методов и средств адекватного математического моделирования характеристик ЭТС становится все более актуальной, а вследствие своей трудоемкости обуславливает необходимость создания в системах автоматизированного проектирования ЭУ подсистем автоматизированного математического моделирования электромагнитных свойств ЭТС (ПАМЭСЭС). Основное назначение такой подсистемы - автоматизированное формирование на основе каталожных данных (или данных производителя) и хранение в базе данных адекватных математических моделей характеристик ЭТС, описывающих их необходимые электромагнитные и другие свойства ЭТС при различной частоте магнитного поля, с учетом физических свойств стали, состава, технологических особенностей производства, процесса обработки, прокатки, направления намагниченности, вида и способа покрытия и т.д.
Задача формирования математического описания характеристик ЭТС в настоящее время решается на основе применения различных аппроксимационных и интерполяционных методов для аппроксимации пре-^де всего основной кривой намагничивания В(Н). Для этого предлагаются (и реализуются) два подхода:
• интерполяционный [1, 2, 5], при котором для каждого сегмента напряженности Не[Н,-, Н+1] кривую В(Н) интерполируют кубическим сплайном. Этот метод отличается простотой и вычислительной эффективностью, позволяет вполне удовлетворительно решать задачу аппроксимации исходной характеристики В(Н) и ее первой производной, но возможности кубического сплайна не позволяют смоделировать ее вторую производную как гладкую непрерывную функцию и производные высших порядков;
• аналитический, с помощью одной формулы [3,
4, 9], когда аппроксимацию выполняют на всем диапазоне напряженности Не [Нтш, Нтах].
В работе [4] решение выполняется на основе предложенной гиперболической функции:
В = (А + В(1 + Б ■ Н )с )Ш(Б • Н), (1)
где А, В, С, Б - коэффициенты, которые определяются методом наименьших квадратов. Этот подход дает хороший результат для изотропных сталей, но для анизотропных ЭТС погрешности существенно возрастают.
В работе [9] вместо кривой В(Н) предложена универсальная аппроксимация зависимости магнитной проницаемости от напряженности в безразмерном виде ц* = у(Н):
М* (H) =
(
^ 1 ^max _ 1
1
l-^min - 0,17H
*(1-1,5# )
BSH
4b—l)
--1
У1 “ -1
H *-1 b 1 +
У ^max
(2)
B, к (H *)
max
где ц тах - максимальное значение относительной магнитной проницаемости для данного материала; Ц тш - минимальное значение относительной магнитной проницаемости для данного материала; Вх - индукция насыщения; Вц*тах - значение индукции, при котором ц максимальна; Ь - безразмерный коэффициент, обусловленный технологическими особенностями процесса прокатки; к(Н) - поправочный коэффициент, который равен
к (H *) = (і + 400exp(- 9 H *2 )^1 + -
1,2
2 + H
0,9
(3)
Достоинством формулы (2) является хорошая точность и то, что для ее использования необходимы только значения ц тш, ц тах, В„ Вц*тах, которые вычисляются по таблице исходных данных для В(Н). Однако при Н = 1 при определенных сочетаниях исходных данных производная от (2) терпит разрыв второго рода, т.е. не позволяет вычислить производные второго рода. Указанные проблемы исключаются при использовании аппроксимации с помощью пологих функций с экспоненциальным ядром, что позволяет получать семейство непрерывных дифференцируемых функций В(Н) и Н(В).
Идея такой аппроксимации состоит в замене переменной искомой функции, (в данном случае Н(В)) на новую у = у1(Я), которая, во-первых, нивелировала бы крутизну исходной кривой В(Н), во-вторых нормировала бы диапазон изменения аргумента к единичному отрезку, т.е. ує[0, 1].
В итоге, задача поиска приближения кривой В(Н), Нє[Нтіп, Нтах], сводится к поиску приближения кривой В(у) = ф(у), ує[0, 1].
Такую переменную логично назвать функцией
нормирования аргумента. Простейшая функция, отвечающая требованию нормирования, это линейная функция вида:
^лин(Я) = / ~ЯГ . (4)
Н max ^ Н min
Но линейная функция хоть и преобразует интервал He[Hmm, Hmax] к единичному виду уе[0, 1] но никак не влияет на форму исходной кривой -В(у), а просто на отрезке уе[0, 1] воспроизводит ее форму. Для изменения же формы (с целью снижения ее крутизны), целесообразно воспользоваться логарифмической функцией нормирования:
^ (Я) = 1П(1 + Н ~Нmin ) . (5)
^ЛОГ' ' 1n(1 + Н max - ffmin)
Функция (5) нивелирует крутизну практически любой кривой намагничивания, делая ее более пригодной для визуального анализа и подбора аппроксимирующих функций. К тому же очень хорошо отображаются оба колена кривой намагничивания (рис. 1).
Рис. 1. Характеристика намагничивания для стали 3404 0,35 мм: а) исходной В(Н), б) преобразованной В(у)
По подобию преобразованной кривой 5(у), была предложена пологая экспоненциальная функция с постоянным масштабом по оси У, следующего вида
ф(х) = Му ■ ехр[-(Мх(х -хзт))2п]; X £ [0,1], (6)
где Му - параметр, который отвечает за масштаб (высоту) функции по оси У; Мх - параметр, который отвечает за масштаб по оси X, с уменьшением которого график функции растягивается вдоль оси X; хт -смещение кривой по оси X (для подбора колена кривизны под исходный график); п - степень крутизны кривой, (целочисленный параметр, при увеличении которого возрастает и крутизна графика функции на обоих коленах).
Функция (6) может быть представлена в более универсальном виде
ф(х) = Мух (х) ■ Му ■ Е(х); х е[0,1] , (7)
где Мух(х) - функция искажения амплитуды по оси У; Е(х) = ехр[-(Мх(х - х5т))2п] - ядро аппроксимирующей функции, определяющее пологость, крутизну и форму кривой в целом.
4
Процесс формирования математической модели (ММ) по таблично заданной функции В(Н) представлен виде блок-схемы на рис. 2.
Как видно из графика функции ошибок Еь для подотрезка уе[0; ~0,23], ее значения близки нулю. Таким образом, полученная функция на данном по-дотрезке приближает исходную кривую с абсолютной погрешностью не более 4 %.
2.2. Поскольку функция ошибок Е1 при у>0.23 не равна нулю, и погрешность аппроксимации составляет десятки процентов, то необходимо подбирать вторую аппроксимирующую функцию ф2(у). Аппроксимирующая функция ф2(у) также будет подбираться, под график функции ошибок Еь Для данного шаблона, ориентируясь по графику Е^у), (рис. 4,а) подбираем соответствующие параметры Мх = 0.75,
1.394, х,т = 1.637, п = 14, М = 0.5189, М2 = 1.
Му
Рис. 2. Блок-схема процесса формирования ММ по исходной зависимости В(Н)
Пример аппроксимации кривой намагничивания для стали 3404 толщиной 0.35 мм
1. Замена переменной.
На рис. 1,а изображена таблично заданная кривая В(Н), подлежащая аппроксимации. Ввиду значительной нелинейности исходной кривой произведем замену переменной аргумента Н. В результате получим более удобную для аппроксимации кривую В(у) (рис.
1.6), на которой хорошо отображены как начальный участок зависимости В(Н), так и участок насыщения.
2. Подбор аппроксимирующих функций.
Аппроксимирующие функции фг(у) подбираются
не для всей зависимости, а лишь для ее конкретного, нелинейного, участка, например, начального. Следует отметить, что результирующая аппроксимирующая функция определена по всему рабочему диапазону аргумента у. Подбор аппроксимирующих функций выполняется в специальном редакторе. Число необходимых функций зависит от характера кривой и требований точности.
2.1. Подбор первой аппроксимирующей функции выполняется под начальный участок исходной преобразованной кривой В(у) для подотрезка уе[0; ~0,23] (рис. 3,а). Для этого в редакторе выбирается экспоненциальный шаблон с постоянным коэффициентом амплитудной коррекции и его параметры Мх = 1, Му = 0.6, xsm = 1.03, п = 2 настраиваются на совпадение кривой В(у) на подотрезке уе[0; ~0,23]. При настройке данных параметров на отдельном графике (рис. 3,а) отображаются как исходная кривая В(у) так и ее аппроксимирующая функция ф1(у). При необходимости аппроксимирующий участок может быть увеличен на все поле просмотра. В том же редакторе предусмотрено поле просмотра функции ошибок (рис.
3.6), в данном случае это функция Е1 = В(у) - ф1(у).
а)
б)
Рис. 3. Подбор первой аппроксимирующей функции: а) поле просмотра функций, б) поле просмотра функции ошибок Е1 = В(^) - ф1(^)
а)
Рис. 4. Подбор второй аппроксимирующей функции для Е1: а) поле просмотра функций, б) поле просмотра функции ошибок Е2
В результате мы получим функцию
(р2 = (0.7233 + 0.6707^)е“3 1750'10-4(^“1 6370)28 , (8)
которая приближает функцию ошибок Е1 на всем рабочем отрезке с относительной погрешностью не хуже
0.025 (<0.5 %). Новая функция ошибок Е2 = Е1 - ф2 представлена на рис. 4,6, из которого видно, что она содержит только ошибки измерений и дальнейшая аппроксимация не имеет смысла.
Результирующая аппроксимирующая кривая выразится суммой аппроксимирующих функций ф1 и ф2:
ф(^) = В(^) = д>1 +Ф2 = 0,6e_(^_1,030)4 + (9)
+ (0,7233 + 0,6707^)e_3,1750'10_4(^_1,6370)28.
В редакторе аппроксимирующих функций, для полученной ММ В(у) можно также просмотреть график ошибок в процентах от таблично заданных точек (рис. 5), для всего отрезка аргумента у, а следовательно, и Н.
і 1 -1.1
ТЯ "
І. .
□ п Iff г - - ■■■ 0 QD 1 ■ * ■' "
G
1 щ
1
о
і
0,2 0,4 0,6 0,8
Рис. 5. График распределения погрешностей (в %) по исходным точкам для функции ф(у)
В результате изначально таблично заданную кривую B(H) для любой точки, рабочего отрезка напряженности магнитного поля Hє[7;7000] можно построить по двум формулам
B(^) = 0,60е_(^“1,030)4 +
(10)
+ (0,7233 + 0,6707^)є-3,175°'1(г4^-1,6370)28 , где у = 0,1130-^^ - 6,0) - функция логарифмического нормирования.
Соответственно по этой ММ, могут быть рассчитаны, все остальные зависимости, такие как B'(H), B"(H), H(B), H,(B), H"(B).
СПИСОК ЛИТЕРАТУРЫ
1. Борзік В.Л. Апроксимація кривої намагнічування сталі в суміщених електричних машинах / В.Л. Борзік // Електротехніка і електромеханіка. - 2003. - № 2. - С. 11-12.
2. Гололобов В.В. Представление кривых намагничивания в численных расчетах магнитного поля / В.В. Гололобов, В.В. Рымша, И.В. Меркулов, А.С. Порайко // Електромашинобудування та електрообладнання. - 1999. - № 5. - С.81-85.
3. Гольдштейн В.Г. Функциональные гиперболические и полиномиальные аппроксимации кривых намагничивания / В.Г. Гольдштейн, В.М. Мякишев, М.С. Жеваев // Моделирование и оптимизация динамических систем и систем с распределёнными параметрами, матем. моделирование и краев. задачи: труды IV Всерос. науч. конф. с междунар. участием, 29-31 мая 2007 г.: тезисы докл. 4.2. - 2007. - С. 44-48.
4. Зубков Ю.В. Аппроксимация кривых намагничивания электротехнической стали при проектировании совмещенных электромеханических преобразователей / Ю. В. Зубков, В. Е. Верещагин // Известия вузов. Электромеханика. -2007. - № 6. - С. 19-22.
5. Зубрицкий С.М. Аппроксимация основной кривой намагничивания магнигомягких материалов / С.М. Зубрицкий // -Спец. вып.: Материалы ежегодн. науч.-теоретич. конф. молодых ученых. - Вестн. Иркут. ун-та. - Иркутск, 2002. - С. 28-29.
6. Інноваційні пріоритети паливно-енергетичного комплексу України // під заг. ред. А.К. Шидловського. - К.: Укр. енциклопедичні знання, 2005. -512 с.
7. Магнитные свойства электротехнических сталей и сплавов / [Казаджан Л.Б.]; под ред. В. Д. Дурнева. - М.: Наука и технологии, 2000. - 224 с.
В. Молотилов Б.В. Холоднокатаные электротехнические стали / Молотилов Б.В. - М.: Металлургия, 19В9. - 16В с.
9. Пентегов И.В. Универсальная аппроксимация кривых намагничивания электротехнических сталей / И.В. Пентегов, А.В. Красножон // Електротехніка і електромеханіка. -2006. - № 1. - С. 66-70.
10. Пуйло Г. В. Современные тенденции совершенствования силовых трансформаторов / Г. В. Пуйло, И. С. Кузьменко, В. В. Тонгалюк // Електротехніка та електромеханіка. -200В. -№ 2.-С. 4В-З2.
11. Прокат холоднокатаной тонколистовой электротехнической стали. Каталог продукции. ОАО "Новолипецкий металлургический комбинат" [Электронный ресурс] // Режим доступа: http://www.nlmk.ru.
Bibliography (transliterated): 1. Borzik V.L. Aproksimaciya krivoi namagnichuvannya stali v sumischenih elektrichnih mashinah / V.L. Borzik / Elektrotehnika і elektromehanika. - 2003. - № 2. - S. 11-12. 2. Gololobov V.V. Predstavlenie krivyh namagnichivaniya v chislennyh raschetah magnitnogo polya / V.V. Gololobov, V.V. Rymsha, I.V. Merku-lov, A.S. Porajko // Elektromashinobuduvannya ta elektroobladnannya. -1999. - № З. - S^l^. 3. Gol'dshtejn V.G. Funkcional'nye giperboliches-kie i polinomial'nye approksimacii krivyh namagnichivaniya / V.G. Gol'd-shtejn, V.M. Myakishev, M.S. Zhevaev // Modelirovanie i optimizaciya dinamicheskih sistem i sistem s raspredelennymi parametrami, matem. modelirovanie i kraev. zadachi: trudy IV Vseros. nauch. konf. s mezhdunar. uchastiem, 29-31 maya 2007 g.: tezisy dokl. Ch.2. - 2007. - S. 44-4В. 4. Zubkov Yu.V. Approksimaciya krivyh namagnichivaniya 'elektrotehnicheskoj stali pri proektirovanii sovmeschennyh 'elektrome-hanicheskih preobrazovatelej / Yu. V. Zubkov, V. E. Vereschagin // Izvestiya vuzov. 'Elektromehanika. - 2007. - № 6. - S. 19-22. 5. Zu-brickij S.M. Approksimaciya osnovnoj krivoj namagnichivaniya magni-tomyagkih materialov / S.M. Zubrickij / - Spec. vyp.: Materialy ezhe-godn. nauch.-teoretich. konf. molodyh uchenyh. - Vestn. Irkut. un-ta. -Irkutsk, 2002. - S. 2В-29. 6. Innovacijni prioriteti palivno-
energetichnogo kompleksu Ukraini / pid zag. red. A.K. Shidlovskogo. -K.: Ukr. enciklopedichni znannya, 200З. -З12 s. 7. Magnitnye svojstva 'elektrotehnicheskih stalej i splavov / [Kazadzhan L.B.]; pod red. V. D. Durneva. - M.: Nauka i tehnologii, 2000. - 224 s. 8. Molotilov B.V. Holodnokatanye 'elektrotehnicheskie stali / Molotilov B.V. M.: Metal -lurgiya, 19В9. - 16В s. 9. Pentegov I.V. Universal'naya approksimaciya krivyh namagnichivaniya 'elektrotehnicheskih stalej / I.V. Pentegov, A.V. Krasnozhon / - Elektrotehnika і elektromehanika. 2006. № 1. <http:lllibrary.kpi.kharkov.ualJUR/eie_2006-1.pdf>- S. 66-70. 10. Pujlo G.V. Sovremennye tendencii sovershenstvovaniya silovyh transformatorov / G. V. Pujlo, I. S. Kuz'menko, V. V. Tongalyuk / Elektrotehnika ta elektromehanika. - 200В. - № 2.- S. 4В-З2. 11. Prokat holodnokatanoj tonkolistovoj 'elektrotehnicheskoj stali. Katalog produkcii. OAO "Novolipeckij metallurgicheskij kombinat" ['Elektronnyj resurs] // Rezhim dostupa: <http://www.nlmk.ru>.
Поступила 03.09.2011
Пуйло Глеб Васильевич, д.т.н., проф.,
Насыпаная Елена Петровна,
ЛевинДмитрий Михайлович
Одесский национальный политехнический университет
кафедра электрических машин
6З044, Одесса, пр. Шевченко, 1
тел. (04В2) 734В494, e-mail: KeM@opu.ua
Puilo G. V., Nasypanaya E.P., Levin DM.
Mathematical simulation of electromagnetic characteristics of electric steels.
The necessity of mathematical simulation and computer-assisted forming of analytical models for electric steels electromagnetic characteristics with the purpose of research designing, synthesis, analysis and optimization of electrical devices is grounded. Requirements for these mathematical models are formulated, an algorithm of their computer-assisted forming on the basis of flat approximating functions with an exponential kernel developed, which allows obtaining mathematical models in the form of repeatedly differentiable functions that satisfy specified requirements for the approximation accuracy.
Key words - subsystem, mathematical models, electric steel, approximating function.