COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6

6. AUTOCAD SCRIPT : HYDROSTATICS

This program reads offset file and generates script files for hydrostatics and an output file

INPUT
OFFSET FILE

HYDROSTATICS OUTPUT
SCRIPT FILE FILE

OFFSET FILE
142.000000 19.086000 6.160000 13.838712
23 no of stations
5 no of offsets
0.000000 0.000000 5.550000 station 0
0.000000 2.100000 5.650000
0.000000 3.800000 5.750000
0.000000 4.700000 6.000000
0.000000 5.066500 6.160000
12
7.101625 0.000000 4.638187 station 1
7.101625 0.425926 4.688814
7.101625 1.880353 4.902474
7.101625 2.977689 5.080120
7.101625 3.881452 5.240041
7.101625 4.643124 5.395631
7.101625 4.984827 5.487795
7.101625 5.307661 5.600055
7.101625 5.615051 5.739171
7.101625 5.902462 5.907385
7.101625 6.163370 6.105809
7.101625 6.223200 6.160000
MAIN PROGRAM
***************************************************************
***** E X A M P L E 6 . F O R *****
***************************************************************

IMPLICIT DOUBLE PRECISION (A-H,O-Z)

DIMENSION XOFS(99),YOFS(99,99),ZOFS(99,99),ZOFS1(99)
DIMENSION SC1(99,99),SM1(99,99),SC(99),SM(99),HB(99)
DIMENSION SCPRN(99,99),SMPRN(99,99),HBPRN(99,99)
DIMENSION X0(99),Y0(99),Z0(99)
DIMENSION VOL(99),DISP(99)
DIMENSION CB(99),CBA(99),CBF(99)
DIMENSION CM(99),CP(99),CPA(99),CPF(99)
DIMENSION XLCB(99),XLCF(99)
DIMENSION CWP(99),CWPF(99),CWPA(99),AWP(99)
DIMENSION CVP(99),CVPA(99),CVPF(99)
DIMENSION CVB(99),CVBA(99),CVBF(99)
DIMENSION TIL(99),TIT(99),BMT(99),BML(99)
DIMENSION TPCM(99),TMTCM(99),BKB(99)
DIMENSION BKMT(99),BKML(99),GMT(99),GML(99),WETALAN(99)

* OFFSET FILE NAME

OPEN(UNIT=1,FILE='EXAMPLE6.DAT')
* OUTPUT FILE NAME
OPEN(UNIT=7,FILE='EXAMPLE6.OUT')

6.1
 COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6
* SCRIPT FILE NAME
OPEN(UNIT=8,FILE='EXAMPLE6.SCR')

* VERTICAL CENTER OF GRAVITY [M]

CG=0.0
* DENSITY OF WATER [T/M3]
RHO=1.025
* NO OF SECTION CORRESPONDİNG TO MIDSHIPS
NMID = 11

* READ MAIN DIMENSIONS AND OFFSETS

READ(1,*) PL,B,T,D
READ(1,*) NS,NW
DO I=1,NS
DO J=1,NW
READ(1,*) XOFS(I),YOFS(I,J),ZOFS(I,J)
ENDDO
ENDDO

* SPECIFY WATERLINE HEIGHTS

DO J=1,NW
ZOFS1(J)=ZOFS(NMID,J)
ENDDO

* CALCULATE SECTIONAL AREA DISTRIBUTION

DO I=1,NS
DO M=1,NW
DO J=1,M
Y0(J) = YOFS(I,J)
X0(J) = ZOFS(I,J)
ENDDO
CALL SPLINE (X0,Y0,M,AREA,CEN, RAD)
SC1(I,M)=AREA*2.
IF(M.EQ.1) SC1(I,M)=0.0
ENDDO
ENDDO

* CALCULATE SECTIONAL MOMENTS

DO I=1,NS
DO M=1,NW
DO J=1,M
Y0(J) = YOFS(I,J)*ZOFS(I,J)
X0(J) = ZOFS(I,J)
ENDDO
CALL SPLINE (X0,Y0,M,AREA,CEN,RAD)
SM1(I,M)=AREA*2.
IF(M.EQ.1) SM1(I,M)=0.0
ENDDO
ENDDO

* CALCULATE SECTIONAL AREA, MOMENT AND HALF BREADTH AT SPECIFIED WATERLINES

DO M=1,NW

ZZZ = ZOFS1(M)
DO I=1,NS
SC(I)=SC1(I,M)
SCPRN(I,M)=SC(I)
SM(I)=SM1(I,M)
SMPRN(I,M)=SM(I)
HB(I)=YOFS(I,M)
HBPRN(I,M)=HB(I)
ENDDO

* WATERLINE BREADTH
BWL = HB(1)
DO I=2,NS
BWL=MAX(BWL,HB(I))
ENDDO
BWL=BWL*2.

* CALCULATE DISPLACEMENT
CALL SPLINE(XOFS,SC,NS,AREA3,XCEN,RAD)

VOL(M) = AREA3
DISP(M) = VOL(M) * RHO

* CB, CM, CP, LCB

XLCB(M) = ((XCEN-PL/2.)/PL)*100.
CB (M) = AREA3/(PL*BWL*ZZZ)
CM (M) = SC(NMID)/(BWL*ZZZ)

IF(CM(M). EQ.0.0) CP(M) = 0.0

IF(CM(M). NE.0.0) CP(M) = CB(M)/CM(M)

* CWP AND LCF

