Bulkcarrier ship preliminary design with the aid of optimization method

0

TECHNISCHE HOGESCHOOL DELFT

Afdeling der Scheepsbouw- en Scheepvaartkunde

Vakgroep ontwerpen van schepen

BULKCARRIER-SHIP PRELIMINARY DESIGN

WITH THE AID OF OPTIMIZATION METHOD

BULKCARRIER-SHIP PRELIMINARY DESIGN WITH THE AID OF OPTIMIZATION METHOD

by Dr.Ing. L.K. Kupras,

Delft University of Technology,

Department of Shipbuilding and Shipping,

Ship Design Panel.

SUMMARY

The computer program for bulkcarrier-ship preliminary design is described.

It offers limited flexibility due to the selection of free variables,

pa-rameters and constraints. The optimization procedure searches for

optimal

solution according to the selected object function. The program can

ana-lyze data of the "similar ship", which will be used later for designed

ship calculations. The main flow charts, program listing in ALGOL-60 (T.H.

Delft-compiler) as well as samples of input/output are shown and explained.

Contents

page

Introduction

1

Free variables, parameters, constraints

and object functions

2

Method of solution

3

Program construction

3

4.1.

Analysis

3

4.2.

Data for optimization procedure

3

4.3.

Optimization procedure

4

4.3.1.

Initial and final step widths,

search cycles

4

4.3.2.

Tolerances and accuracies

4

4.4.

Design model

4

Ship's subdivision

6

Numerical procedures

6

Description of the program listing

7

Sample of input and output

Termination condition

9

Recommendations

9

Job cards

10

Literature

11

Figures

12

Appendix 1

: Computer program listing

16

Appendix 2

:

Input data sheets

31

1. Introduction

The present report describes the program of preliminary bulkcarrier ship

design by the aid of optimization method.

The program allows to find the optimal set of main ships particulars due

to the selected object function, constraints and limits.

The program can execute fresh calculations or can make an use of similar

ship data analysis.

The program is flexible due to:

the range of data concerning the similar ship;

the number of free variables and parameters;

the number and qualification of the constraints;

free selection of the object function;

in such sense that maximal number of similar ship data is limited as

well as number of free variables, parameters, constraints and object

fnnct ions.

The program can be used when the number of equality constraints is lower

than the number of free variables. Minimal number of free variables is

one. If other conditions then see the program for

parametric study [16].

2. Free variables, parameters, constraints and object functions

The following tables give a look in the possible selection of free

va-riables, parameters, constraints and object function .

In most cases

the equality constraints represent owner requirements which must

strict-ly be satisfied (within assumed tolerances), and inequality constraints

represent either boundary limits caused by ship environment, or have

the task to delimite the range for search where optimal solution is

ex-pected.

Table 1: Possible Combination of Free Variables and Parameters

1

1

V

either free variable

or parameter

2

CB

either free variable

or as parameter K;

K = CB + 0.5 x

V/117

3 L

either free variable

or parameter

4 B

either free variable

or parameter

5 T

either free variable

or parameter

-

2-Table 2: Possible Combinations of Constraints

I

Speed, Block Coefficient,

Main Dimensions,

if considered as free variable then

upper and lower limits as inequ.

constraints

2

_{DEADWEIGHT, STOW.FACTOR,}

CARGO CAPACITY:

DEADWEIGHT, STOW.FACTOR

both as equality constraints (exact

values)

DEADWEIGHT, C.HOLDS

CAPA-CITY

both as equality constraints (exact

values)

STOW.FACTOR, C.HOLDS

CA-PACITY

both as equal. constraints (exact

values)

DEADWEIGHT

_{equality constraint (exact value)}

C.HOLDS CAPACITY

_{inequal. constraint (upper limit)}

STOWAGE FACTOR

_{equality constraint (exact value)}

C.HOLDS CAPACITY

_{inequal. constraint (upper limit)}

C. HOLDS CAPACITY

equal. constraint (exact value)

DEADWEIGHT

_{inequal. constraint (upper limit)}

DEADWEIGHT

C. HOLDS CAPACITY

-both as inequal. constraints

(upper limits)

3

i cl__

INITIAL STABILITY IN FULL

LOAD CONDITION

either equality constraint

(exact

value) or inequality constraints

(lower and upper limits)

4

STATICAL STABILITY RANGE

IN FULL LOAD CONDITION

if required then inequality

con-straint (lower limit)

5

FREEBOARD

_{either equality constraint (exact}

value of minimal freeboard required)

