FINISH
/SOL
*SET,LSIZE,0.0002 !
*SET,V,V1/60000 !焊接速度,单位m/s
*SET,TINC,LSIZE/V !
*SET,rb,0.00045 !
*SET,pi,3.1415926
*SET,Qmax,0.25*Q/(pi*rb**2) !
!*************
*SET,MAX_TIME,51 !激光是沿着X方向移动的,每格一步,所以时间步总共51
*SET,MAX_X,51 !X方向上一共51格,所以X最大值为51
*SET,MAX_Y,81 !Y方向上一共17格,所以Y最大值为17
*DIM,arc,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME !
*DO,K,1,MAX_TIME,1 !K:循环控制变量 1:循环变量的初始值,MAX ,TIME:结束值 1
*DO,I,1,MAX_X,1
*SET,arc(I,0,K),(I-1)*LSIZE
*ENDDO
*DO,J,1,MAX_Y,1
*SET,arc(0,J,K),(J-1)*LSIZE
*ENDDO
*DO,I,1,MAX_X,1
*DO,J,1,MAX_Y,1
*SET,XCENTER,V*(K-1)*TINC
*SET,YCENTER,0
*SET,DISTANCE,SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)
*IF,DISTANCE,LE,rb,THEN
*SET,arc(I,J,K),Qmax
*ELSE
*SET,arc(I,J,K),0
*ENDIF
*ENDDO
*ENDDO
*SET,arc(0,0,K),(K-1)*TINC
*ENDDO
!************
ANTYPE,TRANS
TUNIF,20
TRNOPT,FULL !
KBC,0 !指定载荷是递增载荷或阶跃载荷0 递增 1
TIMINT,ON !
TSRES,ERASE
TINTP,0.005, , ,1, , , !采用后Euler
SOLCONTROL,ON,1, !
LNSRCH,0 !
AUTOTS,ON !
NROPT,FULL,,ON !完全的N-R
*SET,TT,50*LSIZE/V+0.0001 !
*DO,TM,1E-6,TT,TINC
TIME,TM !
DELTIM, , , ,1 !
FLST,2,1,5,ORDE,1
FITEM,2,1
/GO
!*
SFA,P51X,1,HFLUX, %ARC% !在所选择的面上施加表面载荷APPLY THERMAL
SOLVE
OUTRES,ALL,ALL !
*ENDDO
FINISH
SAVE
这是热源加载部分………… |