6.2
 COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6
CALL SPLINE (XOFS,HB,NS,AREA1,XCEN,RAD)
ALCF = XCEN
XLCF(M) = ((XCEN-PL/2.)/PL)*100.
CWP (M) = 2.*AREA1/(PL*BWL)
AWP (M) = 2.*AREA1

* CVP
IF(CWP(M).EQ.0.0) CVB(M)= 0.0
IF(CWP(M).NE.0.0) CVB(M)= CB(M)/CWP(M)
IF(CWP(M).EQ.0.0) CVP(M)= 0.0
IF(CWP(M).NE.0.0) CVP(M)= CP(M)/CWP(M)

* KB
CALL SPLINE (XOFS,SM,NS,AREA1,CEN,XRAD)
BKB(M) = AREA1/VOL(M)

* IT
DO I=1,NS
X0(I) = 2.*HB(I)**3.
ENDDO

CALL SPLINE (XOFS,X0,NS,AREA,CEN,RAD)

TIT(M) = AREA/3.

* IL
DO I=1,NS
X0(I) = HB(I)*XOFS(I)*XOFS(I)
ENDDO

CALL SPLINE (XOFS,X0,NS,AREA,CEN,RAD)

TIL(M) = 2.*AREA - AWP(M)*ALCF*ALCF

* BMT , BML
BMT(M) = (TIT(M)/VOL(M))
BML(M) = TIL(M)/VOL(M)

* KMT , KML
BKMT(M) = BKB(M) + BMT(M)
BKML(M) = BKB(M) + BML(M)

* GMT , GML
GMT(M) = BKMT(M) - CG
GML(M) = BKML(M) - CG

* TPcm , MTcm
TPCM(M) = AWP(M)*RHO/100.
TMTCM(M) = DISP(M)*GML(M)/(100.*PL)

ENDDO

* PRINT OFFSETS SECTIONAL AREAS AND SECTIONAL MOMENTS

DO I=1,NS
WRITE(7,500)I,XOFS(I)
DO J=1,NW
WRITE(7,510) ZOFS1(J),HBPRN(I,J),SCPRN(I,J),SMPRN(I,J)
ENDDO
WRITE(7,520)
ENDDO

* PRINT HYDROSTATIC CHARACTERISTICS

WRITE (7,600)
DO M=2,NW
WRITE(7,610)M,ZOFS1(M),VOL(M),DISP(M),CB(M),CM(M),CP(M),CWP(M),
* CVP(M)
ENDDO

WRITE (7,620)
DO M=2,NW
WRITE(7,630)M,ZOFS1(M),XLCB(M),XLCF(M),AWP(M)
ENDDO

WRITE (7,660)
DO M=2,NW
WRITE(7,670)M,ZOFS1(M),TIT(M),TIL(M),BKB(M),BMT(M),BML(M)
ENDDO

WRITE (7,680)
DO M=2,NW
WRITE(7,690)M,ZOFS1(M),BKMT(M),BKML(M),GMT(M),GML(M)
ENDDO

* AUTOCAD PLOTTING ROUTINE

WRITE(8,107)'LIMITS 0,0 285,198 GRID 10'
WRITE(8,106)'ZOOM A'
WRITE(8,109)'PLINE'

6.3
 COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6
WRITE(8,103) 10.,10.
WRITE(8,103) 280.,10.
WRITE(8,103) 280.,190.
WRITE(8,103) 10.,190.
WRITE(8,103) 10.,10.

AX = 15.
AY = 15.
UX = 275.
UY = 185.

WRITE(8,*)'PLINE'
WRITE(8,103)AX,AY
WRITE(8,103)UX,AY
WRITE(8,103)UX,UY
WRITE(8,103)AX,UY
WRITE(8,103)AX,AY

* AG CIZGILERI - LTSCALE = 10 LTYPE = DOT

WRITE(8,*)'LTSCALE'
WRITE(8,108)4.
WRITE(8,109)'LAYER'
WRITE(8,111)'M'
WRITE(8,111)'1'
WRITE(8,111)'L'
WRITE(8,114)'DOT'
WRITE(8,111)'1'

SCALE1= (UX-AX)/PL
WRITE(8,*)'COLOR'
WRITE(8,106)'YELLOW'
DO I=1,NS
WRITE(8,109)'PLINE'
X1 = XOFS(I)*SCALE1+AX
WRITE(8,103) X1,AY
WRITE(8,103) X1,UY
WRITE(8,103)
ENDDO

SCALE2= (UY-AY)/ZOFS1(NW)

DO J=1,NW
WRITE(8,109)'PLINE'
Y1 = ZOFS1(J)*SCALE2+AY
WRITE(8,103) AX,Y1
WRITE(8,103) UX,Y1
WRITE(8,103)
ENDDO

* PLOT BONJEAN SECTIONAL AREA CURVES

WRITE(8,109)'LAYER'
WRITE(8,111)'S'
WRITE(8,111)'0'

SCMAX = 0.0
DO I=1,NS
DO J=1,NW
SCMAX=MAX(SCMAX,SCPRN(I,J))
ENDDO
ENDDO

SCALE3=1.75*((UX-AX)/FLOAT(NS-1))/SCMAX
WRITE(8,*)'COLOR'
WRITE(8,109)'WHITE'
DO I=1,NS
WRITE(8,109)'PLINE'
DO J=1,NW
Z1=ZOFS1(J)*SCALE2+AY
Y1=SCPRN(I,J)*SCALE3+AX+XOFS(I)*SCALE1
WRITE(8,103)Y1,Z1
ENDDO
WRITE(8,*)'PEDIT'
WRITE(8,111)'L'
WRITE(8,111)'F'
WRITE(8,103)
ENDDO

