这是我做的移动温度场焊接加热阶段的命令流,请各位多多指教!
/CLEAR,NOSTART
/TITLE,TemFieldANALYSIS
/FILNAM,TemField
/UNITS,SI
/CONFIG,NRES,10000
/PREP7
ET,1,SOLID70
TOFFST,273
mptemp,1,20,250,500,750,1000,1500
mptemp,7,1700,2500
mpdata,c,1,1,460,480,530,675,670,660
mpdata,c,1,7,780,820
mpdata,kxx,1,1,50,47,40,27,30,35
mpdata,kxx,1,7,140,142
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350
mpdata,dens,1,7,7300,7090
mpdata,hf,1,1,100,350,520,1000,1500,3000
mpdata,hf,1,7,3100,3500
/VIEW,1,1,1,1
BLOCK,0,0.300,0,0.010,0,0.060
SAVE
LESIZE,1,,,5
LESIZE,4,,,150
LESIZE,9,,,30
MSHKEY,1
VATT,1,1,1
VMESH,1
/PNUM,DEFA
EPLOT
FINISH
/SOLU
LSIZE=0.002
V=0.010
TINC=LSIZE/V
KK=3
Qmax=156668147
R=0.004
MAX_TIME=151
MAX_X=151
MAX_Y=6
*DIM,FLUX2,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME
*DO,K,1,MAX_TIME,1
*DO,I,1,MAX_X,1
FLUX2(I,0,K)=(I-1)*LSIZE
*ENDDO
*DO,J,1,MAX_Y,1
FLUX2(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.005
DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)
*IF,DISTANCE,LE,R,THEN
FLUX2(I,J,K)=Qmax/EXP(KK*DISTANCE**2/(0.004**2))
*ELSE
FLUX2(I,J,K)=0
*ENDIF
*ENDDO
*ENDDO
FLUX2(0,0,K)=(K-1)*TINC
*ENDDO
!开始计算
ANTYPE,TRANS
TUNIF,20 !初始温度20度
KBC,1
NEQIT,100
*SET,DELT,1
TIMINT,ON
SOLCONTROL,ON
AUTOTS,ON
TT=30
*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,%FLUX2%
/STATUS,SOLU
SOLVE
OUTRES,ALL,ALL,
/PSF,HFLUX,,2
/REPLOT
ALLSEL,ALL
SFDELE,ALL,HFLUX
*ENDDO
OUTRES,ALL,ALL, |