or, inequality constraint (minimal

Table 3: Possible Object Functions

The choice of the object function influences the proper selection

of constraints. There is no doubt about the importance of this

re-mark. For example you can not select the deadweight as an object

function to be maximized and in the same time ask to design the ship

for exact given value of the deadweight.

METHOD OF SOLUTION

The problem of preliminary bulkcarrier ship design is solved by the

aid of optimization technique, which works on the principle of "better

point algorithm" [14]. This method allows to start with the search

from any point either feasible or infeasible. Equality constraints are

solved every time during the search, by the aid of linearization and

iteration process.

The search proceeds in the direction of feasible

and later when feasible is reached in the direction of the optimum.

PROGRAM CONSTRUCTION

4.1. Analysis

The main flow diagram is shown in Fig.

1. First of all, the program

checks if only fresh calculations are required. If not, then input

data concerning the similar ship are read and analysed in such a way,

that more important technical characteristics are re-calculated. For

such an analysis use is made of the sane procedures and formulae as

will be used later for new design. The experience coefficients are

calculated, which represents ratios of actually to the calculated

values. Designing new ships, the calculated values (e.g. for weights

or capacities, etc.) will be multiplied by proper experience

coeffi-cients.

4.2. Data for optimization procedure

In the next step the program reads designed ship data, counts number

of free variables and of equality and inequality constraints. These

-

3-minimum of building costs

2

minimum of required freight rate

3

maximum of deadweight

maximum of stowage factor

5

maximum of cargo holds capacity

minimum of building costs per

1

ton of deadweight

minimum of building costs per

1

cu.m. of cargo holds capacity

S

minimum stow. factor

4

numbers are required to run the optimization procedure. The start

value of free variables, needed for initiation of the search, are

the mean values of upper and lower limits, as given in input data.

4.3. Optimization procedure

The optimization procedure initiates, governs and terminates the

whole design process. The real accuracy and credibility of the final

solution depends on the accuracy of single procedures or formulae

used, as well as on the final step widths of free variables.

4.3.1. Initial and final step widths, search cycles

The following step widths were assumed, as for:

Three search cycles were assumed. In some cases the method can

re-quire one additional search cycle.

4.3.2. Tolerances and accuracies

The equality requirements will be solved with following tolerances,

assumed, as for:

deadweight

:

+ 20 ton

- stowage factor

:

+ 0.03 cu.ft/ton

cargo holds capacity

: -T 30 cu.m.

initial stability

: 71. 0.01 m

freeboard

: -T 0.01 m

The inequality requirements will be satisfied within assumed accuracies,

calculated by the optimization procedure. For more explanation see

literature

[14 ].

4.4. Design model

We assume the following definition for the "design model". One cycle

of calculation which for one set of free variable and parameter values

allows to compute constraints and object function values will be named

a design model.

The flow diagram of bulkcarrier ship design model is shown on Fig. 2-5.

In the first step the model checks whether ship's speed V is considered

to be free variable or parameter (constant value

). If V will be a

free variable then two inequality constraints for upper and lower

limits are created:

V - V

<0;

and

V.

- V < 0;

max =

min

=

where: V

= free variable, speed

V

= upper limit for speed

vmax

= lower limit for speed.

min

initial

final

- length

1.00

m

0.20 m

- breadth

0.10

m

0.02 m

- draught

0.05

m

0.01

m

- depth

0.10

m

0.02

m

- speed

0.10

kn

0.02

kn

- block coefficient

0.005

0.001

In the next step the model checks which of the main dimensions L, B,

T, D will be a parameter, and which a free variable. If any of them

is

free variable then two inequality constraints are created for

upper and lower limits, e.g.:.

- L

< 0 ;

and L

- L < 0

max =

min

- B

< 0 ;and B

- L < 0.

max =

min

The upper and lower limits are given by the designer as input data.

They are derived either from design requirements or from the range of

validity Of the calculation methods.

The block coefficient can be a free variable, or a parameter, and will

be expressed by a factor K from formula: CEF. K - 0.5 m V/V L/0.304Er

Constraints:

either OF - CB

<0 ;

and CB

.

- CBs0 ;

max=

mmn

or

CB- - 0.84550 ;

and 0.535 - CB 5.0

;

In the next step the model checks if the ship will be designed for

required deadweight and cargo capacity or stowage factor. The model

creates the following constraints,

either:

DWR

- DW

= 0;

and VOLTR - VOLT = 0;

or

DWR

- DW

= 0;

and SFR

- SF

= 0; or

VOLTR - VOLT = 0;

and SFR

- SF

= 0

where: DWR, DW

=

deadweight required and available

VOLTR, VOLT

=

cargo capacity required and available

SRF, SF

=

stowage factor required and available.

The next set of combinations of equality and inequality constraints

due to the deadweight, cargo capacity and stowage factor, allows to

search for optimal deadweight or optimal cargo capacity. In such a

way the deadweight and cargo capacity are treated as indirect

vari-ables.

Initial stability for full load condition can be considered as an

equality constraint:

GM - GMR = 0;

or as two inequality constraints:

GM-GM

<0;

and

GM.

- GM < 0;

max =

min

where: GM

= initial stability available

GM

= upper limit for GM

GMmax

= lower limit for GM.

min

The model allows to create the inequality constraints for minimal

range of positive value of righting arm curve (statical stability):

a .

- a

<0;

min

h=0 =

'

where: amin

= required minimal range of positive value of righting arm

ah=0

= available heeling angle, at which righting arm is equal

to zero.

-5-6

Requirements for freeboard can be expressed by the aid of equality

con-straint:

FE

.

- FB = 0;

min

or inequality constraint:

FB

.

FB < 0

min

where: gmin =

minimal freeboard as required from regulations

=

freeboard available.

In the first case the ship will be designed exact for minimal

free-board, and in the second case the freeboard will not be legs than minimal.

In the last phase the program checks which from the offered object function

has been

selected already by a designer, and then creates and

calculates the value of such a function.

In the proper places of the program there are blocks for resistance,

propulsion, weights, capacities, C. of G., building and operating

costs calculations.

SHIP'S SUBDIVISION

The subdivision is assumed as shown in Figures 6a, 6b and 6c. The

en-gine room (diesel main enen-gine) is positioned aft, and the double

bottom is of the same height on the whole length. The symbols describing

dimensions as for double bottom, sheer, hatch coamings, wing and hopper

tanks, are the same as used in computer program.

NUMERICAL PROCEDURES

Cargo Hold Capacity

Volume coefficients of the whole ship, of the engine room and of double

bottom are derived making use of nondimensional values of hydrostatics,

section area and moment curves, developed for ship forms of Series 60

[24 ]. The proper volumes are calculated by multiplication of volume

co-efficients by a product of main ship's parameters. Corrections for camr

ber, sheer, forepeak, cargo hatch coamings, wing and hopper tanks are

introduced using simple geometrical relationships.

Light Ship Weight

For steel weight calculation, the

method

of Hagen Johnsen and Ovrebo

[11 ]

has been applied. For weight of forecastle and deckhouses

for-mulae of H. Nowacki are adapted. The weight of engine room equipment

is calculated according 'to the diagrams developed at T.H. Delft by

J.B. Polko and A. Groeneweg.

Centres of gravity ab. B.L. were estimated according to the method of

H. Schneekluth [22

1.

Resistance and Propulsion

For main engine power calculation the method of Auf'm Keller has been

applied. The residuary resistance coefficients are calculated for the

Stability Cross Curves

On the basis of the regression analysis carried out for ship forms of

Series 60, a set of linear equations is developed. With the aid of them

the stability cross curves can be estimated [ 15 ].

Freeboard for Ship's Type B-60

The procedure calculates minimal freeboard required, under the

assump-tions:

- the length of forecastle = 0.07 m L;

depth of forecastle - standard;

no sheer in the range of cargo holds, small sheer in the forward and

afterward part;

- no other superstructures included to the calculations.

Stability

Stability has been calculated for ship fully loaded (homogeneous cargo).

Free surfaces corrections are included. KB and MB values are estimated

by formulae derived on the basis of ship form Series 60 data analysis.

Building and Operating Costs

Building costs procedure has been established with the help and advice

of one of the Dutch shipyards (level 1974).

Operational costs procedure follows the model of Tsuneo Kuniyasu[ 13].

It is rather simple and cannot be used for serious economical

con-siderations.

For more explanations and computer listings of the numerical procedures

see separate report [15 ].

7. DESCRIPTION OF THE PROGRAM LISTING

Computer program listing can be found in App.

1.

The sequence is as follows:

5 - 6

- declarations for constraint arrays of the standard design

model; equality - HHH and inequality - GGG

6 - 7

- declaration for basic ship data

8

- declaration for designed ship data

9

- declaration for stability cross curves

10

- declaration for experience coefficients

11

- declaration for output

32 - 132

- calculation of number of equality and inequality constraints

36 -

71

- calculation of start values of free variables

72 - 128

- estimation of final tolerances of equality constraints

134 - 169

- estimation of initial and final step widths of free variables

170 - 173

- declaration for optimization procedure data

175 - 177

- procedure for freeboard calculation

178 - 186

- procedure for propulsion calculation

187 - 189

- procedure for weights calculation

180 - 192

- procedure for fuel oil, diesel oil calculation

193 - 196

- procedure for C.G. calculation

197 - 200

- procedure for stability cross curves calculation

201 - 203

- procedure for BM and KB calculation

204 - 206

- procedure for initial stability calculation

207 - 210

- procedure for righting arm curve calculation

-7-- 7--

8-211 - 239

- procedure for building costs calculation

240 - 244

- procedure for operational costs calculation

244 - 275

_{- procedure for the design model}

276 - 422

- procedure for experience coefficient calculation

422 - 741

- procedure of bulkcarrier ship design model

436 - 509

- free variables and their limits creation

594

- 680

- constraints creation for deadweight, cargo capacity,

stowage factor, stability and freeboard

693 - 724

- output - main ship's data

725 - 730

- output - stability cross curves

731 - 736

- output - righting arm curves including free surfaces

cor-rection (full load condition)

737 - 739

- output - economical results

741 - 748

- optimization procedure

749

- program execution, experience coeff..calculation

1550 - 1551 - program execution, search for optimal solution

752 - 1744 - output.

8. SAMPLE OF INPUT AND OUTPUT

Five different test cases are carried out. Data of the basic ship as

well as Of designed ship requirements are the same in all test cases.

For sample of input sheets and input data (Test Case No.5 ) see

Appendix No. 2.

Test Case No.

1

Free variables

Parameters

Constraints

inequality

equality

Object function

Test Case No. 2

Free variables,

Test Case No.

1.

Object function:

Test Case No. 3

Free variables

Parameters

Constraints

inequality

equality

Object function

:

L, B, T, D

:

V = 16.00 kn.

K = CB + 0.5m

V/117

= 1.107

:

200.0 < L < 245.0;

25.0 < B < 32.25;

10.0 <T 4.-14.0;

18.0 .7 D

25.0;

=

FB

< D- T + s;

2.0 < GM < 3.5

mi

=

n =

=

s - stringer (deck) plate thickness

:

DWR = 70000.0

SFR =

46.0

:

minimum building costs

as well as parameters and constraints the same af for

minimum required freight rate.

:

V, L, B, T, D

:

K= 1.107

:

200.0 < L < 245.0;

25.0 < B < 32.25

10.0

T <

14.0; 18.0

₌

D .Z 23.0

2.0 < GM < 3.5;

14.0 < V < 16.0

FB

.

<D-T+ s

nun =

:

DWR = 70000.0; SFR = 46.0

:

minimum R.F.R

Test Case No. 4

Free variables

:

CB, L, B, T, D

Parameters

:

V = 16.0

Constraints

inequality

:

200.0 < L < 245.0;

28.0 < B < 32.25;

=

13.0< T <

14.0;

18.0 < D <25.0 ;

2.0 < GM <

3.5;

0.80 < CB < 0.83

FB

D --T + s

min =

equality

:

DWR = 70000.0;

SFR = 46.0

Object function

:

minimum building costs

Test Case No. 5

Free variables, parameters and constraints are the same as for test case

no. 4.

Object function: minimum R.F.R.

For output samples of all test cases see Appendix No. 3.

9. TERMINATION CONDITION

The program terminates the calculations when either:

optimal solution is

found. In such case all output are printed out:

basic ship data, experience coefficients (basic ship data analysis),

design requirements and data of the "optimum" ship;

or:

number of search evaluation exceeds the maximal number allowed (ab. 900),

or calculation time is too big. Such case can occur when the start

point lies too far from the optimum.

Advice: after the investigation of between results change start values,

by changing of upper and lower limits for free variables;

does not exist any solution in the space delimited by equality and

in-equality constraints (requirements). It can occur when the ship has to

be designed either for given deadweight or/and stowage factor.

Advice: run the program 4 times, searching for maximal and minimal

values of deadweight and stowage factor (see sample, Fig. 7);

procedure LINORT (solution of set of linear equations) DOES NOT WORK.

The proper warning will be printed out. It can occur when search

ceeds in the range of out of validity of any of the calculation

pro-cedure (method), e.g. resistance calculation for blockcoefficient

value greater than 0.845.

Advice: investigate the between results and change limits.

10. RECOMMENDATIONS

Before starting of the program, following proceedings will be recommended:

select free variables,

- check if some of these variables can be replaced by parameters and if

yes, then do it (e.g. restricted draught and/or breadth),

-9-- 9--

10-- upper and lower limits for blockcoefficient

can't exceed the values

of 0.845 and 0.535 (0.535 < CB <0.845);

- parameter K from formula for blockcoefficient

can't exceed the value:

K.( 0.845 + 0.5 m V .

/VL

/0.3048

min

max

Where: V

.

.= lower limit for speed

min

= upper limit for length.

max

The following differences between upper and lower limits for free

variables are recommended:

: 2

knots

CB

:

0.02

:

20.00 m

:

2.00 m

T

:

1.00 m

:

2.00 m

In such way the "length" of the search space in the first cycle amounts

20 steps in each direction, so makes it quite roomy for the search.

Range of deadweight: 30000 - 80000 T.

11. JOB CARDS

The program is written in ALGOL-60 language,

T.H. Delft compiler.

The following JOB CARDS are required, when program

is executed from

the disc:

//

job statement

/ROUTE

PRINT LOCAL

/mJOBPARM

LINES=5,CARDS=0

//

EXEC

PGM=membername,REGION=256K,PARM=ISIZE=MAX-2K'

//STEBLIB

DD

DISP=SHP,DSN=libname

1/ALGLDD01

DD

SYSOUT=A

//ALGLDD02

DD

SYSOUT=A

//ALGLDD03

DD

SYSOUT=A

//SYSPRINT

DD

SYSOUT=A

//SYSIN

DD

m

"DATA"

For more explanation see T.H.BleIft

Report (Rekencentrum): "ICL-RECEPTEN

Literatiare

Auf'm Keller, W.H., "Extended diagrams for determining the resistance

and required power for single-screw ships", I.S.P. 1973.

Benford, H., "Principles of engineering economy in ship design",

Trans. SNAME, 1963.

Benford, H., "The practical application of economics to merchant ship

design", De Ingenieur, Delft, 1966.

Deetman, E., "The design of a displacementship", T.H. Report No. OvS-2,

1976.

Erichsen, S., "Optimizing containerships and their terminals",

Trans. SNAME, 1972.

Fischer, W., "Procedures in preliminary ship design" (Applied to the

Australian ore trade), RINA, 1971.

Fischer, W., "The inclusion of IMCO tanker design constraints in

general optimization procedures", Trans. SNAME, 1973.

Gallin, C., "Entwurf wirtschaftlicher Schiffe mittels

Elektronenrech-ner", Jahrbuch der STG, 1967.

Gallin, C., "Which way computer aided preliminary ship design and

op-timization", ICCAS-Papers, Tokyo, Japan, August 28-30, 1973.

Heinecke, N., "Rechnereinsatz in der Entwurfspraxis", Kontaktstudium

Institut filr Schiffbau der Universitat Hamburg, 1976.

Hagen, E., Johnson, I., Ovrebo, B., "Hull steel weights of large oil

tankers and bulkcarriers", European Shipbuilding, No. 6, 1967.

Holtrop, J., "Computer programs for the design and analysis of general

cargo ships", NSRC, TNO Report No. 1575, 1971.

Kuniyasu, T., "Application of computer to optimization of principal

dimensions of ships by parametric study", Japan Shipbuilding and Marine

Engineering, 1968.

Kupras, L.K., "Further development of better point algorithm and its

application in preliminary ship design", T.H. Delft Report, 1975.

Kupras, L.K., "Procedures in preliminary ship design", T.H. Delft

Report, 1976.

Kupras, L.K., de Zwaan, A.P., "Preliminary ship design with

interac-tive graphical aids", The Naval Architect

,

July 1977

Murphy, R., Sabat, D.J., Taylor, R.J., "Least cost ship characteristics

by computer techniques", Marine Technology, 1965.

Mandel, P., Leopold, R., "Optimization methods applied to ship

de-sign", Trans. SNAME, 1966.

Nowacki, H., Brusis, F., Swift, P.M., "Tanker preliminary design

-an optimization problem with constraints", Tr-ans. SNAME, 1970.

Nowacki, H., "Modern approach to integrated ship design", Symposium

on 'Development in Merchant Shipbuilding', Delft, 1972.

Schnell, G., "Contribution to optimization methods with penalty functions

for use in ship design", Doctor Thesis in preparation, 1976.

Schneekluth, H., "Zur Frage des Rumpfstahlgewichtes mid des]

Rumpf-stahlschwerpunkts von Handelsschiffen", HANSA No. 18, 1972.

Sading, H., Poulsen, I., "Methoden der Programmierungvon Aufgaben des

Schiffsentwurfs", Jahrbuch der STG, 1975.

Todd, F.H., "Some further experiments of single-screw merchant ship

START

BASIC SHIP DATA REQUIREMENTS TO THE DESIGNED SHIP

ANALYSIS OF BASIC SHIP DATA, CALCULATION OF "CORR.COEFF." CALCULATION OF DATA REQUIRED BY OPTIMIZATION METHOD BLOCK OF OPTIMIZATION

OPTIMIZATION PROCEDURE

OUTPUT: OPTIMUM SHIP DATA

SHIP DESIGN MODEL: (SYNTHESIS1 CREATION OF: FREE VARIABLES, CONSTRAINTS, OBJECT FUNCTION

AND CALCULATION OF THEIR VALUES

DESIGN

IS

ERVICE SPEED V

ONSIDERED AS FREE

TABLE?

KEEP CONSTANT VALUE FOR SPEED

IS

ANY OF THESE

IMENSIONS CONSIDERED AS

E VARIABLE?

KEEP CONSTANT VALUE FOR SUCH DIMENSION

IS

BLOCKCOEFFICIEN

CONSIDERED AS FREE

ARIABLE?

CREATE & CALCULATE INEQUAL. CONSTR. FOR UPPER & LOWER LIMIT OF THE SPEED V CREATE & CALCULATE INEQUAL. CONSTR. FOR UPPER & LOWER LIMITS OF THAT VARIABLE CREATE & CALCULATE INEQUAL. CONSTRAINTS FOR UPPER & LOWER LIMITS OF CB

NJ

KEEP CONSTANT VALUE FOR PARAM. K

CALCULATE

CONSTRAINTS FOR CB

0.535<GB <0.845

CALCULATIONS OF:

RESIST. & PROPULSION, KEIGHTS, STORES, SUBDIVISION, CARGO CAPACITY, DEADWEIGHT, ... ETC.

YES YES

T.H.

Fig. 2

SAMPLE OF BULKCARRIER SHIP DESIGN MODEL

DELFT 1976

END

FIG.

I

SHIP DESIGN OPTIMIZATION PROGRAM

STOW.FACTOR-SFR

CARGO CAPACITY-VOLTR7

THE SHIP WILL BE DESIGNED FOR RE- QUIRED:

CREATE & CALCULATE INEQUAL. CONSTRAINT FOR UPPER LIMIT OF DEADWEIGHT

UPPER LIMIT OF

CARGO CAPACITY

CREATE & CALCULATE EQUAL. CONSTRAINTS FOR REQUIR. DEADWEIGHT & STOWAGE FACTOR CREATE & CALCULATE EQUAL. CONSTRAINTS FOR REQUIRED DEADWEIGHT & CARGO CAPA- CITY CREATE & CALCULATE EQUAL. CONSTRAINTS FOR REQUIRED STOW. FACTOR & CARGO CAPACITY CREATE & CALCULATE EQUAL. CONSTRAINT FOR REQUIRED DEADWEIGHT & INEgUAL. CON- STR. FOR UPPER LIMIT OF CARGO CAPACITY CREATE & CALCULATE EQUAL. CONSTRAINT FOR STOW. FACTOR & INEQUAL, CONSTR. FOR UPPER LIMIT OF CARGO CAP. CREATE & CALCULATE EQUAL. CONSTRAINT FOR CARGO CAPA- CITY & INEQUAL. CONSTRAINT FOR UPPER LIMIT OF DEAD- WEIGHT CALCULATION: INITIAL STABILITY - GM

XACT

VALUE FOR INT

-IAL STABILITY IS RE- QUIRED - GMR?

CALCULATION:

CREATE & CALCULATE INEQUAL, CONSTRAINTS FOR UPPER & LOWER LIMITS OF GM STABILITY CROSS CURVES

INIMAL

RANGE OF RIGHTIN

ARM POSITIVE VALUE

RE-QUIRED?

CALCULATION: MINIMAL FREEBOARD AS RE- QUIRED BY REGULATIONS

MINIMAL

FREEBOARD REQUIRED

CREATE & CALCULATE EQUAL.

FOR INITIAL

STABILITY CREATE & CALCULATE INEQUAL. CONSTRAINT FOR MINIMAL RANGE OF POSITIVE STABI- LITY CREATE & CALCULATE EQUAL. CONSTRAINT FOR MINIMAL FREEBOARD

T.H. DELFT 1976

Fig. 3

T.H.

SAMPLE OF BULKCARRIER SHIP DESIGN MODEL

DELFT

(CONTIN.)

1976

Fig. 4

CREATION & CALCUL. OF INEQUAL.CONSTRAINT FOR LOWER LIMIT, EQUAL TO THE MIN. FREEBOARD REQUIRED CALCULATION OF BUILD. & OPER.COSTS, & R.F.R. OBJECT

FUNCIT/CW

SELECTION

CREATION & CALCULATION OF THE OBJECT FUNCTION:

MIN. . BUILDING COSTS? MIN.

R.F.R. 7.

MAX. DEADWEIGHT ? MAX. STOW.FACTOR ?

.WIN.

BUILD.COSTS/DWT ?

MIN.-ILD:COSTS/R.CAP STOW.FAC7OR ? _{OFV= DWT}

OFV=BUILD. COSTS OFV= R.F.R. OFV= 1/ DWT OFV= 1

SF

MAX.

GO HOLDS CAPAC.

OFV= 1/ C.H.CAPACITY

.

OFV= BUILD.COSTS/DWT OFV=BUILD.COSTS/C.H.CAPAC.

Fig. 5

SAMPLE OF BUIKCARRIER SHIP DESIGN MODEL (CONTIN.)

T.H.DELFT 1976

Fig. 6a

LONGITUDINAL SUBDIVISION B.L. LUM .11

Fig. 6b

DOUBLE BOTTOM, WING TANKS

LWA LUB 1.1BUHP B.L. LBLHP

Fig. 6c

HOPPER TANKS

Fig. 6a,

GENERAL SUBDIVISION OF A BULKCARRIER SHIP

T.H. DELFT

6b, 6c

SF PARAMETERS: V = 14.5 T-= 12.55 K = 1.114 CONSTRAINTS: 210..s.L4230

28B

.s 32.21

17D

18.50 113 = FE . min 1.5 <Mc 3.5 0.535 ,.c..CB..g.0.845 SFmax "--4-

....\,.7T,

. Gm ,.f.. .4,..4.:.. '

,

--___....

s-=

_{_,}c.8 ....,

,

DWTmin---,...

...____ -. "---SFmin DWT 55000 60000 65000 = 32.21 DwT -

15 -

.

Fig. 7 SAMPLE OF FEASIBLE RANGE FOR DEADWEIGHT AND T.H. DELFT

STOWAGE FACTOR 1976

40

-APPENDIX 1

: COMPUTER PROGRAM LISTING

41 42 IIA(/8/):=IIA1/8/1.2; 42 43 44 0IF'DESB1/1/1>0.50THEN, 'BEGIN, IIA(/6/):=IIA(/6/).1; 0 'BEGIN' 45 IRAMIA(/6/)/1:=I3ES13(/2/)+DESB(/3/))*0.5; 1

'COMMENT' STANDARD DECLARATIONS AND DATA FOR OPT. METHOD;

46 I1A( /8/) :=1 IA1/8/1+2; 1 IINTEGERIDSN1OSN2OSN3; 47 'END': 2 eINTEGER0IARAAY'IlAI/1:10/101Al/1:3/1; 48 'IF"DES81/7/1>0.5,THEN0 3 0ARRAYIIRA(/1:18/).GAMO(/1:61); 49 'BEGIN' 4 0REAL0A.OFV; 50 IIA(/6/):=IIA1/6/1.1; 5 51 IR4I/I1A(/6/1/1:=(DE5B1/8/)+DESBI/911)*0.5; 5 'COMMENT' 52 IIA(/8/):=IIA(/8/)+2; 5

STANDARD DECLARATIONS FOR DESIGN MODEL;

53 'ENV; 5 'ARRAY0HHH.GGG(/1:100/); 54 IIFIIDESB(/1.4)>0.51THEN" 6

,COMMENTIDECLARATIONS FOR CARGO SHIP;

55 'BEGIN' 6 56 IIA(/6/):=IIAI/6/1+1; 6 IREAL,BAS; 57 IRA(/1141/6/1/):.(DES81/11/)+DESB(/12/))*0.5; 7 'REALIIIARRAWBASB1/1:60/1; 58 IIA(/8/):=IIAI/8/)+2; 8 'REALloARRAVIDES81/1:71/1; 59 'END"; 9 IREALIIARRAY0KNI/1;6/1; 60 ,IFIDESB(/13/)>C.51THEN, 10 'REALs'ARRAY0EXCO(/1:40/); 61 'BEGIN' 11 'REAL0IARRAY'FIN(/1:100/); 62 IIA(16/):=IIA(161)+1; 12 IINTEGER11; 63 IRAMIAl/6/1/1:=IDESB(/14/)+DESB(/15/))*0.5; 13 64 IIA(/8/):=IIA(/8/)+2; 13 IlAl/1/1:=1; 65 'END"; 14 11A(121):=2; 66 eIFIDESR1/16/)>0.50THEN, 15 11A1/3/1:=3; 67 'BEGIN, 16 IIA(/4/):=1; 68 IIAI/6/1:=IIAl/6/1+1; 17 IIA(/5/);=950; 69 IRAMIA(/6/)/1:=IDESB1/1711*DESBI/18/1)*0.5; 18 IIAI/10/1:=3; 73 IIAl/8/1:=IIA(/8/)4.2; 19 A:=0.1; 71 'END,; 20 OSN1:=IIAI/1/); 72 . 21 05N2:=IIA(12/); 72 22 DSN3:=IIA(/3/); 72 "IF0DESBC/19/))1.58DESBI/20/1>1.51THEN, 23 SETTING(OSN1,132,60); 73 'BEGIN, 24 ,IFI-IDSN2=DSNIITHEN1 74 IIA(/7/):=IIAl/7/)+1; 25 SETTING(DSN2,132,60); 75 GAMOI/IIA(/71)/1:=25.0; 26 '1F'-.0SN3=DSN21...DSN3=DSNIITHENI 76 IIAC1711:=IIA(17/).1; 27 SETTING(DSN3,132,60); 77 GAMOUIIA(/7/)/1:=C.D3; 28 78 'END': 28 INREAL((j1BAS); 79 11F0DESB(/191)>1.5AOESB(/21/1>1.50THEN, 29 IF0BAS>0.50THEN'INARRAYIO.BASB); 80 'BEGIN, 31 INARRAY(O,DES8); 81 IIA(17/):=IIA(/71)+1; 32 82 GAMOC/IIAI/7/1/1:=25.0; 32 IIAI/6/1:=0; 83 IIA117/1:=IIAI/7/).1; 33 IIAl/7/1:=0; 84 GAMOMIA(17/)/1:=30.0; 34 IIA1/8/1:=0; 85 'END'; 35 I1A(/9/):=0; 86 '1F1DES81/2C/I>1.58DESB(/21/),1.50THEN' 36 36 'IFIDESS(/4/)>0.50THEN' 87 88 "BEGIN, IIA(/7/):=11A(/7/)+1; 37 'BEGIN' 89 GAMOUI1A(/7/)/1:=0.03; 38 IIA1/6/1:=IIAI/6/1+1; 90 IIA1/7/1:=IIA(/7/)4.1; 39 IRAl/11A1/6/11):=IDES81/5/I+DES81/6/1)*0.5; 91 GAMOUTIAI/7/1/1:=30.); 40 'END': 92 'ENO':

