ANSYS电磁场分析指南磁宏.docx

ANSYS电磁场分析指南磁宏.docx

《ANSYS电磁场分析指南磁宏.docx》由会员分享，可在线阅读，更多相关《ANSYS电磁场分析指南磁宏.docx（28页珍藏版）》请在冰豆网上搜索。

ANSYS电磁场分析指南磁宏

第十一章磁宏

11.1什么是电磁宏

电磁宏是ANSYS宏命令，其主要功能是帮助用户方便地建立分析模型、方便地获取想要观察的分析结果。

目前，ANSYS提供了下列宏命令，可用于电磁场分析：

·CMATRIX：

计算导体间自有和共有电容系数

·CURR2D：

计算二维导电体内电流

·EMAGERR：

计算在静电或电磁场分析中的相对误差

·EMF：

沿预定路径计算电动力（emf）或电压降

·FLUXV：

计算通过闭合回路的通量

·FMAGBC：

对一个单元组件加力边界条件

·FMAGSUM：

对单元组件进行电磁力求和计算

·FOR2D：

计算一个体上的磁力

·HFSWEEP：

在一个频率范围内对高频电磁波导进行时谐响应分析，并进行

相应的后处理计算

·HMAGSOLV：

定义2-D谐波电磁求解选项并进行谐波求解

·IMPD：

计算同轴电磁设备在一个特定参考面上的阻抗

·LMATRIX：

计算任意一组导体间的电感矩阵

·MAGSOLV：

对静态分析定义磁分析选项并开始求解

·MMF：

沿一条路径计算磁动力

·PERBC2D：

对2—D平面分析施加周期性约束

·PLF2D：

生成等势的等值线图

·PMGTRAN：

对瞬态分析的电磁结果求和

·POWERH：

在导体内计算均方根（RMS）能量损失

·QFACT：

根据高频模态分析结果计算高频电磁谐振器件的品质因子

·RACE：

定义一个“跑道形”电流源

·REFLCOEF：

计算同轴电磁设备的电压反射系数、驻波比、和回波损失

·SENERGY：

计算单元中储存的磁能或共能

·SPARM：

计算同轴波导或TE10模式矩形波导两个端口间的反射参数

·TORQ2D：

计算在磁场中物体上的力矩

·TORQC2D：

基于一个圆形环路计算在磁场中物体上的力矩

·TORQSUM：

对2-D平面问题中单元部件上的Maxwell力矩和虚功力矩求和

本章对这些宏有详细描述。

在ANSYS命令手册和理论手册对这些宏有更详细的描述。

下面的表格列出了这些电磁宏的使用范畴。

电磁宏

前处理器

求解

后处理器

时间历程

MVP域

MSP域

2D平面

2D轴对称

3D

棱边方法

高频

静态

谐波

瞬态

耦合

CMATRIX

Y

Y

-

-

-

-

Y

Y

Y

-

-

-

-

-

Y

CURR2D

-

-

Y

-

Y

-

Y

Y

Y

-

-

Y

Y

Y

-

EMAGERR

-

-

Y

-

Y

Y

Y

Y

Y

Y

-

Y

Y

Y

-

EMF[1]

-

-

Y

-

-

-

-

-

-

-

Y

-

Y

-

Y

FLUXV

-

-

Y

-

Y

-

Y

Y

Y

-

-

Y

Y

Y

-

FMAGBC

Y

Y

-

-

Y

Y

Y

Y

Y

Y

-

Y

-

Y

Y

FMAGSUM[1]

-

Y

Y

-

Y

Y

Y

Y

Y

-

-

Y

-

Y

Y

FOR2D

-

-

Y

-

Y

-

Y

Y

-

-

-

Y

-

Y

-

HFSWEEP

-

Y

-

-

-

-

-

-

-

-

Y

-

-

-

-

HMAGSOLV

-

Y

-

-

Y

-

Y

Y

-

-

-

-

Y

-

-

IMPD

-

