Теоретические аспекты конечно-элементного моделирования состояний и функциональности магнитоуправляемых датчиков. Часть 1

№ 1’2009
PDF версия
В статье рассматриваются аспекты математического моделирования магнитостатических, магнитодинамических и магнитотепловых (мультифизических) процессов, ответственных за функциональность магнитоуправляемых датчиков скорости и положения, и параллельно раскрываются моделирующие возможности и особенности конечно-элементного CAD/CAE пакета FLUX Cedrat.

Все статьи цикла:

Введение

Магнитоуправляемые датчики скорости/положения представляют собой широко распространенный и перспективный класс сенсорных устройств в составе автомобильных, промышленных, аэрокосмических и других систем. Эти устройства отличаются бесконтактностью, готовым типом выхода, часто цифровым исполнением, достаточно высокой чувствительностью, стабильностью, надежностью и точностью. Чтобы разработать и гарантировать работу такого устройства в полном диапазоне рабочих условий, включая в анализ материалы и геометрию магнитной системы, допуски, рабочие температуры и другие параметры, стоит воспользоваться компьютерным конечно-элементным моделированием, которое способно дать хорошие результаты. Пример такого специализированного программного обеспечения — Flux Cedrat [1–9] (рис. 1).

Рис. 1. Типичная задача, решаемая методом конечных элементов, -  моделирование датчика ферромагнитного зубчатого ротора: а) постановка задачи: 1 - ферромагнитный зубчатый ротор; 2 - основные компоненты статора;  3 - обратно-смещающий магнит; 4 - ИС датчика (дифференциальная ИС Холла TLE4921 Infineon); б–г) решение в рамках реализации конечно-элементного проекта во Flux v.10.1.2: б) результат генерирования сетки для линий и контуров; в) результат генерирования сетки для элементов второго порядка Generate Second Order Elements; г) постобработка: изображение векторных результатов решения образцовой модели
Рис. 1. Типичная задача, решаемая методом конечных элементов, —
моделирование датчика ферромагнитного зубчатого ротора:
а) постановка задачи: 1 — ферромагнитный зубчатый ротор; 2 — основные компоненты статора;
3 — обратно-смещающий магнит; 4 — ИС датчика (дифференциальная ИС Холла TLE4921 Infineon);
б–г) решение в рамках реализации конечно-элементного проекта во Flux v.10.1.2:
б) результат генерирования сетки для линий и контуров;
в) результат генерирования сетки для элементов второго порядка Generate Second Order Elements;
г) постобработка: изображение векторных результатов решения образцовой модели

Существуют три основных модели и/или состояния, которые характеризуют функциональность магнитоуправляемых датчиков, соответствующих приложениям (модулям) Flux:

  1. Магнитостатическое (магнитное поле постоянное, источники поля — стационарные токи проводимости или постоянные стационарные магниты, в системе отсутствуют механическое движение источника или модулятора поля и наведенные токи), изучается в Flux MagnetoStatic.
  2. Магнитогармоническое (магнитное поле — установившееся переменное синусоидальное, источники поля — переменные токи, подвижные постоянные магниты и/или неподвижные постоянные магниты, стационарное поле которых синусоидально модулируется подвижным ферромагнитным/проводящим модулятором, в расчет принимаются наведенные токи, скин-эффект и эффект близости), соответствует приложению Steady state AC Magnetic Flux, но данный модуль не рассматривает магниты в качестве источников.
  3. Переходные, соответствующие переходу системы из одного установившегося магнитостатического или магнитодинамического состояния в другое и характеризующиеся переменным несинусоидальным состоянием источников и полей. Во Flux переходные процессы изучаются посредством специализированного модуля Transient Magnetic, допускающего задание в качестве источников и переменных (несинусоидальных) токов, и магнитов, анализ наведенных токов, скин-эффекта и эффекта близости.

Магнитотепловое моделирование, также доступное во Flux, основано на существовании взаимодействия между электромагнитными и тепловыми аспектами процессов нагрева — вследствие вихревых токов.

Теоретические аспекты моделирования основных функциональных процессов, протекающих в датчиках (магнитных и магнитнотепловых), достаточно подробно рассматриваются в рамках данной статьи. Но нельзя не отметить, что возможности моделирования, предлагаемые Flux, шире. Так, возможность объединения с электрической схемой «circuit coupling» доступна пользователю в приложениях Transient Magnetic и Steady State AC Magnetic. Кинематический модуль Flux также допускает исследование подвижной части устройства вследствие механических сил (действия пружин, трения, гравитации).

Магнитостатическое
моделирование

Магнитостатическое моделирование основывается на фундаментальных законах электромагнетизма — уравнениях Максвелла и конституционных соотношениях среды. В самом общем виде уравнения Максвелла (в дифференциальной форме) записываются как:

rotvec{E} = − ∂vec{B}/∂t,
rotvec{H} = vec{j}+(∂vec{D}/∂t), (1)
divvec{B} = 0,
divvec{D} = q,

где vec{B} — плотность магнитного потока, vec{j} — плотность тока, vec{H} — напряженность магнитного поля, vec{D} и vec{E} — электрическое смещение и напряженность электрического поля соответственно, q — плотность статического заряда.

Соотношения для материалов записываются следующим образом:

vec{j} = σ vec{E},
vec{B} = μ vec{H}, (2)
vec{D} = ε vec{E},

где σ — проводимость; μ — магнитная проницаемость; ε — диэлектрическая проницаемость.

В случае низкочастотных квазистатических электромагнитных полей выражения для электрических полей vec{D} и vec{E}, vec{B} и vec{H} могут быть разъединены и рассмотрены отдельно.

Для чисто магнитной системы принимается гипотеза квазистатического состояния, вариацией плотности электрического потока vec{D} можно пренебречь. Уравнения Максвелла тогда запишутся как

rotvec{E} = – ∂vec{B}/∂t,
rotvec{H} = vec{j}, (3)
divvec{B} = 0,

а соотношения материалов, которые принимаются в расчет, так:

vec{j} = σ vec{E},
vec{B} = μ vec{H}. (4)

Условия вычисления в магнитостатическом исследовании

Переменные состояния не зависят от времени, то есть d/dt = 0, вычисляются только B— и H-поля, D— и E-поля не вычисляются.

Наведенные токи не учитываются, и уравнения Максвелла для задач магнитостатики тогда записываются как

