编辑代码

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