588
16 Computer Programs in Vibration Analysis
Table 16.6 (continued)
DK=DIF(I,I)
DO 901 K=J,L
DIF(I,K)=DIF(I,I)
901 DIF(K,I)=DIF(I,I)
IF(J-L)903,904,904
903 I=I+1
GO TO 902
904 DIF(J,J)=DK+(2.6*D)/(E*ASEC(J+1))
C 904 DIF(J,J)=DK+(2.6*D)/(ASEC(J+1))
DO 421 I=1,L
DO 421 J=1,L
421 AC(I,J)=AC(I,J)+DIF(I,J)
410 DO 407 I=1,N
ADDC(I)=ADDM(I)*YZ(MODE)
*
AM(I)=(AMM(I)+ADDM(I))/G
AM(I)=(AMM(I)+ADDC(I))/G
407 CONTINUE
BM=AM(1)
S=0.0
BI=0.0
DO 30 I=2,N
BM=BM+AM(I)
S=S+AM(I)*XX(I)
BI=BI+AM(I)*XX(I)**2.0
K=I-1
X(K)=XX(I)
AM(K)=AM(I)
30 CONTINUE
N=N-1
P=1.0/(BM*BI-S**2.0)
DO 31 I=1,N
A1(I)=0.0
B1(I)=0.0
DO 32 J=1,N
W=AM(J)*AC(J,I)
A1(I)=A1(I)+W
B1(I)=B1(I)+W*X(J)
32 CONTINUE
AA(I)=P*(S*A1(I)*AM(I)-BM*B1(I)*AM(I))
BB(I)=P*(S*B1(I)*AM(I)-BI*A1(I)*AM(I))
31 CONTINUE
DO 101 I=1,N
DO 101 J=1,N
A(I,J)=AC(I,J)*AM(J)
101 CONTINUE
DO 34 I=1,N
DO 34 J=1,N
A(I,J)=A(I,J)+X(I)*AA(J)+BB(J)
34 CONTINUE
CALL GENSS(N,1)
DO 805 I=1,N
805 A(I,I)=ABS(A(I,I))
FREQ=1.0/SQRT(A(MODE,MODE))
Z=0.0
DO 906 I=1,N
906 Z=Z+BB(I)*T(I,MODE)
(continued)