rotvec{H} = vec{j},
divvec{B} = 0, (5)

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

vec{B} = μ vec{H}. (6)

Это принципиальное соотношение для магнитных материалов может быть помещено в расчетную модель в виде B(H)- или H(B)-соотношений:

vec{B} = μrμ0vec{H},
vec{B} = μrμ0vec{H} + vec{B}r, (7)

где μr — относительная магнитная проницаемость среды; μ0 — абсолютная магнитная проницаемость вакуума; vec{B}r — остаточный магнитный поток постоянных магнитов, или

vec{H} = vrv0vec{B},
vec{H} = vrv0vec{B}vec{H}c, (8)

где vr — относительное удельное магнитное сопротивление среды; v0 — абсолютное удельное магнитное сопротивление вакуума; vec{H}c — коэрцитивная сила постоянных магнитов.

Закон Гаусса для магнетизма divvec{B} = 0 также допускает введение магнитного векторного потенциала vec{A}:

vec{B} = rotvec{A}. (9)

Выражение для закона Ампера rotvec{H} = vec{j} позволяет дать выражение для магнитного скалярного потенциала:

vec{H} = − vec{grad}phi + vec{T}, (10)

где phi — магнитный скалярный потенциал; vec{T} — вращательный или нулевой компонент поля источника или электрический векторный потенциал в А/м.

Некоторые геометрические аспекты конструирования магнитных систем

В соответствии с законом Ампера в электромагнитной схеме магнитное поле электрически возбуждаемой магнитной системы (рис. 2) определяется соотношением

ointvec{H}dl = I × n, (11)

где l — длина магнитного пути; n — число витков возбуждающей обмотки.

Рис. 2. Замкнутые электромагнитные схемы: а) электромагнитно возбуждаемая;  б) магнитная система  (возбуждаемая постоянным магнитом)
Рис. 2. Замкнутые электромагнитные схемы:
а) электромагнитно возбуждаемая;
б) магнитная система
(возбуждаемая постоянным магнитом)

Для магнитной системы, возбуждаемой постоянным магнитом:

ointvec{H}× dl = 0 (12)

или

vec{H}B lB+ vec{H}Ф lФ+ vec{H}M lM= 0, (13)

где vec{H}B — напряженность магнитного поля в воздухе; lВ — длина пути линий поля в воздухе; vec{H}Ф — магнитное поле в ферромагнитном материале; lФ — длина пути линий поля в ферромагнитном материале; vec{H}M — магнитное поле в магнитном материале; lМ — длина пути линий поля в магнитном материале.

При этом напряженности поля в ферромагнитной и воздушной частях замкнутой схемы будут характеризоваться тем же самым направлением, в магнитной — противоположным.

Согласно этой записи, увеличения магнитного поля в воздушном зазоре можно достичь, уменьшая воздушный зазор и увеличивая длину магнитно-ферромагнитного пути.

В пределах пути lМпостоянного магнита напряженность и индукция магнитного поля связаны соотношением:

vec{B}M = μ0vec{H}+ vec{M}M, (14)

где vec{M}M — намагничивание постоянного магнита, направление которого определяет направление поля в контуре из ферромагнитной и воздушной частей схемы.

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

vec{B}M =
vec{B}B =
vec{B}Ф,

переписав (13) как

(BBlB0)+(BФlФrμ0)+
+((BMMM)/μ0) lM= 0

и, далее, как

BBlB+(BBlФr)+(BB– MM) lM= 0
и lB+(lФr)+ lM–(MMlM/BB) = 0,

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

BB= MM/(1+(lB/lM)+(lФrlM)). (15)

Из этой формулы явствует, что для увеличения поля в воздухе в схеме с постоянным магнитом должно выполняться соотношение:

lМ>> lВ. (16)

Соотношение lМ> lФ в датчиках можно на практике не принимать в расчет, так как lФ/(μr× lМ) — всегда малое значение из-за того, что μr, магнитная проницаемость стали, обычно >1000. Тем не менее, ферромагнитную часть все же не следует делать превышающей по длине магнитную.

Однако длина ферромагнитного пути в сравнении с длиной воздушного пути — и в том случае, если магнитный участок отсутствует, и во всех остальных — должна быть значительно увеличенной, то есть

lФ>> lВ, lФ+ lМ>> lВ.

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

vec{H}BlB+ vec{H}ФlФ= nI (17)

или сведенным к виду

(BBlB0)+(BФlФrμ0) = nI, (18)

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

BB= nI μ0/(lB+(lФr)). (19)

В то же время, если удается выполнять условие (16), достаточно большие зазоры — порядка нескольких мм — могут не оказывать влияния на величину поля в воздухе (хотя она и уменьшается по сравнению с реальным значением BФ) и не влиять на точность измерений.

Еще один важный практический аспект для геометрического конструирования датчиков — это возможность расширения воздушного зазора за счет увеличения поперечного сечения магнита или электромагнитно возбуждаемой ферромагнитной части.

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

BMSM= BBSB, (20)

это далее дает следующее выражение для индукции в воздушном зазоре (при идеальных условиях):

BBlB+(BBlФr)+(BB(SB/SM)− MM) lM= 0,
BB= MM/((SB/SM)+(lB/lM)+(lФrlM)), (21)

и при соблюдении условия (16) и большом значении μrвыражение (21) аппроксимируется к виду:

BB= MM/(SB/SM) = (MMSM)/SB, (22)

то есть увеличение индукции определяется также фактором SM/SB.

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

Магнитостатическое моделирование во Flux

Во Flux существует приложение Magneto Static, позволяющее получать следующие принципиальные результаты:

  • распределение напряженности магнитного поля H, плотности магнитного потока B, плотности магнитной энергии dWm;
  • расчет магнитной энергии в регионе Wm во всем домене исследования, силы Fm и магнитного крутящего момента Tm;
  • магнитный поток, пересекающий поверхность или навивку.

Основные этапы магнитостатического исследования включают:

  • магнитостатическое исследование образцовой модели, решение опорного сценария;
  • физическое параметризованное исследование;
  • геометрическое параметризованное исследование.

Для решения используются две модели:

  • векторная модель с вычислением переменной состояния магнитного векторного потенциала vec{A}, введение которого допускает закон Гаусса для магнетизма;
  • скалярная модель с вычислением магнитного скалярного потенциала, записываемого как phitot или phired.