93

'IF'DE581/19/1>1.5ADES8(/23/))1.51THENI

150

I1(A(/2*11A(/6/)+I/):=0.20:

94

'BEGIN'

151

'END':

95 96

IIA(/7/):=IIAI/7/1+1:

GAMOUIIA(/7/)/1:=25.0:

152 153 01FIDESR4/1C/I>C.56THEN,

'BEGIN'

97

IIA(/8/):=IIAI/8/1+1;

154

1:=I+1;

98

'END";

155

IRAI/IIA1/6/1+1/1:=0.10:

99

01F0DESB(120/)>1.5&DESB(/23/)>1.5,THENI

156

IRAI/2*11A(/6/)+1/):=3.02;

100

'BEGIN'

157

'END";

101

IIA(/7/):=IIA(/7/)+1:

158

'IF'DES81/13/1>0.5'THEN'

102

GAMOI/IIA1/7/1/1:=0.03;

159

'BEGIN'

113

IIA(/8/):=IIA(/8/)+1;

160 104

'END':

161

IRWIIA(/6/1+I/):=0.051

105

01F'DES3(/21/)>1.58.DESB(/22/)>1.50THEN°

162

!=1ff*IIA(/6/)+1/1:=0.01:

106

'BEGIN'

163 107

IIA1/7/1:=IIA(/7/)+1;

