Weerstandsberekening van schepen volgens de gegevens van H.E. Guldhammer en Sv.Aa. Harvald

(1)

TECHNISCHE HOGESCHOCL DELFT

AFDELING DER SCHEEPSBOUW- EN SCHEEPVAARTKUNDE CENTRALE WRKGROEP WISKUNOF.

Rapport No. 14

WEERSTANDSEEREKENING VAN SCHEPEN VOLGENS DE GEGEVENS VAN H.E.GULDHAMMER EN Sv.Aa.HARVALD

Ing. A.P.de Zwaan

september 1979

DeUI Uriivrsity o Technology

Ship Hydrorn.chanics Laboratory Mekelweg 2

DeUt 2208 Netheçanth

(2)

Inleiding i

1, Algemeen. 2

1.1. Programma gegevens 2

1.2. Doel 2

Organisatie van het programma. 2

2.1. Gebruiksmogelijkheden 2

2.2. Gebruikte methoden en korrekties 3

2.3. Verkiaring van de gebruikte symbolen 5

2.4. Stroomdiagrammen 6

2.4.1. Programma "GULM" Polynoombenadering 6

2 . 4 . 2. Programma HGULWU Wee rstandsberekening ₈

Gebruik van het programma. 10

3.1. Programma 'GULW" 10

3.2. Invoer ₁₀

3.3. Uitvoer ₁₁

Literatuuriijst 12

Append. ix

Listings van de programma's.

(3)

Inleiding

In het rapport wordt de beschrijving gegeven van een korruouter

programma met behuip waarvan de weers band van schepen wordt

be-rekend voor dienstkondities.

Er wordt gebruik gemaakt van de gegevens uib het rapport:

Ship Resistance

Effct of form and Principal Dimensions (Revised)

by

H.E.Gulcìharnmer and Sv.Aa.Harvald

1974

Akademisk Forlag

Copenhagen

Het programma is geschikt voor schepen waarvoor geldt:

a)4

L/1/3<

8 en Prismatische coefficient

0.5CP0.8

b) 4. 5

6 en Prismatische coefficient O.

5 CP 0.85

Op ze t

in[ij zijn grafieken

gegeven voor de weerstandscoCfficient CR

op basis van het froude getal F=V/'jvoor de prismatische

coefficient

CP=O.5,0.55, 0 85

en

L/''3=4,4.5,..

8

De krornmen zijn opgemeten voor de froude getallen 0.14,0.16,..

o 44

Met behuip van deze meetwaarden wordt een polynoom benadering

gegeven volgens 0Theilheimer"

[2]

(programma "GULM")

De berekende coCfficien'cen worden weggeschreven op schijf.

In het orogramma "GULM' zijn de fig. 10 en 11 [i] voor de

kor-rektie van XE ingevoerd.

Het programma UGULWU berekent da weerstarid voor dienstkondities

met de korrekties als gegeven in

[1]

(4)

1 .Algemeen

1. 1 .Prcqramma gegevens

Taal Algol-60

Geheugen 256 K (GULM) ,256 K (GULW)

Rekentijd 12 sec(GULM) ,6 sec (G[JLW)

Naam van de programma1s :GULM en GULW

1.2. Doel

Berekenen van de weerstand voor dienstkondities voor een zo groat mogelijk scala van scheoen.

2.Organisatie van het programma

2. 1. Gehruiksmogeiijkheden

Ret irocJrarnma !GULMH is eenmalig gebruikt voor de polynoom be-rekening van de fig.1 t/m il uit i]

De gebruiker rekent alleen met het programma "GULW".

Ret orogranima berekent voor de opgegeven dienst sneiheden de weerstand van het schio met de voiqende korrekties voor: a) CR

B/T - Breedte diepgang verhouding.

XE - Langsscheepse ligging van het drukkingspunt. Scheepsvorm. Buibsteven. Schroefas en schroefasuithouders. b) _CF 1) Huid en aanhangsels. c) CA Windweerstand. Stuurweerstand.

Ret programma produceert een lijst van weerstanden voor

