一个完整的三维移动高司分布表面热源焊接例子.txt

ansys.com上的ansys命令流,主要为热,流体,接触分析方面的内容

/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