164

'IF'DES8I/16/1)0.51THEN0

108

GAMO(/IIA1/7/1/):=30.01

165 'BEGIN' 109

II4I/8/1:=IIA(/8/)+1;

166 110

'END':

167

IRAI/IIA1/6/1+I//:=0.10:

111

0IF0DES8(/22/1>1.51iDESBI/23/1)1.51THENI

168

IRA1/2*11Al/6/1+I/):=0.02;

112

'BEGIN'

169

'END':

113

IIA(/8/):=IIA(/8/)+1;

170 114

IIA(/8/):=IIA(181)+1;

170 115

'END";

170 'BEGIN' 116 119 120

01FIDES6(/25/)>-C.51THEN'IIA(/8/):=IIA(/8/)+21ELSE.

'BEGIN'

IIA(/7/):=IIA(/71)+1:

171 172 172

'ARRAYI3RA(/I:IIA(/6/)/),GST(/1:IIA(/9/)/), Hl/1;'IF'IIA(/7/)=0,THENI1'ELSEIIIA1/7/1/),

G1/1:0IEIIIAl/8/1=0,THEN010ELSE'IIA(/8/)//:

121

GAM01/IIA1/7///):=0.C1:

172

'ARRAYIX(/1:1)0/);

122 "END": 173

0ARRAY.ACCUR(/0:IIA(/9/)/);

123

'IFICIES81/26/1>0.01THEN.IIAl/8/1:=IIAl/8/1+1;

174 'BEGIN' 125

*IF'DES51/2//,(-0.51THEN

175 126

'BEGIN'

175 127

II41/7/):=IIAI/7/1+1;

175 128

GAM0(IIA(/7/)/):=0.01:

175 129 IEND00ELSE' 175 131

IIA1/8/):=IIA(/8/)+1;

175 0PROCEDURE'FRE8601L,D,T,C804,FREEBI; 132

IIA1/9/1:=IIA(/7/1+IIAl/8/);

176

'REALILO/T.CB,W,FREEB;

133

I:=0;

177 'CODE'; 134 178 134

IIFIDESB(/4/)>0.5,THEN0

178 135

'BEGIN'

178 'PROCEOUREIPOWKEL(L,3.79C80,SERCONIONETAMOB): 136

1:=I+1:

179 IREAL01.03,T,CB,N,SERCON,VIETAM,PBOCODE'; 137

IRAMIA(/6/)+I/1:=0.005;

181 'PR3CEDURE140LBUL(L,8,T,D,CBIWA,WF,CAM,HDB,LMA,LFP,INI,UW890MAIIIN2rL1011 138 139

IRAl/2*11A1/6/1+I/):=0.001:

'END':

182 182 LWA,LUL.IN3oLUB,LUD,NOUHP,UBUNP,LBUHP,OUHP,NOLHP,UBLNPILBLHPOLNP, VOLHol'OLLU,VOLTIIVOLDB,VOLUVIT,VOLLVIT,VOLUHP,VOLLHP); 140 140

11F'DESB(/1/)>0.51THENI

182 183 'REAL'L.811790.C4.WA.WF,CAM,H081LMAI,LFP,UWB,UNA,LWB,LWAgLUL.LUBILUO, VOLN,VULLUOIOLT,VOLDB,VOLUNT,VOLLWT; 141

'BEGIN'

183 'REAL'UOUHP,LBUHPOUHP,UBLHP.LBLNP,DLHP.VOLUMP.VOLLHP: 142

I:=I+1:

184

'INITEGERIINIgIN2,IN3:

143

IRA(/IIA1/6/)+I/I:=0.10;

185 'INTtG-ia0NOUHP,NOLHP; 144

IRA(/2*I1A1/6/1+1/):=0.02;

106 'COOE'; 145

'END':

187 I 146 147

0IFIDESBI/7/1>0.50THEN'

'BEGIN'

187 168 'PRUCEDUREIFIEIBUL(LO,TO.CB,PB,S7EilltSTER2,LO.SIG,WSTE.WOWNWENGINSM, _{KGSTE,KGOUTAGENG,KGSM);} 148

I:=I+1;

188 'aEAL°1.0317.D.C110b.STER19STE42.1.1(F,5IGIWSTE,WOUT,WENGOISM.KGSTE,KGOUT, -4 149

IRAI/IIA(/6/)+1/1:=1.00;

109 KGENG.KGSM;

189 190 190 191

'CODE';

'PROCEDURE'STOBULIL,B.T.D.C.HDB.PB.RAD,PORTIM.V.DOCTIM.WHOFVHOINDO, VD0114LO.VLO.WFW,WSTOR.KGST0R.TRIPTI.YWHOIYWDO,YRLOU

222 223 224 225

1(01/7/):=450000; 1(01/8/) :=o;

IFORII:=11STEPU'UNTIL17.00s

1(01/8/):=M01/8/)+KOUIM

191 REAL61.80.90.CB.HDB.PB.RAD.PORTIM,VIDOCTIM.WHOINH00100.V0094LO,VLO. 226

K01/9/):=245*(IRSTE*1100/88))**0.835)*25;

192 WFWIWSTOR.KGSTOR,TRIPTI.YWHO.YWOOtYWL0; 227

KOl/10/):=383*INSTE*(1C0/88))**0.56)*26;

192 'CODE*; 228