2D-модуль Flux обеспечивает решение только векторной модели, 3D-модуль — и векторной, и скалярной. По умолчанию 3Dрешатель использует векторную модель для 2D-приложений и скалярную для 3D-приложений. Но в дальнейшем на деталях работы с Flux в рамках данной статьи внимание не акцентируется, отчасти в связи с предшествующей публикацией ряда практических руководств именно по магнитостатическому моделированию магнитоуправляемых датчиков во Flux [2–7].

Для векторной модели решаемое во Flux уравнение записывается как

rot (v0[vr] rot ( vec{A}) − vec{H}c ) = vec{j}s, (23)

где [vr] — тензор относительного удельного сопротивления среды; vec{j}s — поверхностная плотность источника тока (в А/м2).

Для скалярной модели vec{T} выбирается различным образом, и этот выбор приводит к различным скалярным потенциалам и применению различных магнитных формулировок во Flux 3D. Каждая область задачи связывается со своей формулировкой, но в общем виде форма решаемого в магнитостатическом приложении уравнения записывается как

div ([μr0(− vec{grad}( phi )+ vec{T})+ vec{B}r ) = 0, (24)

где [μr] — тензор относительной магнитной проницаемости среды; phi — магнитный скалярный потенциал ( phitot или phired); vec{T} — векторная характеристика источника поля; vec{B}r — остаточный магнитный поток магнитов.

Существуют три основных ситуации, относящиеся к источникам:

  • источников тока нет;
  • источники не покрываются сеткой;
  • источники покрываются сеткой.

Соответствующие им уравнения представлены далее.

Тотальный скалярный потенциал phitot используется при отсутствии токовых источников. Плотность тока j равна нулю. Напряженность поля выражается как

vec{H}= − vec{grad}( phitot), (25)

и решаемое уравнение представляет собой

div ([μr0(− vec{grad}( phitot)− vec{B}r ) = 0, (26)

phitot является переменной состояния, относительно которой решается уравнение (26).

Редукционный магнитный скалярный потенциал, соответствующий H или phiredHj, используется, когда источники — не покрываемого сеткой типа. Во Flux это единичные объекты (entity), которые могут быть круговой, прямоугольной или более сложной формы. Это случай, когда
vec{T} = vec{H}j, и Hj вычисляется аналитически с использованием закона Био-Савара-Лапласа. Закон позволяет рассчитать магнитное поле в точке P, созданное проводником с током I, без присутствия магнитных материалов как:

  • элементарное
    dvec{H}(P) = (1/4 π) Ivec{dl}×( vec{r}/r3); (27)
  • тотальное
    vec{H}(P) = (1/4 π) ointIvec{dl}×( vec{r}/r3). (28)

В приведенных выражениях vec{r} — это вектор между элементом vec{dl} проводника и точкой P.

Для поверхностного или объемного вычисления тока выражения (27) и (28) должны быть переписаны: Ivec{dl} заменяется на vec{j}SdS и vec{j}VdV соответственно.

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

vec{H}= − vec{grad}( phiredHj ) + vec{H}j, (29)

а решенное уравнение представляет собой

div ([μr0(− vec{grad}( phiredHj ) + vec{H}j ) − vec{B}r ) = 0, (30)

где phiredHj — переменная состояния.

Кроме того, магнитное поле может быть наложенным (например, магнитное поле Земли).

Редукционный магнитный скалярный потенциал, соответствующий T0 или phiredT0, используется, когда источники — сетчатого типа, для которых

rotvec{T}0 = vec{j}s, (31)

но может использоваться и с несетчатыми источниками, для которых

rotvec{T}0 = vec{j}s + rotvec{H}j. (32)

Магнитное поле выражается как

vec{H} = − vec{grad}( phiredT0 ) + vec{T}0, (33)

а решенное уравнение представляет собой

div ([μr0(− vec{grad}( phiredT0 ) + vec{T}0 ) − vec{B}r ) = 0, (34)

где phiredT0 и vec{T}0 — переменные состояния.

Описание модели во Flux включает:

  • определение границ и граничных условий;
  • описание материальных сред (воздушных, магнитных/ферромагнитных материалов, источников тока);
  • описание источников поля (постоянных магнитов, DC тока, наложенного поля).

Граничные условия задаются в большинстве случаев методом бесконечного бокса: программой автоматически назначается условие граничного нулевого магнитного потенциала. Альтернатива — задание нематериальных областей пользователем и определение направления поля на границах как нормального или тангенциального (по умолчанию поле назначается тангенциальным к границе). Если в модели присутствуют симметрии и периодичности, пользователь также определяет направление поля как тангенциальное или нормальное по отношению к плоскостям симметрии, циклические или антициклические условия по отношению к плоскостям периодичности, определенным в геометрическом контексте.

Объемы, поверхности, линии материальных областей позволяют моделировать во Flux материальные среды — воздушные (с μr= 1), магнитные (магнитомягкие и магнитотвердые — c μr) непроводящие, проводящие немагнитные (μr= 1). Для моделирования доступны тонкие области и филиформы (области с малым поперечным сечением: когда два размера значительно меньше третьего, магнитное поле является тангенциальным к линии моделирования филиформы).

Постоянные магниты как источники поля описываются посредством назначения областей магнитных материалов. При определении материала вводятся характеристики магнита, остаточная индукция, указывается направление намагниченности. Электрические токи как источники в магнитостатическом применении считаются непрерывными, независимыми от времени, с однородно распределенной плотностью тока в поперечном сечении проводника. Электрические источники описываются либо посредством областей материала типа coil conductor, либо посредством entities coils. Дополнительной возможностью является задание наложенного поля как создаваемого магнитным диполем или известным магнитным полем — однородным или пространственно распределенным.

Практическое магнитостатическое моделирование во Flux на примере датчика ферромагнитного зубчатого ротора было показано в публикациях [2–7], к которым (а также к технической документации Cedrat) рекомендуем обращаться за выяснением деталей задания свойств материалов, указания направления осуществления постобработки результатов и т. п.

В случае замкнутой магнитной системы с высокой магнитной проницаемостью возникает также задача соединения различных магнитных областей, если интеграл в выражении (11) не равен нулю. Решением является искусственное моделирование зазора посредством единичного объекта Magnetic Cut, но в том случае, если замкнутая система рассечена плоскостью симметрии и периодичности, введение зазора не требуется.

Анализ результатов предполагает постобработку локальных и глобальных параметров, различных для 2D- и 3D-применений (табл. 1).

Таблица 1. Локальные и глобальные параметры постобработки магнитостатических результатов решения во Flux 2D/3D
Параметры Единицы
измерения
Объяснение
Локальные 2D-параметры
Векторный потенциал:
•для плоской задачи — AZ,
•для осесимметричной задачи — rAZ
Вб/м
Вб
 
Напряженность магнитного поля, vec{H} А/м vec{H} = ν×rotAZ
Плотность магнитного потока, vec{B} Тл vec{B} = [μ]vec{H}
Магнитный поток, Ф
(phiP1P2 — магнитный поток через поверхность,
неограниченную линиями, нормальными к домену
исследования, передаваемый посредством точек P1 и P2)
Вб phiP1P2= (AZP1AZP2)p,
где p — толщина домена плоской задачи
Проницаемость, μ Гн/м  
Плотность тока источника, JZ А/мм2  
Лапласова сила (локальное значение), vec{F}L Н/м3 vec{F}L = vec{J}×vec{B}
Глобальные 2D-параметры
Магнитная сила, vec{F}m(Force) Н vec{F}m=∂Wm/∂l, l — элементарное линейное перемещение
Магнитная сила, vec{F}m(Magnetic Pressure) Н Интегрирование поверхностных сил
(тензорный метод Максвелла)
Магнитный крутящий момент, vec{Gamma}m(Torque) Н/м vec{Gamma}m=∂Wm/∂θ, θ — элементарное угловое перемещение
Лапласова сила, vec{F}L(stress) Н vec{F}L=VJ×vec{B}dV
Магнитная энергия, Wm Дж Wm=VdWmdV=V∫(0vec{B}vec{H}dvec{B})dV.
Для линейной изотропной гомогенной среды:
dWm= (1/2)vec{H}vec{B}= (1/2)μH2= (1/2)(B2/μ)
Магнитная энергия Energy/Co-Energy Дж V∫(0vec{H}vec{B}dvec{H})dV
Ток (источника), I(Current) А I=SJZdS
Магнитный поток в области,
Ф(Inductance/Flux)
Вб Ф = lAZdl
Индуктивность, L(Inductance/Inductance) Гн L = (n2Ф)/I
Локальные 3D-параметры
Скалярные магнитные потенциалы, phitot,phiredHj,phiredT0 А/м  
Напряженность магнитного поля, vec{H} А/м vec{H} = −vec{grad}(phitot), vec{H} = −vec{grad}(phiredHj) + vec{H}j, vec{H} = −vec{grad}(phiredT0) + vec{T}0
Плотность магнитного потока, vec{B} Тл vec{B} = [μ]vec{H}
Магнитное поле источников,
не покрываемых сеткой, vec{H}j
А/м Аналитическое или полуаналитическое вычисление
с использованием закона Био-Савара-Лапласа
Остаточная плотность магнитного
потока постоянных магнитов, vec{B}r
Тл  
Коэрцитивная сила постоянных магнитов, vec{H}c А/м  
Проницаемость, μ Гн/м  
Относительная проницаемость, μr    
Плотность поверхностной силы
между материалами, dvec{F}s
Н/м2 Вычисление на границах сред
с различной магнитной проницаемостью
Плотность источника тока JS в сетчатых обмотках А/м2 Верификация данных
Плотность лапласовой силы dvec{F}L в сетчатых обмотках Н/м3 dvec{F}L = vec{J}S×vec{B}
Объемная плотность магнитной энергии, dWm Дж/м3 dWm=0vec{B}vec{H}dvec{B}
Глобальные 3D-параметры
Магнитная сила, vec{F}m(Force) Н vec{F}m=∂Wm/∂l
Магнитный крутящий момент, vec{Gamma}m(Torque) Н/м vec{Gamma}m=∂Wm/∂θ
Магнитная энергия, Wm Дж Wm=VdWmdv
Магнитный поток через обмотку-проводник (stranded-типа), ФВ Вб  

Моделирование магнитогармонических процессов

Магнитные процессы, достигшие установившегося гармонического состояния, характеризуются одной и той же частотой, но синусоидальным изменением во времени всех физических параметров в уравнениях Максвелла.

Для изучения магнитных процессов, находящихся в установившемся гармоническом состоянии с постоянной частотой во Flux существует применение Steady State AC Magnetic. Магнитное поле в этом применении анализируется вместе с вихревыми токами, наведенными в проводящих областях.

Учитываются скин-эффект и эффект близости. Скин-эффект состоит в том, что электрические переменные токи неоднородно распределяются в поперечном сечении проводника — с преобладанием концентрации токовых линий в прилегающей к внешней стороне поверхности проводника, и этот эффект становится более значительным по мере того, как частота f тока увеличивается. Также скин-эффект становится более заметным с высокой магнитной проницаемостью μ материала и электрической проводимостью σ.

Параметры, ассоциируемые с электромагнитным полем, уменьшаются по мере увеличения глубины z проводника по экспоненциальному закону:

U (z) = US× ez, (35)

где US — значение на поверхности; ω = 2 πf — круговая частота; глубина скин-эффекта:

Эффект близости относится к явлению, которое возникает в системе из нескольких токонесущих проводников, когда магнитное поле, сгенерированное проводником, разрушает токовое распределение в других проводниках.

Результаты магнитогармонического моделирования — те же, что и для магнитостатического приложения, плюс расчет параметров вихревых токов и потерь мощности вследствие эффекта Джоуля. Все вычисленные для данной частоты параметры — синусоидальные.

В контексте Steady State AC Magnetic все параметры выражаются посредством их комплексных изображений. Такое представление действительно при соблюдении следующих гипотез:

  • все материалы линейные;
  • все источники синусоидальные.

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

  • для нелинейных материалов — параметры являются несинусоидальными, но предусмотрена возможность задания так называемой эквивалентной B(H)-характеристики, допускающей принимать в расчет факт синусоидальности параметров;
  • для линейных материалов, но несинусоидальных источников может применяться разложение в ряд Фурье — для получения суперпозиции нескольких гармонических задач (также этот метод применим и для нелинейных материалов).

Важной особенностью приложения является отсутствие в моделях магнитов. Вихревые токи наводятся и в пассивных, и в активных (питаемых внешней схемой) проводниках. Возможно объединение с внешней схемой.

Уравнения и условия вычислений следующие. В вычисления включаются уравнения Максвелла для магнитной системы (3), а соотношения материалов, которые принимаются в расчет, имеют вид (4).

Вычисляются B-, H-, E-поля. D-поле не вычисляется. Переменные состояния зависимы от времени (то есть d/dt ≠ 0).

Для решения используются две модели:

  • векторная, с вычислением переменной состояния магнитного векторного потенциала vec{A};
  • скалярная, с вычислением магнитного скалярного потенциала, записываемого как phitot или phired.

Уравнение Максвелла-Фарадея допускает присутствие электрического скалярного потенциала φ:

vec{E} = −(∂vec{A}/∂t)− vec{grad}φ. (37)

Уравнение для переходных процессов, решаемое в применении Transient Magnetic, записывается как

rot (v0[ vr] rot ( vec{A}))+[ σ ](∂vec{A}/∂t + vec{grad}φ) = 0. (38)

Комплексное изображение этого уравнения, в котором d/dt = jω, имеет следующий вид:

rot (v0[ vr] rot (vec{underline{A}})) + [ σ ](j ω vec{underline{A}} + vec{grad}φ) = 0. (39)

Переменные состояния включают:

  • магнитный векторный комплексный потенциал vec{A}( underline{A}n для плоских задач и r underline{A}n для осесимметричных);
  • электрический скалярный потенциал.

Flux 3D для скалярной модели Steady State AC Magnetic использует также различные автоматические формулировки с различными скалярными потенциалами — тех же, что рассматривались ранее в описании магнитостатических приложений.

Описание задачи в Steady State AC Magnetic включает определение частоты (при выборе применения). Для изучения низкочастотных проблем уравнения электрических полей E и D и магнитных полей B и H также могут быть разъединены. Это разделение зависит от частоты, материалов, размера домена исследования. Обычно разделение действительно от 0 до 10 ГГц.

О моделировании скин-эффекта необходимо сказать еще несколько слов.

Если материал характеризуется линейными, гомогенными и изотропными свойствами, среднее значение глубины скин-эффекта можно вычислять согласно формуле (36).

Физические параметры электрического тока или магнитного поля экспоненциально уменьшаются при увеличении глубины скинэффекта, если область материала плоская. Для того чтобы точно оценить физические параметры в скин-регионе, по крайней мере два слоя элементов (второго порядка) в толщине необходимы. Некоторые проблемы возникают, когда глубина δ << L или δ >> W — длины или ширины материала. Но во Flux существует возможность задания областей материала с поверхностным импедансом.

Гипотеза основывается на существовании нулевого вектора электрического поля в проводнике vec{E} — за исключением поверхности проводящего региона, где присутствует поверхностная плотность тока vec{j}S.

Поверхностный импеданс ZS соотносит тангенциальный компонент магнитного поля с тангенциальным компонентом электрического поля на поверхности проводящего региона с использованием следующего соотношения:

vec{E}S = ZS( vec{n}vec{H}S ), c ZS= (1 + j)/σδ, (40)

где vec{n} — единичный вектор, направленный нормально к поверхности, исходящий из проводящего региона.

Ограничения по определению поверхностного импеданса задаются двумя условиями: δ < L/5 и δ < r/5, где L — толщина, r — радиус кривизны.

Поле должно быть преимущественно тангенциальным (допускаются только незначительные нормальные составляющие). Эти ограничения записываются в виде: δ/L > 10(μ0/μ) и μ > 10μ0.

Проводящий твердотельный регион, описываемый условием поверхностного импеданса (Region of solid conductor type described by surface impedance), представляет собой комбинацию:

  • неактивного объемного региона,
  • контура региона с поверхностным импедансом.

Неактивный регион можно не покрывать сеткой, поскольку переменные состояния на внутренних узлах не вычисляются.

Flux 3D допускает также возможность моделирования тонких проводящих регионов пассивных проводников, в которых плотность тока моделируется постоянной (2D) или гиперболически распределенной по толщине (3D), а также областей с малым поперечным сечением — филиформ.

Источники в применении Steady State AC Magnetic бывают следующих типов:

  • электрический ток (синусоидальный);
  • наложенное значение напряженности H0 магнитного поля (3D-применение) однородное или распределенное, синусоидально изменяющееся во времени;
  • наложенное значение плотности φ0 магнитного потока (2D-применение), характеризующееся теми же свойствами, что и H0.

Параметры постобработки в Steady State AC Magnetic — скалярные или векторные, также синусоидальные.

Синусоидальный скалярный параметр угловой частоты пульсаций (переменного тока) ω может быть геометрически представлен как вектор, вращающийся с угловой скоростью ω (рис. 3) [8]. В комплексном домене вращающийся вектор для данного отрезка времени является представлением комплексного числа. Значение вращающегося вектора для t = 0 представляет собой комплексное изображение синусоидально изменяющегося скалярного параметра.

Рис. 3. Синусоидальный ток и его комплексное изображение
Рис. 3. Синусоидальный ток и его комплексное изображение

Мгновенное значение синусоидального тока равно изображению комплексного числа underline{I}ejωt:

i (t) = Î sin(ω t + β), (41)

где underline{I} = Îejβ (42) — это так называемое комплексное изображение тока, то есть i(t) = Imag( underline{I}ejβ) (43).

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

  • в декартовой форме:
  • i(t) = Î cos(ω t + β) + jÎ sin(ω t + β); (44)
  • в экспоненциальной форме:
  • underline{i} = Iejt + β), (45)

    где I — модуль комплексного значения underline{i}; (ω t + β) — аргумент или фаза комплексного значения underline{i}, β — начальная фаза (для t = 0).

Вращающийся вектор, назначенный синусоидальному параметру

A(t) = Â sin(ω t + v0), (46)

выражается в форме комплексной функции времени A(t). Эта комплексная функция времени может быть записана:

  • в декартовой форме:
  • A(t) =  cos(ω t + v0) + j  sin(ω t + v0); (47)
  • в экспоненциальной форме:
  • underline{A}(t) = Â ejt + v0), (48)

    где далее она может быть разложена на два фактора:

    underline{A}(t) = Â ejt + v0) = Â ejv0ejωt, (49)
    underline{A} =  ejv0 =  cos(v0) + j  sin(v0), (50)

    фазор (phasor) содержит информацию об амплитуде и начальной фазе синусоидально изменяющегося параметра — это и есть комплексное изображение параметра; ejωt = cos ω t + j sin ω t (51) — информацию о временной вариации параметра.

Синусоидальная функция во временном домене может быть представлена в комплексном домене с помощью первого фактора — комплексного изображения функции. Фазор (phasor) определяется как вращающийся вектор в момент t = 0 в комплексном домене, назначенный параметру A(t).

Постобработка скалярных результатов осуществляется различными путями:

  • мгновенная постобработка — для момента времени ω t : значение параметра A(t); значение комплексного вращающегося вектора (модуль/аргумент (фаза) или действительная/мнимая часть);
  • не зависимая во времени постобработка (во Flux 2D): фазор или комплексное изображение (модуль/аргумент (фаза) или действительная/мнимая часть).

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

BX(t) = hat{B}X sin(ω t + v0X), (52)

или, в комплексной системе изображения, как вращающийся вектор:

underline{B}X(t) = hat{B}X × ejt + v0X),
underline{B}X(t) = hat{B}X × cos(ω t + v0X) + jhat{B}X ×
× sin(ω t + v0X), (53)

где hat{B}X — модуль комплексного компонента underline{B}X; (ω t + v0X) — аргумент или фаза underline{B}X.

В общем случае векторные параметры описывают эллиптические пути периода в 2D-вычислительном домене (X, Y). Выражения двух компонентов вектора плотности магнитного потока B в точке вычислительного домена записываются как

BX(t) = hat{B}X ×cos(ω t + v0X),
BY(t) = hat{B}X ×cos(ω t + v0Y) (54)

и представляют собой не что иное, как параметрические уравнения эллипса — эллиптического изменения параметра во времени t (рис. 4) [8].

Рис. 4. Эллиптическое представление  двух компонентов вектора плотности  магнитного потока в 2D-вычислительном домене
Рис. 4. Эллиптическое представление
двух компонентов вектора плотности
магнитного потока в 2D-вычислительном домене

В зависимости от их эллиптической характеристики два компонента плотности B магнитного потока выражаются как следующее:

Постобработка векторных результатов возможна в следующих вариантах:

  • мгновенная постобработка — для данного момента времени ω(t):
    • значение векторного параметра vec{B}(t);
    • значение комплексного вращающегося вектора vec{underline{B}}(t) (модуль/аргумент (фаза)) или действительная/мнимая часть);
  • не зависимая от времени постобработка:
    • электрическое представление (в 2D): эллиптические характеристики (a, b и другие.)

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

