编辑代码

 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
     IF (TIME(2)<=1.0) THEN
         F = 10*(15.0 - COORDS(2))
     ELSE IF (COORDS(2)<=(15.0 - 2.5*(TIME(2)-1.0)/24.0)) THEN
         F = 10*((15.0 - 2.5*(TIME(2)-1.0)/24.0) - COORDS(2))
     END IF
     RETURN
     END
     SUBROUTINE DISP(U, KSTEP, KINC, TIME, NODE, NOEL, JDOF, COORDS)
C
     INCLUDE 'ABA_PARAM.INC'
C
     DIMENSION U(3), TIME(2), COORDS(3)
C
     IF (TIME(2)<=1.0) THEN
         U(1) = 10*(15.0 - COORDS(2))
     ELSE IF (COORDS(2)<=(15.0 - 2.5*(TIME(2)-1.0)/24.0)) THEN
         U(1) = 10*((15.0 - 2.5*(TIME(2)-1.0)/24.0) - COORDS(2))
     END IF
     RETURN
     END