模拟永磁体的实例.docx

模拟永磁体的实例 看到论坛上很多讨论永磁体的，但是有些问题总是不够清楚终于有一天我下决心好好看了 ansys 的相关帮助，自己做了一个最简单，也是最容易验证的例子：模拟一条磁铁的磁力线 下面就是我的输入文件QUOTE:!说明：该例子演示一个永磁体的磁场（使用了 infin9单元）/TITLE, alextest, Test for Permanent Magnet*go,:start:start !利用这个可以让ansys有选择性的读取输入文件!JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS /PREP7emunit,mks !定义电磁单位为国际标准单位，即p0=4 Pi e-7 henries/meter! 定义单元类型*ET,53,PLANE53 ! Define PLANE 53 as element typeET,9,INFIN9 !无限外界（注意：系统原点一定不能在infin9类型的节点上）! 定义材料*MP,MURX,1,1 !Define material properties （permeability） !定义空气（磁导率=1）HC=895000 ! Coercive force!表示矫顽力 有的地方为895000TB,BH,2,,30TBPT,, 130.000000 , 0.100000000TBPT,, 170.000000 , 0.200000000TBPT,, 197.000000 , 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.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,, 52200.0000 , 2.37000000TBPLOT,BH,2,,,/IMAGE,SAVE,BH2,JPEG !将材料 2 的 B-H 曲线存储成 bh2.jpg mp,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/2 rectng,-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,53 lsel,s,line,,1,4 !选择所有的无限外边界type,9 !设定为 infin9 单元lesize,all,,,30 !划分为30 等份lmesh,all !开始划分lsel,s,line,,5,8lesize,all,,,20 !将磁铁边界的每条边分成20 等份asel,allamesh,all!建立载荷ESEL,ALLNSEL,EXT D,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轴上。

是第一幅图）■3._Lcxte=：t,A1TSYS 8 ・ 0 :MAR Q 2 005 19:53:33NODAL SbLUTION 3TEP=2 SVB =1TIME=2A3RS ^3=^05MN ~-2.262 5MK =2,2 E2-2.01 -Lr0^3 -1^508 -L® -Lr 113 -^331609 -,6700S7 -.502565 -.1E7521 .L67E-Q5 ・335D46 .5DZS68 .67DD91 l.DDS 1.173 1.34 1.67S0 432）如果输入文件写成mp,mgxx,2,0mp,mgyy,2,hc那么，得到的磁力线为：（也就是磁极在y轴上是第二幅图）AW5TS 0 ,□ MXR 9 2005 19:54:2 EITODJIL SOLUTIOn STEP=2SUB =1 TIME-2 員是R51：S=0=-1.546 SMS =1・546-1.AB9 一1・375 -1,Z6 -「□J 1 一・ 91640^ -.80185^ 572755 450205 -r343 556 一.1145S6 637E-D5 .229092・ 3 436-lZ .458192 ■&87291 aaoia4i・ 9163 91451.261.4393) 如果输入文件写成mp,mgxx,2,hc*s in( 45)mp,mgyy,2,hc*cos(45)2 ・ 17E那么，得到的磁力线为：（也就是磁极在对角线上。

是第3幅图）AWSTTS Q . Q MAR 9 2DO5 za;31;0ZHODAJC SQI.UTIOM3TEP-2SUB =1TIMH=2 AZRST3-0 3MM =-L ・ 67@ 5MK =1.6^6—1 a b L G -1.^92 一：L・ 360 一1・119 -u994G15 -.a^oaea 621634 -.49730? -□3^29^9 一 ” 1 i 4 3 Z 5 -2USE-05 .Z 4E G36 ・ 37290-9 u49T31i .7^5955 -STDZ 32 r 99^3 6 19 lu2^3 1.2^0 1,616继续试验:如果输入文件写成mp,mgxx,2,0mp,mgyy,2,hc并且，将加载的命令注释掉：!d,all,az,0那么，得到的磁力线为：（也就是磁极在y轴上，并且没有加边界条件从磁力线的形状来看，似乎这个是最正确的？ ？ ？ ？ ？（高手来给解释一下！）AJTSZS 6 ・□MAR ID ZODS Zl:5B:3DNODAL SQIiUTION STEP=2SUB =1TIME=2A2alextest^ Test for Permanent MagrietP3TS-0 SMET =-1.6 69 SP-1Z =lu 669-la 607 -1.^03 一1・3石 -1・1L2 --988867 --36S259 61BMZ-,370826 -B12361-»166E-05 ・24T215 .370023 .43-5431 -7416-17 -8 652 55 ・988863 1-2361.3G1. GO 7接着实验：如果输入文件写成mp,mgxx,2,0mp,mgyy,2,hc并且，不使用infin9单元，而且也注释掉加载的命令：!d,all,az,0那么，得到的磁力线为：（也就是磁极在y轴上，并且没有加边界条件，也没有使用infin9 单元。

从磁力线情况来看，这种情况应该是不对的也就是说，应该要使用nfin9单元，结果才能 够比较正确alje-xtest. Test far Perraanen七 Magne七日・口 MAR 11 2005 io:zosaa MODAL SCLLUriO^ STEP-2 SUB =1 m