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='[Link]')
* OUTPUT FILE NAME
OPEN(UNIT=7,FILE='[Link]')
6.1
COMPUTER AIDED SHIP DESIGN LECTURE NOTES – CHAPTER 6
* SCRIPT FILE NAME
OPEN(UNIT=8,FILE='[Link]')
* 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([Link].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([Link].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([Link].2) AREA = 0.0
IF([Link].2) XCEN = 0.0
IF([Link].2) RETURN
IF([Link].2) AREA = (T(2)-T(1))*(Y(2)+Y(1))/2.
IF([Link].2) XCEN = (T(2)+T(1))/2.
IF([Link].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([Link].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([Link].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([Link].0.0) AMOM = 0.0
IF ([Link].0.0) XCEN=0.0
IF ([Link].0.0) RAD =0.0
IF ([Link].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([Link].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([Link].0.0) ARAD = 0.0
RAD=ARAD/AREA
RETURN
END
Example [Link]
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
[Link]
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