Program Tranship; A computerprogram for the design of the midship section of a transversely framed dry cargo ship Part two

(1)

SSC-2 16

1971

This document has been approved

for public release and sale;

its

distribution is unlimited.

SHIP STRUCTURE COMMITTEE

PROGRAM "TRANSHIP"

A COMPUTER PROGRAM FOR THE DESIGN OF

THE MIDSHIP SECTION OF A

TRANSVERSELY-FRAMED DRY CARGO SHIP

PART TWO

(2)

SHIP STRUCTURE COMMITTEE

AN INTERAGENCY ADVISORY COMMITTEE DEDICATED TO 1MPROVING

THE STRUCTURE OF SHIPS

SR- 175

1972

Dear Sir:

The Ship Structure Committee has undertaken a series of

research projects to develop analytical methods and computer

programs which will apply modern high-speed electronic

com-putational techniques to ship hull structures.

This report contains the details of the computer program

discussed in SSC-215, A Guide for the Synthesis of Ship

Structure--Part One--The Mid-Ship Hold of a Transversely Framed Dry Cargo

Ship.

Comments on this report would be welcomed. Sincerely,

W. F. REA, III

Rear Admiral, U. S. Coast Guard

Chairman, Ship Structure Committee

MEMBER AGENCIES: ADDRESS CORRESPONDENCE TO:

Nl TED STATES COAST GUARD SECRETARY

NAVAl SHIP SYSTEMS COMMAND SHIP STRUCTURE COMMITTEE MII IIARY SEALIFT COMMAND U.S. COAST GUARD HEADQUARTERS

MARITIME ADMINISTRATION WASHINGTON. D.C. 20591

(3)

SSC-216

Final Technical Report

on

Project SR-175, "Rational Ship Structural Design"

PROGRAM 'TRANSHIP'

A COMPUTER PROGRAM FOR THE DESIGN OF THE

MIDSHIP SECTION OF A TRANSVERSELY-FRAMED

DRY CARGO SHIP PART TWO

by

Manley St. Denis University of Hawaii

under

Department of the Navy Naval Ship Engineering Center Contract No. N00024-68-C-54O3

This document has been approved for public release and sale; its distribution is unlimited.

U.S. Coast Guard Headquarters Washington, D.C.

(4)

ABSTRACT

This report presents the computer program corresponding to the

method of design expounded in the Ship Structure Committee Report

SSC-215, A Guide for the Synthesis of Ship Structures - Part One - The

Midship Hold of a Transversely-Framed, Dry Cargo Ship. The program

consists in an executive routine, called TRANSHIP, and twenty seven

(5)

INDEX Page No. Introduction i Description of Input 2 Program TRANSHIP Description 3

Notation - Mathematical Symbols to FORTRAN 4

FORTRAN Symbols to Mathematical 7

Listing 15

Macro Flow Diagram

17

Detailed Flow Diagram 19

Subroutine ASPECT

Description and Notation 24

Listing 25

Flow Diagram 26

Subroutine BNMAT

Listing 30

Flow Diagram 30

Subroutine COSTING

Listing 32

Flow Diagram 33

Subroutine DOCUMENT

Listing 37

Flow Diagram 38

Subroutine EV

Listing 40

Flow Diagram 40

Subroutine FRAME

Listing 43

Flow Diagram 46

Subroutine FSHAPE

Listing 57

Flow Diagram 57

(6)

