- 积分
- 0
- 注册时间
- 2004-10-31
- 仿真币
-
- 最后登录
- 1970-1-1
|
/CLEAR,NOSTART
/FILNAM,sample
/TITLE,RESIDUAL STRESS ANALYSIS
/UNITS,SI
/CONFIG,NRES,10000
/PREP7
ET,1,SOLID70
MP,DENS,1,7850
MP,C,1,438
MP,KXX,1,53
/VIEW,1,1,1,1
BLOCK,0,0.04,0,0.04,0,0.004
SAVE
LESIZE,1,,,80
LESIZE,4,,,80
LESIZE,9,,,4
MSHKEY,1
VATT,1,1,1
VMESH,1
/PNUM,DEFA
EPLOT
FINISH
/SOLU
LSIZE=0.0005
V=0.04 !速度
TINC=LSIZE/V !单步时间长度
KK=1 !能量集中系数
P=1400
R=0.002
r=0.01
MAX_TIME=41
MAX_X=41
MAX_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,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.02+(V*(K-1)*TINC)
DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)
*IF,DISTANCE,LE,R,THEN
FLUX2(I,J,K)=(2*P*0.4)/(3.14*R*R)/EXP(2*kk*DISTANCE**2/(0.002**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
*SET,DELT,1
TIMINT,ON
AUTOTS,OFF
TT=0.5
*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
*ENDDO
OUTRES,ALL,ALL,
FINISH
上面是激光沿斜线扫描的温度场,算出的温度有问题,只有50-60度。
在同样参数下,直线扫描温度场很正常,最高温度为900度左右。 |
|