-

Y

-

-

-

-

-

-

-

Y

-

-

-

-

LMATRIX

Y

Y

-

-

Y

Y

Y

Y

Y

Y

-

Y

-

-

-

MAGSOLV

Y

Y

-

-

Y

Y

Y

Y

Y

Y

-

Y

-

Y

-

MMF

-

-

Y

-

Y

Y

Y

Y

Y

Y

Y

Y

Y

Y

-

PERBC2D

Y

-

-

-

Y

-

Y

-

-

-

-

Y

Y

Y

-

PLF2D

-

-

Y

-

Y

-

Y

Y

-

-

-

Y

Y

Y

-

PMGTRAN

-

-

-

Y

Y

Y

Y

Y

Y

Y

-

-

-

Y

-

POWERH

-

-

Y

-

Y

-

Y

Y

Y

Y

Y

-

Y

-

-

QFACT

-

-

Y

-

-

-

-

-

-

-

Y

-

-

-

-

RACE

Y

-

-

-

-

Y

-

-

Y

-

-

Y

[2]

[2]

-

REFLCOEF

-

-

Y

-

-

-

-

-

-

-

Y

-

-

-

-

SENERGY

-

-

Y

-

Y

Y

Y

Y

Y

-

-

Y

Y

Y

-

SPARM

-

-

Y

-

-

-

-

-

-

-

Y

-

-

-

-

TORQC2D

-

-

Y

-

Y

-

Y

-

-

-

-

Y

Y

Y

-

TORQ2D

-

-

Y

-

Y

-

Y

-

-

-

-

Y

-

Y

-

TORQSUM

-

-

Y

-

Y

-

Y

-

-

-

-

Y

Y

Y

-

1这些宏也应用于静电场问题

2能用在通过界面单元INTER115连接的MVP区域

11.2使用电磁宏

电磁宏根据其实现的功能，可以分为如下四类：

·建模类·求解类

·后处理类·高频分析类磁宏

11.2.1建模类

有三个宏可用作帮助建模：

RACE、PERBC2D和FMAGBC

1）RACE产生一个由条形和弧形基元（SOURCE36单元）组成的“跑道”形电流源

命令：

RACE

GUI：

MainMenu>Preprocessor>-Modeling-Create>RacetrackCoil

MainMenu>Preprocessor>Loads>-Loads-Apply>-Magnetic-Excitation>RacetrackCoil

RACE宏要求的参数如前面图1所示。

“跑道”由二个参数XC和YC定位，这些值是在工作平面内分别沿X和Y轴到线圈厚度中点的距离。

执行该宏时，可以把构成线圈的这些SOURCE36单元定义为一个部件，将部件名作为该宏的一个输入参数即可。

2）PERBC2D宏通过生成两个周期性对称面所必须的约束方程或节点耦合来施加周期对称边界条件，调用该宏的方式如下：

命令：

PERBC2D

GUI：

MainMenu>Preprocessor>Loads>-Loads-Apply>-Magnetic-Boundary>-VectorPoten-PeriodicBCs

下面的图形描述了该宏的三种选项型式：

奇对称选项表示一个半周期对称条件，偶对称条件表示全周期对称条件（重复结构）。

3）FMAGBC用于对单元组件施加Maxwell面标志和虚功边界条件：

命令：

FMAGBC

GUI：

MainMenu>Preprocessor>Loads>-Loads-Apply>-Magnetic-Flag>CompForce/Torq

MainMenu>Solution>Apply>Comp.Force/Torq

11.2.2求解类

有四个宏可帮助求解：

MAGSOLV,HMAGSOLV,CMATRIX和LMATRIX

11.2.2.1MAGSOLV

MAGSOLV宏对大多数静磁分析问题能很快地定义求解选项并开始求解。

它可应用于2D和3D模型，标量势法、矢量势法和棱边单元求解方法，以及线性和非线性分析。

