大家帮我看看命令流有什么问题

dannyboy

做辐射分析，可得到的温度很高，与不加表面单元相比差距太大了，大家帮忙看看原因是什么？谢谢了。还想请教一下，如果加入辐射分析后，在后处理中看到的节点结果的单位是K吗？还是温度偏移前使用的单位制？
/prep7  
et,1,solid70
et,2,surf152,,,,1,1
keyopt,2,9,1
r,2,1,5.67e-8   
mp,kxx,1,157
mp,c,1,700  
mp,dens,1,3000  
mp,emis,2,0.8   
BLOCK,0,0.025,0,0.025,0,0.007   
esize,0.001,
type,1  
mat,1   
mshkey,0
mshape,1,3d
vmesh,1
type,2  
real,2  
mat,2   
n,8000,0.0125,0.0125,0.020  
esurf,8000  
*SET,aa,0.005   
*SET,bb,0.002   
*SET,cc,0.005   
*SET,u,150  
*SET,i,12   
*SET,ef1,0.75   
*SET,pi,3.1415926   
*SET,xcenter,0  
*SET,ycenter,0  
*SET,zcenter,0  
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNC_C6
*DEL,_FNC_C7
*DEL,_FNC_C8
*DEL,_FNC_C9
*DEL,_FNC_C10   
*DEL,_FNCCSYS   
*SET,_FNCNAME,'hgen2'   
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*DIM,_FNC_C6,,1
*DIM,_FNC_C7,,1
*DIM,_FNC_C8,,1
*DIM,_FNC_C9,,1
*DIM,_FNC_C10,,1
*SET,_FNC_C1(1),ef1
*SET,_FNC_C2(1),u   
*SET,_FNC_C3(1),i   
*SET,_FNC_C4(1),xcenter
*SET,_FNC_C5(1),aa  
*SET,_FNC_C6(1),ycenter
*SET,_FNC_C7(1),cc  
*SET,_FNC_C8(1),zcenter
*SET,_FNC_C9(1),bb  
*SET,_FNC_C10(1),pi
*SET,_FNCCSYS,0
! /INPUT,hgen2.func
*DIM,%_FNCNAME%,TABLE,12,35,1,,,,%_FNCCSYS%
!   
! Begin of equation: 6*sqrt(3)*ef1*u*i*exp((-3)*(((({X}-xcenter)/aa)^2+(({Y}-   
! ycenter)/cc)^2+(({Z}-zcenter)/bb)^2)))/(aa*bb*cc*pi*sqrt(pi))
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)%
*SET,%_FNCNAME%(9,0,1), %_FNC_C7(1)%
*SET,%_FNCNAME%(10,0,1), %_FNC_C8(1)%   
*SET,%_FNCNAME%(11,0,1), %_FNC_C9(1)%   
*SET,%_FNCNAME%(12,0,1), %_FNC_C10(1)%  
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 1, -2, 3, 19
*SET,%_FNCNAME%(0,7,1), 0.0, -2, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,8,1), 0.0, -3, 0, 1, 0, 0, -2
*SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, -2, 2, -3
*SET,%_FNCNAME%(0,10,1), 0.0, -2, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,11,1), 0.0, -3, 0, 1, -4, 3, -2   
*SET,%_FNCNAME%(0,12,1), 0.0, -2, 0, 1, 2, 2, 20
*SET,%_FNCNAME%(0,13,1), 0.0, -4, 0, 1, -2, 4, 21   
*SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,15,1), 0.0, -5, 0, 1, -4, 17, -2  
*SET,%_FNCNAME%(0,16,1), 0.0, -2, 0, 1, 3, 2, 22
*SET,%_FNCNAME%(0,17,1), 0.0, -4, 0, 1, -2, 4, 23   
*SET,%_FNCNAME%(0,18,1), 0.0, -2, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,19,1), 0.0, -6, 0, 1, -4, 17, -2  
*SET,%_FNCNAME%(0,20,1), 0.0, -2, 0, 1, -5, 1, -6   
*SET,%_FNCNAME%(0,21,1), 0.0, -4, 0, 1, 4, 2, 24
*SET,%_FNCNAME%(0,22,1), 0.0, -5, 0, 1, -4, 4, 25   
*SET,%_FNCNAME%(0,23,1), 0.0, -4, 0, 2, 0, 0, -5
*SET,%_FNCNAME%(0,24,1), 0.0, -6, 0, 1, -5, 17, -4  
*SET,%_FNCNAME%(0,25,1), 0.0, -4, 0, 1, -2, 1, -6   
*SET,%_FNCNAME%(0,26,1), 0.0, -2, 0, 1, -3, 3, -4   
*SET,%_FNCNAME%(0,27,1), 0.0, -2, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,28,1), 0.0, -3, 0, 1, -1, 3, -2   
*SET,%_FNCNAME%(0,29,1), 0.0, -1, 0, 1, 21, 3, 25   
*SET,%_FNCNAME%(0,30,1), 0.0, -2, 0, 1, -1, 3, 23   
*SET,%_FNCNAME%(0,31,1), 0.0, -1, 0, 1, -2, 3, 26   
*SET,%_FNCNAME%(0,32,1), 0.0, -2, 16, 1, 26, 0, 0   
*SET,%_FNCNAME%(0,33,1), 0.0, -4, 0, 1, -1, 3, -2   
*SET,%_FNCNAME%(0,34,1), 0.0, -1, 0, 1, -3, 4, -4   
*SET,%_FNCNAME%(0,35,1), 0.0, 99, 0, 1, -1, 0, 0
! End of equation: 6*sqrt(3)*ef1*u*i*exp((-3)*(((({X}-xcenter)/aa)^2+(({Y}-
! ycenter)/cc)^2+(({Z}-zcenter)/bb)^2)))/(aa*bb*cc*pi*sqrt(pi))
!-->
allsel,all  
/solu   
alls
antype,trans
toffst,273  
d,8000,temp,20  
tunif,20
autots,on   
kbc,1   
outres,,all
outpr,nsol,last
*SET,ycenter,0.005  
*SET,hgen2(8,0,1),ycenter   
LOCAL,16,2,0,0.005,0,,, ,1.25,1,
nsel,s,loc,x,0,0.005
nsel,r,loc,y,0,180  
nsel,r,loc,z,0,90   
nplot   
cm,cm3,node
bf,cm3,hgen,%hgen2%
*SET,t,0.2   
time,t  
cmdele,cm3  
allsel,all  
solve