K01/11/1:=1(0(/8/)+KO(/9/)4.K0(/10/);

193 229

K0(/12/):=3000004.3#DW+36*PB;

193 230

1(0(/13/):=4500004.280PPB;

193 PROCEDURE'RGBUL 231

(01/14/):=1611000+105*Pd4.23*Dw;

194 IL.B.T.0.CB.HOB9IN1.UWB,UWA.IN2gLWB.LWA.NOUHP.DUHP.NOLHP, 232

K01/15/1:=RO(/13/)+KO(/14/);

194 DIHNIUD,WOLH.VOLLU,VOLTINOLDB,103LUWTOFOLLWT,YOLUHP.VOLLHP.WSM,KGSM, 233

1(01/16/1:=(5*PB+1*DW)*26;

194 WSTORIAGSTOW.WCAR.KGCAROWOW); 234

KJI/17/1:=(1.5*Pd4.0.28*D4)*26;

194

'REAL'

235

K01/18//1-403(/11/)+KO(/12/).K0(/15/)+KO(/16/)+KO(/17/);

194 1.0.7.0.CB.HDB.U43.UWA.LWB.LWAIIDUHP, 236

'COMMENT.BUILD.COSTS FOR 1973 WERE CALCULATED;

195

OLHP,LUD.VOLH,VOLLU.VOLT.YOLDB.VOLUWT.VOLLWT.VOLUHP.VOLLHP,WSM.KGSM, 236

1COMMENTiNOW 4E INTRODUCE THE CORR. FOR INFLATION;

195

4ST0RIAGSTORIWCAR.KGCAROW,Cil

236

1COMMENTiTAKE MEAN YEARLY INFL.FACTOR WHICH WAS ESTABLISHED;

195

l'INTEGEMIN1.1N2,NOUHP.NOLHP;

236

'COMMENT' FOR PERIOD 1973NOW;

196

'CODE';

236

'FOR°I:=OSTEPl'UNTIOYEAR-1973001

197 237

1(01/18/1:=K31/18/)*INFL;

197 238

8101:=K0l/18/;;

197 9PROCEDURE'PANKUP(LOIT.D.C3.WAIWF,RN;; 239 IENDCOSBUL; 198 'REAL1L.B.T,D,CB,WA.WF; 240 199

'REAL"ARRAY'RN;

240 1PROCEDURE'RFRBUL(L.BoTICB.V.PB.DWIWSTUR.YWHOgY400.YWLOgRAD.PORTIM. 200

'CODE';

241 DOCTIMOCAL.NCANgbALTRI,YIR,YDEP.BKB.YRTRItYCAL,YPTIM.ANCAROCOSTS.Y); 201 241 1REAL'IL,B.T.CB.YIPBOW.WSTORgYWHO,YWDO.Y4LO.RAD,PORTIMIDOCT1M,NCALOCAN 201 242 BALTRI.YIROIDEPOKB,YRTRI.YLALtYPTIM.ANCAR.UCOSTS; 201 'PROCEDURE'BMKB(B.T.C81K8.13M); 242 °REAL"ARRAY6Y; 202

iREAL°897.03.1(8.BM;

243 'CODE'; 203 IICODEs; 244 204 244 'COMMENT' DLSOPT; 204 244 204 IPROCEDUREUNSTAB(B.T.CB.KBOMIKG.GMI; 244 'PROCEDURE'OESOPT(HIGOFV,X); 205 REAL6B.T.CEI.KBOM.KG.GM; 245 0ARRAY8H,G,X; 206

'CODE':

246 IREALIOFV1 207 247

'BEGIN'

207 248

'INTEGER'!;

207 sPROCEDUROAAMSTA13.T.C8,1(G,KN.RANSTA.RARM;; 249 DES3021EXCO,DES6.FIN.HHH.GGG.X,FREB6000W1(EL,WEIBUL,VOL8UL, 208 'REAL.613.7,CB.KG,RANSTAIRARM; 250 STOBUL.KGBULORKLI.INSTAdtPANKUP.ARMSTA,COSBULAFRBULU 209

REALitARRAY'RN;

250

'IWIIIA(/7/)=OITHEN"GOTO'DES1;

210

'CODE';

252

'FOR*1:=1STEPil'UNIIL'IIAI/7WOO.

211 253

HUI/I:=HHHUI/l;

211 254 DES1: 211 254

sIF6IIAI/8/1=01THEN"GOTO'DES2;

211

'PROCEDURE' OSBULIL.B.D.D408,WSTEILUL,LUB.YEAR.INFLOKB;;

257

'FORi1:=1'STEPtPUNTILIII4(/8/)100'

212 gREAOL.B.D.DW.PB.WSTEILULILUB,YEAR.INFLOKB; 258

GUI/J:=GGG(/I/);

213

'BEGIN'

2,9

DES2: 214

REAL'IARRAWKOl/1:18/;;

259

9IF'DES6(/71/)(1.58THEN8OFV:=FIN(/78/)/1000000;

215

'INTEGER' I;

262

gIFI3ES3I/71/)(2.530ESB1/71/)>1.511THEWOFV:=FIN(/80/)1

216

1(0(11/):=5550*(3W**0.546);

264

01FIDES13(/71/)(3.530ESB1/7111>2.51THENIOFV:=100000/FIN(/7/):

217

1(01/2/):=4S1E*(100/88)*1.35*650;

2o6

'IF'0ESB(/71/)<4.5S0ESBI/71/1>3.51THEN'OFY:=10/FINI/59/1;

218

KO(/3/1:=4STE*(100188)*604AL*13*B+4*0)+10*D*8)*40;

263

91FIDESM/71/)(5.5&DES131/71/1>4.51THEN'OFV:=100000/FIN1/8/);

219

1(01/4/1:=LUL*(LUB**2)4.50;

270

elFIDES3l/71/1<6.53DES3(/71/)>5.51THENe

220

KO(/5/):=7*L*84,041400*45*SORTIL;;

271

OFif:=J.0C1*FIl1/70/1/FIN1/7/11

221 K01/6/1:=25O..)30+1500:1*45; 272

IFDES8(/71/)<7.5,30ES8(171/)>6.5'THcNI

273

OFV:=0.001*FINI/78/1/FINI/8/1;

303

LKF:=mASB1/54/1;SIG:=BASB1/55/1;

274 275 275 ₂₇₅

'IFIDES131/71/)(8.51LOES8f/71/7.51THEN.

OFV:=FIN(/59/);

0IF'DESIII/71/1<9.54DE56U71/16.8.5,INEN0

OFV:=FINE/7/1/100000.0;

305

3)6

307 307

LKF:=LKF*L;

WEIBUL(L,B,T,),CF,P6,STERI,STER2,LKF,SIG,

WSTETWOUT,WENG,WSM,KGSTE,KGOUI,KGENG,KGSM);

LKF:=LKF/L;

275 275

'END':

'COMMENT' END OF DESOPT;

308 310

'IFI8ASBI/9/1>10.01THEN'EXCO(/2/):=BA5131/9/1/W5TE

IELSE'EXCOI/2/1:=1.00;

276

'COMMENT' PROCEDURE COML.;

312

0IF'8AS61/1,1/160.11THEN'EXCOl/3/):=BAS8(/10/1/KGSTE

277 0PROCEDURE'COEBUL(BASB,EXCO); 314

'ELSE'EXC0(/3/):=1.00;

278 'ARRAYI8AS8,EXCO; 316

IIFIBASIA/12/)>10.01THEN'EXCOl/5/1:=BASB(/12/)/WOUT

278 'BEGIN' 318

'ELSE'EXC04/5/1:=1.00;

278 279 'REAL' 320 322

IIFIBASB1/13/1)C.I'INEN'EXCOE/6/1:=BAS8(/13/)/KGOUT

'ELSE'EXCOI16/1:=1.00;

279 279 279 279 279 279 279 279 279 279 279 279 279 279 279 279 279 279 279 279 279 279 279 279

1.03,70101,C8,

DWOWROWMA, VOLTR,VOLTMA,SF,SFR, WSTE,KGSTE,XGSTE. WOUTOGOUTOCGOUTt WENG,KGENG,XGENG, WSMIKGSM,XGSM, WSTOR,KGSTORIKGSTOR, WCAR,KGCAROMCAR, KG,XG,GMIKBAM,KM, 1.108,WAIWF,W,CAMOWBOWA,LWB,LWA, LULILUB,LUD,LMA,LFP, USUMP,LESUNP,DUMP, UBLHPILBLHPOLHP, VOLH/VOLLUOIOLTINOLOB, VOLUNTIFVOLLWTIVOLUHP,VOLLHP, STER1ySTER2p1AFtSIG. NtSERCONgETAM,ETABUL,PB, FREE3, RAD,PORTIM,DOCTIM, WHOO/HO,WOO,VDO,WLO,VLOTWFWI, TRIPTI,YWHO,YWDO,YWLOACALOCANOALTRI. YIR,YDEPOKB,YRTU,YCALOPTIM,ANCAR, UCOSTS,

324 326 328 330 332 333 334 335 336 337 338 339 340 341 342 343 344 345 345 346 348 350 351

IIFIBASB1/15/1610.00THEN'EXCO(/8/):=8ASa(/15/)/WENG

'ELSE'EXCO(/8/):=1.00;

IF0BASW16/1>C.IITHEN'EXCO(/9/1:=BASBI/16/)/KGENG

'ELSE'EXC0C/9/1:=1.'50;

WSTEW5TE*EXC0(/2/);

KGSTE:=KGSTE*EXCOI/3/1; WOUT:=WOUT*EXC0(/5/); KGOUT:=KGOUt*EXCN/6/1; WENG:=WENG*EXC01/8/1; KGENG:=KGENG4EXCOI/9/1;

0W:=BAS8(/7/);

WSM:=WSTE+WOUT+WENG;

EXCO(/11/):=(1.03*L*i*T*C8DW)/WSM;

WSM:=WSM*EXC01/11/1; WSTE:=MISTE*EXC01/11/1; WENG:=WENG*EXCOI/11/1; W3UT:=WOUT*EXC01/11/1; 'COMMENT' W6M IS CALCULATED MORE PRECISELY 3UT WSTE 4ENG WOUT APPROX.; KGSM:=IWSTE*KGSTE.WOUT*KGUUT+WENG*KGENGI/WSM;

01F,SASBI/18/)>3.11THEN'EXCOI/12/):=BAS81/18/1/KGSM

0ELSE'EXCOI/12/I:=1.00;

KGSM:=KGSM*EXCOl/12/1;

LMA:=0.:67*L+3.U)136*P3;

280 RANSTA,RARM; 352

EAC01/14/1:=6ASul/41/1/L4A;

280 0REAL'YEAR,INFL; 353

LMA:=EXC0(/14/)*LMA;

281 'INTEGER' 354

WA:=8AS8l/31/);wF:=BAS31/32/);CAM:=3ASBI/34/1;

281 INI,IN2,I143,NOUHP,NDLHP; 357

HD8:=BAS8l/30/1;

282 358

'IF0BAS9(/35/)<0.11THEN'INI:=09ELSE'INI:=2;

283

'IF0BAS<0.5"THEN,

362

UW8:=RASi1/35/1;JWA:=RAS1(/36/);

284

'BEGIN'

364

0IFIBAS8(/37/)<C.I.T1EN'IN2:=00EL5E'IN2:=2;

285

'FOR'I:=1'STEPII'UNTIL143°000

368

LW8:=8AS8(/37/);LWA:=BAS3(/38/);LFP:=BAS8(/43/);

286 EXCOM/I:=1.00;

371

LUL;=8AS61/39/1;LUL:=LUL/(LLMALFP*L);

287 IGOTO0C01; 373

IN3:=2;LUB:=8AS31/4C/MUD:=8AS(/41/);

287

'END':

376

NOUHP:=8ASW44/1;(J8JHP:=8AS81/45/);

287

LIIII8AS8(/1/1:8:=1ASR1/2/1;T:=RASR(/3/);

378

LBUHP:=8ASB(/46/);DUHP:=BASSI/47/1;

290

CB:=BAS8C/6/1;N:=8AS8I/56/1;SERCON:=8ASS(/57/);

380

NOLHP:=8AS8(/48/);UBLHP:=BAS3(/49/);

293

V:=BAS8I/5/1;ETAM:=BAS8(/58/);ET48UL:=BASB1/59/1;

382 LEIL1IP:=BASI3(/5)/1;OLHP:=BAS81/51/1; 296 POWKEL(L,B,T,CB,N,SERCON,V,ETAM,PB); 384 297 P81=P8/ETABUL; 384 VA3UL(L13,10,C3,WA,WF,:AM,HDBpLMA,LFPIIN1,UWBOWAgIN20.WB,LWAt 298

EXCO(/1/):=8ASBI/60/1/Pe; P8:=EXCO(/1/)*PB;

385 LUL,IN3,LUB,LJD,q0UHP,UdUHP,LBUHPODUNPIN3LHP,U8LMP,L8LHP,0LHP, 300

0:=8ASBI/4/1;STERI:=8AS81/52/1;STER2:=8AS8(/53/);

385 VOLN,dOLLU,VOLT,V3L013,VJLUWT,VOLLWT,VOLUHP,VOLLHP): 385

385

EXC04/15/1:=BASBI/8/1/VOLTI

428 DNOWR.DWMA, 386

JOLT:=EXC9(/15/)*VOLT;

428 VOLTRIVULTMA,SFISFR, 387

EXCOl/16/1:=1.00:

428 WSTE,KGSTEOGSTE, 388

EXC01/17//:=1.00;

428 MOUTIKGDUTIXGOUT, 389

NST0R:=BASB(/20/):

428 WENG,KGENGOMENG, 390 391 391 392 392 392 392 394

KGST0RBASIA/21/):

KG8UL(L.R.T.DICRO3R,INI.UWBOWAFIN2ILMB,LWArNOUHP,DUHP,NOLHP. DOPpLUD,VOL4.VOLLU,VOLT,VOLD8,VOLUMTIVOLLMT,VOLUMPIVOLLHP. WSM,KGSM,WSTOR.KGSTUR,WCAR,KGCAR,DW,KG):

IIF'dASI5I/23/1>C.1'THEN'EKCO(/19/)3ASB(/23/)/WCAR

'=LSE'EXC0(/19/):=1.00:

428 428 428 428 428 428 428 428 428

WSM.KGSm,XGSM, WSTOR,KGSTOR,XGSTOR, WCAR,KGLAR,XGCAR, KG0(G.GM,K8,84,KM, 1-108,WAIMF,WitCAMOWB,UWA,LWB,LWA, LUL,LUB.LUDgLMAILFP, UBLIFIPILBUNP,DIMP, UBLMP,LBLIO.DLNPJ VOLN,VOLLU,VOLT,VOLOB,

396

'IFIBASIA/24/)>C.18T4EN'EXC0(/20/):=BASdI/24/1/KGCAR

428 VOLUNI,VOLLWT,VOLUNP,VOLLNP, 398

'ELSE'EXC01/20/1:=1.00;

428 STER1gSTER2ILKF,SIG, 400 WCAR:=WCAR*s_XC0(/19/): 428 N9SERCON'ETAM'ETABUL,P8, 431 KGCAR:=KGCA4*EXC0(/2U): 428 FREEB, 402

'IF'dAS6(/26/)>0.1"THEN°

428 RADOORTIM,30CTIM, 403 EXC0I/22/1:=8ASB(/26/)/(IWCAR*KGCAR+WSTOR*KGSTOR+WSM*KGSM) 428 WN011040,WDOIVD0OILD'VLO,WFW, 404

/11.03*L*134.14C8IIIELSE'EXCOl/22/1:=1.00:

428 TRIPTItYWHO,YWDO,YWLO,NCAL,NCAN,BALTRI, 406 KG:=EXCU(/22/)*(WCAR*KGCAR+WSTOR*KGSTOR+WSM*KGSM) 428 YIR.YDEPOKB,YRTRI,YCAL,YPTIM,ANCAR, 407 407

/(1.03*L*B*T*CB):

428 428 OCOSTS, RANSTAtRARM; 407 13MKB(BI,T,C13,K3,5M); 428 IREALIYEARtINFL: 408 429

'ARRAI'lf(/0:30/);

408

IIF*BAS8l/29/)>C.10THEN'