该宏不需要用户使用MAGOPT命令和二步或三步求解顺序（一定情况所要求），它也允许定义非线性收敛标准，并提供选项来控制电流源Biot-Savart积分的重新计算。

命令：

MAGSOLV

GUI：

MainMenu>Solution>-Solve-Electromagnet>-StaticAnalysis-Opt&Solv

11.2.2.2HMAGSOLV

HMAGSOLV宏对谐波分析能很快定义求解选项和开始求解。

对2-D模型，它使用磁矢势（MVP）方法求解。

它可用于线性和非线性分析。

对于非线性分析，此宏不需要定义二步求解过程，并允许用户自己定义收敛标准。

用下列方法之一，调用HMAGSOLV宏：

命令：

HMAGSOLV

GUI：

MainMenu>Solution>-Solve-Electromagnet>-HarmonicAnalys-Opt&Solv

11.2.2.3CMATRIX

CMATRIX可计算“对地”和“集总”电容矩阵。

“对地”电容值表示一个导体的电荷与导体对地电压之比。

“集总”电容值表示二个导体之间的电容值。

实例详见多导体系统求取电容和本手册电容计算实例（命令方法）。

详细情况见ANSYS理论手册第五章。

用下面方法之一，调用CMATRIX宏命令：

命令：

CMATRIX

GUI：

MainMenu>Solution>Solve-Electromagnet>CapacMatrix

11.2.2.4LMATRIX

LMATRIX宏可以计算任意线圈组中每个线圈的微分电感矩阵和总磁链。

参见《ANSYS理论手册》第5章。

LMATRIX宏用于在静磁场分析的一个“工作点”上计算任意一组导体间的微分电感矩阵和磁链。

“工作点”被定义为在系统上加工作（名义）电流所得到的解，该宏命令既可用于线性求解也可用于非线性求解。

必须用波前求解器来计算“工作点”的解。

LMATRIX宏的计算依赖于对工作点进行求解的过程中建立的多个文件。

该宏在执行求解之前在这些文件前面加一个前缀OPER来重命名文件，并在完成求解后自动保存这些文件。

用户自己也可以保存这些文件的拷贝以进行备份。

该宏命令返回一个N×N＋1矩阵参数，N×N部分表示N-绕组系统的微分电感值，此处N表示系统中的线圈数。

N＋1列表示总磁链。

第I行表示第I个线圈。

另外，电感矩阵的值还以文本文件的格式输出，以供外部使用。

文件中第一个列表表示每个线圈的磁链。

第二个列表表示微分电感矩阵的上三角部分。

命令：

LMATRIX

GUI：

MainMenu>Solution>-Solve-Electromagnet>-StaticAnalysis-InductMatrix

在调用LMATRIX宏之前，还需要给线圈单元赋一个名义电流值。

对于使用磁矢势（MVP）法或基于棱边元方法进行求解的静磁分析，可以使用BFV、BFA或BFE命令来给线圈单元赋名义电流（以电流密度的方式）。

对于使用简化标势法（RSP）、差分标势法（DSP）和通用标势法（GSP）的静磁分析，可以使用SOURCE36单元的实常数来给线圈单元赋名义电流。

为了使用LMATRIX宏，必须事先用*DIM命令定义一个N阶数组，N为线圈数，数组的每行都表示一个线圈。

数组的值等于线圈在工作点时每匝的名义电流值，且电流值不能为零，当确实有零电流时，可以用一个很小的电流值来近似。

另外，还需用CM命令把每个线圈的单元组合成一个部件。

每组独立线圈单元的部件名必须是用一个前缀后面再加线圈号来定义。

一个线圈部件可由标量（RSP/DSP/GSP）或矢量单元（MVP）混合组成，最重要的一点是这些单元的激励电流与前面数组中所描述的电流相同。

在LMATRIX宏中需定义一个用于保存电感矩阵的数组名，用LMATRIX宏的对称系数（symfac）来定义对称性。

如果由于对称性而只建了n分之一部分模型，则计算出的电感乘以n就得到总的电感值。

