Advanced methods for ship motion and wave load prediction

SSC-333

ADVANCED METHODS FOR

SHIP MOTION AND WAVE

LOAD PREDICTION

This &cumcnt has been approved for public release d sale; its

distribeticst is unlimited

SHIP STRUCTURE COMMITTEE

RADM J. D. Sipes, USCG, (Chairman) Chief, Office of Marine Safety,

Security and Environmental Protection U. S. Coast Guard

Mr. Alexander Malakhoff Director, Structural Integrity

Subgroup (SEA 55V) Naval Sea Systems Command

Dr. Donald Uu Senior Vice President American Bureau of Shipping

CONTRACTING OFFICER TECHNICAL REPRESENTATIVES Mr. William J. Siekierka Mr. Greg D. Woods

SEA 55Y3 SEA 55Y3

Naval Sea Systems Command Naval Sea Systems Command

SHIP STRUCTURE SUBCOMMITTEE

THE SHIP STRUCTURE SUBCOMMITTEE acts for the Ship Structure Committee on technical matters by providing technical coordinating for the determination of goals and objectives of the program, and by evaluating and interpreting the results in terms of structural design, construction and operation.

U.S. COAST GUARD

Dr. John S. Spencer (Chairman) CAPT T. E. Thompson

Mr. David L. Motherway CDR Mark E. NOII

NAVAL SEA SYSTEMS COMMAND Mr. Robert A. Sielski

Mr. Charles L. Null Mr. W. Thomas Packard Mr. Allen H. Engle

MARITIME ADMIN IStBAIIQN

Mr. Frederick Seibold Mr. Norman O. Hammer Mr. Chao H. Lin Dr. Walter M. Maclean

U.S COAST GUARDACADEMY

LT Bruce Mustain

U.S. MERCHANT MARINE ACADEMY Dr. C. B. Fm

U. S. NAVAL ACADEMY Dr. Ramswar Bhattacharyya

SIIJJJ4MSlTYOF NEW YORK

MARITIME COLLEGE Dr. W. R. Porter

WELDING RESEARCH COUNCIL Dr. Glen W. Oyler

SHIP STRUCTURE COMMITTEE

THE SHIP STRUCTURE COMMITTEE is constituted to prosecute a research program to improve the hull structure of ships and other marine structures by an extension of knowledge pertaining to design, materials and methods of construction

Mr. H. T. Haller

Associate Administrator for Ship-building and Ship Operations Maritime Administration Mr. Thomas W. Allen Engineering Officer (N7) Military Sealift Command

CDR Michael K. Parmelee, USCG, Secretary, Ship Structure Committee U. S. Coast Guard

MILITARY SEALIFT COMMAND Mr. Glenn M. Ashe

Mr. Michael W. Tourna Mr. Albert J. Attermeyer Mr. Jeffery E. Beach

AMERICAN BUREAU OF SHIPPING Mr. John F. Conlon

Mr. Stephen G. Arntson Mr. William M. Han2alek Mr. Philip G. Rynn

SHIP STRUCTURE SUBCOMMITTEE LIAISON MEMBERS

Tl. A ..C. ii F CENCES

MARINE BOARD Mr. Alexander B. Stavovy

NATIONAL ACADEMY OF SCIENCES COMMITTEE ON MARINE STRUCTURES Mr. Stanley G. Stiansen

SOCIETY OF NAVAÇARCHITECTS AND MARINE

ENGINEERS-HYDRODYNAMICS COMMITTEE Dr. William Sandberg

AMERICAN IRON AND STEEL INSTITUTE Mr. Alexander D. Wilson

Member Agencies:

United States Coast Guard Naval Sea Systems Command Maritime Administration American Bureau of Shipping Military Sealift Command

Ship

Structure

Committee

An Interagency Advisory Committee Dedicated to the Improvement of Manne Structures

August 2, 1990

ADVANCED METHODS FOR SHIP

MOTION AND WAVE LOAD PREDICTION

Advanced numerical methods are needed by ship designers to better

predict and simulate ship motions and hull girder loads.

Co:plex

structural loading problems such as bottom slamming,

bow fiare

impact, and green water on deck cannot be satisfactorily analyzed

using linear strip theory.

This report provides a numerical method for predicting transient

three-dimensional

hydrodynamic pressures

and

resulting

loads.

This work is based on an initial level of investigation and

development,

and will require further testing, validation,

and

refinement of the numerical methods and computer programs.

SIP

Rear Admiral, U.

S. Coast Guaid

Chairman, Ship Structure Committee

-

33

Address Correspondence to:

Secretary, Ship Structure Committee U.S. Coast Guard (G-Mm)

2100 Second Street SW. Washington, D.C. 20593-0001 PH: (202) 267-0003 FAX: (202) 267-0025 SSC- 333 SR- 1277

I C

SUBROUTINE FOIST

COMMON/BD/XPAN(120),YPAN(120),ZPAN(120),AREA(120),ST(j.20),

*

_{ACN(120),ACNW(120),AN(120,3),E(120),P(120,6L.PRFS(12o),}

* STOLD(120), PX(120, 6) COMMON/BD2/XPT( 150 ) . YPT( 150), ZPT( 150 ) WRF( 150 ) WRFR ( 150) *

KK(150,4)

COMMON/A/NPAN. NPT, CEE, RHO, NKX, NKY8 EYE, DT8 TIM, UFWD COMPLEX A8 B, EYE

DIMENSION XPSL(3,4), XPSLR(3,4).PEB(120, 120)

COMMON/PTST/ARE4(2001 4), X4(200, 4), Y4(200, 4), Z4(2001 4) *

,SEL(200,4)

DO 1500 .J1,NPAN

ARE4(J, 4)-1. O JT=4 IF(RK(J14). EQ. 0) JT=3 DO 1500 %J'J1lJT %J2= i IF(JJ. LT. siT) J2=JJ+1 KF=KK(J, JJ) KG=KK(J, U2)

X4(J,JJ)=(XPT(KF)+XPT(KC)+XPAN(J))/3.0

Y4(J,JJ)=(YPT(KF)+YPT(KC)+YPAN(J))/3.0

Z4(J, JJ)=(ZPT(KF)+ZPT(KC)+ZPAN(.J) )/3. O

AF=XPT(KF)XPAN(J)

BF=YPT(KF)YPAN(.J)

CF=ZPT(KF)ZPAN(J)

AG=XPT(KG)XPAN(J)

BQ=YPT(KO)YPAN(J)

C0=ZPT(RQ)ZPAN(U)

CALL SELF(AF1 BE, CF, AG, BG, CG, FEE)

SEL(J, JJ)=FEE

CR=AF*BGBF*AQ

AR=BF*C0CF*BG

BR=CF*AGAF*CQ

ARE4(J, JU)=0. 5*SQRT(AR*AR+BR*BR+CR*CR)

1500

CONTINUE

DO 127 NJ18NPAN

DO 1277 MJ=1,NPAN

1277 PBB(NJ, MU)=0. 00 P(NJ. i )=0. 00 P(NJ, 2)=0. 00 P (NJ1 3)=0. 00 P(NJ1 4)0. 00 F(NU, 5)=0. 00 P(NJ, 6)0. 00

DO 128 NK=114

ARN=ARE4(NJ, NK) IF(ARN. LT. 0. 0) 00 TO 129 P1=0. 0 P2=0. 0 P3=0. 0 P4=0. 0 P5=0. 0 P6=0. 0 X=X4(NJ, NR)

YY4(NJ, NR)

Z=Z4(NJ, NK)

DO 138 MJ=1,NPAN

DO 138 MK1,4

XF=X4(MJ, MR) YF=Y4(MU, MR) ZF=Z4(MJ, MR) ARM=ARE4(MJ, MR) IF(ARM. LT. 0. 00) 00 TO 138 IF(NJ. NE. MJ) CO TO 140 IF(MK. NE. NR) CO TO 140 FRA=SEL (ft.), MR) /ARM

GO TO 1380

C

SUBROUTINE GE(XF, VF, ZF, J, Vi, V21 V3 V1R1 V2R, V3R, NBT)

COMMON/BDIXPAN(120)1 YPAN(120)1 ZPAN(120)1 AREA(120)1 ST(120),

*

ACN(120),ACNW(120),AN(i20,3),E(1.20),P(120,6),PRFS(120),

*

STDLD(120),PX(12016)

COMMONfBD2fXPT(15O)1VPT(150) ZPT(.50)WRF(150),WRFR(150)

*

,KK(150,4)

COMMON/ARE/RR(500), XZJ(200), YXJ(200), ZYJ(200)

DIMENSION XSA(314),XFA(3)1XSAR(314)

J4=J*4

V10. 00

V2=0. 00 V3=0. 00 V1R=O. 00 V2R=O. 00 V3R=O. 00 XN.J=AN(J, i) YNJ=AN(J, 2) ZNU=AN(J, 3)

NSIDE=4

IF(KK(J,4).EG.0) NSIDE=3

DO 20 UJ=1,NSIDE

J2= i

IFLJJ. LT. NSIDE) J2=JU+1

..J4=J4+i KF=KK CU1 JJ) AF= XP T C KF B F=YPT C KF C F= ZP T C KF R=RR (U4) K0=KK(J, J2)

ANXCAF-XPT(KQ) ) /R

ANY=(BF-YPT(KC) ) IR

ANZ=(CF-ZPT(KG) )/R

A=AF-XF

B=BF-VF

CCF-ZF

TX=XZJ( U)*ANZ-YXU(J)*ANY

TV=VXU(J)*ANX-ZYU(J)*ANZ

TZ=ZVJ(J )*ANY-XZJ( J) *ANX

EX I=A*ANX+B*ANY+C*ANZ

CALL 0O(EX1, R FF, WRF(KF), WRFCKO))

V1=V1+FF*TX

V2=V2+FF*TY

V3=V3+FF*TZ

XSA( 1

UU)=-A/WRF(KF)

XSA(2, JJ)=-BIWRF(KF)

XSA(3, JJ)=-C/WRF(KF)

EX1-R=EX1+2. 0*ZF*ANZ

CR=-CF-ZF

CALL QOCEX1R, R, FR1 WRFR(KF), WRFR(KG)) VI R=V1R-FR*TX

V2R=V2R-FR*TY

V3R=V3R+FR*TZ

XSAR(1, JJ)=-A/WRFR(KF)

XSAR(2, JJ)=-B/WRFR(KF)

XSAR(3, JJ)=CR/WRFR(KF)

20

CONTINUE

0=6. 283185307

IF(J. EQ. NBT) GO TO 84 CALL SDLID(XSA, G, NSIDE)

AGQ=A* X NJ+B* YNJ+ C* Z NJ

G=-SIGN(G, AGO)

84

CONTINUE

CALL SOLID(XSAR1 GR, NSIDE)

AGOR=A*XNJ+B*YNJ-CR*ZNU

GR=SIGN(0R1 AGOR) 85

CONTINUE

7371

FORMAT(' G,GR=', 2F15.5)

V1=V1XNJ*Q

V2=V2+YNJ*Q

V3=V3+ Z NJ*Q

V1R=V1R+XNJ*QR

V2R=V2R+YNJ*OR

V3R=V3R-ZNJ*QR

5590

FORMAT(' V1,V2,V3=',3F15.5)

5591 FORMAT(' V1R1 V2R, V3R', 3F15. 5)

RETURN

END

DO 500 I13