O.14<FO.44 en een lijst van weerstanden ';oor een door de

gebruiker opgegeven aan bai sne lheden.

Op deze waarden geeft de gebruiker een extra toesiag voor dienst kcndities,volgens een lijst die het programma uitvoert.

(5)

-2--Voorschip Achterschip

Extreme U

2.2 Gebruikte methoden en korrekties

Het Drogramrna maakt gebruik van de krommen voor de weerstands-coefficient CR zoals deze gepubliceerd ziln in de publikatie van Guidhanmier en Harvald

In dit rapport zijn tevens korrekies gegeven voor CRICF en CA. De korrekties voor CR zijn:

B/T-korrektie 1O3CR=1O3CR(B/T_25)+0.16(B/T_2.5)

XB korrektie 1O3CR

_{+10CR/XBtXß-XB}

I

-

standaard standaard

zie fig.10 en fig.11 i]

Geen korrektie wordt toegepast ais XB4XB Korrektie voor de scheepsvorm

Extreme V

De korrektie geldt voor 0.2OF 0.25

Korrektie voor de buibsteven

F 0.15 0.18 0.21 0.24 0.27 0.30 0.33 0.36 C n p 0.2 0 -0.2 -0.4 -0.4 -0.4 0.2 0 -0.2 -0.3 -0.3 0.2 0 -0.2 --0.3 -0.3 0.1 0

-0.2

De korrektie geidt ais

Als 0<ABT/Ay0.10 wordt de korrektie evenredig aangepast. De bovenstaande korrekties zijn in poiynoomvorm opgenornen in het programma.

Korrektie voor de schroefas en schroefasuithouders

- n o. -

van

0.50

0.60

O . 70

0.80 standaard

3 -O . 1

+0.1

+0. 1

-o

. i

(6)

De korrekties voor CF zijn: Huid en aanhangsels

L loo m io3 _CA=O.4

L= 150m

=0.2

L = 200 in =0

L = 250 m =-O.2

L > 300 m =-O.3

Deze korrekties zijn in polynoomvorm opgenomen in het programma.

Weerstands verhoging voor de aanhangsels

CF'=CF x S1/S met S is het nat opperviak van de huid en S1 het nat opper-viak van de huid + aanhangsels.

Het nat opperviak wordt berekend m.b.v. de methode voigens Holtrop en

Iennen [3].

Wind en Stuurweerstand

Windweerstand

_103CA=0.07.

Stuurweerstancl 1O3C=0 .04

Toesiag voor dienst kondities

Noord Atlantische route,00stwaarts resp. in zomer en winter Ncord Atlantische route,westwaarts resp. in zomer en winter Pacific route

Zuid-Atlantic en Australische routes Oost-Aziatische route

Deze waarden worden uitgevoerd door het programma.

-4-15-20 % 20-30 % 15-30 % 12-18 % 15-20 %

(7)

2.3 Verklarinq van de gobruikte symbolen S ymbool SAPP AB T S RUD s AX B CA CB CF CM CP CR Cwp DCR FNV L LD RT T vs LCB DELTA RHO

Programma Dirnens.ie Omschrijving ¿ m 2 ti 2 nl 2 rn 2 m ra 3 m 3 kg/rn

Nat opperviak van de aanhangsels. Opperviak van de bulb op de

voor-loodlijn.

Opperviak van het roer. Nat oppervlak van de huid. Grootspant opperviak.

Breedte.

Weerstandscoëfficient voor het schaaleffekt.

Blokcoëfficient.

Coefficient voor de wrijvings-weerstand.

GrootspantcoCfficient. Prismatische coefficient.

Coefficient van de restweerstand. Waterlijncoëfficient.

Korrektie van CR. Fraude getal VA/ x L

punt t.o.v. L /2 in van L. no

Waterverpiaatsing.

Soortelijke massa van water.

"5-rcL/sec2

ri-i

Versneliing van de zwaarLekracht. Lengte op de waterlijn.

m Lengte tussen de loodlijnen.

Verhouding

L/13.

N Totale weerstand.

ra Diepgang.

