j=1                ！定義求解時(shí)間
timint,on

*do,i,1,20

esel,s,mat,,1     ！只選擇材料1的單元
nsle,s,1     !選擇當(dāng)前單元的節(jié)點(diǎn)
nsel,r,loc,z,(zmin+zmax)/2-dz,(zmin+zmax)/2+dz   ！在已選節(jié)點(diǎn)中重新去選
nsel,r,loc,y,2.5
nsel,r,loc,x,-dx,dx       ！在大正方形中選擇出中心的正方形
nplot
sf,all,hflux,qm    ！中心高斯熱源


r=1                    !給個(gè)初值
*do,m,1,3
qr=qm*exp(-3*r*r/16)           !高斯熱源計(jì)算公式，16為等離子體云根部半徑的平方(有效加熱面)。

esel,s,mat,,1
nsle,s,1
nsel,r,loc,y,2.5
nsel,r,loc,x,-(m+1)*dx,(m+1)*dx
nsel,r,loc,z,(zmin+zmax)/2-(m+1)*dz,(zmin+zmax)/2-m*dz
nplot
sf,all,hflux,qr

esel,s,mat,,1
nsle,s,1
nsel,r,loc,y,2.5
nsel,r,loc,x,-(m+1)*dx,(m+1)*dx
nsel,r,loc,z,(zmin+zmax)/2+m*dz,(zmin+zmax)/2+(m+1)*dz
nplot
sf,all,hflux,qr

esel,s,mat,,1
nsle,s,1
nsel,r,loc,y,2.5
nsel,r,loc,z,(zmin+zmax)/2-(m+1)*dz,(zmin+zmax)/2+(m+1)*dz
nsel,r,loc,x,-(m+1)*dx,-m*dx
nplot
sf,all,hflux,qr

esel,s,mat,,1
nsle,s,1
nsel,r,loc,y,2.5
nsel,r,loc,z,(zmin+zmax)/2-(m+1)*dz,(zmin+zmax)/2+(m+1)*dz
nsel,r,loc,x,m*dx,(m+1)*dx
nplot
sf,all,hflux,qr

r=r+1

*enddo   


time,j                     
deltim,j,j,j         
alls
solve

esel,s,mat,,1
nsle,s,1               
sfdele,all,hflux   


j=j+1
zmin=zmin+zd
zmax=zmax+zd
*enddo
alls
eplot 返回小木蟲(chóng)查看更多