*SET,ycenter,0.01  
*SET,hgen2(8,0,1),ycenter   
LOCAL,16,2,0,0.01,0,,, ,1.25,1,
nsel,s,loc,x,0,0.005
nsel,r,loc,y,0,180  
nsel,r,loc,z,0,90   
nplot   
cm,cm3,node
bf,cm3,hgen,%hgen2%
*SET,t,t+0.2   
time,t  
cmdele,cm3  
allsel,all  
solve

*SET,ycenter,0.015  
*SET,hgen2(8,0,1),ycenter   
LOCAL,16,2,0,0.015,0,,, ,1.25,1,
nsel,s,loc,x,0,0.005
nsel,r,loc,y,0,180  
nsel,r,loc,z,0,90   
nplot   
cm,cm3,node
bf,cm3,hgen,%hgen2%
*SET,t,t+0.2   
time,t  
cmdele,cm3  
allsel,all  
solve

*SET,ycenter,0.02  
*SET,hgen2(8,0,1),ycenter   
LOCAL,16,2,0,0.02,0,,, ,1.25,1,
nsel,s,loc,x,0,0.005
nsel,r,loc,y,0,180  
nsel,r,loc,z,0,90   
nplot   
cm,cm3,node
bf,cm3,hgen,%hgen2%
*SET,t,t+0.2   
time,t  
cmdele,cm3  
allsel,all  
solve

*SET,ycenter,0.025  
*SET,hgen2(8,0,1),ycenter   
LOCAL,16,2,0,0.025,0,,, ,1.25,1,
nsel,s,loc,x,0,0.005
nsel,r,loc,y,0,180  
nsel,r,loc,z,0,90   
nplot   
cm,cm3,node
bf,cm3,hgen,%hgen2%
*SET,t,t+0.2   
time,t  
cmdele,cm3  
allsel,all  
solve

[ 本帖最后由 dannyboy 于 2007-1-30 10:23 编辑 ]

wxpand2004

顶一下！！！