U型磁轭仿真，求大神看看
浏览：97996 回答：3
ANSYS新手，想对工件的磁感应强度进行分析，但是结果有问题，求大神改进。finish
/clear
keyw,magedg,1
/prep7
ET,1,SOLID97
ET,2,SOLID97
ET,3,SOLID97,2 !the coil
ET,4,SOLID97
mp,rsvx,1,3.00e-8 !the iron
mp,murx,1,1000
mp,murx,2,1 !the air
mp,murx,3,1 !the coil
mp,rsvx,3,1e-8
mp,murx,4,200
mp,rsvx,4,2e-7
!
!define the paremeters
!
iw=0.03 !the width of the iron frame
ih=0.12 !the height of the iron frame
id=0.15 !the depth of the iron frame
ix=0.007 !the coil
il=0.5 !the length of the base
ig=0.05 !the height of the base
im=0.6 !the length of the air
gap=0.003
!plot the frame
/pnum,volu,1
block,0,iw,0,ih,0,iw
block,iw+id,iw+id+iw,0,ih,0,iw
block,0,iw,ih,ih+iw,0,iw
block,iw,iw+id,ih,ih+iw,0,iw
block,iw+id,iw+id+iw,ih,ih+iw,0,iw
vglue,all
numcmp,all
vsel,s,volu,,1,5
cm,ironmodel,volu
!plot the coil
block,iw+gap,iw+id-gap,ih,ih+iw,-gap,-ix-gap
block,iw+gap,iw+id-gap,ih+iw+gap,ih+iw+ix+gap,0,iw
block,iw+gap,iw+id-gap,ih,ih+iw,iw+gap,iw+ix+gap
block,iw+gap,iw+id-gap,ih-gap,ih-ix-gap,0,iw
local,12,0,iw+gap,ih+iw,iw,,,90
wpcsys,-1,12
cyl4,0,0,gap,90,ix+gap,180,id-2*gap
cyl4,iw,0,gap,0,ix+gap,90,id-2*gap
cyl4,iw,-iw,gap,-90,ix+gap,0,id-2*gap
cyl4,0,-iw,gap,180,ix+gap,270,id-2*gap
vsel,s,volu,,6,13
vglue,all
numcmp,all
vsel,s,volu,,6,13
cm,coilmodel,volu
allsel,all
!plot the base
wpcsys,-1,0
block,-il/2+id/2+iw,id/2+il/2+iw,-gap,-ig-gap,-il/2+iw/2,il/2+iw/2
vsel,s,volu,,14
cm,basemodel,volu
!plot the air
block,-im/2+id/2+iw,id/2+im/2+iw,-ig-0.01,ih+iw+ix+0.08,-im/2+iw/2,im/2+iw/2
allsel,all
vplot,all
vovlap,all
numcmp,all
!mesh
r,1,1.008e-3,100,1.390556e-4,0,1,0
RMORE,1,0.8
LOCAL,11,1,iw+id/2,ih+iw/2,iw/2
cmsel,s,ironmodel
vatt,1,,1
cmsel,s,coilmodel
vatt,3,1,3,11
vsel,s,volu,,15
vatt,2,,2
cmsel,s,basemodel
vatt,4,,4
cmsel,s,ironmodel
cmsel,a,coilmodel
cmsel,a,basemodel
vmesh,all
vsel,s,volu,,15
esize,0.08
mshape,1,3d
mshkey,0
vmesh,all
allsel,all
!自由度
ESEL,S,TYPE,,3
nsle,s
CP,21,curr,ALL !激励线圈自由度约束
!边界条件
allsel,all
ASEL,S,,,75,76
NSLA,S,1
D,ALL,AX,0
D,ALL,AZ,0
ASEL,S,,,73,74
NSLA,S,1
D,ALL,AX,0
D,ALL,AY,0
ASEL,S,,,77,78
NSLA,S,1
D,ALL,AY,0
D,ALL,AZ,0
/solu
antyp,harm
esel,s,mat,,3
bfe,all,vltg,,20
esel,all
harfrq,10
allsel,all
eqslv,sparse
solve
/clear
keyw,magedg,1
/prep7
ET,1,SOLID97
ET,2,SOLID97
ET,3,SOLID97,2 !the coil
ET,4,SOLID97
mp,rsvx,1,3.00e-8 !the iron
mp,murx,1,1000
mp,murx,2,1 !the air
mp,murx,3,1 !the coil
mp,rsvx,3,1e-8
mp,murx,4,200
mp,rsvx,4,2e-7
!
!define the paremeters
!
iw=0.03 !the width of the iron frame
ih=0.12 !the height of the iron frame
id=0.15 !the depth of the iron frame
ix=0.007 !the coil
il=0.5 !the length of the base
ig=0.05 !the height of the base
im=0.6 !the length of the air
gap=0.003
!plot the frame
/pnum,volu,1
block,0,iw,0,ih,0,iw
block,iw+id,iw+id+iw,0,ih,0,iw
block,0,iw,ih,ih+iw,0,iw
block,iw,iw+id,ih,ih+iw,0,iw
block,iw+id,iw+id+iw,ih,ih+iw,0,iw
vglue,all
numcmp,all
vsel,s,volu,,1,5
cm,ironmodel,volu
!plot the coil
block,iw+gap,iw+id-gap,ih,ih+iw,-gap,-ix-gap
block,iw+gap,iw+id-gap,ih+iw+gap,ih+iw+ix+gap,0,iw
block,iw+gap,iw+id-gap,ih,ih+iw,iw+gap,iw+ix+gap
block,iw+gap,iw+id-gap,ih-gap,ih-ix-gap,0,iw
local,12,0,iw+gap,ih+iw,iw,,,90
wpcsys,-1,12
cyl4,0,0,gap,90,ix+gap,180,id-2*gap
cyl4,iw,0,gap,0,ix+gap,90,id-2*gap
cyl4,iw,-iw,gap,-90,ix+gap,0,id-2*gap
cyl4,0,-iw,gap,180,ix+gap,270,id-2*gap
vsel,s,volu,,6,13
vglue,all
numcmp,all
vsel,s,volu,,6,13
cm,coilmodel,volu
allsel,all
!plot the base
wpcsys,-1,0
block,-il/2+id/2+iw,id/2+il/2+iw,-gap,-ig-gap,-il/2+iw/2,il/2+iw/2
vsel,s,volu,,14
cm,basemodel,volu
!plot the air
block,-im/2+id/2+iw,id/2+im/2+iw,-ig-0.01,ih+iw+ix+0.08,-im/2+iw/2,im/2+iw/2
allsel,all
vplot,all
vovlap,all
numcmp,all
!mesh
r,1,1.008e-3,100,1.390556e-4,0,1,0
RMORE,1,0.8
LOCAL,11,1,iw+id/2,ih+iw/2,iw/2
cmsel,s,ironmodel
vatt,1,,1
cmsel,s,coilmodel
vatt,3,1,3,11
vsel,s,volu,,15
vatt,2,,2
cmsel,s,basemodel
vatt,4,,4
cmsel,s,ironmodel
cmsel,a,coilmodel
cmsel,a,basemodel
vmesh,all
vsel,s,volu,,15
esize,0.08
mshape,1,3d
mshkey,0
vmesh,all
allsel,all
!自由度
ESEL,S,TYPE,,3
nsle,s
CP,21,curr,ALL !激励线圈自由度约束
!边界条件
allsel,all
ASEL,S,,,75,76
NSLA,S,1
D,ALL,AX,0
D,ALL,AZ,0
ASEL,S,,,73,74
NSLA,S,1
D,ALL,AX,0
D,ALL,AY,0
ASEL,S,,,77,78
NSLA,S,1
D,ALL,AY,0
D,ALL,AZ,0
/solu
antyp,harm
esel,s,mat,,3
bfe,all,vltg,,20
esel,all
harfrq,10
allsel,all
eqslv,sparse
solve
yyyywin 
