vdload子程序，可模拟列车移动荷载

      SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,
     1 COORDS,JLTYP,SNAME)
C
      INCLUDE 'ABA_PARAM.INC'
C
      DIMENSION TIME(2),COORDS(3)
      CHARACTER*80 SNAME
C     ****** standard fatigue car JTG D64-2015 modle III ******
C     ****** vel = 1m/s EQ2 36km/h ******
      VEL=1
      PLENGTH=0.2
      ZSTART=0
C     ****** ZF = Z front & ZB = Z BEHIND ******
      ZF1=ZSTART+VEL*TIME(1)
      ZB1=ZF1-PLENGTH
      ZF2=ZF1-1.200
      ZB2=ZB1-1.200
      ZF3=ZF2-6.000
      ZB3=ZB2-6.000
      ZF4=ZF3-1.200
      ZB4=ZB3-1.200
C    
      F=0
C     ****** unit of F is N/m^2 or Pa ******
C    
    IF(COORDS(1).LE.ZF1.AND.COORDS(1).GE.ZB1)THEN
      F=1
    END IF
C    
    IF(COORDS(1).LE.ZF2.AND.COORDS(1).GE.ZB2)THEN
      F=1
    END IF
C    
    IF(COORDS(1).LE.ZF3.AND.COORDS(1).GE.ZB3)THEN
      F=1
    END IF
C    
    IF(COORDS(1).LE.ZF4.AND.COORDS(1).GE.ZB4)THEN
      F=1
    END IF
    
    RETURN
    END