* LABEL BONJEAN SECTIONAL AREA CURVES

WRITE(8,109)'STYLE'
WRITE(8,104)'COMPLEX'
WRITE(8,104)'COMPLEX'
WRITE(8,108)1.50
WRITE(8,105)1
WRITE(8,105)0
WRITE(8,111)'N'

6.4
 COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6
WRITE(8,111)'N'
WRITE(8,111)'N'

NW23 = NW*2/3

WRITE(8,112)'TEXT'
WRITE(8,111)'S'
WRITE(8,104)'COMPLEX'
X1=(SCPRN(NMID,NW23)*SCALE3+AX+XOFS(NMID)*SCALE1)*0.99
Y1=ZOFS1(NW23)*SCALE2+AY
X2=(SCPRN(NMID,NW)*SCALE3+AX+XOFS(NMID)*SCALE1)*0.99
Y2=ZOFS1(NW)*SCALE2+AY
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
XSCALE1=SCMAX/(1.75*26./FLOAT(NS-1))
WRITE(8,113)'SECTIONAL AREA 1CM=',XSCALE1,'M^2'

* PLOT BONJEAN SECTIONAL MOMENT CURVES

WRITE(8,109)'LAYER'
WRITE(8,111)'M'
WRITE(8,111)'2'
WRITE(8,111)'L'
WRITE(8,106)'DASHED'
WRITE(8,111)'2'

SCMAX = 0.0
DO I=1,NS
DO J=1,NW
SCMAX=MAX(SCMAX,SMPRN(I,J))
ENDDO
ENDDO

SCALE3=1.25*((UX-AX)/FLOAT(NS-1))/SCMAX

DO I=1,NS
WRITE(8,*)'PLINE'
DO J=1,NW
Z1=ZOFS1(J)*SCALE2+AY
Y1=SMPRN(I,J)*SCALE3+AX+XOFS(I)*SCALE1
WRITE(8,103)Y1,Z1
ENDDO
WRITE(8,*)'PEDIT'
WRITE(8,111)'L'
WRITE(8,111)'F'
ENDDO

* LABEL BONJEAN SECTIONAL MOMENT CURVES

WRITE(8,*)'TEXT'
WRITE(8,111)'S'
WRITE(8,104)'COMPLEX'
X1=(SMPRN(NMID,NW23)*SCALE3+AX+XOFS(NMID)*SCALE1)*0.99
Y1=ZOFS1(NW23)*SCALE2+AY
X2=(SMPRN(NMID,NW)*SCALE3+AX+XOFS(NMID)*SCALE1)*0.99
Y2=ZOFS1(NW)*SCALE2+AY
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
XSCALE2=SCMAX/(1.25*26./FLOAT(NS-1))
WRITE(8,115)'SECTIONAL MOMENT 1CM=',XSCALE2,'M^3'

WRITE(8,109)'LAYER'
WRITE(8,111)'S'
WRITE(8,111)'0'

* PLOT VOLUME AND DISPLACEMENT CURVES

SCALE4=0.95*(UX-AX)/VOL(NW)
WRITE(8,*)'PLINE'
DO M=1,NW
Z1=ZOFS1(M)*SCALE2+AY
Y1=VOL(M)*SCALE4+AX
WRITE(8,103)Y1,Z1
ENDDO

WRITE(8,*)'PEDIT'
WRITE(8,103)Y1,Z1
WRITE(8,111)'F'

WRITE(8,*)'PLINE'
DO M=1,NW
Z1=ZOFS1(M)*SCALE2+AY
Y1=DISP(M)*SCALE4+AX
WRITE(8,103)Y1,Z1
ENDDO

WRITE(8,*)'PEDIT'
WRITE(8,103)Y1,Z1

6.5
 COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6
WRITE(8,111)'F'

* LABEL VOLUME AND DISPLACEMENT CURVES

NWA = NW*5/6
NWB = NW

WRITE(8,*)'TEXT'
X1=VOL(NWA)*SCALE4+AX
Y1=(ZOFS1(NWA)*SCALE2+AY)*1.01
X2=VOL(NWB)*SCALE4+AX
Y2=(ZOFS1(NWB)*SCALE2+AY)*1.01
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
XSCALE3=VOL(NW)/(0.95*26.)
WRITE(8,116)'VOLUME 1CM=',XSCALE3,'M^3'

WRITE(8,112)'TEXT'
X1=DISP(NWA)*SCALE4+AX
Y1=(ZOFS1(NWA)*SCALE2+AY)*0.98
X2=DISP(NWB)*SCALE4+AX
Y2=(ZOFS1(NWB)*SCALE2+AY)*0.98
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
XSCALE4=XSCALE3
WRITE(8,117)'DISPLACEMENT 1CM=',XSCALE4,'T'

* PLOT AWP AND TPCM CURVES

SCALE5=0.7*(UX-AX)/AWP(NW)
WRITE(8,109)'PLINE'
DO M=1,NW
Z1=ZOFS1(M)*SCALE2+AY
Y1=AWP(M)*SCALE5+AX
WRITE(8,103)Y1,Z1
ENDDO

WRITE(8,*)'PEDIT'
WRITE(8,103)Y1,Z1
WRITE(8,111)'F'

