[转帖]载流铁磁导体非线性静态磁场分析
载流铁磁导体非线性静态磁场分析<BR>截面为园环的一段铸钢载有均匀分布的电流I,求导体内不同点处的磁感应强度B.模型的几何尺寸为:ri=.3 inro=.45 inr1=.75 in负载为: I=100A铸钢的BH曲线为:B.21.55.80.951.01.11.151.251.40H150300460640720890102012801900假定导体无限长,因此边端效应可以忽略,采用二维平面分析.由于分析区域轴对称,所以只选择5度来建模.外部空气取半径r1=.75in,并且采用无限边界单元INFIN9.求解时第一个载荷步作5次迭代,收敛标准为1*10.单位制采用MKS系统.转换因数为1in=.0254m. 电流密度j=i/a=438559(A/m2) <BR><BR>命令流：<BR>/VERIFY,VM165<BR>JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS<BR>/SHOW,JPEG<BR>/PREP7<BR>SMRT,OFF ?1<BR>/TITLE, VM165, CURRENT CARRYING FERROMAGNETIC CONDUCTOR<BR>C*** PRINCIPLES OF ELECTRIC AND MAGNETIC FIELDS, BOAST, PAGE 225<BR>ET,1,PLANE13 ! 2-D COUPLED FIELD SOLID<BR>ET,2,INFIN9 ! 2-D INFINITE BOUNDARY ELEMENT<BR>EMUNIT,MKS ! MKS UNITS<BR>MP,MURX,1,1 ! SET RELATIVE PERMEABILITY FOR AIR TO 1<BR>TB,BH,2 ! B-H CURVE FOR MATERIAL 2<BR>TBPT,,150,.21 ! H AND B RESPECTIVELY<BR>TBPT,,300,.55<BR>TBPT,,460,.80<BR>TBPT,,640,.95<BR>TBPT,,720,1.0<BR>TBPT,,890,1.1<BR>TBPT,,1020,1.15<BR>TBPT,,1280,1.25<BR>TBPT,,1900,1.40<BR>TBPLOT,BH,2<BR>TBPLOT,NB,2<BR>/wind,1,top<BR>/wind,2,bottom<BR>/gtype,1,grph,1<BR>/gtype,2,grph,1<BR>/gcmd,1,tbplot,bh,2<BR>/gcmd,2,tbplot,nb,2<BR>gplot<BR>/wind,2,off<BR>/wind,1,full ! 将BH曲线输出到jpeg图形中<BR>SF=.0254 ! SET CONVERSION (INCHES TO METERS)<BR>CSYS,1<BR>K,1<BR>K,2,.3,-2.5 ! INNER RADIUS OF RING<BR>K,3,.45,-2.5 ! OUTER RADIUS OF RING<BR>K,4,.75,-2.5 ! OUTER RADIUS OF SURROUNDING AIR<BR>KPSCALE,ALL,,,SF,,,,,1 ! MOVE ORIGINAL KEYPOINTS TO NEW POSITION<BR>CSYS,0<BR>L,1,2<BR>LESIZE,1,,,5 ! 分为5段<BR>L,2,3<BR>LESIZE,2,,,6<BR>L,3,4<BR>LSYMM,2,ALL<BR>NUMMRG,KP ! MERGE KEYPOINTS<BR>L,4,8 ! 定义最外面的圆弧<BR>LESIZE,7,,,1 ! 划分网格时只分为一段<BR>TYPE,2 ! 单元类型为INFIN9<BR>ESIZE,,,1<BR>LMESH,7 ! MESH LINE SEGMENT (WITH BOUNDARY ELEMENT)<BR>A,2,6,1,1<BR>A,2,3,7,6<BR>A,3,4,8,7<BR>ASEL,S,AREA,,2<BR>AATT,2 ! ASSIGN MATERIAL 2 TO STEEL AREA<BR>ASEL,ALL<BR>ESIZE,,1 ! DEFAULT ELEMENT DIVISIONS=1<BR>TYPE,1 ! 面的单元类型都设置为PLANE13<BR>AMESH,1,2<BR>LDVA ?2<BR>AMESH,3<BR>FINISH<BR>/SOLU ! ENTER SOLVER<BR>BFA,2,JS,,,438559 ! APPLY CURRENT DENSITY JS(Z) <BR>MAGSOLV<BR>FINISH<BR>/POST1<BR>CSYS,1<BR>NSEL,S,LOC,X,.325*SF<BR>NSEL,A,LOC,X,.375*SF<BR>NSEL,A,LOC,X,.425*SF<BR>RSYS,1 ! SET RESULTS C.S. TO CYLINDRICAL<BR>PRNSOL,B,COMP ! PRINT NODAL FLUX DENSITY<BR>*GET,B1,NODE,21,B,SUM<BR>*GET,B2,NODE,23,B,SUM<BR>*GET,B3,NODE,25,B,SUM<BR>*DIM,LABEL,CHAR,3,2<BR>*DIM,VALUE,,3,3<BR>LABEL(1,1) = \'@ R=\',\'@ R=\',\'@ R=\'<BR>LABEL(1,2) = \'.325\',\'.375\',\'.425\'<BR>*VFILL,VALUE(1,1),DATA,.48,1.03,1.22<BR>*VFILL,VALUE(1,2),DATA,B1,B2,B3<BR>*VFILL,VALUE(1,3),DATA,ABS(B1/.48),ABS(B2/1.03),ABS(B3/1.22)<BR>/COM<BR>/OUT,vm165,vrt<BR>/COM,------------------- VM165 RESULTS COMPARISON --------------<BR>/COM,<BR>/COM,B, TESLA<BR>/COM, | TARGET | ANSYS | RATIO<BR>/COM,<BR>*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)<BR>(1X,A8,A8,\' \',F10.2,\'\',F10.2,\' \',1F5.3)<BR>/COM,<BR>/COM,-----------------------------------------------------------<BR>/OUT<BR>FINISH<BR>*LIST,vm165,vrt<BR>
页:
[1]