高手帮忙看看

ecjtuljsh

麻烦大哥帮帮看看

我找了一个高斯热源加载的APDL，但是中间有些地方不是很明白，望高手教教我，谢谢！

三维移动高斯分布热源表面焊命令流

/prep7  
/UNITS,si
   
et,1,solid70
mp,c,1,544.28
mp,kxx,1,15.06
mp,dens,1,4500
  
block,0,30,0,15,0,3
  
esize,1  
vmesh,all  
FINISH  
  
/SOLU  
!PARAMETER TABLE TO INPUT HEAT FLUX AS A MOVING SOURCE  
!SPEED OF HEAT SOURCE IN X DIRECTION v=(200/60)m/s  
!GAUSS VARIATION ASSUMED IN X & Y DIRECTION AS q=qmax/EXP(K*R**2)  
!TOTAL HEAT INPUT Q=3.14*qmax/K  
!DIMENSIONING TABLE NAMED "FLUX2"  
v=1  
L=30  
NUM_TIME=L/v+1
NUM_X=L/v+1  
   
*DIM,FLUX2,TABLE,NUM_TIME,NUM_X,7,TIME,X,Y  
   
*DO,K,1,7,1  
  
*DO,I,1,NUM_TIME,1  
FLUX2(I,0,K)=(I-1)  
*ENDDO  
  
*DO,J,1,NUM_X,1  
FLUX2(0,J,K)=(J*v-v)  
*ENDDO  
  
Q=1000  
KAPA=1/4  
qmax=((Q*KAPA)/3.14159)
PERCENT=0.98
RADIUS=ABS(SQRT(1/KAPA*LOG(1/(1-PERCENT))))  
  
!N=[RADIUS/v]  
  
*DO,ANTIME,1,NUM_TIME,1  
*DO,N,0,15,1  
XDIST=(N*v)  
*IF,(ANTIME-N),GT,0,THEN  
*IF,SQRT((XDIST)**2+((RADIUS/6)*(K-1))**2),LE,RADIUS,THEN  
!  
GAUSS=qmax/EXP(KAPA*((XDIST)**2+((RADIUS/6)*(K-1))**2))  
FLUX2(ANTIME,ANTIME-N,K)=GAUSS*10e5
!  
*ELSE  
*ENDIF  
*ENDIF  
!  
*IF,(ANTIME+N),LE,NUM_TIME,THEN  
*IF,SQRT((XDIST)**2+((RADIUS/6)*(K-1))**2),LE,RADIUS,THEN  
!  
GAUSS=qmax/EXP(KAPA*((XDIST)**2+((RADIUS/6)*(K-1))**2))  
FLUX2(ANTIME,ANTIME+N,K)=GAUSS*10e5  
!  
*ELSE  
*ENDIF  
*ENDIF  
*ENDDO  
*ENDDO  
!  
!INPUT Y DIRECTION DIVISIONS  
!---------------------------  
FLUX2(0,0,K)=(K-1)*(RADIUS/6)
*ENDDO  
  
antype,trans
tunif,20
kbc,1
   
*SET,delt,1  
  timint,on  
  autots,off
  
*do,tm,1e-6,30,0.5   
/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, %FLUX2%   
/STATUS,SOLU
SOLVE
OUTRES,ALL,ALL,
/PSF,hflux,,2
  
/replot  
     
*enddo  
OUTRES,ALL,ALL,   
  
FINISH

问题（1）是为什么要在循环加载中用两次
FLST,2,1,5,ORDE,1   
FITEM,2,2   
/GO  
!*   
!*   
（2）在加载过程中为什么是AUTOTS,0 ？
       AUTOTS，0是否就是AUTOTS，OFF？
此处用ＡＵＴＯＴＳ，０是有什么原因吗？
还有DELTIM, , , ,1　中省略了几个参数是该做如何理解，该把他们理解为多少？
（３）TSRES,ERASE　是什么意思？好像用到不是很多。

ecjtuljsh

怎么没人啊