Мгновенная электрическая мощность как сумма двух компонентов P(t) — абсорбируемая или поставляемая посредством диполя — представляет собой сумму двух компонентов:

  • не зависимого во времени:
  • Pm= U × I × cosφ;(56)
  • периодического компонента с двойной частотой:
  • Pp= U × I × cos(2ωt + φ). (57)

Мгновенное значение силы на лопастях контактора — это также сумма двух компонентов (компонента, не зависимого от времени Fm, и периодического):

Fm× cos(2ωt + φ). (58)

Для постобработки скалярного параметра, такого как мощность, главный интерес представляет вычисление постоянного компонента (среднего значения). Условие пульсации малоинтересно. Вычисление мгновенного значения мощности Pm(t) не выполняется ни во Flux 2D, ни во Flux 3D. Еще одна возможность для вычислений — временное интегрирование:

Для векторных параметров, таких как сила, представляют интерес именно два компонента — среднее значение и параметр пульсаций. Вычисляются как мгновенные значения, например F(t), так и среднее значение — посредством временной интеграции:

производится постобработка условия пульсации vec{F}P: стандартная постобработка, эллиптическое 2D-представление.

Локальные и глобальные параметры, доступные для постобработки/анализа результатов решения задачи в приложении Steady State AC Magnetic Flux, включают представленные в таблице 2.