rn/sec Sclieepssnelheid.

iS Lengte ligging van het

drukkings-AA ABT AR AT Ax B CA CB CF CM Cp CR Cwp dCR F n g L L pp RT T V XB V (o

(8)

2. 4 Stroomdiagramrnen

2 . 4. 1 Programma HGTJLMU polynoombenadering

flT7e k laratiesj L LTMAX=0

.J

P

-6-(.5

) B LCP=O.L [1FN0.14 (.02

3TT7'i1

Lees de meetwaarde van CR in,die behoord

bij L/13,FN en OP

nee /

,,

ja

De coëfficienten COP van dc polynomen

CR zijn bepaald voor L/v1/3 _{en FN met OP als} variabele

FNMAXL/1/3Cp=FN

Aantal meetpunLen NL/l/3FN=l

Y(1)=CRL/,1/3NCp;X(1)=CP

L CP=X ( 1) +0 . 05 (0 . 05)

E

Lees de meetwaarde van CR in.die behoort bij L/Vh/3,FN en CP

L/V1/3

n ee

CR/71/3Fp bestaat niet;verhoog CP meL

LQ

3 nec

j a

De laatste meetwaarde was CRJ./,l/3

FN CPU 05

FN AXL/vl/3 _FN

_cP=;NL/l/3

_{FN_NL/y} 1/3

(9)

C'-, U) (1) 1/3 L,V ,LN T CCPL/ 1/3 _FN ;XCPL/91/3 _FN =x I

=o:--:

-j j / Schrijfopschijf:N(1:9,1;16),XCP(1:9,i:16,i:8) / en FNMAX(1;9,1:8) (9 is het aantal L/9"3116 is

/

het aantai froude getailen FN's en 8 het maximum

/ aantal CPs).

/ Lees de cofficienten CLCB van de líjnen uit fig.

u

[i) voor de korrektie van CR voor LCB

Bepaal de cofficienten van cte polynomen op basis van CP: CDLCBFN

_i8

Zet cleze waarden op schijt.

-7-U) U) :j U) r C') -.

Bepaal de coëfficienten C(1:NL/1/3FN)van het polynoom op basis van X en Y

(10)

2.4.2 Programma "GULW" ieerstandsberekening

V

Dekiara aies

L

Lees de coëfficienten van schijf. Lees de gegevens in voor het schip.

Bereken:CB,AX,AET,L/3

en nat opperviak

schip S van de huid..

Zet

de

waarden voor FN=O.1470.16 0 44

in

de array Xi.

L/34 (.5

) 8

i

IFNiL/1/3=minFNiAXL/V1/3CprFNMAXL/9i/3Cp+OOS)

n e

j L/913 geen waarde

van CR

voor gegeven CP

bekend.

n

FN=0. 14 (0.02) FNi1 1

rn1 :i) aantai CP's beh.

bij

en FNí

XXCPL/ 1/3

_,FN m;Cnh CCPL/71/3 _,FN,m Y1T/?1/3 (CP)

79

I EET= 014(0 Q2) O