WRITE(3) (AN(J IL. J=1, NPAN)

500

CONTINUE

WRITE(3) (XPAN(J)1U=1,NPAN)

WRITE(3) (YPAN(J).U=1,NPAN)

WRITE(3) (ZPAN(J)1 J=1, NPAN)

WRITE(3) (AREA(J)1J=1,NPAN)

DO 309 JL=1,NPAN

JL..J=JL

CALL QE(XFI YF1 ZF1 JLJ, VX, VY, VZ,VXRJ VYR VZRI JJU)

VX=VX+VXR

VY=VY+VYR

VZ=VZ+VZR

C

COMPUTE NORMAL VELOCITY AT PANEL J DUE TO PANEL iL

E (JL ) =AX*VX+AY*VY+AZ*VZ

C INCREMENT PX MATIRX

FR1=-AREA(J)*VX*AN(J, 1)

FR2=-AREA(J)*VX*AN(J, 2)

FR3=-AREA(J)*VX*AN(J, 3)

PX(JL, 1)=PX(UL1 1)+FR1 PX(JL, 2)=PX(JL, 2)+FR2 PX(JL, 3)=PX(JL, 3)+FR3

PX(UL 4)=PX(JL, 4)+YF*FR3-ZF*FR2

PX(JL, 5)=PX(JL, 5)+ZF*FR1-XF*FR3

PX(,JL, 6)=PX(JL, ¿,)+XF*FR2-YF*FR1

309

CONTINUE

WRITE(99) (E(JL)1JL=11NPAN)

308

CONTINUE

DO 2424 IÇ1,6

2424

WRITE(3) (PX(JL.K),JL=1,NPAN)

CLOSE (UNIT=99)

C

INVERT E MATRIX

CALL MATIN(NPAN)

RETURN

END

SUBROUTINE MATIN(NPAN)

C INVERSI MATRIX

DIMENSION E(120 120),BB(120),EST(120)

OPEN(UNIT=99,FILE='SCR',FORtI'UNFORMATTED',TYPE ='OLD')

DO 120 J=1,NPAN

120

READ(99) (E(J,IL.11,NPAN)

DO 130 J=1,NPAN

DO 11 MM=1,NPAN

EST(Mrl)=0. 00 11

BB(MM)=0.00

BB(U)1. O

EST(J)1. 0/E(J J)

DO 17 NIT1,6

DO 17 K=1NPAN

B=BB(K)

DO 15 I=1,NPAN

15

IF(I.NE.K) B=B-E(K,I)*EST(I)

EST(R)B/E(K. K)

17

CONTINUE

WRITE(3) (EST(R). K=1, NPAN)

130

CONTINUE

RETURN

SUBROUTINE EBD

C INITIALIZE PANELS AND COMPUTE BODY MATRIX

C

COMMON/BD/XPAN(120),YPAN(120),ZPAN(120),AREAU2O),ST(120),

*

ACN(120), ACNW(120), AN(120 3), E(1.0), P(120, 6), PRFS(120),

*

STOLD(120) PX(120,

)

COMMON/BD2/XPT( 150 ) YPT( 150 ), ZFT( 150 ) WRF( 150 ), WRFR C 150 h

*

KK(l5O,4)

COMtIONIA/NPAN, NPT, GEE, RHO, NVX, NKY, EYE, DI, TIM, UFWD

DIMENSION EF (120), EPP (120)

COMPLEX EYE

C READ IN BODY PANEL PARAMETERS

101

FORMAT(415)

100 FORMAT(3F10. 0) 103 FORMAT(3F10. 2>

104

FORMAT(13)

C NUMBER OF POINTS AND PANELS

READ(2,101) NPTINPAN

TYPE 101, NPTSNPAN

C COORDIANTES OF POINTS

READ(2, 100) (XPT(N)YPT(N).ZPT(N),N=11NPT)

DO 7777 N=1,NPT

TYPE 103e XPT(NL'YPT(N)1ZPT(N)

7777

CONTINUE

C DEFINE CORNER PINTS OF EACH PANEL

READ(2, 101) (KK(N, 1),KK(N,2),KK(N,3),KK(N,4),N=1,NPAN)

C COMPUTE PANEL AREAS

DO 150 U=1,NPAN

K1=KK(J 1)

K2KK(), 2)

K3=KK(J, 3)

K4KK(J, 4)

IF(K4.EG.0) GO TO 8

XPAN(J)(XPT(K1)+XPT(K2)+XPT(K3)+XPT(K4))*0.25

YPAN(J)(YPT(K1)+YPT(K2)+YPT(K3)+YPT(K4))*0.25

ZPAN(U)=(ZPT(K1)+ZPT(K2)+ZPT(K3)+ZPT(K4))*0.25

GO TO 9

C TRIANGULAR PANELS

B

XPAN(J)=(XPT(K1)+XPT(K2)+XPT(K3))/3.0

YPANLJ)=(YPT(K1 )+YPT(K2)+YPT(K3) )/3. O

ZPAN(J)=(ZPT(K1)+ZPT(K2)+ZPT(K3))13.0

K4=K3

9

XA=XPT(K3)XPT(K1)

XB=XPT(K4)XPT(K2)

YA=YPT(K3)YPT(K1)

YB=YPT(R4)YPT(K2)

ZA=ZPTK3)ZPT(K1)

ZB=ZPT(K4)ZPT-(K2).

- - - -

-C -COMPUTE PANEL AREAS

AZ=XA*YBYA*XB

AX=YA*ZBZA*YB

AY=ZA*XBXA*ZB

ARE=SQRT (AX*AX+AY*AY+AZ*AZ)

AREA(J)=ARE*0. 50

AN(J, 1)=AXIARE

AN(i, 2)=AY/ARE

AN(U, 3)=AZ/ARE

150

CONTINUE

808

FORMAT(1X, 15, 5F11. 4)

DO 1308 J=1,NPAN

J J J=

CALL PREP(JJJ)

ST(%J)=0. 00

DO 1308 K1,6

PX(J, K)=0. 00 130E

CONTINUE

DO 308 J=1NPAN

J J j= J

AXAN(J, 1)

AY=AN(J, 2)

AZ=AN(J, 3)

XF=XPAN(U)

YF=YPAN(J)

Z F= ZP AN C

DO 157 L=i,NPT

WRF(L)=SQRT( (XPT(L)XF)**2+(YPT(L)YF)**2+(ZPT(L)ZF**2

157

WRFR(L)=SGRT((XPT(L)XF)**2+(ypl(L)YF)**2(ZpTL)+zF)**2)

SUBROUTINE SOLID(XPN, C, NSIDE) DIMENSION CS(4), SN(4). Z(4) XPN(3, 4) Q=-6. 283185308

ACR12=XPN(1, 1)*XPN(1,2)+XPN(2, 1)*XPN(2,2)+XPN(3, 1)*XPN(3,2)

ACR13=XPN(1,1)*XPN(1,3)+XPN(2,1)*XPN(2,3)+XPN(3,1)*XPN(3,3)

ACR23=XPN(1,2)*XPN(1,3)+XPN(2,2)*XPN(2,3)+XPN(3,2)*XPN(3,3)

IF(NSIDE. EQ. 4) 00 TO 40 0=-3. 141592659 CS( i )=ACR23-ACR13*ACR12

CS(2)=ACR13-ACR12*ACR23

Cs (3) =ACR 12-ACR23*ACR 13 SN( l)=XPN( 1. 1)*( XPN(2, 2)*XPN(3, 3)-XPN(3, 2)*XPN(2, 3) +

XPN(2,1)*(XPN(3,2)*XPN(1,3)-XPN(i,2)*XPN(3,3))+

+

XPN(3, 1)*(XPN(1,2)*XPN(2,3)-XPN(2,2)*XPN(1,3))

SN(2)=SN( 1) SN(3)=SN( 1) SN(4)=0.

CO TO 50

40

ACR14=XPN(1, i)*XPN(i,4)+XPN(2, 1)*XPN(2,4)+XPN(3, 1)*XPN(3,4)

ACR24=XPN(1, 2)*XPN(1, 4)+XPN(2, 2)*XPN(2, 4)+XPN(3,2)*XPN(3, 4) ACR34XPN(1, 3)*XPN(1, 4)+XPN(2, 3)*XPN(2, 4)+XPN(3, 3)*XPN(3, 4) CS( i )=ACR24-ACR14*ACR12 CS (2) =ACR 13-ACR23*ACR 12 CS (3) =ACR24-ACR34*ACR23 CS (4) =ACR 13-ACR34*ACR 14

B241=XPN(2,2)*XPN(3,4)-XPN(3,2)*XPN(2,4)

B242=XPN(31 2)*XPN(1, 4)-XPN(11 2)*XPN(3, 4)

13243=XPN(1,2)*XPN(214)-XPN(2,2)*XPN(1,4)

8131=XPN(2, i )*XPN(3, 3)-XPN(31 i )*XPN(2 3) B132=XPN(3, I )*XPN( 1 3)-XPN( 1, 1 )*XPN(3, 3)

13133=XPN(1,I)*XPN(2,3)-XPN(2,1)*XPN(1,3)

SN(1)=XPN(1. 1)*B241+XPN(2, 1)*B242+XPN(3, 1)*B243

SN(2)=-(XPN(1,2)*B131+XPN(2,2)*B132+XPN(3,2)*3133)

SN(3)=-(XPN(i,3)*B241+XPN(2,3)*13242+XPN(3,3)*B243)

SN(4)=XPN(1,4)*B131+XPN(2,4)*B132+XPN(3,4)*B133

50

CONTINUE

D TYPE 8844, SN(1), CS(i) D

TYPE 8844, SN(2),CS(2)

D TYPE -8844. SN(3)--CS(3)- - --- - -. -D TYPE 8844, SN(4L. CS(4)

8844

FORMAT('

SN, CS='.2F15. 8) SUM=SN(1 )-i-SN(2)+SN(3)+SN(4)

IF(ABS(SUMLCT. 0. 01) 00 TO 25

IF(ABS(CS(1)).CT.ABS(SN(1))) 00 1025

IF(ABS(CS(2)).GT.ABS(SN(2))) COTO 25

IF(ABS(CS(3)).QT.ABS(SN(3))) GO 1025

1090

Q=SUM*. 25

IF(NSIDE. EQ. 3) 0=SUPI*. 1666666667

RETURN

25

ST=SN(NSIDE)

DC 30 I=1,NSIDE

IF( (ABS(CS( I)). LT. 9E-8). AND. (ABS(SN(I)). LT. .9E-05))

+

CD TO 1090

IF(ST*SN(I). LT. 0.) 00 TO 1090 ST=St'4( I) C2=CS( I) /SQRT(SN( I )**2+CS( I )**2)

C=Q+ACOS(C2)

30

CONTINUE

RETURN

END

SUBROUTINE PREP(J)

COMMON/BD/XPAN(120),YPAN(120),ZFAN(120),AREA(120)1ST(120),

*

ACN(120)1ACNW(120),AN(120,3L.E(120).P(12O16),PRFS(120)

*

STOLD(120),PX(120,6)

COMMONIBD2/XPT(150)1 YPT(150), ZPT(150).

WRF(150), WRFR(150)1 *

KK(15014)

COMMON/ARE/RR(500)1 XZJ(200). YXJ(200)

ZYJ(200)

ZYT=0. O

YXT=O. 00

XZT=0. 00

,J4=J*4

JT=4

IF(KK(.J,4). EQ. 0) UT=3

DO 20 JJ=11UT

,J4=J4+j .12=1 IF(JJ. LT. JT) J2=JU+1 KF=KK(J, JJ) KC=KK(i .12)

AGXPT(VG)

B G=YP T ( KG

CG=ZPT(KG)

AF=XPT(KF)

B F=YP T ( KF

CF=ZPT(KF)

R=SGRT( (AF-AG)**2+(BF-BQ)**2+(CF-CG)**2)

XT=AF-XPAN ( U) YT=BF-VPAN (U) ZT=CF-ZPAN( U)

ANX=(AF-AG)/R

ANY C SE-BG) IR ANZ=(CF-CG) IR

DOT=ANX*XT+ANY*YT+ANZ*ZT

XT=XT-DOT*ANX

YT=YT-DOT*ANY

ZT=ZT-DOT*ANZ

ZYT=ZVT+ZT*ANY-ANZ*YT

YXTYXT+YT*ANX-ANY*XT

XZT=XZT+XT*ANZ-ANX*ZT

RR(U4)=R

20

CONTINUE

XZU(U)=SIGN(AN(U12)XZT)

YXJ(J)=SIGN(AN(J,3),YXT)

ZYU(U)=SIGN(AN(J,1),ZYT)

RETURN

SUBROUTINE SELF(AF, BE, CF, AG, BG, CG, FEE)

REAL LB21,LA21

ASG=AF*AF+BF*BF+CF*CF

BSQAG*AO+BO*BO+CQ*CO

ADBAF*AG+BF*BO+CF*CG

ADB2=ADB+DB

ASAS= (AF*BG-BF*A0 ) **2+ (CF*BG-BF*CG ) **2+

(AF*CG-BF*AG)**2

FF=0. 00 DO 15 MK=1, 10

DO 15 NK=1,MK

LA2 I =FLOAT C NK-MK) A2SG=ASG*LA2 1 *LA2 i

DO 15 ML=1111-MK

DO 15 NL=1, 11-ML

- LB21=FL0AT(NL-ML) - --- -- - -IF(LA2I. NE. 0. )

GO TO 5

IF(LB21. LT. 0. ) GO TO 5

GO TO 15

5

R=SGRT(A2SQ+A13D2*LA21*L521+BSQ*L521*L521)

FF=FF+1. O/R 15

CONTINUE

FEE=FF*ASAS*O. 002

RETURN

END

C

PROGRAM HYDREX3

CHARACTER*25 PANFIL, MATFIL

COMMON/BD/XPAN ( 120 ), YPAN ( 120), ZPAN C 120 ), AREA ( 120 ) * ST(120), ACN(120)1 ACNW(120), AN(120, 3) E(120)1 P(120, 6),

* PRFS(120)1 STOLD( 120), PX(1201 6)

COMMON/FS/AKZ(100, 100), SS(100, 100), CC(100, 100)

*

DKX(100),DKY(100),AKX(100),AKY(100)

COMMON/2D2/XPT( 150 ) YPT( 150 ) ZPT( 150 ), WRF( 150 )

*

WRF(150),KK(15014)

COMMON/A/NPAN, NPT, GEE, RHO1 NKX, NKY, EYE, DT1 TIM, UFWD

DIMENSION PF(6), PB(6), P1(6) CONtION/WAVEX ¡OMEGA COMMON/BP/BPRES( 120), TPRES 120)

COMPLEX EYE

EYE=(0. 0, 1. 0) C

TYPE 1

ACCEPT 4,PANFIL

TYPE 2

ACCEPT 4,MATFIL

OPEN(UNIT=2, FILE=PANFIL, TYPE='OLD")

OPEN(UNIT=3. FILE=MATFIL. FORM='UNFORMATTED', TYPE'NEW')

OPEN(UNIT=991 FILE='X. DAT' FORM='L)NFORMATTED'1 TYPE ='NEW)

CALL EBD

CALL POTSI

WRITE(6, 6)

WRITE(6,7) J,AN(J, 1).AN(J12)1AN(J3)1

*

XPAN(.J),YPAN(U),ZPAN(J),AREA(J)

150

CONTINUE

STOP

FORMAT(' Input name of EPAN] file >'$)

2

FORMAT(' Input name of [MAT] file >'$)

4

FORMAI(A)

6 FORMAT( 'J',X. 'NX',9X, 'NY',9X, 'NZ',9X, 'XP',9X,

* 'YP'..9X, 'ZP',9X1 'AREA')

7

FORMAT(1X,1517F11.4)

180

CONTINUE

190

CONTINUE

200 CONTINUE

210 D=A(NN)*D

IF (A(N1N).EQ. 0. 0) GO TO 270 Z(N)=1. /A(N.N) C

C..

OBTAIN SOLUTION 13V BACK SUBSTITUTION.

C

DO 220 L1,N1

B(N L)=Z(N)*B(N, L)

220 CONTINUE

IF (NM1. EQ. 0) 00 TO 260

DO 250 K=1NM1

J=N-K

j1='.j+1

DO 240 L=1,N1

w=0.

DO 230 I=J1N

W=A(J, I)*B(I,L)+W

230

CONTINUE

B C J, L (B C J, L )-W )*Z (J)

240

CONTINUE

250 CONTINUE

260 LNEGT=0

IF (ABS(D).GE. ERROR) RETURN

LNEQT= i

WRITE (OUTPUT,280) D,ERROR

RETURN

C

C.. SINGULAR MATRIX--MAXIMUM ELEMENT

IN COLUMN IS ZERO.

C

270 LNEQT=3

WRITE (OUTPUT, 290)

RETURN

C

280 FORMAT (1BHO*** DETERMINANT =.1PE13.5,24H1

ERROR SPECIFICATION

=1

1E13. 5)

290 FORMAT (20H0*** SINGULAR MATRIX)

C

FUNCTION LNEQT(MI N, Nl, A1 B, ERROR 2)

C

C.. SOLVES SIMULTANEOUS LINEAR EQUATIONS BY GAUSSIAN REDUCTION.

C.. SOLVES

A*X5

FOR

X ¡

AND STORES THE

X

VECTOR(S) IN

B

C

REAL ACM, M), B(M, M) ,Z(M), ERROR RMAXI RNEXT, W

C

CS

C5

INPUT AND OUTPUT LOGICAL UNITS..

COMMON/ID/INPUT1 OUTPUT, BIF. OFF, COF INTEGER OUTPUT, BhF, OFF, COF

CS C D=1. O

NM1=N-1

IF (NMI. EQ. 0) 00 TO 210

DO 200 J=1,NM1

J1=J+1

C

C.. FIND ELEMENT OF COL J, ROWS J-N, WHICH HAS MAX ABSOLUTE VALUE.

C

LMAX=J

RMAX=ABS(A(J, J))

DO 110 K)1,N

RNEXT=ABS(A(K, J))

IF (RMAX. QE. RNEXT) GO TO 110

R MA X =RNE XT

LMAXK

110

CONTINUE

IF (LMAX.NE.J) GO TO 120

C

C..

MAX ELEMENT IN COLUMN IS ON DIAGONAL

C

IF (A(J,J)) 150,270,150

C

C.. MAX ELEMENT IS NOT ON DIAGONAL. EXCHANGE ROWS J AND LMAX.

C 120

DO 130 L=UN

W=A(J, L) A(J, L)=A(LMAX, L) A(LMAX, L)=W 130

CONTINUE

DO 140 L=1,N1

W=E(J, L) B (J, L)=B (LMAX, L) B(LMAX, L)=W 140

CONTINUE

D=-D

C

C..

ZERO COLUMN J BELOW THE DIAGONAL.

C 150

D=A(J,J)*D

Z(J)=1. 0/A(J, J)

DO 190 k=J1,N

IF (A(K,J)) 160,19O160

160

W=-Z(J)*A(K,J)

DO 170 L=J1N

A ( K, L ) =W*A ( J, L) +A (K, L) 170

CONTINUE

DO 180 L=i.,N1 B C K, L ) W*B C J, L )+E (K, L)

C..

ZERO ROW J TO RIGHT OF DIAGONAL.

C 170

D=A(J,J)*D

V.1=1.

0/A(J J)

Z(J)=V

DO 200 K=J1NOE

IF (A(J,K).EG. 0. 0) 00 TO 200 W=-V*A(J, K)

DO 180 L=J1,NOE

L, K) =W*A( L1 J) +A (L, K) 180

