Scribd
Upload
51 views9 pages

Function of Longitudinal Acceleration

The document describes equations and calculations for modeling ship dynamics, including longitudinal acceleration, wind forces, ship resistance, propeller thrust, and other factors. It includes equations with coefficients for variables like density, drag, lift, wind speed, ship dimensions, resistance components, and propeller specifications. The goal is to calculate the longitudinal acceleration of a ship based on its operating parameters and environmental conditions.

Uploaded by

Sugiati Widodo
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
51 views9 pages

Function of Longitudinal Acceleration

The document describes equations and calculations for modeling ship dynamics, including longitudinal acceleration, wind forces, ship resistance, propeller thrust, and other factors. It includes equations with coefficients for variables like density, drag, lift, wind speed, ship dimensions, resistance components, and propeller specifications. The goal is to calculate the longitudinal acceleration of a ship based on its operating parameters and environmental conditions.

Uploaded by

Sugiati Widodo
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 9

%% Function of longitudinal acceleration

function [fX1] = functionla_no_tug(u,v,r,n,L,B,d,Cb,Cp,delta,Uw,H,Attack)

%% wind

CD=0,75; % f (B/2T ) - Leite, 1998)

At=1224.8; %

Al=4682; %

Hbr=30,.81;

C=-0.1899;

Hc=11.904

Aod=411.66;

Cbr=-9.977;

%%

%% Calculation of Wind Force

airDensity=1.1644; % at 30 C

X00=-0.330;

X01=0.293;

X02=0.193;

X03= 0.682;

...

X10=-1.353;

X11=1.700;
X12=2.87;

X13=-0.463;

X14=-0.570;

X15=-6.640

X16=-0.0123;

X17=0.0202

....

X30=0.830;

X31=-0.413;

X32=-0.0827;

X33=-0.563;

X34=0.804;

X35=-5.67;

X36=0.0401;

X37=-0.132;

.....

X50=0.0372;

X51=-0.0075;

X52=-0.103;

X53=0.0921;

X0=X00+(X01*B*Hbr/At)+(X02*C/Hc)+(X03*Aod/L^2);

X1=X10+(X11*Al/(L*B) )+(X12*L*Hc/Al)+(X13*L*Hbr/Al)+(X14*Aod/Al)+(X15*At/(L*B) ) +
(X16*(At/L^2)^-1)+ (X17*(Hc/L)^-1)

X3=X30+(X31*(L*Hbr/Al)^-1)+(X32*Al/At)+(X33*L*Hc/Al)+(X34*Aod/Al)+(X35*Aod/L^2)+(X36*C/Hc)+
(X37*Cbr/L);
X5=X50+(X51*(Aod/Al)^-1)+(X52*Cbr/L)+(X53*Al/(L*B) );

Xw=0,5*airDensity*Uw.^2*(X0+X1*cos(Attack)+X3*cos(3*Attack)+X5*cos(3*Attack) );

%% OK

% calculation of ship resistance

xG=-0.0139*L; % buku Jepang

% xG=-0.0905*35;

deltaCf=(105*(125*10^-6/L)^1.3-0.64)*10^-3;

visco=1.05/10^6;

g=9.81;

rho=1025;

UU=sqrt(u^2+(v-xG*r)^2);

Fn=UU/sqrt(g*L);

slx=3.36*Fn+Cp-0.7;

c=18.32*slx^3-27*slx^2+14.2*slx-1.83;

Cr=(c+012*(B/d-3)+50*( (L*B*d*cb) /L^3-0.007) )/1000;

Rn=UU*L/visco;

logRn=log10 (Rn) ;

Cf=0.455*logRn^-2.58;

Ct=Cf+deltaCf+Cr;

WSA=(1.7*d+Cb*B) *L;

X0=-(Ct*0.5*rho*WSA*UU*UU) ;

%% Ship resistance Holltrop and Mennen Method


%% Calculation of friction resistance

% Cstern=-10; % V-shaped section

% Cstern=0; % normal section shape

Cstern=10; % U-shaped section with Hogner stern

c13=1+0.003*Cstern;

dPerL=d/L;

if d/L > 0.05

c12=dperL^0.2228466;

elseif 0.05>=dPerL && dPerL>=0.02

c12=48.20*(dPerL-0.02)^2.078+0.479948;

else

c12=0.479948;

end

lcb=1.81/100*L; % SR108, JASNAOE 176, Matsumoto N, Suemitsu K, The Prediction of Manuevering

% Cb=0.5617; %SR108, JASNAOE 176, Matsumoto N, Suemitsu K, The Prediction of Manuevering

% Cp=0.5821; %SR108, JASNAOE 176, Matsumoto N, Suemitsu K, The Prediction of Manuevering

% Cm=0.977+0.085*(cb-0.60); %Preliminary ship design parameter estimation

Vs=19.5*.5144;

FN=Vs/sqrt(9.81*L);

Cb = 1.08-1.68*Fn; % H. Scheekluth, V. Bertram, Ship Design for Efficiency

Cm=0.9+0.1*Cb; % H. Scheekluth, V. Bertram, Ship Design for Efficiency

CP=Cb/Cm;

Cwp=0.95*Cp+0.17*(1-Cp)^(1/3) ; % U section form, H. Scheekluth, V. Bertram, Ship Design

% Cwp=(1+2*Cb/Cm^0.5)/3; % average section, H. Scheekluth, V. Bertram, Ship Design


