SFA-Vessels Guide E-Oct17
SFA-Vessels Guide E-Oct17
GUIDE FOR
This Guide becomes effective on the first day of the month of publication.
Users are advised to check periodically on the ABS website www.eagle.org to verify that this version of
this Guide is the most current.
We welcome your feedback. Comments or suggestions can be sent electronically by email to [email protected].
ABS GUIDE FOR SPECTRAL-BASED FATIGUE ANALYSIS FOR VESSELS . 2016 iii
Table of Contents
GUIDE FOR
CONTENTS
SECTION 1 Introduction ............................................................................................ 1
1 Purpose and Applicability....................................................................1
3 Background .........................................................................................2
5 Areas for Fatigue Strength Evaluation ................................................2
7 Methodology and Assumptions ...........................................................2
SECTION 1 Introduction
The main underlying assumptions of the Spectral-based Fatigue Analysis method are:
i) Ocean waves are the source of the fatigue stress range acting on the structural system.
ii) The load and structural analyses are assumed to be linear as required for the frequency domain
formulation and the associated probabilistic analysis to be valid. As such, scaling and the
superposition of stress transfer functions from unit amplitude waves are considered valid.
iii) Non-linearities due to non-linear roll motion and intermittent loads, such as wetting of the side
shell in the splash zone, can be effectively accounted for using correction factors.
iv) Due to their insignificant contributions in typical load cases it is appropriate to disregard structural
dynamic amplification, transient loads and effects such as springing. This allows for the use of
quasi-static finite element analysis.
For the specific SFA method presented in Appendix 2, it is assumed that the short-term stress variation for
a given sea-state is a random, narrow-banded, stationary process. Therefore a Rayleigh distribution can be
used to represent the short-term stress range distribution.
The Spectral-based Fatigue Analysis method is applied to each of the selected structural locations by
implementing the following process:
• Determine the stress range distributions.
• Determine fatigue capacity (S-N data).
• Calculate fatigue strength (fatigue life or damage).
Section 2 provides an overview of the spectral-based fatigue analysis procedure. A schematic representation
of the SFA procedure can be found in Section 2, Figure 1.
An effort is made in this Guide to avoid the discussion of complicated formulae and instead to focus on the
concepts being presented. It is necessary to present the complex formulae used in the calculation of fatigue
damage resulting from the predicted stress range distributions, which are presented in Appendix 2. It should
be noted that the contents of Appendix 2 are intended to serve as an example of a valid SFA method. It is
not necessary that the process be followed exactly; however, any method used should adhere to the same
principles.
1 General
This Section provides an overview of the spectral-based fatigue analysis procedure for sea-going vessels.
The section in which a detailed discussion of each step can be found is included in parentheses.
The spectral fatigue analysis for ocean-going vessels covers a wide scope of technical disciplines ranging
from the analysis of ship’s Response Amplitude Operators (RAOs) to the assessment of fatigue damages of
structural details. Four main steps are typically involved:
• Seakeeping Analysis: to determine ship motion and wave pressure RAOs.
• Structural Analysis: to determine stress transfer functions.
• Spectral Analysis: to determine stress range distributions.
• Fatigue Strength Assessment: to determine fatigue life or damage.
The tasks that comprise each step are described below.
FIGURE 1
Schematic Spectral-based Fatigue Analysis Procedure
SEAKEEPING ANALYSIS
(Section 3)
Develop Hydrodynamic Model
Subsection 3/5
(Section 3)
Determine Stress
Transfer Functions
Subsection 4/1
5 Hydrodynamic Models
Boundary element methods generally require that the wetted surface of the vessel be discretized into many
three dimensional panels. The panel mesh should be fine enough to resolve the radiation and diffraction
waves with reasonable accuracy. A general rule is that there are at least four panels within the shortest
wave length being analyzed. The plots of the RAOs versus frequency should be reviewed for any abrupt
changes; they are indicative of an insufficient number of panels to fully capture the wave characteristics.
The hull offset data can be used to facilitate the generation of the panel model. For each loading condition,
the draft at the forward and aft perpendiculars, the location of the center of gravity, the radii of gyration
and the sectional mass distribution along the ship length are to be in accordance with the Loading Manual.
In addition, there should be sufficient compatibility between the hydrodynamic and structural models such
that the mapping of fluid pressures onto the structural model can be done appropriately.
FIGURE 1
Ship Coordinate System
A structural analysis is to be conducted to calculate the stress transfer functions necessary for spectral-
based fatigue analysis. The stress transfer functions are used in the spectral response analysis to calculate
the stress range distributions as discussed in Section 5.
The wave induced motion and load components described in this Section are RAOs, and the stress responses
are stress transfer functions. The following loads are to be calculated via a structural analysis: liquid tank
pressures (4/3.1), dry bulk cargo loads (4/3.3) and container loads (4/3.5). In each Paragraph, the static,
quasi-static, and inertial loads are discussed. In addition, the inertial loads of discrete masses, which are
applicable to significant items of equipment, are discussed (4/3.7).
The effect of sloshing loads due to fluid in partially filled tanks is not within the scope of the SFA (years)
or SFA (years, WWT) classification notations. However, the designer is encouraged to perform and
submit such calculations, if deemed important.
There may be situations where partial models of the structural system are used. In such instances, hull
girder shear and axial forces, and bending and torsional moments should be determined to appropriately
represent the boundary conditions at the ends of the partial model.
The general approach used in the calculation methods described below is to calculate the total stress response
considering both static and dynamic loads. Subsequently, the static stress is deducted from the total,
leaving the pure wave-induced stress response. Alternative methods and formulations that directly produce
the dynamic stress response may also be used.
3 Internal Load
The procedures to calculate liquid tank pressure, bulk cargo pressure, container loads, and loads on discrete
masses are presented in this subsection. Each of these loads is composed of static and dynamic components.
The static load component results from gravity. The dynamic load component can be further decomposed
into quasi-static and inertial components. The quasi-static load component results from gravity, considering
the roll and pitch inclinations of the vessel. The direction of gravitational forces in the ship-fixed coordinate
system varies with roll and pitch motions, resulting in the quasi-static loads. The inertial load component
results from the local accelerations caused by the ship motions in six degrees-of-freedom. It is assumed
that the ship motions are calculated in the ship-fixed coordinate system. Therefore, coordinate transformations
are not needed for the local accelerations.
The following definitions are used in this Subsection:
θ = roll angle, positive starboard down
φ = pitch angle, positive bow down
gV = acceleration of gravity. It can be considered as the vertical component of gravitational
acceleration relative to the ship-fixed coordinate system to the accuracy of the first
order.
gV ≈ g
gT, gL = transverse and longitudinal components of gravitational acceleration relative to the
ship-fixed coordinate system due to roll and pitch inclinations
= (–g sin φ, g sin θ)
aV, aT, aL = vertical, transverse and longitudinal components of local accelerations caused by ship
motions relative to the ship-fixed coordinate system at a calculation point
= A + Θ× R
A = surge, sway and heave acceleration vector at vessel’s center of gravity
Θ = roll, pitch and yaw acceleration vector
R = distance vector from the vessel’s center of gravity to the calculation point
hi = internal pressure head defined by the height of projected liquid column in the direction
of a resultant acceleration vector. For a completely filled tank, the pressure head is to
be measured from the highest point of the tank to the load point (see Section 4, Figure 2).
For a partially filled tank, the pressure head is to be measured from the free surface
level to the load point (see Section 4, Figure 3). The free surface is defined as the liquid
surface normal to the resultant acceleration vector. In the aforementioned figures,
only vertical and transverse accelerations are indicated.
aV, aT, aL = local accelerations at the center of gravity (CG) of tank contents
The liquid tank pressure p in the above formula becomes the static pressure if aV, aT, aL, gT and gL are all
set to zero, and hi is calculated in the direction of g.
FIGURE 2
Liquid Pressure on a Completely Filled Tank
z z
y aV y
θe θe
aT
hi
hi
hi
hi
FIGURE 3
Liquid Pressure on a Partially Filled Tank
z aV z
θe
y y
aT
hi
hi
hi hi hi
The following definitions are used in the formulas for the static and dynamic load calculations in this Paragraph.
αo = angle of repose for the bulk cargo considered
= 30 deg. in general, 35 deg. for iron ore, 25 deg. for cement
ρ = density of the bulk cargo
aV, aT, aL = local accelerations at the CG of cargo contents
n = surface unit normal vector in the ship coordinate system, towards the inside of the
cargo hold
= nx i + ny j + nz k
a = (aL + gL) i +(aT + gT) j + (aV + g) k
a = magnitude of a
= |a |
α = the angle between the surface unit normal vector, n , and the acceleration vector, a
(see Section 4, Figure 4)
a⋅n
= cos −1
| a || n |
Both static and dynamic bulk cargo pressures can be decomposed into normal and tangential components
relative to the surface of cargo hold. The positive normal component is defined in the opposite direction of
n , and the positive tangential component is defined in the direction of n × ( n × a ).
FIGURE 4
Depiction of Surface Normal Vector, n , Acceleration Vector, a , and α
n
α
The normal and tangential components of bulk cargo pressure can be calculated using the following equations,
each of which is expressed as a combination of the static and dynamic pressure components:
pn = ρ hTa{cos2α + (1 − sin αo) sin2α} 0° ≤ α ≤ 90°
pn = 0 otherwise
where
pn = pressure component normal to the hold surface
pt = pressure component tangential to the hold surface
hT = bulk cargo pressure head defined by the height of projected bulk cargo column in the
direction of the acceleration vector. The pressure head is to be measured from the top
of the cargo surface to the load point.
The total pressure components pn and pt in the above formulae become the static pressure components if
aV, aT , aL , gT and gL are all set to zero, and hT is calculated in the direction of g.
FIGURE 5
Vertical and Transverse Container Load Components
∇ FT ∇
FV
aV
aT
Section 4, Figure 6 depicts an acceptable method that can be used to calculate the hot spot stress for a weld
toe. Element sizes near the detail of interest are to be approximately equal to the plating thickness. The
hot spot stress is found via linear extrapolation of the calculated stress results at distances of t/2 and 3t/2
from the detail. For stresses obtained using this approach, the use of E class S-N data (see Appendix 1) is
considered to be most appropriate. The surface stresses (considering a “bending plate” element type) should
be used to determine the hot spot stress. A detailed description of the numerical extrapolation procedure
can be found in 5C-1-A1/13.7 of the Steel Vessel Rules.
FIGURE 6
Definition of Hot Spot Stress
Peak Stress
t Weld Toe
≈t
Weld Toe Location
t/2
3t/2
1 General
The spectral analysis is used to calculate the stress range distributions using the stress transfer functions
calculated via structural analysis and the selected wave spectra and a wave scatter diagram. Mathematically,
Spectral-based Fatigue Analysis begins after the determination of the stress transfer functions. A detailed
presentation of an example method is presented in Appendix 2.
H2 2
4
5 1 2 4
4
S() = s exp
4 Tz Tz
where
Hs = significant wave height, in m
TABLE 1
ABS Wave Scatter Diagram for Unrestricted Service Classification
(Representing the North Atlantic)
* Wave heights taken as significant values, Hs
** Wave periods taken as zero crossing values, Tz
Wave Period (sec)** Sum Over
3.50 4.50 5.50 6.50 7.50 8.50 9.50 10.50 11.50 12.50 13.50 All Periods
0.5 8 260 1344 2149 1349 413 76 10 1 5610
1.5 55 1223 5349 7569 4788 1698 397 69 9 1 21158
2.5 9 406 3245 7844 7977 4305 1458 351 65 10 25670
3.5 2 113 1332 4599 6488 4716 2092 642 149 28 20161
4.5 30 469 2101 3779 3439 1876 696 192 43 12625
Wave Height (m)*
TABLE 2
ABS Wave Scatter Diagram for World-Wide Trading Service
Classification (1 October 2017)
* Wave heights taken as significant values, Hs
** Wave periods taken as zero crossing values, Tz
Wave Period (sec)** Sum Over
3.50 4.50 5.50 6.50 7.50 8.50 9.50 10.50 11.50 12.50 13.50 All Periods
0.5 1888 5417 5836 3392 1260 313 55 8 1 18170
1.5 493 3739 8374 9339 6517 2985 940 219 41 6 1 32654
2.5 96 1128 3682 5934 6159 4216 1930 629 158 33 6 23971
3.5 18 294 1217 2439 3269 2959 1769 734 229 57 12 12997
4.5 4 79 388 893 1417 1568 1142 568 208 60 15 6342
Wave Height (m)*
1 General
Section 4 of this Guide addressed calculating the stress transfer functions for locations in the structure for
which the adequacy of fatigue strength is to be evaluated. Appendix 2 presents a spectral-based analysis
procedure that may be used for the calculation of the stress range distributions. The capacity of a location
to resist fatigue damage is characterized by the use of S-N Data, which are described below. Refer to
Appendix 1 of this Guide and Part 5C of the Steel Vessel Rules concerning the S-N Data recommended by
ABS.
S-N data is typically used in two ways. In the first, the nominal stress approach, it is assumed that the acting
variable stress range can be adequately determined from the nominal stress distribution (which may include
geometric stress concentration effects) in the area surrounding the location for which the fatigue life is to
be evaluated. The second method, the hot spot approach (see 4/7.3), is used for locations in which complicated
geometry or relatively steep local stress gradients invalidate the use of the nominal stress approach.
Reference should be made to Part 5C of the Steel Vessel Rules for further explanation and application of
these two approaches and for guidance on the categorization of structural details into the various S-N data
classes.
Where an improvement technique is applied, full details of the technique are to be submitted for review
along with supporting calculations indicating the proposed fatigue life improvement factor. If grinding is
used, the full details of the grinding standard, including the extents, profile smoothness particulars, final
weld profile, improved workmanship technique and quality acceptance criteria are to be clearly indicated
on the applicable drawings.
Grinding is to produce a smooth concave profile at the weld toe with a penetration depth of at least 0.5 mm
into the plate surface below the bottom of any visible undercut. It is preferable that a rotary burr be used
and that the grinding extend below the plate surface in order to remove toe defects. The ground area is to
have effective corrosion protection. The groove depth is to be minimized and generally should not exceed
1 mm. In no circumstances may the grinding depth exceed the lesser of 2 mm or 7% of the plate gross
thickness. The grinding area is to extend well beyond the high stress region.
A weld surface treated using ultrasonic peening should have a smooth finished shape and all traces of the
weld toe are to be removed. The minimum peening depth is to be 0.2 mm below the original surface.
Generally, the maximum peening depth is 0.5 mm.
A maximum fatigue life improvement of 2 times may be granted provided the above recommendations are
followed.
5 Fatigue Strength
An appropriate method is to be employed to calculate the fatigue damage resulting from each considered
sea state. The damage resulting from individual sea states is referred to as short-term. The total fatigue
damage resulting from combining the damage from each of the short-term conditions can be calculated
using a weighted linear summation technique (e.g., Miner’s Rule).
Appendix 2 contains a detailed description of the steps involved in a suggested Spectral-based Fatigue Analysis
method that follows the basic elements mentioned above. ABS should be provided with background and
verification information that demonstrates the suitability of the analytical method employed.
where
LC = combined fatigue life
Li = fatigue life for the ith loading condition (i = 1 to n) prior to accounting for non-sailing
time
n = number of loading conditions.
If the exposure time ratios have not been considered in the fatigue life calculations for each loading condition,
the combined fatigue life is given by:
1
LC =
p1 p 2 p
αs + ++ n
L1 L2 Ln
where
pi = exposure time ratio of the ith loading condition (i =1 to n) satisfying:
n
∑p
i =1
i = 1.0
FIGURE 1
S-N Curves
TABLE 1
Parameters for Basic S-N Design Curves
N ≤ 107 N > 107
A m C r
Class (For MPa units) (For MPa units)
B 1.013 × 1015 4 1.020 × 1019 6
C 4.227 × 1013 3.5 2.584 × 1017 5.5
D 1.519 × 1012 3 4.239 × 1015 5
E 1.035 × 1012 3 2.300 × 1015 5
F 6.315 × 1011 3 9.975 × 1014 5
F2 4.307 × 1011 3 5.278 × 1014 5
G 2.477 × 1011 3 2.138 × 1014 5
W 1.574 × 1011 3 1.016 × 1014 5
Refer to Part 5C of the Steel Vessel Rules for guidance on the categorization of structural details
into the indicated classes.
Notes for Application of Classes:
Class B: Parent material with automatic flame-cut edges ground to remove flame cutting drag line.
Class C: Parent material with automatic flame-cut edges and full penetration butt welds ground flush in way of
hatch corners in container carriers or similar deck areas in other vessel types.
Class D: Full penetration butt welds in way of hatch corners in container carriers or similar deck areas in other
vessel types.
1 General
In the “short-term closed form” approach described below, the stress range is normally expressed in terms
of probability density functions for different short-term intervals corresponding to the individual cells or
bins of the wave scatter diagram. These short-term probability density functions are derived by a spectral
approach based on the Rayleigh distribution method, whereby, it is assumed that the variation of stress is a
narrow-banded random Gaussian process. To take into account effects of swell, which are not accounted
for when the wave environment is represented by the scatter diagram, Wirsching’s “rainflow correction”
factor is applied in the calculation of short-term fatigue damage. Having calculated the short-term damage,
the total fatigue damage is calculated through their weighted linear summation (using Miner’s rule).
Mathematical representations of the steps of the Spectral-based Fatigue Analysis approach just described
are given next.
mn ( V 2 cos θ / g ) n S ( | H s , Tz , θ)d .............................................................. (2)
0
Most fatigue damage is associated with low or moderate seas, hence, confused short-crested sea
conditions must be allowed. Confused short-crested seas result in a kinetic energy spread, which
is modeled using the cosine-squared approach, (2/) cos2. Generally, cosine-squared spreading is
assumed from +90 to –90 degrees on either side of the selected wave heading (refer to Appendix 2,
Figure 1). Applying the wave spreading function modifies the spectral moment as follows:
θ 90
2
cos θ ( V cos / g ) S ( | H s , Tz , )d d .................. (3)
2 2 n
mn
θ 90 0
FIGURE 1
Spreading Angles Definition
0
ve
Wa
θ
α
θ - 90°
θ + 90°
4. Using the spectral moments, the Rayleigh probability density function (pdf) describing the short term
stress-range distribution, the zero up-crossing frequency of the stress response and the bandwidth
parameter used in calculating Wirsching’s “rainflow correction” are calculated as follows:
Rayleigh pdf:
s s
2
g (s) = exp − ......................................................................................... (4)
4σ 2 2 2σ
Zero-up crossing frequency, in Hz:
1 m2
f = ................................................................................................................... (5)
2π m0
Bandwidth Parameter:
m 22
ε = 1− ................................................................................................................ (6)
m0 m 4
where
s = stress range (twice the stress amplitude)
σ = m0
5. Calculate cumulative fatigue damage based on Palmgren-Miner’s rule, which assumes that the
cumulative fatigue damage (D) inflicted by a group of variable amplitude stress cycles is the sum
of the damage inflicted by each stress range (di), independent of the sequence in which the stress
cycles occur:
J J
∑ ∑N
ni
D= di = ......................................................................................................... (7)
i =1 i =1 i
where
ni = number of stress cycles of a particular stress range
Ni = average number of loading cycles to failure under constant amplitude loading
at that stress range according to the relevant S-N curve
J = number of considered stress range intervals
Failure is predicted to occur when the cumulative damage (D) over J exceeds a critical value equal
to unity.
The short term damage incurred in the i-th sea-state, assuming an S-N curve of the form N = AS-m,
is given by:
∞
T
Di = ∫ (kt kms kh s ) m f 0i pi g ( s )i ds ........................................................................... (8)
A 0
where
Di = damage incurred in the i-th sea-state
kh = a factor for high tensile steel, applicable to parent material only
= 1.000 for mild steel or welded connections
= 0.926 for H32 steel
= 0.885 for H36 steel
= 0.870 for H40 steel
kt = a factor for thickness effect, which is not applicable to longitudinal stiffeners
which are flat bars or bulb flats
n
t
= for t ≥ 22 mm
22
= 1.0 for t < 22 mm
n = 0.25 for cruciform joints, transverse T joints and plates with
transverse attachments
= 0.20 for transverse butt welds
= 0.10 for butt welds ground flush, base metal, longitudinal welds or
attachments
If it can be conclusively established that the detail under consideration is always subject
to a mean stress of σms, D is to be adjusted by a factor kms.
kms = a factor for mean stress effect, which is
= 1.0 for σms > s4/2
= 0.85 + 0.3σms/s4 for –s4/2 ≤ σms ≤ s4/2
= 0.7 for σms < –s4/2
where
D = total cumulative damage
f0 = “average” frequency of s over the design fatigue life
∑f pg 0i i i
g (s) = i
......................................................................................................... (10)
∑f p i
0i i
and
NT = total number of cycles in design fatigue life = f0T
the expression for total cumulative damage, D, can be rewritten as:
∞
D = (k h k t k ms )
NT
∫s
m m
g ( s )ds ..................................................................................... (11)
A 0
6. If the total number of cycles, NT, corresponds to a required minimum design fatigue life of 20
years, the calculated fatigue life would then be equal to 20/D. Increasing the design fatigue life to
higher values can be done accordingly. The fatigue safety check is to be done in accordance with
the applicable Rules where factors of safety (or Fatigue Design Factors) are specified.
∑ λ(m, ε ) f
T
D= (2 2 ) m Γ(m/2 + 1) i 0i pi (k h k t k msσ i ) m ........................................................ (12)
A i =1
where
Γ = complete gamma function with the argument (m/2 + 1)
∑ λ(m, ε )µ f
T
D= (2 2 ) m Γ(m/2 + 1) i i 0i pi (k h k t k msσ i ) m ................................................... (14)
A i =1
where µi is the endurance factor having its value between 0 and 1 and measuring the contribution of the
lower branch to the damage. It is defined as:
[ ]∫ s
SQ SQ
A
∫ s g i ds − (k h k t k ms )
m ∆m m + ∆m
g i ds
C
µi = 1 − 0
∞
0
.................................................................... (15)
∫s
m
g i ds
0
µi = 1 −
[ ]
Γ0 (m / 2 + 1, vi ) − (1 / vi ) ∆m / 2 (k h k t k ms )∆m Γ0 (r / 2 + 1, vi )
................................................ (16)
Γ(m / 2 + 1)
where
2
SQ
vi =
2 2σ
i
Γ0 = incomplete gamma function and is
x
∫
Γ0 (a, x) = u a −1 exp(−u )du
0