3d静态磁场分析（标量法）

yu2008fish

在《ansys电磁场分析指南》一书中，第五章3d静态磁场分析（标量法）里用gui方式和apdl方式做同一个例子，我照书上做了，按照书上的apdl我得出来的结果却和gui得出的不一样，主要是使用apdl文件在执行过程中出现了error提醒，最后计算出的电感值也差别比较大。有没有大虾能够帮我看一下，到底是怎么回事呢？
我把gui做法的log文件和apdl文件都传上来了。
希望能有人帮我弄明白。

yu2008fish

/TITLE,3d static force problem - tetrahedral

/COM,   
/COM,Preferences for GUI filtering have been set to display:
/COM,  Magnetic-Nodal   

*SET,n,500  
*SET,i,6
/PREP7  

ET,1,SOLID96

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MURX,1,,1
tb,bh,2,1,30

tbpt,,355,0.7   
tbpt,,405,0.8   
tbpt,,470,0.9   
tbpt,,555,1
tbpt,,673,1.1   
tbpt,,836,1.2   
tbpt,,1065,1.3  
tbpt,,1220,1.35
tbpt,,1420,1.4  
tbpt,,1720,1.45
tbpt,,2130,1.5  
tbpt,,2670,1.55
tbpt,,3480,1.6  
tbpt,,4500,1.65
tbpt,,5950,1.7  
tbpt,,7650,1.75
tbpt,,10100,1.8
tbpt,,13000,1.85
tbpt,,15900,1.9
tbpt,,21100,1.95
tbpt,,26300,2   
tbpt,,32900,2.05
tbpt,,42700,2.1
tbpt,,61700,2.15
tbpt,,64300,2.2
tbpt,,110000,2.25   
tbpt,,135000,2.3
tbpt,,200000,2.41   
tbpt,,400000,2.69   
tbpt,,800000,3.22   

tbcopy,all,2,3  
!*  
SAVE
/PNUM,KP,0  
/PNUM,LINE,0
/PNUM,AREA,0
/PNUM,VOLU,1
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0   
!*  
/PNUM,ELEM,0
/REPLOT
!*  
block,0,63.5,0,12.5,0,25
/VIEW, 1 ,1,1,1
/ANG, 1
/REP,FAST   
block,38.5,63.5,0,12.5,25,125   
block,13.5,63.5,0,12.5,125,150  
vglue,all   
block,0,12.5,0,5,26.5,125   
block,0,13,0,5.5,26,125.5   
vovlap,1,2  
numcmp,volu
cyl4,,,0,0,100,90,175   
vovlap,all  
numcmp,volu
SAVE
CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       3, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,3,6,ORDE,2   
FITEM,5,3   
FITEM,5,-5  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       2, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
ALLSEL,ALL  
SMRT,6  
SMRT,7  
SMRT,8  
MSHAPE,1,3D
MSHKEY,0
!*  
FLST,5,6,6,ORDE,2   
FITEM,5,1   
FITEM,5,-6  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
/UI,MESH,OFF
/PNUM,KP,0  
/PNUM,LINE,0
/PNUM,AREA,0
/PNUM,VOLU,1
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,1   
!*  
/PNUM,MAT,1
/REPLOT
!*  
EPLOT   
ESEL,S,MAT,,3   
CM,arm,ELEM
!*  
!*  
FMAGBC,'ARM'
SAVE
ALLSEL,ALL  
FLST,2,6,6,ORDE,2   
FITEM,2,1   
FITEM,2,-6  
VLSCAL,P51X, , ,0.001,0.001,0.001, ,0,1
!*  
LOCAL,12,0,0,0,75/1000, , , ,1,1,   
WPCSYS,-1,12,   
RACE,0.0285,0.0285,0.014,n*i,0.018,0.0966, , , 'coil1'  
!*  
/SHRINK,0   
/ESHAPE,1   
/EFACET,1   
/RATIO,1,1,1
/CFORMAT,32,0   
/REPLOT
!*  
EPLOT   
FINISH  
/SOL
FLST,2,1,1,ORDE,1   
FITEM,2,2   
D,P51X,MAG,0,   
ALLSEL,ALL  
ALLSEL,ALL  
MAGSOLV,3, , ,0.001,25,1
FINISH  
/POST1  
!*  
!*  
FMAGSUM,'ARM'   
*DIM,cur,ARRAY,1,1,1, , ,   
!*  
*SET,CUR(1,1,1) , 6
FINISH  
/SOL
!*  
!*  
LMATRIX,1, 'coil' , 'CUR' , ' '
/COM,ANSYS RELEASE 10.0    UP20050718       17:06:51    05/24/2011
/COM,ANSYS RELEASE 10.0    UP20050718       17:06:58    05/24/2011

yu2008fish

finish
/clear
/prep7
/title,3-d static force problem-tetrahedral
/com,magnetic-edge

n=500
i=6

et,1,96
mp,murx,1,1
tb,bh,2,1,30
!tbtemp,0
tbpt,,355,0.7
tbpt,,405,0.8
tbpt,,470,0.9
tbpt,,555,1
tbpt,,673,1.1
tbpt,,836,1.2
tbpt,,1065,1.3
tbpt,,1220,1.35
tbpt,,1420,1.4
tbpt,,1720,1.45
tbpt,,2130,1.5
tbpt,,2670,1.55
tbpt,,3480,1.6
tbpt,,4500,1.65
tbpt,,5950,1.7
tbpt,,7650,1.75
tbpt,,10100,1.8
tbpt,,13000,1.85
tbpt,,15900,1.9
tbpt,,21100,1.95
tbpt,,26300,2
tbpt,,32900,2.05
tbpt,,42700,2.1
tbpt,,61700,2.15
tbpt,,64300,2.2
tbpt,,110000,2.25
tbpt,,135000,2.3
tbpt,,200000,2.41
tbpt,,400000,2.69
tbpt,,800000,3.22
!mpcopy,,2,3
tbcopy,all,2,3

/pnum,volu
block,0,63.5,0,12.5,0,25
/view,1,1,1,1
/replot
block,38.5,63.5,0,12.5,25,125
block,13.5,63.5,0,12.5,125,150
vglue,all

block,0,12.5,0,5,26.5,125
block,0,13,0,5.5,26,125.5
vovlap,1,2
numcmp,volu
cyl4,,,0,0,100,90,175
vovlap,all
numcmp,volu

vsel,s,volu,,1
vatt,3,1,1
vsel,s,volu,,3,5
vatt,2,1,1

allsel,all
smrt,8
mshape,1,3d
mshkey,0
vmesh,all
/pnum,mat,1
/number,1
eplot

esel,s,mat,,3
cm,arm,elem
fmagbc,'arm'

allsel,all
vlscale,all,,,0.001,0.001,0.001,,0,1

local,12,0,0,0,75/1000, , , ,1,1,   
wpcsys,-1,12,
race,0.0285,0.0285,0.014,n*i,0.018,0.0966,,,'coil1'
/eshape,1
eplot
save
finish

/solu
d,2,mag,0

allsel,all
magsolv,3,,,,,1
finish

/post1
/com
frmagsum,'arm'
finish
save

/solu
*dim,cur,array,1
cur(1)=i
lmatrix,1,'coil','cur',
finish

jindan1224

我也在看这个例题啊~~ 也遇到很多不理解的地方啊~
error有没有提醒你哪里出问题了啊？