ANSYS电磁场分析指南磁宏.docx
- 文档编号:14470471
- 上传时间:2023-06-23
- 格式:DOCX
- 页数:27
- 大小:1,006.66KB
ANSYS电磁场分析指南磁宏.docx
《ANSYS电磁场分析指南磁宏.docx》由会员分享,可在线阅读,更多相关《ANSYS电磁场分析指南磁宏.docx(27页珍藏版)》请在冰点文库上搜索。
ANSYS电磁场分析指南磁宏
第十一章磁宏
11.1什么是电磁宏
电磁宏是ANSY宏命令,其主要功能是帮助用户方便地建立分析模型、方便地获取想要观察的分析结果。
目前,ANSYSI供了下列宏命令,可用于电磁场分析:
•CMATRIX计算导体间自有和共有电容系数
•CURR2:
计算二维导电体内电流
•EMAGERF计算在静电或电磁场分析中的相对误差
•EMF沿预定路径计算电动力(emf)或电压降
•FLUXV计算通过闭合回路的通量
•FMAGBC对一个单元组件加力边界条件
•FMAGSUM对单元组件进行电磁力求和计算
•FOR2D计算一个体上的磁力
•HFSWEEP在一个频率范围内对高频电磁波导进行时谐响应分析,并进行
相应的后处理计算
•HMAGSOLN定义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力矩和虚功力矩求和
本章对这些宏有详细描述。
在ANSY命令手册和理论手册对这些宏有更详细的描述。
下面的表格列出了这些电磁宏的使用范畴。
电磁宏
刖
处
理
器
求
解
后处理器
时间历程
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
-
EMF1]
-
-
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
-
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
-
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连接的MVF区域
11.2使用电磁宏
电磁宏根据其实现的功能,可以分为如下四类:
•建模类•求解类
Thickness=DZ
圏1跑道型銭圈电流源
•后处理类•高频分析类磁宏
11.2.1
建模类
有三个宏可用作帮助建模:
RACEPERBC2和FMAGBC
1)RACE^生一个由条形和弧形基元(S0URCE3单元)组成的“跑道”形电流源
命令:
RACE
GUI:
MainMenu>Preprocessor>-Modeling-Create>RacetrackCoil
Main
Menu>Preprocessor>Loads>-Loads-Apply>-Magnetic-Excitation>RacetrackCoil
RACE宏要求的参数如前面图1所示。
“跑道”由二个参数XC和YC定位,这些值是在工作平面内分别沿X和丫轴到线圈厚度中点的距离。
执行该宏时,可以把构成线圈的这些SOURCE366元定义为一个部件,将部件名作为该宏的一个输入参数即可。
2)PERBC2宏通过生成两个周期性对称面所必须的约束方程或节点耦合来施加周期对称边界条件,调用该宏的方式如下:
命令:
PERBC2D
GUI:
Main
Menu>Preprocessor>Loads>-Loads-Apply>-Magnetic-Boundary>-Vector
Poten-PeriodicBCs
MS
R1
R2
AZfx.y1)=~AZ(xjy2):
奇衬称AZ(x,y1)=AZ(xjy2).1对称
■H/ors
图4两个周期性对称面平行于X轴—*
奇对称选项表示一个半周期对称条件,偶对称条件表示全周期对称条件(重复结构)。
3)FMAGB用于对单元组件施加Maxwell面标志和虚功边界条件:
命令:
FMAGBC
GUI:
MainMenu>Preprocessor>Loads>-Loads-Apply>-Magnetic-Flag>CompForce/Torq
MainMenu>Solution>Apply>Comp.Force/Torq
11.2.2求解类
有四个宏可帮助求解:
MAGSOLV,HMAGSOLV,CMAT和IXMATRIX
11.2.2.1MAGSOLV
MAGSOL宏对大多数静磁分析问题能很快地定义求解选项并开始求解。
它可应用于2D和3D模型,标量势法、矢量势法和棱边单元求解方法,以及线性和非线性分析。
该宏不需要用户使用MAGOP命令和二步或三步求解顺序(一定情况所要求),它也允许定义非线性收敛标准,并提供选项来控制电流源Biot-Savart积分的重新计算。
命令:
MAGSOLV
GUI:
MainMenu>Solution>-Solve-Electromagnet>-Static
Analysis-Opt&Solv
11.2.2.2HMAGSOLV
HMAGSO宏对谐波分析能很快定义求解选项和开始求解。
对2-D模型,它使
用磁矢势(MVP方法求解。
它可用于线性和非线性分析。
对于非线性分析,此宏不需要定义二步求解过程,并允许用户自己定义收敛标准。
用下列方法之一,调用HMAGSOL^:
命令:
HMAGSOLV
GUI:
MainMenu>Solution>-Solve-Electromagnet>-HarmonicAnalys-Opt&Solv
11.2.2.3CMATRIX
CMATRI河计算“对地”和“集总”电容矩阵。
“对地”电容值表示一个导体的电荷与导体对地电压之比。
“集总”电容值表示二个导体之间的电容值。
实例详见多导体系统求取电容和本手册电容计算实例(命令方法)。
详细情况见ANSYSS论手册第五章。
用下面方法之一,调用CMATRIX宏命令:
命令:
CMATRIX
GUI:
MainMenu>Solution>Solve-Electromagnet>CapacMatrix
11.2.2.4LMATRIX
LMATRIX宏可以计算任意线圈组中每个线圈的微分电感矩阵和总磁链。
参见
《ANSYS1论手册》第5章。
LMATRIX宏用于在静磁场分析的一个“工作点”上计算任意一组导体间的微分电感矩阵和磁链。
“工作点”被定义为在系统上加工作(名义)电流所得到的解,该宏命令既可用于线性求解也可用于非线性求解。
必须用波前求解器来计算“工作点”的解。
LMATRIX宏的计算依赖于对工作点进行求解的过程中建立的多个文件。
该宏在执行求解之前在这些文件前面加一个前缀OPE眯重命名文件,并在完成求解后自动保存这些文件。
用户自己也可以保存这些文件的拷贝以进行备份。
该宏命令返回一个NXN+1矩阵参数,NXN部分表示N-绕组系统的微分电感值,此处N表示系统中的线圈数。
N+1列表示总磁链。
第I行表示第I个线圈。
另外,电感矩阵的值还以文本文件的格式输出,以供外部使用。
文件中第一个列表表示每个线圈的磁链。
第二个列表表示微分电感矩阵的上三角部分。
命令:
LMATRIX
GUI:
MainMenu>Solution>-Solve-Electromagnet>-StaticAnalysis-InductMatrix
在调用LMATRIX宏之前,还需要给线圈单元赋一个名义电流值。
对于使用磁矢势(MVP法或基于棱边元方法进行求解的静磁分析,可以使用BFVBFA或
BFE命令来给线圈单元赋名义电流(以电流密度的方式)。
对于使用简化标势法
(RSP、差分标势法(DSP和通用标势法(GSP的静磁分析,可以使用SOURCE36单元的实常数来给线圈单元赋名义电流。
为了使用LMATRIX宏,必须事先用*DIM命令定义一个N阶数组,N为线圈数,数组的每行都表示一个线圈。
数组的值等于线圈在工作点时每匝的名义电流值,且电流值不能为零,当确实有零电流时,可以用一个很小的电流值来近似。
另外,还需用CM命令把每个线圈的单元组合成一个部件。
每组独立线圈单元的部件名必须是用一个前缀后面再加线圈号来定义。
一个线圈部件可由标量
(RSP/DSP/GSP或矢量单元(MVP混合组成,最重要的一点是这些单元的激励电流与前面数组中所描述的电流相同。
在LMATRIX宏中需定义一个用于保存电感矩阵的数组名,用LMATRIX宏的对称系数(symfac)来定义对称性。
如果由于对称性而只建了n分之一部分模型,则计算出的电感乘以n就得到总的电感值。
当工作点位于BH曲线的弯点处时,切向磁导率变化最快,会导致计算的感应系数随收敛标准而变化。
为了获得更加准确的解,收敛标准要定义得更加严格一些,不仅仅是缺省值1.0X10一3。
一般在执行MAGSOL命令时,选择1.0X10一4或1.0X10-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的单元
cmwindl,elem!
给选出的单元赋予部件名windl
esel,s,,!
选择线圈2的单元
cmwind2,elem!
给选出的单元赋予部件名wind2
esel,s,,!
选择线圈3的单元
cmwind3,elem!
给选出的单元赋予部件名wind3
symfac=2!
对称系数
Imaxtrix,symfac,'wind','curr','ind'!
计算微分电感矩阵和总磁链
*stat,ind!
列出ind电感矩阵
11.2.2.5下面是以命令流方式进行的一个计算电感矩阵的例子
该例计算一个二线圈系统(永磁电感器件)在非线性工作点下的微分电感矩阵和总磁链,其示意图如下:
几何性质:
x仁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!
symmetricfactorforinductancecomputationnc=2!
numberofcoils
*dim,cur,array,nc!
nominalcurrentsofcoils
*dim,coils,char,nc!
namesofcoilcomponents
!
cur
(1)=0.25!
nominalcurrentof1stcoilcoils
(1)="wind1"!
nameofcoil1component!
cur
(2)=-0.125!
nominalcurrentof2ndcoilcoils
(2)="wind2"!
nameofcoil2component
!
auxiliaryparameters
mu0=3.1415926*4.0e-7
x3=x1+x2!
xcoordinaterighttocoil2leftx4=x3+2*x!
xcoordinaterighttocorex5=x4+x2!
xcoordinaterighttocoil2rightx6=x5+x1!
xcoordinaterighttocoil1rightjs1=cur
(1)*n1/(x1*y)!
nominalcurrentdensityofcoil1js2=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)A2;BS=2T;HS=100A/m
*do,qqq,1,20
B=qqq/10*Bstbpt,,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,1asel,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=0d,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,!
Storedelectromagneticenergysavelen=S_ENG
senergy,1!
Co-energy
savelce=C_ENG
fini
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
computeinductance
andflux
lmatrix,symfac,"wind","cur","ind",!
computeinductancematrix!
/com
finish
你将得到如下结果:
SUMMARYOFSTOREDENERGYCALCULATIONLoadStepNumber:
1.
SubstepNumber:
1.
Time:
0.1000E+01
MaterialNumberofStoredEnergyMaterialDescriptionNumberElements(J/m)
1.4.0.52360E-05LinearIsotrp...
2.1.-0.33314E+00Nonlin.MagnetIsotrp.
TOTAL5.-0.33313E+00
Note:
Theenergydensityfortheactiveelementsusedintheenergy
calculationisstoredintheelementitem"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-energydensityfortheactiveelements
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ANSYS 电磁场 分析 指南