4fl-L/74(0.

[FNiL/î/ 3-FN

n e e

nee

-

_L/r1-ENT(L/V

)=O

---

---r:

(11)

z

k1=k1+1;X(ki)L/Ç13;Y(k1)Y1T, 1/3

i ,FN

Bepaal de korrektie voor CR,CF,CAA en CAS. Bepaal de weerstand van het schip voor

FN=O.14,O.16, 0 44

Bepaal m.b.v. interpolatie de weerstand van het schip voor de opgegeverì srielheden.

9

k3=k3+1;Bepaai de coëfficienten van het o1ynoorn voor CR,die behoren hij FN op basis van

L/13.

Bepaal de waarde CR(k3) ,die behoort. bij de waar-de

L/3;kiO;FNV(k3)FN.

(12)

3. Gebruik van hot orogramma

3.1. Procrramraa "GULW"

Het orogramma staat in de hibliotheek SBAL.LIBCWW van de "Centrale Werkgroep Wiskunde".

Het programma berekent de weerstand voor opqegeven sneiheden en de standaard fraude getalien 0.14,0.16 ...0.44 en maakt gebruik van de coëfficienl:en,herekend door het progranuna "GULM", van de polynornen voor CR.

Deze coëfficien.ten staan in de hibliotheek _{SBALRES.SBALCWW.CNTL,}

onder de naam GULCS van de "Centrale Werkgroep Wiskunde".

3.2. Invoer

Benaming Dimensie Oprnerking

CP - Prismatische cofficient.

L rn Lengte op de konstruktie waterlijn.

B rn Breedte.

T rn Diegang.

DELTA rn Waterverolaatsing.

CM - Grootspant cofficient.

CWP - Waterlijn cofficient.

BT m2 Opperviak van de bulb op de voorloodlijn

S m2 Nat opperviak van de huid.Als S=0 wordt ingo-voerd,wordt doze berekend door het programma.

SRUD m2 Opperviak van het roer.

2

SAPP in Opperviak van de aanhangsels.

LCB Lengte ligging drukkingsunt t,o.v. het grootspant(L/2) in procenten van L.

(voor 4- achter -)

RHO _kg/rn3 _{Soortelijke massa van water.}

OP EEN NIEUWE KAART INVOEREN.

EFJFB _{Extreme U-vorrn vcorschip:'TRUE' anders 'FALSE'}

EVFB

-

Extreme V-vorrn voorschip:'TRUE' anders

'FALSE'

EUAB - Extreme U-vorm achterschip;'TRUE' anders

FALSE'.

EVAB

-

Extreme V-vorm achterschip:'TRUE' anders

'FALSE'.

(13)

-VS (1) VS (2)

VS( NV )

rn/sec

of

knopen

Scheepssnelheid.

DE VLG4DC T.ELADEN 'YR DI

TKOJJITIZS 3SD2

ND3RD ATLANTI3CHE

0'JTE,3STJATS RES?.

I'1 Z0ME

NOORD ATL4NTISCr1 RCUT STWÀART.S RESP. IN ZOMER.

PAÇIIC. 0UT

ZLID 4TLNTIC EN AUST&LI3CH

ROUTES

OÚST AZIATISCE ROUT3

Krntrole Invopr,

.5'5,i9G.7 30, . 8,-, 31232,.965,.708,11.ti.,0,30,0,-.03,1025

'FALSE',

'FALSE','

FM.SE' 'FALSE'

₁₁

--3,1),20,

21,22,23,2!i.,25,26

ANGE ÁDN EN WINTE.:1S-?C EN WINTER:2C-3C :15-3C' 12-16 : 15-2O 3.3.

Uitvoer

+I1.9554

l3.633

15.3712

+17.0791 :ERSTArIO N +312356 + 406s 7 s +525062

+534915

FN

.1400

.1600

13t3

2CO3

SNELHETD N/SFC

+5.149

+7.0334

7cC7

+,7355

127C

+9.b4-1

+fl3.7870 + 770 C 25 +.2'+0C

+1Q.54

+2C.453

;- 93 331

+.26C0

+114212

+22.2329

129 C34 S.

+2900

+12.2997

+3.91CS

+1 3690

+.3000

13.17E3

+25.4137

+.320

14.0S6

+27.326G

+363 958 3ri

+.340C

+14.S354 .C345

+4919509

4.3600

+15.8139

+O.7424

+6347O5B

-.38C0

1ó.492E

+32.4303

+94544

+.4i00

+17.5710

+34.153

+141 54183 +..42C +13.44-5 +35.352

+..225

+9.7733 .2342 ±1C.2S8 .,_. + 3 Ç 033

+2459

+10.8024

4i.0000

+103 19 i-.2576 1.3153 + 123 959 5

+.293

11.3312

+23 .0303 +152 1044 4431O

1Z.34E6

+ 133 9416

.228

+12.3603 +25.0303

+.345

+13 .3744- +26.CC.CC

2339533

Vervoig van do invoercjegevens

Benarning

Dimensie

Oprnerking

OP EEN NIEUWE KAAIT INVOEPEN.

NV

-

Aantai snelheden waarvoor de weerstand meet

worden

berekend.

NV<0 Sneiheid in knopen invoeren. NV>0 Sneiheid in rn/sec invoeren.

(14)

Literatuuriij st

Ship Resistance

EffecL cf form and Principal Dimensions (Revised) by

H.E.Guldharrimer and Sv.Aa.Harvaid, 1974

Akademisk Forlag, Copenhagen

2 _{Theilheimer F. and W. Starkweabher,}

The fairing of ship lines on a high-speed Electronic Computer, D.T.M.B. Report no. 1474, 1961.

3 _{A Statistical Power prediction Method,}

by

J.Holtrop and G.G.J.Mennen,

International Shiohuilding Progress,

Sliipbui iding and Marine engineering i1cnthly,

t--'

Vol 25,no. 290,oktober 1978.

(15)

-Programma "GULM" -. 2 2 -, t) 1T:G ' IL ) ,IFf C,P , , i , ;