-111-INDEX (Cont'd)

Page No. Subroutine GEOM

Listing 59

Flow Diagram 60

Subroutine GRILLAGE

Listing 62

Flow Diagram 64

Subroutine INPUT

Listing 69

Flow Diagram 71

Subroutine INTERMED

Listing 75

Flow Diagram 76

Subroutine LONGIT

Listing 78

Flow

Diagram 79

Subroutine LONGMAT

Listing 81

Flow Diagram 82

Subroutine MATINV

Listing 86

Flow

Diagram 87

Subroutine NTPLATE

Listing 90

Flow

Diagram 91

Subroutine PLATING

Listing 94

Flow

Diagram 96

Subroutine REACT

Listing 102

Flow

Diagram 102

(7)

-iV-Subroutine RTPLSUB

Listing 104

Flow Diagram 105

Subroutine SECTION

Listing 111

Flow Diagram 112

Subroutine SHAPES

Listing 115

Flow Diagram 115

Subroutine SOLVE

Listing 118

Flow Diagram 120

Subroutine STEP

Listing 126

Flow Diagram 126

Subroutine STRESS

Listing 129

Flow Diagram 130

Function T

Listing 135

Flow Diagram 135

Function THETA

Listing 137

Flow Diagram 137

Subroutine TRANSV

Listing 139

Flow Diagram 140

Subroutine XLOAD

Listing 143

Flow Diagram 143

-V-INDEX (Cont'd)

(8)

The SHIP STRUCTURE COMMITTEE is constituted _{to prosecute a research}

program to improve the hull structures of ships by an extension of knowledge

pertaining to design, materials and methods of fabrication.

RADM W. F. Rea, III, USCG, Chairman Chief, Office of Merchant Marine Safety

U.S. Coast Guard Headquarters

Capt. J. E. Rasmussen, USN

Head, Ship Systems Engineering and Design Department Naval Ship Engineering Center Naval Ship Systems Command Mr. K. Morland, Vice President American Bureau of Shipping

U.S. COAST GUARD

LCDR C. S. Loosmore, USCG - Secretary

CDR C. R. Thompson, USCG - Member

CDR J. W. Kime, USCG - Alternate

CDR J. L. Coburn, USCG - Alternate

MARITIME ADMINISTRATION

Mr. F. Dashnaw - Member Mr. A. Maillar - Member Mr. R. Falls - Alternate Mr. R. F. Coombs - Alternate MILITARY SEALIFT COMMAND

Mr. R. R. Askren - Member

LTJG E. T. Powers, USNR - Member

AMERICAN BUREAU OF SHIPPING Mr. S. G. Stiansen - Member

Mr. F. J. Crum - Member

SHIP STRUCTURE COMMITTEE

SHIP STRUCTURE SUBCOMMITTEE

The SHIP STRUCTURE SUBCOMMITTEE acts for the Ship Structure Committee on technical matters by providing technical coordination for the determination of

goals and objectives of the program, and by evaluating and interpreting the

re-sults in terms of ship structural design, construction and operation.

NAVAL SHIP ENGINEERING CENTER OFFICE OF NAVAL RESEARCH

Mr. E. S. Dillon Chief

Office of Ship Construction Maritime Administration

Capt. L. L. Jackson, USN

Maintenance and Repair Officer Military Sealift Command

Mr. P. M. Palermo - Chairman Mr. J. M. Crowley - Member

Mr. J. B. O'Brien - Contract Administrator Dr. W. G. Rauch - Alternate

Mr. G. Sorkin - Member Mr. H. S. Sayre - Alternate

Mr. I. Fioriti - Alternate

NAVAL SHIP RESEARCH & DEVELOPMENT CENTER

Mr. A. B. Stavovy - Alternate NATIONAL ACADEMY OF SCIENCES

-Ship Research Committee

Mr. R. W. Rumke, Liaison

Prof. R. A. Yagle, Liaison

SOCIETY OF NAVAL ARCHITECTS & MARINE ENGINEERS

Mr. T. M. Buermann, Liaison

BRITISH NAVY STAFF Dr. V. Flint, Laision

CDR P. H. H. Ablett, RCNC, Liaison

WELDING RESEARCH COUNCIL

Mr. K. H. Kooprnan, Liaison

(9)

-1-INTRODUCTION

This report presents the computer program

corresponding to

the method of design expounded in the Ship Structure Committee

Report SSC-215, 'A GUIDE FOR THE SYNTHESIS OF SHIP STRUCTURES

PART ONE - THE MIDSHIP HOLD OF A

TRANSVERSELY FRAMED,

DRY CARGO SHIP". The program consists in an executive routine,

called TRANSHIP, and twenty seven subroutines

arranged in alphabetical

order.

These have the following names:

ASPECT

2

BNMAT

3 COSTING 4

DOCUMENT

5

EV

6

FRAME

7

FSHAPE

8

GEOM

9

GRILLAGE

10

INPUT

11

INTERMED

12

LONGIT

13

LONGMAT

14

MATINV

15

NTPLATE

16

PLATING

17

REACT

18

RTPLSUB

19

SECTION

20

SHAPES

21

SOLVE

22

STEP

¿3

STRESS

24 T 25

THETA

26

TRANSV

27 XLOAD

For each routine, the following items are given:

Description

Notation

Listing

dl

Flow Chart

Appreciation is expressed for the help received by Mrs. Linda

(10)

-2-DESCRIPTION OF INPUT

Input data are introduced as follows:

ist Card

Format (9 I 5)

1K = i

Center Keel only

1K = Z

Z Side Keelsons plus Center Keel

1K = 3

4 Side Keelsons plus Center Keel

1K = 4

6 Side Keelsons plus Center Keel

IPILL =

O

No pillars

IPILL =

i

Pillars included

[HOG =

O

Sag case

IHOG = i

Hog case

NI

Number of iterations

2nd Card

Format (212, 8F8, 2)

NDECKS

Number of decks (double bottom excluded)

KPANELS

Number of bays

TLENGTH

Ship length(ft)

BEAM

Ship beam (ft)

DRAFT

Ship draft (ft)

HMAIN

Ship depth (ft)

XLHOLD

Length of hold (ft)

XLHATCH

Length of hatchway (ft)

HATCH

Width of hatchway (ft)

HFLOOR

Height of floor (ft)

3rd Card

Format (7F lO. 3)

TWEENH

'Tween deck height (ft)

4th Card

Format (7F 10. 3)

DKLO

Uniformly distributed deck ioading(lb 1n2)

5th Card

Format (3F 15.2)

YIELD

Elastic limit of material (lb

in2)

_-2

E

Young s modulus of material (lb in

)

Z

PRM

Criterion of primary stress intensity (lb in

6th Card

Format (7F 10.3)

HRATE

Labor rate (dollars

hrj')

DOLLPP

Plate cost (dollars lb

WPRICE

Welding cost (dollars lb

)

TWF

Time-weight factor (dollars lb

SG

Material density (lb in

3)

7th Card

Format (F 10.2)

(11)

(2

0 1 lO

(4

1

496.

71.5

30.

43.5

lOO.

33.

20.

5.

(9. 62

. . 63 . O.

17.76

15.54

2.1

1.78 35000.

30000000.

19000.

0.065

3.

0.02

0.283

0.

Program TRANSHIP

Tranship is the executive routine. It calls all the primary

sub-routines as needed and makes any computation

_{necessary for transition}

between them.

A general Flow diagram of all routines is given starting

on page

17.

A detailed flow chart of TRANSHIP begins

_{on page}

_19.

A listing of terms common to all the routines precedes the

_{flow diagrams}

for TRANSHIP, while those specific to a given subroutine are included

with the description of that subroutine.

_{In all cases, the appropriate}

mathematical symbol from the basic report SSC-215

_{are given when}

applicable.

-3-Sample input values used in the computer runs for the WOLVERINE

STATE are as follows:

(12)

-4-Notati on

Mathematical Symbols to FORTRAN

Mathematical

Definition

FORTRAN

Symbol

Term

a

Spacing of frames and floors, or side

of plate rectangles (input)

A

Frame spacing

AAA

b

Distance between longitudinais

B

Effective width of plating

EFFW

b

Effective breadth or width of plating

EFFB1, EFFWI,

e

_EFFBR

d

Height of floors (input)

HFLØØR

f(A.)

Time - size factor

F4(H)

f1 (h)

Weight per inch of vee weld

Fi (H)

Weight per inch of continuous double

fillet weld

F2(H)

f3(h)

Weight per inch of intermittent

double fillet weld

F3(H)

f(w.)

Time - weight factor

TWr

h

Thickness of plating

THIKK, THIKK1

Wave height

WAVEH

r

Hourly rate of labor (input)

HRATE

A.

Grillage matrix

AX

Ape

Effective area of plating

APE

B

Beam of ship (input)

BEAM

B (x)

Boundary matrix used in grillage

p

_calculation

B

C

Combined cost of plating, labor

and welding per inch of ship length

CcZSTS

Cf

Total labor cost of erection

e

_{and fabrication}

_XLABOR

C

Total material cost

PLCST

(13)

-5-Mathematical

Definition

FORTRAN

Symbol

Term

C

Total cost of welding

TWELD

w

D

Depth of ship (input)

HMALN

Flexural rigidity of plating

FLEXR

E

Youngs modulus (input)

E

H

Draft of ship (input)

DRAFT

I

Second central moment of area

SMAIF

I

Second central moment of effective

pe

plating (about the faying surface)

XIPE

Second central moment of plate-frame

p

_combination

_XIPFX

Second central moment of

transference

XIT

K

Stiffness factor

XJ

L

Ship length (input)

TLENGTH

N

Number of weld seams in inner

bottom

NSEAM

S

Total weight of hull structure

per inch

WBAY TOT

Influence coefficients

ALFA

1,3

Aspect ratio parameter

BE, B2

y

Specific weight of material (input)

SG

e

Joint rotation

T

X

Normalized effective breadth

El, E2

IT

Circle ratio

PIE

Criterion for primary stress intensity

(input)

PRM

Primary stress intensity (design stress

-

intensity)

PSL

o

Secondary stress intensity in longitudinal

2x

_direction

SESL

o

Secondary stress intensity in transverse

(14)

cb' aXb

px

cpy

X

Maximum bending stress intensities

in the longitudinal direction (at the middle

of the sides) when the compressive

stresses are of such magnitude as to

cause buckling

PHIL

PHIS

SIXE

-6-Mathematical

De finition

FORTRAN

Symbol

Trm

a

_{Secondary stress intensity in :he}

Zx

longitudinal direction in the flange

FSESLI

a2

_{Secondary stress intensity in the}

y

_{transverse direction in the flange}

FSEST1

a3

_{Tertiary stress intensity in longitudinal}

X

direction

_TESL

a3

_{Tertiary stress intensity in transverse}

y

direction

_TEST

Gcr

Critical buckling stress intensity

_CRITB

Xb

_{Maximum bending stress intensities in}

the longitudinal direction when there is

no compressive stress

SXB P

a

,

a

As SIXB, but in transverse direction

_SIYB

yb

a

b

Maximum bending stress intensities in

the transverse direction when there is

no

compressive stress

SYBP

Factor depending on aspect ratio in

equation for SXBP

Factor depending on aspect ratio in

equation for SYBP

Distance of neutral axis of plate-frame

(15)

FORTRAN

Term

-7-Notation

-FORTRAN Symbols to Mathematical

Definition

Mathematical

Symbol

A

Spacing of frames and floors, or side

of plate rectangles (input)

a

AA

Storage location for A

AAA

Frame spacing

a

AFAC

Scale factor = 1.0

ALFA

Influence coefficients

ALAD

Loading (life and/or hydrostatic load)

AREAT1

Cross -sectional area of transverse material

AX

Grillage matrix

A

i, i

B

_{Distance between longitudinals}

b

B

Boundary matrix used in grillage calculation Bk)

BB

_{Storage location for B}

BEAM

Beam of ship (input)

B

CC

Coefficients of characteristic equation in

descending order

CN

Determinant of grillage rrtrix

CØDE

Code =

i fixity of plating

Code =

-

I simple support

DEFL

Deflection

DKLØ

_{Uniformly distributed deck loading (input)}

DQLLPP

Dollar per pound of plating material (input)

DRAFT

Draft of ship (input)

H

E

Young's modulus (input)

_E

EFFB1

_{Effective breadth of plating}

_b

e

EFFW1

_{Effective width of plating}

_b

(16)

FORTRAN

Term

FIX CDE

FSECLØ

FSECMØD

F SE C TR

FSESL

FSESL1

FSEST

FSEST1

-8-Definition

Flag denoting ideal fixities of the

inner-bottom structure

= -

1 fixed support

1 simple support

Frame section modulus in longitudinal

direcLion

Required section modulus

Mathematical

Symbol

GP

Grillage pressure load

GNT

Number of oil tight double bottom

longitudinale

GNNT

_{Number of non-tight double bottom}

longitudinal s

HATCH

Width of hatchway (input)

HF1(I)

_{Plate thickness if the plating is fixed along}

the edges of the plate rectangle (I

= 1

indicates bottom shell,

I =

2 indicates

inner bottom)

HFLØR

Height of floors (input)

d

HMAIN

Depth of ship (input)

D

HNEUT

Height of neutral axis of the midship

section

HNTX

Height of neutral axis of beams

HRATE

Hourly rate of labor (input)

r

Frame section modulus in transverse

direction

Maximum allowable secondary stress

intensity in longitudinal direction in flange

Secondary stress intensity in the longitudinal

direction in the flange

f

Maximum allowable secondary stress

intensity in transverse direction in flange

Secondary stress intensity in the transverse

direction in the flange

(17)

-9-FORTRAN

Definition

Term

HS1 (I)

Plate thickness if the plating is 3imply

supported along the edges of the plate

rectangle (1 =

i indicates bottom shell,

I

2 indicates inner bottom

HWF(L)

Thickness of oil tight floors (L = 1)

or oil tight double bottom longitudinals

(L =

2), fixed along edges

HWFS(L)

_{Thickness of oil tight floors (L =}

₁₎

or oil tight double bottom longitudinals

(L = 2), simply supported along

one long

edge fixed along the other

HWS(L)

_{Thickness of oil tight floors (L}

₁₎

or oil tight double bottom longitudinais

(L =

2), simply supported along edges

IßT

Limiting index on innerbottom

IHcG

Code = O Sag case

Code =

i Hog case

1K

Code =

i Center keel only

(input)

Code = 2 2 Side keelsons plus center keel

Code = 3 4 Side keelsons plus center keel

Code = 4 6 Side keelsons plus center keel

ILG

_{Limiting index on side shells}

IPILL

Code =0 No pillars

Code =

i Pillars included

IPLATE

Number of plating items

IPLT

Code =

i

Bilge

Code = 2 0. T. Floors

Code = 3

N. T. Floors

Code = 4 O. T. Longitudinals

Code

5 N.T.Longitudinals

ITR

_{Limiting index on decks}

J

Code =

i

outboard

Code = 2

inboard

(18)

-lo-FORTRAN

Definition

Ma thematj

Term

Symbol

LLGX

Code =

i

denotes longitudinals

Code =

O

denotes transverses

LONGLAB

Label For Longitudinale

i

Keel

Z = Side girders

3 = Fourth deck

4 = Third deck

5 = Second deck

6 = Main deck

NBBG

Code =

i

corresponds to

i

double bottom longitudina

Z corresponds to

3

double bottom longitudina

3

corresponds to

5

double bottom longitudina

4 corresponds to 7 double bottom longitudina

NDECKS

Number of decks (input)

NG

Number of longitudinal girders

NNT

Number of non-tight longitudinals

NQC

Code =

i

Unrestrained frame deflections

at girder intersections

Z

Influence coefficients

3 Frame bending moment distribution

NØT

Number of oil tight longitudinals

NTS

Code =

i Bulkhead Bay

= Z Hatch Bay

NSB

Number of unsupported frames in

bulkhead bay

NSH

Number of frames in hatch bay

NSTT

Number of non-tight floors

PIE

Circle ratio

PLATLAB

Label For Plate Elements

= Double bottom

Z = Inner bottom

3 = Fourth deck

4 = Third deck

5 = Second deck

6 = Main deck

7 = First side sirake

S = Second side strake

(19)

10 = Top side strake

11 = Bilge

12

_{Oil-tight floor}

13 = Non-tight floor

14 = Oil-tight longitudinal

15 = Non-tight longitudinal

PPP

Concentrated loading

PRHEAD

_{Pressure head}

*

PRM

Criterion for primary stress intensity (input)

o,

PS L

_{Primary stress intensity (design stress}

Q,

intensity)

RABI

Bilge radius

RL

_{Output argument from EV for adjoint of matrix}

RT

Array of real roots of characteristic equation

SESL

Sec9ndary stress intensity in longitudinal

direction

X

SEST

Secondary stress intensity in transverse

direction

y

SG

Specific weight of material (input)

SHMAX

Critical buckling shear stress

SMAF

Second central moment of area

I

SMALl

Second central moment of area of

longitudinal girders

SMAT1

Second central moment of area of

transverse framing

SMDL1

Section modulus to the plate of longitudinal

girders

SMØDL2

Section modulus to the flange of

longitudinal girders

SMØDT1

Section modulus to the plate of transverse

framing

SMDT2

Section modulus to the flange of transverse

framing

-11-Mathemaica1

FORTRAN

Definition

_Symbol

Term

(20)

-42-FOR TRAN

Definition

Term

SS

Spacing of stiffeners

SUMAREA

Total cross-sectional area

of plating

SUMAREAL

Total cross-sectional area of

all webs of

longitudinals

SUMAREAP

Total cross-sectional area of

plating of a

ship section

SIJMML

Total first moment of area of longitudinals

about neutral axis of ship section (neutral axis

from previous iteration or assumed)

SUMMM

Total first central moment of area of ship section

SUMMP

SUMSMA SU MSMA 1

SUMSMAL

SUMSMAP

TESL

TES T

Total first moment of area of plating material

about neutral axis of cross section

Total second central moment of area of ship

section

Correct second central moment of

area of ship section

Total second central moment of area of

longitudinals of ship section

Total second central moment of area of

plating of ship section

Tertiary stress intensity in longitudinal

direction

Tertiary stress intensity in transverse

direction

Ma thematic al

Symbol

3x

3y

THIKK1

Thickness of plating

h

TLENGTH

Ship length (input)

L

TRANLAB

Label For Transverse Elements

1

Oil-tight floor

2

Non-tight floor

3 z Fourth deck beam

4 z Third deck beam

5 z Second deck beam

6

Main deck beam

7

= Innerbottom to fourth deck frame

8 = Fourth deck to third deck frame

9 = Third deck to second deck frame

10 = Second deck o main deck frame

(21)

-13-FORTRAN

Definition

Mathematical

Term

Symbol

TWEENH

Tween deck height (input)

TWF

Time-weight factor (input)

f( w.)

V

Array of shear coefficients

VS

Array of slope coefficients

VV

Array of moment coefficients

VY

Array of deflection coefficients

WAVEH

Wave height

h

WEBL1

Web heights of longitudinals

WPRICE

Cost of weld material per pound (input)

XFL

Length of stiffener (between joints)

XII

Maximum field bending on longitudinal

elements

XKHcZG

Hog coefficient in equation for wave

bending moment

XKSAG

Sag coeffícient in equation for wave

bending moment

XLHATCH

Length of hatchway (input)

XLHØLD

Length of hold between two bulkheads (input)

XLL

Length of longitudinal girder

XLPANEL

Length of bay

XMANH

Number of man hours per square foot of

the 'equivalent surface' of the structure

XMMT

Maximum field bending moment on transverse

elements

XNG

Number of double bottom longitudinals

XNST

Number of non-tight floors or frames

per bay

(22)

-14-FORTRAN

Definition

Term

ZL

Vertical coordinate of centroid of

longitudinals

(23)

PPOG8l TkANSHIP CÙMMON/TLEN, TLUTH

,IMUt

,RAU),XKSAj,XKHUtj

DIMENSION HV1(2),S1(2,Hw (2),MWS(2),cR(2),cs(2),hs(n,

_{1CR1 (2), 411 (2) ,14C2. (2), w I I () ,H12 (2) * (2 (2), Mw) ¿ (2)} COMMON/PK/ DKLU(3,fl)

-COMMON/LANEL, TA'JLA(1U),ILA1LA)3(j),LON(LA(Ö) ÇfMMON/SH/

SHMAX

MMflJiflII fl1(*fl%

lJ*j

EF8R( 12* 2) , AREATj( 12, 2, 51,WEbl( 12. 2.5) ,SMAL1(1/,2,) ,SMO)JL1(17,2 4.3).WFbLi(7*2.SS),A14EAL1(7,*3).SUMAHAp(S),SUMAREAL),SUl1Mp(5;,SU 6(7),TIKKj(j7,.Ç3UE(1/),ALOAD(j7)*suM5Ma*(5),suMSMAL(,),suM5MA(,)

OMMON/GR/ Bl(9.9)a32(9,9),M.S(9,9)iXL1(ç).xL(9),Zj),ULt17),

i

ITEST.IPILL,NOUS*N,NSR,NSS,D.WI.MM*NUTS,NM,Nt,NU,NJ,XFL(1,), 2 SMAV(15) DATA(LOAHa4kKEEL.

S '.IRDER4TII UECX3kU DECIÇ2NI) ¡JECKMAIN DK

)

OATA( TRANLA:8UH0T Fj..UONN1 FLOOR4TN UECK3HI) UECK2ND

ClÇIiAI

UK

I

1,R.-4TN4T14.3H1) 3MO-2MO 2ND-MAIN)

-OAIA( PLATLA):12UNIL SPIELLIN BUTT,4rM DEcK3RD )Et,k2NU LJICKMAIN OK _{ISTRAKE 4SIRAKE 3STNAKE 2S StRAK}

STNAIÇEUT r(.OURNT

LUOI4OT LUND,NT

2 LONG.)

9991 FORMAT(*

TEST ¿UNE a, F1Q,1/)

N OJ T: II RABIa 120 0 REPP 921b, II, IPILL, IHIE, (II

I 9218 ORMAT(9I5 I -CALL INPUT(3E)UIN0,NG, 1K .NJ,XFL.XII ,NSH,NSb,N( MPM: N B IO G IF(NIOPG.EO.1) NRMa7 VIXCOUFa-i.fl COMMON/COST/ HNATE,OULLPI.',WPNIUE,TwF SD -NX2:j CflMMDN/LO/FSECLOIi2),rESLI17),FSEsLi(i7),SuoLL2(t/,2,),i.LL(l7, C00F7

lin

-

--C(3MMON/TR/ 3411 CONTINUI: 1 17.2,5) iQ-JEIIT a O

COMMON IC, TWELD

347 CALL SECTIUN( (MEUT, SMUNEC,KUUNT,BENUING,NX2) COMMON/rI HF1,HS1.H,$Mb,14WIS,H12,Hw12,IfL,IIL CALL DOCUMENT COMMON/NI XMMT(17) PRINT 701 DIMENSION NSA(,),PPI'(V),wL(9),Ax8(i5,i5) 7111 rITRMAI(IHO,* T4ANS, SMAT1 SìIOUT1 MOUT ARATt

COMMON lEI IBT,E

IFFIOl ZP ALUAU*,//) COMMON /0/ VS(iS),VY(5) )O 702 1:1,10 DIMEtSION SEST1(17),StSL1(17)

7112 PRINT 705 ,TNANLAB( I).

SMAT1(I,1,1),SMOUT1(I*1.1),SM0012CO,1P1), COMMON /8/ A,R,NUEG'(S.IJLAIE,PSL,PKM,I.MAIN,NNEUT.TESL.ILST,SESL, i AREAT1(1.1.1),EFFIOI(I,i,1 .2

I!_!

AO(I) 1 SESTp KPALS,IOEAM,XNG,HATCN,TH1Ì;ai.F-BN;LFF-Bi. PRINT 703 2 EFVW*EFFW1.I*LU)N* IEENH,THIKX,bBk.WEM.WEB1,XII.*X*Cfl),

71i3 FORM&T(1IlQ.* LONU.

SMALl SMOUL1

SMOUL2AI4AL1

E 3 XIS, AREAT1,S'IAL1,UDL1*WEBL1,SUMAREAAWEAL,SUMMP,SUMML, lFFwl

iL *,//)

4 AR,SUMAREA,SM3M,WPL*XPANEL,WBAYLP,XNST.NST,IIN, DII 734

Il.ô

5 WTR.BAYTI4,SJMWT.SUMWTR,LLGX, YIELU,WIOAYTN,XLAdOH,PLCUSTP

-7P4 PRINT 735 ,LUNGLAB(I ),SMAL1II,i,1)

,SMUUL1II,j,1),ST400L2TI.1,1), 6 COST,COSTMIN,GNT.IiNNT,KNEUT,SUMSMAP,SUMSHAL.SUMMA.COUE,

i ARF4L1(I,1il),E1h1CIhl,tI*ZL(I)

lWAVEH,PRHEAD,DRAT,ALOAO,XLHULU,XLpANE1j,A,AA,LjELrES1, 735 cOMMAT(* a, A8,7F1O,21 R SUMSMA1,LPLLL, IL.,AMEAL13ASX, SMATX,5MQDI1 NA20 CDMMON/C,/ FIxCODE PRM1:PRM COMMON

ALFAÇ1,i5),AXI 9,9 ),CC(1u),LN,UEI-L( 9),I.-',Ii,.j,J,K,NIO,N5,

II(HNEUT,LE, fl, * 4l'111.N) (iO TU 11 1PO(15,4fl), PRM1:PRM (PIMAINHNEJI) /lINEUT LXMIj15)*XM(15>,XA(4U).1-R*SUAL,AFAÇ,NNN

-11009547 ¡:1,10

DIMENSION lP(17),SESL(i7).SEST(j7).TbsL(j/),TT(j7),oN,),THjKK(

PSLL(1):PRM1*(HU1_2Il)/(NMAI)__

P.SLL( I ):APSF(PSLL( i) I

IF( ¿P(I) ,LI,

INEUT

I

LL(I)P5LL(I)*XKNUG/SAG

9547 CONTINUE

C

--GRILLAGE JOiNT CALCULATIN

3000 CONTINUE 351 CONI ¡NuE -NSS3 9871 CONTINUE CALL GRILLAGE i ALOAD,SMATOISMAL1 1 DO 3050 LLL:1,NG IF(ROOT(LLLH 2400,-S0U,3U50

2400 WRITE OUTPUT TAPE b,9144,J,HUOT(J)

GO TO

9999

3050

CONTINUE

30O PRINT 9139, cRuT(J),J:i,Nur

(NOECKS, HATCH.PRl4AD,BAM,L1)UR,AAA,

SECTION MODULUS C11LCULAT)DN

pP MC a P141-4 P014:19000, SMONEC a BEND INIO / PM NOUS: O K DONT : Q C X a XLL 0110 CALL 1-INMAT(AX,

I,

1)

(24)

CALL (4NMAI(AX, 2, ff2T -CALL(4NHAT(AX. 3, dl) N1:NM - --1120 CALL

LflA[)(AX,2XL2,XLL)

-I(Nd(G.GT1

U ID 9/lb

CALL LÚAD (Ax, 1. 4L1,XLL) PPP(6):PPP(5) NSSNSs_.j. PPP(5)tPPP(4)

JF(WSS,LT.t) nÑ 36u6

3I9ÇÙT.ÑF

PPP(3):PPP(2)

1410 CALL SOLVE (bi,

,2.thG XL1, XL2, NUUS,NSS.,II1ILL)

I(NSS-1)36fl1,91,U,j6u1 NItNG*1

__

--. S950 CONTINUE

9736 CALL FRAME(SMAV,NJ, XEL,OL,W1 ,O,Z,HthPP,U,,N1,NUTS,NOC,WL,KK,NM

DO 1?t(1, I:iNU

36(15 CONTINUE -3O'1 CONTINUE

vsi-vscii

IV(NSS,LE,fl) U TD 3615 V (I>=V (I)/E/XILJ(1) DO 831J. I:1,NG --17(11 CONTINUE

V (1):V (I)**XID(1)

GO TO 146fl VV(I)VV(I)**X1UC1) 36(16 DO 36fl7 itj,rj 83(11 CONTINUE

V (lIt-V Ì)//icÎu(í)

r, ro 3601.

S(I)z-VS(j)

3611' COTINUF

VV(l): VV(I)/E/XID(1)

-XJA:SS.(NSHG1) 36(17 CONTINUE Ir(NSS,EO.0 AND, .UT,b) CALL STP(AX,N5,NS2,xJA,M) 36(11 Ir(NSS.Gr,o) GO TU 94/1 C _14A0 ÇONTjNUE C STRESS CALCULATION C CALL STRESSI N,N) C GPjLL4GMMbEH GALCUÌATIUN Ir(NN) B4fl2,b402,33 NS3:4 --8402 IF(N) 84fl0,b4flU,3S LF(NSS,FD,0,ANU,NG.JTIØ) NSJt 123 KKM t (KM +

i

,1S2:NS,7 <((MM t4 I(K((M - KKMM)324,J24.6400 NSA(2):2 -4PRINT 326, K((M NSA '4S2-2 SESL(2) t SEL(1) lr(NSS,Ew,n) NSA(3).'12/2.T SFST(2) ¡ SEST(-1) 'ASA(4 ):NS2 NUITtNOIT.1 00 36115 Jti,NSJ Ir(NOIT.GE,N1 uo Tu d400 KK:NSA(J) GO TO 34? 8ÚCONT1SE

--XJJ:FLOATE(JU)eSS PRINT 326, KM PRINT 76,XJJ 1< t 76 FORNAT(1i.41,* STEI LOCATION

t e

L = 1 lF(NS.EQ,fl) P'IÑF /8

'

-' -78 rO(MAT(/,. 181CM 38V

e,j

PRINT 2221 ,frIXCU0 Ic(NSS,EU,1) P1NT 727 FORMAÎ( e FIXCOL1 s 77 FÙNMT(/, EULKM0U VAY

*./)

2222 FORMATI * e,A8.F4,1/) CALL STEP(AX,K((,NS.XJJ,N(1M) DO 9632 11,15 lr(Kç-NS2( Jbfl4,360,,3bU1' 9632 !RINT ¿222,

I),(I)

361(4 CONTINUE --(F(F1<C0E.C9. I) G1 T1 801 00 9U27 I:i,NG lF1(t) THlK1(1,1..K) 009127 Kti,NG

4r1(2) t TH!lcK1(.1,(()

AXI(I,K):ALFA(I,Ic) IrL TIR + 2 902? CONTINUE -IlL = !TI- + 4

CALL REACT tNG,AXd,P0I,W,UFL)

HWI(1) : TI.4lK1(1rL,1,r) S'-4114X50, HW(2) = T1-1IKI<l(OIL,1,K). 00 8319 I:1,NN8G FINCTIDE = 1.0 8319 IF(SHMAX,LT, ASCP-'(1?)) SHMAX:Ads(ppp(I)) NOC 3

(25)

00 8502 1 1 IPLAIL R52 CODOCI) --GO TO 344 8501 HS1(1)

THIKKj(i,j,Ç)

RS1(2) THIçK1(2,l,K) NW5(1) TH!r.K1cjrL.i.,Iç) = THI;K1(I11.,1,K)

RJS(L) = 0.SC-)IF(L) I- HIS(L)) _{CALL I1TZLE0}

9999 CALL DOCUMENT

CA L_ÇQSJ_

FORMAT (2)15,2) 326 FORMAT( * KM*,i3/) i9 FORMAT (1Hfl,j0X13r uu, // 51, ii REAL 9144 FORMAT (1Hfl,lflX, 6h R3UT(I,4H)

:1J5)

END Macr:

'ow Diagram - No Loops Sho

GEM

PLATING

ASPECT

NT PLATE

(26)

FRAME - MATINV EV RTPLSUB

[EV

LT - THETA

EV

- THETA

rEv

[T - THETA

EV _{T - THETA} EV _{T - THETA}

IN TERMED

FRAME

MATINV

STRESS

(27)

Prorarn TRANSHIP

(detailed flow chart) NØIT = O

RABI = l0.O

R EAD JE, ITEST, IPILL, IHlG, NI

CALL INPUT (BENDING, NBBG, 1K, NJ, XFL, XII, NSII, NSB, NG)

NBM = NBBG NBM = Z Section Modul Calculation PRMC = PRM PRM = 19000 SMDNEC BENDING/PRM

NIØS = o KØIJNT r 0 FIXCDE =

-1 NXZ r i CØDE(7) = KNEUT = O

CALL SEC TIIN(KNEUT, SMDNEC, KØUNT, BENDING, NXZ)

CALL DØCUMENT

Print Values Calculated in SECTIQN for Transverses

(28)

NXZ = O PRM1 = PRM Yes PSLL(I) = PRM1*(HNEUT -ZP(I))/ (HMAIN - HNEUT) PSLL(I) = ABSF(PSLL(I)) HNEUT g O. 5'HMAIN Dl I = 1,10

--No

PRM1 = PRM*(HMAIN - HNEIJ T) /HNEIJ T

PSLL(I)

PSLL(I)*

XKH'G/X1ÇSAG

Grillage Joi Calculation

CALL GRILLAGE (NDECKS, HATCH, PRHEAD, BEAM, HFLR, AAA, ALØAD, SMAT1, SMALl)

Grillage Matrices

[ X

ILL

_t

CALL BNMAT (AX, 1, Bl)

t

CALL BNMAT (AX,2, BZ)

t

CALL BNMAT (Ax, 3, B3)

f

CALL XLAD (AX, Z, xLZ, XLL)

t

CALL XLAD(AX.1,XL1,XLL)

(29)

I

L

n I = l,NG

t

V(I) = -V(l)/E/XID(J) VS(I) = -VS(I) vv(I) = VV(1)/E/XID(I) Grillage Member Calco ation

CALL SILVE (BI, BZ, BZ, B3, XL!,

xLZ, NIS, NSS, NBM, IPILL)

--D I = 1,NG VVU) = - VV(I)/E/XID(l) vS(l) = - VS(l) V(I) = V(l)/ E/XID(I) NSA(4) = NSZ DØ J = !,NS3

t

KK = NSA(J)

JJ r KK

-i XJJ = FLATF(JJ)*SS Print XJJ

CALL STEP (AX, IcK, NSZ, XJJ, NBM)

t

Yes

_{\ Print 'BULKHEAD BAY/}

NSZ = NS + Z No NSA(I) S i NSA(Z) = Z NSA(3) = NSZ -Z No

(30)

CALL REACT (NG,AXB,PPP,W,DEFL)

r

130 I = 1,NBBG NOC = 3 Nl = NM Yes

L)

Yes SHMAX = pPp(I) Nl = NG +

CALL FRAME (SMAF,NJ,XFL,DL, WI, D, Z, HB, PPP,Q, E, Nl, NOTS, NOC, WZ, KK, NBM)

CON TI N L ¡ E No DO I = 1,NG

t

V(I) V(I)*E*XID(I) VV(I) = VV(I)*E*XID(I) PPP(6) = PPP(5)

PPP() = PPP(4) PPP(4) = PPp3)

PPP(3) PPF(2) PPP(Z) O XJA = SSs(NSH -f 1)

(31)

Stress Calculation _{CALL STRESS (N, NN)} No Yes KKM = KKM + 1 KKMM 4 No Yes o

SESL(Z) = SESL(1) SEST(2) = SEST(1)

Intermediate Fixity Calculation

\ Print

DO I

=

1,15

PLATLAB(I), CODE(l/

HF1(1) = THIKK1(1,1,K) HF1(Z) = THIKK1(2, 1, K) {FL = ITR + 2 IlL = ITR + 4 HWF(l) = THIKK1(IFL,1,K) HWF(2) = THIKK1(IIL,l,K)

FIXCODE = 1.0 I DO I = l,ÍPLATE CODE(I) = -1.0

L)

HSI(1) = TFIIKKÌ(1,l,K) HS1 (2) = Tl-IIKKI(2, 1, K)

HWS(1) = THIKK1(IFL,1,K) HWS(2) = THIKK1 (IlL, 1, K) HWFS(L) =

0. 5*(HWF(L) + HWS(L))

t

CALL INTERMED

t

CALL DOCUMENT CALL COSTING

(32)

i.

Subroutine ASPECT

Abstract:

Subroutine ASPECT is called from subroutine SECTION.

It calculates

the aspect ratio of all plating and calls subroutine PLATING or NTPLATE

as appropriate.

Terms specific to this subroutine:

-24.-Definition

Side of plate rectangle used to

determine aspect ratio

Length of plate rectangle used to

determine aspect ratio

Code =

i

denotes non-tight floors

= O denotes non-tight longitudinals

FORTRAN

Term

Al

Bi

(33)

SHdWcJUJIN 01 COMMON/ILEN,

TL%TH

1ti'tU ,R,XÇSAG,X*HD

cOMMON,nEr

ori

-COMMON /E/ IbI,E

GO TO 919

i

SEST,

KPAL,dE4M.XNG, HATCH, THItçç, THIK(1,

R,EFFB1, 59 IP.T:l-phNDELXS.91 2 TK<(IÌHI)U((i) 3 XIS, AREAT1,bMAL1,UD1,WELi,SUMAAP,SUMAREAL,SUMMP,SUMML, IF(J,FU,2) UO TU 9111 ---GO TO (67,68,69,5.B)LPLI 5 ,TRAyTH,SUMJT,5JM,1R,LLGX, YltLjJ,WHAYTH,XLAUN,PLCUST, C dILGE h 67 CONT IUE 7HAVEH,PRHEAD,DRA T, ALOAO,XLMULU,XLPANL.AAA,Bdb,AA,dd,ULTSL, C CALL RILGE 8 SUMSHAj,ZP,LL,ILG,AiE,AS), SMAT[IffTT EFERR(1) r AA DIMENSION ¿P(17),StSL7),SSTI17),1L(17),TESTC17),XNG(/),1HLKK( EFEWIl) = 117),

TH1rr1T,2yFTt1ErrTÏnfTtT2,ä).nÑ1(i7

SO TO 119 2,2.5),AREA(17,5),T(),SMATi(1,2,MODi(t2,2,5),THAKXCi2,2 C O,T,FLOORS

[12

EArii122

111 25TS

CTf.2 SM(L .2 Ad A1rHELOON

4,3), wI-EIL1(7,.3), AMEALXC 7,,3) ,SUMANCAP(5),SuMARAL(5) ,SUMMP(5),SU

PSI. = 6(7),TIKK1(i7,S),CODEU/),ALUAD(i7),SUMSMAP(5),SUMSMAL(5),UMSMA(5) TESL(I) YIELL) -TFST(I) OftLil - sESTU) A1:AA SESTIl> SET (1)

IV(HHEUT,LE,Th.5 HMA1uÏ

-SESLII) fl,fl PRM PRH * (HMAIN-MNJJ) /HNUI GO 10 919 11 Dû 111 I=1,IPLAT

PSL =PWM *(HNEUT - 1PII))/IHMAIN -HNEUT) PSLABSFCPSL)

-.

Ir(

ZP(I> ,Lî.HNEUT I SSL( I)AHSF(SESL(

¡Y

SESTIl )AE4SF (SbSTC II) TESLI1) r VI¼LO - PSI.

SEL(!) TESTI!> YftLU SSTI)) IF(TEST(I).L.0b,>iESTU) 3000,

--IF(TESL(l),LT. lEU,) tESL(I)= 3000,

DO 111

'1

DO 9111 J:l, Air AA 3j HU lEI I-2NUECKS.l0T)90,9U,S9

90 jF( I.NDECKS-TISE, 55.1U 58 1F I.IBT)51,S1,S0 51 CONTINUE GO TO 919 C OECKS__ 50

CONTINUE GO TO (55,56J ElI: (NEAM.HY GO 10 ji)

56 BjIIATCH 10 CONTINUE GO TO 919 C WEBS C SIL)E SHELL i12 Ifl:I-NUECIS-IbT I1(J,EIJ,2) c,ff TU 9jj 01= TWENC!LJ) N T, F LO O R S 69 R1:N0 NIE LOON L000E1 CALL GO TO 119 C O,T,LONGITUGINALS 75 IO1HFLJOR A1;8A TESL( Il r YIELU-SES.111

TESI(I) = YIL[) SEST(I) SFSL(I)

r SESL(1.I uO TO 919 C N,T,LONOITUDINALS C TLONG 11M LCODEO CALL NTPLAT(I , J,A4,

LOO,PSL, INlIcK, YIELU,E , LCUÊTh, ONG, BEAM, CODE)

GO TO 119

919 CALL PLATING

(At,i,>,J)

119 TWI(pÇ1(l,J.K)

:TMÍ,Çjj)

EFFPi(I,.J,p() EFFW1(I,J.K> r EF(I)

-Riti CONTINUE iii CONTINUE RETURN EJD

(34)

.

Subroutine ASPECT Enter PRM = PRM (HMAIN - HNEUT) IINEUT No 1)0 I = 1,IPLATE

Calculate primary, secondary and tertiary stress intensities frr each member

TEST(I) = 3000 Side Shell No No TESL(l) < 100 Yes TESL(I) = 3000

Bitge and Web Platings of Floors and Longitudinals

IDO K = I,KPANELS DO J 1,2

t

SESL(l) = ISESL(I) SEST(11 = ISEST(I) I Al = AA

TESL(I) = Yield - PSL - SESL(I)

BI = NB

(35)

IPLT =

I

Bilge Plating

O. T. Floors

Thickness as for bottom shell

t

EFFBR(I) = AA EFFW(I) = 80

PSL = 0.0 TESL(I) = YIELD TEST(I) = YIELD

- SEST)I)

SEST(I) = SEST(I) SESL(I) = 0.0

t

CALL PLATING (Al, Bi, I, J)

T. Floors

t

o.

BI = BB LCODE = _CALL

N TPLA TE (I, J, AA, HFLOOR PSL, T}U XE, YIELD, E, LCODE, XNG, BEAM, CODE)

1<

4 _{N. T. Longitudinais}

T. Longitudinais

BI = HFLØOR Al = AA

CALL PLATING (Al, Bi, I, J)

o

LCODE = O CALL NTPLA TE (I, J, AA, HFLQJOPI, PSL, TI-lIKE, YIELD, E, LCØDF, XNG, BEAM, CODE)

Yes IPLT = I -Z NDECKS IßT THIKK(I) = THIKK(l)

f2

N.

t

TESL = YIELD - SESL(i)

(36)

E ID = I - NDECKS _{- IBT} Bi = TWEENH(ID)

f

CALL PLATING (Ai,Bi, I, J)

TF1IKK1 (I, J, K)

THIKK(l)

EFFB1(I,J,K) = EFFBR(I) EFFWI (r, J, K)

EFFW(T)

CALL PLATING

4

Deck plating

(Al, Bi, I,J)

Bottom shell and inner bottom

Deck Panel Outboard

Deck Panel Inboard

of Hatch Side Girder

Bi = (BEAM

- HATCH)

* 0. 5

CALL PLATING(Al,Bi,I,Jl

t

TI-IIKK1 (I, J, K) = THIKK)I) EFFB 1(1, J, K) = EFFBR (J) EFFWI)I,J,K)

EFFW(I)

CcNTINUE CØNTINUE _(RETURN

Bi = HATCH

j=i

(37)

.

Subroutine BNMAT (A,B,M)

Abstract:

Subroutine BNMAT calculates the boundary matrices used in the

grillage calculation called from subroutine GRILLAGE, SOLVE

and STEP.

Terms specific to this subroutine:

FORTRAN

Definition

Term

A

Grillage matrix

AL

Characteristic matrix

M

_{Index of Nielsen function}

(38)

-29-SUBROUTINE BNMAT(A,M,0) COMMON lEI IUT,E COMMON ALFA, X, CC. CN, UEFL, P, N. J..), (, NG, NS, PD, 1ROOT, S, SS, V, VV, ), X, XID, Xjj, XL, XLL, XML, XMl), XX,NR,SC 2 ,AFAC,NNN DIMENSION

ALFA(15,15),AX( 9,9 ),Ce(jQ),DEL( 9>,PU(15,4fl),

1RQOT( 9) ,S(1),V(1),VV(j5),W(j5) ,XID(15).XML(15),XMu(15),XX(40), 2X1140),A(9,9),X(9,9),AL(9,9),NL(9,9,l)(9,9)

98 FORMAl' (13) 91 FORMAT (Fi5,) Subroutine BNMAT(A, M, D)

n

Enter f DO

J = ING

G)I,J) = 0.0

r

Go, J) = 6(0, J) TT''RL(I, J)

(RETURN)

DO 1.000 I1,NG DO 1.000 Jj,NG

LL

G) I,.j)fl,0 1.000 CONTINUE CONTINUE DO 5000 Ii,NG oo 4000 I1,NG DO 3000 .J1,NG DO JR = 1,NG f AL) I, J) $A( I, J) /ROUT( IM) DO I = 3000

CONTINUE ALCI , I) AL (I, I) '1,0

DO I = 1,NG 4000 CONTINUE J = 1,NG CALL EV(AL,NG,Rt,,CCO

TTT( X, IR,M, NG, ROOT, CN, AFAC)

I,NG DO J DO 4600 I1,NG DO 4500 ..i1,NG G) I ,..J)'G( I ,J).TTRL( L ,,j) AL(I,J) -A)), J)/ROØT(IR)

B(I,J) = 0.0

45fl0 CONTINUE 4600 CONTINUE 5000 CONTINUE

f DO 6000 jj,NG

L

ALU, I) = ALU, J) -- 1.0 DO IR = 1,NG DO 6000 .Jj,NG

B( I,J)0,fl

DO 6000 IR:t.NG

t

CALL EV)AL, NG, RL, CC) B(I, J) = B)1, J) -B (I,J) -B (I J) -G t I IX) A oX, J) 6000 CONTINUE RETURN f

H

6(1, IR)*A(IR, J) TT = T(X,IR,M.NG, END ROOT. CN, AFAC) DO I = 1,NG

(39)

3.

Subroutine COSTING

Abstract:

Subroutine COSTING is called from program TRANSHIP.

It calculates

the weight and cost of all transverse and longitudinal members, as well

as all welding and labor costs.

The overall weight and cost of hull

structure per inch of length are then determined.

Terms specific to this subroutine:

FORTRAN

Definition

Ma thematic al

Term

Symbol

COSTS

Combined cost of platin, labor

and welding per inch of

ip length

C

F1(H)

Weight per inch of vee weld

f1 (h)

F2(H)

Weight per inch of continuous

double fillet weld

f (h)

F3(H)

Weight per inch of intermittent

double fillet weld

f3(h)

F4(H)

Time - size factor

f(A)

NSEAM

Number of weld seams in inner

N

bottom

S

PLCOST

Total material cost

C

m

TWELD

Total cost of welding

C

w

WBAYTOT

Total weight of hull structure

per inch

S

XLABOR

Total labor cost of erection

and fabrication

_Cfe

(40)

-31-SIJ[O!lTINE COSTING flM'IO!'!/LANEL/ T4ANLu),PLATLABC15),LO4GLAH(b)

O11NSI)N CSTA(jn) ,S444IlO) ,SAALI1fl) ,WTL(17)

C)Mù,COST/ kR4TE,DOLLPPPICE,T4F SO A1E4 (7) *31, SAR5 (j) 4fl . CST4'J(I)rI1IEtNS (ILJ)*4,. SI /XLHOLO*SAN-A(I)