430 409

EXC0(/251):=3AS8(/29/)/(KB+BM)6ELSE'

430

'INTEGER'IFJ.KK:

411

M01125/1:4.0C:

431

'INTEGER'INI,IN2,IN3;

412

K8:=EXC0(/25/)*KB:

432 'INTEGER'NOUNP,NOLNP; 413

3M:=EXCO(/25/)*BMI

433 414 433 'COMMENT'START; 414

INSTAB(8,TICB,K8,8M,KGIGM):

433

I:=0: J:=0; KK:=0;

415 415

'IF'13458(/28/)>G.G0THEN0

436 436 416 418 419 420 420 422

EXCOl/2411:=BA581/28/1/GMIELSEI

EXC0(/24/):=1.00:

GM:=EXCO(/24/)*GM: COI:

'END':

ICOMMENTIEND OF COEBOL:

436 437 438 440 441 442 443

'IF'DE58I/4/1>3.5°THEN'

"BEGIN'

KK:=KK+1;J:=J+I:

CB:=X(/KK/):

GGGI/J/II=DES8(/5/)C13: J:=J+1; GGG(/J/):=Ce-3ESB1/6/1;

422 444

'END':

422 445 422 422

'COMMENT' Pi0C. DESCO2;

446 447 01FIDES81/1/I>0.5"THEN, 'BEGIN' 422 PROCEDURE'DES002(EXCD,DES8,FINgHHH,GGG,X,FREB6000MKEL,WEIBULgVOLBUL, 448 KK:=KK+1;J:=J4.1; 423 STOBUL,KGBULOMKB,INSTABOANKUP,ARMSTA,COSBUL,RFRBUL): 450

V:=X(/KK/);

423 'REAL01ARRAY'EXCOOESB,FINI 452

GGG(/J/):=DESB(/2/1V:

424 'ARRAYIHHIgeGGG,X; 453

J:=J+1:

425 'PROCEDURE'

453 GGGI/J/I:=VDESB(/3/);

426 FREB609POWKEL,WEIBUL,V0LdULISTOEWLOWBUL.BMKB,INSTA8,PANKOP, 454 'END' 426 ARMSTA,COSBULAFRWL: A55

'ELSE'V:=DESBI/2/I;

426

'BEGIN'

457

0IF'DESB(/7/)>0.50THEN'

427 458 'BEGIN' 427

'REAL'

427

L'13010,V,C4

521 LKF:=LKF/1.; 459 Kitt.KK+1;JJ+1; 522 WSTE:=WSTE*EXCO(/2/): 460 L:=XI/KK/); 523 KGSTE:=KGsTE*EXC3(/3/); 461 GGG(/J/):=0ES9(18/)L; 524 4OUT:=WOUT*EXCOI/5/1; 462 J:=J+I; 525 KGOUT:=KGOUT*EXCO(/6/); 464 465 GGG(/J/):=L-0ESB(/91); 'END' 526 527 WENG:=WENG*EXC01/8/1; KGENG:=KGENG*EXC:A/9/); 4(36

'ELSEIL:=DEsmism

528 WSM:=WSTE+WOUT+WENG: 468 459 470 471 472 473 475 476 477

,IF,DESB10/)>0.50T1ENI 'BEGIN' KK:=Ki+I;J:=J+1; B:=Xl/KK/); GGGI/J/1:41ES3(/11/)-3; SGGI/J/1:=8-0ESB(/12/); 'END' IELSEIB:=DES41/11/1;

529 530 531 532 533 533 533 533

WSM:=WSM*EXC01/11/); KGSM:=(NSTE*KGSTE+WOUT*KGOUDWENG*KGENG)/WSM; KGSM:=KGSM*EXC01/12/1; LMA:=DESES(/43/1; 11F'LMA<1.5.7HEW

479 IIFIDESB(/13/150.59THEN0 534 'BEGIN' 490 'BEGIV 535 LMA:=0.057*L+0.00136*PB: 481 KIC:=10(+1;J:=J+1; 536 LMA:=LMA*EXC0(/14/); 482 7:=X(/K/); 537 'END'; 493 GGGI/J/I:=DES3l/14/1-4; 538 WA:=DES3(/29/1;s4F:=DES81/30/MAM:=DESBI/32/1;

484 JJ+1;

541 H08:=5ES8(/26/); 446 GGC(IM):=TOESS(/15/1; 542 ,IFIH3B>-1.5,111136<0.C,THEN, 487 'ENS' 543 HOB:=0.C31*(600.9*B*(T**0.5)); 488 IELSPIT:=JES3(/14/); 544 01F,H08<-1.5,THLN8 490 'IF'DES8l/16/1>C.5,74ENI 545 HDB:=3.0014,1600+9*8*(7**0.5))*(1+1L-150)/300); 491 'BEGIN' 546 LFP:=DES8(/44/); 492 KK:=KK+1;J:=J+1; 547 IlF00ES8(/33/1(0.5*THEN'IN1:=3; 493 0:=Xf/KK/I: 549 0IFIDES8l/33/))6.58DES8(/33/)<1.50THEWIN1:=1; 494 GGG(/J/):=DESB(/17/)-0; 551 ,IFIDES81/33/1>1.50THENIIN1:=2; 495 J:=J+1; 553 UW3:=5ES3(/34/); UWA:=DESB(/35/); 497 GGGC/J/1:=D-0ES8l/18/1; 555 eIF°0ESEI(/36/)<.5°THEN'IN2:=0; 498 'END' 557 IIF,DES8(/36/)>k,.58,DES8(/36/)<1.52THEN'IN2:=1; 499 IELSE1D:=5ES41/17/1; 559 IIFI0ES31/3b/i>1.51THEVIN2:=2; 501 5)2 534 IIFIDE5IA/4/1<:,5,THENe 'BEGIN' CR:=DES8l/5/1-0.5*V/S0RT(L/0.3'4.3); 561 563

LWB:=DES8(/37/); LWA:=DES13(/38/); LUL:=DES8(/39/);

564

,IFIDES61/40/1<1.507HEN'IN3:=10ELSE'IN3:=2;

506 507 508

GGGl/J/I:=0.535C8: J:z.14.1; G63l/J/I:=CBC.645;

568 573

LUB:=DESB(/41/); LUD:=DESB(/42/); NOUHP:=DESB(/45/); UBUHP:=DESEil/46/1;

59

'ENO;

572

LBUHP:=DES151/47/1; DUHP:=DESB(/48/);

574

NOLHP:=DESBI/49/1; UBLHP:=DESIll/55/1; LeLHP:=DES5I/51/1;

577

DLHP:=DES81/52/I;

510 511 512 513 514 515

N:=DESS(/57/); SERC3I:=DES8t/58/1; ETAM:=DES9(/59/); ETABUL:=DES5(/60/); PUWKEL(LO,T,C30,SERCONIVtETAM,P8); P8:=PB/ETABUL: PB:=EXCO(/1/)*P8; STER2:=DES8(/54/);

578 578 579 579 579 579

V0LBULIL,d17.10,C6,4A,WF,CAM,HDB,LMA,LFP,IN101.18,UWA,I42,LMI,LWAT LULON3,LUB,LUD,NOUHPOWHPILBUMPIDUHPOOLHPOBLHP,L8LHP,DLHP, VOLH,VOLLU,VOLT,VOLD3IVOLUNT,VOLLWTIVOLUHP,VOLLHP); VOLT:=E4C0(/15/)*VOLT;

516 STER1:=0ESE11/53/1; 530 RAD:=GE5B4/61/1; P0WIIME:=DESB(/62/); 517 LKF:=DES8(/55/); 582 DUCTIME:=DES5(/63/); 518 LKF:=LKF*1.; 593 519 SIG:=DE513(/56/); 583 5104UL(LtagT,3,C30-104,P4gRA)IP0RTIM,VO0CIIM.WH0,VH00100,VDO,WLO, 520 WEI8UL(1.0.70,090A,STER1ISTER2,LKF.SIG,WSTEIW0U1.WENG, 584 VLO,WFW,WSTDA,KGSTOR,TRIPTI,YWHDIYWDO.YWL0/: 521 WS4,KGSTE,KGOUT,KGENG.KGS4); 584 584 KG8UL(Lv6,T,0,Ci,HD3,IN1gUw3,UWA,IN2gL14130.WA.NOUHP,DUHP,NOLHP,

585 DLHP,LUDeVOLHIVOLLU'VOLTINOLDB,VOLUNT,VOLLWT,VOLUHP,VOLLHP, 640 J:=J+1; 585 4SM.KGSM,W5TOR,KGSTO,WCAR.KGCARIOW,KG); 641 GGGI/J/I:=VOLTVOLTMA; 585 642 'END': 585 WCAR:=EXC01,19/1*WCAR; 643 BMKBIB,T,CB,KBOMI; 586 KGCAR:=E4C04/20/)*KGCAR; 644 K8:=K8*EXCO(/25/); 587 KG:=EXCO(/22/)*(WCAR*KGCAR*WSTOR*KGSTOR+WSM*KG5M)/(1.03*L*8*T*C8); 645 BM:=BM*EXC0(/25/); 588

DWR:=DESBI/19/l; SFR:=DES8I/20/1; VOLTR:=DESB(/21/);

646 KM:=KB+8M; 591 D4MA:=DES8(/22/); VULTMA:=DES81/23/1; 647 INSTABIB,T.CBIKB,BM,KG,GM); 593 5F:=VOLT*35.314/WCAR; 648 GM:=GM*EXCO(/24/); 594 'IF'DMR)1.5A5FR>1.50THEN' 649 111FIDES81/25/1>-0.5,THEN. 595 'BEGIN' 650 'BEGIN' 596 I:=I+1; 651 J:=J+1; 597 HHMI/1/1:=DwROW; 652 GGG(/J/):=GMDES3(/25/); 598 I:=I+1; 653 J:=J+1; 599

MHMUI/I:=SFRSF;