(JCI00004G

CCCCC.

OCCLCGu

CCCC,0070 COCCCCEC

7

-,t*1\i1_ :,

t,

'\/} ,'INT,'

4DRAYPK,.(/1.N/) ,('CCCC(9C 11

ucccc1c

it; LcGcc?11c. l&) ;

cT''FY'J -1'S1.1J\iIL-.?.'(/I,J/) =IIFt4R.(X(/I/})>I.

CCOC1:C 1 [H\h1tLLStO CCOC11C 19

cCCC1C

2i

'THEN' (X(/ I/)-X(/J-I) }3'LSEC ;

1 CCCCC1( 24 25

I\)*?'<cJ,Jj

crH,

CCCLCI7v CO,C1bC 2

'G1N''INTEG'1:'k:ÁL'YY;

oc,ccc1s 2 33

Y\. =(/1/)+(/f)X)+.(/i)' \X>+2(/'-/)XI*XX*XX,

F<'I ='TP'1'U1IL1

xx-x(,I-_,)>'-:'f-s:' :5

yy:=yy+(fI/)(XXX(/I_/)):.;ThEiL?:=yy;

ccccc::o t bJ T1 ILr,

'IL

'T'1'UNT1L'15'J'

Lt. LCCL C C C CC Z SC CC CCC2 SC 41 'FCR'ILE;'STEP'1'JNT1L'16't,' )CCCC27C 42

'C1i'I:1:

44 2tIF'!:=7tSTPt1uN1TLZZtDOt

(CCCO2O

45 * t 'tS Li;

R'ADL ,.( /IL, iF',T/) }

; IF 'CR

(fILE) ,I FN,I/)=-1 'TN'

t

:î1

'r i+1 (CCCC1C

CLCCC

'T ''T...'L3,F1AX(/1LD,I/) .«''CLCC )?3. 55 JFN;':(/1L)1IF\/):=t;Y(/1,'):=.Co(/ILD,IFN,I/);II:I+i; CCCCC34 59

X(/1/}:=.G3-1;

CCCCC3SC 60

II:cï):_ITscTpti?uTIL;i7p:cI

f