CST(114(AA11(I,1,11)T4EENII(Iù)2. /oA

C1IO' /.' 4,&.NDEC(S,

'Lf)TE,P5L,PHM,HMAIN,HNEIJT. 1ESL,TEST,SESL,

i

SST. ASELS.FAM,XNG,HATC,THIK1HIKK1,EHR,EFB1, 7

EIrWFr41,LOJR.14E-,THIKx,bFhR,WEd.wEbi,xIx,CH!.

i

XIS, AI4EAT1 , SlALj., S'400L1, WHLj,SIJMAfEAP,SUMANEAL .SUMMP, SIJMML,

' ANEA,5UAA,SM3M,4PL,ANEL,BAP.4NST,NSThTR, 5 WTB,BAYIN, SIMJT, SJM,JT,LLGX, YIf)L[,W3AYT, XLAOk,PLCOST, 2 AT(T)rf3(WEB1(I,1,j.)*I).036)*TWEE4M(I))*2, 414(1) r 44&AT1U,1.'C) TrEENM(ID) SG

* 2,0

/44 SIJIIWTP('C) r SIJM4TR(()

* 4T'( I)

SUMCTrS(JMCS1*CST(I)/A4 SIM5TS(JMSTAN*CS14N(!I/AA 14F StliIWMATrSUMWMAT*W'IAT(I)/AA C'NT1'IJF 7AVF.PPEAfl,I)RAFT,ALDAO,XLHOLU,XLPANEL,AAA,H,AA,Ab,DLTESL, SUMSMA1,lP,LL,IL,AREAL1,ASX, SMAT1,SHOOTI C 4AVI4 .,4AyTR SJM4T9() L4GIIUTIINAL SEAMS C')M')Q /E/ INI, E DT

(1SIflJ ¿P(17),SESL(17),SEST(I7),TESL(i7).flST(17),XN((7) ,THIK(

1j7>, 7,5) APEA()7,5), THEE () , SMJj(j 25) ,SM000j( j2 2 .5) ,1H1K4 (j2,2 WILflj SELf1LI Burri 0.0 lT ISr 0, _{TY1 2j} I j. IPLATE

Ir(I.T.ILG) 001022

3).ErBB(12,2),ONEAT1?,2,5),WEdI(12.2.5),SNALI(1/,2,5),S?IODL1(17,2 AMr4FATM/36fl 4,i),4L1(7,2,31,#REAL1(1.2,3),SUMAP(5),SLiNANEAL(5).SUNMP(5),SU RUITS(IlF1fT-4I('C1( 1,1,0) )XN(-4rl

cSL(T)FA!*F5(19I<IU(I,t,1))

(7),TI1(17,5),C00E(i?),ALOAD(17),Ut'SMA'(5),SUMSMAL(5),SUMSMA(5)

_{Dr4SI'J ,LD(17),-1jT I S)) /),wI.LUL(7)}

4(L) I ) M.TlI(K1( I.j,Ic)SG D!M(SION CSL(j7) ,CST(i7. ,9M41(l7).CSLL(17)

fi (H)r (fl,2HH.,04*N)

3

lLo(I)

NSEA'I* Fi) 14151Cl) 1,1,1))

O Ti 2

VN)r

2f(3*.4

2

IF(I,I,ITR) 00 Ii) 24

1(1.4)r i.07*H*N IT I -!LG

F4(M)r0,fl4RIr(I-)

TSL(I)I4EENH(IL))f5(1lI'Ckl(I,i.l))

2 V5(I):fl.fl4*SORIF(H) IItTSIl)rIWFfrbi1lIr))*F1(T41.thi.(I,1.i0)/360

*.

NT:1)4r( i) 0N01 p T )rT,.EN'l( ¡1)1*

)-lIK'Cj( I 1,5) ?,ti

.50 4FL1j( I )rF1(T4!KSj(1,1,1) 1 *7. S IM S T IN r O 4 ¡P(J r I-ITS 51-41 :SUMS8RE O. ;ri

ro (7o.27,2a.29,300 ILr

P)IIT 998'I C AILSE 9955 FIRMAT(IHI,*

Cisl IRAût1.N ,i .

PLATbS *,/)

fl45

25 e.JTTS(I)rkFLUON

F1(T4IIcI5j(1,j,l))/340

3.ii

CSL(I)*

F5(I-4I'C'cl(I,l.i))

.$.1415*HfLUOR 4'l' T rSUMCSTr 50MG SL rS JMCSLLrO

4IL(I)

3,1416. 5.I,*IHTSK1(1.1,5 .50 W)4YTi1 rfl,fl WELD( II rFj(T41<Pcj( 1,1,1) 1)0 1 1< :j,KP4NELS G) TO ¿5 IMSIM ('C I :fl I C 01 FLOOPS C O 2 1 = 1,114 2 4PLt)(I) F2(THI)(51(I,j,1)

) 2,

(HEAM+XG(l)*HfLO04)/XL4OLU CSL(l)r7, JXLHOLG*F5(1411(l,j,1)).HFLUOH.01-4M

ICT,r,T.I11) oü ro5

, TL ( ¡ ) r I i)

4 Jrj

414(1) r 4REATI(TI,<)*NEAM*S1. C GO 1(3 25 NT FLI1014S

I( 1.10.1)

14( I )=4T41 I )/XLHOLI)*2 ¿q

4SL)(I)2.n

