DFLUX子程序如何控制热源移动路径啊?

czy2019

      程序如下，如果R=sqrt((x-x0-d)**2+(y-y0)**2)，则热源沿x方向移动，如果 R=sqrt((x-x0)**2+(y-y0-d)**2)，则热源沿y方向移动，但是我想根据时间控制移动路径成回字型，比如总时长0.3秒，前0.1秒沿x正方向移动，然后0.1秒沿y正方向移动，再0.1秒沿x负方向移动，不知道怎么编程？求大佬帮忙




SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,
     1 JLTYP,TEMP,PRESS,SNAME)


      INCLUDE 'ABA_PARAM.INC'

      parameter(one=1.d0)

      DIMENSION COORDS(3),FLUX(2),TIME(2)

      CHARACTER*80 SNAME
      
      Q=70
      n=1
      Rh=0.001
      PI=3.1415
      
      v=0.05
      x=COORDS(1)
      y=COORDS(2)
      z=COORDS(3)

      x0=0
      y0=-4e-3
      d=v*TIME(1)

C     JLTYP＝0，表示为面热源      
      JLTYP=0
      
      R=sqrt((x-x0-d)**2+(y-y0)**2)
      heat=3*Q*n/(PI*Rh**2)
      shape1=exp(-3*R**2/Rh**2)
      FLUX(1)=heat*shape1

      RETURN
      END

1099

楼主问题解决了吗？

aqustut

学习一下学习一下