ABT=0.25*pi*(0.6*d)^2; % etimation, transverse section area of

LR = L*(1-Cp+0.06*Cp*lcb/(4*Cp-1) );

SatuPLUSk1=c13*(0.93+c12*(B/LR)^0.92497*(0.95-Cp)^-0.521448*(1-Cp+0.0225*lcb)^0.6906);

CF=0.075/(logRn-2)^2;

WSA=L*(2*d+B)*Cm^0.5*(0.453+0.4425*Cb-0.2862*Cm-0.003467*B/d+0.3696*Cwp)+2.38ABT/Cb;

RF=0.5*rho*UU^2*WSA*CF

RAPP1=0.5*rho*UU^2*Srudder*CF*(SatuPLUSk2rudderbehindskeg) ;

RAPP2=0.5*rho*UU^2*Sskeg*CF*(SatuPLUSk2skeg);

RAPP3=0.5*rho*UU^2*Sbossing*CF*(SatuPLUSk2hullbossings);

RAPP4=0.5*rho*UU^2*Sbilgakeel*CF*(SatuPLUSk2bilgakeels);

RAPP=RAPP1+RAPP2+RAPP3+RAPP4;

%% Ship resistance Holltropand Mennen Method

%% Calculation of wave resistance

vol=L*B*d*Cb;

lcb=L*(8.80-38.9*Fn)/100; % H. Scheekluth, V. Bertram, Ship Design for Efficiency and Ecomnomy,


second

iE=1+89*exp(-(L/B)^0.80856*(1-Cwp)^030484*(1-Cp-0.0225*lcb)^0.6367*(LR/B)^0.34574*(00*vol/
L^3)^0.16302);

if (L/B) <12

lamdawave=1.446*Cp-0.03*L/B;

else

lamdawave=1.446*Cp-0.36;

end
if Cp<0.80

c16=9.07981*Cp-13.8673*Cp^2+6.984388*Cp^3;

else

c16= 1.73014-0.7067*Cp;

end

m1=0.0140407*L/d-1.75254*vol^(1/3)/L-4.79323*B/L-c16;

if (L^3/vol)<512

c15=-1.69385;

else if 512<=(L^3/vol) && (L^3/vol) <=1727

c15=-1.69385+ (L/vol^(1/3)-8.0)/2.36;

else

c15=0.0;

end

m2=c15*Cp^2*exp(-0.1*Fn^-2);

AT=0.5*B*0.025*d; % immerse part of transverse area of transom part a

c5=1-0.8&*AT/(B*d*Cm) ;

dF=d; % draft at FP

hB=0.25*d; % estimation, center of ABT from kell line

c3=0.56*ABT^1.5/(B*d*(0.31*ABT^0.5+dF-hB) );

c2=exp(-1.89*sqrt(c3) );
if B/L<0.11

c7=0.229577* (B/L) ^ 0.33333;

elseif 0.11<=B/L && B/L<=0.25

c7=B/L;

else

c7=0.5-0.0625*L/B;

end

cl=2223105*c7^3.78613*(d/B)^1.07961*(90-iE)^-1.37565;

dwave=-0.9;

% RW=c1*c2*c5*vol*rho*g*exp(m1*Fn^dwave+m2*cos(lamdawave*Fn^-2) )

RW=0;

%% Ship resistance Holltrop and Mennen Method

%% Resistance of bulbous and immerese transom

FNT=UU/sqrt (2*g*AT/ (B+B*Cwp) );

if FNT<5

c6=0.2*(1-0.2*FNT);

else

c6=0.0;

end

RTR=0.5*rho*UU*UU*AT*c6;

FNI=vol/sqrt (g*(dF-hB-0.25*ABT^0.5)+0.15*UU^2);

PB=0.56*ABT^0.5/(dF-1.5*hB);

RB=0.11*exp(-3*PB^-2)*FNI^3*ABT^1.5*rho*g/(1+FNI^2);
%% Ship resistance Holltrop and Mennen Method

%% Resistance of model-ship corelation

if dF/L<=0.04

c4=dF/L;

else

c4=0.04;

end

CA= 0.006* (L+100) ^-0.16-0.00205+0.003*sqrt(L/7.5)*Cb^4*c2*(0.04-c4);

% RA=0.5*rho*UU^2*WSA*CA

RA=0;

%% Ship resistance Holltrop and Mennen Method

X0=-(RF*SatuPLUSk1+RAPP+RW+RB+RTR+RA);

% Ship resistance Holltrop and Mennen Methodot finished

%% Tahanan dari Buku Jepang

xG=-0.0139*L; % buku jepang

UU=sqrt(u^2+(v-xG*r)^2);

% buku Jepang

rho=1025;

X0D=-0.0062;

X0=x0D* (rho/2)*L*d*UU*UU;

% calculation of propeller trust based of Kt-J B5-0.73

falsekeel=0.0;
dem=d+falsekeel;

% Ar=1/45.8*L*dem;

% Dp=0.1857*35;

% buku Jepang

Ar=0.02466*L*dem;

% Dp=5.7;

Dp=0.6*d;

% global~n;

rho=1025;

% buku Jepang

J0=0.33;

J1=0.22;

J2=0.16;

wp = 1-0.780;

tp = 1-0.82;

if n==0

Xp=0;

Xr=0;

else

You might also like