F3CTHI1CC1(I,j.j) )*X5Sl*(iAM*XNG(j)*SFLOO4)/XLHOLD 11(1,1 Q2) 414f! )r,4T9( t ),xLNOLD* 4NST C5L(I)rXN5T/XL13LD*F5(I4I5h1(I,1,iI(fL03R*BEAM

I)

f, 5 I'( I .T, ILG) GO Ti 12 11 qT1i( I)r(A4EAri(I,1,4)*(34M -HATC4))bU/AA C 4TL( 1):Ofl SO T) ¿s 01 LUD CT I rF4(AE4T1 1,1,1)).BEAMNAtCN) /44 ¿9

4L1(I)

f2ITHIKS1(I.j,1)

I AT(1)rF3(NF41(!,j,1)*4,Û6)*(8bAM_HATCk( CSL(l)r oNoT HILJOR*F5ItH)KlCi(I.1,1) ) . T) b i? 1 -1L1 C ASSL)'1I D ANEAS

NILfI) r IHISK)(I,1,i). 4FL0014 * G4OI.S,i

_4TL(II

r T.j51Cj(I,l.l).

ifLl(04

cio ri 2"

S 454 C 4 ) ZiO C NI LONG SARFA(3)20 31

W)LO(I)rf3(THI<51(I,),I)

