模拟永磁体的实例

1、模拟永磁体的实例看到论坛上很多讨论永磁体的，但是有些问题总是不够清楚。终于有一天我下决心好好看了ansys的相关帮助，自己做了一个最简单，也是最容易验证的例子：模拟一条磁铁的磁力线。下面就是我的输入文件。QUOTE:!*!说明：该例子演示一个永磁体的磁场（使用了infin9单元）!*/TITLE,alextest,TestforPermanentMagnet*go,:start:start!利用这个可以让ansys有选择性的读取输入文件!JPGPRF,500,100,1!MACROTOSETPREFSFORJPEGPLOTS/PREP7emunit,mks!定义电磁单位为国际标准单位，即科0=

2、4Pie-7henries/meter!*!定义单元类型*!*ET,53,PLANE53!DefinePLANE53aselementtypeET,9,INFIN9!无限外界（注意：系统原点一定不能在infin9类型的节点上）!*!定义材料*!*MP,MURX,1,1!Definematerialproperties（permeability）!定义空气（磁导率=1）HC=895000!Coerciveforce!表示矫顽力有的地方为895000TB,BH,2,30TBPT,130.000000,0.100000000TBPT,170.000000,0.200000000TBPT,197.00

3、0000,0.300000000TBPT,218.000000,0.400000000TBPT,250.000000,0.500000000TBPT,290.000000,0.600000000TBPT,338.000000,0.700000000TBPT,400.000000,0.800000000TBPT,472.000000,0.900000000TBPT,570.000000,1.00000000TBPT,682.000000,1.10000000TBPT,810.000000,1.20000000TBPT,975.000000,1.30000000TBPT,1600.00000,1.

4、40000000TBPT,2520.00000,1.50000000TBPT,3520.00000,1.60000000TBPT,4760.00000,1.70000000TBPT,8300.00000,1.80000000TBPT,12000.0000,1.90000000TBPT,17000.0000,2.00000000TBPT,23000.0000,2.10000000TBPT,32000.0000,2.20000000TBPT,42500.0000,2.30000000TBPT,44500.0000,2.32000000TBPT,48200.0000,2.35000000TBPT,5

5、2200.0000,2.37000000TBPLOT,BH,2,/IMAGE,SAVE,BH2,JPEG!将材料2的B-H曲线存储成bh2.jpgmp,mgxx,2,0!对于永磁铁，必须定义mgxx(或mgyy)!mp,mgyy,2,hc!*!*建立模型*!*/PNUM,AREA,1!定义显示模式wall=20hall=10w1=4h1=2!x1=(wall-w1)/2!x2=x1+w1!y1=(hall-h1)/2!y2=y1+h1!rectng,0,wall,0,hall!rectng,x1,x2,y1,y2rectng,-wall/2,wall/2,-hall/2,hall/2rectn

6、g,-w1/2,w1/2,-h1/2,h1/2aovlap,allnumcmp,area!将生成的面重新编号aplot!/eof!配合前面的:start使用!*!建立材料属性!*asel,s,area,1!选择中间的磁铁aatt,2asel,s,area,2!选择周围的空气aatt,1!*!建立单元类型，并划分网格!*asel,alltype,53lsel,s,line,1,4!选择所有的无限外边界type,9!设定为infin9单元lesize,all,30!划分为30等份lmesh,all!开始划分lsel,s,line,5,8lesize,all,20!将磁铁边界的每条边分成20等份as

7、el,allamesh,all*!建立载荷!*ESEL,ALLNSEL,EXTD,ALL,AZ,0!*!求解!*ALLSEL,ALLMAGSOLV!*!后处理!*FINISH/POST1PLF2D,27,0,10,1!显示磁力线/IMAGE,SAVE,mf,JPEG!将磁力线保存成jpg文件说明：1）如果输入文件写成：mp,mgxx,2,hcmp,mgyy,2,0那么，得到的磁力线为：（也就是磁极在x轴上。是第一幅图）ANSUSfl.0MAP9200519:53:33NODALSOIiUTIQM5TEP=ZSITE-1TIME=2AZRS¥S=0SMN=7.2625MX=2.Z62-2.17ST.Q1