ansys求电感的方法总结.doc

11.2.2.4 lmatrixlmatrix宏可以计算任意线圈组中每个线圈的微分电感矩阵和总磁链。参见ansys理论手册第5章。lmatrix宏用于在静磁场分析的一个“工作点”上计算任意一组导体间的微分电感矩阵和磁链。“工作点”被定义为在系统上加工作（名义）电流所得到的解，该宏命令既可用于线性求解也可用于非线性求解。必须用波前求解器来计算“工作点”的解。lmatrix宏的计算依赖于对工作点进行求解的过程中建立的多个文件。该宏在执行求解之前在这些文件前面加一个前缀oper来重命名文件，并在完成求解后自动保存这些文件。用户自己也可以保存这些文件的拷贝以进行备份。该宏命令返回一个nn1矩阵参数，nn部分表示n-绕组系统的微分电感值，此处n表示系统中的线圈数。n1列表示总磁链。第i行表示第i个线圈。另外，电感矩阵的值还以文本文件的格式输出，以供外部使用。文件中第一个列表表示每个线圈的磁链。第二个列表表示微分电感矩阵的上三角部分。命令：lmatrixgui：main menusolution-solve-electromagnet-static analysis-induct matrix在调用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.0103。一般在执行magsolv命令时，选择1.0104或1.0105。在使用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！给选出的单元赋予部件名wind1esel，s！选择线圈2的单元cm，wind2，elem！给选出的单元赋予部件名wind2esel，s！选择线圈3的单元cm，wind3，elem！给选出的单元赋予部件名wind3symfac=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-coil inductor with a permanent magnet/nopr!geometry data!n=1! meshing parameterx=0.1! width (x size) of corey=0.1!hight of core, y size of windowz=1! thickness of iron in z directionx1=0.1! width (x size) of coil 1x2=0.1! width (x size) of coil 2hcy=25! coercive magnetic field in y directionn1=10! number of turns in coil1n2=20! number of turns in coil2!excitation data used by lmatrix.mac!symfac=1! symmetric factor for inductance computationnc=2! number of coils*dim,cur,array,nc! nominal currents of coils*dim,coils,char,nc! names of coil components!cur(1)=0.25! nominal current of 1st coilcoils(1)=wind1! name of coil 1 component!cur(2)=-0.125! nominal current of 2nd coilcoils(2)=wind2! name of coil 2 component!auxiliary parameters!mu0=3.1415926*4.0e-7x3=x1+x2! x coordinate right to coil2 leftx4=x3+2*x! x coordinate right to corex5=x4+x2! x coordinate right to coil2 rightx6=x5+x1! x coordinate right to coil1 rightjs1=cur(1)*n1/(x1*y)! nominal current density of coil1js2=cur(2)*n2/(x2*y)! nominal current density of coil2!/prep7et,1,53!mp,murx,1,1! air/coilmp,mgyy,2,hcy! coercive termbs=2! saturation flux densityhs=100! saturation magnetic fieldtb,bh,2! core: h = hs (b/bs)2; bs=2t;hs=100a/m*do,qqq,1,20b=qqq/10*bstbpt,hs*(b/bs)*2,b*enddo!rect, 0,x1,0,y! coil1 leftrect,x1,x3,0,y! coil2 leftrect,x3,x4,0,y! corerect,x4,x5,0,y! coil2 rightrect,x5,x6,0,y! coil1 right!aglue,all!asel,s,loc,x,x1/2! coil 1 volume attributeaatt,1,1,1asel,s,loc,x,x5+x1/2aatt,1,2,1asel,s,loc,x,x1+x2/2! coil 2 volume attributeaatt,1,3,1asel,s,loc,x,x4+x2/2aatt,1,4,1asel,s,loc,x,x3+x! iron volume attributeaatt,2,5,1asel,all!esize,namesh,all!nsel,s,loc,x,x6!flux parallel dirichlet at symmetry plain, x=x6!homogeneous neumann flux normal at yoke, x=0d,all,az,0nsel,all!esel,s,real,1! coil 1 left componentbfe,all,js,js1! unite current density in coil 1!esel,s,real,2! coil 1 right componentbfe,all,js,-js1!return unite current density in coil 1!esel,s,real,1,2cm,coils(1),elem!esel,s,real,3! coil 2 left componentbfe,all,js,js2! unite current density in coil 2!esel,s,real,4! coil 2 right componentbfe,all,js,-js2!return unite current density in coil 2!esel,s,real,3,4cm,coils(2),elem!allsel!fini!/com/com obtain operating solution/com!/solucnvtol,csg,1.0e-4/out,scratchsolvefini!/post1!/out!/com,/com,senergy,! stored electromagnetic energysavelen=s_engsenergy,1! co-energysavelce=c_eng!fini! compute inductancelmatrix,symfac,wind,cur,ind, ! compute inductance matrix and flux!/comfinish你将得到如下结果：summary of stored energy calculationload step number:1.substep number:1.time:0.1000e+01material number ofstored energymaterial descriptionnumberelements(j/m) 1.4.0.52360e-05linearisotrp.2.1.-0.33314e+00nonlin.magnetisotrp._t o t a l5.-0.33313e+00note:the energy density forthe active elements used in the energycalculation is storedin the element item mg_eng for displayand printing. the totalstored energy is saved as parameter (s_eng)summary of coenergy calculationload step number:1.substep number:1.time:0.1000e+01material number ofcoenergymaterial descriptionnumberelements(j/m) 1.4.0.52360e-05linearisotrp.2.1.0.33314e+00nonlin.magnetisotrp._t o t a l5.0.33314e+00note:the co-energy density forthe active elements used in the co-energycalculation is storedin the element item mg_coeng for displayand printing. thetotal coenergy is saved as parameter (c_eng)_ lmatrix solution summary _flux linkage of coi