一个焊接模拟高斯热源的实例

一个实例（摘自ANSYS百思论坛）原型为高斯热源在两个相互接触的平板上移动希望对大家有点用计算移动热源的温度分布，代码如下/CLEAR,NOSTART/FILNAM,2211/TITLE,RESIDUALSTRESSANALYSIS/UNITS,SI/CONFIG,NRES,10000/PREP7ET,1,SOLID70MP,DENS,1,5250MP,C,1,650MP,KXX,1,5/VIEW,1,1,1,1BLOCK,0,0.04,0.02,0.04,0,0.004block,0,0.04,0,0.02,0,0.004esize,0.001vsel,allvmesh,allallsel!建立接触单元ET,2,170ET,3,174KEYOPT,3,9,0R,3,RMORE,RMORE,,150RMORE,0KEYOPT,3,1,2!Generatethetargetsurfacevsel,s,volu,,1,1aslv,sASEL,S,loc,y,0.02,0.02,0CM,_TARGET,AREATYPE,2NSLA,S,1ESLN,S,0ESURF,ALLallsel,all!Generatethecontactsurfacevsel,s,volu,,2,2aslv,sASEL,S,loc,y,0.02,0.02,0CM,_CONTACT,AREATYPE,3NSLA,S,1ESLN,S,0ESURF,ALLRMORE,,150allsel,all!接触完毕FINISH/SOLULSIZE=0.001V=0.4!速度TINC=LSIZE/V!单步时间长度KK=1!能量集中系数Qmax=1800000R=0.0005MAX_TIME=41MAX_X=41MAX_Y=41*DIM,FLUX2,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME*DO,K,1,MAX_TIME,1*DO,I,1,MAX_X,1FLUX2(I,0,K)=(I-1)*LSIZE*ENDDO*DO,J,1,MAX_Y,1FLUX2(0,J,K)=(J-1)*LSIZE*ENDDO*DO,I,1,MAX_X,1*DO,J,1,MAX_Y,1XCENTER=V*(K-1)*TINCYCENTER=0.02DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)*IF,DISTANCE,LE,R,THENFLUX2(I,J,K)=Qmax/EXP(KK*DISTANCE**2/(3*r**2))*ELSEFLUX2(I,J,K)=0*ENDIF*ENDDO*ENDDOFLUX2(0,0,K)=(K-1)*TINC*ENDDO!开始计算ANTYPE,TRANSTUNIF,20!初始温度20度KBC,1*SET,DELT,1TIMINT,ONAUTOTS,OFFTT=0.1*DO,TM,1E-6,TT,TINC/VIEW,1,1,1,1/ANG,1/REP,FASTFLST,2,1,5,ORDE,1FITEM,2,2/GO!*TIME,TMAUTOTS,0DELTIM,,,,1KBC,0!*TSRES,ERASE!FLST,2,1,5,ORDE,1!FITEM,2,2/GOasel,s,loc,z,0.004,0.004SFA,all,1,HFLUX,%FLUX2%/STATUS,SOLUSOLVEallsel,allOUTRES,ALL,ALL,/PSF,HFLUX,,2/REPLOT*ENDDOOUTRES,ALL,ALL,FINISH