16.5 Computer Program for the Free Vibration Analysis of Ships …
589
Table 16.6 (continued)
Z=Z*FREQ**2.0
FREQ=9.545*FREQ
WRITE(12,420)MODE,FREQ
420
FORMAT(/'NATURAL FREQUENCY OF MODE ',i4,' IS : ',G12.4)
430
FORMAT(/'
MODE SHAPES ',/(5G12.4))
WRITE(12,430)(T(I,MODE),I=1,N)
N=N+1
IF(MODE-MOD)417,419,419
417 MODE=MODE+1
GO TO 410
419 END
*******************************************************************
SUBROUTINE INCOF(AI,X,N,I,S,B)
*
Subroutine to determine the influence coefficients
DIMENSION B(50),X(50),C(50),CG(50),AI(50)
M=N-1
DO 30 J=1,N
B(J)=0.
C(J)=0.
CG(J)=0.
30 CONTINUE
L=I-1
DO 31 J=1,L
Y1=(X(I)-X(J))/AI(J)
Y2=(X(I)-X(J+1))/AI(J+1)
C(J)=C(J)+(Y1+Y2)*S/2.0
CG(J)=CG(J)+((2.0*Y1+Y2)*S)/(3.0*(Y1+Y2))
31 CONTINUE
DO 22 J=1,L
DO 22 K=1,M
IF(J-K)5,6,22
6 B(K)=B(K)+C(J)*CG(J)
GO TO 22
5 B(K)=B(K)+C(J)*(CG(J)+X(K+1)-X(J+1))
22 CONTINUE
RETURN
END
********************************************
SUBROUTINE GENSS(N,ISW)
*
Subroutine to solve eigenvalue problem by QR Transformation
Method
*
The subroutine deals with solid and unsymmetric matrices
DIMENSION AA(50,50),A(50,50),T(50,50)
COMMON A,T
IF(ISW)100,106,100
100 DO 104 I=1,N
DO 102 J=1,N
102 T(I,J)=0.0
104 T(I,I)=1.0
C
SET UP CONSTANTS, CALCULATE EN
106 YR=0.1E-07
EPS=50.0*YR
EP=YR*1.0E-02
MARK=0
NM1=N-1
C
ITM LOOP
(continued)