Hull Deformation Korbetis
Hull Deformation Korbetis
Hull Deformation Korbetis
Abstract
A typical VLCC driven by a two-stroke Diesel engine is studied. A detailed finite element model of the
hull structure is generated. The propulsion shaft of the ship is modelled as a statically indeterminate
multi-supported beam. First, a reference shaft alignment plan is assumed, and the static equilibrium
of the shaft is calculated using matrix analysis. Next, different loading conditions (laden/ballast) of
the ship are assessed. For each loading condition (a) hydrostatic equilibrium of the ship is computed,
(b) the corresponding hull deformations are calculated using Finite Element Analysis, (c) the relative
vertical displacements at the bearing locations are determined and (d) the static shaft equilibrium is
re-evaluated. The computed bearing loads are compared to those of the reference case.
1. Introduction
The propulsion system of conventional cargo ships typically consists of a two-stroke Diesel engine,
and a shafting system, which transmits the engine power to the propeller, Fig.1. Radial shaft loads
(propeller/shaft/engine weights) are supported by journal bearings (stern tube bearings, line bearings,
crankshaft bearings). Proper design, installation and alignment of the shafting system of a ship is cru-
cial for stable, efficient and reliable operation, ABS (2004), NKK (2006). Primarily, shaft alignment is
concerned with the determination of proper longitudinal and vertical bearing positions, aiming at
equi-distribution of bearing loads. The successful application of a static shaft alignment plan is essen-
tial for trouble-free dynamic operation of the propulsion system, aiding in decreasing bearing wear,
increasing bearing expected lifetime and decreasing maintenance and replacement costs.
In operation, shaft alignment may be considerably influenced by hull deflections, due to different
loading and environmental conditions. The effect of hull deflections on shaft alignment is more pro-
nounced in very long ships, with relatively flexible hulls and stiff shafts. In such cases, the robustness
of shaft alignment at different loading conditions of the ship, taking hull deflections into account,
should be carefully assessed. In this respect, the use of detailed Finite Element Analyses for the calcu-
lation of hull deflections is imperative.
Recently, the subject of shaft alignment has gained increasing attention. Devanney and Kennedy
(2003) underlined the drastic deterioration of tanker newbuilding standards in the last decade, and the
corresponding effect on the reliability of the shafting system. Specifically, emphasis was put on the
severity of stern tube bearing failures in modern VLCCs and ULCCs, which may lead to loss of pro-
pulsion and vessel immobilization. The authors claimed that the main reason of this failure is the de-
sign of propulsion shafts with decreased diameters, followed by improper shafting alignment. They
546
suggested that (a) hull deflections should be thoroughly taken into account for a range of loading con-
ditions of the ship, (b) the engine room structure should be reinforced, to minimize additional offset of
the bearings, and (c) time varying loads on the stern tube bearing and heat dissipation in the lubricant
domain should be taken into account.
Šverko (2003) highlighted several design concerns in propulsion shafting, especially for VLCC and
large bulk carrier vessels. In such vessels, shaft alignment is very sensitive to hull deflections; this
behavior was attributed to the increased hull flexibility of such ships (due to scantling optimization
and increased ship lengths) and to the increased stiffness of the propulsion shaft (due to the demand
for higher propulsion power and, consequently, larger shaft diameters). If the hull deformations can be
predicted accurately, an optimal set of bearing offsets for the vessel on even keel may exhibit a rea-
sonably good performance at other loading conditions of the vessel; however, since hull deflections
cannot be easily calculated accurately, a practical solution could be to complete the alignment at dry
dock conditions, and make provisions to correct (if needed) bearing vertical offsets when the reactions
are verified afloat. Šverko (2006) addressed the problem of predicting hull deflections through analy-
sis of series of collected real life data. Hull deflections were estimated by measurement of shaft de-
flections using bending gauges. The goal of this study was to find appropriate dry dock bearing offsets
that will result in acceptable alignment performance over a wide range of vessel loading conditions.
Murawski (2005) also utilized a FEM model of a large containership, and introduced a new parameter
to be considered: the stiffness characteristics of the bearing foundations. He concluded that, in a holis-
tic approach to the shaft alignment problem, bearing stiffness and oil film characteristics of each bear-
ing should be taken into account in the design stage. Dahler et al. (2004) reported the results of an
joint industrial project between DNV, MAN B&W and DAEWOO concerned with the numerical and
experimental study of shaft deflections and bearing loads in large ships propelled by two-stroke Die-
sel engines. They utilized a complete FEM model of the ship, which exhibited a fine mesh at the aft
end of the ship hull (engine room). Focus was given on engine and crankshaft deflections and on the
corresponding bearing loads. To this end, FEM analyses were performed taking into account the real
crankshaft geometry, and the results were compared with simulations using simplified crankshaft
models. Simulation results were also compared to experimental measurements. They concluded that
FEM-hull analyses can capture the general trend of hull deflections reasonably well, but fail to ac-
count for local variations in the curvature of the shaft, leading to inaccurate predictions of bearing
loads. Finally, they suggested that by applying the final shafting plan after vessel launch, possible er-
rors due to wrong estimation of hull deflections could be avoided. BV (2013) released Rule Note NR
592, concerned with Elastic Shaft Alignment (ESA) of ships. The proposed methodology of shafting
alignment calculations takes into account hull deformations, oil film characteristics and stiffness of
the bearings’ foundation. The rule is mainly applicable to ships characterized by a propeller shaft di-
ameter greater than 750 mm, or between 600 mm and 750 mm, but with propeller weight greater than
30 tones or a prime mover with power output greater than 20 MW.
In the present work, a typical VLCC vessel, driven by a two-stroke Diesel engine, is studied. The ves-
sel has a propeller shaft diameter of 815 mm; therefore it is within the scope of the ESA Rule of BV.
Here, a detailed finite element model of the hull structure of the ship, complying with the meshing
requirements set by Classification Societies, is generated with the use of the ANSA pre-processor.
The propulsion shaft of the ship is modeled as a statically indeterminate multi-supported beam; the
bearing stiffness and clearance are taken into account, and the static equilibrium of the shaft is calcu-
lated using matrix analysis. Considering the undeformed hull of the vessel, a reference shaft align-
ment plan is assumed, and the static equilibrium of the shaft is calculated. Next, different loading
conditions (laden/ballast) of the ship are assessed. For each loading condition (a) hydrostatic equilib-
rium of the ship is computed, (b) the corresponding hull deformations are calculated, (c) the relative
vertical displacements at the bearing locations are determined and (d) the static shaft equilibrium is re-
evaluated. The computed bearing loads (reaction forces) are compared to those of the reference case.
547
2. Problem Definition
The main characteristics of the studied VLCC are presented in Table I. Finite Element Analysis is
performed to calculate the hull deformations of the vessel, at different loading conditions. Of
particular importance are the deformations at the bearing locations of the propulsion shafting system.
The static analyses are conducted with the aid of the ANSA pre-processor the MSC/NASTRAN
solver. Here, thermal loads from the engine or the environment are not taken into consideration. First,
a FEM model of the ship structure is generated. The whole structure of the ship is represented by first-
order shell elements; at the stern tube region, solid tetrahedral elements are used. A coarse mesh is
generated for the whole structure (element length of 0.95 m), except from the engine room floor,
where finer mesh (element length of 0.2 m) ensures better accuracy results, Figs.2 and 3.
Fig.3: Detail of the generated FEM mesh at the engine room region of the vessel
548
Fig.4: Detail of the generated FE mesh at the stern tube region of the vessel
The mesh generation is an automated process performed by the ANSA Batch Meshing Tool. Meshing
parameters and quality criteria are defined in two meshing scenarios (fine mesh for the engine room
floor and coarse mesh for the rest of the structure), Table II. Re-meshing algorithms act on areas with
poor mesh quality until the predefined quality criteria are fulfilled. The final model comprises of
about 402.000 shell elements, 143.000 beams and 17.000 solid tetrahedrals, Figs.2 to 4.
Stiffeners are represented by beam elements pasted on the shells. This method simplifies the model by
avoiding the generation of very small shell elements. The properties of the beam elements are calcu-
lated in accordance with the cross section of each stiffener.
Machinery, auxiliary structures and small constructions that do not contribute to ship strength are not
modeled in the present FEM model. Their mass is applied to the model as non-structural mass. This
mass is appropriately distributed over the FEM model, so as to reach the prescribed lightship weight
and the corresponding center of gravity. The mass of the present structural model is 34442 t, while the
lightship weight is 43938.7 t and its center of gravity L.C.G. at 151.338 m. Thus, 9496.7 t of lumped
masses are appropriately distributed in holds, stern and bow by the automatic process of the ANSA
Mass Balance Tool, Fig.5. The engine mass is represented by a lumped mass of 990 t distributed to
the engine foundation positions by RBE3 elements. Bearing positions where measurements will take
place are represented by single nodes on the bearing axis, connected to the engine room floor with
RB2 elements, Fig.6.
1918 tons 1360 tons 1205 tons 1334 tons 911 tons 358 tons
167 tons
168 tons 532 tons 208 tons 24 tons 130 tons 179 tons
549
Fig.6: Engine and bearings representation in the present FEM model
Three representative loading conditions of the vessel, namely full-load departure, ballast arrival and
departure with partial load, are considered in the present analysis, Fig.7. The contents of the tanks are
represented by lumped mass connected to the each hold bottom with RBE3 elements. The ship is
positioned on steel water considering the vessel’s total displacement and center of gravity. Buoyancy
is applied as pressure at the hull underneath the waterline using PLOAD4 entities, Fig.8. Finally, the
vessel is trimmed in order to achieve static equilibrium between weight and buoyancy, which makes
the model able to run without the need of displacement constraints (SPCs), which would lead to high
local stresses. A NASTRAN keyword for inertia relief (INREL) is added for this solution.
550
2.2 Calculation of Static Shaft Alignment
As noted in the introduction, a successful application of a static shaft alignment plan is important for
trouble-free operation of the ship in the anticipated service conditions. The propulsion shaft of the
ship is supported by the stern tube bearings, the line bearing(s) and the crankshaft bearings of the
main engine. At first, a reference line can be defined as the one passing through the centers of the aft
and fore stern tube bearings. Shaft alignment is concerned with the determination of the proper verti-
cal offset of the center of the remaining bearings from the reference line, that result in even pressure
distribution amongst all of the bearings of the system. This should stand both for static and dynamic
conditions of the vessel.
In the present paper, the propulsion shaft is represented by an assembly of two-node beam elements
subjected to purely flexural deformations. External loads and deformations are applied at the beam
nodes. All internal loads (e.g. the distributed weight of a beam) can be expressed in terms of equiva-
lent nodal generalized loads, through the application of basic principles of mechanics, Hughes and
Paik (2010). The degrees of freedom allowed for each node are three rotations about each axis of a
Cartesian 3D coordinate system and three displacements along each axis of the same system. We
shall, from now on, refer to all parameters related to each simple beam element as “local” parameters
and similarly, we shall denote all parameters related to the whole assembly of beams as “global” pa-
rameters.
For a single beam element, we may consider a vector f, containing the values of external and internal
nodal loads of each of the six degrees of freedom (DOFs) of each node, and a vector u with the corre-
sponding generalized displacements (i.e. displacements and rotations). A linear relationship between
the nodal generalized displacements and nodal forces is assumed, namely f = ku. Matrix k represents
the stiffness of each beam; the elements of matrix k are a function of the geometric and material prop-
erties of the beam (length, moment of inertia, Young’s modulus).
In Fig.9(a), a simple shaft consisting of four beam elements is presented. Using vectors f, u, and ma-
trix k of each beam element, a global linear relationship between generalized forces and displace-
ments of the system can be defined. To this end, vectors F and U are defined, which hold the values of
all nodal DOFs, the total number of which evaluates to six times the number of the system nodes. The
corresponding stiffness matrix of the system K (global stiffness matrix) is produced by appropriately
utilizing the local stiffness matrix of each beam. The global problem can now be defined as F=KU,
and can be solved for F or U, Hughes and Paik (2010).
The propulsion shaft of a ship can be modeled as a multi supported beam. A simple type of support is
that presented in Fig.9(a), denoted as a small triangle below the constrained node of the beam. Those
idealized supports allow zero displacements of the shaft in the radial direction. In practice, the shaft is
supported by hydrodynamically lubricated journal bearings. Geometrically, a journal bearing is a hol-
low cylinder, which encloses a solid shaft that rotates about its axis. The radius of the bearing is
slightly larger than that of the shaft; the difference between the bearing and the shaft radius is called
clearance. Therefore, the shaft may undergo a small displacement before contact with the bearing sur-
face. Further, the journal bearing foundation is also deformable, therefore it will elastically deform
551
when a load is applied. The elasticity of the bearing foundation can be taken into account by introduc-
ing an appropriate bearing stiffness coefficient. In practice, several of the bearings of propulsion
shafts will be deliberately shifted in the y direction by an appropriate vertical offset, described in the
shafting plan of the vessel. Therefore, a predefined displacement of the shaft at the bearing position
should be taken into account. The above concepts are presented in Fig.9(b).
In summary, at bearing of the shafting system, an initial y-offset may be imposed (as displayed in
Fig.9(b)); at equilibrium, the shaft will either “float” within the boundaries of the bearing clearance,
without interacting with the bearing, or it will come into contact with the bearing upper or lower inner
surface. In the latter case, the “spring-like” behavior of the bearing foundation will cause an additional
deformation of the bearing support, to such a degree, that the resulting reaction force will balance the
shaft weight that this bearing was meant to support.
(a) (b)
Fig.9: (a) Simplified model of a shafting system. (b) Sketch of a bearing; clearance, offset and
foundation stiffness.
The above considerations become even more important as we examine the overall behavior of a given
shafting system under different ship loading conditions. The deflections of the ship hull due to the
action of load and buoyancy, directly affect the vertical position of the bearings, Fig.10. This addi-
tional disturbance can be taken into account through the application of an additional vertical offset to
each bearing, relative to the reference line.
Fig.11 shows a model of the shafting system studied. The propeller shaft, the intermediate shaft and
part of the crankshaft of the main engine are considered. The propeller shaft is supported by two stern
tube bearings, the intermediate shaft by a line shaft bearing and the engine crankshaft by the
crankshaft bearings (five of them are included in the present calculations). Bearing details are
presented hereinafter:
552
Stern tube bearings: Aft bearing: L/D =2.22, foundation stiffness of 4x109 N/m. Fore bearing:
L/D=0.53, foundation stiffness of 5x109 N/m. Both bearings have a radial clearance of
0.55 mm.
Line shaft bearing: L/D =0.78, radial clearance of 0.425 mm, foundation stiffness of 109 N/m.
M/E crankshaft bearings: Foundation stiffness of 6x109 N/m, radial clearance of 0.345 mm.
Density of the shaft material: 7850 kg/m3, Young’s modulus of elasticity: 2.06x1011 N/m2.
The shaft is discretized with 39 beam elements. The geometry of each beam and load details
are presented in Table III.
A reference line of the shafting system is defined as that passing through the centers of the aft and
fore stern tube bearings. Initially, no hull deformations are considered (this case resembles dry-
docking conditions of the ship). The line shaft bearing and the engine crankshaft bearings are
appropriately offset from the reference line. In Table IV, the corresponding vertical offsets of each
bearing are presented, accompanied by the properties of each bearing and by the calculations of
reaction forces at each bearing support location.
Table IV: Initial shaft alignment plan: Bearing properties, vertical offsets and reaction forces.
Radial Mean
Bearing Stiffness Offsets Reactions
Bearing Clearance L/D Pressure
No. (N/m) (mm) (kN)
(mm) (MPa)
1 Aft S/T 0.550 4.0x109 2.221 -0.06 997 0.676
9
2 For S/T 0.550 5.0 x 10 0.528 0.00 84.4 0.241
3 Intermediate 0.425 1.0 x 109 0.780 -3.90 165 0.426
4 M/E 1 0.345 6.0 x 109 - -6.60 181 -
5 M/E 2 0.345 6.0 x 109 - -6.60 301 -
9
6 M/E 3 0.345 6.0 x 10 - -6.60 406 -
7 M/E 4 0.345 6.0 x 109 - -6.60 396 -
8 M/E 5 0.345 6.0 x 109 - -6.60 161 -
553
Table III: Discretization details of the shafting system of the present study.
Element Dist. to Diameter (m)
right External
Length
end of Load (N)
No. Type (m) Left Right
element
(m)
1 0.050 0.050 0.650 0.650
2 Load at right end 0.375 0.325 1.035 1.035 6.99E+06
3 0.405 0.030 0.650 0.650
4 Load at right end 1.372 0.967 0.726 0.775 6.14E+05
5 2.175 0.803 0.775 0.815
6 2.505 0.330 0.815 0.815
7 2.635 0.130 0.815 0.815
Bearing at right
8 2.907 0.272 0.815 0.815
end
9 4.445 1.538 0.815 0.815
10 7.835 3.390 0.815 0.815
11 7.895 0.060 0.817 0.817
Bearing at right
12 8.110 0.215 0.817 0.817
end
13 8.325 0.215 0.817 0.817
14 8.505 0.180 0.817 0.817
15 9.020 0.515 0.817 0.817
16 9.120 0.100 0.817 0.817
17 9.970 0.850 0.817 0.705
18 10.105 0.135 1.320 1.320
19 10.240 0.135 1.320 1.320
20 12.555 2.315 0.705 0.705
21 12.955 0.400 0.705 0.705
22 13.130 0.175 0.710 0.710
Bearing at right
23 13.405 0.275 0.710 0.710
end
24 13.680 0.275 0.710 0.710
25 13.855 0.175 0.710 0.710
26 17.655 3.800 0.710 0.710
27 19.200 1.545 0.705 0.705
28 19.335 0.135 1.458 1.458
29 Load at right end 19.336 0.001 1.458 1.458 5.19E+04
30 19.555 0.219 1.458 1.458
Bearing at right
31 20.205 0.650 0.980 0.980
end
32 Load at right end 20.840 0.635 0.980 0.980 5.93E+04
Bearing at right
33 21.205 0.365 0.980 0.980
end
34 Load at right end 21.955 0.750 0.552 0.552 3.51E+05
Bearing at right
35 22.705 0.750 0.552 0.552
end
36 Load at right end 23.455 0.750 0.552 0.552 3.51E+05
Bearing at right
37 24.205 0.750 0.552 0.552
end
38 Load at right end 24.955 0.750 0.552 0.552 3.51E+05
Bearing at right
39 25.705 0.750 0.552 0.552
end
554
3. Computational Results
Hull deformations have been computed for the three different loading conditions presented in Section
2.1. For loading condition 1 (ballast arrival condition), the hull exhibits a hogging behaviour, which
causes considerable displacements at the bearing positions. Figs 12 and 13 show distributions of Von
Misses stresses at the ship hull and at the engine room region. For loading condition 2, the hull is
bending towards the opposite direction (sagging), whereas for loading condition 3 a hogging at aft and
sagging at fore behaviour is exhibited. In Table V and Fig.14, the bearing offsets from the reference
line are presented for all loading conditions considered in the present study.
Fig.12: Loading condition 1 (ballast arrival condition): Distribution of Von Misses stresses on hull.
Fig.13: Loading condition 1 (ballast arrival condition): Distribution of Von Misses stresses at engine
room region.
Table V: Bearing vertical offsets at different loading conditions of the ship (distance from a reference
line passing through the centers of the aft and fore stern tube bearings)
Bearing Initial case Loading Loading Loading
(even keel) Condition 1 Condition 2 Condition 3
Aft S/T -0.06 -0.06 -0.06 -0.06
For S/T 0.00 0.00 0.00 0.00
Intermediate -3.9 -1.89 -4.66 -3.12
M/E 1 -6.60 -2.27 -8.99 -2.44
M/E 2 -6.60 -1.96 -9.32 -2.37
M/E 3 -6.60 -1.50 -9.85 -2.28
M/E 4 -6.60 -1.06 -10.44 -2.20
M/E 5 -6.60 -0.63 -11.07 -2.14
555
Fig.14: Bearing vertical offsets at different loading conditions of the ship.
In Table VI and Fig.15, the calculated bearing reaction forces are presented for the three different
loading conditions studied in the present work.
Table VI: Calculated bearing reaction forces for different loading conditions of the vessel
Bearin Loading Condition 1 Loading Condition 2 Loading Condition 3
g No. Differ- Differ- Differ-
Bear- Bear- Bear-
Reac- ence Reac- ence Reac- ence
ing ing ing
tions from tions from tions from
Offsets Offsets Offsets
(kN) Initial (kN) Static (kN) Static
(mm) (mm) (mm)
(%) (%) (%)
1 -0.06 1030 3% -0.06 976 -2% -0.06 995 0%
2 0 25.6 -70% 0 132 56% 0 138 64%
3 -1.89 193 17% -4.66 127 -23% -3.12 28.3 -83%
4 -2.27 185 2% -8.99 204 13% -2.44 468 159%
5 -1.96 300 0% -9.32 295 -2% -2.37 180 -40%
6 -1.50 407 0% -9.85 406 0% -2.28 326 -20%
7 -1.06 402 2% -10.44 416 5% -2.20 397 0%
8 -0.63 153 -5% -11.07 137 -15% -2.14 159 -1%
Reference Condition Loading Condition No. 1 Loading Condition No.2 Loading Condition No.3
1030000
997000
995000
976000
468000
416000
407000
406000
406000
402000
397000
396000
326000
301000
300000
295000
204000
193000
185000
181000
180000
165000
161000
159000
153000
138000
137000
132000
127000
84400
28300
25600
1 2 3 4 5 6 7 8
BEARING NUMBER
Fig.15: Calculated bearing reaction forces for different loading conditions of the vessel
556
Overall, although the vertical offsets of certain bearings are of the order of 10 mm, the differences in
bearing reaction forces are not very pronounced. In particular, the reaction force of the aft stern tube
bearing ranges from 976 kN to 1030 kN (maximum deviation of approximately 5%). Bearing 4 (aft
engine bearing) exhibits the maximum deviations in reaction forces, ranging from 181 kN (even keel
condition) to 468 kN (L.C. 3). Bearings 7 and 8 display the least amount of deviation.
4. Conclusions
A preliminary study of shaft alignment in a typical VLCC vessel was conducted. A detailed finite
element model of the hull structure of the ship was generated; a very fine mesh was utilized at the
engine room region of the ship. The propulsion shaft of the ship was modeled as a statically
indeterminate multi-supported beam and solved using matrix analysis. Bearing clearance and the
stiffness of the bearing foundation were taken into account. First, considering the undeformed (even-
keel) hull of the vessel, a reference shaft alignment plan was assumed, and the static equilibrium of
the shaft was calculated, yielding the reaction forces at the shaft bearings. Next, three representative
loading conditions of the vessel, corresponding to full-load, partial load and ballast conditions were
simulated. The corresponding hull deflections were computed, the offset of the bearings due to hull
deflections were determined, and the bearing reaction forces were calculated.
In general, the differences in bearing reaction forces at different loading conditions are not very
pronounced. At the aft stern tube bearing, the reaction force exhibits a maximum deviation of
approximately 5%. The bearing 4 (aft engine bearing) exhibits the most pronounced deviations in
reaction forces. The results support, for this specific case and vessel, conclusions drawn by other
researchers in recent literature: An appropriate even-keel shaft alignment plan exhibits reasonably
good performance at other loading conditions of the vessel. This study could be further extended to
account for (a) hot conditions of the engine / application of eccentric thrust loads, (b) the full range of
loading conditions of the ship and (c) detailed behaviour of the oil film at each bearing (solution of
the Reynolds equation in the lubricant domain).
References
ABS (2004), Guidance notes on propulsion shafting alignment, American Bureau of Shipping
DAHLER G.; BRODIN, E.; VARTDAL, B.J.; CHISTENSEN, H.W.; JAKOBSEN, S.B.; OK, Y.K.;
HEO, J.H.; PARK, K.R. (2004), A study on flexible hulls, flexible engines, crank shaft deflections and
engine bearing loads for VLCC propulsion machinery, CIMAC Congress, Kyoto
DEVANNEY, J.; KENNEDY, M. (2003), The down ratchet and the deterioration of tanker
newbuilding standards, Center for Tankship Excellence
HUGHES, F.O.; PAIK J.K. (2010), Ship structural analysis and design, SNAME
MURAWSKI, L. (2005), Shaft line alignment analysis taking ship construction flexibility and
deformations into consideration, Marine Structures 18, pp. 62–84
SVERKO, D. (2003), Design concerns in propulsion shafting alignment, ABS Technical Papers
SVERKO, D. (2006), A solution to robust shaft alignment design, ABS Technical Papers
557