654 GGG(/J/1:=DFSB(/241)GM; 600 'END': 655 'END'IEL5E0 601 01F0DWR)1.511V0LTA>1.5,THEN' 657 'BEGIN' 602 'BEGIN' 658 I:=I+11 603 I:=I+1; 659 HHHI/I/):=GM-0E58(/24/); 604 HHH1/1/):=DOR-13W; 660 'END'; 6a5 I:=1.1q; 661 PANKUP(L097.01C804A9WF,KN); 606 MHHI/I/I:=VJLTRVOLT; 662 RANSTA:=DES81/26/1; 607 'END': 663 "IFIRANSTA)0.)1THENI 608 IIF'SFR>1.5.1VOLTR51.51THEN. 664 'BEGIN' 609 'BEGIN' 665 ARMSTA(8,T,C30(G,KN,RANSTA,RAAM); 610 I:=I*1; 666 J:=J+1; 611 1-11.011/I/J:=SFRSF; 667 GGGI/J/I:=RARM; 612 1:=I+1; 668 'END': 613 HHHI/I/I:=VOLTRVOLT; 669 W:=DESB(/311); 614 'END': 670 FREB60(1.,0,T,CB,W,FREE9); 615 0IF'DWR>1.56MOLTmA>1.50THEN, 671 0IFI0ES8(/27/)<;.50THEN' 616 'BEGIN' 672 'BEGIN' 617 I:=I+1; 673 I:=1.0.1; 618

HHMI/I//:=DADM;

674 HHHI/I/I:=0TFREE8; 619 J:=.141; 675 'ENE:119E1SE' 620 GGGI/J/I:=VOLTVOLTMA; 677 'BEGIN' 621

'ENV;

678 J:=J+1; 622 IIFISFR>1.5aVOLTMA)1.5*THEN, 679 GGG(/J/):=FREEB()I); 623 "BEGIN' 680 'END': 624 I:=I+1; 681 YEAR:=DES8l/69/1; INFL:=DESB(/70/); 625 HHHI/I/I:=5FRSF; 683 COSBUL(L,B,DOW,P6,W5TE.LUL,LUBIYEAR.INFLOKB); 626 J:=J+1; 684 RAD:=0E5B(/61/);110RTIM:=DES8(/62/); 627 GGGI/J//:=VOLTVOLTMA; 666 00CTIM:=DESE(/63/);NCAL:=DESBI/64/1; 628 'END'; 688 NCAN:=DESBI/65/1;BALTRI:=DES8(/66/); 629 gIF'V3LTR,I.51WWM4>1.50THENI 690 YIR:=DESB(/67/);YDEP:=DESBI/68/1; 630 'BEGIN' 692

RFRBUL( LIB, ItCB irVirP15,13W9WSTORgYWHO,YWDOpYWI-OtRAD,PORTIM gDOCTIMOCAL,

631 1:=I4.1; 693 NCANOALT(lgYIRtYDEPOKEItYRTRI,YCALtYPTIM,ANCARgUCOSTS,Y); 632 HMMI/1/1:=V0LTRVOLT; 693

'COMMENT' OUTPUT PREPARATION

633 J:=J+1; 693 FIN(/1/):=L; 634 GGGI/J/1:=DW-0WMA; 694 F1N(/2/):=B; 635 'END"; 695 FIN1/3/1:=7; 636 0IFIDWMA>1.58NOLTMA>1.51THENI 696 FIN(/4/):=D; 637 'BEGIN' 697 FIN1/5/1;=V; 638 J:=J+1; 698 FIN(/6/):=CB; 639 GGGI/J/):=DWDWMA; 699 FIN(17/):=DW;

700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715

FIN(/8/):=VOLT; FINI/9/I:=i4STE; FINC/10/1:=KGSTE; FINI/12/1:=MOUT; FIN(/13/):=KGOUT; FINI/15/):=MENG; FIN(/16/):=KGENG; FINC/17/1:=WISM; FIN(/18/):=KGSM; FINI/20/):=WSTOR; FINI/21/1:=KGSTOB; FIN(/23/):=MCAR; FINI/24/1:=KGCAM; FIN(/26/):=G; FIN(/28/):=GM; FIN(/29/):=KM;

749 749 749 750 750 750 750 750 750 752 752 752 752 753 754

COM.IL(PASB.EXCD): "IFI3AS<0.5"THEN"GOT0'3A1: 'COMMENT' BASIC SHIP OUTPUT; IbEGIN"

SYSAMOS11,14,3);

ODTSTRINGIDS410("8ASIC SNIP DATA')':: SYSACT(DSNI,14,3):

716 717 FIN(/30/1:=M08; FIN(/42/):=LMA; 756 758 OUTSTRINGIOS410('L(M) =')"): FIX(DSN1,3113,BASSI/1/:); SYSACT(DS511,1491); 718 719 FIN(/59/):=SF; FIN(/60/):=PS; 759 761 OUTSTRIWilDS410(13(N) .1)'): FIXIDSN1.3,308AS81/2/1/I SYSACT(DSN1,14,1): 720 721 FIN1/61/):=FREE8; FIN(/62/1:=M-10; 762 764 OUTSTRINGIDSNIO(iT(A)

=1)i):

FIX(DSN1,3,303AS81/3/1): SYSACT(DS41,14,1); 722 723 724 725 726 727 728 729 730 731 732 733 ₇₃₄ 735 736 737 738 739 740 740 741 741

FINI/63/1:=NDO; FIN(/64i):=Wl0; FINI/65/1;=mFM; FINI/66/):=KNI/1/I; FINI/67/11=KM/2/1; FIN1/68/):=K4(/3/); FINI/69/I:=KNI/4/1; FIN(/70/):=KM/5/I; FIN(/71/):=KM/6/I; FIN(/72/):=KM/1/1(KG+B*8/(500*T*CB))*SIN(10/57.297); FIN(/73/):=KW/2/1-1KG+B*8/(500*T*CM*SIN(20/57.297); FINI/74/):=KNI/3/)(KG+B*8/(500*T*C8))*5IN(30/57.297); FIN(/75/):=KNI/4/1-1KG+8*B/I500*T*C811*SIN(40/57.297); FIN(/76/):=KNi/5/1(KG+8*8/(500*T*CE0)*SIN(50/57.297); FIN(/77/):=KN(/6/)-1KG*8*8/(500*T*CB))*SIN(60/57.297); FIN(/78/):=8K8; FIN(/79/):=ANCAR; FINI/80/1:=UCOSTS; 'COMMENT' END OF

OUTPUT PREPARATION;

'END': 'COMMENT' ENO OF JES002;

765 767 768 770 771 773 774 776 777 779 780 781 782 784 785 786 787 788 793 791 792 793 OuTSTRING(DS11Ol'O(M)

.1m;

FIXIDSN1,3.3.8ASBI/4/11; SYSACT(OSN1,1411); OuTSTRING(OSNIO(IV(KM) =1)'); FIX(DSN1.3,30ASBI15/1); SYSACTOSN1,14,1); OUTSTAINGIOSN10(0C8() =0)0); EIX(DSN1,3,3,BASBI/6/:): SYSACI(DSN1,14t1): DOSTRINGIDSN101"DN(TON) ="16): FIzIOSN1,7.0,BASB(/7/));

SYSACT(OSN1,14,1); OUTST4ING(DS1I,IleV0LT(CuM)=1,11); FIX(OSN1,7,0,BASB(/8/)); SYSACT(DSN1,14/3): "IF'BASB(/9/)5.10THEN" "BEGIN' 3UTSTRING(OSN1.0("WSTE(T0N)

1.1)0); FIX(DSN1,6,0,3ASB(/9/));

SYSACT(OSN1,14.1); 'ENO"; "IF'8ASi-1/1J/1>C.56THEN° 'BEGIN' OUTSTR/NG(DSN1,001KGSTEIMI

=11"); FIX(DSN1,6,3,BASBI/10/1);

SYSACI(DSN1,14,1): 'END': "IF'8AS(/11/)>C.51THEN' I6EGINI

741 741 'PROCEDURE'DOPT10 774 796 797 OUTSTRINGIDA1,'('XGSIE(N)

=1)1):

FIXIDSNIt6,301ASB(/11/):: SYSACT(DSN1,14,1): 'END': I 742 742 743 744 745 746 747 748 (IIA,IRA,0IA.ORA,A,(JFvgx,H,GAM0,G,GST,ACCUR,DES3PT,RESULTS); 'INTEGERI'ARRAY011A,JIA; 'ARAAYeIRA,ORA.X,H,GAMO.G,G5T; "ARRAY'ACCUR; 0REAL'AOFV; ILABEORESULTs; 'PROCEDUREmDESOPT; 'CODE'; 798 799 800 8J2 a03 804 505

11FsBAS:A/12/1>10,THEN" '5EGIN" OUTSIRINGlOsNI,01"WUUT(TON)

=110);

EIXIDSN1116,018A53(/121)/I

SYSACT(OSN1,14,1): 'END': 'IF"BASi.11/13/1>C.5"THEN' '6EGIN'

749

896 8313

DUTSTa1iGIOSN1040KG3UT(m)

.'10); FIXIDSN1,6113.8AS81/13//1:

809 810 811 812 814 815

'END':

'IPIBASS(/14/)>G.591HENi

'BEGIN'

OUTSTRING(OSN10(1)XGOUTCM) SYSACT(DSN1,14,1);

'END';

=i/i); FIX1DSNI,6,3,BAS4(/14/1);

876 878 879 880 881 882

OUTSTRINGIDSN10IIKAW) SYSACTOSN1g14,11;

'END';

ilFi3AS8(/21/)>4.5*THEN°

ibEGINs

DUTSTRING(DSNIO('XAM)

=°)'); FIX(DSN1,6,3,BASB(/26/));

=t)e); FIX(DSN1,6,3,BAS8(/27/));

816

iIF.BAS8(/15/)>10.09THENI

884 SYSACTIOSN1,14,1); 817

'BEGIN'

885

'END':

818 OUTSTRING(DSN1O(NENG(TON)

=9)'); FIA(DSN1,690,BAS8(/15/));

886

'IF'BASB(/26/)>J.51THEN'

820 SYSACTIOSN1,14.1); 887

'BEGIN'

821

'END';

888 OUTSTRINGIDSN1019GM(M)

FIXOSN1,6939BASB(/28/));

822

'1FIBASB(/16/)>0.5THEN'

890 SYSACT(OSN1,14,1); 823

'BEGIN'

891

'END';

824 OUTSTRINGIDSN1011KGENG(M)

=1)'); FIX(DSN196,3,BAS3(/16/));

892

11F.BAS8I129/1>3.5iTHEN°

826 827 SYSACTIDSN1,14,1);

'END';

893- 894

'BEGIN'

OUTSTRING(DSN10(6KM(M)

=9)°); FIX(DSN1,693,BASB(/29/));

828

1IFBASH1/17/1>0.59THEN6

896 SYSACT(JSN1,1491); 629

'BEGIN'

897

'END';

830 OUTSTRINGIOSNIOMWENG(M)

21)'); FIX(D5N1,6,39BASB(/17/));

898

SYSACT(DSNI.14,3);

832 SYSACTICISN1.1411/: 899 OUTSTRINGIDSN1011A0S(M)

=1)9); FIXIDSN1,3,3,BASB(/30/));

833 834

'END'; gIFBASB(/18/)>0.56THEN°

901 932 SY$ACT(DSN1,14,1); OUTSTRINGIDSNIO(NA(M)

--111,1; Fix(DSNI,318.BASM/31/1);

835 836

'BEGIN'

OUTSTRINGIDSN1.1(1KGSM(M)

=').); FIX(CISN196,31BASB(/18/));

904 905 SYSACTIDSN1,14.1); OUTSTRING(DALIO('WF(M)

7-0).); Fix(DSNI,30,BAS81/32/1);

838 839 SYSACT(OSN1.14,1);

'END':

907 938 SYSACT(OSN1,14,1);

0UTSTAING(DSN10(1W()

=e),); Fix(DSN1,3,3,8A5B(/38/));

840 841

sIP/BAS1/19/)>0.51THEN'