CONTINUE

DO 190 L=1,N1

B ( K, L ) =W*B J, L ) +B (K, L) 190

CONTINUE

200

CONTINUE

210 CONTINUE

D=A(NOE, NOE)*D

IF (A(NOEINOE).EQ. 0. 0) 00 TO 260

Z(NOE)=1. 0/A(NOE1 NOE)

C

C.. OBTAIN SOLUTION ¡3V BACK SUBSTITUTION.

C

DO 220 L=1,N1

B(NOE, L)=Z(NOE)*B(NOE, L)

220 CONTINUE

DO 250 K=i,NM1

J=NOE-K

JIJ+1

DO 240 L1,N1

W=0. O

DO 230 I=J1,NOE

W=A( I, J)*B( I, L)+W

230

CONTINUE

B(J, L)=(B(J,L)-W)*Z(J)

240

CONTINUE

250 CONTINUE

IF (ABS(D).LT. ERR) GO TO 270

C

JULIET=0

C

C.. NO PROBLEMS DURRING THIS EXECUTION.

C

RETURN

C

C.. SINGULAR MATRIX--MAXIMUM ELEMENT IN ROW IS ZERQ.

C

260 JULIET=3

WRITE (OUTPUT, 280)

RETURN

C

C.. ABSOLUTE DETERMINANT VALUE LESS THAN ERROR VALUE..

C

270 JULIET1

WRITE (OUTPUT290) DIERR

RETURN

C

280 FORMAT (20H0*** SINGULAR MATRIX)

290 FORMAT (1BHO*** DETERMINANT =,1PE13.5,24H,

_{ERROR SPECIFICATION}

1E13. 5) C

COMMON /

/

