- 积分
- 62
- 回帖
- 0
- 西莫币
-
- 贡献
-
- 威望
-
- 存款
-
- 阅读权限
- 5
- 最后登录
- 1970-1-1
该用户从未签到
|
!********************************************
! 说明:该例子演示一个永磁体的磁场(使用了infin9单元)
!********************************************
/TITLE, alextest, Test for Permanent Magnet
*go,:start !:start 是*go的行标签 即批处理要执行的下一行
:start !利用这个可以让ansys有选择性的读取输入文件 :start 行标签
!JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/PREP7
emunit,mks !定义电磁单位为国际标准单位,即μ0=4 Pi e-7 henries/meter
!********************************************
! 定义单元类型*
!********************************************
ET,53,PLANE53 ! Define PLANE 53 as element type
ET,9,INFIN9 !无限外界(注意:系统原点一定不能在infin9类型的节点上)
!********************************************
! 定义材料*
!********************************************
MP,MURX,1,1 !Define material properties (permeability) !定义空气(磁导率=1)
HC=895000 ! Coercive force!表示矫顽力 有的地方为895000
TB,BH,2,,30
TBPT,, 130.000000 , 0.100000000
TBPT,, 170.000000 , 0.200000000
TBPT,, 197.000000 , 0.300000000
TBPT,, 218.000000 , 0.400000000
TBPT,, 250.000000 , 0.500000000
TBPT,, 290.000000 , 0.600000000
TBPT,, 338.000000 , 0.700000000
TBPT,, 400.000000 , 0.800000000
TBPT,, 472.000000 , 0.900000000
TBPT,, 570.000000 , 1.00000000
TBPT,, 682.000000 , 1.10000000
TBPT,, 810.000000 , 1.20000000
TBPT,, 975.000000 , 1.30000000
TBPT,, 1600.00000 , 1.40000000
TBPT,, 2520.00000 , 1.50000000
TBPT,, 3520.00000 , 1.60000000
TBPT,, 4760.00000 , 1.70000000
TBPT,, 8300.00000 , 1.80000000
TBPT,, 12000.0000 , 1.90000000
TBPT,, 17000.0000 , 2.00000000
TBPT,, 23000.0000 , 2.10000000
TBPT,, 32000.0000 , 2.20000000
TBPT,, 42500.0000 , 2.30000000
TBPT,, 44500.0000 , 2.32000000
TBPT,, 48200.0000 , 2.35000000
TBPT,, 52200.0000 , 2.37000000
TBPLOT,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=20
hall=10
w1=4
h1=2
!x1=(wall-w1)/2
!x2=x1+w1
!y1=(hall-h1)/2
!y2=y1+h1
!rectng,0,wall,0,hall
!rectng,x1,x2,y1,y2
rectng,-wall/2,wall/2,-hall/2,hall/2
rectng,-w1/2,w1/2,-h1/2,h1/2
aovlap,all
numcmp,area !将生成的面重新编号
aplot
!/eof !应该是关闭 aplot 语句的
!建立材料属性
!********************************************
asel,s,area,,1 !选择中间的磁铁
aatt,2
asel,s,area,,2 !选择周围的空气
aatt,1
!********************************************
!建立单元类型,并划分网格
!********************************************
asel,all
type,53
lsel,s,line,,1,4 !选择所有的无限外边界
type,9 !设定为infin9单元,因为是开域问题则可以加个无限外边界来描述,当然也可以用无限单元infin10,则需要远场区域描述
lesize,all,,,30 !划分为30等份
lmesh,all !开始划分
lsel,s,line,,5,8
lesize,all,,,20 !将磁铁边界的每条边分成20等份
asel,all
amesh,all
!********************************************
!建立载荷
!********************************************
ESEL,ALL
NSEL,EXT
D,ALL,AZ,0
!********************************************
!求解
!********************************************
ALLSEL,ALL
MAGSOLV
!********************************************
!后处理
!********************************************
FINISH
/POST1
PLF2D,27,0,10,1 !显示磁力线
/IMAGE,SAVE,mf,JPEG !将磁力线保存成jpg文件
-------------------------------------------------------------
我对于开域问题的描述理解不是很到位,
lsel,s,line,,1,4 !选择所有的无限外边界
type,9
仅仅加个无限外边界能很好的描述开域问题吗?
难道不需要远场区域描述开域问题? |
|