Calculation Methods For Shaft Alignment
Calculation Methods For Shaft Alignment
The content of this service document is the subject of intellectual property rights reserved by DNV GL AS ("DNV GL"). The user
accepts that it is prohibited by anyone else but DNV GL and/or its licensees to offer and/or perform classification, certification
and/or verification services, including the issuance of certificates and/or declarations of conformity, wholly or partly, on the
basis of and/or pursuant to this document whether free of charge or chargeable, without DNV GL's prior written consent.
DNV GL is not responsible for the consequences arising from any use of this document by others.
DNV GL AS
FOREWORD
DNV GL class guidelines contain methods, technical requirements, principles and acceptance
criteria related to classed objects as referred to from the rules.
If any person suffers loss or damage which is proved to have been caused by any negligent act or omission of DNV GL, then DNV GL shall
pay compensation to such person for his proved direct loss or damage. However, the compensation shall not exceed an amount equal to ten
times the fee charged for the service in question, provided that the maximum compensation shall never exceed USD 2 million.
In this provision "DNV GL" shall mean DNV GL AS, its direct and indirect owners as well as all its affiliates, subsidiaries, directors, officers,
employees, agents and any other acting on behalf of DNV GL.
CHANGES – CURRENT
Changes - current
This is a new document.
DNV GL AS
CONTENTS
Contents
Changes – current.................................................................................................. 3
Section 1 General.................................................................................................... 5
1 Objective..............................................................................................5
2 Scope................................................................................................... 5
3 Application........................................................................................... 5
4 Abbreviation list.................................................................................. 5
Changes – historic................................................................................................ 25
DNV GL AS
SECTION 1 GENERAL
Section 1
1 Objective
The propulsion line aftmost bearing life can be extended and maintenance cost reduced by using more
comprehensive shaft alignment calculation methods for the load on the bearing. This guideline is a
supporting document for the DNV GL Shaft align(2) class notation.
The objective of the computational fluid dynamic (CFD) calculation is to estimate and assess the maximum
shaft forces and moments induced by hydrodynamic propeller loads which shall be used as input to the shaft
alignment analysis at design and off-design conditions.
2 Scope
The class guideline provide details of methods used to carry out viscous CFD calculations to predict the
hydrodynamic propeller loads and finite element analysis (FEA) to estimate the propeller shaft aft bearing
contact pressure. In addition, an alternative for aft bearing lubrication criteria for exceptional aft bearing
designs and logic for laser based measurement are provided.
3 Application
This guideline may be applied to calculate propeller forces and aftmost bearing pressures for the DNV GL
rules for shaft align class notation, Shaft align(2), see DNVGL-RU-SHIP Pt.6 Ch.2 Sec.10.
4 Abbreviation list
Table 1 Abbreviations
Abbreviation Description
CG class guidelines
DNV GL AS
SECTION 2 CFD ANALYSIS
Section 2
1 General
The geometry and CFD model shall be of sufficient detail level and quality to represent the dominating
physical effects needed to obtain shaft forces and moments of satisfying accuracy, considering the objective
of the analysis, which is to predict the bearing loads.
DNV GL has a considerable amount of ship scale measurement data from various ship types and
corresponding CFD computations. This will be used in the approval of submitted measurement and analysis
reports.
2 Geometry
The geometry shall represent the as-built design of the hull, appendages and other geometry features that
may influence the findings of the analysis. Any assumptions to geometry features that may and may not be
considered of importance, as well as the reason why, should be included in the analysis report.
Both port and starboard side of the ship shall be modelled.
Typically, features that shall be included are:
— hull geometry
— propeller with hub
— nozzles
— gondola
— shaft
— bracket and struts
— rudder
— tunnel thrusters
— stabilizer fins
— bilge keels
— any propulsion improvement devices that may influence the results such as:
— ducts
— pre-swirl stator
— vortex generators
— propeller boss cap (PBCF)
— etc.
3 Modelling
The model shall be sufficiently large in extent to avoid influence from the boundaries on the results of
interest. Boundary conditions shall be selected so that the physics in the far field is representative to the
objective of the calculation. At least one pressure outlet condition should be included to ensure correct static
pressure in the model.
The mesh shall be sufficiently refined to resolve expected flow features and the propagation of these features
to a level where the results no longer respond to refinement.
A grid sensitivity study should either be made, or justified through reference to previous similar studies.
Typical flow features in need of special attention is large scale separation and vortices from any appendages
or hull details that may be convected to locations at or close to locations within the reach of the propeller
blades. To resolve the boundary layer, a prism layer of sufficient resolution and thickness shall be used
around all no-slip surfaces.
The mesh and boundary conditions shall be documented.
DNV GL AS
Turbulence models and other models directly affecting the field variables should be justified and referred to in
Section 2
the report.
The propeller shall be modelled as physically rotating in the frame of reference under consideration, and a
time-step sufficiently small to capture the static loads should be applied.
4 Acceptable assumptions
Reasonable simplifications and assumptions may be applied and stated in the report. Typical examples are:
— Gravity has no or negligible effect on the resulting velocity field.
— Deformations of shaft and propeller has negligible impact on the forces and moments.
— Turbulence may be represented by a RANS closure model.
— Cavitation is assumed not to influence the results.
5 Definitions
Table 1 Symbol list
n order of harmonic
an amplitude of harmonic N
Trim Tr-Ta m
DNV GL AS
Section 2
5.1 Coordinate system
Coordinate system shall be defined and reported, an example of coordinate system is given in Figure 1:
DNV GL AS
Section 2
5.2 Moments and forces
The results from the load calculations should be given as mean values and harmonic components according
to the following convention:
(1)
where:
5.3 Draught
The vessel draft to be used in the calculations are:
Trim = Tf - Ta
The turning characteristics can be directly computed in the propeller load computation or imposed as motions
taken from either ship trial, models tests or separate manoeuvring computation. The background of the
imposed motions should be justified and presented in the report.
DNV GL AS
Section 2
Figure 2 Basic ship's motions
DNV GL AS
The report should contain at least:
Section 2
— main particulars of ship and propeller including definition of propeller coordinate system
— operating conditions and definitions including possible limitations
— torque limit curve as a function of rpm of the propulsion motor
— geometry definitions
— description of physical models and justification of their appliance
— mesh characteristics including snapshots of the most important areas
— boundary and initial conditions
— assumptions
— results for all operating conditions:
— forces and moments for all 6 degrees of freedom in the origin of the propeller coordinate system
— time series for all components for one revolution where the peak loads occur
— thrust variation per blade as a function of circumferential position for the same revolution as above
— pressure distribution, excluding the effect of gravity, on the entire propeller from back and face side
when one blade is located at 12 o'clock position
— pressure distribution, excluding the effect of gravity, on the aft ship at the same time instant as above
— total wake velocity field including the propeller action in a plane upstream of the propeller , typically a
distance of 20% of the propeller diameter
— pressure and velocity magnitude should be presented on all surfaces in the domain, including bounding
surfaces
— conclusions including a short justification of the results.
See App.A for illustration of an example and recommended format.
DNV GL AS
SECTION 3 STERN TUBE BEARING CONTACT MODELLING
Section 3
1 General
Calculation of hydrodynamic propeller loads by use of CFD gives forces and moments in the plane of rotation,
not only in vertical direction. By using contact elements between the shaft and bearing surface in a model
of the shaft system, it is possible to calculate the gap in 3D. The gap is defined as the clearance between
the shaft surface and the bearing surface, and is a measure of the shaft vs. bearing misalignment (not the
oil film thickness). It is possible to visualize this gap with contour plots for the bearing surface. The gap and
gradient of gap is important for evaluation of the bearing load capacity. Both relative slope and area with a
gap less than a given value will be used. A larger load carrying bearing area gives less edge pressure and
thereby better bearing performance and reduced risk of premature bearing failure.
The aft stern tube bearing contact modelling, as presented in this guideline, is an advanced methodology
used to optimize designs considering extreme propeller forces. The methodology has been developed over
several years by DNV GL based on research and numerous troubleshooting cases on actual installations.
2 Modelling requirements
The model of the propulsion shaft system shall be a complete model of the low speed shaft line, either
including the engine bearings or the output gear shaft bearings.
The minimum requirement is to make a 3D model of the two aft most bearings. For a single stern tube
bearing configurations, this means that the aft intermediate bearing should be modelled in 3D. For
conventional ship designs it is sufficient to model the two stern tube bearings only in 3D. The bearing support
structure can be assumed to be infinitely stiff. The hull structure is by far more stiff than the deflection of the
shaft, and modelling the complete hull will cause an unnecessary large model. Therefore, it is acceptable to
do this simplification. An exemption from this could be when the aft bearing is supported by struts or the aft
stern tube hull structure is particularly slender. Then these structures should be included in the model.
All bearing offsets and sloping (vertical and horizontal) in accordance with the shaft alignment calculation
shall be included in the model. It is not necessary to model the bearing white metal layer. However the
sloped bearing surface and the bearing oil grooves shall be included. The shaft outside the 3D bearings can
be beam elements if a proper connection between beam and solid elements are made. The representation
of the other bearings shall include the bearing clearance. Normally the intermediate bearing, the gear
output shaft bearing and the engine bearings have downward acting forces from the shafts and the bearing
clearance can be represented by a deflection. The figure below shows a typical model of a shaft line system.
In a simplified model as shown in Figure 1 the propeller is represented with a mass element in the COG
of the propeller. The correct propeller mass adjusted for buoyancy shall be applied. The propeller loads
(moments and forces) shall be applied in the COG of the propeller, see Figure 2.
DNV GL AS
Section 3
Figure 2 Simplified FE-model of propulsion shaft system, stern tube
3 Mesh requirements
The mesh shall be made sufficiently fine to capture pressure gradients between the shaft and the bearing.
The same mesh size should be used in the axial direction of the bearing giving square elements. In areas
were contact is not expected a less dense mesh can be used. It is recommended to do a check of the
element size by refinement. The refinement should have small influence on the results. The same mesh
density shall be used for the 3D shaft inside the bearing. Preferably there should be a node to node match for
the initial model. A typical mesh is shown in Figure 3.
On an indicative basis, an initial mesh division of 100 elements around the circumference is expected to be
suitable.
20-node brick elements should be used. This is normally the best type of elements to be used for curved
surfaces. 10-node tetrahedral element may also be used, but dependent of the contact element formulation
this may lead to a requirement for significantly finer mesh to obtain the required quality of results.
DNV GL AS
Section 3
4 Load cases
In addition to the load cases defined in the rules DNVGL-RU-SHIP Pt.6 Ch.2 Sec.10 [7.2.2], a load case with
a downward bending moment of 30% of nominal torque should be calculated. This will be a reference case
for the CFD load cases, and a reference to the shaft alignment calculation with 30% downward propeller
moment and minimum speed from the oil film criteria given in DNVGL-RU-SHIP Pt.4 Ch.2 Sec.4 [2.1.6]. The
calculated maximum contact pressure in the CFD load cases will be compared to the 30% downward case.
5 Results
Contour plots with contact pressure, penetration and gap (clearance between shaft and bearing surface) for
all load conditions shall be reported. It is recommended to use 40 µm as a maximum contour limit to be able
to evaluate the gap pattern. Typical plots for a comparison between a single slope and a double slope bearing
is shown below.
The reason for setting limits on the penetration and having a maximum value for the contour gap is based
on a typical minimum oil film thickness of 20 µm. However, this size should be seen together with the
penetration values. More penetration will give larger contour gap areas. The size of the gap area will show
how the shaft surface is aligned with the bearing surface.
Figure 4 Contact pressure, typical difference between single and double slope
DNV GL AS
Section 3
Figure 5 Contour plot corresponding to Figure 4, gap between shaft and bearing surface, gap
contour < 40 µm
Calculated maximum contact pressure for the extreme propeller load conditions from CFD predictions shall
be compared to the contact pressure for the 30% downward condition. It is required that the 30% downward
condition fulfils the bearing lubrication criteria given in DNVGL-RU-SHIP Pt.4 Ch.2 Sec.4 [2.1.6]. Also, the
gap contour plots for the extreme load conditions should be carefully evaluated. In these conditions the
shaft could be in contact in both rear and front end of the bearing and the gap contours will reveal the shaft
position in the bearing. If the shaft contact is in diametrical opposite positions in the rear and front end of the
bearing there will be a significant increase in contact pressure and risk of bearing damage during extreme
propeller load conditions.
An aft stern tube bearing, designed for extreme load conditions, will require some iterations before an
optimized design is achieved. Results from two different designs shall be reported to document the sensitivity
to design changes in way of change in maximum contact pressure.
DNV GL AS
Section 3
Figure 6 Contour plot with gap < 40 µm, MCR load condition
6 Evaluation criteria
The maximum calculated contact pressure is dependent of several parameters such as contact element
formulation/settings, shaft material, bearing material and mesh size. The mesh size effect should be
eliminated by sufficient fine mesh, but the other parameters could be different between finite elements
solvers and choice of material. Therefore, this guidance does not give an absolute value for the maximum
contact pressure. However, the maximum contact pressure shall be checked against the 30% downward
load case (where an oil film is retained). Increase in maximum contact pressure during transient turning
conditions should preferably not exceed 1.5 times the results from the 30% downward load case, however
maximum contact pressure up to 2.0 times may be acceptable.
The results for maximum pressure should be seen in comparison with the gap contour plots. The gap
gradient (rate of gap vs distance) will be a measure of the relative slope between the shaft and the bearing.
DNV GL AS
SECTION 4 ALTERNATIVE METHOD FOR CALCULATION OF OIL FILM
Section 4
THICKNESS
1 General
The DNV GL rules DNVGL-RU-SHIP Pt.4 Ch.2 Sec.4 [2.1.6] specify an analytic formulation for calculation of
oil film thickness. For complex bearing geometries, such as trumpet shape, the analytic formulation is not
applicable. Further, we have also experienced that the model becomes too conservative in case the bearing
geometry consists of three sloped segments or more. In such cases, more sophisticated models for prediction
of oil film thickness should be used with a more direct use of the Reynolds equation (less simplifications
than in the analytical formula). The applied theory with simplifications should be described and justified, and
simplifications vs. the Reynolds equation should be specified.
DNV GL AS
APPENDIX A CFD REPORT RECOMMENDED FORMAT
Appendix A
1 General
The below example is given as guidance and inspiration and is not meant to be a complete report.
Justification of the chosen methodology, assumptions, definitions and solution should be presented and
discussed in detail in the final report.
2 Main particulars
General
Hull
CB 0.78
Main engine
Number of sets 1
Number of cylinders 6
Propeller
Number of sets 1
Number of blades 5
Diameter 5.0 m
P0.7R/D 0.70
Ae/A0 0.57
DNV GL AS
Appendix A
3 Operating conditions
3.2 Turning
Entrance Speed Power
Conditions Draught aft Entrance RPM
[kn] [kW]
4 Geometry definitions
Description of geometry and justification of exclusion of geometry features.
5 RANS/CFD model
Documentation of the domain, mesh, boundary and initial conditions and justification of chosen physical
models and assumptions.
6 Results
DNV GL AS
Appendix A
6.2 Harmonic analysis
Table 3 Forces and moments at MCR, design draught
Mean n=1xz
an θn [deg]
Thrust [kN] 1200 120 -60
DNV GL AS
Appendix A
6.3 Time series at peak loads
DNV GL AS
Appendix A
Figure 2 Forces during one revolution
DNV GL AS
Appendix A
Figure 3 Moments during one revolution
DNV GL AS
APPENDIX B LASER MEASUREMENT
Appendix B
1 Laser measurement layout
Figure 1 A typical laser sighting process with 0-0 offsets measured in way of the aft and forward
seal flanges
Note:
The number of measurement points applicable inside bearings and stern tube housing are detailed in DNVGL-RU-SHIP Pt.6 Ch.2
Sec.10 [9.2.1]
---e-n-d---o-f---n-o-t-e---
DNV GL AS
CHANGES – HISTORIC
Changes – historic
There are currently no historical changes for this document.
DNV GL AS
About DNV GL
DNV GL is a global quality assurance and risk management company. Driven by our purpose of
safeguarding life, property and the environment, we enable our customers to advance the safety
and sustainability of their business. We provide classification, technical assurance, software and
independent expert advisory services to the maritime, oil & gas, power and renewables industries.
We also provide certification, supply chain and data management services to customers across a
wide range of industries. Operating in more than 100 countries, our experts are dedicated to helping
customers make the world safer, smarter and greener.