Таблица 2. Локальные и глобальные параметры постобработки магнитогармонических результатов во Flux 2D/3D
Параметры Единицы
измерения
Объяснение
Локальные комплексные 2D-параметры
Векторный потенциал:
•для плоской задачи — underline{A}Z,
•для осесимметричной задачи — r underline{A}Z
Вб/м
Вб
 
Напряженность магнитного поля, bar{underline{H}} А/м bar{underline{H}} = v×rotunderline{A}Z
Плотность магнитного потока, bar{underline{B}} Тл bar{underline{B}} = [μ]bar{underline{H}}
Магнитный поток, Ф(underline{phi}P1P2 — магнитный поток через
поверхность, не ограниченную линиями, нормальными к домену
исследования, передаваемый посредством точек PP2)
Вб underline{phi}P1P2= (underline{A}ZP1underline{A}ZP2)p,
p — толщина домена плоской задачи
Плотность тока источника, underline{J}Z А/мм2 bar{underline{J}}Z= rot(bar{underline{H}})
Лапласова сила (локальное значение), bar{underline{F}}L Н/м3 bar{underline{F}}L= bar{underline{J}} × bar{underline{B}}
Локальные средние 2D-параметры
Объемная плотность активной мощности, dP Вт/м3 dP= (1/2)ρunderline{J}Z × underline{J}Z*
Другие локальные 2D-параметры
Проницаемость, μ Гн/м  
Удельное сопротивление, ρ Ом/м  
Глобальные комплексные 2D-параметры
Магнитная сила, vec{F}m(Magnetic Pressure) Н Интегрирование поверхностных сил
(тензорный метод Максвелла)
Магнитный крутящий момент, vec{Gamma}m(Torque) Н/м
Ток (источника), IZ(Current) А underline{I}Z=Sunderline{J}ZdS
Магнитный поток в области, Ф(Inductance/Flux) Вб Ф= lunderline{A}Zdl
Глобальные действительные 2D-параметры
Индуктивность, L(Inductance/Inductance) Гн L= (n2Ф)/I
Глобальные средние действительные 2D-параметры
Магнитная сила, bar{F}m Н bar{F}m= ∂Wm/∂l
Магнитный крутящий момент, bar{Gamma}m Н/м bar{Gamma}m= ∂Wm/∂θ
Активная мощность, P Вт P=V∫(dP)dV
Реактивная мощность, Q Вар Q=V∫(1/2)(ω/μ)|B|2dV
Лапласова сила, bar{underline{F}}L= V∫(bar{underline{J}} × bar{underline{B}})dV:
•средний действительный вектор;
•пульсирующий компонент.
Н  
Локальные комплексные 3D-параметры
Скалярные магнитные потенциалы, phitot, phiredHj, phiredT0 А/м  
Электрический векторный потенциал, bar{underline{T}} А/м Непрямой доступ
Напряженность магнитного поля, bar{underline{H}} А/м bar{underline{H}} = −