I GNN r

(41)

CL 1 t s NrT. 4FL JQT.F (1H ¡ (Kl (i,) 1) t'-j

LTTr

J1T1 JUT1S( 1) 4Lflj. r -LrIl RLfl( I) S'Mc3I rjMCSL*CSL( i S41=ÇAREA1+WTL(I( 71

C)Tiu-

riT 999

99 FTR4T(/, *

SOAP WELL) dUTT 'LJ AN(4OUHS Lit

l/tNp *,j)

¿Y

47'9 ¡rj,l5

A79 PiT 953

PLATLAR(t),WELD(H, pUTTs(() ,CLH),WTL{i)

-'T 946

,WELUI,4UTTI,SUMCSL AREAl

'I9T oi

3

FTlMat(//.*

TRANSVEI4SS

.,/

Ld/IN WELt) MATER tIUL MAN4OLU-S ,/) jrt

¡i.lfl

H55 P9191 953, TRANLA4(I), TR(i),WMAT(I),CST(i)

P4iiT 946,SUMWTR(t)jSJM4M4T,SUMCST UFC( LOG(TUDINALS 3)

5J,

r 3

¡LU CLL C t) rFl C AREAL1 I

l'i.)

S IM r St Lr S UPIC S L.0 SLL(1i HL'L( ¡t r rNr,( i t T.114K5( ¡.1.1> SARP4L o)

AALt( I l5U( ¡1 *S9

SIIMSAME SUMSAtE .SAMEALOI) 9ELOI i IELOL1 4ELL?L(I)

i

()ONIThJIJE i 3iN1i"iJF kt9T 993

)fl

Rr'Mt5T(/* f1ECi

LONiTUDiNALS //

WELD MAN IW')U9S ç)r 33,ir3,ILU 53 PIJ1 93, LUNGLAtt(i)W5L1)L(I),CSLL(i) SA1PAL(i) P1NT 9r6 ,SfLOL1,SJMZSLL ,SUMSA4E YL.0 r SAtEA1+ SUMArE wMAYT1 r AYLP + JEJAYTR PLC)ST r 'OLLP W-(AYT3T = (rELU1*JELEJL1 + dUITi+SU'MAT).P4ICE

T4ELD1 53.TitELI) St PAL:SUMCST+SUMCSL4SUMCSLL 5LAt)WrSiG'lAL1RATT SJMS IAN OSTS

PLCOST + XLA3O

*

Pi

2, IdAYLP, W44YTR, SIM4L

P9iT 9951 ,ATE,DDLLPP,WPRICE,SG

'l1 FWPAT(

LAIIOM RATE r*,f55,.3,* L)OLLAbS PER