WRITE(8,*)'PLINE'
DO M=1,NW
Z1=ZOFS1(M)*SCALE2+AY
Y1=100.*TPCM(M)*SCALE5+AX
WRITE(8,103)Y1,Z1
ENDDO

WRITE(8,*)'PEDIT'
WRITE(8,103)Y1,Z1
WRITE(8,111)'F'

* LABEL AWP AND TPCM CURVES

NWA = NW/2-2
NWB = NW/2

WRITE(8,*)'TEXT'
X1=AWP(NWA)*SCALE5+AX
Y1=(ZOFS1(NWA)*SCALE2+AY)*1.015
X2=AWP(NWB)*SCALE5+AX
Y2=(ZOFS1(NWB)*SCALE2+AY)*1.015
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
XSCALE5=AWP(NW)/(0.7*26.)
WRITE(8,118)'AWP 1CM=',XSCALE5,'M^2'

WRITE(8,112)'TEXT'
X1=(100.*TPCM(NWA)*SCALE5+AX)*1.01
Y1=(ZOFS1(NWA)*SCALE2+AY)*0.95
X2=(100.*TPCM(NWB)*SCALE5+AX)*1.01
Y2=(ZOFS1(NWB)*SCALE2+AY)*0.95
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
XSCALE6=TPCM(NW)/(0.7*26.)
WRITE(8,119)'TPCM 1CM=',XSCALE6,'T/CM'

* PLOT KB AND BMT CURVES

BKMTMAX=BKMT(2)
DO I=2,NW
BKMTMAX=MAX(BKMT(I),BKMTMAX)
ENDDO

SCALE7=0.8*(UX-AX)/BKMTMAX
WRITE(8,109)'PLINE'
DO M=2,NW
Z1=ZOFS1(M)*SCALE2+AY
Y1=BKB(M)*SCALE7+AX

6.6
 COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6
WRITE(8,103)Y1,Z1
ENDDO

WRITE(8,*)'PEDIT'
WRITE(8,103)Y1,Z1
WRITE(8,111)'F'

WRITE(8,*)'PLINE'
DO M=2,NW
Z1=ZOFS1(M)*SCALE2+AY
Y1=BKMT(M)*SCALE7+AX
WRITE(8,103)Y1,Z1
ENDDO

WRITE(8,*)'PEDIT'
WRITE(8,103)Y1,Z1
WRITE(8,111)'F'

* LABEL KB AND BMT CURVES

NWA = NW*3/4
NWB = NW

WRITE(8,*)'TEXT'
X1=BKB(NWA)*SCALE7+AX
Y1=(ZOFS1(NWA)*SCALE2+AY)*1.01
X2=BKB(NWB)*SCALE7+AX
Y2=(ZOFS1(NWB)*SCALE2+AY)*1.01
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
XSCALE7=BKMTMAX/(0.8*26.)
WRITE(8,121)'KB 1CM=',XSCALE7,'M'

NWA=2
NWB=3
WRITE(8,112)'TEXT'
X1=BKMT(NWA)*SCALE7+AX
Y1=(ZOFS1(NWA)*SCALE2+AY)*1.01
X2=BKMT(NWB)*SCALE7+AX
Y2=(ZOFS1(NWB)*SCALE2+AY)*1.01
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
XSCALE8=XSCALE7
WRITE(8,122)'TRANS. BM 1CM=',XSCALE8,'M'

* PLOT LCB AND LCF CURVES

SCALE9=0.1*(UX-AX)/5.
WRITE(8,109)'PLINE'
DO M=2,NW
Z1=ZOFS1(M)*SCALE2+AY
Y1=XLCB(M)*SCALE9+AX+(UX-AX)/2.
WRITE(8,103)Y1,Z1
ENDDO

WRITE(8,*)'PEDIT'
WRITE(8,103)Y1,Z1
WRITE(8,111)'F'

WRITE(8,*)'PLINE'
DO M=2,NW
Z1=ZOFS1(M)*SCALE2+AY
Y1=XLCF(M)*SCALE9+AX+(UX-AX)/2.
WRITE(8,103)Y1,Z1
ENDDO

WRITE(8,*)'PEDIT'
WRITE(8,103)Y1,Z1
WRITE(8,111)'F'

* LABEL LCB AND LCF CURVES

NWA = NW*4/5
NWB = NW

WRITE(8,*)'TEXT'
X1=(XLCB(NWA)*SCALE9+AX+(UX-AX)/2.)*1.02
Y1=ZOFS1(NWA)*SCALE2+AY
X2=(XLCB(NWB)*SCALE9+AX+(UX-AX)/2.)*1.02
Y2=ZOFS1(NWB)*SCALE2+AY
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
XSCALE9=5./2.6
WRITE(8,123)'LCB 1CM=',XSCALE9

WRITE(8,112)'TEXT'
X1=(XLCF(NWA)*SCALE9+AX+(UX-AX)/2.)*0.99
Y1=ZOFS1(NWA)*SCALE2+AY

6.7
 COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6
X2=(XLCF(NWB)*SCALE9+AX+(UX-AX)/2.)*0.99
Y2=ZOFS1(NWB)*SCALE2+AY
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
XSCALE10=XSCALE9
WRITE(8,124)'LCF 1CM=',XSCALE10

* PLOT BML CURVE

BMLMAX=BML(2)
DO I=2,NW
BMLMAX=MAX(BML(I),BMLMAX)
ENDDO

