RUNT VERSICN FEB 7^ B 13«0* 0«/09/76 



SUBROUTINE SPRINGdC) 



5 C0HM0N/TW0/WY(2/3»,EX(2»3»20)»FX(2,3,2P),FY(2»3,20)*TENS(2»3»20)» 

 2 FEXT(9),SEP(9),TENS1(9),TENS2(9),NANGLE»NC0ND,TITLE(8) 



6 CCf*H0N/THREE/X(2)»F(2) 



6 C0HM0N/F0UR/YCG,XCAB(2).YCAB(2)»TTDE(5),ITIDF 



6 DIMENSION OF X0X(2) ,DFYOX( 2) t DFXOY( 2 ),DFYOY( 2) »0 t 3) # FXX (2# 3 )» 



2 FYX«2,3),FVt2)»DTDX(2).DT0Y(2)»0TDR(2)»T(2»3> 

 6 REAL KMll,KH12»Km3»KM21»KH22»KM23,KH3I,KM32»KM33 



C»***SUBR0UT1NE SPRING COMPUTES THE BREAKWATER SPRING CONSTANTS 



C***COMPUTE THE SPRING CONSTANTS FOR EACH CABLE 



C HORIZONTAL FORCE AT EQUIL I9RIUM«F ( I ) FOR CABLE I 



C VERT FORCE AT EOUILIBRIUM»FVt I » FOR CABLE I 



6 DO 1« I«l»2 



7 00 12 J«l,3 



10 CALL LTERPS (I, J,NAN6LE» EX, FX» X ( I) , FXX < I, J ) ,DF ) 



27 IF(J .EQ. 1) DFXDX(I)— OF 



35 CALL LTERPS ( I» J^NANGLE. EX, TENS,X( I ),T( I, J ) »0T) 



54 IF(J .EQ. 1) DTDX(I)— DT 



62 12 CALL LTERPS ( I, J,NANGLE, EX, FY, X( I) , FYX ( I, J ) ,D ( J ) ) 



107 0TDY(I)"<T(I,3)-T<I,2))/(WY(I,2)-WY(I,3)) 



133 DFYDX(I)--D<1) 



140 FV(I)-FYX(I,1) 



147 0FYDY(I)-(FYX(I,2)-FYX(I,3)>/(WY(I,2)-WY(I,3))*(-1.) 



174 0FXDY(I)»(FXX(I,2)-FXX(I,3))/(HY(I,2)-WY(I,3))*(-1.> 



221 14 CONTINUE 



223 PRINT 16, TITLE 



230 1€ F0RMAT(1H1,8A10) 



230 PRINT 15,TIDE(ITIDE),FEXT(IC),SEP(IC).TEN?1(IC),TENS2(IC) 



261 15 FORMAT(///X*FOR THE CONDITIONS — */ 



1 5X+TI0E » ♦F5.2/ 



2 5X*EXTERNALLY APPLIED HORIZONTAL FORCE, FEXT" ♦F1C.3*LB.*/ 



3 5X*H0RIZ0NTAL ANCHOR SEPERATION, SEP" *F10.3*FEET*/ 



4 5X*N0MINAL TENSION IN CABLE 1 •♦F10.3* LB.*/ 



5 5X*N0MINAL TENSION IN CABLE 2 ■♦F10.3* LB.*//) 

 261 PRINT 18 



265 18 F0RMAT(/5X*CABLE NO. FX*10X*F Y*11X, 1HX,11X* Y*9X*Df XDX*7X, 



2 ♦DFXDY*7X*DFY0X*7X*DFYDY*,5X, ♦TENSION*//) 

 265 DO 20 1-1,2 



270 20 PRINT 25, I, F( I ) , F V( I ), X ( I ),UY( 1,1 ) ,DFXDX( I ) ,DFYDX( I ) , DFXDY ( I ) , 



2 DFYDY(I),T(I,1) 

 337 25 F0RMATt9X,Il,4(2X,F10.3),5(X,E11.3)) 



C***NOW CALCULATE FORCES AND SPRING COKSTANTS FOR THE BREAKWATER 



C S-SWAY MOTION •fX DIRECTION FE^T 



C H-HEAVE MOTION +Y DIRECTION 



C R-ROLL MOTION COUNTERCLOCKWISE RADIANS 



C FS-FORCES CAUSING SWAY DUE TO THE MnPING LINES 



C FH-FORCES CAUSING HEAVE DUF TO MOORING LINFS 



C EMR«MOMENTS CAUSING ROLL DUE TO MOORING LINES 



C CHANGE YCAB TO BE DIST TO C6 IN Y DIRECTION 

 337 YCAB(1)-YCAB(1)-YCG 



345 YCAB(2)"YCAB(2)-YCG 



352 FS-F(1)*F(2) 



360 FH-FV(1)*FV(2) 

 365 EMR»FV(l)*XCAB(l)+FV(2)*XCAB(2)-F(n*YCAB(l)-F{2)*YCAB(2) 



C***CALCULATE CHANGE IN TENSIONS WITH BREAKWATER MOTIONS 



Table B-3. Continued 



101 



