对于一个永磁铁的三维分析,分别采用标量法和棱边法计算,但差别很大,请问到底是什么原因?是因为电磁场的理论,还是我的命令流有错误,请大家讨论一下。
模型简介:
永磁铁:0.1mX0.1mX0.1m,剩磁1.12T,矫顽力780000,计算得到的相对磁导率为1.143
空气范围10mX10mX10m,相对磁导率为1
标量法:永磁铁和空气采用solid96,无限边界采用infin47
棱边法:永磁铁和空气采用solid117,无限边界采用infin47
提取坐标(3,3,3)的磁感应强度:
标量法:Bx=0.49831E-06
By=0.27287E-06
Bz=0.49831E-06
棱边法:Bx=0.54270E-06
By=-0.59086E-06
Bz= 0.54270E-06
!****棱边法****
/prep7
emunit,mks
!空气边界
b1=10
!磁铁尺寸
b0=0.1
!1个磁铁
!三维分析
!**************
!设置单元属性
!**************
!磁铁和空气
!低频棱边单元
et,1,solid117
keyopt,1,1,5
!远场边界
et,2,infin47
!**************
!设置材料属性
!**************
!空气
mp,murx,1,1.0
!永磁铁
mp,murx,2,1.143
mp,mgxx,2,0
mp,mgyy,2,780000
mp,mgzz,2,0
!**************
!建立几何模型
!**************
block,-b1/2,b1/2,-b1/2,b1/2,-b1/2,b1/2
wpoffs,,,-b0/2
vsbw,all,
wpcsys,-1
wpoffs,,,b0/2
vsbw,all,
wpcsys,-1
wpoffs,-b0/2
wprota,,,90
vsbw,all,
wpcsys,-1
wpoffs,b0/2
wprota,,,90
vsbw,all,
wpcsys,-1
wpoffs,,-b0/2
wprota,,90
vsbw,all,
wpcsys,-1
wpoffs,,b0/2
wprota,,90
vsbw,all,
wpcsys,-1
!**************
!建立有限元模型
!**************
!磁铁
csys,0
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
aslv,r
lslar,r
lesize,all,,,4
type,1
mat,2
vmesh,all
!6个小条空气
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,y,b0/2+0.01,b1/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,y,b1/2-0.01,b1/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,y,b0/2-0.01,b0/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,b0/2,b1/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,-b1/2,-b0/2,
aslv,r
lslar,r
lsel,r,loc,y,-b1/2+0.01,-b0/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,y,-b1/2+0.01,-b1/2-0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,y,-b0/2-0.01,-b0/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,-b1/2,-b0/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b1/2,-b0/2,
aslv,r
lslar,r
lsel,r,loc,z,-b1/2+0.01,-b0/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,z,-b1/2+0.01,-b1/2-0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,z,-b0/2-0.01,-b0/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b1/2,-b0/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,z,b0/2+0.01,b1/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,z,b0/2-0.01,b0/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,z,b1/2-0.01,b1/2-0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,b0/2,b1/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,x,-b1/2+0.01,-b0/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,x,-b0/2-0.01,-b0/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,x,-b1/2-0.01,-b1/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,-b1/2,-b0/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,x,b0/2+0.01,b1/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,x,b0/2-0.01,b0/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,x,b1/2-0.01,b1/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,b0/2,b1/2
type,1
mat,1
vmesh,all
!12个薄板空气
allsel
vsel,s,loc,z,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,y,-b0/2,b0/2
aslv,r
lslar,r
lsel,u,loc,z,-b0/2+0.01,b0/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,z,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,y,-b0/2,b0/2
aslv,r
lslar,r
lsel,r,loc,z,-b0/2+0.01,b0/2-0.01
lesize,all,,,4
allsel
vsel,s,loc,z,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,y,-b0/2,b0/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,u,loc,y,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
aslv,r
lslar,r
lsel,u,loc,x,-b0/2+0.01,b0/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,u,loc,y,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
aslv,r
lslar,r
lsel,r,loc,x,-b0/2+0.01,b0/2-0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,u,loc,y,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
aslv,r
lslar,r
lsel,u,loc,y,-b0/2+0.01,b0/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
aslv,r
lslar,r
lsel,r,loc,y,-b0/2+0.01,b0/2-0.01
lesize,all,,,4
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
type,1
mat,1
vmesh,all
!大部分空气
allsel
vsel,u,loc,y,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
aslv,r
lslar,r
lesize,all,,,20
allsel
vsel,u,loc,y,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
type,1
mat,1
vmesh,all
!远场边界
allsel
asel,s,loc,x,-b1/2-0.01,-b1/2+0.01
asel,a,loc,x,b1/2-0.01,b1/2+0.01
asel,a,loc,y,-b1/2-0.01,-b1/2+0.01
asel,a,loc,y,b1/2-0.01,b1/2+0.01
asel,a,loc,z,-b1/2-0.01,-b1/2+0.01
asel,a,loc,z,b1/2-0.01,b1/2+0.01
type,2
amesh,all
!**************
!求解
!**************
/solu
csys,0
allsel
!静力分析
antype,static
!施加边界条件
!棱边法
!垂直边界自动满足
!施加平行边界
allsel
nsel,s,loc,x,-b1/2-0.01,-b1/2+0.01
nsel,a,loc,x,b1/2-0.01,b1/2+0.01
nsel,a,loc,y,-b1/2-0.01,-b1/2+0.01
nsel,a,loc,y,b1/2-0.01,b1/2+0.01
nsel,a,loc,z,-b1/2-0.01,-b1/2+0.01
nsel,a,loc,z,b1/2-0.01,b1/2+0.01
d,all,az,0
!选择求解器
allsel
!矢量法或棱边法
magsolv,0
!**************
!后处理
!**************
/post1
!提取节点编号
csys,0
rsys,0
nx=3
ny=3
nz=3
noden=node(nx,ny,nz)
!提取磁感应强度
*get,bx1,node,noden,b,x
*get,by1,node,noden,b,y
*get,bz1,node,noden,b,z
*get,bs1,node,noden,b,sum
!提取磁场强度
*get,hx1,node,noden,h,x
*get,hy1,node,noden,h,y
*get,hz1,node,noden,h,z
*get,hs1,node,noden,h,sum
!写入到文件
*create,ansuitmp
*cfopen,d:\ansys\data,txt
*vwrite,noden,nx,ny,nz,bx1,by1,bz1,bs1,hx1,hy1,hz1,hs1
(f20.0,/f20.5,/f20.5,/f20.5,/e25.5,/e25.5,/e25.5,/e25.5,/e25.5,/e25.5,/e25.5,/e25.5)
*cfclos
*end
/input,ansuitmp
!****end****
!****标量法****
/prep7
emunit,mks
!空气边界
b1=10
!磁铁尺寸
b0=0.1
!1个磁铁
!三维分析
!**************
!设置单元属性
!**************
!磁铁和空气
!磁实体标量
et,1,solid96
!远场边界
et,2,infin47
!**************
!设置材料属性
!**************
!空气
mp,murx,1,1.0
!永磁铁
mp,murx,2,1.143
mp,mgxx,2,0
mp,mgyy,2,780000
mp,mgzz,2,0
!**************
!建立几何模型
!**************
block,-b1/2,b1/2,-b1/2,b1/2,-b1/2,b1/2
wpoffs,,,-b0/2
vsbw,all,
wpcsys,-1
wpoffs,,,b0/2
vsbw,all,
wpcsys,-1
wpoffs,-b0/2
wprota,,,90
vsbw,all,
wpcsys,-1
wpoffs,b0/2
wprota,,,90
vsbw,all,
wpcsys,-1
wpoffs,,-b0/2
wprota,,90
vsbw,all,
wpcsys,-1
wpoffs,,b0/2
wprota,,90
vsbw,all,
wpcsys,-1
!**************
!建立有限元模型
!**************
!磁铁
csys,0
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
aslv,r
lslar,r
lesize,all,,,4
type,1
mat,2
vmesh,all
!6个小条空气
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,y,b0/2+0.01,b1/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,y,b1/2-0.01,b1/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,y,b0/2-0.01,b0/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,b0/2,b1/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,-b1/2,-b0/2,
aslv,r
lslar,r
lsel,r,loc,y,-b1/2+0.01,-b0/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,y,-b1/2+0.01,-b1/2-0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,y,-b0/2-0.01,-b0/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,y,-b1/2,-b0/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b1/2,-b0/2,
aslv,r
lslar,r
lsel,r,loc,z,-b1/2+0.01,-b0/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,z,-b1/2+0.01,-b1/2-0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,z,-b0/2-0.01,-b0/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b1/2,-b0/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,z,b0/2+0.01,b1/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,z,b0/2-0.01,b0/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,z,b1/2-0.01,b1/2-0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,r,loc,y,-b0/2,b0/2
vsel,r,loc,z,b0/2,b1/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,x,-b1/2+0.01,-b0/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,x,-b0/2-0.01,-b0/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,-b1/2,-b0/2
aslv,r
lslar,r
lsel,r,loc,x,-b1/2-0.01,-b1/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,-b1/2,-b0/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,x,b0/2+0.01,b1/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,x,b0/2-0.01,b0/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,b0/2,b1/2
aslv,r
lslar,r
lsel,r,loc,x,b1/2-0.01,b1/2+0.01
lesize,all,,,4
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,r,loc,z,-b0/2,b0/2
vsel,r,loc,x,b0/2,b1/2
type,1
mat,1
vmesh,all
!12个薄板空气
allsel
vsel,s,loc,z,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,y,-b0/2,b0/2
aslv,r
lslar,r
lsel,u,loc,z,-b0/2+0.01,b0/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,z,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,y,-b0/2,b0/2
aslv,r
lslar,r
lsel,r,loc,z,-b0/2+0.01,b0/2-0.01
lesize,all,,,4
allsel
vsel,s,loc,z,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,y,-b0/2,b0/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,u,loc,y,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
aslv,r
lslar,r
lsel,u,loc,x,-b0/2+0.01,b0/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,u,loc,y,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
aslv,r
lslar,r
lsel,r,loc,x,-b0/2+0.01,b0/2-0.01
lesize,all,,,4
allsel
vsel,s,loc,x,-b0/2,b0/2
vsel,u,loc,y,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
type,1
mat,1
vmesh,all
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
aslv,r
lslar,r
lsel,u,loc,y,-b0/2+0.01,b0/2-0.01
lesize,all,,,20
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
aslv,r
lslar,r
lsel,r,loc,y,-b0/2+0.01,b0/2-0.01
lesize,all,,,4
allsel
vsel,s,loc,y,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
type,1
mat,1
vmesh,all
!大部分空气
allsel
vsel,u,loc,y,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
aslv,r
lslar,r
lesize,all,,,20
allsel
vsel,u,loc,y,-b0/2,b0/2
vsel,u,loc,x,-b0/2,b0/2
vsel,u,loc,z,-b0/2,b0/2
type,1
mat,1
vmesh,all
!远场边界
allsel
asel,s,loc,x,-b1/2-0.01,-b1/2+0.01
asel,a,loc,x,b1/2-0.01,b1/2+0.01
asel,a,loc,y,-b1/2-0.01,-b1/2+0.01
asel,a,loc,y,b1/2-0.01,b1/2+0.01
asel,a,loc,z,-b1/2-0.01,-b1/2+0.01
asel,a,loc,z,b1/2-0.01,b1/2+0.01
type,2
amesh,all
!**************
!求解
!**************
/solu
csys,0
allsel
!静力分析
antype,static
!施加边界条件
!标量法
!平行边界自动满足
!施加垂直边界
allsel
nsel,s,loc,x,-b1/2-0.01,-b1/2+0.01
nsel,a,loc,x,b1/2-0.01,b1/2+0.01
nsel,a,loc,y,-b1/2-0.01,-b1/2+0.01
nsel,a,loc,y,b1/2-0.01,b1/2+0.01
nsel,a,loc,z,-b1/2-0.01,-b1/2+0.01
nsel,a,loc,z,b1/2-0.01,b1/2+0.01
d,all,mag,0
!选择求解器
allsel
!简化标量法
magsolv,2
!**************
!后处理
!**************
/post1
!提取节点编号
csys,0
rsys,0
nx=3
ny=3
nz=3
noden=node(nx,ny,nz)
!提取磁感应强度
*get,bx1,node,noden,b,x
*get,by1,node,noden,b,y
*get,bz1,node,noden,b,z
*get,bs1,node,noden,b,sum
!提取磁场强度
*get,hx1,node,noden,h,x
*get,hy1,node,noden,h,y
*get,hz1,node,noden,h,z
*get,hs1,node,noden,h,sum
!写入到文件
*create,ansuitmp
*cfopen,d:\ansys\data,txt
*vwrite,noden,nx,ny,nz,bx1,by1,bz1,bs1,hx1,hy1,hz1,hs1
(f20.0,/f20.5,/f20.5,/f20.5,/e25.5,/e25.5,/e25.5,/e25.5,/e25.5,/e25.5,/e25.5,/e25.5)
*cfclos
*end
/input,ansuitmp
!****end**** |