SCALE10=0.9*(UX-AX)/BMLMAX
WRITE(8,109)'PLINE'
DO M=2,NW
Z1=ZOFS1(M)*SCALE2+AY
Y1=BML(M)*SCALE10+AX
WRITE(8,103)Y1,Z1
ENDDO

WRITE(8,*)'PEDIT'
WRITE(8,103)Y1,Z1
WRITE(8,111)'F'

* LABEL BML CURVE

NWA = NW*3/4
NWB = NW

WRITE(8,*)'TEXT'
X1=(BML(NWA)*SCALE10+AX)*0.98
Y1=(ZOFS1(NWA)*SCALE2+AY)
X2=(BML(NWB)*SCALE10+AX)*0.98
Y2=(ZOFS1(NWB)*SCALE2+AY)
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
XSCALE11=BMLMAX/(0.9*26.)
WRITE(8,125)'BML 1CM=',XSCALE11,'M'

* PLOT CB,CP AND CWP CURVES

SCALE11=0.2*(UX-AX)
WRITE(8,109)'PLINE'
DO M=2,NW
Z1=ZOFS1(M)*SCALE2+AY
Y1=CB(M)*SCALE11+AX
WRITE(8,103)Y1,Z1
ENDDO

WRITE(8,*)'PEDIT'
WRITE(8,103)Y1,Z1
WRITE(8,111)'F'

WRITE(8,*)'PLINE'
DO M=2,NW
Z1=ZOFS1(M)*SCALE2+AY
Y1=CP(M)*SCALE11+AX
WRITE(8,103)Y1,Z1
ENDDO

WRITE(8,*)'PEDIT'
WRITE(8,103)Y1,Z1
WRITE(8,111)'F'

WRITE(8,*)'PLINE'
DO M=2,NW
Z1=ZOFS1(M)*SCALE2+AY
Y1=CWP(M)*SCALE11+AX
WRITE(8,103)Y1,Z1
ENDDO

WRITE(8,*)'PEDIT'
WRITE(8,103)Y1,Z1
WRITE(8,111)'F'

* LABEL CB, CP AND CWP CURVES

NWA = NW*5/6
NWB = NW

WRITE(8,*)'TEXT'
X1=(CB(NWA)*SCALE11+AX)*0.99
Y1=ZOFS1(NWA)*SCALE2+AY
X2=(CB(NWB)*SCALE11+AX)*0.99
Y2=ZOFS1(NWB)*SCALE2+AY
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2

6.8
 COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6
XSCALE12=1./(0.2*26.)
WRITE(8,126)'CB 1CM=',XSCALE12

WRITE(8,112)'TEXT'
X1=(CP(NWA)*SCALE11+AX)*1.06
Y1=ZOFS1(NWA)*SCALE2+AY
X2=(CP(NWB)*SCALE11+AX)*1.06
Y2=ZOFS1(NWB)*SCALE2+AY
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
WRITE(8,126)'CP 1CM=',XSCALE12

NWA=NW/2

WRITE(8,112)'TEXT'
X1=(CWP(NWA)*SCALE11+AX)*1.05
Y1=ZOFS1(NWA)*SCALE2+AY
X2=(CWP(NWB)*SCALE11+AX)*1.05
Y2=ZOFS1(NWB)*SCALE2+AY
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
WRITE(8,127)'CWP 1CM=',XSCALE12

* PLOT TMTCM CURVE

SCALE13=(UX-AX)*0.95/TMTCM(NW)
WRITE(8,109)'PLINE'
DO M=2,NW
Z1=ZOFS1(M)*SCALE2+AY
Y1=TMTCM(M)*SCALE13+AX
WRITE(8,103)Y1,Z1
ENDDO

WRITE(8,*)'PEDIT'
WRITE(8,103)Y1,Z1
WRITE(8,111)'F'

* LABEL TMTCM CURVE

NWA = NW*3/5
NWB = NW

WRITE(8,*)'TEXT'
X1=TMTCM(NWA)*SCALE13+AX
Y1=(ZOFS1(NWA)*SCALE2+AY)*1.01
X2=TMTCM(NWB)*SCALE13+AX
Y2=(ZOFS1(NWB)*SCALE2+AY)*1.01
WRITE(8,103)X1,Y1
WRITE(8,103)X2,Y2
XSCALE13=TMTCM(NW)/(0.95*26.)
WRITE(8,128)'MOMENT TO TRIM BY 1CM=',XSCALE13,'TCM'

* LABEL SHIP NAME

X1=(UX-AX)/2.
Y1=AY+(ZOFS1(NW)*SCALE2)*1.01
WRITE(8,112)'TEXT'
WRITE(8,103)X1,Y1
WRITE(8,108)0.0
WRITE(8,107)SNAME

Y1=11.0
WRITE(8,112)'TEXT'
WRITE(8,103)X1,Y1
WRITE(8,108)0.0
WRITE(8,*)'HYDROSTATICS
* K.S'

103 FORMAT(E10.5,',',E10.5)
104 FORMAT(A7)
105 FORMAT(I1)
106 FORMAT(A6)
107 FORMAT(A26)
108 FORMAT(E10.5)
109 FORMAT(A5)
111 FORMAT(A1)
112 FORMAT(A4)
113 FORMAT(A19,F9.3,A4)
114 FORMAT(A3)
115 FORMAT(A21,F9.3,A4)
116 FORMAT(A11,F9.3,A4)
117 FORMAT(A17,F9.3,A2)
118 FORMAT(A8,F9.3,A4)
119 FORMAT(A9,F9.3,A5)
121 FORMAT(A7,F9.3,A2)
122 FORMAT(A14,F9.3,A2)
123 FORMAT(A8,F9.3)
124 FORMAT(A8,F9.3)

