forums.silverfrost.com

[shah]

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 = RHOUR*3600.0+ RMINUT*60.0 + RSECON + RHUND/100.0
SECNDS = X1-X
RETURN
END