_{HAO(24,8),BLOGP(20,20.4),CONH(40,3),CROLL(40}

1.40, 2) HEAVI(20, 20, 4>, EJI. CZRI, CZLI, SZRI, SZLII RARI1 RALI, RBRI, RELI. 2CLI, CRO, SLI, SRI1 Z(40)

N2NOE/2

D=1. O C

C.. COMPLETE THE MATRIX A.

C

DO 120 J=1,N2

L=N2+J

DO 110 I=1,N2

K=N2+I

A( I, L)=-A(K, J) A(K, L)=A(11 J) 110

CONTINUE

120 CONTINUE

C

NM1=NOE-1

C

C.. SOLVE

AT*XB

FOR X ,

WHERE

AT

IS THE TRANSPOSE OF THE

C

MATRIX

A -

STORE THE

X

VECTOR(S) IN

B

C

DO 210 J=1NM1

Jj =J+ j

C

C...

FIND ELEMENT OF ROW J,

COLS J--N, WHICH HAS MAX ABSOLUTE VALUE.

C

LMAX=J

RMAX=ABS(A(J1 J))

DO 130 K=U1,NOE

RNEXT=ABS(A(J, K))

IF (RMAX. QE. RNEXT) GO TO 130

RMAX=RNEXT

LMAX=K

130

CONTINUE

IF (LMAX. NE. J) GO TO 140

C

C..

MAX ELEMENT IS ON DIAGONAL.

C

IF (A(J.U)) 170,260,170

C

C.. MAX ELEMENT IS NOT ON DIAGONAL.

C.. EXCHANGE COLUMNS J AND LMAX.

C

140

DO 150 L=J,NOE

W=A(L, J)

A(L, J)=A(L1 LMAX)

A(L. LMAX)W

150

CONTINUE

C

C.. EXCHANGE ROWS U AND LMAX.

C

DO 160 L=1.N1

W=B (J, L) B(J, L)=B(LMAX, L)

B(LMAX, L)W

160

CONTINUE

D=-D

C

C IF (R. LT. 1. )

GO TO 130

TEST=. 1*TEST IF (R. LT. 2. )

00 TO 130

TEST. 1*TEST

IF (R.LT.4. )

GO TO 130

TEST. 1*TEST

130 SUMC=QAMMA+AL+Y

SUMS AT+ X

TC=Y

TS=X

COX=1. DO 140 K=2, 501

TOTC

FACT=COXfFLOAT (K) **2

COX=K

TC=FACT* (Y*TCX*TS)

TS=FACT* (y*TS+X*TO) SUMC=SUMC-i-TC SUMS=SUMS+TS

IF ((ABS(TC)+ABS(TS)).LE. TEST) 00 10 150

140 CONTINUE

WRITE (OUTPUT, 190) XXSTA,X.Y,WN

STOP 12

150 C IN=E* ( C*SUMC+S*SUMS)

SON=E* ( S*SUMCC*SUMS)

RA=ALC IN

RB=ARQ+SON

RETURN

C

190 FORMAT (59H0*** NONCONVERGENT EXPONENTIAL INTEGRAL FOR STATION AT

I

X =,F13.5/22H *** PARAMETERS-- X =,1PE13.515H, Y =E13.3, 15H, WA

2VE NUMBER

.E13.5)

C

END

FUNCTION JULIET (N1,A,B)

C

C. -

SOLVES SIMULTANEOUS LINEAR EQUATIONS BY GAUSSIAN REDUCTION.

C.. FOR THE SPECIALIZED MATRICES IN THE SUBROUTINE WINE.

C

REAL A(40, 1),13(40, 1)

C2

C2

THE A MATRIX MUST BE DIMENSIONED WITH EXACTLY 40 ROWS AND

C2

AT LEAST 40 COLUMNS.

THE B MATRIX MUST ALSO BE DIMENSIONED

C2

WITH EXACTLY 40 ROWS AND AT LEAST Nl. COLUMNS.

C2 C5

C5

INPUT AND OUTPUT LOGICAL UNITS..

COMMON/IO/INPUT1 OUTPUT, 131F. OFF, COF INTEGER OUTPUT, 131F, OFF, COF

CS CS

COMMON/GEOMETRY/MSTA, LPTS(25)s YOFF(25, 25). NAFT1 XAFT(25),

*

_{YAFT(25),NFWD,XFWD(25),YFWD(25),XQFF(25,,}

*

_{ZOFF(25,25), XFPERP, XAPERP,SHIPL,SHIPB,SHIPT,}

Yl (21, 25). ZL.JL(25). WL(25), INPTS(25), XWLF. XWLA. XXFI

* _{XXAJ TAN, NON, NOE. NWLI CR, XXFWD. XXSTA. !AFT,DX, Dxl,}

*

_{DX2,Z2(21),Y2(21),Zz(2o),yy(2o),SNE(2o),csE(20),}

* _{DEL(20), ROL(20), ADJUST. WMAX, YMAX, ZMAX, AREA. VERT}

COMMON 1/ HA1,SA1,RA1,CA1,HV1,SV1,RV1,CV1,RHO2,RSIQ,WN,W1,W2,ERR,Z

1RI. YRI, EUT

SAS=SAS+PAS*DS I

R AR =R AR +P AR*DFR

C6

CCAPAR*DS I+PAS*DFR+CCA

C6

INTEGRATION TO OBTAIN FORCE VELOCITY COEFFICIENTS..

C6 HVH=HVH+PVH*DC I SVS=SVS+PVS*DS I R VR R VR +P VR*DFR

CCVPVR*DSI+PVS*DFR+CCV

C6

Có

AT THIS POINT THE PRESSURES REQUIRE THE MODIFICATIONS NOTED

C6

ABOVE TO GIVE THE DIMENSIONAL VALUES.

C6

THE INTEGRATIONS OF THE PRESSURES ARE COMPLETED ELSEWHERE

C6

IN THE PROGRAM TO GIVE DIMENSIONAL FORCE COEFFICIENTS.

C6

RETURN

C

END

f

SUBROUTINE ROMEO (C, S, RA1 RB, CIN1 SON)

C

Cl

EXPONENTIAL INTEGRAL WITH COMPLEX ARGUMENT.

Cl

C2

THE ARGUMENT IS SUPPLIED THROUGH BLANK COMMON AS THE VARIABLES

C2 X

AND

Y.

C2

C3

PARAMETERS AND VARIABLES..

C3 X

-- REAL PART OF ARGUMENT

C3 Y

-- IMAGINARY PART OF ARGUMENT

C3 E

---EXP( -Y)

C3 C

-- COS( X

C3 S

-- SIN( X

C3

CIN

-- REAL RESULT

C3

SON

-- IMAGINARY RESULT

C3 RA

-- ALOG( X**2 + y**2 )/2.0 - GIN

C3 RB -- ATAN C X/V ) - P 1/2. 0 + SON

C3 CS

CQMMON/GEOMETRY/MSTA, LPTS(25), YOFF(251 25), NAFTS XAFTC25), *

YAFT(25),NFWD,XFWD(25),YFWD(25),XOFF(25),

*

ZOFF(25,25)1XFPERP1 XAPERPSHIPL1SHIPB,SHIPT,

* Vi (21, 25), ZWL(25),WL(25)1 INPTS(25), XWLF, XWLAS XXF,

* XXA, TAN, NON, NOE,NWL, CR, XXFWDI XXSTA. XXAFT, DX, Dxl,

*

DX2,Z2(21),Y2(21),ZZ(20)1YY(20).SNE(20),CSE(20),

* DEL(20)5 ROL(20), ADJUSTS WMAX, YMAX, ZMAX, AREA, VERT

COMMON f

/ HAl, SAl, RAi, CAl, HV1, SV1, RV1, Cvi,

1

RH025 RSIG, WN,

Wi, W2, ERR,

X,

Y, E

CS

CS

INPUT AND OUTPUT LOGICAL UNITS..

COMMON/rn/INPUT1 OUTPUT, BIF, OFF, COF INTEGER OUTPUT, ElF, OFF, COF

DATA

GAMMA /0. 5772 1564 90153 286O 0512/

DATA

HALFPI/1. 5707 9326 79489 66192 31322/

C AT=ATAN2(X, Y) ARQ=AT-HALFP I

C=COS(X)

SSIN(X)

R=X*X+Y*Y

AL=0. 5*ALQG(R) TEST=0. 00001

SZLI=SZLT

110 CONTINUE

RETURN

END

SUBROUTINE SONG (HAI, HOT, RAI, ROT1 I)

DIMENSION HAI(i), HOT(1), RAI(I), ROT(1), PP(6)

COMMON/GEOMETRY/MSTA1 LPTS(25)1 YOFF(25. 25), NAFT, XAFT(25), *

YAFT(25),NFWD,XFWD(25),YFWD(25),XOFF(25),

* ZOFF(25. 25) XFPERP, XAPERP, SHIPL, SHIPB, SHIPT1

* Yl (21, 25) ZWL(25) WL(25), INPTS(25), XWLF, XWLAI XXF,

* XXA. TAN1 NON, NOE, NWL, CR, XXFWD, XXSTAI XXAFT,DX, Dxl,

*

DX2,Z2(21),Y2(21),ZZ(20),YY(20)1SNE(20),CSE(20)

* DEL(20). ROL(20), ADJUST, WMAX, YMAX, ZMAX, AREA, VERT

COMMON II HAH, SAS. RAR, CCA, HVH, SVS, RVR, CCV, RHO2, RSIG, WN, Wi, W2. ERR, Z 1RI. YRI, EJT

COMMON /

f

HAO(24,8),BLOGP(20,20).YLOGP(20,20),BLOGM(20

120) YLOQM(20, 20) CONH(40). CONR(40, 2)

PAHO. O

PAS=0. O

PARO. O

PVHO. O

PVS=0. O

PVRO. O

DO 110 J11NON

NJ=NON+J

PAH=PAH+CONH(J)*HOT(J)-CONH(NJ)*HAI(J)

PAS=PAS+CONR(U, i)*ROT(J)-CONR(NJ, 1)*RAI(J)

PAR=PAR+CONR(J, 2)*ROT(J)-CONR(NJ, 2)*RAI(J)

PVH=PVH+CONH(J)*I-1A1 (J)+CONH(NJ)*HOT(J)

PVS=PVS+CONR(J,1)*RAI(J)+CONR(NJ,1)*ROT(J)

PVR=PVR+CONR(J,2)*RAI(J)+CONR(NJ,2)*ROT(J)

110 CONTINUE

DDDDEL ( I) DCI=CSE( I )*DDD DSI=-SNE( I )*DD DFR=ROL( I )*DDD

C6

C6

THE PRESSURES ON THIS SEGMENT OF THE CYLINDER MAY BE CALCULATED.

ce

Cé

THE PRESSURES IN PHASE WITH THE SINUSOIDAL DISPLACEMENT ARE..

Ce

C6

HEAVE -- PAH = PAH*RHO*ESIG*ESIO

SWAY

-- PAS = PAS*RHO*ESIQ*ESIG

C6

ROLL

-- PAR = PAR*RHO*ESIG*ESIQ

Ce

C6

OF COURSE THE ACCELERATION COMPONENTS OF THE FORCE ARE EQUAL

Ce

IN MAGNITUDE TO THE ABOVE, BUT HAVE THE OPPOSITE SIGN.

Ce

C6

THE PRESSURES IN PHASE WITH THE SINUSOIDAL VELOCITY ARE..

Co

CO

HEAVE -- PVE-{ = PVH*RHO*ESIO*ESIG

C6

SWAY

-- PVS = PVS*RHO*ESIC*ESIG

ce

ROLL

-- PVR = PVR*RHO*ESIQ*ESIG

ce

ce

ce

INTEGRATION TO OBTAIN FORCE ACCELERATION COEFFICIENTS..

CO

HAHHAH+PAH*DC I

C

C

1RT, YRT, EJT

COMMON /

f HAO(24, 8), BP(20, 20), YP(20, 20), BM(20, 20), YM(2

10,20) CONH(40), CONR (40, 2), CROLL(40, 40), CHEAV(40, 40), HEAVI (20e 20)1H 2EAVT(20, 20), ROLLI (20, 20) ROLLT(20, 20), EJI, CZRI, CZLI, SZRI, SZLI, RARI

3, RALI, RBRIS RI3LI, CLI, CRO, SLI SRI

DATA TP 1/6. 28318530717958/

YYI=YY( I)

ZZI=ZZ(I)

SISNE( I)

CI=CSE( I) DO 110 J=l, NON

XRT=WN*( ZZIZ2(J+l))

YRT=WN*(YYI+Y2(J+1))

EJT=EXP (YRT)

CALL ROMEO (CZRT, SZRTJ RART, RBRTICRT, SRI)

XRT=WN*(ZZI+Z2(J+1))

CALL ROMEO (CZLT. SZLTI RALT, RBLTS CLT, SLT)

C J=C SE ( J

SJ=SNE(J)

SSS=SI*C')

TTT=SJ*CI

ULJU=C I*CJ VVV=S I *SJ

C IPJ=UUUVVV

S IPJ=SSS+TTT

SI MJ=SSSTTT

C I MJ=UUU+VVV

SSS=SIMJ*(CLICLT)CIMJ* (SLISLT)

TTT=SIPJ*(CRQCRT)C LPJ*(SRISRT)

UUU=SJ*(RALIRALT)+CJ*(RBLTRBLI)

VVV=SJ* (RARIRARI> +CJ* (RBR 1RERT)

WWW=EJT*(SZRT*CIPJCZRT*SIPJ)EJI*(SZRI*CIPUCZRI*SIPJ)

RRR=EJT*(SZLT*CIMJCZLT*SIMJ)EJI*(SZLI*CtMJCZLI*SIMi)

OQQ=EJI*(SZRI*CJCZRI*SJ)EJT*(SZRT*CJCZRT*SJ)

PPP=EJI*(SZLI*CJ+CZLI*SU)EJT*(SZLT*CJ+CZLT*SJ)

CHI(J)=BLOOP(J)+2.0*(TTTSSS)

CRI (J)BLOQM(J)+2. 0*(TTT+SSS)

HIl (J)=YLOQP (J)+UU*(VVV+UJU)

RI I (U>=YLOQM(J)+UU*(VVVUUU)

C

CHN(NJ) = CHI(J)

C

CRN(NJ) = CRI(J)

NJ=NON+J

CHI (NJ)TPI*(4WWRRR)

CRI (NJ)=TPI*(WWW+RRR)

HTI (J)=W1*(000PPP)

RTI (J)=W1*(GGG+PPP)

C

CHN(J> = CHI (NJ)

C

CRN(J) = CRI(NJ)

IF (J. EQ. NON) GO TO 110

EJI=EJT

CR Q=CR T

SRI=SRT

CLI=CLT

SLI=SLT

RARI=RART

RBRI=RBRT

RALI=RALT

RBLI=RBLT

CZRI=CZRT

SZR I=SZRT

CZLI=CZLT

SAS=0. O

RARO. O

CCAO. O

HVH=O. O SVS=o. O RVR=O. O CCV=o. o C9

C9

SLIGHT INCREASE IN SPEED IF THE FINAL INTEGRATION

_{AVOIDS THE}

C9

INTERIOR SURFACE SEGMENTS..

C9

NI = NON - NUL

C9

DO *** I=1,NI

C9

DO 170 1=1, NON

CALL SONG (HEAVI(1, I), HEAVT(11 I), ROLLI(1, I), ROLLT(1 I), I)

170 CONTINUE

C6

C6

FORCE COEFFICIENTS.

C6

_{FORCE IS THAT WHICH MUST BE APPLIEP TO THE CYLINDER (PER UNIT}

C6

_{LENGTH) TO CAUSE SINUSOIDAL OSCILLATIONS AT THE GIVEN FREQUENCY}

C6

AND UNIT AMPLITUDE.

C6

_{COEFFICIENTS ARE THE PARTIAL DERIVATIVES OF THE}

_{FORCE BY THE}

C6

_{ACCELERATION OR VELOCITY COMPONENT OF THE GIVEN}

_MOTION.

C6

C6

ACCELERATION TERMS..

C6 HAH=HAH*R H02

SASSAS*RHO2

RAR=RAR*RHO2

CCACCA*RHO2/2. O

C6

C6

VELOCITY TERMS..

C6

HVH=HVH*RSIG

SVS=SVS*RSIG

RVRRVR*RSIG

CCV=CCV*RSIG/2. O C6

RETURN

C3

leo FORMAT (36H *** HEAVE MATRIX, FREQUENCY

_{INDEX =, 13, 15H, WAVE NUMEE}

IR =,IPE13.5/31H *** COEFFS. FOR

_{STATION AT X =,OPF13.5)}

190 FORMAT (40H *** SWAY-ROLL MATRIX,

_{FREQUENCY INDEX =, 13, 15H, WAVE N}

lUMBER

_{, 1PE13. 5/31H *** COEFFS. FOR STATION AT X}

_{, OPFI3, 5)}

200 FORMAT (26H *** EXECUTION TERMINATED.

C

END

SUBROUTINE WOMEN

i

_{(I,BLOOP,YLQQP,BLOQM,YLOQM,CHI,CRI,HII,HTL,RII,RTI)}

C

REAL BLOQP(I)1

_{BLOQM(l)YLDCP(1)YLOQM(l)1CHI(1)HII(1),RII(1)CRI(}

11) NIl (1), Rh (1)

Cs

COtIMON/GEOMETRY/MSTA, LPTS(25), YOFF(25, 25), NAFTA XAFT(25),

* _{YAFT(25), NFWDS XFWD(25), YFWD(25), XOFF(25),}

*

_{ZOFF(25,25),XFPERP,XAPERP,SHIPLFS'PBSSHIPT,}

*

Y1(21,25),ZWL(25)1WL(25),INPTS(25),Xt.JLF,XWLA,XXF,

* _{XXAS TAN, NON, NOEI NWL CR, XXFWDI XXSTA, XXAFTI}

DX, Dxl,

*

DX21Z2(21)1Y2(21),ZZ(20),YY(2Q),SNE(20),CSE(20),

* _{DEL(20), ROL(20), ADUUST1 UMAX, YMAX1 ZMAX1}

AREA, VERT

C

160 HAH=0. O C6

CS

CS

INPUT AND OUTPUT LOGICAL uNiTS..

COMMON/ID/INPUT9 OUTPUT. BIF, OFF, COF INTEGER OUTPUT, BhF, OFF. COF

C5

C5

COMMDN/QEDMETRY/MSTA, LPTS(25), YOFF(25, 25), NAFT, XAFT(25),

*

YAFT(25),NFWD,XFWD(25).YFWD(25),XOFF(25),

* ZOFF(25, 25), XFPERP, XAPERP, SHIPL, SHIPB, SHIPT,

* Vi (21, 25). ZWL(25). WL(25), INPTS(25), XWLF, XWLAI XXFJ

* XXA, TAN, NON, NUE, NWL, CR, XXFWDI XXSTAI XXAFT, DX, DX1,

*

DX2,Z2(21),Y2(21),ZZ(20),YY(20),SNE(20),CSE(20),

* DEL(20), ROL(20), ADJUST, L4MAX, YMAX, ZMAX, AREA, VERT

COMMON /1 NAH. SAS, RAR, CCA, HVH1 SVS, RVR, CCV, RHO2, RSIQ. WN, Wi, W2.ERR. Z 1RI. YRI, EJT

COMMON /

/

_{HAO(24,8),BLOGP(20,20),YLDGP(20,20),BLOQM(20}

1,20). YLOQM(20, 20). CONH(40), CONR(40, 2>, CROLL(40, 40), CHEAV(40, 40), HE 2AVI(20. 20), HEAVT(20, 20), ROLLI(20, 20), ROLLT(20, 20), EJI, CZRI, CZLI, SZ 3RI SZLI, RARI, RALI, RBRI, RBLI, CLI. CRI, SLI, SRI, I, IPESO, NI

DO 110 I=1,NON

NI=NON+I

CONH(I)=0. O CONR(I, 1)=0. O CONR(I, 2)=0. O CONH(NI )=CSE( I) CONR(NI, 1)=-SNE(I) CONR(NI, 2)=ROL(I)

ZRIWN*ZZ( I)

YRI=-WN*(YY( I )+Y2(1)) EJT=EXP (-YR I)

EJI=EJT

CALL ROMEO (CZRI, SZRI, RARI1 RBRI, CRI, SRI) C ZL I =C ZR I

SZLISZR I

RALIRAR I

RBLI=RI3RI

CLI=CRI

SLISR I

CALL WOMEN (I, BLOOP (1, I) YLDOP (1. I), BLIJOM( 1, I), YLOQM( 1, 1). CHEAV 1

(i,I),CROLL(1,I),HEAVI(1,I),HEAVT(1,I),ROLLI(1,I),ROLLT(1,I))

110 CONTINUE

IF (NWL. EQ. 0) GO TO 130 I =NOE-NWL+ i

DO 120 I=I,NOE

CONH( I )0. O CONR(I, 1)O. O CONR(I, 2)=0. O

120 CONTINUE

130 IT=JULIET(1, CREAVI CONk) IF (IT. EQ. 0) 00 TO 140

WRITE (OUTPUT, 180) K,WN,XXSTA

IF (IT.NE. 1) GO TO 150

140 IT=JULIET(2. CROLL, CONR) IF (IT. EQ. 0) 00 TO 160

WRITE (OUTPUT, 190) K, WN, XXSTA IF (IT. EQ. 1) 00 TO 160

150 WRITE (OUTPUT, 200)

STOP 11

APRT=ATAN2(YMT, ZMT) IF (ZNT. 0E. 0. 0) GO TO 130

IF (J1.GT.I) GO TO 110

IF ('t'MT. LT. 0. 0) APRT=APRT+TPI

GO TO 120

110 IF ('t'MT. GE. 0.0) APRT=APRTTPI 120 IF ('(PT. LT. 0. 0) 00 TO 130

ACRTPIN

GO TO 140

130 ACRT=ATAN2(YPT, ZMT) 140

ACLT=ATAN2(VPT, ZPT)

APLT=ATAN2(VMT1 ZPT) FPRT=ALOG ( ZMT*ZMT+YMT*YMT) /2. 0 FPLT=ALOG ( ZPT*ZPT+YMT*VMT) /2. 0 FCRT=ALOG( ZMT*ZMT+VPT*YPT) /2. 0 FCLT=ALOQ C ZPT*ZPT+YPT*VPT) /2. 0

SIMJ=SNE( I )*CSE(J)SNE(J)*CSE( I)

CIMJ=CSE( I )*CSE(J)+SNE( I )*5NE(J)

SIPJ=SNE( I )*CSE(J)+SNE(J)*CSE I

s

CIPJ=CSE( I )*CSE(J)SNE( I )*5NE(J)

DPNR=SIFIJ*(FPRIFPRfl+CIMJ*(APRIAPRT)

PPR=CSE(J)*(ZMI*FPRIYMI*APRIZMIZMT*FPRT+YMT*APRT+ZMT)+SNE

(J)*(YMI*FPRI+ZMI*APRIYMIYMT*FPRTZMT*APRT+YMT)

DPNL=SIPJ*(FPLTFPLI)+CIPU*(APLTAPLI)

PPLCSE(U)*CZPT*FPLTYMT*APLTZPTZPI*FPLI+YMI*APLI+ZPI)+SNE

(J)*('(MI*FPLIZpI*APLI+VMTYMT*FPLTZPT*APLTYMI)

DCNR=SIPJ* (FCR IFCRT) +C IPJ* (ACRIACRI)

PCR=CSECJ)*(ZMI*FCRIYPI*ACRIZMIZrIT*FCRT+VPT*ACRT+ZNT)+SNE

j

(J)

*

(YPT*FCRT+ZMT*ACRT+YP IVP I*FCRIZMI*ACR 1YPT)

DCNL=SIMJ*(FCLTFCLI)+CIMJ*(ACLTACLI)

PCL=CSE(J)*(ZPT*FCLTYPT*ACLTZPTZPI*FCLI+YPIsACLI+ZPI)+SNE

(J)*(YPT*FCLT+ZPT*ACLTVPTypI*FCLIZpI*ACLI+ypI)

BLOGP(JI I )=DPNR+DPNLDCNRDCNL

YLOGP(J, I )=PPR+PPLPCRFCL

BLOGM(J, I )=DPNRDPNLDCNR+DCNL

VLOGrl(J, I )=PPRPPLPCR+PCL

IF (J. EQ. NON) GO TO 150

FPRI=FPRT

FPLI=FPLT

FCRI=FCRT

FCLI=FCLT

APRI=APRT

APLI=APLT

ACR I=ACRT

ACLI=ACLT

ZMI=ZMT

YMI=YMT

ZP I=ZPT VP I=YPT 150

CONTINUE

160 CONTINUE

RETURN

C

END

StJBROUTINE WINE (K)

C

Cl

TWODIMENSIONAL HYDRODYNAMIC CALCULATION FOR NONZERO FREQUENCIES.

cl

C6

THIS SUBROUTINE IS CALLED FOR EACH STATION AND ALL NONZERO

co

FREQUENCIES WHEN THE HYDRODYNAMIC COEFFICIENTS ARE BEING

C6

C

END

SUBROUTINE GIRL

C

Ci

CALCULATION OF FREQUENCY INDEPENDENT TERMS TO BE USED IN THE

Cl

TWO-DIMENSIONAL HYDRODYNAMIC CALCULATIONS.

Ci

C6

THIS SUBROUTINE IS CALLED ONCE FOR EACH STATION OF THE SHIP

Có

WHEN THE HYDRODYNAMIC COEFFICIENTS ARE BEING GENERATED.

C6 C5

COMMON/QEOMETRY/MSTA. LPTS(25), YOFF(25, 25). NAFT. XAFT(25),

*

YAFT(25).NFWD,XFWD(25),YFWD(25),XQFF(25),

*

ZOFF(25,25), XFPERP. XAPERP,SHIPL,SHIPB,SHIPT,

*

_{Y1(21,25).ZWL(25).WL(25),INPTS(25).XWLF,XWLA,xxE,}

* XXA, TAN, NON. NOE. NWL, CR, XXFWDS XXSTA, XXAFT.DX, Dxl,

*

_{DX2.Z2(21)1Y2(21).ZZ(20),YY(20),SNE(20),CSE(20),}

* DEL(20), ROL(20). ADJUST. WMAX, YMAX. ZMAX, AREA, VERT

COMMON If HAl, SAl. RAi. CAl, HV1. SV1, RV1. Cvi, RHO2, RSIQ. WN. Wi, W2. ERR. Z 1RI. YRI. EJT

COMMON /

/

HAO(24,8).BLOGP(20,20),YLOGP(20,20).BLOQM(20,

120), YLOGM(20. 20), I. J1 ACLI, ACLTI ACRI. ACRT, APLI. APLT. APRI, APRT, CIMJ. 2CIPU, DCNL, DCNR, DPNL, DPNR, ECLI. FCLT, FCRI, FCRT, FPLI. FPLT. FPRI, FPRT. P

3CL, PCR. PPL, PPR, SIMJ. SIPJ, ZMI. ZMT, ZPI, ZPT, YMI, YMT, YPI, YPT

DATA PIN/-3. 14159265358979/

DATA TP 1/6. 28318530717958/

DO 160 1=1 NON

ZMI=ZZ (I) ZP I=ZMI YMIYY( I )-Y2( 1) YPI=YY( I )+Y2( 1)

FPRI=ALOGZMI*ZMI+YMI*YrlI)/2. O

FPLI=FPRI

FCRI=ALOG( ZMI*ZMI+YPI*YPI ) /2. 0 ECL I=FCR I

APRIATAN2(YMI, ZMI)

APLI=APR I

ACRI=ATAN2(YPI, ZMI)

ACLI=ACRI

DO 150 J=1,NON

J1=J+1

YMT=YY( I )-Y2(Ji) YPT=YY( I )+Y2(Ji) ZMT=ZZ( I )-Z2(J1) ZPT=ZZ( I )+Z2(Ji)

C

CALCULATE ANGLES (MEASURED OUTSIDE SECTION)..

220 CONTINUE

HAH=HAH*RHO2

SASSAS*RHO2

RAR=RAR*RHO2

CCA=CCA*RHO2/2. O

IF(WN. EQ. 0. ) HAH=99.

INFINITE AT ZERO FREG

RETURN

C

230 FORMAT (43H *** HEAVE MATRIX, ZERO ENCOUNTER FREQUENCY)

240 FORMAT (37H *** HEAVE MATRIX, INFINITE FREQUENCY)

20 FORMAT (47H *** SWAY-ROLL MATRIX, ZERO ENCOUNTER FREQUENCY)

260 FORMAT (41H *** SWAY-ROLL MATRIX, INFINITE FREQUENCY)

270 FORMAT (31H *** COEFFS. FOR STATION AT X =,F13.5)

280 FORMAT (26H *** EXECUTION TERMINATED.

ROLLT(t,J)=-YLOOM(U1 I)

HEAVT(I1J)=-YLOQP(U1 I)

140

CONTINUE

10 CONTINUE

C3

C3

SOLUTION FOR EITHER THE ZERO OR INFINITE FREQUENCY CASE..

C3

160 CONTINUE

DO 170 I=1NN

CONH( I )CSE( I)

CONR( I 1 )=-SNE( I)

CONR (I, 2)=ROL( I)

170 CONTINUE

IT=LNEGT(401 NN, 1, CREAVI CONH, ERR, HEAVI)

IF (IT. EQ. 0) QQ TO 180

IF (WN.EQ.0.0) WRITE (OUTPUT1 230)

IF (WN.NE.0.0) WRITE (OUTPUT124O)

WRITE (OUTPUT,270) XXSTA

IF (IT. NE. 0) GO TO 190

180 CONTINUE

IT=LNEQT(40, NN, 2 CROLL, CONR1 ERR1 ROLL I)

IF (IT. EQ. 0) GO TO 200

IF (WN.EQ.0.0) WRITE (OUTPUT125O)

IF (WN.NE.0.0) WRITE (OUTPUT126O)

WRITE (OUTPLIT,270) XXSTA

IF (IT. EQ. 1) GO TO 200 190 WRITE (OUTPIJT128O)

STOP 10

C3

C3

EVALUATE VELOCITY POTENTIALS AND FORCE COEFFICIENTS,.

C3

200 DO 220 I=1,NN

PAHO. O

PAS=0. O

PARO. O

DO 210 U=1NN

PAHPAH+CONI-1(U)*HEAVT(J, I) PAS=PAS+CONR(U, 1)*ROLLT(J, I) PAR=PAR+CONR (U, 2)*ROLLT(U, I) 210

CONTINUE

C6

Cá

_{THE PRESSURES IN PHASE WITH THE SINUSOIDAL DISPLACEMENT ARE..}

C6

C6

HEAVE -- PAR

=

PAH*RHO*ESIQ*ESIG

C6

SWAY

-- PAS

=

PAS*RHO*ESIG*ESIG

C6

ROLL

-- PAR

=

PAR*RHO*ESIQ*ESIG

C6

C6

THE ACCELERATION COMPONENTS OF THE FORCE ARE EQUAL

C6

_{IN MAGNITUDE TO THE ABOVE1 BUT HAVE THE OPPOSITE SIGN.}

C6 DDD=DEL ( I) DC I=CSE ( I) *DDD DSI=-SNE(I )*DDD DFR=ROL( I )*DDD C6

C6

INTEGRATION TO OBTAIN FORCE ACCELERATION COEFFICIENTS..

C6

HAH=HAR+PAH*DC I

SAS=SASPAS*DS I

R AR =R AR +P AR *DFR

CCA=0. O HVH=O. O SvS=0. o RVR=0. O cCv=0. o

NN=NON-NWL

C3

C C IF (WN. NE. 0. )

GO TO 130

C3

ZERO FREQUENCY CASE..

C3

DO 120 I=1,NN

XM1=ZZ (I )-Z2( 1) XP1=ZZ( I )+Z2( 1) YP1=YY( I )+y2( 1)

FCR1=. 5*ALOQ(XM1**2+Yp1**2)

FCL1=. 5*ALQQ( XP1**2+YP1**2) ACR1=ATAN2(YPI. XMl) ACL1=ATAN2(YP1, XP1)

DO 110 U=1,NN

XM2ZZ( I )-Z2(J+1) XP2=ZZ( I )+Z2(J+1) YP2YY( I )+Y2(J+1) FCR2=. 5*ALQQ ( XM2**2+YP2**2) FCL2=. 5*ALQQ ( XP2**2+YP2**2) ACR2=ATAN2 ( YP2, XM2)

ACL2=ATAN2(YP2, XP2)

SIMJ=SNE( I )*CSE(J)-SNE(J)*CSE( I) CIMJ=CSE( I )*CSE(J)+SNE( I )*SNE(J)

SIPU=SNE(I)*CSE(J)+SNE(.J)*CSE(I)

CIPJCSE(I)*CSE(U)-SNE(I)*SNE(J)

DCNR=SIPJ*(FCR1-FCR2)+CIPJ*(ACR1-ACR2)

PCR=CSE(J)*(XM1*FCR1-YP1*ACR1-XM1-XM2*FCR24-YP2*ACR2+XM2)+SNE

(J)*(yp2*FCR2+XM2*ACR2+Yp1-Yp1*FCR1-X*ACR1-yp2)

DCNL=SIMJ*(FCL2-FCL1 )+CIMU*(ACL2-ACL1)

PCL=CSE(U)*(XP2*FCL2-YP2*ACL2-XP2-XP1*FCL1+YP1*ACL1+XP1)+SNE

(J)*(YP2*FCL2+XP2*ACL2-YP2-YP1*FCL1-XP1*ACL1+YP1)

CROLL(IJ J)BLOQM(J, I)+2. O*(DCNR-DCNL)

CHEAV(I J)=BLOGP(J, I)+2. O*(DCNR+DCNL)

ROLLT(I, J)=-YLOGM(J, 1)-2. O*(PCR-PCL) HEAVT( I, J)=-YLOGP

(J1 1)-2.

O*(PCR-4-PCL)

IF (J. EQ. NN) GO TO 110 X M 1= X M2

XP1=XP2

YP1YP2

FCR 1=FCR2

FCL1=FCL2

ACR i ACR2

ACL1=ACL2

110

CONTINUE

120 CONTINUE

CO TO 160

C3

C3

INFINITE FREQUENCY CASE..

C3

130 CONTINUE

DO 150 I=1NN

DO 140 J=i.,NN

CROLL(I,J)=BLOQM(J, I)

CHEAV(II J)=BLOGP(J1 I)

C

INTERPOLATE FOR WATERLINE..

C

X0=XAFT C Ui)

XX (XXX0) / (Y0YNEXT) *yo+Xo

267

IF (XX.GE.XWLA) GO TO 269

XWLAXX

C

268

XXA=AMIN1(XX, XXA)

269

Y0=VNEXT

270 CONTINUE

290 IF (ISWL.EG.0) 1SWL26

C

420 FORMAT (32H0S T A T I O N

C E O M E T R Y)

480 FORMAT (22H0

DRAFT FWD (AT X =,F10.3,3H)

_{,F10.3/5X17HDRAFT AF}

lT (AT X = F10. 3 3H) =1 F10. 3)

570 FORMAT

(61H0*** TWO

DRAFTS SPECIFIED, AND LENGTH BETWEEN PERPE.

_IS

i

ZERO. /31H *** BOTH PERPENDICULARS AT X =F12.4/16H *** DRAFT

FWD

2=,F12.4/16H *** DRAFT AFT =,F12.4)

575 FORMAT

(12H0***

_{STATION, 13, 5H (X =1F12.41 18H) IS OUT OF ORDER. /29H}

i *** PREVIOUS STATION HAS X =F12.4)

580 FORMAT (25H0*** SHIP IS ABOVE WATER. )

C

END

SUBROUTINE ERROR(NO, IDUM,RDUM)

COMMON/IO/INPIJT, OUTPUT, BIF, 0FF. COF

INTEGER OUTPUT, BIF, OFF, COF

WRITE(OLJTPUT, 10) NO

WRITE(OUTPLIT, il) IDUM,RDUM

10

FORMAT(' STOPPED DUE TO ERROR NO.

',12,//)

11

FORMAT(1X. 13, SX,F10. 3)

STOP END

SUBROUTINE BEER (K)

C

Cl

_{TWODIMENSIONAL HYDRODYNAMIC CALCULATION FOR THE SPECIAL CASE}

Cl

OF ZERO OR INFINITE FREGUENCY.

Cl

C5

COMMON/GEOMETRY/MSTA, LPTS(25), YOFF(25, 25), NAFT, XAFT(25),

*

YAFT(25), NFWDI XFWD(25), YFWD(25), XOFF(25),

*

ZOFF(25, 25)1 XFPERP1XAPERP1SHIPL,SHIPB,SHIPT,

*

_{Y1(21125)1ZWL(25)..WL(25)1INPTS(25),XULF,XWLA,XxF,}

*

XXA1 TAN, NON, NOE, NUL, CR, XXFWDI XXTAI XXAFT, DX, DXI,

*

_{DX21Z2(21),Y2(21),ZZ(20)1YY(20),SNE(20),CSE(20),}

*

DEL(20),ROL(20), ADJUST, WMAXI '(MAX1 ZMAXI AREA, VERT

COMMON // HAH, SAS, RAR, OCA1 HVH,SVSI RVR, CCV, RHO2, RSIQ UN, Wi, W2,

ERRI Z

1RI, '(RI1 EUT

COMMON /

/

_{HAO(24,8),BLOGP(20,20),YLOQP(20,20),BLOGrI(20}

1,20), YLOGM(201 20), CONH(40), CONR(40, 2), CROLL(40, 40), CHEAV(40,

_{40) HE}

2AVI (20, 20), HEAVT(20, 20), ROLLI (20e 20), ROLLT(20, 20), EJI, CZRI,

_{CZLI1 SZ}

3R1, SZLI, RARI,

RALI, RBRI, RBLI, CLI, CRI, SLI, SRI, I

IPESO, J, NU

C5

C5

INPUT AND OUTPUT LOGICAL UNITS..

COMMON/I0/INPUT1 OUTPUT, ElF, OFF, COF

INTEGER OUTPUT, ElF, OFF, COF

CS

CS

_{OUTPUT LISTING PAGE HEADING DATA..}

C5

HAH=0. O

SAS=0. O RAR=0. O

C

C

FIND FORWARD AND AFTER ENDS OF WETTED HULL..

C 200 XWLA=1. 0E32

XWLF=-XWLA

IF (ISWL.EG.0) CO TO 210

XWLF=XOFF( ISUL) XWLA=XOFF C LSWL)

210 XXFXOFF(ISTA)

IF (NFWD. EQ. 0) CO TO 250

IF (NFWD.OT.1) GO TO 220

C

C

XFWD(1) DEFINED AS FORWARD END OF WATERLINE..

C

XWLF=AMAX1 (XFWD( 1) XWLF) XXF=AMAX1 (XWLF1

XXF)

GO TO 250

C

C

FIND FORWARD END

OF

WATERLINE AND FORWARD END OF WETTED HULL..

C

220 Y0=(XFWD(1 )-XFPERP)*TAN+YFWD(1)-TF

DO 230 U=21 NFWD

X X=XFWD C J)

YNEXT= CXX-XFPERP )*TAN+YFWD (J) -TF

IF (YNEXT. LT. 0. 0) CD TO 228 IF (YO. CT. 0. 0) 00 TO 229 IF (YNEXT. EQ. 0. 0) 00 TO 227

C

C

INTERPOLATE FOR WATERLINE..

C XO=XFWD(J-1.) XX=(XX-X0) / (YO-YNEXT)*Y0+XO

227

IF (XX. LE. XWLF) GO TO 229

XWLF=XX

C

228

XXF=AMAX1(XX1 XXF)

229

YO=YNEXT

230 CONTINUE

C

C

FIND AFTER END OF WATERLINE AND AFTER END OF WETTED HULL..

C

250 XXA=XDFF(LSTA)

IF (NAFT. EQ. 0) 00 TO 290

IF (NAFT.GT.1) GO TO 260

C

C

XAFT(1) DEFINED AS AFTER END OF WATER LINE..

C

XWLA=AMIN1(XAFT(1), XWLA)

XXA=AMIN1 (XWLA, XXA)

CO TO 290

C

C

FIND AFTER END OF WATERLINE AND AFTER END OF WETTED HULL..

C

260 YO=(XAFT(1 )-XFPERP)*TAN+YAFT(1)-TF

DO 270 J=2,NAFT

X X=XAFT C J)

YNEXT=(XX-XFPERP )*TAN+YAFT(J)-TF

IF (YNEXT. LT. 0. 0) 00 TO 268 IF (YO. CT. 0. 0) GD TO 268 IF (YNEXT. EQ. 0. 0) GO TO 267 C

(

L. C C C

XFPERPXOFF( 1)

XAP ERP = X OFF ( MSTA

IF (XFPERP.NE.XAPERP) QQ

TO 130

TAN=0. O IF (TF. EQ. TA) GO TO 140

WRITE(OUTPUT 570) XFPERP, 1F. TA

STOP 5

130 TAN(TATF)/(XFPERPXAPERP)

140 ISTAO

ISWLO

LSWLO

XXXOFF( 1)

DO 190 U=1,MSTA

IF (XX. QE. XOFF(J)) QQ TO 142

WRITE (OUTPUTs 575) J. XOFF(J)1XX

STOP 6

C

C

TU IS DRAFT OF STATION U.

C

142

XX=XOFF(J)

TU= ( XFPERPXX ) *TAN+TF

YNEXT=YOFF( 1 J)TU

N0

IF (YNEXT. GT. 0. 0) 00 TO 152

LSTA=J

IF (ISTA.EQ.0) ISTA=J

Y1(1. J)=YNEXT

N=LPTS(U)

DO 150 I=2.N

YNEXT=YOFF( I, J)TJ

IF (YNEXT.QT.0.0) GO TO 160

Yl (I, .J)=YNEXT 150

CONTINUE

I=N+1 IF (YNEXT. EQ. 0. 0) GO TO 170 C

C

SECTION IS NOT SURFACE PIERCING.

C

152

INPTS(J)=N

WL(J)=. FALSE.

GO TO 190

C

C

SECTION IS SURFACE PIERCING.

FIND WATERLINE COORDINATES..

C

160

YOYOFF( I-1, J)TU

IF (YO. EQ. 0. 0) GO TO 170

INPTS(J)=I-1

ZO=ZOFF( I-1, U)

ZWL(J)=Z0Y0*ZOFF(I,J)Z0)/(YNEXTY0)

GO TO 180

170

INPTS(J)=I-2

ZWL(J)=ZOFF(I-1. J) 180 WL(J)=. TRUE.

LSWLJ

IF (ISWL. EQ. 0) ISWL=J

190 CONTINUE

IF (ISlA. NE. 0) GO TO 200 WRITE (OUTPUT, 580)

STOP 77

(

C

C

C

C

CNEW(K)=C

250

SNEW(K»S

ZZ (I )=ZZNEW(N1) VY( I )YYNEW(N1)

ROL(I)=RNEW(N1)

DEL(I )D

NUTNON+ i

CALL INSERT (Z2, ZNEW(2), Il, NUT.

NUMBER)

CALL INSERT (Y2,YNEW(2), Il. NUT, NUMBER) CALL INSERT (ZZ, ZZNEWI I, NON, NUMBER) CALL INSERT (VV. YYNEW, I, NON. NUMBER)

CALL INSERT (DEL,DNEW, I,NON, NUMBER)

CALL INSERT (ROL, RNEW, L NON. NUMI3ER) CALL INSERT (CSE. CNEW, L NON, NUMBER) CALL INSERT (SNE,SNEW, I, NON1 NUMBER)

260

I=I1-fNUMBER

NON=NON+NUMBER

270 CONTINUE

280 NOE=NON+NON

RETURN

C C

ERROR DIAGNOSTICS..

C

C

TOO MANY WET SEGMENTS..

C

290 WRITE (OUTPL!T,300) XXSTA

STOP 7

C

300 FORMAT(' More than 20 wet segments for station at X='.F13.5)

C

END

SUBROUTINE FLOAT

C

C

THIS SUBROUTINE APPLIES THE GIVEN DRAFT TO THE ORIGINAL TABLE

C OF OFFSETS.

COMMON/IO/INPLJT, OUTPUT, BIF, OFF, COF

INTEGER OUTPUT, ElF, OFF. COF

C

COMMON/SHIP/ISTA. LSTA, ISWL, LSWLI 1F, TA, XCQ.YCO, DISPL

COMMON/DRFT12/ DRAFT(6, 2), IDRAFT

COMMON/GEOMETRY/MSTA. LPTS(25), VOFF(25, 25). NAFT, XAFT(25),

*

YAFT(25),NFWD.XFWD(25).YFWD(25).XOFF(25),

* ZOFF(25. 25), XFPERP. XAPERPI SHIPL, SHIPE, SHIPT,

*

\'i(21,25),ZWL(25).WL(25),INPTS(25)XWLF,XWLA.XXF,

* XXA. TAN, NON, NOE. NWL, CR, XXFWD, XXSTAI XXAFT, DX, DX1.

*

DX2,Z2(21).V2(21),ZZ(20),YY(20),SNE(20),CSE(20),

* DEL(20). ROL(20), ADJUST, WMAX, YMAX, ZMAX. AREA. VERT

LOGICAL UL, ADJUST

C

C

PLACE SHIP AT GIVEN DRAFT..

C C C TF=DRAFT( IDRAFT. 1) TA=DRAFT( IDRAFT, 2) C C C C

IF (XFPERP.NE.XAPERP) GO TO

130

C

r-CSE(K)=-1. O SNE(K)=O. O YY(K)=O. O ZZ (K)=Z2(K)-D ROL (K) =ZZ (K)

NON=K

K=K+1

Y2(K)=O.O

Z2(K)=Z2(K-1 )-ZINT

190 CONTINUE

Z2(K)=0. O NON =K- 1 C

C

END OF FIRST PASS.

ADD ADDITIONAL SEGMENTS IF REQUIRED..

C

200 IF (NON. OT. LIMIT) GO TO 290

IF (. NOT. ADJUST) GO TO 280 IF (NON. QT. MAXPTS) GO TO 280 IF (MTOT. EQ. 0) GD TO 280

MTOT=MTOT+NON

M1=NON-NWL

IF (MTOT. LE. MAXPTS) GO TO 230

C

C

DECREASE

MORE

UNTIL

MTOT

IS EQUAL

MAXPTS

C

210 DO 220 K=1,M1

IF (FIORE(K).LE.0) GO TO 220

MORE(K)=MORE(K)-1

MTQT=MTOT-1

IF (MTOT.LE.MAXPTS) GO TO 230

220 CONTINUE

00 TO 210

C

C

INSERT ADDITIONAL SEGMENTS AS INDICATED BY

MORE

C

230 1=1

DO 270 M=1,M1

11=1+1 NUMBER=MORE (M) IF (NUMI3ER. LE. 0) GO TO 260 Nl =NUMBER+ 1. Z0=Z2 C I) Y0=Y2 C I) ZINT=Z2( Il )-Z0 YINT=Y2( 11)-YO

ZINT=ZINT/FLOAT(2*N1)

YINT=YINT/FLOAT(2*N1)

D=DEL( I) ÌFLOAT(N1) CCSE C I) S=SNE( I)

DO 240 Ki.,N1

P(K-1 )*2

Q=P+1. O ZNEW (K) =Z0+P*ZINT YNEL.J ( K)=Y0+P*YINT ZZNEW(K ) =ZO+Q*ZINT YYNEW(K ) =YQ+Q*VINT

RNEW(K)=(CR-YYNEW(K) )*S-ZZNEW(K)*C

240

CONTINUE

DO 250 K=1, NUMBER

DNEW(K)=D

C

SUBROUTINE POTST

COMMON/BD/XPAN( 120) YPAN(120), ZPAN(120), AREA(120), ST(120),

*

_{ACN(120),ACNW(120),AN(120,3),E(120),P(120,6),PRFS(120),}

*

_{STOLD(12O)PX(12016)}

COMMON/BD2/XPT(150). YPT(150), ZPT(150), WRF(150). WRFR(j.50),

*

KR(150,4)

COMMONIA/NPAN1 NPT GEE, RHO, NKX, NKY, EYE1 DL TIM, UFWD

COMPLEX

AI B3 EYE DIMENSION XPSL(3, 4), XPSLR(31 4) PBB(120, 120) COMMON/PTST/ARE4(2001 4), X4(200, 4), Y4(2001 4), Z4(2001 4) *

,SEL(200,4)

DO 1500 U1,NPAN

ARE4(J, 4)-1. O

JT=4

IF(KK(.J,4).EG.0) JT=3

DO 1500 UU=1,JT

J21

IF(JJ. LT. UT) J2=JJ+1 KF=KK(31 JJ) KG=KK(J, U2)

X4(U,JJ)(XPT(KF)+XPT(KQ)+XPAN(J))/3. O

Y4(J,JJ)=(YPT(KF)+YPT(KQ)+YPAN(U))/3.O

Z4(J, UJ)=(ZPT(KF)+ZPT(KG)+ZPAN<J) )/3. O

AF=XPT(RF)-XPAN(U)

BF=YPT(KF)-YPAN(U)

CF=ZPT(KF)-ZPAN(J)

AG=XPT(KQ)-XPAN( J)

EG=YPT(KG)-YPAN(J)

CG=ZPT(KC )-ZPAN( J)

CALL SELF(AF, EF, CF, AQ, BG, CG, FEE)

SEL(J, UJ)=FEE

CR=AF*BG-BF*AG

AR=BF*CG-CF*BG

B R C F* AG-AF*C G

ARE4 (U, JJ ) =0. 5*SQRT (AR*AR+ER*BR+CR*CR)

1500

CONTINUE

DO 127 N=1,NPAN

DO 1277 MU=1,NPAN

1277 PBE(NJ, MJ)=0. 00 P(NU, 1)=O.00 P(NJ, 2)=0. 00 P(NJ1 3)=O. 00 P(NJ, 4)=0. 00 P(NJ, 5)=0. 00 P(NJ, 6)=0. 00

DO 123 NK114

ARNARE4(NJ NR)

IF(ARN. LT. 0. 0) GO TO 128 P 1=0. 0 P2=0. 0 P3=0. 0 P4=0. 0 P5=0. 0 P6=0. 0 X=X4(NJ, NR) Y=Y4(NJ, NR)

ZZ4(NJ, NK)

DO 138 MJ=1,NPAN

DO 138 MR1,4

XF=X4(MJ, MR) YF=Y4(MJ, MR) ZF=Z4(MU, MR) ARM=ARE4(MJ, MR)

IF(ARM.LT.0.00) GO TO 138

IF(NU. NE. MU) GO TO 140 IF(MK. NE. NR) GO TO 140

FRA=SEL(MU MK)/ARM

GD TO 1380

C

CONTINUE STANDARD PROCEDURE..

C

140 DEL(K)=ABS(Z0)

CSE(K)=-1. O

SNE(K)0. O

ZZ(K)=O. 5*Z0 VV (K) =VO ROL ( K) =ZZ (K)

NON=K

K=K1

Z2(K)=0. O

Y2(K)Y0

GO TO 200

C

C

ADD SEGMENT UP TO WATERLINE..

C

150 ZINT=ZW-Z0

YINT-Y0

IF (NUL. LT. 0) NWLO

IF (ZW. LE. 0.0) NWL=0

D=SGRT( ZINT*ZINT+VINT*YINT)

IF (D. EQ. 0. 0) 00 TO 170

IF (.NOT. ADJUST) GO TO 160

C

C

CODE INSTRUCTIONS FOR THE ADDITION OF POINTS..

C

NUMBER=MAXO(IABS(IFIX(ZINT/ZMAX))1 IABS(IFIX(YINT/YMAX)))

MORE (K) =NUMBER

MTOT=MTOT+NUMBER

C

C

CONTINUE STANDARD PROCEDURE..

- C

160 ZS=Z0+ZW

YS=Y0

AREA=AREA+YINT*ZS

VERTVERT-YINT*(Z0*(CR+VINT/1. 5)+ZW*(CR+YINT/3. 0))

CZ INT/D

S=YINT/D

CSE(K)=C

SNE(K)S

DEL(K)=D

ZZ(K)=0. 5*ZS

VV(K)=0. 5*5

ROL ( K ) (CR-VY ( K) ) *S-ZZ (K)

NON=K

KK+1

Z2(K)=ZW

Y2(K)=0.0

C

C

ADD DECK AT WATERLINE..

r

170 CONTINUE

IF (NUL. EQ. 0) 00 TO 200 Z I NT=ZW/FLOAT ( NUL) IF (UMAX. EQ.0. 0) 00 TO 180 IF (ZINT. LT. UMAX) GO TO 180 NWLIFI X (ZW/UMAX )+1 Z INT=ZW/FLOAT (NUL) 1BO D=ZINT*0. S DO 190 1=1, NUL

DEL(K)=ZINT

MTOT=O

K 1 Il =K

IF (Z(1).EQ.O.0) 11=2

IF (Il. Ql. NPTS) 00 TO 130 C

C

CALCULATION LOOP FOR SL)BMERQED OFFSET POINTS..

C

DO 120 I=I1,NPTS

ZINT=Z( I )-ZO YINT=Y( 1)-YO

D=SGRT( ZINT*ZINT+YINT*YINT)

IF (D. EQ. 0. 0) 00 TO 120

IF (.NOT.ADJUST) CO TO 110

C

C

CODE INSTRUCTIONS FOR THE ADDITION OF POINTS..

C

NUMBER=MAXO(IABS(IFIX(ZINTÍZMAX)), IABS(IFIX(YINT/YMAX)))

MORE (K) =NUMBER

MTOT=MTOT+NUMBER

C

C

CONTINUE STANDARD PROCEDURE..

C 110

ZS=Z0+Z(I)

YS=Y0+Y( I) AREA=AR EM-Y I NT* ZS

VERT=VERTYINT*(Z0*(Y0-CR+YINT/3. 0)+Z(I)*(Y0-CR+YINT/1. 5))

C=Z I NT ÍD

S=YINT/D

CSE(K)=C

SNE(K)=S

DEL C K ) =D ZZ(K)=0. 5*ZS YY(K)=0. 5*YS

ROL(K)=(CR-YY(K) )*S-ZZ(K)*C

Z2(K)=Z0

Y2(K)Y0

Z0=Z C I)

Y0=Y(I)

K=K+1

120 CONTINUE

C

C

END OF CALCULATION LOOP FOR SUEMERCED POINTS.

C

130 Z2(K)=Z0

Y2(K)Y0

NON=K-1

C

C ADD UPPERMOST SEGMENT.. C

IF

(NUL. NE.

0) 00 TO 150

C

SECTION IS SUEMEROED.

IF (ZO. EG. 0. 0) 00 TO 200

IF (.NOT. AD.JUST) 00 TO 140

C

C

CODE INSTRUCTIONS FOR THE ADDITION OF POINTS..

C

NUMEER=IABS( IFIX(Z0/ZMAX))

MORE (K) =NUMEER

MTOT=MTOT+NUMB ER

C

END

SUBROUTINE INSERT(Al, A2, Ji Li, L2)

C

_{Purpose: Inserts array A2 into arra4 Al at}

location Ji.

REAL

Al(i)1 A2(l)

IF (Li. LT. Ji) GO TO 120

M=Li+L2

I=L1

K=LlJl+i

DO 110 J=I,K

Al (M)=Al (I)

M=M1

1=11

110 CONTINUE

120 I=Ji1

DO 130 K=1,L2

M=K+I

A 1 (M ) =A2 ( K)

130 CONTINUE

RETURN

END

SUBROUTINE STATN (Z,Y,ZW,NPTS)

C

C

_{CALCULATION OF DATA CONCERNING STATION}

_GEOMETRY. C

_{REVISION OF OFFSETS FOR GOOD RESULTS}

MAY SE PERFORMED.

C

C

NON = NUMBER OF CALCULATED MIDPOINTS.

C

_{NUL = NUMBER OF WATERLINE MIDPOINTS.}

C Z

_{= HORIZONTAL COORDINATE OF SEGMENT ENDPOINT.}

C Y

_{= VERTICAL COORDINATE OF SEGMENT ENDPOINT.}

C ZZ

_{= HORIZONTAL COORDINATE OF SEGMENT MIDPOINT.}

C 's'Y

= VERTICAL COORDINATE OF SEGMENT MIDPOINT.

C

_{SNE = HORIZONTAL COMPONENT OF UNIT}

NORMAL TO SEGMENT.

C

_{CSE = VERTICAL COMPONENT OF UNIT}

_NORMAL.

C

_{DEL = LENGTH OF SEGMENT.}

C

_{ROL = MOMENT OF UNIT NORMAL ABOUT}

CENTER OF ROLL (CG).

C

REAL Y(l),Z(1)

COMMON/ID/INPUT1 OUTPUT1 51F, OFF, COF INTEGER OUTPUT, BIF, OFF, COF

COMMON/GEOMETRY/MSTA,LPTS(25), YOFF(25, 25),

NAFTI XAFT(25),

* _{YAFT(25), NFWD, XFWD(25), YFWD(25),}

XOFF(25),

* _{ZOFF(25, 25)1 XFPERPI XAPERP.SHIPL,SHIPB,SHIPT,}

*

Y1(21,25)1ZWL(25),WL(25),INPTS(25)1XWLF,xWLA,XXF,

* _{XXAI TAN, NON1 NOE, NUL, CR, XXFWD, XXSTAS}

XXAFTI DX, Dxl,

*

DX2,Z2(21),Y2(21),ZZ(20),YY(20),SNE(20),CSE(20),

* _{DEL(20)I ROL(20)1 ADJUST, UMAX, YMAX,}

ZMAXS AREA, VERT

COMMON /

_{/ HAl, SAI, RAi, CAl, HVI, SV1,}

_{RV1, Cvi,}

1

_{RHO2, RSIG, UN,}

_{WI, W2,}

ERR,

_{XRII YRI, EUT}

COMMON /

/

_{HAO(24), SAO(24), RAO(24), CAO(24),}

i

_{HVO(24), SVO(24), RVO(24)1 CVO(24)}

COMMON 1/ ZNEW(20)s YNEW(20)1

_{ZZNEW(20), YYNEW (20), CNEW(20), SNEW(20),}

1DNEW(20)1 RNEW(20), MORE(20)1 I K1 ZO, YO1 Ml, MTOT, ZINT, YINT, NUMBER1 _{M1 Nl}

2, I1 JOB, ZS, YS, D C1 S, P,G, NUT

LOGICAL POJUST

DATA LIM1T/20/

DATA MAXPTS/20/

AREA=0. O VERT=0. O Z0=0. O YO=Y( 1)

I'-S-.

NFWD=0

NAFT=0

C

C *** CARD TYPE D

30 N=1 MSTA=MSTA+ i

IF

(MSTA. GT. 25) CALL ERROR(10 IDUM, RDUM) READ (OFF, 416) STATNOS Yll, Zi, JTEST

WRITE(OUTPUT 417) STATNO, Yll, Z1 JTEST

XOFF (MSTA ) =STATNO*SPACE

GO TO 50

40 CONTINUE

1oop within each station

N=N+1

IF

(N. GT. 25) CALL ERROR(11, MSTA, RDUM)

READ COFF,416) S,Y11.Z1,JTEST

WRITE(OUTPUT, 417) S Vil, Zi UTEST

IF (S. NE. STATNO) CALL ERROR( 12, MSTAI RDUM)

50 YOFF(N1 MSTA)=Z1*ZSCAL

ZOFF(N1 MSTA)=Y11*YSCAL

IF

(.JTEST. EQ. 0 . OR. JTEST. EQ. 77777) GO TO 40

LPTS(MSTA)N

!No.

of points- MSTA

IF

(N. LT. 2) CALL ERROR(131 MSTAI RDUM)

IF (.JTEST. EQ. 88888) GO TO 30

!Qo onto next station

IF (JTEST. NE. 99999) CALL ERROR(14 UTEST RDUM)

C C

DO 220 U=1MSTA

XDFF(J)=XOFF(J)

220 CONTINUE

IF (NFWD. EQ. 0) GO TO 240 X=XOFF( 1.)

DO 230 L=1,NFWD

XFWD( I )XFWD( I )+X

230 CONTINUE

240 CONTINUE

X=XOFF(MSTA)

DO 250 I=1NAFT

XAFT( I )XXAFT( I)

250 CONTINUE

RETURN

180 FORMAT (5X, 15)

190 FORMAT (15)

197 FORrIAT(1HI/,81(1H*)/, '

INPUT DATA ECHO ',T64,

*'PROQRAM HYDREX'/,Bl(IH*)//,33(1H),

*'EBIF] DATA FILE',32(1H)/)

198 FORMAT(1H1/,81(IH*)/, '

INPUT DATA ECHO ',T64,

*'PROQRAM HYDREX'11 81(1H*)//, 33(1H),

196 199

* '[OFF] DATA FILE', 32( 1H)!)

FORMAT(1X,A)

FORMAT(A)

200 FORMAT (6F10.2)

201 FORMAT (3Fb. 2, 15)

210 FORMAT (F10.2, IS, 5X,Fi0.2)

C

410 FORMAT (A)

412 FORMAT (4F10.3,13XI24X1I1)

414 FORMAT (SX,I5,5X. 'INPUT OF

SHCP TYPE D OFFSET DATA')

416 FORMAT (F6. 3 2F7. 0, 16)

417 FORMAT (F7.3,2F10.21I6)

420 FORMAT (215,F10.2)

430 FORMAT (2Fb. 2)

'S.

r

_S.. C C.

*

DX2,Z2(21),,V2(21),ZZ(20),YY(20),SNE(20),CSE(20),

*

DEL(20), ROL(20) AD.JUST. WMAXS YMAX, ZtIAX1 AREA, VERI

LOQICAL ADJUSTIWL

C

CHARACTER*81 CARDID

DATA

DEOREE/0.01745 32925 19943/

DATA

0

/32. 17/

DATA NWL/1/

C

WRITE(OUTPUT1 197)

C *** TITLE

READ (B IF1 199) TITLE

WRITE(OUTPUT. 196) TITLE

C *** DRAFT (fwd),DRAFT (aft)

long.

icc's of DRAFT marks

READ (ElF1 200) TE, TA, XFPERPP XAPERP

WRITE(OUTPUT1 200) IF, TA, XFPERPI XAPERP

C *** Center of Qravitq (XCG aft of FP, YCO above BL)

READ (EIF,200) XCQ,YCQ,ZCQ

WRITE(OUTP(JT, 200) XCQ, YCO, ZCC

C *** Six DRAFTS at which hudro. coeffs are computed

READ(B1F1200) (DRAFT(I11),I16)

WRITE(OUTPIjT, 200) (DRAFT( I, 1), 1=1,6)

C *** Minimum segment lengths for Frank Close Fit

READ (BIF,201) YMAXIZMAXIWMAXINWL

WRITE(OUTPUT, 201) YMAX, ZMAXI WMAX, NWL

ADJUST=ZMAX. 01. 0. 0. AND. YMAX. 01. 0. 0

C *** Number of forware profile points

READ (RIF1 190) NFWD

WRITE(OLJTPUT, 190) NFWD

IF (NFWD. CT. 25) CALL ERROR( 15, IDUM, RDUM)

C *** Coordinates of forward profile points

IF (NFWD. CT. 0) READ (RIF, 430) (YFWD( I), XFWD(I)1 1=1, NFWD)

WRITE(OtJTPIJT,430) (VFWD(I),XFWD(I),I=1.NFWD)

C *** Number of aft profile points

READ (B IF4 190) NAFT

WRITE(OUTPUT1 190) NAFT

IF (NAFT. CT. 25) CALL ERROR(161 IDUM, RDUM)

C *** Coordinates of aft profile points

IF (NAFT. Ql. 0) READ (B IF, 430) (YAFT( I), XAFT( I), 1=1, NAFT)

WRITE(OUTPUT, 430) (YAFT(I), XAFT(I)1 11, NAFT)

C C C C

C

Section 2.0 - READ OFFSET file

C

The offset file can be an actual SHCP DATA File

C

WRITE(OUTPUT, 198)

C

C *** CARD TYPE A

READ (OFF,410) CARDID

WRITE(OUTPUT, 410) CARDID

C *** CARD TYPE B

READ (OFF,410)

C *** CARD TYPE C

READ (OFF, 412) SPACES ZSCAL, YSCAL, SHIPL, NAPNI KINDO

WRITE(OUTPUT, 412)SPACE5 ZSCAL, YSCAL, SHIPL, NAPNI KINDO

IF (SPACE. EQ. 0. 0) SPACE=1.0

ZSCAL=1. O

YSCAL1. O

*

Tb, 'Vert. Moment

'T40E13.7,T60, '

units ',1)

292

FORMAT(T1O

'

METACENTRIC HEIGHTS

'/

*

Tb, 'BM (longitudinal)

'T4O. F13. 2, T60, '

units ',/

* TiO, 'BM (transverse)

'T4O1F13.2,T60, '

units ',/

*

TiO, 'GM (longitudinal)

'1T40F13.2,T60, '

units '.1

* TiO, 'GM (transverse)

'T40,F13.2,T60, '

units '.1)

293

FORMAT(T10. '

HYDROSTATIC FORCES

*

Tb, 'Roll Restoring Moment

'.T401E13.7,T60, '

units ',/

*

Tb, 'Pitch Restoring Moment

'T40Eb3.7,T6O, '

units

,/

*

Tb, 'Heave Restoring Force

'1T4OE13.7,T60, '

units

'si

*

Tb, 'Pitch Induced FReave Force

'1T40,E13.7,T60, '

units ',1)

294

FORMAT(/,T55 '

LWL begins at ',F10.2,/

* ,T5, ' LWL ends at

',F10.2,/)

300 FORMAT( 1H11, 81 ( 1H=), I lx, A30, 'ADDED MASS/DAMPING COEFFICIENTS'

*,T66, 'PROGRAM HYDREX'/,Bl(bH=))

301 FORMAT (T36, 'Station '12/,T36, '

'I,

* T5, '

Dist. from F.P.

'.F8.2,T50, 'Area ', F11. 3,!

* T5, ' DRAFT (fwd)

',F9.2,T50, 'Roll Ctr abv WL'.F11.3/

* T5, ' DRAFT (aft)

',F8.2.!1)

302 FORMAT(1OX, '

HEAVE----',

*4X, '

SWAY

',4X, ' ROLL.

*4X, '--SWAYROLL---', f. 14X, 'A22', 5X, 'B22'.

*7X, 'A33',5X, '533'1 GX, 'A44',ÓX. 'B44',6X, 'A34',6X, '534',/

*,2X, 'Freg. '1)

309 FORMAT ((lXi 0PF22. 4, 3(5X, 1P2E1O. 2)))

310 FORMAT(1X, F5. 2, 3X. F8. 4, F8. 4, 2X. F8. 4 F8. 4,

* 2X, F9. 2 F9. 2 F9. 1, F9. 1)

C

420 FORMAT (1H1/, 33(1H*), ' N O T E S ',33(1H*)/)

430 FORMAT(T5, 'Generate additional offset points. '1,

*

T5, 'Maximum segment heiqht=',F1O.3/,

*

T5, 'Maximum segment width

', F10. 3/)

440 FORMAT(/T5, 'Use segments as defined bij table of offsets. ')

450 FORMAT(T51 'No internal freesurface segments are used. ')

460 FORMAT(T5, 13,2X, 'nodes for internal freesurface. ')

470 FORMAT(SX, 'Add internal nodes if surface segment lengths',

*1

_{exceed',F10.3)}

706 FORMAT (6X, 'Height above', 5X. 'Half', 15X,

*'Submerged Offsets'!, lOX, 'Baseline',5X, 'BREADth',

.*14X, 'Y',9X. 'Z'!)

705 FORMAT (f/I/I, 33(1H), 'STATION OFFSETS', 32(1H), I)

710 FORMAT (6X,2F12. 3, F18. 3,F12. 3) 720 FORMAT (39X, 9HWATERLINE5 F12. 3)

730 FORMAT (6X,2F12.3)

END

SUBROUTINE INDATA

COMMON/lO/INPUT, OUTPUT, ElF1 OFF, COF INTEGER OUTPUT, ElF, OFF, COF

COMMON /SIGMA /

NV.,

SIGMA(24), SIGFIAO, ERRO, OM(12)

COMMON / DRFT12/ DRAFT(6,2),IDRAFT

COMMON/IOFILE/ OFFIL,BIFIL,COFIL

COMMON/HEAD/TI TLE

CHARACTER*30 TITLE

COMMON ¡U

/

RHO, Q

COMMON/SHIP/ISTA, LSTA, ISWL1 LSWL, TF5 TA, XCG, YCG, DISPL

COMMONIOEOMETRY/MSTA1 LPTS(25). YOFF(25. 25), NAFT, XAFT(25),

*

YAFT(25),NFWD.XFWD(25),YFWD(25),XOFF(25),

*

ZOFF(25,25),XFPERP,XAPERP,SHIPL,SI-11P5,SHIPT,

*

Y1(21,25),ZWL(25),WL(25),INPTS(25),XWLF,XWLA,XXF,

IF (WL(J)) WRITE (OUTPUT1 720) ZWL(J)

670

M=M+1

IF (M. LE. N) WRITE (OUTPUT1 730) (YOFF(I. J), ZOFF(I, U), I=M, N)

C

220

CONTINUE

240 CONTINUE

QAMMA=HO*0

D ISP L=VOLO*GAMMA C

WPT=WPT/3. o

YFYTOWP 0*QAMMA YFZRO=WP 1*QAMMA

XUAO. O

IF (WPO. NE. 0. 0) XWA=WP1/WPO

XCB=VDL1 /VOLO YCB=VDLVI VOLO XBM=WPT/ VOLO ZBM= (WP2WP 1*XWA) /VOLO X GM=X B M+YC B ZQM=ZBM+YCB XWA=XWA+XCG XCB=XCB+XCG YCB=YCB+YCC XMXRO=XOM*DISPL ZMZRO=ZGM*DISPL

WRITE (OUTPUT, 270) TITLE

WRITE (OUTPUTS 281) DRAFT(IDRAFT, 1). DRAFT(IDRAFT1 2)

WRITE(OUTPUT, 290) WPO, XWA1 WP1, WP2, WPT

WRITE(OUTPUT, 291) DISPL, VOLO1 XCB1 YCB, VOLl, VOLV

WRITE(OUTPUT, 292) ZEN, XBMP ZOM, XQM

WRITE(OUTPUT, 293) XMXRO. ZtIZRO, YFYTO, YFZRO

WRITE(OUTPUT, 294) XXF, XXA

WRITE (OUTPUT,280) XCC,YCG

C *** OUTPUT OF NOTES

WRITE(OUTPLJT, 420)

IF (ADJUST) WRITE( OUTPUT,430) ZtIAXIYMAX

IF (.NOT.AD%JUST) WRITE (OUTPUT144O)

IF (NUL. CT. O) GO TO 110

CO TO 120

110 WRITE (OUTPUT, 460) NWL

IF (UMAX. QT.O. O) WRITE (OUTPUT147O) WMAX

120 CONTINUE

C

RETURN

270 FORMAT (1H1/1 Bi (IH=) / iX, A30, 'HYDROSTATIC COEFFICIENTS'

*, 164, 'PROGRAM HYDREX', /,81(1H)/)

280 FORMAT (//S1(1H)f, ' NOTE: All moments are about',

*'center of gravit' /,7X, 'XCG=',F1i.3, ' YCG=',F11.3,/)

281 FORMAT(T5, ' DRAFT (f wd) =',F8.2,/,T5, ¡ DRAFT (aft) ='F8.2/)

290 FORFIAT(T10, '

_*

WATERPLANE

'1

110, 'Area

',T40, F13. 2. T60, '

units '/

*

TiO, 'LCF

_{',T40,F13. 2,160,' units ',/}

*

Tb, 'ist Long. Moment

',140, E13. 7,160, '

units

l

_f

*

_{110, '2nd Long. Moment}

_{',T40.E13.7,T60, '}

_{units "f}

*

Tb, '2nd Transy. Moment

'.T4O.E137,T60, '

units '.1)

291

FORMAT(T10, '

VOLUME

'.1

*

110, 'Displacement

'1T4O,'13.2,T6O, '

units '/

*

_{110, 'Volume of Displacement}

_{',T401E13. 7 160, '}

_{units ',1}

*

TiO, 'LCB

',T40, F13. 2, 160, '

units

'

/

*

110, 'VCB

_{',T40.F13.2.T6O, ' units ',1}