/units,si
/prep7
et,1,70
mp,dens,1,1190
mp,kxx,1,0.2
mp,c,1,1466
block,0,0.02,0,0.03,0,0.002
block,0,0.02,0,0.03,0,-0.002
vglue,all
numcmp,all
allsel,all
lsel,s,,,1,3,2
lsel,a,,,6,8,2
lsel,a,,,13,15,2
lesize,all,0.0005,,,-0.1
lsel,inve
lesize,all,0.0005
mshkey,1
mshape,0,3d
allsel
vmesh,all
*SET,lsize,5e-4
*SET,Q,13
*SET,v,5e-3
*SET,tinc,lsize/v
*SET,pi,3.14159
*SET,rb,5e-4
*SET,Qmax,Q/(pi*rb*rb*10)
*SET,R,5E-4
*SET,MAX_TIME,41
*SET,MAX_X,41
*SET,MAX_Y,63
*DIM,laser,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME
*DO,K,1,MAX_TIME,1
*DO,I,1,MAX_X,1
laser(I,0,K)=(I-1)*LSIZE
*ENDDO
*DO,J,1,MAX_Y,1
laser(0,J,K)=(J-1)*LSIZE
*ENDDO
*DO,I,1,MAX_X,1
*DO,J,1,MAX_Y,1
XCENTER=V*(K-1)*TINC
YCENTER=0.015
DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)
*IF,DISTANCE,LE,R,THEN
laser(I,J,K)=Qmax/EXP(2*DISTANCE**2/(rb**2))
*ELSE
laser(I,J,K)=0
*ENDIF
*ENDDO
*ENDDO
laser(0,0,K)=(K-1)*TINC
*ENDDO
/solu
antype,trans
tunif,293
TIMINT,ON
*SET,TT,0.02/v
*DO,TM,1E-6,TT,TINC
TIME,TM
KBC,1
AUTOTS,on
Nropt,auto
tsres,erase
asel,s,,,1
sfa,all,,hflux,%laser%
nsel,s,loc,y,0
sf,all,conv,100,293
allsel,all
nsel,s,loc,y,0.03
sf,all,conv,100,293
allsel,all
nsel,s,loc,x,0
sf,all,conv,100,293
allsel,all
nsel,s,loc,x,0.02
sf,all,conv,100,293
allsel,all
nsel,s,loc,z,0.002
sf,all,conv,100,293
allsel,all
nsel,s,loc,z,-0.002
sf,all,conv,100,293
allsel,all
outres,all,all
solve
outres,all,all
*ENDDO
/solu
antype,,rest
sfadele,all,,hflux 去掉热流
time,10
solve
可是我的温度居然是升高的。。。。