关于移动热源放电通道变化的问题

ttuanjui75

以下是一个移动高斯热源的例子，命令的移动热源放电通道半径是一个确定的常量R=2.4e-5。但实际上，这个半径不是阶跃性的达到的，而是随时间扩张达到的。这个扩张时间是2us，变化权且先当作是线性的。。如果想实现这个，在命令中该怎样体现呢？
请高手具体指教！！
/UNITS,SI
/CONFIG,NRES,10000
/title,Velocity=10m/s
/PREP7
ET,1,SOLID70
MPTEMP   
MPTEMP,1,20,100,200,300,400,500
MPTEMP,7,600,800,1000,1200
MPDATA,c,1,1,245,260,264,268,272,281
MPDATA,c,1,7,289,297,306, 310
MPDATA,Kxx,1,1,144.6,139,135,131,127,123  
MPDATA,Kxx,1,7,116,109,103,101
mp,dens,1,10220
/VIEW,1,1,1,1
BLOCK,0,0.0010,0,0.0001,0,0.0001        
SAVE

LESIZE,1,,,20
LESIZE,4,,,100
LESIZE,9,,,20
MSHKEY,1
VATT,1,1,1
VMESH,1
/PNUM,DEFA
EPLOT
FINISH

LSIZE=5e-6                              
V=5
TINC=LSIZE/V
pi=3.1415926
R=2.4e-5
Q=105
Qmax=Q/(pi*R*R)
MAX_TIME=101                              
MAX_X=101
MAX_Y=21
*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=5e-5
DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)
*IF,DISTANCE,LE,R,THEN
laser(I,J,K)=Qmax/EXP(3*DISTANCE**2/(R**2))
*ELSE
laser(I,J,K)=0
*ENDIF
*ENDDO
*ENDDO
laser(0,0,K)=(K-1)*TINC
*ENDDO                                    

/SOLU
ANTYPE,TRANS
TUNIF,25
KBC,1
*SET,DELT,1                              
TIMINT,ON
AUTOTS,OFF
TT=60e-6                                   
*DO,TM,1E-6,TT,TINC                        
/VIEW,1,1,1,1
/ANG,1                                       
/REP,FAST
FLST,2,1,5,ORDE,1                           
FITEM,2,2
/GO
!*
TIME,TM                                    
AUTOTS,0
DELTIM, , , ,1
KBC,0
!*
TSRES,ERASE
FLST,2,1,5,ORDE,1                           
FITEM,2,2
/GO
!*
!*
SFA,P51X,1,HFLUX,%laser%                     
SOLVE
OUTRES,ALL,ALL,
/PSF,HFLUX,,2                                   
/REPLOT
*ENDDO
OUTRES,ALL,ALL,
FINISH

ttuanjui75

哪位大侠能关注一下啊！！期待中……

ttuanjui75

没有人关注吗？

守护蓝天

我来看看怎么样！不过我也不是很精通！

ttuanjui75

先谢过哈

jx04319

先看看 不过我也不知道呀