C REDUCE RIGHT HAND SIDE LOAD VECTOR 150 DO 180 N=1,NN KL=NDS(N)+1 KU=NDS(N+1)-1 IF(KU-KL) 180,160,160 160 K=N C=0.0 DO 170 KK=KL,KU K=K-1 170 C=C+SK(KK)*P(K) P(N)=P(N)-C 180 CONTINUE C BACK SUBSTITUTION DO 200 N=1,NN K=NDS(N) 200 P(N)=P(N)/SK(K) IF(NN.EQ.1) RETURN N=NN DO 230 L=2,NN KL=NDS(N)+1 KU=NDS(N+1)-1 IF(KU-KL) 230,210,210 210 K=N DO 220 KK=KL,KU K=K-1 220 P(K)=P(K)-SK(KK)*P(N) 230 N=N-1 RETURN 222 FORMAT(//20X,'STOP-STIFFNESS MATRIX NOT POSITIVE DEFINITE', *'NON POSITIVE PIVOT FOR EQUATION',I4,//10X,'PIVOT=',E20.12) END
SUBROUTINE PASLOD(P,NDF,NC,NSN,NEQ)
IMPLICIT REAL*8(A-H,O-Z)
COMMON/TAPES/ISTRES,NDARAY,IPR
DIMENSION P(NEQ),NDF(6,NSN)
DIMENSION CNL(6)
IF(IPR.EQ.0) WRITE (6,30)
DO 20 J=1,NC
READ(5,11) NODE,(CNL(I),I=1,6)
IF (IPR.EQ.0) WRITE (6,40) NODE,(CNL(I),I=1,6)
DO 20 I=1,6
II=NDF(I,NODE)
IF(II) 20,20,10
10 P(II)=P(II)+CNL(I) 20 CONTINUE RETURN 11 FORMAT(I10,6F10.4) 30 FORMAT(10X,'NODE',10X,'X-AXIS',10X,'Y-AXIS',10X,'Z-AXIS', *10X,'X-AXIS',10X,'Y-AXIS',10X,'Z-AXIS'/25X,'FORCE',10X,'FORCE', *10X,'FORCE',10X,'MOMENT',10X,'MOMENT',10X,'MOMENT'/) 40 FORMAT(10X,I4,6F16.3) END
SUBROUTINE DISP(D,NDF,NSN,NEQ)
IMPLICIT REAL*8(A-H,O-Z)
DIMENSION D(NEQ),NDF(6,NSN),DISPV(6)
DO 30 J=1,NSN
DO 10 I=1,6
10 DISPV(I)=0.0 DO 20 I=1,6 KK=NDF(I,J) 20 IF(KK.NE.0) DISPV(I)=D(KK) 30 WRITE(6,22) J,(DISPV(I),I=1,6) 22 FORMAT(5X,I3,4X,E19.7,2E19.7,1X,E20.6,2E18.6/) RETURN END C INTEFRACE ROUTINE FROM DEC TOPC FOR GETTING TIME FUNCTION SECNDS(X) C RETURN CURRENT TIME MIDNIGHT- X REAL*4 SECNDS,X
INTEGER*2 IHOUR, IMINUT, ISECON, IHUND
REAL*4 RHOUR,RMINUT, RSECON, RHUND
REAL*4 X1
C CALL GETTIM(IHOUR, IMINUT , ISECON, IHUND) RHOUR = FLOAT( IHOUR ) RMINUT =FLOAT ( IMINUT ) RSECON = FLOAT ( ISECON) RHUND = FLOAT( IHUND ) X1 = RHOUR3600.0+ RMINUT60.0 + RSECON + RHUND/100.0 SECNDS = X1-X RETURN END