grad

(−phitot), bar{underline{H}} = −

grad

(−phiredHj) + bar{underline{H}}j, bar{underline{H}} = −

grad

(−phiredT0) + bar{underline{T}}0

Плотность магнитного потока, bar{underline{B}} Тл bar{underline{B}} = [μ]bar{underline{H}}
Электрическое поле, bar{underline{E}}(в твердотельном проводнике) В/м bar{underline{E}} = ((1/σ)rotbar{underline{T}})
Магнитное поле источников, не покрываемых сеткой, bar{underline{H}}j А/м Аналитическое или полуаналитическое вычисление
с использованием закона Био-Савара-Лапласа
Плотность источника тока в сетчатых обмотках, bar{underline{J}}S А/м2 Верификация данных
Плотность вихревых токов, bar{underline{J}}(в твердотельном проводнике) А/м2 bar{underline{J}} = rot(bar{underline{T}})
Другие локальные 3D-параметры
Проницаемость, μ Гн/м  
Проводимость, σ См  
Удельное сопротивление, ρ Ом·м  
Относительная проницаемость, μr    
Средние комплексные локальные 3D-параметры
Комплексное значение плотности мощности, dP+jdQ:
•действительное значение, dP;
•мнимое значение, jdQ
Вт/м3
Вар/м3
 