( Il':\(:,O 1/IL. ,IN,1»/)), '.

C.CCCt36C

i(/IL),iCN/):=(/ïLD,IrN/)1;FIX(1,2,C,(JrLJ,IFN ¡H;

C C C C 02 CCCCO39C CCC C. +C C 7

1(!(/ILD,IFN1/)/)..C(/ILJ,IFN,ICP/);

CCCCC41(, C0CCA2C 71 L2: CÜCCC43C COCCCi't4C CC C. C 77

'iN';1:iC

(OCCC4ÒC 81 'ENr)'; CC'CCO47C 82

L: '[ '

CCCCC4eC CCCCC4'L' C CC CCC C O? '-JII-_7':TEP'1'UFTIL'2'DO' CCC(C510 _.cL>IF, CL¼LLS2C

(16)

i

"4

--'I'»'t/T--i I)

5,Y(/ICP-1.,/)>;

'.V' cj:Rr-'(y,y,7J)

CCQCO53C (C(C54C CCCCC55C CCG55C 1C)2

ccccci

13 'r';JTARAY(3,DLcA);:L'TJRAY(Z,CLC1;

CCCCC50 106 CCCCO590 1:7

':iG1N4turSTING(i, (!L/(3LTÁ)**1/;)t);FI(ii1y1,1LD.);LINL.(1,1 CCCC6CC

111 ) ; CQCCO61C

ill

GUTT-\ING(i,(' ...0

J.55 C.6C C.u5 ...70

0.75

C.00

COCCGÔ2C 11

.95')LI:(1,1H

CC.c63C

113

rCtTFN:=7TD1SUNT1Lt2t,Dü1

CCCO64c 114

.r,I

' 'F' 1

T:r* i VU TrL*N(/LL ,IFN/ 'j

L1 CÇCCC6ÌL

1i'

OCCCC67C

10 CCC6C

i'3

14 'F'I

_LrttL(T1L,17,0,

CCCC7CC

16 >.G1UTHNV

CC.CCO71 126 LCOCOI2O 1 ,C V CCC t. i C 13C

FIX(1,?,1,1HEIL3(C,X,(/ILD,IFN/),..5*IC)/.3625);

CCCCO74Ü L

i

L'

1,

'ENJ' ;

CCCCC7EÇ 1?.;

'ND ; 'COTí' LÇ

CCCG77C 137

SNG:uU1STI

(1(FCUT'}');

OCrCCC78C

19 L:'NL

CCCC73

(17)

Programma "GULW" Pf ),Y1,CLO(/3:1',1:16/),X,Y,NV,DO,Ci/1i7/),C,X1(/11óI),X2(/19/), CCC0005L _{CL,LCB(/72,1:7/ )} -1EsOL5AN'JFB,EVF,EUAe,EVl.F3 :3 I INT ' IFNi (/5:16/); 4 S

sRuD,CM,C,CP,R,I-tJ,LC,51;

5 'ITEG:s'FN,ILD,,K1,K:,:NTLo,ICP,NV; b 9

$EGISIJTGRhI,J;IARR4Y*A(/1:N,1:N!)

13

PR0cUREP1VLU;

'COD ; 'PROCEDU'R 17

*Ffl11:1STcPP1UNT1 ''_n'

18 21 21 23 23 26 'END'2RflCDUR C.OTH; REAL' "'RO

DU'THEILP (O,XN, XX) ; ' ASSA

31 -' ' INTGER' 1: .EAL'YY; -¿4 VY: (/1/ )+C(

/ /)i>XX+C(/3/ )*XXX+C(

35 37 39 'END'PR0CEDUR THEILP: 40 N,xCP,ccP,N AX1

LO);ADR,CP

S°, LOB Rdu); 42

-- LINE(C,1);AD(O,EUF.,EVEUAEJAJ

IF' C>. t LK'+. 5 LD> ) 'T 'GOTO LF 'FÛR'IFN:75TEP'l'UNTIL'2'DO'Xi(/IFN-

_{Â:=ENTIR(P/..CCC1);31:=.+1;}

-FOR 'ILD :8' STP l'UNI IL' 16

'00'

'BEGI(J

'Nl: F;iAX (II

LD

AI) ; 'IF' F'1>FM,A

'1 I F'iL( /ILD/

)=e 'T;L

0TC'Li;

'FJR IFN:

7 'STEP 'i

''iNri

L 'I FN1(

/ I LD/)

'I'F.JM:1'ST.2l'UNTIL'N(/ILD

'GIN'X(/i/):=XCP(/ILD,IFN,/);C

'';yi(/IL0,IN/:THEILP(C,X

'ENO

Li:

ND';

<1: K3 :O ; ; IFN : 7 t s ' i' UNTIL' 22 'DC -3EGIN' 'FO

ILT:=3'ST ?'

1' UNTiL '16 'Dci

'B'GIN' 'REAL' T;T:=LFNI(IILD/)IFN; I

IF ILD'NTL C' THEN ' GCJTO

'ND';K1:K1+1;

X( /K1/).5*ILD; O 2 000001 (0 j/):='ÌF'X(/I/)X(IJ-2/)>'--5 300Cc-180 CO CCC 190 00000200 CL L C 1 C

y,X;'REALXX;'iNTLGN;

-CCCCi2C C0000230

/;/)XXXX*XX;

CCtjCC

(/1-2/

)>'S 'THEN' CCCCC2. COCO 026 CO 0002 70 C C O G C 280

CCCC290

C C C 003C C X*.O1; OCCCO3IQ CC001); OCOCC320 cul; (CCC34C- CCC CC 3 50

-CCCO03O

-00000370 x(fILEE1f)'THE'i' CCCCO36C rI:R(1/.3Z+.GC00C1); C000039C CC 0CC 420

'00'

000(0453

IFN/i'DO'

COCCC4C

t/-/):=CP(/ILD,IFN,M/);

CO300470 ,N(/ILD,1FN/),CP); 000004EC

000510

C 30CL E) 3 L CCC CO 4C 00 0CC 5 50 CC 060560 F"í<O&ABS(T)>.000001'THEN' Cc-000570 'L3;'GCTO'L2; COCOLSECì CCC 00590 OC 0Z06 10

;'

iTE(.E.'N;

Cc-000120

;'ITER''ARÂY'PK,R(/1:N/);0CCCC13C

_T.UB;'CCDE'; 0000014G CCC C C 150 (/I,J/):='IF'A3S(X(/I/))>'-6 CCcCcl&C Th{'EN'i 'ELSE 'O; D: L,B,T,CLTA,CM,CWP, -C-COCO-060

-CCC0CG7i CC COLO 8L 00000090 OÇCC01CO 000001 10

(18)

8g Y(/U/):.yl(/1LD,1FN.'/); CCCCC62C 90 L2:

'fl;K3:K3+l;

COCCC63U 93 COOCO&4C

FNV(/K3fl:=x1(/IFN-6/);

OOCOCi6óO 97 CCC(C670 98

1(0TO'EINDi;

CCCCC7IC 99 CCCC3720 103

.03*ICP;

cCGw07C

13 FAK:i/4,;'1F'FAK>.l'T'K:1'ELS('AK=FAK'1C'

00000740 13

L/1.1',S

T(CH( .4+0GCL

C CCCC07.0 110 COCCO7Oe

ll

00000810

lic

CGC0C2C 122 124

0000C350

130 CC0C'87G 13&

00000910

140

C0000940

141 000LC95C 147

i1

00CCC99C::.

-U0CCiC30:1T 153 GCOO1CIC

00 oc io o

C C C 10 30

00001040

172

GCGC150

17 CCC1c,6

0000109G

IS?

00001100

00001130

204 030CI1ÔC 217 000012CC 112 .4EETÀ'flP)');LINE(l,l); C'CC0C79C 114 ')');LI.\(l,l); CCCCCSCC 11% CCC0C.S3C 12S 00000SÔC 131 00000E8C 133

00000850.

1?4

C0000903

140 C000093C

13

(CCCC»7ôC

S:=S+SRUD;H:'S:.RHJS*.)l;

0UTSt1NG(1,

FN NLH'TL)

DUTSTIGC1,'(t

M/S.0 KNOOP

b0IN

A'AY'

T ( /1:K/ ; FOR

:'11

STtP' 1. W'JT IL iO

'j

f(/1

K/ },

; I1F'LCB>DLCFTH'4t DLC3: =LCBDL0b'EL

ErD,

L5:

'IN''IN/(/K1/)<-.27'Tr4EN'

tLSE1FtCP<=.7tT.N

' ßEGIN ' I F NV 1K 1/) .2 7 T HN tBEOINIV 178 OCOC 1070-130 C0001O8C 138 0000.1.110

11 C0CC112C...

201 CCcc'114G 202

C00C11t

20S

00001I7C

209 C00011EC 2.14

00001150

218

C00012IC

tNDIIFLttIFFN1\f(/K1/)<.3tTHcNt

?'flt -.!G1NttIFtFNV(/Ki/)>.18&FNV(/K1/)<=.2L1THNS

'8G1Y:=-20FNV(/<1 /)/3+2 .8-202;

''y<IT'

END'; FNV 1(1/ ) < 1S FNV /K1 I) >=

.15

'THeN 'Y = 1 10*FNV (1 Ki/ ) /2. 5) +

(10CP4E)*( 10FN'.'C/K1/) /3.4) 'EL

'Y: =.9CP;

'END' 'LSC'

'ßcGINi'IF'FNV(/K1/)>=.18»N1(/K1/)<=.2i'THEN'Y:=-2C*FNV(/K1fl/3

&FNv(/Kl/)>=.15tTHEY:=-10*FNv(/K1/)/3+.6tELSE'Y:=.1;

'END';

_{DCR(/K1/):=DOR(/K1/)+YFAK;CR(/K1/):=CR(/Kl/)1.045;}

(19)

220 Cr(/K1/}:=(/<1/)+DCR(/Kl/) CCCCl22C 221 cF:7/(L!.s(F4V(/K1Í)*WGL)/LrtU)-2)**2; CCCC1Z3C. ¿22 C1 000C12'tL 223 (CGC12L 223

C:CF(S+S4?P)/;

CCOO12C

224 CR(/K1/):=R(/1/)Ct+CF+.C7+.O4V(/K1/)FNV(fK1/)*SQRT(9.81*L) CCCG1270 226 RT(/K1/):=RROSV(/Kl/)V(/K1/)»CR(/K1í) CCOCLZ8O 227 FIX(1 ,l,+,FNV(/Kl/) ;FIX(1,2,4,J(/K1/) );FIX(1,2,4,V(/K1f}/.5144) CCGC12C 230 FIX(1,l5,C,RT(/K1/));LIN(i,1); COCC13CC-232 t;ND; CCCC13IO 2

LTN(L,5) LI-(,1),A(,\'V),

(C132

CCC.C13U 239 -RE(O,VS);GL?l/SRT(9.F31»L);CUEFTH(FNV,RT,K3,C); - --COCc1340 242 c0CC135C 24

GCCC1ôC

245 UGOC137C RTvTHILP(C,F'!V,K3,VS1*WGL -cGCO138c .47 IX(i,1,4,V1*WGL);FIX(l,2,4,Vi);biX(l,2,4,VS2)FIX(1,lC,RlV000L1î90 251 );LINE(11);

cccc'+co

252 'ENo' ; CCc(141C 253

(CCC12C

254 '[ND'; cocci OL'TST I"JG(1,'('D VOLGI

1U-SLN

JOc'R

ûIFSTKCNDITIS

J'rN

AANcC_C141

256

R4DcNt));LINE(1,1); -00001432 257

OTST1NG(l,'('NO) tTLA4TISC-t-- RJut:,co:T4ARTs

i ZJ'ER EN WLCtiC143 258

INTE1-2OU)'};LINf1,1);

T T

CCCC144

259

OUTTI(1('NJJRD 41LMNTI3Cr4

JTE,ESTtAATS RES.

I ZOMER EN 3CCi4 ¿60

INTEP:2-3')t)LINC1,l)

26i OUTSTING(i,'('PACIIC ROUTE

-CCCCi47

262

:1-3c.V));LIN(1,1);

(CC14h

263 OUTSTRINC(i,'('ZUID TLANTIC EN AUST.f4LTSCHE ROUTSS

CCC143

¿54

:12H'));LIN(l,i);

CCCCI4AO 265 JJTST'JNG(1,'('00ST ÁZITISCH ROUTE - --GCC1441 :13-2C')');LiNE(1,3);

CCCC1'L

267

'OTC'EIND;

CCCC144 Sj\' CCCC144e 2(2 LFJUT:JUTSTRING(I ( EU1TEN

ÙE GNZN' )

) ;LINE(1, 1); £C001449 275

I4D::

CCCC140

275 'END' CC0C14O