(OUSt

*/

j

* PLA'd COST r*,IIU,3,* ÙOLLAS PEW

POUND *1

2

*

*fj)

COST r*,fjO,3,* IJOLLAbS PL)R POUND *1

3 * JE9SITY r*,I1J,3,* POUNDS PE Cubic 1NCH*/) i FORIAT(//. WI, LONG,

.,Fi!j,2/,*

91, TRANS r*,Fjfl,2/, 1* SUM OF MR

.,F1O,2//)

933 flRPAT(+

*,AA ,Ej3,3 )

9q6 rOI1AT(* TOTAL *,5013,3)

Pill '(1 413, cOSTS, LCUS T

XLAEQR, TItEL))

FORMATO i41, TOTAL COSI

*,F1O,2/.CLJ5T OF PLATING r.,Flfl,2/ 1*COST OF LAROR *,F1O,2/*CIJST 0F 9OLIJiNG

,flO,'/)

343 RETURN END Subroutw G(ESIING Enter Fi(E) = O. Z*H*H -f 0.4tH F2(H) = O. 283*H*H

F3(H) = O. O7*Hp3.j F4))-!) = O. F5(H) = 0. 04*/if GNNT = XNG(J) - GNØT ONT = z SIJMSTAN = O SAREA1 = SUMSARE

O

Print

COST BREAKDOWN"

XMANK = 0.45 SUMWMAT = SIJMCST = SUMCSL = SUMCSLL = O WBAYTR = O

DO K = O,KPANELS SUMWTR(K) Trans ve r s e s _{DO I}

= lITRi

Yes

J=l

WTR(I) = AREATI(I,1,K) BEAM *SG

(42)

*0. 036)*(BEAM -lIATCH) CONTINUE ID = I - ILG

_f

Assumed Areas

