Modelling reasonable operation regimes of the main propulsion plant main diesel engine - propeller - hull on the general cargo ship

Modelling reasonable operation regimes of the main propulsion plant main  
diesel engine - propeller - hull on the general cargo ship  
Prof. DSc. Do Duc Luu, Prof. Dr. Luong Cong Nho, Dr. Tran Ngoc Tu  
Vietnam Maritime University  
luudoduc@gmail.com; luongcongnho@vimaru.edu.com;tutn.dt@vimaru.edu.vn  
Abstract: The article presents the modeling results of the characteristics of main propulsion plant Diesel  
Engine Propeller Hull, based on the combination of the analytical method and regressive analysis of  
the empirical data. The authors create the algorithms on the base of the created models to select the  
reasonable operating regime of the main propulsion plant (MPP) according to the different operating  
conditions (in still water, in sea wave, in the canal, etc.). This paper provides the numerical simulation  
results of the drawn combination characteristics of MPP on the 34000 DWT general cargo ship built in  
Vietnam with the main engine MAN 6S 46 MCC, and of the selected regime respectively the desired  
operated conditions. The research results would be used for the simulation of monitoring and controlling  
the MPP on the general cargo ship, in order to improve the education quality of the electrical and marine  
engineering field at Vietnam Maritime University.  
Keywords: Main Engine, Propeller, Hull resistance, Main Propulsion Plant.  
1. Basis of selecting reasonable exploitation regime  
Choosing reasonable working regime of the combination of hull propeller - main engine under  
different exploitation regimes ensures the safety of diesel main engine (DME), so that mechanical and  
thermal overload may not be occurred. In addition, the proposing voyage efficiency of (time, economy)  
for a cargo ship and sea-going cargo fleet can be achieved by choosing reasonable exploitation regime  
of engine’s revolutions and ship’s speeds in different environmental conditions, in which MPP works.  
Further more, the efficiency of reasonably exploiting working regime of MPP is not only suitable for  
ship’s working purpose from the concepts of safety, high reliability and economy, but also ensures  
environmental norms set by International Association of Classification Societies (IACS).  
The frequent regimes in MPP’s exploitation of a marine vessel: working in still water, carrying cargo  
(Cargo Load Index, CLI) from ballast (CLI = 0) to full load (CLI = 100%); working in sea wave and  
wind with given CLI; working in canal; vessel in shallow water. In working process, the resistance states  
of ship hull and propeller are deteriorated. In this case, the reasonable exploitation regime of MPP may  
varies from normal exploitation regime (normal regime).  
During operating the marine vessel’s MPP, the captain chooses the revolution regime of the DME (from  
the Control Bridge) or gives a command to the engine control room (ECR), which is recommended by  
the Chief Engineer in order to have the desired velocity in practice. The engineer- officer staff, who  
chooses the exploitation revolution regime of main propulsion shaft (MPS), firstly follows the captain  
bidding, but need to ensure the working safety of MPP: DME is not overloaded, MPS does not fall into  
the area of dangerous revolutions; secondly ensures the Specific Fuel Oil Consumption (S.F.O.C) with  
possible minimum, attached to supervision of environmental norms in chosen DME regime.  
To have necessary knowledge and skills to manage and operate the MPP in different operating regimes,  
the education and training electro-mechanical crew, as well as deck officer with the management level  
is only deployed by independent simulation software for self-studying, self-improving ability or by  
intergraded software’s module built in the simulation system of ECR. The essence of building imitation  
software is to establish mathematical model and algorithm to select reasonable exploitation regime for  
MPP in all situations that may happen in reality. In the process of researching this complicated issue,  
we can choose one of mathematical models with different accuracy and coherence, which is illustrated  
360  
by general cargo ship 34000 DWT, used DME of brand MAN 6S 46MCC built in 2013 (in Pha Rung  
Shipyard, Hai Phong city, Vietnam).  
Choosing reasonable exploitation regime for MPP of marine vessel, which is set up for educating and  
training engineering cadets, and officer training courses, is a fundamental problem. In different levels  
of training programme, the principle of diesel operation follows the guide of diesel operation provided  
by the engine manufacturer [1, 2, 3] to ensure the DME the thermo-mechanical safety, the working  
regime with the minimum S.F.O.C (ge =min). In technical file provided by the ship builder or diessel  
manufacturer as well as the data of sea trial test and the certification of the diesel, there are characteristic  
curves that allows to exploit engine, by which we may assess the area of the MPP exploitation. For  
example, for the ship 34000 DWT mentioned above, following the engine document of MAN 6S  
45MCC, engine factory [6] provides the table of test data in 4 loading regimes with power 25-50-75-  
100% of the normal power and the measured parameters: engine revolution (nE, rpm); fuel pump index  
(FPI, hp), maximum combustion pressure (Pz, Bar), maximum compression pressure (Pc, Bar), mean  
effective pressure (MEP, Pme, Bar), the Specific Fuel Oil Consumption (S.F.O.C, ge, g/kWh), air intake  
pressure (Pscav, Bar), speed of turbo-changer (nT, rpm). In the ship sea trial test document [7], the data  
table of measured test regimes and parameters: mean revolution of propeller (np, rpm), ship speed (V,  
mile/h, or Vp, m/s), torque (M, kNm), relatively the 4 test regimes of consumption power (Np%) is given.  
In addition, the operation guidance of DME engine MAN 6S 46 MCC [5] shows the characteristics of  
DME working limit and characteristics of propeller  
limited by operating revolution as figure 1.  
The limited power characteristics, given in the figure  
1, provide us rather full information to maintain safely  
exploitation DME and propeller. However,  
characteristic line ge(n) and speed area of dangerous  
resonance are still not given.  
It is essential for exploitation to get the information  
and maintain ship speed according to the propeller’s  
revolution (through engine revolution and the  
transmission ratio of gear-box). Therefore, the graph  
of characteristic V-n combined with power  
characteristics (NE(n), Np(n)) [1,2,4] is often used.  
Figure 1 Power Diagram for MPP, used  
In essence, the problem of choosing reasonable  
exploitation regime is optimal control problem, the  
main control parameter is revolution of main engine  
=/n (relative revolution by normal revolution  
n (radian/s), n rpm) in the reality of exploitation  
conditions (X- vector of parameters affecting the  
choice of control law), consists of choosing  
parameters in order to set up initial control regime for  
MAN 6S 46 MCC [5]  
Line 1- Propeller curve through point A; Line 2 –  
Propeller curve, heavy running, recommended limit for  
founled hull at calm weather condition; Line 3 Speed  
limit; Line 4- Torque /speed limit; Line 5 Mean  
Effective pressure limit; Line 6 Propeller curve, light  
running (rang 2.5 -5.0%) for clean hull and calm  
weather condition; Line 7 Power limit for continuous  
running; Line 8 Overload limit.  
the voyage (slow change): x1CLI (Cargo Load Index,  
from 0 -100%); x2 and x3surface roughness state of  
hull and propeller relatively (changing by the number  
of working months after the ship last dry dock); x4vessel in shallow water; x5vessel in canal with  
limited depth and width and fast changing parameters (or stochastic) in the voyage; x6 wave and wind,  
X=[x1, x2, ..., x6]. Two parameters, which are affected directly to control the engine revolution in the  
MPP with DME (with turbo-changer), are the fuel pump index, hp (or relative position =hp/hpn) and  
scavenge air pressure (Pk (Bar), or = Pk/Pkn) to ensure working process of DME in the revolution  
regime n (rpm), or relative revolution (=n/nn, or = /n), at the same time ensure appropriate  
economical, technical and environmental parameters. Model of control exploitation revolution is chosen  
as one of these following formal equations:  
361  
= ( X);  
(1)  
(2)  
= (, , X),  
The next model depends on rotational current velocity state c(rad/s):  
= (, , X, c)  
(3)  
The system of essential models in controlling exploitation regime consists of:  
- Limiting characteristic models in exploiting engine, propeller and ship hull relatively on the  
coordinates of power- revolution- velocity;  
- Current state models of engine power, propeller’s consumption power, ship speed following real  
exploitation condition X=[x1, x2, ..., x6];  
- Models that give the solution for reasonable operation regime.  
In this article, the authors present the establishing the basic models as above and some results received  
for MPP of general cargo ship 34000 DWT used DME 6S46 MCC, which was built in Vietnam in 2013.  
2. Mathematical models of characteristics on the coordinates N n V  
2.1 Limited characteristics  
Input database (DB) in building model is limited characteristic curves, given by engine manufacturer  
and shipyard. This DB is relative and contains noises, thus the data, obtained from the graph and used  
to input into the DB input data table, does not require high accuracy. It is important to investigate the  
reliability of the model according to the appropriate standard statistic, such as Fisher standards. The  
building method of regression models bases on the Least Square of Errors Method (LSEM) .  
Mathematical basis of building regression model of experimental data in function ym=f(a, x), contains  
required coefficient vector a =[a1, a2,…, an], at a time x1, x2,…, xN and measuring value relatively y1,  
y2,…, yN is paraphrased shortly as below:  
Error between required regression model and measuring values at a certain measured time k: k= f(a, xk)  
yk. The coefficient vector a is found in order to minimize the sum of squares of the errors so we have  
N
function: J ( f (a, x ) y )2 min  
k
k
k1  
Since then, we take the partial derivative of each coefficient am, with m =1n,  
and  
n : J / am 0  
receive the system of m simple linear equations with m unknowns of elements in vector a. Resolve this  
equation system we get vector with solution a. Received solution of regressive equation is checked  
reliability (accuracy) by statistical method, Fisher function.  
Limited characteristic in figure 1, although graphs’ figure can be seen as straight lines, but in fact they  
could be high-order polynomial functions, for example, characteristic lines of propeller are often  
exponential function Np = c.n3. As researching results of choosing model, the authors realized that power  
characteristics of diesel propullsion plant - propeller have general figure and the best suitable is third-  
order polynomial [2]:  
y =(a1 + a2.x + a3.x2)x  
(4)  
in which, coefficients aj (j=1,2, 3) are found by input data and algorithm use LSEM  
The authors create file writing in MATLAB for statistical processing and establishing regressive  
functions.The result of models applied to ship 34000 DWT is shown in the table 1 by input data  
corresponding to lines (Line 1 8) in the figure 1.  
362  
Table 1 Regression model of limiting characteristics for MPP used MAN 6S46MCC  
Model’s coefficient (2)  
a1 a2 a3  
Fref.(4,ns,0.99) Fcalculate  
Name of characteristic lines  
Line 1- Propeller curve through point A  
Line 2 Propeller curve, heavy running  
Line 2A Total Load curve, heavy running,  
NL = NP + NGen.  
0.0000 0.0350 0.9661 12.06; ns=5  
0.0001 -0.1493 1.1058 12.06; ns=5  
7.11e+03  
1.51e+03  
0.0001 0.2077 0.7913 12.06; ns=5  
5.27e+03  
-
-
-
-
-
Line 3 Speed limit, = 1.05  
Line 4- Torque /speed limit  
Line 5 Mean Effective pressure limit  
-0.0002 1.0707 -0.0420 6.22; ns=10  
1.0000  
1.02e+03  
Line 6 Propeller curve, light running  
(range 2.5 -5.0%) for clean hull and calm 0.0000 0.0037 0.8320 12.06; ns=4  
weather condition  
4.10e+03  
Line 7 Power limit for continuous  
running, 1.0 <<1.05  
Line 8 Overload limit.  
1.1  
0
0
-
-
-0.0001 1.0979 0.0064 5.95; ns=11  
1.97e+03  
Utilization characteristic when exploit MPP of vessel 34000 DWT with DME MAN 6S46MCC is the  
revolution’s area of torsional vibration with dangerous resonance BARR(1)=(33 43) (rpm), or BARR  
=[0.26, 0.33] in case of all cylinders work normally and only one misfire cylinder (the worse is cylinder  
number 5), the area of torsional vibration with additional dangerous resonance: BARR(2) > 66 (rpm),  
orBARR(2)> 0,51 (by calculating table of shaft’s torsional vibration, approved by Register DNV [7]).  
Choosing operation revolution regime following the fuel-saving norm expressed by the S.F.O.C. ge  
(g/kWh, SFOC). Hence, models ge() are reference characteristics. According to document [5], we build  
model in the form of second-order polynomial with input data of 16 samples from graphs and get results  
of relative quantities SFOCrel. = SFOCabs./SFOCmin:  
2
SFOCrel. = 1.1425 0.3982+ 0.2781,  
with high reliability, because: Ft = 5.98e+03 >> F(3, 11, 0.99) = 6.2.  
According to received model, we have values in the table below  
Table 2 Calculated parameters ge.rel.( )  
50  
55  
60  
65  
70  
75  
80  
85  
90  
95  
100  
%  
ge.rel. 1.0129 1.0076 1.0037 1.0011 1.0000 1.0003 1.0019 1.0049 1.0094 1.0152 1.0224  
In order to select the operation regime, we perform chracteristics of engine power, propeller and ship  
speed by rotational velocity (relative ) of engine [1, 3].  
Limited characteristic of ship speed according to revolution of engine determined through revolution of  
propeller (by the ratio of transmiting velocity iEp =E/P, which is always determined in each regime of  
gearbox, iEp =1 without gearbox using clutch- direct propeller). Ship limit speed [V()] (knots)  
P
received when propeller’s revolution chosen by changing band, in accordance with propeller’s  
revolution band , in the condition of exploiting that relative to building limited power characteristics  
P
(in figure 1).  
It is seen that when the engine (propeller) revolution is less than some value , ship speed does not  
P0  
exist (V=0) because the propeller thrust is not dominant the hull resistance. From the value > ,  
P0  
ship speed is determined by the relationship between propeller power and hull:  
363  
1t  
RV HR M pp ;H  
,
(1wt )  
where: V ship speed (m/s);Mpmoment of propeller (Nm); R ship hull resistance (N); H Hull  
hydrodynamic efficiency; R Relative rotative efficiency; t Thrust deduction; wt wake fraction.  
Studying propeller’s power characteristic Np(p) for general cargo ship, we can choose the model Np(p)  
= C(p)3, where: C- coefficient, which represents ship resistance-propeller, depends on: their roughness  
state, cargo loading and marine condition. We build the model: C =C(X).  
The total ship resistance is summation of component resistances and each component is proportional to  
square of ship speed. We have relationship between ship speed and propeller revolution:  
, in which, hull resistance’s coefficient kH depends on density of water (kg/m3),  
kH1/3V (H p C )1/3p  
wetted surface S (m2) and received from experiment or analytical model of hull total resistance in  
exploitation conditions. This coefficient also depends on X, so we have: kH=kH(X).  
Then, we have first-order relation between ship speed and propeller revolution:  
V=  p; v = ,  
(5)  
p
Concerning the hysteresis of system, we build the model in the following equations:  
V = (- ) =   - 0  
(6)  
(7)  
p
p0  
p
v = - 0; or v = X - X0  
p
p
where: perfomance operating characteristics curve = X = (X) may be derived by experimental or  
analytical methods; V(m/s) and v =V/Vn are absolute and relative ship speed, Vn - speed of ship at  
nominal operating condition Vn = V(X=Xn).  
We use the limited resistance coefficient [n] to describe the ship limited speed when the propeller works  
on the maximum propeller characteristics with nominal operating condition Xn (curve 2, fig. 1). Likely  
to equation (7), we change to equation for the limited speed:  
[vn] = [n] - [n0]  
(8)  
p
The significant factor is that, in the most severe, the main engine can sustain upper DME limited power  
characteristics. The method to establish this perfomance operating characteristics curve is as follows:  
Propeller torque M = D5n2KQ, KQ - Propeller torque coefficient. This coeficient can be described  
approximately with a line form. KQ = a2 b2(v/n) therefore, we obtain the computational models of the  
power of main engine, which Prof. Emil Stanchev defined [2, page 13]:  
M = A20n2 B20nV, (Nm)  
where: A20, B20 - coefficients depending on propeller parameters.  
A20 = D5a2, (kg.m2); B20 = D4b2, (kg.m)  
(9)  
(10)  
Where - density of seawater, kg/m3 (=1025 kg/m3); D diameter of propeller (m); a2 and b2 can be  
calculated from propeller torque coefficient.  
From (9), we have the limited speed equals the power of main engine with limited torque characteristics  
M0() = N0/; = n/30:  
364  
[V] = (A20/B20)- M0()/B20, (m/s)  
(11)  
Hence:  
A
[v] k1 20   k2  
B20  
M 0  
B20  
k1  
k1  
can be changed to dimensionless.  
The coefficients  
and  
The curve of limited speed based on torque characteristics of engine is a hyperbol curve.  
Then, we define A20 and B20 from equation (9) base on sea trial data with LSEM.  
With 34000DWT ship using DME, MAN 6S46 MCC, base on the sea trial data [7], we obtain A20 =  
6.1374 and B20 = 6.0041 with the following model:  
M = 6.1374 2 6.0041V, (kN.m);  
Fisher standard: Fr = 878.67 >> F(2,4,0.99) = 18, prove that the reliability of model is 99%.  
The limited speed equation (11) with limited torque of engine:  
[V] = 1.0222- 0,1666 M0()/, (m/s)  
(12)  
(13)  
If "Fuel Pump Index" unchanges at the maximum value, the limited speed curve is hyperbol form.  
2.2. The current operation characteristics  
To establish the mathematical model of MPP in real operation condition by assuming X vector and the  
operation conditions are determined (or set up the concrete operation), we will establish the operating  
characteristics of engine power, propeller and ship speed.  
NE(, XE), NP(, XP), and V(, XH),  
(14)  
Where: XE,XP and XH - the concrete operation condition vectors for engine, propeller and ship hull,  
corresponding to the factor condition of X vector which was mentioned above.  
Models are established (14) by analytical methods or experimental method, and by combination of  
analytical and experimental methods.  
2.2.1. Establishing the charateristics of MPP using analytical method  
The advantage of analytical method is that, it is more flexible and has the sense of initiative to describe  
the different conditions from X operation condition to obtain the relationships (4). The only difficulty  
of analytical method is that, we need experimental data with a large of range fluctuation, or the reliability  
is low. To overcome that disadvantage, we need to take test in some favourable conditions, therefore,  
the analytical data can be confirmed.  
Ship resistance at various operation conditions  
The total resistance in the still water. The total resistance in the still water includes the following  
components:  
(15)  
R RF RV RW RAPP RA RAA , kN  
Where: RF frictional resistance; RV form resistance; RAPP appendage resistance; RA roughness  
resistance; RAA Air resistance.  
365  
The discussion of selecting the formulas to determine the resistance components in (15) is represented  
concretely in [4].  
Added resistance in waves: when ship operates in the condition of wave and wind, ship resistance  
increases. The increase of resistance will depends on height of wave, the direction of wave propagation  
in comparison with the movement of ship. According to [10], the added resistance due to wave influence  
is determined as follows:  
RAW RA0W  
(16)  
where: β - is a function taking into account operation angle compared to propagation direction of wave  
0
RAW  
(), determined by using Fig. 2;  
- is added resistance of ship in the irregular incoming wave (that  
means head-on incoming wave 180o compared with movement direction). These quantities are  
determined by formula (17) and (18) as follows:  
RA0W g(1t)(L /100)3(1/ qw2 ){A 100(CB 0,5)2 A A2 (L / B)   
0
1
(17)  
A (L /10T) A4 (xB /10L) 10A (ry / L) M[A A (L /10T)]10A Fr}  
3
5
6
7
8
(1,2 0,1 CWP ) (2,9 CWP 1,9)Fr.  
(18)  
R
Coefficients in (17) and (18) are determined as follows:  
2
2
m
D2  
10h  
5
L C  
J y  
3%   
qw ks2 ; ry   
; A Bim qw  
(i 0,1,...,8).  
WL  
i
  