6.9
 COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6
125 FORMAT(A8,F9.3,A2)
126 FORMAT(A8,F9.3)
127 FORMAT(A9,F9.3)
128 FORMAT(A22,F9.3,A4)

500 FORMAT(//1X,'KESIT NO:',I2,7X,'X =',F8.3/

* 1X,'--------------------------------------------------'/
* 1X,'WL [M]',5X,'OFSET[M]',5X,'ALAN[M^2]',5X,'MOM[M3]',1X,/
* 1X,'--------------------------------------------------')
510 FORMAT(F7.4,F13.4,F13.4,F13.4)
520 FORMAT(1X, '--------------------------------------------------')

600 FORMAT(//'========================================================
*================='/
* ' WL DRFT[M] VLM[M3] DISP[T] CB CM CP
* CWP CVP '/
* '========================================================
*=================')
610 FORMAT(I4,F8.4,1X,2F9.2,5F8.4)
620 FORMAT(//'========================================================
*================='/
* ' WL DRFT[M] LCB [+FWD] LCF [+FWD] AWP [M2
*] '/
* ' [%L] [%L] '/
* '========================================================
*=================')
630 FORMAT(I4,F8.4,5X,2F12.4,F17.4)
640 FORMAT(//'========================================================
*================='/
* ' WL DRFT[M] CBA CPA CWPA CVPA CBF CPF
* CWPF CVPF '/
* '========================================================
*=================')
650 FORMAT(I4,2X,9F7.4)
660 FORMAT(//'========================================================
*================='/
* ' WL DRFT[M] IT[M4] IL[M4] KB[M] BMT[
*M] BML[M] '/
* '========================================================
*=================')
670 FORMAT(I4,F8.4,F14.4,F15.4,2F8.4,F11.4)
680 FORMAT(//'========================================================
*================='/
* ' WL DRFT[M] KMT[M] KML[M] GMT[M] GML[M]'/
* '========================================================
*=================')
690 FORMAT(I4,4X,2F8.4,F10.4,2X,F8.4,F12.4)

STOP
END

******************************************
***** S U B R O U T I N E T R I D *****
******************************************

SUBROUTINE TRID(A1,A2,A3,Z,N)

IMPLICIT DOUBLE PRECISION (A-H,O-Z)

DIMENSION A1(399),A2(399),A3(399),Z(399)

DO I=2,N
A1(I)=A1(I)/A2(I-1)
A2(I)=A2(I)-A1(I)*A3(I-1)
Z(I)=Z(I)-A1(I)*Z(I-1)
ENDDO

Z(N)=Z(N)/A2(N)
NM1=N-1

DO I=NM1,1,-1
Z(I)=(Z(I)-A3(I)*Z(I+1))/A2(I)
ENDDO

RETURN
END

******************************************
***** S U B R O U T I N E SPLINE *****
******************************************

SUBROUTINE SPLINE (T,Y,N,AREA,XCEN,RAD)

IMPLICIT DOUBLE PRECISION (A-H,O-Z)

6.10
 COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6
DIMENSION T(399),Y(399),D(399),Z(399),C(399)
DIMENSION A1(399),A2(399),A3(399)
DIMENSION AA(399),BB(399),CC(399),DD(399),S(399)

IF(N.LT.2) AREA = 0.0

IF(N.LT.2) XCEN = 0.0
IF(N.LT.2) RETURN
IF(N.EQ.2) AREA = (T(2)-T(1))*(Y(2)+Y(1))/2.
IF(N.EQ.2) XCEN = (T(2)+T(1))/2.
IF(N.EQ.2) RETURN

DO I=1,N
A1(I)=0.0
A2(I)=1.0
A3(I)=0.0
Z(I)=0.0
ENDDO

D(1)=1.0
C(1)=0.0
Z(1)=0.0
NM1=N-1

DO I=2,NM1
D(I)=2.0*(T(I+1)-T(I-1))
C(I)=T(I+1)-T(I)
TEMP=(Y(I+1)-Y(I))/(T(I+1)-T(I))
Z(I)=6.*(TEMP-(Y(I)-Y(I-1))/(T(I)-T(I-1)))
ENDDO

D(N)=1.0
C(N)=0.0
Z(N)=0.0

DO I=1,N
A2(I)=D(I)
ENDDO

DO I=2,NM1
A1(I)=C(I-1)
ENDDO

DO I=2,NM1
A3(I)=C(I)
ENDDO

CALL TRID(A1,A2,A3,Z,N)

AREA = 0.0
AMOM = 0.0

DO I=2,NM1
AA(I)=(Z(I+1)-Z(I))/(6.*C(I))
BB(I)=(-T(I)*Z(I+1)+T(I+1)*Z(I))/(2.*C(I))
CC(I)=(-3.*Z(I)*T(I+1)**2+Z(I)*C(I)**2+3.*Z(I+1)*T(I)**2-
%6.*Y(I)+6.*Y(I+1)-Z(I+1)*C(I)**2)/(6.*C(I))
DD(I)=-(T(I)**3*Z(I+1))/(6.*C(I))+(T(I+1)**3*Z(I))/(6.*C(I))-
%(Z(I)*C(I)*T(I+1))/6.+(Z(I+1)*C(I)*T(I))/6.+(Y(I)*T(I+1))/C(I)
%-(Y(I+1)*T(I))/C(I)
S(I)=AA(I)*T(I)**3+BB(I)*T(I)**2+CC(I)*T(I)+DD(I)
ENDDO

AREA=AA(2)*(T(2)**4/4.-T(1)**4/4.)
% +BB(2)*(T(2)**3/3.-T(1)**3/3.)
% +CC(2)*(T(2)**2/2.-T(1)**2/2.)
% +DD(2)*(T(2)-T(1))

DO I=2,NM1
AREA=AREA+AA(I)*(T(I+1)**4/4.-T(I)**4/4.)
% +BB(I)*(T(I+1)**3/3.-T(I)**3/3.)
% +CC(I)*(T(I+1)**2/2.-T(I)**2/2.)
% +DD(I)*(T(I+1)-T(I))
ENDDO
IF(AREA.LT.0.0) AREA = 0.0

AMOM=AA(2)*(T(2)**5/5.-T(1)**5/5.)
% +BB(2)*(T(2)**4/4.-T(1)**4/4.)
% +CC(2)*(T(2)**3/3.-T(1)**3/3.)
% +DD(2)*(T(2)**2/2.-T(1)**2/2.)
IF(AMOM.LT.0.0) AMOM = 0.0

DO I=2,NM1
AMOM=AMOM+AA(I)*(T(I+1)**5/5.-T(I)**5/5.)
% +BB(I)*(T(I+1)**4/4.-T(I)**4/4.)
% +CC(I)*(T(I+1)**3/3.-T(I)**3/3.)

6.11
 COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6
% +DD(I)*(T(I+1)**2/2.-T(I)**2/2.)
ENDDO
IF(AMOM.LT.0.0) AMOM = 0.0

IF (AREA.EQ.0.0) XCEN=0.0
IF (AREA.EQ.0.0) RAD =0.0
IF (AREA.EQ.0.0) RETURN

XCEN=AMOM/AREA

ARAD=AA(2)*(T(2)**6/6.-T(1)**6/6.)
% +BB(2)*(T(2)**5/5.-T(1)**5/5.)
% +CC(2)*(T(2)**4/4.-T(1)**4/4.)
% +DD(2)*(T(2)**3/3.-T(1)**3/3.)
IF(ARAD.LT.0.0) ARAD = 0.0

DO I=2,NM1
ARAD=ARAD+AA(I)*(T(I+1)**6/6.-T(I)**6/6.)
% +BB(I)*(T(I+1)**5/5.-T(I)**5/5.)
% +CC(I)*(T(I+1)**4/4.-T(I)**4/4.)
% +DD(I)*(T(I+1)**3/3.-T(I)**3/3.)
ENDDO
IF(ARAD.LT.0.0) ARAD = 0.0

RAD=ARAD/AREA

RETURN
END

Example 6.dat
16.0000 1.6000 1.0000 1.0000
21 20
0.0000 0.0000 0.0000 1
0.0000 0.0000 0.0526
0.0000 0.0000 0.1053
0.0000 0.0000 0.1579
0.0000 0.0000 0.2105
0.0000 0.0000 0.2632
0.0000 0.0000 0.3158
0.0000 0.0000 0.3684
0.0000 0.0000 0.4211
0.0000 0.0000 0.4737
0.0000 0.0000 0.5263
0.0000 0.0000 0.5789
0.0000 0.0000 0.6316
0.0000 0.0000 0.6842
0.0000 0.0000 0.7368
0.0000 0.0000 0.7895
0.0000 0.0000 0.8421
0.0000 0.0000 0.8947
0.0000 0.0000 0.9474
0.0000 0.0000 1.0000
0.8000 0.0000 0.0000 2
0.8000 0.0156 0.0526
0.8000 0.0303 0.1053
0.8000 0.0442 0.1579
0.8000 0.0573 0.2105
0.8000 0.0695 0.2632
0.8000 0.0808 0.3158
0.8000 0.0914 0.3684
0.8000 0.1011 0.4211
0.8000 0.1099 0.4737
0.8000 0.1179 0.5263
0.8000 0.1251 0.5789
0.8000 0.1314 0.6316
0.8000 0.1368 0.6842
0.8000 0.1415 0.7368
0.8000 0.1453 0.7895
0.8000 0.1482 0.8421
0.8000 0.1503 0.8947
0.8000 0.1516 0.9474
0.8000 0.1520 1.0000
1.6000 0.0000 0.0000 3
1.6000 0.0295 0.0526
1.6000 0.0574 0.1053
1.6000 0.0838 0.1579
1.6000 0.1085 0.2105
1.6000 0.1316 0.2632
1.6000 0.1532 0.3158
1.6000 0.1731 0.3684
1.6000 0.1915 0.4211

6.12
 COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6
1.6000 0.2082 0.4737
1.6000 0.2234 0.5263
1.6000 0.2369 0.5789
1.6000 0.2498 0.6316
1.6000 0.2593 0.6842
1.6000 0.2681 0.7368
1.6000 0.2752 0.7895
1.6000 0.2808 0.8421
1.6000 0.2848 0.8947
1.6000 0.2872 0.9474
1.6000 0.2880 1.0000
2.4000 0.0000 0.0000 4
2.4000 0.0418 0.0526
2.4000 0.0814 0.1053
2.4000 0.1187 0.1579
2.4000 0.1537 0.2105
2.4000 0.1865 0.2632
2.4000 0.2170 0.3158
2.4000 0.2453 0.3684
2.4000 0.2712 0.4211
2.4000 0.2950 0.4737
2.4000 0.3165 0.5263
2.4000 0.3357 0.5789
2.4000 0.3526 0.6316
2.4000 0.3673 0.6842
2.4000 0.3797 0.7368
2.4000 0.3899 0.7895
2.4000 0.3978 0.8421
2.4000 0.4035 0.8947
2.4000 0.4069 0.9474
2.4000 0.4080 1.0000
3.2000 0.0000 0.0000 5
3.2000 0.0525 0.0526
3.2000

Example6.scr
LIMITS 0,0 285,198 GRID 10
ZOOM A
PLINE
.10000E+02,.10000E+02
.28000E+03,.10000E+02
.28000E+03,.19000E+03
.10000E+02,.19000E+03
.10000E+02,.10000E+02
PLINE
.15000E+02,.15000E+02
.27500E+03,.15000E+02
.27500E+03,.18500E+03
.15000E+02,.18500E+03
.15000E+02,.15000E+02
LTSCALE
.40000E+01
LAYER
M
1
L
DOT
1
COLOR
YELLOW
PLINE
.15000E+02,.15000E+02
.15000E+02,.18500E+03

PLINE
.28000E+02,.15000E+02
.28000E+02,.18500E+03

6.13
COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6

6.14
 COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6

Homework Assignment No 6
Prepare a script file for hydrostatic curves of a fishing vessel with following nondimensional offsets.

Sta WL WL WL WL WL WL WL WL WL WL WL
0.0 0.1 0.2 0.3 0.4 0.6 0.8 1.0 1.2 1.4 1.6
0.5 0.000 0.022 0.028 0.032 0.050 0.060 0.187 0.505 0.741 0.860 0.913
0.75 0.000 0.040 0.059 0.071 0.093 0.137 0.298 0.592 0.790 0.889 0.933
1.0 0.000 0.060 0.092 0.117 0.143 0.222 0.400 0.663 0.835 0.915 0.951
1.5 0.000 0.113 0.173 0.225 0.272 0.406 0.592 0.786 0.901 0.951 0.976
2.0 0.000 0.178 0.276 0.351 0.422 0.575 0.744 0.875 0.943 0.975 0.990
2.5 0.000 0.256 0.395 0.487 0.568 0.716 0.852 0.933 0.974 0.991 0.998
3.0 0.000 0.348 0.509 0.612 0.694 0.823 0.919 0.971 0.992 0.999 1.000
4.0 0.000 0.534 0.725 0.821 0.886 0.956 0.988 1.000 1.000 1.000 1.000
5.0 0.000 0.658 0.853 0.930 0.970 0.998 1.000 1.000 1.000 1.000 1.000
6.0 0.000 0.576 0.766 0.847 0.890 0.940 0.966 0.979 0.988 0.995 1.000
7.0 0.000 0.362 0.499 0.587 0.651 0.738 0.802 0.854 0.900 0.938 0.970
7.5 0.000 0.259 0.374 0.452 0.520 0.611 0.680 0.750 0.815 0.871 0.926
8.0 0.000 0.173 0.262 0.326 0.381 0.467 0.546 0.624 0.700 0.776 0.852
8.5 0.000 0.105 0.169 0.216 0.258 0.332 0.407 0.480 0.562 0.650 0.742
9.0 0.000 0.049 0.089 0.123 0.152 0.205 0.262 0.326 0.401 0.486 0.583
9.25 0.000 0.025 0.054 0.080 0.104 0.146 0.192 0.248 0.318 0.393 0.490

Main dimensions of the vessel are given as follows:

No Name Length Breadth Draught Curve

1 080190017 Ata Emir Erdoğdu 20 6 3 KB - BMT
2 080200003 Ayşegül Özşen 22 6 3 KB - BML
3 080200019 Cemal Yasin Atasoy 24 6 3 CB - CP
4 080200020 Alper Toprak 26 7 3 CP - CM
5 080200034 Yiğitcan Tilaver 28 7 3 AWP - CWP
6 080200040 Sefa Ramazan Nacakcı 30 7 3 LCB - LCF
7 080200052 İlker Aksoy 32 8 4 Volume - LCF
8 080200059 Emre Berber 34 8 4 Volume - Diplacement
9 080200065 M. Mustafa Karakaya 36 8 4 Displacement - CB
10 080200070 Ekin Saraçoğlu 38 9 4 KB - BML
11 080200071 Burak General 40 9 4 BMT - LCB
12 080210057 Yağmur Tuncay 42 9 4 BML - LCF
13 080210069 Emre Altun 44 10 5 AWP - CWP
14 080220026 İsranur Bülbül 46 10 5 BMT - BML
15 080220713 Ceylan Şavluk 48 10 5 Volume - CWP
16 50 11 5 Displacement - LCB
17 52 11 5 CWP - LCB
18 54 11 5 BMT – KMT
19 56 12 6 BML – KML
20 58 12 6 IT – IL
21 60 12 6 AWP – TPcm
22 62 12 6 CB – KMT
23 64 12 6 Displacement – AWP
24 66 12 6 LCB - IT

Plot Bonjean curves and ONLY the specified hydrostatic curves. Attach your code, Autocad drawing and
the script file.

6.15