论坛 › C01：ANSYS--低频软件(Maxwell \PExprt\Rmxprt\Simplorer) › 永磁铁的三维磁场分析：标量法与棱边法差别巨大？

gigolo 发表于 2011-8-30 11:05:10

永磁铁的三维磁场分析：标量法与棱边法差别巨大？

对于一个永磁铁的三维分析，分别采用标量法和棱边法计算，但差别很大，请问到底是什么原因？是因为电磁场的理论，还是我的命令流有错误，请大家讨论一下。


模型简介：
永磁铁：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****

gigolo 发表于 2011-8-30 11:12:01

另外，将永磁铁视为磁偶极子，采用磁矩来描述永磁铁，将剩磁和矫顽力转化为等效磁矩，上述模型中的磁矩为111，计算得到坐标(3,3,3)处的磁感应强度：

Bx=7.91184E-08
By=0
Bz=7.91184E-08

与有限元计算结果差别巨大。


有没有人出来讨论一下啊

ourxialing 发表于 2011-9-5 20:32:07

在理论上，当存在非均匀介质时，用基于节点的连续矢量位A来进行有限元计算会产生不精确的解，所以我感觉，用棱边法会精确一些

gigolo 发表于 2011-9-6 09:19:08

ourxialing 发表于 2011-9-5 20:32 static/image/common/back.gif
在理论上，当存在非均匀介质时，用基于节点的连续矢量位A来进行有限元计算会产生不精确的解，所以我感觉， ...

谢谢你的回复。
根据计算结果，Bx和Bz相差比较小，在10%~50%；但是标量法计算的By为正，棱边法计算的By为负，差别太大了。


另外，根据磁偶极子理论计算的磁感应强度，和有限元结果，差2个数量级。这个也很难解释。

frankygtd 发表于 2012-3-17 15:18:01

很好的例子，学习学习

panda_neiep 发表于 2012-3-21 10:23:16

正有此问！有没有详细的理论分析啊？

186hesheng 发表于 2012-6-12 15:39:18

请问楼主这个问题找到原因了吗？这是一个很好的对比标量法和单元法的例子啊

gigolo 发表于 2012-6-20 21:52:41

186hesheng 发表于 2012-6-12 15:39 static/image/common/back.gif
请问楼主这个问题找到原因了吗？这是一个很好的对比标量法和单元法的例子啊 ...

实在惭愧，土木工程专业的，电磁方面水平有限。没有搞定，放弃了有限元分析，最后只考虑了磁偶极子的理论解。

gigolo 发表于 2012-6-20 21:53:53

但是如果采用线圈模拟磁偶极子，有限元计算结果和磁偶极子的计算结果，两者比较接近，因为本身，磁偶极子将磁性物质视为环形电流。

查看完整版本: 永磁铁的三维磁场分析：标量法与棱边法差别巨大？