当工作点位于BH曲线的弯点处时，切向磁导率变化最快，会导致计算的感应系数随收敛标准而变化。

为了获得更加准确的解，收敛标准要定义得更加严格一些，不仅仅是缺省值1.0×10－3。

一般在执行MAGSOLV命令时，选择1.0×10－4或1.0×10－5。

在使用LMATRIX命令前，不要施加（或删除）非均匀加载，非均匀加载由以下原因生成：

·自由度命令（D,DA,等）在节点或者实体模型上定义非0值

·带有非0约束的CE命令

不要在不包含在单元组件中的单元上施加任何载荷（如current）

下面的例子是一个3线圈系统，每个线圈的名义电流分别为1.2、1.5和1.7安/匝，其分析的命令流如下。

在这个例子中，数组名为“curr”，线圈部件名前缀为“wind”，电感矩阵的计算值存贮在名为“ind”数组中。

值得注意的是，在LMATRIX命令行中，这些名字必须用单引号引起来。

*dim，cur，，3！

3个线圈系统数组

cur

（1）=1.2！

线圈1的名义电流为1.2安培/匝

cur

（2）=1.5！

线圈2的名义电流为1.5安培/匝

cur（3）=1.7！

线圈3的名义电流为1.7安培/匝

esel，s……！

选择线圈1的单元

cm，wind1，elem！

给选出的单元赋予部件名wind1

esel，s……！

选择线圈2的单元

cm，wind2，elem！

给选出的单元赋予部件名wind2

esel，s……！

选择线圈3的单元

cm，wind3，elem！

给选出的单元赋予部件名wind3

symfac=2！

对称系数

Imaxtrix，symfac，’wind’，’curr’，’ind’！

计算微分电感矩阵和总磁链

*stat，ind！

列出ind电感矩阵

11.2.2.5下面是以命令流方式进行的一个计算电感矩阵的例子

该例计算一个二线圈系统（永磁电感器件）在非线性工作点下的微分电感矩阵和总磁链，其示意图如下：

几何性质：

x1=0.1,x2=0.1,x=0.1,y=0.1

材料性质：

μr=1.0（空气），Hc=25（永磁体），B-H曲线（永磁体，见输入参数）

线圈1：

名义电流＝0.25安/匝，匝数＝10

线圈2：

名义电流＝0.125安/匝，匝数＝20

目标值：

L11=4,L22=16,L12=8

命令流如下：

/batch,list

/title,Two-coilinductorwithapermanentmagnet

/nopr

!

geometrydata

!

n=1!

meshingparameter

x=0.1!

width（xsize）ofcore

y=0.1!

　hightofcore,ysizeofwindow

z=1!

thicknessofironinzdirection

x1=0.1!

width（xsize）ofcoil1

x2=0.1!

width（xsize）ofcoil2

Hcy=25!

coercivemagneticfieldinydirection

n1=10!

numberofturnsincoil1

n2=20!

numberofturnsincoil2

!

!

excitationdatausedbyLMATRIX.MAC

!

symfac=1!

symmetricfactorforinductancecomputation

nc=2!

numberofcoils

*dim,cur,array,nc!

nominalcurrentsofcoils

*dim,coils,char,nc!

namesofcoilcomponents

!

cur

（1）=0.25!

nominalcurrentof1stcoil

coils

（1）="wind1"!

nameofcoil1component

!

cur

（2）=-0.125!

nominalcurrentof2ndcoil

coils

（2）="wind2"!

nameofcoil2component

!

!

auxiliaryparameters

!

mu0=3.1415926*4.0e-7

x3=x1+x2!

xcoordinaterighttocoil2left

x4=x3+2*x!

xcoordinaterighttocore

x5=x4+x2!

xcoordinaterighttocoil2right

x6=x5+x1!

xcoordinaterighttocoil1right

js1=cur

（1）*n1/（x1*y）!

nominalcurrentdensityofcoil1