11,4CB 12  
L
m1  
where: ks coefficient depending on the wave type (ks = 1 with developed wave, ks = 0.6 with developing  
wave); Jy mass inertia moment for the axis Oy; L, B, D, T is length, width, depth and draft in turn;  
ρ density of water; t thrust deduction; M effective coefficient of bulbous bow (M=1 if a ship has  
bulbous bow; conversely, M = 0); Bimcoefficient, determined by Table 3.  
Figure 2 The curve determining the influence of angle of attack [10]  
366  
Table 3 The value of Bim [10]  
m
3
i
1
2
4
5
0
1
2
3
4
5
6
7
8
13.1320 -103.71 222.0600 -170.78 44.981  
-1.6109  
1.5981  
0.1151  
33.442 -49.2910 35.4960 -9.9726  
-14.900  
6.5228  
7.8786  
-0.6628 -0.0248  
5.0745 -1.5133  
-7.6263  
-5.5428 48.5920 -74.5600 52.5470 -14.0310  
-11.1200 98.9780 -152.070 105.25 -27.5260  
-7.5139  
2.0858  
-8.843  
51.670 -58.0050 28.5910 -5.1947  
-16.263 -17.7170 -9.6519 2.2724  
62.3330 -73.3620 -37.075 -6.7179  
Added wind resistance being head-on a ship are determined as follows [10]:  
2
R
C  
a V 2 S ,  
WAA  
WAA  
A A  
where: CWAA added wind resistance coefficient (CWAA = 0,7 for cargo ship); ρa air density; VW –  
wind velocity; SA the projection area of ship above water line to the midship section.  
Added ship resistance in the shallow water. As many previous research suggests that when ship runs  
in the shallow water, there will be an increase in the resistance. The clear effect of flow depth on  
resistance will be seen apparently when  
(
; where: h- water depth [m], V- ship  
Frh V gh  
Frh 0,3  
speed [m/s]) [11].  
Added ship resistance in the shallow water is determined as follows [11]:  
1
Rh (CF CR )V 2S.  
(19)  
h
h
2
where:  
,
CF C R  
- the frictional resistance coefficient and residual added resistance in the shallow water  
h
h
respectively.  
bf  
T
h
CF CF 1  
1; CR (Cw k1CF )(CR2 1),  
(20)  
h
h
h
where: CF coefficient of frictional resistance in deep water; CW coefficient of wave resistance in deep  
water; k1CF coefficient of viscous resistance; CF, CW, k1 refer to [4]  
L
lg Re  
B
b
6,630,884  
.102 ;  
(21)  
(22)  
f
B lg Re3,35 T  
1  
1,644  
2   
T
T
T
B
   
   
   