Средние действительные локальные 3D-параметры
Плотность потерь вследствие эффекта Джоуля
(объемная плотность активной мощности, dP)
Вт/м3 dP= (1/2)vec{E}×vec{J}
Объемная плотность реактивной мощности, dQ Вар/м3 dWm=0&intvec{B}vec{H}dvec{B}
Объемная плотность магнитной энергии, dWm Дж/м3  
Локальные параметры 2ωпульсации
Плотность поверхностной силы между материалами, dbar{F}s:
•среднее значение;
•пульсирующий компонент
Н/м2 Вычисление на границах сред
с различной магнитной проницаемостью
Объемная плотность лапласовой силы, dbar{F}L,
в проводниках stranded-типа/сетчатых обмотках:
•среднее значение;
•пульсирующий компонент
Н/м3 dbar{underline{F}}L = bar{underline{J}}S × bar{underline{B}}
Глобальные 3D-параметры (средние значения)
Магнитная сила, bar{F}m Н bar{F}m = ∂Wm/∂l
Магнитный крутящий момент, bar{Gamma}m Н/м bar{Gamma}m = ∂Wm/∂θ
Магнитная энергия, Wm Дж Wm= VdWmdv
Комплексное значение магнитного потока
через обмотку-проводник (stranded-типа), Фk
Вб Фk = ∂Wm/∂Ik = bar{B} × bar{h}ok,где hok — поле, созданное обмоткой k
в вакууме с током Ikв 1А, Ω — домен исследования

Переходные процессы Transient Magnetic

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

Transient Magnetic допускает схемное и кинематическое объединение.

Объединение с внешней схемой позволяет пользователю принимать в расчет наведенные токи в обмотках типа stranded (многожильных) — постоянные наведенные токи в проводниках без скин-эффекта. Кроме того, применение позволяет исследовать угловое движение и вихревые токи в пассивных твердотельных проводниках, в которых токи наводятся, но не проводятся.

В переходной магнитной модели действует гипотеза квазистатического состояния, согласно которой вариациями во времени электрического потока D можно пренебречь.

Принципиальные результаты, достижимые с этим применением, — те же, что и для магнитостатического применения. А также наведенные токи и рассеиваемая вследствие эффекта Джоуля мощность.

Уравнения, используемые для решения переходных проблем, представляют собой уравнения Максвелла для магнитной системы (3) и конституционные соотношения среды (4) — те же, что и для магнитогармонических процессов.

Условия вычислений следующие:

  • переменные состояния зависимы от времени d/dt ≠ 0;
  • вычисления производятся для B-, H-, E-полей, D-поле не вычисляется.

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

  • векторная (с магнитным векторным потенциалом A);
  • скалярная (с магнитным скалярным потенциалом phitot или phired).

Уравнение, решаемое в векторной модели применения Transient Magnetic методом конечных элементов, записывается в виде (38).

Переменные состояния включают:

  • магнитный векторный потенциал vec{A};
  • электрический скалярный потенциал φ.

Применение Transient Magnetic Flux 3D использует автоматические формулировки с различными скалярными потенциалами (в скалярной модели).

Описание задачи во многом сходно с магнитостатическим и магнитогармоническим применением. Определяются границы; граничные условия либо назначаются автоматически с использованием (в большинстве случаев) метода бесконечного бокса (в бесконечности магнитный потенциал нулевой), либо — без использования метода бесконечного бокса — задаются вручную с указанием нематериальных областей. В отсутствие контурных областей применяются стандартные граничные условия. Магнитное поле является тогда тангенциальным (электрическое — нормальным) к границе.

Если домену исследования назначаются симметрии, пользователь должен указать тип симметрии. Нормальное магнитное и тангенциальное электрическое поля или тангенциальное магнитное и нормальное электрическое поля могут быть наложены на плоскости симметрии, определенные в геометрическом контексте. Если домену исследования назначаются периодичности, пользователь должен указать тип периодичностей. Циклические или антициклические условия также могут накладываться на плоскости периодичности, определенные в геометрическом контексте.

