3d静态磁场分析（标量法）
在《ansys电磁场分析指南》一书中，第五章3d静态磁场分析（标量法）里用gui方式和apdl方式做同一个例子，我照书上做了，按照书上的apdl我得出来的结果却和gui得出的不一样，主要是使用apdl文件在执行过程中出现了error提醒，最后计算出的电感值也差别比较大。有没有大虾能够帮我看一下，到底是怎么回事呢？我把gui做法的log文件和apdl文件都传上来了。
希望能有人帮我弄明白。
gui法的log文件
/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
apdl文件
本帖最后由 yu2008fish 于 2011-5-24 21:48 编辑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 我也在看这个例题啊~~ 也遇到很多不理解的地方啊~
error有没有提醒你哪里出问题了啊？
页:
[1]