整理ANSYS电磁场分析指南第十一章磁宏.docx
《整理ANSYS电磁场分析指南第十一章磁宏.docx》由会员分享,可在线阅读,更多相关《整理ANSYS电磁场分析指南第十一章磁宏.docx(33页珍藏版)》请在冰豆网上搜索。
整理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
意愿调查评估法(简称CV法)是指通过调查等方法,让消费者直接表述出他们对环境物品或服务的支付意愿(或接受赔偿意愿),或者对其价值进行判断。
在很多情形下,它是唯一可用的方法。
如用于评价环境资源的选择价值和存在价值。
-
疾病成本法和人力资本法是用于估算环境变化造成的健康损失成本的主要方法,或者说是通过评价反映在人体健康上的环境价值的方法。
-
(1)建设项目概况。
Y
(5)为保障评价对象建成或实施后能安全运行,应从评价对象的总图布置、功能分布、工艺流程、设施、设备、装置等方面提出安全技术对策措施;从评价对象的组织机构设置、人员管理、物料管理、应急救援管理等方面提出安全管理对策措施;从保证评价对象安全运行的需要提出其他安全对策措施。
对策措施的建议应有针对性、技术可行性和经济合理性,可分为应采纳和宜采纳两种类型。
-
-
Y
2.早期介入原则;[2]
[2]
-
4.建设项目环境影响评价文件的分级审批REFLCOEF
(2)评价范围。
根据评价机构专业特长和工作能力,确定其相应的评价范围。
-
表三:
周围环境概况和工艺流程与污染流程;-
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"ford