C 10,618  
[2,0611,6(CB 0,40)3 ]0,333 10,314 3,59 5,50  
ln  
,
  
Rh  
h
h
h
6T  
Added ship resistance in the canal. When a ship runs in a canal, there has an interaction between the  
ship and the canal because of the limitation of width and depth of the chanel. This impact makes the  
367  
change in water resistance of a ship. As a result, additional resistance causing by the influence of width  
and depth is calculated [11]:  
V 2S  
2
R CK  
;
(23)  
K
Where: CK the coefficient of water resistance of a ship in the canal, determined as follows:  
C
C C  
F R  
K K  
(24)  
K
Where:  
,
- the frictional resistance coefficient and residual added resistance in the canal. They  
C F CR  
K
K
are determined as follows:  
1  
T
0,335  
0,322  
5 3   
1  
CF 0,285CB  
24,4m  
0,54  
C
(25)  
(26)  
F
K
K
h
1,26 (v vKP  
)
lg Re3,16  
  K   
1   
CR (CW k1CF )V 2  
K
These quantities in formula (25) and (26) are determined as follows:  
11,195 mK (10,345mK2 )  
V
VKP  
gh ,  
m
2C BT [h (B B )].  
M K W B  
Fr   
;
K
K
h
1
10,88b  
gh  
K
i
3
mK  
V
V   
1(2,44 4,48mK ) 11  
;
b (B B ) (B B )  
W
B
W
B
i
1mK  
VKP1  
where: bi = (BW BB)/(BW + BB) - Trapezoidal coefficient; BW width of canal’s water surface; BB –  
width of canal’s bottom; hK canal depth; VKP1 first limmited speed of ship in the canal; CM midship  
section coefficient.  
Formula (25) and (26) are correct in the range as follows: 106Re 2.109; V 0,95VKP1; 0,04 mK   
0,30; bi 0,60; 0,450 CB 0,920; 0,2 T/hK 0,90.  
Ship resistance when there is a change in ship deadweight. Throughout the voyage, deadweight of  
ship varies incessantly because of the change of fuel, cargo weight and so on,… It leads to the variation  
in the draft and eventually in ship resistance. In order to assess the change in resistance of a ship due to  
the change draft we use the following formula [9]:  
(27)  
Ri ki R  
where: Ri ship resistance in the random loading case; RShip resistance full loaded condition;  
coefficient ki is determined according to [9, page 201].  
Characteristic model of propeller  
The work of combination “propeller – ship hull – main engine” causes the mechanical and  
hydrodynamic impact between each component.  
Performance of the propulsion combination η is calculated by [12]:  
368  
1t  
R 1wt  
strdp  
,
(28)  
where: ηs shaft performance (ηs = 0,97 ÷ 0,98 when engine room locate at the end of a ship); ηtrd  
actuator performance (ηtrd =1 when engine have no gear box, ηtrd =0,98 ÷ 0,99 when engine have gear  
box); ηp propulsion performance; ηR relative rotative efficiency; t thrust deduction; wt wake  
fraction.  
The discussion of selecting the formulas to determine the thrust deduction, wake fraction, coefficient  
with the effect of erratic wake on hydrodynamic moment Q of propeller and propeller performance of a  
ship in still water is represented concretely in [4]. In other operation condition, such as in ballast  
condition (or a part load), the wake fraction tends to be 515% larger than the wake fraction in the  
loaded condition. The value can be determined by equation of Moor and O’Connor [9].  
Algorithm determining the relationship between velocity - revolution - engine power  
Algorithm determining the relationship between velocity - revolution - engine power in the still water  
is presented in [4]. In other operation condition, algorithm is the same one, therefore, it is no longer  
necessary to discuss about it.  
Essential input parameters for determining the relationship between velocity - revolution - engine power  
include: geometric parameters of ship hull and propeller, wave and wind level, water depth, geometric  
parameters of canal ship pass through.  
Example. Using above theory to calculate the relationship between velocity - revolution - engine power  
of bulk carrier 34000DWT built in PhaRung Shipyard [7].  
Input parameters: L = 179.95 m; Lpp = 176,75 m; B = 30 m; TF = 9,75 m; TA = 9,75 m; T =9,75 m; CB =  
0,8137; CWP = 0,9606; CM = 0,995; XB = 3,849 m; ABT = 14,85 m2; hB = 5,85 m; Cstern = -22; DP = 5,6 m;  
Z = 4; Zp = 1; P/D = 0,73; AE/AO = 0,85.  
Operation condition: when a ship is in the wave and wind h3% = 6 m, wave propagation compared to  
ship movement βw = 1350, ship type ks =1, wind velocity VW = 19 m/s; when a ship is in the shallow  
water: water depth h = 13m; when a ship is in the canal having following parameters: depth hK = 25 m;  
surface width Bw = 50 m; bottom width BB = 40 m; limited velocity in the canal VKP1 = 6 m/s; when a  
ship is in ballast condition: ship draft Ti = 5.51 m.  
The calculation results are shown in fig.3.  
Figure 3 The relationship curves between ship velocity - revolution - engine power in different  
operation condition  
369  
2.2.2. Building the characteristics of MPP by experimental method  
In operation condition, in some allowable regime, we measure three basic parameters: Engine power  
DME, propeller revolution and ship velocity.  
By using measured figures (or indirect calculating), some above parameters allow us to estimate the  
cooperation of components in MPP. Regression model is built based on database measured in  
experimental condition, using LSEM.  
- The propeller characteristics is built by model, formula (4);  
- Engine power can be indirectly calculated or measure by available devices;  
- Character V =V(n) is modelled in the form of straight line, formula (7).  
2.3. Selection the reasonable revolution for operation of MPP  
Figure 4 Principle scheme of the automatic selecting the  
reasonable operation mode of the MPP  
Based on modelling the limited characteristics of ship diesel MPP, working characters of ship hull,  
propeller in the current operation condition (as well as prognose working characters of ship hull,  
propeller in the alleged condition) allow us to give database ZREF (reference) in order to make a decision  
of selecting the reasonable operation regime. In the Fig.4, we suggest a principle scheme of the  
automatic selecting the reasonable operation mode of the MPP, using a electronic governor with the  
regime (revolution) selection. The Fig.4, we show the principle of controlling revolution of the DME by  
speed governor including MPU (Magnet Pulse Unit) and FVC (Frequency to Voltage Convector). The  
control signals DME in the modern diesel engine are usually fuel oil rack and scavenge air pressure  
(relative quantities are abbreviated by and ). Based on the decision of optimal selection signal op,  
the desired signal XD impacts on the PID controller as a revolution error input.  
Block optimum solution (BOS) in above principle diagram plays a key role in selecting the reasonable  
revolution condition for operation. Input data for BOS is information of current operation condition with  
the limited condition of X0 and database as well as character S.F.O.C. of engine. BOS will put forward  
the best reasonable mode according to the optimal standard so that controlling function is established,  
for instance, minimum fuel consumption in the safe range of system, or possible maximum velocity in  
the scope for MPP working safely.  
370  
Standard database includes: The limited characters {[Ne], [Np], [V]} and character SFOC ge, that are  
regressive models, some of them are shown in the table No1 in the paper; The basic operation characters  
are built by the analytical method {Ne(, X0), Np(,X0)and V(,X0)}, according to the above mentional  
bases in the point 2.2.  
3. Conclusion  
This paper builds the system of regression models for limited characteristic curve of DME, limited  
character of the ship speed when engine works above maximum constant moment characteristic curve;  
characteristic curve of the S.F.O.C. of DME and propeller characteristic curve in the full loaded  
condition. Those models have an essential reference data for maintaining the operation mode in a safe  
and economical way.  
Based on analytical models determining ship resistance, moment and thrust of propeller in every  
operation mode MPP will be the prediction of propeller characteristic curve and ship hull so that we are  
able to select a suitable operation mode.  
This article puts forward principle diagram of the automatic selecting the reasonable operation regimes;  
specifically it puts forth round mode of main engine and propeller according to reference data (standard  
characteristic model system) as well as actual operation mode of ship hull, propeller (these cases depend  
on working condition vector X of MPP).  
References  
[1] Bogdanov H.Marine Energy Equipment (Operation). Military Edition. Sofia, 1983. (In  
Bulgarian).Богданов Х.Г.Корабни енергични уредби. София, «Военно издателство», 1983. 244  
с.  
[2] Emil Stanchev and others. Dynamic researching on the marine propullsion plants in the exploitation  
conditions. Military Edition. Sofia, 1982. (In Bulgarian). Емил Станчев и др. Динамични  
изледвания на корабния пропилсивен комплекс в эксплуатационни условия. София, «Военно  
издателство», 1982. – 284 с.  
[3] Ovsianicov M. and Petukhov V. Marine propulsion plants. Leningrad, 1987. (In Russian).  
Овсянников М.К, Петухов В.А. Дизели в пропульсивном комплексе морских судов:  
Справочник. Л.: Судостроение, 1987. 256 с.  
[4] LuuD.D., Nho L.C, Tu Tr.Ng.Modelling characteristics of the hull and propeller of sea-going  
general cargo ship in normal operating conditions. Maritime Science-Technology Journal, N46,  
March 2016 (VietNam). ISBN.  
[5] MAN DIESEL 6S46 MC-C. VOLUME I. Operation.  
[6] Hồ sơ tàu 34000 DWT. Technical File 6S46 MC MK II. HHMShanghai Hujiang Diesel Emission  
Test Scie-Tech Co.Ltd. Approved DNV 20 -08-2007.  
[7] Hồ sơ tàu 34000 DWT. Hồ sơ thử tàu 34000 DWT. NMĐT Phà Rừng.2007.  
[8] Hồ sơ tàu 34000 DWT. Torsional Vibration Calculation Report DE 6S46MC C7860 kW @ 129  
r/min 34,000 DWT. HUDONG HEAVY MACHINERY Co. Ltd. PhaRung Shipyard. DNV  
(Classification Society) Approved 2007 - 06 -06.  
[9] Molland, Anthony F. Ship resistance and propulsion: practical estimation of ship propulsive power  
. ISBN 978-0-521-76052-2 (hardback).  
[10] Под ред. Я.И. Войткунского. Справочник по теории корабля: В трех томах. Том 1.  
Гидромеханика. Сопротивление движению судов. Судовые движители. (1985). Л.:  
Судостроение - 768 с.  
[11] В.Ф. Бавин, В.И. Зайков и др.; Под ред. В.Г. Павленко.Ходкость и управляемость судов  
(1991). Учебник для вузов/ – М.: Транспорт. 1991. 397 с.  
[12] Гайкович А.И. Теория проектирования водоизмещающих коряблей и судов. TOM 2: Изд.  
НИЦ МОРИНТЕХ, 2014.  
371  
pdf 12 trang yennguyen 20/04/2022 1620
Bạn đang xem tài liệu "Modelling reasonable operation regimes of the main propulsion plant main diesel engine - propeller - hull on the general cargo ship", để tải tài liệu gốc về máy hãy click vào nút Download ở trên

File đính kèm:

  • pdfmodelling_reasonable_operation_regimes_of_the_main_propulsio.pdf