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(KSTEP.EQ.1) THEN U(1)=45.*TIME(1) U(2)=0.0 U(3)=0.0 ELSE TMP1=1. TMP4=4. PI=ATAN(TMP1)*TMP4 TMP22=.22222222 TMP45=45. U(1)=2.*SIN(TMP22*PI*TIME(1))+TMP45 U(2)=(2.*TMP22*PI)*COS(TMP22*PI*TIME(1)) U(3)=-(2.*(TMP22*PI)**2)*SIN(TMP22*PI*TIME(1)) END IF RETURN END