js2=cur

（2）*n2/（x2*y）!

nominalcurrentdensityofcoil2

!

/prep7

et,1,53

!

mp,murx,1,1!

air/coil

mp,mgyy,2,Hcy!

coerciveterm

Bs=2!

saturationfluxdensity

Hs=100!

saturationmagneticfield

TB,BH,2!

core:

H=Hs（B/Bs）^2;BS=2T;HS=100A/m

*do,qqq,1,20

B=qqq/10*Bs

tbpt,,Hs*（B/Bs）**2,B

*enddo

!

rect,0,x1,0,y!

coil1left

rect,x1,x3,0,y!

coil2left

rect,x3,x4,0,y!

core

rect,x4,x5,0,y!

coil2right

rect,x5,x6,0,y!

coil1right

!

aglue,all

!

asel,s,loc,x,x1/2!

coil1volumeattribute

aatt,1,1,1

asel,s,loc,x,x5+x1/2

aatt,1,2,1

asel,s,loc,x,x1+x2/2!

coil2volumeattribute

aatt,1,3,1

asel,s,loc,x,x4+x2/2

aatt,1,4,1

asel,s,loc,x,x3+x!

ironvolumeattribute

aatt,2,5,1

asel,all

!

esize,,n

amesh,all

!

nsel,s,loc,x,x6!

　fluxparallelDirichletatsymmetryplain,x=x6

!

!

　homogeneousNeumannfluxnormalatyoke,x=0

d,all,az,0

nsel,all

!

esel,s,real,,1!

coil1leftcomponent

bfe,all,JS,,,,js1!

unitecurrentdensityincoil1

!

esel,s,real,,2!

coil1rightcomponent

bfe,all,JS,,,,-js1!

　returnunitecurrentdensityincoil1

!

esel,s,real,,1,2

cm,coils

（1）,elem

!

esel,s,real,,3!

coil2leftcomponent

bfe,all,JS,,,,js2!

unitecurrentdensityincoil2

!

esel,s,real,,4!

coil2rightcomponent

bfe,all,JS,,,,-js2!

　returnunitecurrentdensityincoil2

!

esel,s,real,,3,4

cm,coils

（2）,elem

!

!

allsel

!

fini

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

/com

/comobtainoperatingsolution

/com

!

/solu

cnvtol,csg,,1.0e-4

/out,scratch

solve

fini

!

/post1

!

/out

!

/com,

/com,

senergy,!

Storedelectromagneticenergy

savelen=S_ENG

senergy,1!

Co-energy

savelce=C_ENG

!

fini

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

!

computeinductance

lmatrix,symfac,"wind","cur","ind",!

computeinductancematrixandflux

!

/com

finish

你将得到如下结果：

SUMMARYOFSTOREDENERGYCALCULATION

LoadStepNumber:

1.

SubstepNumber:

1.

Time:

0.1000E+01

MaterialNumberofStoredEnergyMaterialDescription

NumberElements（J/m）

1.4.0.52360E-05LinearIsotrp...

2.1.-0.33314E+00Nonlin.MagnetIsotrp.

_____________________________________________________________________

TOTAL5.-0.33313E+00

Note:

Theenergydensityfor　theactiveelementsusedintheenergy

calculationisstored　intheelementitem"MG_ENG"fordisplay

andprinting.Thetotal　storedenergyissavedasparameter（S_ENG）

_____________________________________________________________________

本文系e-works专稿，XX严禁转载

SUMMARYOFCOENERGYCALCULATION

LoadStepNumber:

1.

SubstepNumber:

1.

Time:

0.1000E+01

MaterialNumberofCoenergyMaterialDescription

NumberElements（J/m）

1.4.0.52360E-05LinearIsotrp...

2.1.0.33314E+00Nonlin.MagnetIsotrp.

_____________________________________________________________________

TOTAL5.0.33314E+00

Note:

Theco-energydensityfor　theactiveelementsusedinthec