求助

dannyboy

高斯热源和双椭球同时加载，发现高斯热源不起作用，郁闷好几天了，请指教一下吧，多谢了。
t=0
kk=435
ang=135
r=0.008
*do,i,1,2
csys,0  
wpcsys,-1   
*get,xcoor,kp,kk,loc,x  
*get,ycoor,kp,kk,loc,y  
*get,zcoor,kp,kk,loc,z  
LOCAL,13,1,xcoor,ycoor,zcoor,ang, ,-45,1,1,
LOCAL,14,2,XCOOR,YCOOR,ZCOOR,ANG,,-45,1,1,
local,15,2,xcoor,ycoor,zcoor,ang+180,,-45,2.5,1   
allsel,all
csys,13
wpcsys,-1
asel,s,,,364
nsla,s,1
nsel,r,loc,x,0,r
*SET,xcenter,xcoor  
*SET,ycenter,ycoor  
*SET,zcenter,zcoor
*SET,myguass1(5,0,1),xcenter
*SET,myguass1(6,0,1),ycenter
sf,all,hflux,%myguass1%
*SET,hgen1(6,0,1),xcenter   
*SET,hgen1(8,0,1),ycenter   
*SET,hgen1(10,0,1),zcenter
allsel,all
csys,14
wpcsys,-1  
nsel,s,loc,x,0,0.0055   
nsel,r,loc,y,0,180
nsel,r,loc,z,-90,0
cm,cm1,node
bf,all,hgen,%hgen1%
allsel,all  
csys,15
wpcsys,-1
nsel,s,loc,x,0,0.008   
nsel,r,loc,y,0,180
nsel,r,loc,z,-90,90
cm,cm2,node
bf,all,hgen,%hgen1%  
t=t+0.65
time,t
autots,0
nsubst,1
kbc,0
allsel,all
solve
ang=ang+1.2
kk=kk+1
*enddo

加载部分的命令流，热源的设置都做完了。