Объемы, поверхности и линии материальных регионов допускают в Transient Magnetic моделирование следующих материальных сред:

  • воздуха или вакуума (μr= 1);
  • магнитомягких и магнитотвердых непроводящих материалов (μr> 1);
  • обмоток, проводящих ток, немагнитных (μл= 1);
  • твердых проводников (активных или пассивных, магнитных (μr> 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)
Параметры Единицы
измерения
Объяснение
Локальные 2D-параметры
Плотность вихревых токов, vec{J} А/мм2 Компонент JZ
Удельное сопротивление, ρ Ом·м  
Объемная плотность активной мощности, dP Вт/м3 dP= ρ×j2
Глобальные 2D-параметры
Ток, I(Current) — источника coil
и проводника типа Solid Conductor
А I=SJZdS
Активная мощность, P Вт P=V∫(dP)dV
Локальные 3D-параметры
Электрический векторный потенциал, vec{T} А/м Непрямая доступность
Плотность вихревых токов, vec{j} А/м2 vec{j} = rot(vec{T})
Напряженность электрического поля, vec{E} В/м vec{E} = (1/σ)vec{j}
Проводимость, σ См  
Удельное сопротивление, ρ Ом·м  
Объемная плотность потерь вследствие эффекта Джоуля, dP Вт/м3 dP= vec{E} × vec{j}, может также записываться как dP = ρ × j2= σ × E2
Глобальные 3D-параметры
Рассеянная мощность или потери вследствие
эффекта Джоуля в объемной области, Р
Вт P=V∫(dP)dV

Принципы вычисления потерь

Потери мощности в электромеханическом устройстве бывают трех видов:

  • магнитные потери в магнитных схемах;
  • потери вследствие эффекта Джоуля;
  • механические потери (за счет трения и вентиляции).

Потери мощности в магнитных материалах относятся к явлениям, ассоциируемым с временными вариациями магнитного поля.

Они классически подразделяются на гистерезисные потери — микроскопического происхождения, и потери вследствие токов Фуко — макроскопического.

Гистерезисные потери генерируются посредством временной вариации магнитной микроструктуры вещества, главным образом вследствие движения стен магнитных доменов.

Когда магнитное поле характеризуется циклической временной вариацией, магнитные материалы характеризуются геометрическим местом точек (B, H), представляющим собой замкнутую кривую. Для сниженных значений частоты (f < 1 Гц), замкнутая кривая B(H), называемая статическим гистерезисным циклом, не зависит от частоты. Соответствующие магнитные гистерезисные потери пропорциональны области статического гистерезисного цикла и частоте (квадрату частоты возбуждающего магнитного поля).

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

При увеличении частоты область замкнутой кривой B(H), называемая динамическим гистерезисным циклом, увеличивается (рис. 5).

Рис. 5. Магнитные потери и гистерезисный цикл
Рис. 5. Магнитные потери и гистерезисный цикл

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

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

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

  • гистерезис принимается в расчет на уровне магнитного закона B(H);
  • магнитные потери вычисляются непосредственно;

либо:

  • гистерезисом на уровне B(H) пренебрегают;
  • вычисление магнитных потерь выполняется апостериори посредством теоретических и экспериментальных формул.

Моделирование гистерезиса увеличивает время вычислений и требует большего объема памяти. Поэтому во Flux применяется второй подход.

Flux обеспечивает пользователю две модели вычислений магнитных потерь. В обоих случаях вычисление производится апостериори. Гистерезис не моделируется на уровне процесса решения, и B(H)-характеристика представляет собой однозначное соотношение.

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

Ph= kh× f × BηMAX, (61)

где η — коэффициент Штейнметца (1,6–2); f — частота; BMAX — пиковое значение плотности потока магнитного поля.

Потери, сгенерированные токами Фуко, вычисляются как

Pec= kecf2B2MAX. (62)

Магнитные потери во Flux вычисляются средствами формул Бертотти. Общие магнитные потери могут быть разделены на следующие три категории:

  • потери посредством гистерезиса (P1), пропорциональные частоте f, которые являются наиболее значительными компонентами магнитных потерь на низкой частоте;
  • потери вследствие токов Фуко (P2), пропорциональные квадрату частоты;
  • дополнительные потери (P3), пропорциональные f3/2.

Теория Бертотти дает выражение магнитных потерь как функции частоты и пикового значения плотности магнитного потока:

dP = c1B2MAXf + c2(BMAXf)2+ c3(BMAXf)3/2, (63)
P1= c1B2MAXf,
P2= c2(BMAXf)2,
P3= c3(BMAXf)3/2,

где c1, c2, c3 — соответствующие коэффициенты, BMAX — максимальное значение достижимой индукции.

Коэффициенты c2 и c3 выражаются посредством следующих соотношений:

с2 = (π2× σ × d2)/6, (64)

с3= √σ × G × V0× S, (65)

где d — толщина листа (пластины); G — безразмерная константа; S — поперечное сечение листа (пластины); V0 — постоянное поле, которое зависит от разницы коэрцитивной силы между двумя магнитными объектами.

В применении Steady State AC Magnetic используется комплексное представление. Объемная плотность средней мощности dPmoy записывается как

dPmoy= khfB2MAX + ((π2σ d2)/6)(BMAXf)2+
+ ke(BMAXf)3/2× 8,67, (66)

где kh, ke — коэффициенты потерь.

В вычислительном фрейме Flux это выражение записывается как

dPmoy= [khfB2MAX + ((π2σ d2)/6)(BMAXf)2+
+ ke(BMAXf)3/2× 8,67] kf, (67)

где kf — коэффициент филлинга — заполнения (близкий к единице), или электрической изоляции пластин магнитного сердечника.

Средняя мощность, рассеиваемая в объемном регионе, записывается как

Pmoy= ∫∫∫regdPmoydV. (68)

Важно отметить, что BMAX представляет собой пиковую индукцию, и в симуляциях, где вычисленное значение не соответствует пику (ротор и статор моделируются неподвижными, а удельное сопротивление моделируется переменным), вычисления магнитных потерь имеют ограничения. Так же и в том случае, когда нелинейная аппроксимация B(H) и насыщение вводятся средствами эквивалентной модели намагничивания, значения B могут изменяться.

В применении Transient Magnetic изучается переменный режим поля. В зависимости от магнитных потерь объемная плотность мощности будет записываться как

dP(t) = khfB2MAX + σ (d2/12)((dB/dt)(t))2+
+ ke((dB/dt)(t))3/2, (69)

в вычислительном фрейме Flux записываемый с коэффициентом заполнения:

dP(t) = [khfB2MAX + σ (d2/12)((dB/dt)(t))2+
+ ke((dB/dt)(t))3/2] kf. (70)

Объемная плотность средней мощности за период dPmoy записывается как

которая далее становится во Flux

Средняя мощность, рассеиваемая в объемном регионе, записывается в соответствии с (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)=
= Hstatic(B) + Hdynamic(B, dB/dt). (73)

Кривая 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.

Добавить комментарий

Ваш адрес email не будет опубликован. Обязательные поля помечены *