CC BY
50Теоретические аспекты конечно-элементного моделирования состояний
и функциональности магнитоуправляемых датчиков
Светлана СЫСОЕВА
В статье рассматриваются аспекты математического моделирования магнитостатических, магнитодинамических и магнитотепловых (мультифи-зических) процессов, ответственных за функциональность магнитоуправляемых датчиков скорости и положения, и параллельно раскрываются моделирующие возможности и особенности конечно-элементного CAD/CAE-пакета FLUX Cedrat.
Введение
Магнитоуправляемые датчики скорости/положения представляют собой широко распространенный и перспективный класс сенсорных устройств в составе автомобильных, промышленных, аэрокосмических и других систем. Эти устройства отличаются бесконтактностью, готовым типом выхода, часто цифровым исполнением, достаточно высокой чувствительностью, стабильностью, надежностью и точностью. Чтобы разработать и гарантировать работу такого устройства в полном диапазоне рабочих условий, включая в анализ материалы и геометрию магнитной системы, допуски, рабочие температуры и другие параметры, стоит воспользоваться компьютерным конечно-элементным моделированием, которое способно дать хорошие результаты. Пример такого специализированного программного обеспечения— Flux Cedrat [1-9] (рис. 1).
Существуют три основных модели и/или состояния, которые характеризуют функциональность магнитоуправляемых датчиков, соответствующих приложениям (модулям) Flux:
1) Магнитостатическое (магнитное поле постоянное, источники поля — стационарные токи проводимости или постоянные стационарные магниты, в системе отсутствуют механическое движение источника или модулятора поля и наведенные токи), изучается в Flux MagnetoStatic.
2) Магнитогармоническое (магнитное поле — установившееся переменное синусоидальное, источники поля — переменные токи, подвижные постоянные магниты и/или неподвижные постоянные магниты, стационарное поле которых синусоидально модулиру-
ется подвижным ферромагнитным/проводящим модулятором, в расчет принимаются наведенные токи, скин-эффект и эффект близости), соответствует приложению Steady state AC Magnetic Flux, но данный модуль не рассматривает магниты в качестве источников.
3)Переходные, соответствующие переходу системы из одного установившегося магнитостатического или магнитодинамического состояния в другое и характеризующиеся переменным несинусоидальным состоянием источников и полей. Во Flux пе-
Рис. 1. Типичная задача, решаемая методом конечных элементов, — моделирование датчика ферромагнитного зубчатого ротора:
а) постановка задачи: 1 — ферромагнитный зубчатый ротор; 2 — основные компоненты статора;
3 — обратно-смещающий магнит; 4 — ИС датчика (дифференциальная ИС Холла TLE4921 Infineon); б—г) решение в рамках реализации конечно-элементного проекта во Flux v.10.1.2:
б) результат генерирования сетки для линий и контуров;
в) результат генерирования сетки для элементов второго порядка Generate Second Order Elements;
г) постобработка: изображение векторных результатов решения образцовой модели
реходные процессы изучаются посредством специализированного модуля Transient Magnetic, допускающего задание в качестве источников и переменных (несинусоидальных) токов, и магнитов, анализ наведенных токов, скин-эффекта и эффекта близости. Магнитотепловое моделирование, также доступное во Flux, основано на существовании взаимодействия между электромагнитными и тепловыми аспектами процессов нагрева — вследствие вихревых токов.
Теоретические аспекты моделирования основных функциональных процессов, протекающих в датчиках (магнитных и магнитнотепловых), достаточно подробно рассматриваются в рамках данной статьи. Но нельзя не отметить, что возможности моделирования, предлагаемые Flux, шире. Так, возможность объединения с электрической схемой “circuit coupling” доступна пользователю в приложениях Transient Magnetic и Steady State AC Magnetic. Кинематический модуль Flux также допускает исследование подвижной части устройства вследствие механических сил (действия пружин, трения, гравитации).
Магнитостатическое
моделирование
Магнитостатическое моделирование основывается на фундаментальных законах электромагнетизма — уравнениях Максвелла и конституционных соотношениях среды. В самом общем виде уравнения Максвелла (в дифференциальной форме) записываются как:
ro_E =_>-dB(pt, rotH = j+(3D/3t), (1)
divB = 0, divD = q,
где B — плотность магнитного потока, j — плотность тока, H — напряженность магнитного поля, D и E — электрическое смещение и напряженность электрического поля соответственно, q — плотность статического заряда.
Соотношения для материалов записываются следующим образом:
I = ц
B = цЯ, (2)
D = гБ,
где а — проводимость; ц — магнитная проницаемость; s — диэлектрическая проницаемость.
В случае низкочастотных квазистатичес-ких электромагнитных полей выражения для электрических полей D и Б, B и H могут быть разъединены и рассмотрены отдельно.
Для чисто магнитной системы принимается гипотеза квазистатического состояния, вариацией плотности электрического потока D можно пренебречь. Уравнения Максвелла тогда запишутся как
rotE =—dB/dt, rotH = j, divB = 0,
(3)
а соотношения материалов, которые принимаются в расчет, так:
j = aE ’
B = цН.
(4)
rotH = j , divB = 0,
(5)
B = цН.
(6)
_B = цгц0Н, _ B = цг ц0Н +Br,
(7)
Н = -grad ф+Г,
Некоторые геометрические аспекты конструирования магнитных систем В соответствии с законом Ампера в электромагнитной схеме магнитное поле электрически возбуждаемой магнитной системы (рис. 2) определяется соотношением
jHdl = Ixn,
(11)
где I — длина магнитного пути; п -витков возбуждающей обмотки.
Условия вычисления в магнитостатическом исследовании
Переменные состояния не зависят от времени, то есть dl dt = 0, вычисляются только B-и H-поля, D- и Е-поля не вычисляются.
Наведенные токи не учитываются, и уравнения Максвелла для задач магнитостатики тогда записываются как
а соотношение материалов, которое принимается в расчет (только для связи магнитных параметров):
Это принципиальное соотношение для магнитных материалов может быть помещено в расчетную модель в виде В(Н)- или Н(В)-соотношений:
Рис. 2. Замкнутые электромагнитные схемы:
а) электромагнитно возбуждаемая;
б) магнитная система (возбуждаемая постоянным магнитом)
Для магнитной системы, возбуждаемой постоянным магнитом:
jHxdl = 0
где цг — относительная магнитная проницаемость среды; ц0 — абсолютная магнитная проницаемость вакуума; Вг — остаточный магнитный поток постоянных магнитов, или
HBlB+ НФ1Ф+ HMlM = 0
(12)
(13)
_Н = vrv0B __ Н = vrv0B -Hc,
где ут — относительное удельное магнитное сопротивление среды; у0 — абсолютное удельное магнитное сопротивление вакуума; Ис — коэрцитивная сила постоянных магнитов.
Закон Гаусса для магнетизма divB = 0 также допускает введение магнитного векторного потенциала А:
B = rotA.
Выражение для закона Ампера тоШ = ) позволяет дать выражение для магнитного скалярного потенциала:
где НВ — напряженность магнитного поля
(8) в воздухе; 1В — длина пути линий поля в воздухе; НФ — магнитное поле в ферромагнитном материале; 1Ф — длина пути линий поля в ферромагнитном материале; Нм — магнитное поле в магнитном материале; 1М — длина пути линий поля в магнитном материале.
При этом напряженности поля в ферромагнитной и воздушной частях замкнутой схемы будут характеризоваться тем же самым направлением, в магнитной — противоположным.
(9) Согласно этой записи, увеличения магнитного поля в воздушном зазоре можно достичь, уменьшая воздушный зазор и увеличивая длину магнитно-ферромагнитного пути.
В пределах пути 1М постоянного магнита напряженность и индукция магнитного по-(10) ля связаны соотношением:
где ф — магнитный скалярный потенциал; Т — вращательный или нулевой компонент поля источника или электрический векторный потенциал в А/м.
BM = цоН+Мм>
(14)
где Ым — намагничивание постоянного магнита, направление которого определяет на-
или
правление поля в контуре из ферромагнитной и воздушной частей схемы.
Поскольку и электрически, и магнитно возбуждаемые схемы бесконтактных датчиков конструируются так, чтобы выполнялось равенство
ВМ = ВВ = Вф, переписав (13) как
(вв1в/Ио)+(вф1ф/ИгИо)+
+ ((вм-Мм)/Ио) 1м = 0
и, далее, как
Вв1в+(Вв1ф/^г)+(Вв-Мм)1м = 0 и 1в+(1ф/^г) + 1м-(Мм1м/ВВ) = 0,
в итоге получим, что индукция в воздухе магнитной системы может рассчитываться по формуле:
ВВ = ММ/(1 + (1В/1М) + (1ф/Иг1М))- (15)
Из этой формулы явствует, что для увеличения поля в воздухе в схеме с постоянным магнитом должно выполняться соотношение:
1М ^ В (16)
Соотношение М > 1ф в датчиках можно на практике не принимать в расчет, так как 1ф/(ИгхМ) — всегда малое значение из-за того, что цг, магнитная проницаемость стали, обычно >1000. Тем не менее, ферромагнитную часть все же не следует делать превышающей по длине магнитную.
Однако длина ферромагнитного пути в сравнении с длиной воздушного пути — и в том случае, если магнитный участок отсутствует, и во всех остальных — должна быть значительно увеличенной, то есть
1ф ^ В 1ф+ 1М ^ В
В электромагнитной схеме возбуждения, в которой магнитное поле в замкнутом контуре будет определяться выражением
Нв1в+Нф1ф = п1 (17)
или сведенным к виду
(5вВ/Ио)+(ВфУИгИо) = п1, (18)
считая, что поле в каждой точке контура одинаково, получим, что магнитная индукция в зазоре будет также значительно зависеть от величины зазора В:
Вв = п1цо/(В+(УИг)). (19)
В то же время, если удается выполнять условие (16), достаточно большие зазоры — порядка нескольких мм — могут не оказывать
влияния на величину поля в воздухе (хотя она и уменьшается по сравнению с реальным значением ВФ) и не влиять на точность измерений.
Еще один важный практический аспект для геометрического конструирования датчиков — это возможность расширения воздушного зазора за счет увеличения поперечного сечения магнита или электромагнитно возбуждаемой ферромагнитной части.
Поскольку магнитный поток в замкнутом контуре принимается равным, то есть
BMSM = BBSB, (20)
это далее дает следующее выражение для индукции в воздушном зазоре (при идеальных условиях):
Вв1в+(вв1ф/цг) + (вв(^в/^м)-мм)1м = 0,
Вв = ЗД^мЖВМЖУит(21)
и при соблюдении условия (16) и большом значении цг выражение (21) аппроксимируется к виду:
BB = Mm/(SB/ SM) = (MMSM)/SB, (22)
то есть увеличение индукции определяется также фактором SM/SB.
Существуют и могут быть разработаны еще более развернутые аналитические методы технической оценки полей датчиков — с учетом рабочих точек кривых намагничивания магнитов и разнообразных факторов утечки. Но аналитические модели прежде всего полезны для упрощенного решения какой-то частной прикладной задачи. В рамках данной статьи подробно обсуждаются именно универсальные конечно-элементные, а не упрощенные аналитические модели для конкретных и идеальных условий и набора принятых допущений.
Магнитостатическое моделирование во Flux
Во Flux существует приложение Magneto Static, позволяющее получать следующие принципиальные результаты:
• распределение напряженности магнитного поля H, плотности магнитного потока B, плотности магнитной энергии dWm;
• расчет магнитной энергии в регионе Wm во всем домене исследования, силы Fm и магнитного крутящего момента Tm;
• магнитный поток, пересекающий поверхность или навивку.
Основные этапы магнитостатического исследования включают:
• магнитостатическое исследование образцовой модели, решение опорного сценария;
• физическое параметризованное исследование;
• геометрическое параметризованное исследование.
Для решения используются две модели:
• векторная модель с вычислением переменной состояния магнитного векторного потенциала A, введение которого допускает закон Гаусса для магнетизма;
• скалярная модель с вычислением магнитного скалярного потенциала, записываемого как ф^ или фы.
2D-модуль Flux обеспечивает решение только векторной модели, SD-модуль — и векторной, и скалярной. По умолчанию 3D-решатель использует векторную модель для 2D-приложений и скалярную для 3D-прило-жений. Но в дальнейшем на деталях работы с Flux в рамках данной статьи внимание не акцентируется, отчасти в связи с предшествующей публикацией ряда практических руководств именно по магнитостатическому моделированию магнитоуправляемых датчиков во Flux [2-7].
Для векторной модели решаемое во Flux уравнение записывается как
rot(V)[Vr]rot(A)-Hc) = _S, (23)
где [vr] — тензор относительного удельного сопротивления среды; js — поверхностная плотность источника тока (в А/м2).
Для скалярной модели T выбирается различным образом, и этот выбор приводит к различным скалярным потенциалам и применению различных магнитных формулировок во Flux 3D. Каждая область задачи связывается со своей формулировкой, но в общем виде форма решаемого в магнитостатическом приложении уравнения записывается как
div([^]^(-grad (Ф)+Т)+Вг) = 0, (24)
где [Цг] — тензор относительной магнитной проницаемости среды; ф — магнитный скалярный потенциал (ф^ или ф^; T — векторная характеристика источника поля; Br — остаточный магнитный поток магнитов.
Существуют три основных ситуации, относящиеся к источникам:
• источников тока нет;
• источники не покрываются сеткой;
• источники покрываются сеткой. Соответствующие им уравнения представлены далее.
Тотальный скалярный потенциал ф^ используется при отсутствии токовых источников. Плотность тока j равна нулю. Напряженность поля выражается как
H = -grad (ф^), (25)
и решаемое уравнение представляет собой
div([^]Vo(-grad (<U-Br) = 0, (26)
ф0 является переменной состояния, относительно которой решается уравнение (26).
Редукционный магнитный скалярный потенциал, соответствующий Hj или фгащ., ис-
Таблица 1. Локальные и глобальные параметры постобработки магнитостатических результатов решения во Flux 2D/3D
Параметры Единицы измерения Объяснение
Локальные 2Э-параметры
Векторный потенциал: • для плоской задачи — Л2, • для осесимметричной задачи — гЛ2 Вб/м Вб
Напряженность магнитного поля, Н А/м H = vxrotAz
Плотность магнитного потока, В Тл B = MH
Магнитный поток, Ф (ФР1Р2 — магнитный поток через поверхность, неограниченную линиями, нормальными к домену исследования, передаваемый посредством точек Р1 и Р 2) Вб фР1Р2 = (AZP1AZf2)p> где p — толщина домена плоской задачи
Проницаемость, ^ Гн/м
Плотность тока источника, А/мм2
Лапласова сила (локальное значение), Рь Н/м3 Fl = J х B
Глобальные 2Э-параметры
Магнитная сила, Fm (Force) Н Р т = д^т1д1,1 — элементарное линейное перемещение
Магнитная сила, Fm (Magnetic Pressure) Н Интегрирование поверхностных сил (тензорный метод Максвелла)
Магнитный крутящий момент, Гт (Torque) Н/м Гт = 5^/58, 8 — элементарное угловое перемещение
Лапласова сила, FL (stress) Н = |/х ВВМ V
Магнитная энергия, Wm Дж В мт = = |(|НсВ )ВМ. V V 0 Для линейной изотропной гомогенной среды: В^т = (1/2)НВ = (1/2)цН2 = (1/2)(В2/ц)
Магнитная энергия Energy/Co-Energy Дж |(|В вН)вм V 0
Ток (источника), I (Current) А I = ¡;га$ $
Магнитный поток в области, Ф (Inductance/Flux) Вб е
Индуктивность, L (Inductance/Inductance) Гн L = (я‘Ф)/1
Локальные ЗЭ-параметры
Скалярные магнитные потенциалы, ф^, фтщ, Ф^То А/м
Напряженность магнитного поля, Н А/м Н = -рнВ(ф,„,), Н = -¡таВ(ф„dн)+Hi, Н = -¡гоВ (фГест0)+Во
Плотность магнитного потока, В Тл В = [ц]Н
Магнитное поле источнико в, не покрываемых сеткой, Н^ А/м Аналитическое или полуаналитическое вычисление с использованием закона Био-Савара-Лапласа
Остаточная плотность магнитного потока постоянных магнитов, Вг Тл
Коэрцитивная сила постоянных магнитов, Нс А/м
Проницаемость, ^ Гн/м
Относительная проницаемость,
Плотность поверхностной силы между материалами, Вр1 Н/м2 Вычисление на границах сред с различной магнитной проницаемостью
Плотность источника тока ]$ в сетчатых обмотках А/м2 Верификация данных
Плотность лапласовой силы ВРЬ в сетчатых обмотках Н/м3 Вр1 = /$*В
Объемная плотность магнитной энергии, Дж/м3 В сшт=¡НсВ
Глобальные ЗЭ-параметры
Магнитная сила, Fm (Force) Н F щ = dWm/dl
Магнитный крутящий момент, Гт (Torque) Н/м Г m = dWJde
Магнитная энергия, Wm Дж Wm = jdWmdv V
Магнитный поток через обмотку-проводник (stranded-типа), Фв Вб
пользуется, когда источники — не покрываемого сеткой типа. Во Flux это единичные объекты (entity), которые могут быть круговой, прямоугольной или более сложной формы. Это случай, когда T = H, и Hj вычисляется аналитически с использованием закона Био-Савара-Лапласа. Закон позволяет рассчитать магнитное поле в точке P, созданное проводником с током I, без присутствия магнитных материалов как:
• элементарное
dH(P) = (1/4n)Idfx(?/r!); (27)
• тотальное
H(P) = (1/4n)fldfx(?/r!). (28)
В приведенных выражениях j — это вектор между элементом dl проводника и точкой P.
Для поверхностного или объемного вычисления тока выражения (27) и (28) должны быть переписаны: Idl заменяется на jsdS и j VdV соответственно.
Магнитное поле, создаваемое источником, не покрываемым сеткой, записывается как Hj. Магнитное поле в магнитостатической системе тогда выражается как
H = -grnJ (^redH^^ (29)
а решенное уравнение представляет собой div([^r]^o(-grad(фredHj)+H j)—B Г) = 0, (30)
где фтещ — переменная состояния.
Кроме того, магнитное поле может быть наложенным (например, магнитное поле Земли).
Редукционный магнитный скалярный потенциал, соответствующий Т0 или ф^Т), используется, когда источники — сетчатого типа, для которых
rotT 0= j, (31)
но может использоваться и с несетчатыми источниками, для которых
rotT 0 = j s+ rotH j. (32)
Магнитное поле выражается как
H = -grad (ФredT0)+T0, (33)
а решенное уравнение представляет собой
dV([^rK(-grad^redT0)+T0)-Br) = 0, (34)
где ф^о и T0 — переменные состояния. Описание модели во Flux включает:
• определение границ и граничных условий;
• описание материальных сред (воздушных, магнитных/ферромагнитных материалов, источников тока);
• описание источников поля (постоянных магнитов, DC тока, наложенного поля).
Граничные условия задаются в большинстве случаев методом бесконечного бокса: программой автоматически назначается условие граничного нулевого магнитного потенциала. Альтернатива — задание нематериальных областей пользователем и определение направления поля на границах как нормального или тангенциального (по умолчанию поле назначается тангенциальным к границе). Если в модели присутствуют симметрии и периодичности, пользователь также определяет направление поля как тангенциальное или нормальное по отношению к плоскостям симметрии, циклические или антициклические условия по отношению к плоскостям периодичности, определенным в геометрическом контексте.
Объемы, поверхности, линии материальных областей позволяют моделировать во Flux материальные среды — воздушные (с цг = 1), магнитные (магнитомягкие и магнитотвердые — с цг) непроводящие, проводящие немагнитные (цг = 1). Для моделирования доступны тонкие области и филиформы (области с малым поперечным сечением: когда два размера значительно меньше третьего, магнитное поле является тангенциальным к линии моделирования филиформы).
Постоянные магниты как источники поля описываются посредством назначения областей магнитных материалов. При определении материала вводятся характеристики магнита, остаточная индукция, указывается направление намагниченности. Электрические
токи как источники в магнитостатическом применении считаются непрерывными, независимыми от времени, с однородно распределенной плотностью тока в поперечном сечении проводника. Электрические источники описываются либо посредством областей материала типа coil conductor, либо посредством entities coils. Дополнительной возможностью является задание наложенного поля как создаваемого магнитным диполем или известным магнитным полем — однородным или пространственно распределенным.
Практическое магнитостатическое моделирование во Flux на примере датчика ферромагнитного зубчатого ротора было показано в публикациях [2-7], к которым (а также к технической документации Cedrat) рекомендуем обращаться за выяснением деталей задания свойств материалов, указания направления осуществления постобработки результатов и т. п.
В случае замкнутой магнитной системы с высокой магнитной проницаемостью возникает также задача соединения различных магнитных областей, если интеграл в выражении (11) не равен нулю. Решением является искусственное моделирование зазора посредством единичного объекта Magnetic Cut, но в том случае, если замкнутая система рассечена плоскостью симметрии и периодичности, введение зазора не требуется.
Анализ результатов предполагает постобработку локальных и глобальных параметров, различных для 2D- и 3D-применений (табл. 1).
Моделирование магнитогармонических процессов
Магнитные процессы, достигшие установившегося гармонического состояния, характеризуются одной и той же частотой, но синусоидальным изменением во времени всех физических параметров в уравнениях Максвелла.
Для изучения магнитных процессов, находящихся в установившемся гармоническом состоянии с постоянной частотой во Flux существует применение Steady State AC Magnetic. Магнитное поле в этом применении анализируется вместе с вихревыми токами, наведенными в проводящих областях.
Учитываются скин-эффект и эффект близости. Скин-эффект состоит в том, что электрические переменные токи неоднородно распределяются в поперечном сечении проводника — с преобладанием концентрации токовых линий в прилегающей к внешней стороне поверхности проводника, и этот эффект становится более значительным по мере того, как частота f тока увеличивается. Также скин-эффект становится более заметным с высокой магнитной проницаемостью ц материала и электрической проводимостью а.
Параметры, ассоциируемые с электромагнитным полем, уменьшаются по мере увеличения глубины z проводника по экспоненциальному закону:
U(z) = USxe-z/5, (35)
где US — значение на поверхности; ю = 2nf— круговая частота; глубина скин-эффекта:
5 = V2/(®^ct). (36)
Эффект близости относится к явлению, которое возникает в системе из нескольких токонесущих проводников, когда магнитное поле, сгенерированное проводником, разрушает токовое распределение в других проводниках.
Результаты магнитогармонического моделирования — те же, что и для магнитостатического приложения, плюс расчет параметров вихревых токов и потерь мощности вследствие эффекта Джоуля. Все вычисленные для данной частоты параметры — синусоидальные.
В контексте Steady State AC Magnetic все параметры выражаются посредством их комплексных изображений. Такое представление действительно при соблюдении следующих гипотез:
• все материалы линейные;
• все источники синусоидальные.
Тем не менее, представление комплексных изображений может быть расширено для двух случаев:
• для нелинейных материалов — параметры являются несинусоидальными, но предусмотрена возможность задания так называемой эквивалентной В(Н)-характеристи-ки, допускающей принимать в расчет факт синусоидальности параметров;
• для линейных материалов, но несинусоидальных источников может применяться разложение в ряд Фурье — для получения суперпозиции нескольких гармонических задач (также этот метод применим и для нелинейных материалов).
Важной особенностью приложения является отсутствие в моделях магнитов. Вихревые токи наводятся и в пассивных, и в активных (питаемых внешней схемой) проводниках. Возможно объединение с внешней схемой.
Уравнения и условия вычислений следующие. В вычисления включаются уравнения Максвелла для магнитной системы (3), а соотношения материалов, которые принимаются в расчет, имеют вид (4).
Вычисляются B-, H-, E-поля. D-поле не вычисляется. Переменные состояния зависимы от времени (то есть d/dt Ф 0).
Для решения используются две модели:
• векторная, с вычислением переменной состояния магнитного векторного потенциала A;
• скалярная, с вычислением магнитного скалярного потенциала, записываемого как
Фй>(или Фга^
Уравнение Максвелла-Фарадея допускает присутствие электрического скалярного потенциала ф:
E = -(дA/дt)-gradф. (37)
Уравнение для переходных процессов, решаемое в применении Transient Magnetic, записывается как
Tot(v0[vr]rot(A ))+[a](SA/dt+grad<$) = 0.
(38)
Комплексное изображение этого уравнения, в котором d/dt = jro, имеет следующий вид:
rot(v0[vT]rot(A))+[o](j®A+gradq) = 0. (39)
Переменные состояния включают:
• магнитный векторный комплексный потенциал A (An для плоских задач и rAn для
осесимметричных);
• электрический скалярный потенциал.
Flux 3D для скалярной модели Steady State
AC Magnetic использует также различные автоматические формулировки с различными скалярными потенциалами — тех же, что рассматривались ранее в описании магнитостатических приложений.
Описание задачи в Steady State AC Magnetic включает определение частоты (при выборе применения). Для изучения низкочастотных проблем уравнения электрических полей E и D и магнитных полей B и H также могут быть разъединены. Это разделение зависит от частоты, материалов, размера домена исследования. Обычно разделение действительно от 0 до 10 ГГц.
О моделировании скин-эффекта необходимо сказать еще несколько слов.
Если материал характеризуется линейными, гомогенными и изотропными свойствами, среднее значение глубины скин-эффекта можно вычислять согласно формуле (36).
Физические параметры электрического тока или магнитного поля экспоненциально уменьшаются при увеличении глубины скин-эффекта, если область материала плоская. Для того чтобы точно оценить физические параметры в скин-регионе, по крайней мере два слоя элементов (второго порядка) в толщине необходимы. Некоторые проблемы возникают, когда глубина 8 ^ L или 8 ^ W — длины или ширины материала. Но во Flux существует возможность задания областей материала с поверхностным импедансом.
Гипотеза основывается на существовании нулевого вектора электрического поля в проводнике E — за исключением поверхности проводящего региона, где присутствует поверхностная плотность тока jS.
Поверхностный импеданс ZS соотносит тангенциальный компонент магнитного поля с тангенциальным компонентом электрического поля на поверхности проводящего региона с использованием следующего соотношения:
Es = Zs(nHs), c Zs = (1+j)/a8, (40)
где n — единичный вектор, направленный нормально к поверхности, исходящий из проводящего региона.
Рис. 3. Синусоидальный ток и его комплексное изображение
Ограничения по определению поверхностного импеданса задаются двумя условиями: 8 < L/5 и 8 < г/5, где L — толщина, r — радиус кривизны.
Поле должно быть преимущественно тангенциальным (допускаются только незначительные нормальные составляющие). Эти ограничения записываются в виде: 8/L > 10(^0/ц) и ц > 10ц0.
Проводящий твердотельный регион, описываемый условием поверхностного импеданса (Region of solid conductor type described by surface impedance), представляет собой комбинацию:
• неактивного объемного региона,
• контура региона с поверхностным импедансом.
Неактивный регион можно не покрывать сеткой, поскольку переменные состояния на внутренних узлах не вычисляются.
Flux 3D допускает также возможность моделирования тонких проводящих регионов пассивных проводников, в которых плотность тока моделируется постоянной (2D) или гиперболически распределенной по толщине (3D), а также областей с малым поперечным сечением — филиформ.
Источники в применении Steady State AC Magnetic бывают следующих типов:
• электрический ток (синусоидальный);
• наложенное значение напряженности Н0 магнитного поля ^D-применение) однородное или распределенное, синусоидально изменяющееся во времени;
• наложенное значение плотности ф0 магнитного потока (2D-применение), характеризующееся теми же свойствами, что и Но.
Параметры постобработки в Steady State AC Magnetic — скалярные или векторные, также синусоидальные.
Синусоидальный скалярный параметр угловой частоты пульсаций (переменного тока) ю может быть геометрически представлен как вектор, вращающийся с угловой скоростью ю (рис. 3) [8]. В комплексном домене вращающийся вектор для данного отрезка времени является представлением комплексного числа. Значение вращающегося вектора для t = 0 представляет собой комплексное изображение синусоидально изменяющегося скалярного параметра.
Мгновенное значение синусоидального тока равно изображению комплексного числа le’"’:
i(t) = /sin(®t+P),
(41)
где I = Ie№ (42) — это так называемое комплексное изображение тока, то есть i(t) = = Imag(Ie®) (43).
Комплексное мгновенное значение синусоидального тока определяется следующим соотношением:
• в декартовой форме:
i(t) = Icos(rot+P)+jIsin(rot+P); (44)
• в экспоненциальной форме:
-(t) = Iej(rot+e), (45)
где I — модуль комплексного значения i(t); (rot+P) — аргумент или фаза комплексного значения -(t), Р — начальная фаза (для t = 0). Вращающийся вектор, назначенный синусоидальному параметру
A(t) = Âsin(rot+y0), (46)
выражается в форме комплексной функции времени A(t). Эта комплексная функция времени может быть записана:
• в декартовой форме:
A(t) = Âcos(rot+y0)+jÂsin(rot+y0); (47)
• в экспоненциальной форме:
A(t) = Âe1(rot+v0), (48)
где далее она может быть разложена на два фактора:
A(t) = Âei(rot+y0) = Âeyirot, (49)
A = Âe^0 = Âcos(v0)+jÂsin(v0), (50)
фазор (phasor) содержит информацию об амплитуде и начальной фазе синусоидально изменяющегося параметра — это и есть комплексное изображение параметра; e’rot = cosrot+jsinroi (51) — информацию о временной вариации параметра. Синусоидальная функция во временном домене может быть представлена в комплекс-
ном домене с помощью первого фактора — комплексного изображения функции. Фазор (phasor) определяется как вращающийся вектор в момент t = 0 в комплексном домене, назначенный параметру A(t).
Постобработка скалярных результатов осуществляется различными путями:
• мгновенная постобработка — для момента времени rnt. значение параметра A(t); значение комплексного вращающегося вектора (модуль/аргумент (фаза) или действительная/мнимая часть);
• не зависимая во времени постобработка (во Flux 2D): фазор или комплексное изображение (модуль/аргумент (фаза) или действительная/мнимая часть). Векторные параметры раскладываются на
компоненты, каждый из которых записывается как синусоидальный параметр:
Bx(t) = Bxsin{at+vox), (52)
или, в комплексной системе изображения, как вращающийся вектор:
Bx(t) = BXxe*“+4 Bx(t) = Bxxcos(®t+vox)+ jBXx
sin(®t+v0x), (53)
где Bx — модуль комплексного компонента Bx; (rot+vox) — аргумент или фаза Bx В общем случае векторные параметры описывают эллиптические пути периода в 2D-вычислительном домене (x, У). Выражения двух компонентов вектора плотности магнитного потока B в точке вычислительного домена записываются как
BX(t)
BY(t)
- BXxcos(®t+v0X), : BYxcos(®t+v0Y)
(54)
и представляют собой не что иное, как параметрические уравнения эллипса — эллиптического изменения параметра во времени і (рис. 4) [8].
В зависимости от их эллиптической характеристики два компонента плотности В магнитного потока выражаются как следующее:
Bx(t) = V(axcos0)2+(bxsin0)2xcos(®t+vox), tg(vox) = (b/a)tg0,
Ву 4 B(t) /
^2
Рис. 4. Эллиптическое представление двух компонентов вектора плотности магнитного потока в 2D-вычислительном домене
By(t) = V(axsin0)2+(bxcos0)2xcos(at+voy), tg(voy) = -(b/a)ctg0. (55)
Постобработка векторных результатов возможна в следующих вариантах:
• мгновенная постобработка — для данного
момента времени at.
- значение векторного параметра B(t);
- значение комплексного вращающегося вектора B(t) (модуль/аргумент (фаза)) или действительная/мнимая часть);
• не зависимая от времени постобработка:
- электрическое представление (в 2D): эллиптические характеристики (a, b и другие.)
Пульсирующие величины — периодические параметры с ненулевым средним значением — также могут анализироваться с постобработкой.
Мгновенная электрическая мощность как сумма двух компонентов P(t) — абсорбируемая или поставляемая посредством диполя — представляет собой сумму двух компонентов:
• не зависимого во времени:
Pm = UxIxcos9; (56)
• периодического компонента с двойной частотой:
Pp = UxIxcos(2at+9). (57)
Мгновенное значение силы на лопастях контактора — это также сумма двух компонентов (компонента, не зависимого от времени Fm, и периодического):
Fmxcos(2at+9). (58)
Для постобработки скалярного параметра, такого как мощность, главный интерес представляет вычисление постоянного компонента (среднего значения). Условие пульсации малоинтересно. Вычисление мгновенного значения мощности Pm(t) не выполняется ни во Flux 2D, ни во Flux 3D. Еще одна возможность для вычислений — временное интегрирование:
pm=^]mdt. (59)
* о
Для векторных параметров, таких как сила, представляют интерес именно два компонента — среднее значение и параметр пульсаций. Вычисляются как мгновенные значения, например F(t), так и среднее значение — посредством временной интеграции:
-> 1 Т ->
Pm = -jF(i)dt, (60)
1 о
производится постобработка условия пульсации FP: стандартная постобработка, эллиптическое 2D-представление.
Локальные и глобальные параметры, доступные для постобработки/анализа результатов решения задачи в приложении Steady State AC Magnetic Flux, включают представленные в таблице 2.
Таблица 2. Локальные и глобальные параметры постобработки магнитогармонических результатов во Flux 2D/3D
Параметры Единицы измерения Объяснение
Локальные комплексные 2Э-параметры
Векторный потенциал: • для плоской задачи — А2, • для осесимметричной задачи — тА2 Вб/м Вб
Напряженность магнитного поля, Н А/м - = vx rotAz
Плотность магнитного потока, В Тл 1 = [ц)1
Магнитный поток, Ф (фрр — магнитный поток через поверхность, не ограниченную линиями, нормальными к домену исследования, передаваемый посредством точек Р1 и Р2) Вб фР1Р2 = ZPl~-AZP2)P’
p — толщина домена плоской задачи
Плотность тока источника, ]2 А/мм2 -Z = rot (H)
Лапласова сила (локальное значение), Рь Н/м3 - L = - x1
Локальные средние 2D параметры
Объемная плотность активной мощности, dP Вт/м3 dP = (1/2)р-2х-/
Другие локальные 2D-параметры
Проницаемость, ^ Гн/м
Удельное сопротивление, р Ом/м
Глобальные комплексные 2D-параметры
Магнитная сила, Pm (Magnetic Pressure) Н Интегрирование поверхностных сил
Магнитный крутящий момент, Pm (Torque) Н/м (тензорный метод Максвелла)
Ток (источника), Iz (Current) А ¿ = jj2dS S
Магнитный поток в области, Ф (Inductance/Flux) Вб Ф = \jAzdl í
Глобальные действительные 2D-параметры
Индуктивность, L (Inductance/Inductance) Гн L = [пЩії
Глобальные средние действительные 2D-параметры
Магнитная сила, Рт Н Fm = dWJdl
Магнитный крутящий момент, Г т Н/м Г m = dWrn/дв
Активная мощность, Р Вт P = j(dP)dV V
Реактивная мощность, 0 Вар Q = j(1/2)(ra/|i)IBPdV V
Лапласова сила, Рь = |(- хВ )dV: Н
• средний действительный вектор; • пульсирующий компонент.
Локальные комплексные SD-параметры
Скалярные магнитные потенциалы, ф{0{, фгец Фга$То А/м
Электрический векторный потенциал, — А/м Непрямой доступ
Напряженность магнитного поля, Н А/м H = -grad(-tot)'H = -grad(-.r<dH)+Hj. H = -gradMwiiH-o
Плотность магнитного потока, В Тл 1 = №
Электрическое поле, Е (в твердотельном проводнике) В/м E = ((1/a) rot-)
Магнитное поле источников, не покрываемых сеткой, Ц А/м Аналитическое или полуаналитическое вычисление с использованием закона Био-Савара-Лапласа
Плотность источника тока в сетчатых обмотках, А/м2 Верификация данных
Плотность вихревых токов, - в твердотельном проводнике) А/м2 ; = rot (f)
Другие локальные SD-параметры
Проницаемость, ^ Гн/м
Проводимость, а См
Удельное сопротивление, р Ом-м
Относительная проницаемость,
Средние комплексные локальные SD-параметры
Комплексное значение плотности мощности, dP+jdQ: • действительное значение, dP; • мнимое значение, jdQ Вт/м3 Вар/м3
Средние действительные локальные SD-параметры
Плотность потерь вследствие эффекта Джоуля (объемная плотность активной мощности, ¿Р) Вт/м3 dP = (1/2)Ех/
Объемная плотность реактивной мощности, dQ Вар/м3 dWm = jU dB
Объемная плотность магнитной энергии, Дж/м3
Локальные параметры 2ф пульсации
Плотность поверхностной силы между материалами, ¿Р5: • среднее значение; • пульсирующий компонент Н/м2 Вычисление на границах сред с различной магнитной проницаемостью
Объемная плотность лапласовой силы, ¿Рь, в проводниках stranded-типа/сетчатых обмотках: • среднее значение; • пульсирующий компонент Н/м3 dF L = ÍsxB
Глобальные SD-параметры (средние значения)
Магнитная сила, Рт Н Р т = д^тЮХ
Магнитный крутящий момент, Рт Н/м Гт = дWmlдв
Магнитная энергия, Wm Дж ^т = V
Комплексное значение магнитного потока через обмотку-проводник (stranded-типа), Фк Вб Фк = д^т!д!к = Вхк0к^П,где коК — поле, созданное обмоткой к в вакууме с током Iк в1А, П — домен исследования
Переходные процессы Transient Magnetic
Модуль Transient Magnetic позволяет пользователю производить магнитные исследования устройств с принятием в расчет переходных магнитных процессов, если это необходимо. Применение Transient Magnetic делает возможным исследование явлений, обусловленных переменным во времени магнитным полем. Переменное магнитное поле создается переменными токами различной формы и постоянными магнитами. Приложение позволяет исследовать наведенные токи, скин-эффект и эффект близости в проводящих областях.
Transient Magnetic допускает схемное и кинематическое объединение.
Объединение с внешней схемой позволяет пользователю принимать в расчет наведенные токи в обмотках типа stranded (многожильных) — постоянные наведенные токи в проводниках без скин-эффекта. Кроме того, применение позволяет исследовать угловое движение и вихревые токи в пассивных твердотельных проводниках, в которых токи наводятся, но не проводятся.
В переходной магнитной модели действует гипотеза квазистатического состояния, согласно которой вариациями во времени электрического потока D можно пренебречь.
Принципиальные результаты, достижимые с этим применением, — те же, что и для магнитостатического применения. А также наведенные токи и рассеиваемая вследствие эффекта Джоуля мощность.
Уравнения, используемые для решения переходных проблем, представляют собой уравнения Максвелла для магнитной системы (3) и конституционные соотношения среды (4) — те же, что и для магнитогармонических процессов.
Условия вычислений следующие:
• переменные состояния зависимы от времени d/dt Ф 0;
• вычисления производятся для B-, H-, E-полей, D-поле не вычисляется.
Для того чтобы решать эти уравнения, также используются две модели:
• векторная (с магнитным векторным потенциалом A);
• скалярная (с магнитным скалярным потенциалом фш или фг^).
Уравнение, решаемое в векторной модели применения Transient Magnetic методом конечных элементов, записывается в виде (38). Переменные состояния включают:
• магнитный векторный потенциал A;
• электрический скалярный потенциал ф. Применение Transient Magnetic Flux 3D использует автоматические формулировки с различными скалярными потенциалами (в скалярной модели).
Описание задачи во многом сходно с магнитостатическим и магнитогармоническим применением. Определяются границы; гра-
ничные условия либо назначаются автоматически с использованием (в большинстве случаев) метода бесконечного бокса (в бесконечности магнитный потенциал нулевой), либо — без использования метода бесконечного бокса — задаются вручную с указанием нематериальных областей. В отсутствие контурных областей применяются стандартные граничные условия. Магнитное поле является тогда тангенциальным (электрическое — нормальным) к границе.
Если домену исследования назначаются симметрии, пользователь должен указать тип симметрии. Нормальное магнитное и тангенциальное электрическое поля или тангенциальное магнитное и нормальное электрическое поля могут быть наложены на плоскости симметрии, определенные в геометрическом контексте. Если домену исследования назначаются периодичности, пользователь должен указать тип периодичностей. Циклические или антициклические условия также могут накладываться на плоскости периодичности, определенные в геометрическом контексте.
Объемы, поверхности и линии материальных регионов допускают в Transient Magnetic моделирование следующих материальных сред:
• воздуха или вакуума (цг = 1);
• магнитомягких и магнитотвердых непроводящих материалов (цг > 1);
• обмоток, проводящих ток, немагнитных
(Ил = 1);
• твердых проводников (активных или пассивных, магнитных (цг > 1) или проводящих (р), в которых скин-эффект очевиден), а также тонких областей и филиформ.
Различают проводники обмоток (с токонесущим проводником), где скин-эффект пре-небрежим вследствие того, что диаметр провода является малым в сравнении с глубиной скин-эффекта, и плотность тока считается однородной, и проводники типа solid conductor, в которых скин-эффект заметен, и плотность тока неоднородна. Кроме того, различают пассивные проводники (питаемые током вследствие закона Фарадея) и активные (непосредственно питаемые током от источника). Активные проводники могут быть открытой (U, E) или замкнутой (не существующей в действительности, но моделируемой во Flux) формы.
Приложение Transient Magnetic отличается от Magneto Static тем, что в нем присутствуют переменные во времени источники тока и анализируются вихревые токи. В магнитостатическом применении электрические токи в проводниках являются непрерывными токами проводимости, распределение токов в проводящей области является однородным — без пространственных вариаций. В Transient Magnetic электрические токи могут быть синусоидальными, трапецеидальными, наведенными. Распределение тока может быть однородным (пренебрежимый скин-эффект, эффект близости отсутствует)
либо неоднородным (скин-эффект/эффект близости).
Различают следующие типы источников поля:
• постоянные магниты;
• токи (с однородным распределением плотности тока) — постоянные или временнозависимые;
• наложенное электрическое поле (3D).
Привязка электрической схемы осуществляется заданием связей регионов Coil conductor type или Solid conductor type конечноэлементного домена и компонентов Stranded conductor type или Solid conductor type ассоциированной электрической схемы.
Наложенное на домен исследования дополнительное поле H0 может быть однородным или пространственно-распределенным, изменяться во времени. Магнитный поток Ф0, наложенный на линейный регион, может также быть однородным или пространственно-распределенным, но не может изменяться во времени.
В случае моделирования проводников с отверстиями (активных или пассивных), например кольцевых, возникает проблема соединения различных потенциалов (рис. 2). Для решения этой проблемы создается математический срез посредством Entity cut of electric loop.
Переходное приложение представляет собой физическое применение, характеризующееся временно-зависимыми свойствами: переменными источниками, переменными характеристиками материалов, дифференциальным temporal-уравнением первого порядка. Численный метод решения представляет собой получение временных разрешений на каждом шаге. Решения независимы друг от друга, каждое из них связано только с предыдущим temporal-уравнением. В общем случае начальные условия предполагают, что переменные состояния сбрасываются к нулю. Во Flux 2D также существует опция инициализации статического вычисления: initialization by static computation.
Анализ результатов решения в переходном приложении сходен с анализом в магнитостатическом и магнитогармоническом приложениях.
В применении Transient Magnetic параметры постобработки являются мгновенными (вычисляются для мгновенного t, соответствующего шагу времени, или представляют собой временную вариацию параметра).
Локальные и глобальные параметры включают представленные в таблицах 1, 3.
Принципы вычисления потерь
Потери мощности в электромеханическом устройстве бывают трех видов:
• магнитные потери в магнитных схемах;
• потери вследствие эффекта Джоуля;
• механические потери (за счет трения и вентиляции).
Таблица 3. Локальные и глобальные параметры постобработки результатов решения в переходном применении во Flux 2D/3D (в дополнение к представленным в табл. 1)
Параметры Единицы измерения Объяснение
Локальные 2Э-параметры
Плотность вихревых токов, ] А/мм2 Компонент Jz
Удельное сопротивление, р Ом-м
Объемная плотность активной мощности, dР Вт/м3 dP = px_f
Глобальные 2Э-параметры
Ток, I (Current) — источника coil и проводника типа Solid Conductor А I = \jzdS S
Активная мощность, P Вт P = ¡(dP)dV V
Локальные ЗЭ-параметры
Электрический векторный потенциал, Т А/м Непрямая доступность
Плотность вихревых токов, j А/м2 j = rot (T )
Напряженность электрического поля, Е В/м E = (1/a)j
Проводимость, а См
Удельное сопротивление, р Ом-м
Объемная плотность потерь вследствие эффекта Джоуля, dР Вт/м3 dP = Exj, может также записываться как dP = pxj2 = axE2
Глобальные ЗЭ-параметры
Рассеянная мощность или потери вследствие эффекта Джоуля в объемной области, Р Вт P = j(dP)dV V
Потери мощности в магнитных материалах относятся к явлениям, ассоциируемым с временными вариациями магнитного поля.
Они классически подразделяются на гис-терезисные потери — микроскопического происхождения, и потери вследствие токов Фуко — макроскопического.
Гистерезисные потери генерируются посредством временной вариации магнитной микроструктуры вещества, главным образом вследствие движения стен магнитных доменов.
Когда магнитное поле характеризуется циклической временной вариацией, магнитные материалы характеризуются геометрическим местом точек (В, Н), представляющим собой замкнутую кривую. Для сниженных значений частоты ^ < 1 Гц), замкнутая кривая В(Н), называемая статическим гистерезисным циклом, не зависит от частоты. Соответствующие магнитные гистерезисные потери пропорциональны области статического гисте-резисного цикла и частоте (квадрату частоты возбуждающего магнитного поля).
Процесс циклического намагничивания определяет трансформацию части электромагнитной энергии в тепловую (магнитотепловые потери). Циклическая вариация магнитного поля во времени определяет токи Фуко в магнитных материалах с электропроводящими свойствами. Потери, связанные с токами Фуко, вызываются эффектом Джоуля вследствие наведенных токов.
При увеличении частоты область замкнутой кривой В(Н), называемая динамическим гистерезисным циклом, увеличивается (рис. 5).
В этом случае магнитные потери вследствие циклического намагничивания больше гистерезисных потерь. Потери, соответствующие разнице между областью динамического и статического гистерезисного цикла, обычно называются дополнительными магнитными потерями.
Для магнитных сталей этими потерями можно пренебрегать на частотах до 50 Гц, но на частотах килогерцового диапазона они становятся непренебрежимыми.
Вычисление магнитных потерь и моделирование магнитомягких материалов основывается на следующих принципах:
• гистерезис принимается в расчет на уровне магнитного закона B(H);
• магнитные потери вычисляются непосредственно;
либо:
• гистерезисом на уровне B(H) пренебрегают;
• вычисление магнитных потерь выполняется апостериори посредством теоретических и экспериментальных формул. Моделирование гистерезиса увеличивает
время вычислений и требует большего объема памяти. Поэтому во Flux применяется второй подход.
Flux обеспечивает пользователю две модели вычислений магнитных потерь. В обоих случаях вычисление производится апостериори. Гистерезис не моделируется на уровне процесса решения, и В(Н)-характеристика представляет собой однозначное соотношение.
Вычисление мощности, соответствующей магнитным потерям, производится в соответствии со следующим выражением для магнитных потерь в связи с гистерезисом:
B,
/ / S / / \ --fl * / : —f2>1
/ / / 1 ^
/ / / / У н
yfe-
Рис. 5. Магнитные потери и гистерезисный цикл
Ph = khxfxB MAX, (61)
где п — коэффициент Штейнметца (1,6-2); f— частота; BMAX — пиковое значение плотности потока магнитного поля.
Потери, сгенерированные токами Фуко, вычисляются как
Pec = kecf 2BMAX- (62)
Магнитные потери во Flux вычисляются средствами формул Бертотти. Общие магнитные потери могут быть разделены на следующие три категории:
• потери посредством гистерезиса (P1), пропорциональные частоте f, которые являются наиболее значительными компонентами магнитных потерь на низкой частоте;
• потери вследствие токов Фуко (P2), пропорциональные квадрату частоты;
• дополнительные потери (P3), пропорциональные f3/2.
Теория Бертотти дает выражение магнитных потерь как функции частоты и пикового значения плотности магнитного потока:
dP = c1BMAx/+ c2(BMAXf)2+ c3(BMAXf)Ш,(63) P1 = c1BMAXf,
P2 = c2(BMAXf )2,
P3 = c3(BMAXf )3/2,
где c1, c2, c3 — соответствующие коэффициенты, BMAX — максимальное значение достижимой индукции.
Коэффициенты С2 и c3 выражаются посредством следующих соотношений:
С2 = (n2xaxd2)/6, (64)
с3 = 4<jxGxV0xS, (65)
где d — толщина листа (пластины); G — безразмерная константа; S — поперечное сечение листа (пластины); V0 — постоянное поле, которое зависит от разницы коэрцитивной силы между двумя магнитными объектами.
В применении Steady State AC Magnetic используется комплексное представление. Объемная плотность средней мощности dPmoy записывается как
dPmoy = khfBMAX+((n2tf d 2)/6)(BMAXf )2+
+ ke(BMAXf )32x8,67, (66)
где kh, ke — коэффициенты потерь.
В вычислительном фрейме Flux это выражение записывается как
dPmoy = [ khfBMAX+((n2tfd 2)/6)(BMAXf )2+
+ke(BMAXf )3/2x8,67] kf, (67)
где kf — коэффициент филлинга — заполнения (близкий к единице), или электрической изоляции пластин магнитного сердечника.
Средняя мощность, рассеиваемая в объемном регионе, записывается как
Pmoy \\\ reg dPmoydV'
reg moy
(68)
Важно отметить, что BMAX представляет собой пиковую индукцию, и в симуляциях, где вычисленное значение не соответствует пику (ротор и статор моделируются неподвижными, а удельное сопротивление моделируется переменным), вычисления магнитных потерь имеют ограничения. Так же и в том случае, когда нелинейная аппроксимация B(H) и насыщение вводятся средствами эквивалентной модели намагничивания, значения B могут изменяться.
В применении Transient Magnetic изучается переменный режим поля. В зависимости от магнитных потерь объемная плотность мощности будет записываться как
dP(t) = khfB[MAX+a(d2/12)((dB/dt)(t))2+ +ke((dB/dt)(t))3/2, (69)
в вычислительном фрейме Flux записываемый с коэффициентом заполнения:
dP(t) = [khfBjMAX+a(d 2/12)((dB/dt)(t))2+ +ke((dB/dt)(t))3/2]kf. (70)
Объемная плотность средней мощности за период dPmoy записывается как
о
которая далее становится во Flux
dPmoy = К /®МАХ^/ +
if ‘ г J
12
dB
di
(*)
dB
di
4.3/2
/ .
Средняя мощность, рассеиваемая в объемном регионе, записывается в соответствии с (68).
Вычисление магнитных потерь во Flux также производится с помощью модели LS (Loss Surface). Эта модель представляет собой метод оценки магнитных потерь апостериори и основана на модели динамического гистерезиса, ассоциируемого с конечно-элементной симуляцией. В соответствии с LS-моделью магнитное поведение материала определяется экспериментально — посредством характеристической поверхности H(B, dB/dt).
Таким образом, от сигнала B(t) конечной формы и частоты можно через поверхность H(B, dB/dt) прийти к полю H(t), что позволяет реконструировать динамический цикл гистерезиса.
Аналитическая модель позволяет восстановление сигнала H(t) напряженности магнитного поля с начальными данными сигнала B(t) плотности магнитного потока:
H(B, dB/dt) =
dynamic(B dB/dt).
(73)
(71)
kfdt. (72)
Кривая B(H), которая может быть получена, позволяет осуществлять точные вычисления железных потерь. Можно также отметить, что LS-модель выполняется во Flux не для всех материалов. ■
Продолжение следует. Литература
1. Сысоева С. Развитие концепции математического и расчетного моделирования автомобильных датчиков скорости/положения // Компоненты и технологии. 2007. № 12.
2. Сысоева С. Геометрическая и физическая параметризация в конечно-элементном пакете Flux v.10.1.2 Cedrat как подход к оптимизации датчиков скорости/положения // Компоненты и технологии. 2008. № 4.
3. Сысоева С. Геометрическая и физическая параметризация в конечно-элементном пакете Flux v.10.1.2 Cedrat как подход к оптимизации датчиков скорости/положения // Компоненты и технологии. 2008. № 5.
4. Сысоева С. Практический курс моделирования магнитоуправляемых датчиков положения и скорости на базе конечно-элементного программного пакета Flux. Часть 1. Введение в рабочую среду Flux // Современная электроника. 2008. № 5.
5. Сысоева С. Практический курс моделирования магнитоуправляемых датчиков положения и скорости на базе конечно-элементного программного пакета Flux. Часть 2. Магнитостатическое моделирование компонентов рабочей зоны датчика // Современная электроника. 2008. № 7.
6. Сысоева С. Практический курс моделирования магнитоуправляемых датчиков положения и скорости на базе конечно-элементного программного пакета Flux. Часть 2. Магнитостатическое моделирование компонентов рабочей зоны датчика (продолжение) // Современная электроника. 2008. № 8.
7. Сысоева С. Практический курс моделирования магнитоуправляемых датчиков положения и скорости на базе конечно-элементного программного пакета Flux. Часть 2. Магнитостатическое моделирование компонентов рабочей зоны датчика (продолжение) // Современная электроника. 2008. № 9.
8. Техническая документация фирмы Cedrat. Flux 10 User's Guide. Vol.3. Physical applications: Magnetic, Electric, Thermal. Aug 2007.
9. Техническая документация фирмы Cedrat. Flux 9.1 3D application. User's Guide. Vol.5. Physical applications (complements for advanced users). March 2005.