'BEGIN'

910 911 SYSACT(DSN1,14,1);

0UTSTRING(0SN10(ICAM1)

=IP); FIXICISN1.3,3,BASB(/34/));

842 OUTSTRING(OSN101XGSM(M)

=ill); FIX(DSN196,37BASB(/19/));

913 SYSACT(DS11.14,3); 844

SYSACT(D5N1,14,1):

914

11Ft3A56(/35/)>0.5sTHENI

845

'END';

915

'BEGIN'

846 SYSACT(USN1,14,3); 916

0UTSTRINGIDSN1OMPPER WING TANKS')');

847 OUTSTRING(DSN10('WSTOR(TON)

=IP); FIXOSN116.3,BASB(/20/));

917 SYSACT(DSN1,14,3); 849 850 SYSACT(DSN1,14,1); 0UTSTRING1DSN1tilKGST0B(M)

=1)9); 1IX(DSN1,6,318ASB(/21/));

918 919

'END':

OUTSTRING(DGN10(WW3(M)

=8)1); Fix(Dsm1130,8AS81/35/1);

852 SYSACI(DSN1914,1); 921 SYSACTI3SN1914,1); 853 OUTSTRING(DSN10(8XGSTORIM)

=')'); FIA(DSN1,6,3,BASB(/22/));

922 OUTSTRINGIDSNW(WWA(DEG)

=')');FIX(DSN1,3,3,BASB(/36/));

855 SYSACYCIAN1,1491); 924 SYSACT(DSN1,14,1); 856

9IFILIASBC/23//,10.06THEN1

925

'IFIBAS8(/37/)>C.51THENI

857 obEGINg 926 'BEGIN' 858 OUTSTAI4G(1354101eWCAR(TON)

o1)1); FIA(OSNI,6,0,BASS(/23/));

927 SYSACT(OSN1,1413); 860 SYSACTIOSN1,14,1); 928

OUTSTRING(DSN10(8LOWER WING TANKS')');

861

'END':

929 SYSACT(DSN191413); 862

IIFIBASM/24/1,0.56THEN*

930

'ENV;

863

'BEGIN'

931 0UTSTRING(DS41019LWBIN)

i)'); FIX(DSN1,3,3tBASB(/37/));

864

0UTSTRING(35410MGCAR(M)

=1)'); FIKOSN1,6,3,BASbl/24/1);

933 SYSACT(DSN1,14,1); 866 SYSACTIUSN111411); 934

0UTSTRINGI05N10('LWA(0G)

.°)1); FIXIDSN1,3,318AS81/38/));

867

'END':

936 SYSAC1(05N1,1493); 868

11F'BAS8(125/)>C.51THENI

937

OUTSTRINS(0N1016CARGO AATCH COAMINGS').);

869

'BEGIN'

938 SYSACT(DSN1,14,3); 870 OUTSTRING(DhNlyileXGCAR(M)

=1).); FIA(DSN1,60.8A56(/25/));

939 OUTSTRINGIOSN191(°LUL(M)

.1)1); FIAIDSN1.3,3,BASBI/39/1);

872

SYSACTI0SN1,14,1);

941 SYSACT(DSN1,14,1); 873

'END':

942

OUTSTRING(DAIO('LUS(M)

-7-')°); FIXIOSNI,303,BASBI/40/1);

874

11F1BASB(/26/)>0.5ITHENe

944 SYSACT(35N1,1491); 875

'BEGIN'

945

OJTSTRING(OSNIOMUAM)

="16); FIA10SN123,3,BAS8(/41/));

947 948 949 950 952 953 955 956 957 958 959 960 461 963 964 966 967 969 970 972 973 974 975 976 977 978 979 981 982 984 985 987 988 990 991 993 994 996 997 999 1000

132

1003 1005 1006 1008 1009 1011 1012 1314 1015 1017 1318 1319 1020 SYSACT(0SN11114,3): OUTSTRING(OSNIO('.ENGTH OF THE

ENG.ROOM')'):

SYSACT(DSN1,14,3):

OUTSTRING(05N1.1('LMA(M)

=6)8):

FIX(OSN1,313.8AS3l/42/:):

SYSACT(0SN1,14,3):

OUTSTRINGIDSMW('LFP/L()

=1)1):

FIX(OSN1,3.3,80,68(/43/)):

SYSACT(OSN1.1493):

IIFIBASB(/44/))0.5:THEN'

'BEGIN'

_{OUTST1ING(OSN10('UPPER}

HOPPER

TANKS')');

SYSACTIUSN1114.3):

'ENO':

OUTSTRING(DSNIO(INU0HPIPIECW).):

FIK(OSN1,3909BAS8(/44/)):

SYSACT(OSN111491): OUTSTRINGIOSN1.619UBUHP(M) 2211'):

FIX(DSN1.3,3,BAS8(/45/));

SYSACT(OSN1.14,1): OUTSTRING(DSNIOI'LRUHP(M)

.11');

FIX(DSN1,313.8AS3(/46/));

SYSACT(DSN1.14.1): OUTSTRINGIOSNIOI:DUHPIMI

=111):

FIX(OSN1,3,3,8:68(/47/)):

SYSACTIOSN1,14,1): 11F:BASB(/48/)>C.51THEN: 'BEGIN'

SYSACT(05N1,14,3):

OUTSTRINGIOSNI0(1104ER

HOPPER

TANKS')');

SYSACT(OSN1,14.3):

'ENO':

OUTSTRING(DSN10(°N0LHPIPIECW).):

FIK(DSN1p3,001ASB(/48/)):

SYSACTOSN1.14,1):

OUTSTRINGIOSNIpi(lUBLHP(M)

09)1):

FIX(DSNI,3,3,BASB(/49/)):

SYSACTIOSN1.1411):

OUTSTRINGIDSNIO(:OLHP(4)

=i)e):

FIA(DSN1,313,BASa(/50/)):

SYSACT(DSN1,1491): OUTSTRING(0SN1.O(10LHP(M)

=6)61:

FIX(OSN113,3,8AS31/51//):

SYSACTIOSN1.1493):

0UTSTRING(OSN10('S1ERI()

9)1);

FIX(DSNI,3,3,8AS8I/52/1):

SYSACT(O5N1.14,1):

OUTSTRING(DSV10('STER2I)

=1)8):

FIX(3SNIg3.3,8AS3(/53/)):

SYSAMOSNI.14.1): OUTSTRING(DSNW('LKF/L()

=)'):

FIX(DSN1.3.3.8ASR(/54/)i:

SYSACT(DSN1,14.1): OUTSTRING(OSNIO:ISIG(K0/MM2)

=°)°):

FIX(OSN1.3.37BASB(/55/));

SYSACT(OSM1914.3); OUTSTRING(OSNIOlyN11/MINI

=')'):

FIXIOSN113.3.13ASB(/56/)7:

SYSACT(DSN1,14.1): OUTSTRING(OSNIOCISEACUN()

=111):

FIXOSN1,3.3,BAS0(/57/)):

SYSACT(35N19144,1):

OUTSTRING(OSNW('ETAM()

=8)'I:

F1X(DSN1,3,3.8AS3(/58/));

SYSACT(DSN171411): OUTSTRING(DSN11,0('E1AWL() =8)11):

FIAOSN193.39BASB(/59/)::

SYSACT(OSN1.14,1); OUTSTRING(CISNIO(iPilHP)

=1)1):

FIAIDS.N1,6,0,8ASR(/60/)/:

SYSACT(DSN1.14.3): SYSACT(DSN1,14,3): OUTSTRING(OSNIO(IEKPER.

COEFFICIENTS')');

SYSACT(DSN1.1493):

1021 1322 1023 1025 1027 1029 1331 1033 1035 1037 1039 1041 1042 1042 1044 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056

167

1058 1059 1060 1061 1362 1063 1364 1065 1066 1367 1368 1069 1070 1371 1072 1073 1374 1075 1076 1077 1078 1079 1083 1081 1061 1081

IFOR,1:=1'STEP1OUNTIL:25'DO' 'BEGIN'

1IF'I=46THEN"GOTO.Lal:

$1F°I=7°THEN"GOT0:LE1:

'IF'I=10'THEN"GOTO'LE1:

1IF11=13:THEN"GOTOILE1:

81FII=18:THEN"GOTOILEI:

'IF'I.211THEN"GOTO:LE1:

61F11=23'THEN"GOTOLEI: OJTSTRING(DSN1Ol'EXCOMII:

FIX(DSN1.2,0,1):

OUTSTRING(DS410(1)=1)8):

FIK(DSN1,2.3,EXCOUIM:

SYSACTOS41.14.11:

LEI: 'END':

OUTSTRING(OSNW('EXPLA)ATIONO)'): SYSACT(DSNI.14.3):

OUTSTRING(DSNW(I

IF

EXCO(I)=1.0000

THEN FRESH

CALCULATION')');

SYSACT(DSN1.14,3): OUTSTAING(DS41,6(6EXC01COEFF. FOR

POWER')');

SYSACT(OSN1,1491): OUTSTRING(OSNI0(6EXCO2COEFF. FOR STEEL

WEIGHT')');

SYSACTIOSN1214,1:: OUTSTRINGIDSNIO('EXCO3COEFF. FOR KG

STEEL')');

SYSACT(US31,14,1):

OUTSTRIAG(C00410(IEK,05COEFF. FOR OUTFIT

WEIGHT')');

SYSACT(OSN1.14,1); OUTSTRING(OSNIO(IEXC06COEFF. FOR KG

OUTFIT')'):

SYSACT(USN1,14.1): OUTSTRINGIDSNIO('EXC08C3EFF FOR MACH.OltIGHTi)10:

SYSACTOSN1.14.1):

3UTSTAING(OSN1Ol'EA..09COEFF. FOR KG

MACHINERY')');

SYSACT(0SN1114.1): OUTSIRING(DSN1.111EXC011COEFF. FOR

L.S.WEIGHTIA1):

SYSAC1(OSN1.14,1):

OUTSTRING(DS410('EXC012COEFF. FOR KG

L.SHIP').):

SYSAC1IDSN1.14,1); OUTSTRING(DSN10(*EXC014COEFF. FOR ENG.ROOM

LENGTH')');

SYSACT(DSNI,14.1):

OUTSTRIAG(DaN101.EXC015COEFF. FOR GRAIN

CAPACITY')');

SYSACT(USN1.14,1): OUTSTRINGIDSN1.10EAC019COEFF. FOR CARGO

WEIGHT')');

SYSACTIOSNI,14,1);

OUTSTRINGI03N10(1EXCO20COEFF. FOR KG

CARGO')');

SYSACT(USN1,14,11: OUTSTRING(USNIte('tXCO22COLFF. FOR KG

SHIP')');

SYSACTIO.SNI,14,1);

OUTSIKINGID2IOliEC:024COEFF

FOR GM

SHIP')');

SYSACTIGS41,14,11;

3UTSTRING(DSNIO('EXJ25COEFF.

FUR KM

SHIP')');

SYSACTIUSN1,1411); :COMMENT' EN1U OF BASIC SHIP OUTPUT I

BAI: