[转帖]载流铁磁导体非线性静态磁场分析

gghhjj

载流铁磁导体非线性静态磁场分析<BR>截面为园环的一段铸钢载有均匀分布的电流I,求导体内不同点处的磁感应强度B.模型的几何尺寸为:ri=.3 inro=.45 inr1=.75 in负载为: I=100A铸钢的BH曲线为:B  .21  .55  .80  .95  1.0  1.1  1.15  1.25  1.40H  150  300  460  640  720  890  1020  1280  1900假定导体无限长,因此边端效应可以忽略,采用二维平面分析.由于分析区域轴对称,所以只选择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>