CSTAN(I) = TWEENH(ID)* 4*SG/ XLHOLD*SAREA(I) CST(I) = F4(AREAT1(I, 1,1)5 TWEENH(ID)2/AA WMAT(I) = F3(WEBJ(I, 1, 1)5 0. 036*TWEENHaD)*2 WTR(I) = AREATI(I,l,K)* TWEENH(D)*SG2/ AA

L

WBAYTR = WBAYTR + SUMWTR(K)

Longitudinal Seams

WELD1 = WELDL1 = BUTTi = BUTTS = O

DO I

= 1,IPLATE

X BEAM = BEAM/360 BUTTS(I) = F1(TFIIKKI(I, 1.1)

5XBEAM

CSL(I)

BEAM

F5(THIKK1 (1, 1, 1))

NSEAM = BEAM/72 WTL(l) = BEAM*

THIKK1(1, 1, K)*SG

WELD(L) = NSEAM*

Fi (THIKK1 (I, 1, 1))

s

BUTTi + BUTTS(1) WELD1 + WELD(I) SUMCSL + CSL(I)

= SAF(EA1 + WTL(I)

L)

ID =

I

- ILG

CSL(1) = TWEENH(ID)* FS(THIKK1 (I,!, 11)52 BUTTS(I) = TWEENH(ID)* F1 (THIKK1 (1,1, l))/360*2 WTL(I) = TWEENH(ID) THIKK1(I, 1, K)*2*SG WELD(I) =

Fl (THII(K1 (1,1,1)2 IPLT = I -ITR

L)

WTR(I) = (AREAT1(I,1K)* SAREA(4) = 10. (BEAM - HATCHfl*SG/AA CST(I) = F4(AREAT1(I, 1, l))* SAREA(3) = 20. (BEAM - HATCH)! AA SAREA(2) = 30. WMATU) = F3(WEB1(l,l,l) SAREA(I) = 40.

SUMWTR(K) = SUMWTR(K) + WTR(I) SUMCST = SUMCST + CST(I)/AA SUMSTAN = SUMSTAN + CSTAN(I)/

BUTTi =

AA *TWF

SUMWMAT = SUMWMAT + WMAT(I)/AA

(43)

5

BUTTSU) = HFLR*

F1 (T1-IIKK1 (I,!, l))/360n CSL(I)

F5(THIKK1 (1,1, 1)) *u*HFLØØR WTL(I)

r

5 THIKK1 (1, 1 K) *SG

WELD(I) = Fi (TI-IIKK1 (I, 1

1 ))*Z

L)

WELD(I) = FZ(THIKKI (1, 1, 1 ))* ¿*(BEAM + XNG(1)HFLØØR) /XLHlLD CSL(I) = 2/XLHØLD*F5(THIKKI (1,1, 1)* H FLR * BEAM WTL(l) = 0.0

NT F1oors

WELD(l) = ¿*F3(THI}ÇKI (1,1, 1 ))* XNST*(BEAM + XNG(l)*HFLØR) / XLHLD CSL(1) = XNST/XLHLD* F5(THIKK(I, I, i))*HFLØtR* BEAM WTL(T) = 0.0

OT

Longitudinais

WELD(J) = FZ(THIKK1 (I, i

1 ))*

GNIZtT

CSL(I) = GNTHFLR

F5(THIKK1(l, 1,1)) WTL(J) = Ti-lIKE! (I, 1, 1)5 HFLIØR*GNT*SG

NT

Longitudinais

WELD(I)

F3(THIKK1(1, 1, 1))*

GNN T

CSL(I) = GNNT*HFLOR F5(T1-!IKKI(I,1, 1))

Print Pertinent Material

for Piales

Print Pertinent Material

_\

for Transverses

S____________

Deck Longitudinais

WBAYLP = SAREA1 + SUMSARE WBAYTOT = WBAYLP + WBAYTR PLCOST = DØLLPP*WBAYTOT TWELD = (WELD! + WELDL1 + BUTTi + SUMWMAT)WPRICE TWELD = 1.03*TWELD SIGMAL = SUMCST + SUMCSL + SUM CSLL XLABOR = SIGMAL* HRATE + SUMSTAN COSTS = PLCOST + XLABOR + T WELD Print Pertinent Material for Deck Longitudinais /

00 1 = 3,ILG

Print Overall Weight

_{\ and Cost Figures}

cSLL(I) = F4(AREAL1(I, 1,1)) SUMCSLL = SUMCSLL + CSLL(I) WELDL(1 = XNG(1)*THIKK1(1,l,l) SAREAL(I) = AREAL1 (Ii, K)+ XNG(I)SSG SUMSARE = SUMSARE + SAREAL(

(RETURN)

WELDL1 = WELDL1 -f WELDL(I)

L

CONTINUE IPLT = i Bilge

t

O. T. Floors

(44)

4. Subroutine DOCUMENT

a)

Abstract:

Subroutine DOCUMENT is called from program TRANSHIP.

It prints

out tables of scantlings of all structural members, the values of which

have been computed and stored in the memory of the computer.

(45)

-36-S T INI

COM'rN /11/ 1,:,l:;:L,tLATL,

L,(!.CÀIC,-fNLUT,T SFST 1' j1M. \ICR, HI 1C THI THI 1 2 EFrw,Fr1,HFLJ2r.. CHi, 3

XIS, A)-A Ti.

CODL1 WLRL1 ,SU.AP ,SIJM1 I

LCMP, SU14NL. ' AFCEA,SU'ARFA, J' 4,HPL,XIANEL,W.-)YLP,XNST,I,.II. I 5 WTR,DAYT,Si>;1Wi.,IJC1T,LLCI. YIE>.D,WHYTR,XL'(,PLCOST, 6 COsT. CC)'11' IN, G4JT

, K.)C)j, SUlN.C', SIJISCIAL ,IUNSMA ,c°OF.

7AVFl, PR))A1 , J(1jT, ALL1AIJ, XLHC1I J

Xl 'AN), )A , RhE(, A 3, DEL ESL, 8 SUNSMA1,/P,L,iLR,ARALi,ASX. 51iAjj,OL)T1

cni)o> // ÎICT,

DINf-NSjON 71(17) ,SERL(17),SEST)17),TESL(17) TES) 17 ,XNO( 1),t)1IKK( _117>,

II(XR1(J

',i),EFFJ(17>,1FW(1ï) ,EFF)C1(17,2,1)5)W1(li

2,2,5),AREA(1 /,5),T L'CH(5),SHATI (12,2,5),SNODT1(17,2,5),'If41<X(12,2 3). 12' 3) AREATIC j, 2. '1 1(12,2,> SN/L1(11,2. 1) ,SNO[ILI(17.2

4,3>,wERiu.2,3C,II\L1(7,2í>;uNAfiNAp),sUMAPE L).SCjC)iP(5).SI)

5N14L(5) ,SUCIA'LI(5) , -, JiI'IU'1(

) Pt(5) ,1LPArJEL( 5) ,S)(rIJTR(,) W hoi?) ZL

PRINT 20 Rn PRINT 2110,).> THIKKI( 1,J,X).J1,2),Krl.N) PRINT

j,)) THIKK1( 2,J,K),J-J.2),Kni,N)

PRINT 2,(

THIl<1( 3,J»)*Jn1,2,Xn1,N)

PRINT 3,1) TNIKK1( 4,J,1),J1,2>,Xn1,W) PRINT 4,)) THIKKI(

,J,4),Jn1,2),K:1,N)

PRINT

TH)KK1( 6,J),Jnl.,2j,KrC,N)

2RI).T 7 I>NI

i3.(

IHIKK1(j1,J,(),Jrl,2) Krj,N)

PRINT

8,)) THIKICI( /,J,4)i,2),Kn1N)

PIIRT

9,)) THIKKI.( 8,J),Jr1,2),K:1,N)

PRINT

in,)

TH)I<K1( 9,J,K),J1,2),Knl,N) PRINT

il,)) THIKKI (i,J,K) !Jnj.2> ,Krl,N)

PRINT 77 PRINT

i4,(

THIKK1(12,J,4),Jnl,2),K*1,N) PRINT 115,1

THIKKI.(.t3,J,K).Jri,2),Kr!,N)

PRINT

16.)) THIKKI(i1,J,K),.Jnl,2),Kn1,N)

PRiNT

17,)) THIKK1(j,J,K).Jnj,2),Kn1,N)

PRINT 21 PRINT 14, (C

IERi (1,J,'<),J

1,2). K r inN)

PRIllI 15, 1

IERi (i.J.).J r 1.2>. K r j.N)

PRINT 2, (C WEB1(3,J,K)J 1,2) , K

i.N)

PAINT 3, CC

WFfl(4,J,K),J

1,2) , K r 1,N) PRINT 4. (1

WF(i(5,J.K),J r 1,2>

,

K r 1,N)

PRINT .

I)

Ed1(6,J,K),J = 1,2> ,

K r 1,H)

PRINT 7 PRINT f), (C WFI31(7,J,K),.i 1.2> , 1< PRINT 9. (C WEJI108,J,K),J r 1,2) K r 1,N) PRINT 10, ((

WF1(9,J,K),J r 1,2)

, K r 1,N) PRINT 11. o

WEj)10,J,K),Jj,2),Krj,N)

PRINT 22 PRINT 14,

tI

(EOLiO i,J.K),Jr1,2>,K1,N)

PRINT 15.)

WERL1( j,J,K),J=i,2),Krl,N)

PRINT

2,0) WERL1( i,J,K>,Jnl,2),Kni,N)

PRINT 3,

C) WEOLI( 4,J'K),J:C,2),K1,1)

PRINT 4, (C CECILiO 5,J,K),Jrl,2),Knl,N) PRINT 5, (0 WEICLI) 6,J,(Ç),Jr1,2),Knj,) 20 VP-1A ICOHI. .r?írE lOLO *1/. 11E TRI,ES5O:S

//,181,

1* (N 4AY 11V II,ft1 2* .00ATTON PA'>EL I PAN7. 2

*1/)

21

FOR(IAT(1N1,* CLITEK HIJLI)

/1* W-O HELORTS OF TRINSVWS IRA11S OR S

1HAPES *///,:IX,

I'> .Y

f HATCH 2 // LOGATIU" PANEL i PANEL 2

//)

22 FORMAT)

/1* WI] HEIr,HTS OF LONGITUiJINAL IJIRDERS

i

*/J/,1.5X, I4 WAy OF HAtCH 2 *//* LORJ1 Tul) )'A\i!.

j

('AREL 2

*//)

2DO FONCIAT(* HOTTJ(1 3HLL

.,6(Ff'.2.4X)/

> 1. FOR(IAT(. T,', PLATI'JG

*,5(f62,4x),

C 2 )ORM4T( ClERK 4 3 FO4ilAT(* [1E'K 3 *,O(Fó.2,4X>, 4 FORNAI). ICECK 2

ò,6(f6,2.4X)/

FûflfiAI(* 1100K

j

*,6(F6,2,4X(/

7 FORNCAT(* SIRC SWELL./ 8 FOROAT(*

I,r>.TO 1

*,6(F6,?,4X)/

) 9 F09114T(* 4 13 3

*,ö(F6,?,4x/

10 FOR'IAT(* 3 TO 2

*,5(V.7,4X),

11 FI.CRIIAT(* i TO 1

*,ò(F6.2,4X),

13 FÙR1IAI(* BILliE

,6(,2,4X),

14 FORNATC 0,1, FLUORS ,6CF6.?,4X)/ 115 FORMAT(* N,t, FLUORS a,b(f-ó.2,4X)/

16 FORIIAIC* çI G)RflE5

.6(F6.2,4x)J 17 FflOIIAT(. Nt, (iIrOUï'TS

',o(Fh,2,4x)/

P1111)1 18 lA F021111). ji-li,-S)P- ROOPEITIEC *, /1) 77 FORMAT (* UUuLji O 311 -/,') RE TuRN END

(46)

Subroutine DCUMENT

Enter

Plate Thicknesses

T

= 1,15

Print THIKK1 (I, J, K)

J =

1,2

K = 1.N

Web Heights of Transverses

Web Heights of Longitudinal Girders

I = 1,6 Print WEBLI (1, J, K) J 1, z K = 1,N c'i co

I

1=1,10

Print WEBI(I, J, K)

J = 1,2

K = 1,N

(47)

5.

Subroutine EV

Abstract:

-39-This subroutine is called from subroutine

_{GRILLAGE, BNMAT, XLAD.}

It calculates the characteristic equation

_{of the A matrix.}

_The

Cayley-Hamilton theorem is used to reduce the characteristic equation of the

A matrix.

Terms specific to this subroutine:

FORTRAN

_Definition

Term

A

_{A matrix}

CCC

_{Coefficients of characteristic equation}

in descending order

D

_{Identity matrix}

NG

_{Order of matrix (Number of}

longitudinal girders)

(48)

SUBROUTINE EV( A, NG, R, CCC)

DIMENSiON A(9,9)?019) ,U(9,9),AN(9,9),BN(9,9),CCCI1Û),R(9,9)

NGI;NQ-1 DO 90 ¡j,NG DO 90 Jj,NG IF(I.J) 1,2,1 1 Ù(I,J)rO,fl GD TO 90 2 O(I,J):1,0 90 ANII,,.J):A(I,J) DO 200 Lrl,NG 0CL) rO D DO lin (.1,NG 110 R(L Cr8)). )+ANCK,Pc)/FLOATfr CL 00 120 Ir1,1t Do 120 Jrl,NG 120 BN) i ,J)rAN( i ,J)3(L "U) C if (i-NOI) 23CC,210,3U 210 1)0 22E Irl,NG DO 220 .J1,NG 220 Rfl,Jl:BN(I,,j) 230 nO 130 L"l,NG DO j30 ,Jrt,NG AN) i,J)rO,0 DO 130 1,N(i 130

AN) i,J)rAN( I,J).A(I,)')BN(K,,J)

200

CCCCL+t)"lflL)"C-1 ,)

Cccli.) r j

1)

RE T END

Subroutine EV(A,NG, R, CC)

HNo

DII, J) = Enter NG1 NG -D) I = 1,NG D J = AN)I, J) = A(1, J) D L = I,NG B(L) = O I, NG D(I,J) = i D K = i,NG ¡

t

4

L

BN(1, J) = AN)I,J) -B(L)*D(i, J) B(L) = B(L) + AN(K, K)/FLATF(L) I, NG

f

DcJI J

= ING

t

(49)

L

J=ING

AN(J.J) = O

- DO K = 1,NG

r

AN(T,J) = AN(J, T) f A(I, KBN(K, J) CCC(L + 1 = B(L)(-1) DO J I,NG DO J = 1,NG

j

-DO RU, J) = BN(I, J)

(50)

6.

Subroutine FRAME

a) Abstract:

This subroutine is called from program GRILLAGE and TRANSHIP.

It calculates the required frame moments and defleclion, using the

method of slope deflection.

Three options of calculation are provided

as follows:

NOC =

1,

_{Unrestrained frame deflections at girder intersections}

NOC = 2, Influence coefficients

NOC

3,

Frame bending moment distribution

h) Terms specific to this subroutine:

FORTRAN

Definition

Mathematical

Term

Symbol

A

Slope deflection matrix

BB

Fixed bending moment

D

External head at innerbottom

DD

External head at the node

DLL

Height of node above innerbottom

PD

Deflection array

T

Joint rotation

e

TM

Transverse deck beam moment

V

_Loading

WD

_{Uniformly distributed load on a given}

level

WDL

External load at the node

-42-WDLX

WDLXT

WDX

Fixed end moment corresponding to

external head of water

Fixed end moment corresponding to

external head of water

Fixed end moment due to a distributed

load

XJ

Stiffness factor

K

XLF

Length of stiffener

(51)

C9)CX:U*(0

1V

tjT'tI

01109 üG

0U:(T')V

IJrhHX

66690100 (03Sj0N)J1 ZI1

AVEI OV3HW1fld 3

u u(i)UVoÌv

(I)CX4i('?.)Vr(?Y!

1ttilS10N)3l

LIL)l pU

0(1

ccV)rx*(L)rx).]c?1v

oi. oü

(1NaO)3I

(L)rxuc:U

01V u C 1)1X1üM

U(t'L IV

0(TT IXUM

UCt''t)V

C1Ue/)I)UO't/(l)31X)..(C

)Jx/.( )('a(C I

10rs(

1)Ix1UM 99V DIO

1)TlV

(/?.0 0U.9/

I )(11.( )3ÌX-09/o

..(T

1x)i..0 CLflK

T (1)CX.i.O (CTIV

/

4CI)II)IU0(I)XCUM

C (L)cxd4)t)rx)

.(T'T)V

'0T/,(C)31X(I 1i1Mu(1)XO) xI)9W t40I13EiC300 3d01S 3

1LC 9i9

00 flNIJNO3 666

(t)J.v).)93

660)' UflNI1(U3 6 OÛNIJ.N(J3 ThOU (1HX( I )A (NN)31X -I o CI OU CII 11(1)0(1)31 CtW)1X1(1M(o)x1(J'(TW)XUM.CW)XOM T 110-( I )31X+(Jul I IOU (C)UM(I)1UM II0t tUUÇ 01 1,0 9.NW (I IJiXuC I 100 (1)UI- CC) (I) iUM

>1UNN 6 0(1

C0f11(Í1I.Mu(I)uM T ?,OOt' 0.100

(1l(1'J.)U)JI

(I'OU'JHOIJ)lI CI131X'11QV1U (('o3DEJN)3C

1-Irr

(t)x1U* (L)XUM*CT)U(T)Ii8

0.1N'LrI O(P0 OU / C 1H.Bi-9i). I O/I 1)J1X. (lId :( 1)118 U+rNr0.LN AVO HDAVH 1303 No JV1OD1V'l 1NW(j11 013)011813 08X13 3 -IONI !Ñû1 lu/O (I)31X/Ci)3V13S rU)rX OCiU0 OT'T1 0019 ca

dflrh1iu3 69

(080110N)3I

eHX(1)13r(1).f1X c/XuO

(00Fi0N)31

T43(1flN'Çz1 6019 OU

'/Xrj

(t3-'UN)J1 ûI-s-H X HH13NX AVIJ i31V11 3 0-r(T)OVC1V OONI!11U3 6669 (1)lUu(I)090IY T311 UÚÔU Di 05 >h3do,'iuI 10çÇ OU dON Ii NO3 Dt

CIH1;- t+X9-iN >h:IUJN

I 1W)1X1(1M( IX1('r(TW)X(IM()'I)X(JM +

t

'ç,

(N)(1.()O9O1+ (N)1/?/(1(N)U)r(1)(t)5h3

(L)Ud'C1)hX'(/)('/)ç'I I '(ill!! '(f.)TlI'CçT)f'X T

'(/11' (9)OP,( '/10 (/)d'(lT )0V01V'(Ç1).'jX (Ç,T)31U11 NUlSI4Jh)0

(/)i.j NUIS

&:IU

CL)NM(0T)1x'CoTCJM:n1,N N0IS11d,-,IU

)4310oir

CT 00 19'C,II.XX N01b1131'dU

11//(?r4U(

t

Ct)jNX NQI':NJ.JI1

rU)r(',1+

(ttjjiJI3IJN)3I

(L1.))W,X 13,3I,:JII0 (((1301J0N131 1))' /13/ NO-0i93 Ct)x1 (/)X(iO'C1)*Cr(t))5

(dIX' (('10 NOl'-..11l0

(1.1)10 NOI,h,r.i.l1) lof) 1303 11(1V1fl31V iN-3Iilli, ONION-JO 01X13 3 (flTUo'CUTYlfl;-'1:L)11 N0XN-Jh)U 3flí1jNUJ 0009 NI0d T 131X/1 ii 3VlI; C I )CX '3013 'Ç) 0\

'(( '/'flN 1'(tM

31

'--.

(52)

1/ F SINJNoW 1i (or OdVfldli

4',

)iVW'u4 S l(It' ININd i -1 N r N (4'('T9!)J.Yk1.FU. i19P 3nNl!NO3 (Hi.1 CJSRV( I )1X

TTl' 0)00

((l)IHX'j'1(HW1I3SflV)31 Te 0100 raÌN'l'lIo'O'u'r>i).(1 Z'Hk.L' UO0 1N11d (1) i:irk-,J.

(ILflX/

T 9/Ç.Z.(1 )lQM+'?/.IL4Z.( (1 flçiM+ I )0)-(7( IA+

(')i'9>1)FiXX)(1

liii t (I )3lX( 1./( I )(1M+ i c- i 91 9r4 (HN/r)# (1) 3iXZ

durTr

-0 1ó°ç, 1NIld

LSSI 'riI1 itOl: OC

-rN4SsI issI

NO iuNs s I t, k) X: C p 'ç) W XX ç i, k) X I Ç'P) W XX ('vWXÇ.'P)WXX £rNXÇ'Cr)WXX ï W)) I ï ' I k XX ïk1X:(C'1.)WXX -flJt1N(J3 099v :(fli,ljNCO TÇ92 ( KNI ) .1SHV

(Fir) IWA '1ro

tÇ6ô T99S ' T99ç' ((r) LuX- (Nih. I .:,f1V I T h r t99c Ut 0

(U'o3f>)j1

Z' I4)4L'IJÇOr,INI)Jd (1)Fe1kiJ. ¿ç9 (X-Z t (i) d- (i I k.L' (1> ol (T ' ri 1) . I XVI JV3IIrinO 80.1 3 t X7 t t Id-(T )hj C T)kij. ç9ç' 01 00 CX'J.T'Z).jl Trir (i-Z).(Ç),j-()k1('1.)I.,L

(3'1N'Z),I

Z)(P)d-(I)ii(T)i

fl'jN'7'(iNV' t',Ei'OiPI'u)JI

TICt, 01 OD C

i1'0i't'n0't'ir' 1)31

I )flVclyZ

1)A+( I)0X)'( oui Cts9ç

T. dU N r NH-Cl )( CIX

hUX-(I,FZ

TC FiGcç 0). Ut)

(X'il'

2)31 01 0)) (T (SF1 uro S ON) 31 (I) Ei

T-rrN

'':('

99ç OC )WI ).c,IN IIT:dwN TP9 (/'Çi'4 iOoN

/wO83 NOIIVO01 ,',''.ciNikuOuI

,/)1VU80. TT,O), (/'ífF. 1r.u(OC d2'tOl.', NOdi

NOIlVODi .'/'iitriDN

/).j)k)803 O5 (ltd-t I IA") t I/u (T10 1)31

(ç)dITtTA

('(J:]'0ri8N)31 )Z)d-C( r > (t)IVÛ1T,)rflX.(T)A'-IT IA I )UuyuV "I I (A T-008N+ t l'TçOS INIdd

.3O'T:t'099

00

(t,3ç)IYhit)O4

OçU

(I

jNRWQW WY313 )'DRCI 35)43ASNVF,i' '/ )V,.)jQJ 8jE& POOSh 1IlHd 0656 N-]N013 l'Oh(I) KO sNot .LV'(ri1"3 N3WOP( O

ç/c0'çt,'ç,o (-D(jtl)fl

ç/?#()H.o+)

tIiItrxriIu-ux

(Ï03JON)31 X-r) U I iiX Çt,UX-St'I'lX"(t,)WX uhX-'WX-C I (c)i'lX C T ) NX "dt lix T IÇ'),-I-(T)WX' I T)iX (p0301u)JN)31 (u1E3OoflN).11

(t)d(T)JrIïl'Xri r(T)liX

(6)lXlOr)6)XUul+

C (r)Lt(r.u'?)4()rxd.udrrp'lx

(6)X1)JM-(6 )XCM ( (P)j+(l'J

1I 'dIr (t,l'X:UIliX

11.) Xli;) r. I -r -r '' 'r rr"

(L)IX1ONI+(LXUM 'C (T)i+.'C1»(L)i X.h(.d:diX

Nr

I i 'lflD'ivo IN2WOh .ir roc

a r

(1)1 Ti

rf!'Trr

T/.09 0 I (J Tr1 IL1)') 0G NOIlVifl3'1VO ïdOlS 3 (w81EI 'Ti! 'dj.I

'T.j1'uJ''

¿'rN'l)

flNIiPi4lV3

F)oIiVlDulV3 XIHIVii 3SuiiANI 3

C9)rxF3+(Ç)V('V

('(ÚOsJ.CN)iI C 0T)rx3. t'r( S'SIP U" CL S IV

O)

rYS) U" (I '51V C UT) rx.3 d"( Ç, PIP

C ç)rX-ì.( t' 't') t(t''t, Iv

(i

u I (6 )rX.3 'd" t S 'PI V 0" CC P 1V U" I i J') i L I V C 6 1 r X 3 ¿" (t' £ 1 y t t,i r C '