Sfa Vessels Guide July22
Sfa Vessels Guide July22
July 2022
GUIDE FOR
The source document for this Guide is the Guidance Notes on Spectral-Based Fatigue Analysis for Vessels
(2004). Promoting the Guidance Notes to a Guide adheres to ABS standard practice as the source
document was published more than ten years ago. Furthermore, classification notations are customarily
provided by only Rules or Guides. Promotion to a Guide alleviates any confusion caused by deviation from
standard practice.
The technical content in the Guide has not been significantly changed from that in the Guidance Notes.
However, it has been significantly reorganized such that it more closely aligns with the procedure and steps
commonly used to conduct spectral-based fatigue analyses for ship structures. Nomenclature has also been
modified such that it is more consistent throughout the document. The organization of the Guide is detailed
in Section 2 . A flowchart depicting the spectral-based fatigue procedure is given by 2/1 FIGURE 1.
The table below shows how the existing Guidance Notes have been reorganized in the creation of this
Guide.
1 Introduction 1 Introduction
5 Spectral Analysis 9 Fatigue Life (Damage) The discussion of spectral analysis was
moved to its own section and
expanded.
The July 2022 version adds a new Subsection 6/4 to provide S-N data for aluminum by referring to the
ABS Rules for Building and Classing Light Warships, Patrol and High-Speed Naval Vessels.
This Guide becomes effective on the first day of the month of publication.
ABS GUIDE FOR SPECTRAL-BASED FATIGUE ANALYSIS FOR VESSELS • 2022 iii
GUIDE FOR
CONTENTS
SECTION 1 Introduction.......................................................................................... 7
1 Purpose and Applicability (1 October 2017)................................... 7
3 Background (1 October 2017)........................................................ 8
5 Areas for Fatigue Strength Evaluation............................................8
7 Methodology and Assumptions.......................................................8
In addition to the simplified fatigue strength criteria required for classification by ABS, the Owner may
wish to apply more extensive Spectral-based Fatigue Analysis (SFA) techniques to the vessel’s structural
systems. It may be an added objective of these Spectral-based Fatigue Analyses to demonstrate a longer
design fatigue life than that required for classification.
Spectral-based Fatigue Analysis techniques are used in addition to the SafeHull Fatigue Assessment
technique, a Permissible Stress Range method (discussed in Subsection 1/3). The fatigue life of each
critical location in the structural system is assessed for adequacy. The critical locations are to be selected
using the results of the SafeHull Fatigue Assessment technique which is to be employed in the overall
structural design and analysis effort.
The list of critical structural locations which are to be subjected to Spectral-based Fatigue Analysis is to be
submitted to ABS for approval.
Provided that Spectral-based Fatigue Analysis is conducted in accordance with the procedures included in
this Guide or equivalent, ABS will grant the optional classification notation, SFA (years) or SFA
(years, WWT). In either case, as a minimum, a vessel is to meet the fatigue strength criteria in the Marine
Vessel Rules as described above. The SFA (years) notation is granted if the design fatigue life value is
equal to 20 years or greater under the wave conditions of unrestricted service defined in 5/5 TABLE 1. The
SFA (years, WWT) notation is granted if the design fatigue life value is equal to 40 years or greater
under the wave conditions of worldwide trading service defined in 5/5 TABLE 2.
If a vessel is primarily operating in seas that exceed the worldwide trading pattern (i.e., a vessel operating
half or more than half of its life in onerous seas including North Atlantic or North Pacific), SFA (years,
WWT) is not applicable. It should be noted that the calculated fatigue lives for different design wave
conditions may produce significantly different fatigue lives.
The value in parentheses is the design fatigue life in years specified by the applicant in 5-year increments
starting at 20 for SFA (years) and 40 for SFA (years, WWT). The structural system as a whole is
analyzed to verify that the calculated fatigue life values for the entire system meet or exceed the design
fatigue life. The actual service life of a vessel is dependent on many factors such as initial design,
operational and maintenance schemes. The SFA (years) or SFA (years, WWT) notation denotes the
design fatigue life of a vessel and is not a guarantee that the vessel or structure will achieve the design
fatigue life.
In order to obtain the SFA (years) or SFA (years, WWT) notations, SFA methods must be applied to
the entire vessel.
For vessels complying with Part 5A and 5B “Common Structural Rules for Bulk Carriers and Oil Tankers”
of the Marine Vessel Rules, the design fatigue life for Spectral-based Fatigue Analysis is equal to 25 years
or greater in 5-year increments.
Part 5A and 5B “Common Structural Rules for Bulk Carriers and Oil Tankers” of the Marine Vessel Rules
requires a design fatigue life of 25 years with wave loads representing the North Atlantic and multiple
fatigue assessment techniques are employed. A simplified stress analysis method based on beam theory is
presented in Section 5A-9-4 of the Marine Vessel Rules, a finite element stress analysis approach is
described in Section 5A-9-5 of the Marine Vessel Rules and structural design following the detailed design
standards is discussed in Section 5A-9-6 of the Marine Vessel Rules.
Supplementary to the SafeHull Fatigue Assessment technique or methods employed for vessels complying
with Part 5C of the Marine Vessel Rules, ABS may require the use of additional or alternative techniques to
demonstrate the fatigue strength adequacy of structural components. These techniques may include
Spectral-based Fatigue Analysis methods. In many instances the structural details cannot be adequately
analyzed via the permissible stress range fatigue assessment approach. An optional classification notation,
FL (years), may be requested in cases in which the owner or designer increases the target fatigue lives of
some or all of the structural components above the 20 year minimum value. This notation is discussed in
the Marine Vessel Rules and does not require the use of SFA methods.
A typical spectral fatigue analysis for a structural location is to evaluate its fatigue strength by comparing
its stress range distribution against its fatigue capacity. The following definitions are used in the context of
this Guide:
● Stress Range Distribution: Stress range probability density functions calculated per this Guide.
● Fatigue Capacity: S-N data (S-N curves) representing the number of stress cycles at fatigue failure.
● Fatigue Strength: Fatigue life (or damage) calculated per this Guide.
● Fatigue Demand: Design fatigue life.
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 A2, 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:
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 A2. It
should be noted that the contents of Appendix A2 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:
1.1.3 Determination of Ship Motion and Wave Pressure RAOs (Subsection 3/11)
The primary objective of seakeeping analysis is the determination of Response Amplitude
Operators, which are mathematical representations of the vessel responses and load effects to unit
amplitude sinusoidal waves. The ship motion and wave pressure RAOs are to be calculated for
each vessel loading condition.
simplified approach such as beam theory may be used instead. This Guide is written with the
assumption that a finite element model will be created.
The global structural finite element model should adequately represent the stiffness and inertial
properties of the vessel. In order to analyze fatigue critical details, the creation of fine mesh
models is typically necessary.
● Internal tank pressures and cargo hold loads due to fluid and cargo inclinations and
accelerations,
● Inertial forces on the masses of structural components and, as applicable, significant
equipment items.
1.3.3 Determination of Stress Transfer Functions (Subsection 4/1)
The stress transfer function for each location where the fatigue strength is to be evaluated should
be determined via a suitable structural analysis approach. Finite element analysis (FEA) is a
commonly used approach.
The stress transfer functions are required in order to derive the stress range distributions.
Therefore, structural analysis is a fundamental task of spectral-based fatigue analysis.
1.7.2 Calculation of Fatigue Damage and Life (Subsections 6/5 and 6/7)
For a given time span, the short term fatigue damage corresponding to each individual sea state
can be calculated from the stress range distributions, the selected S-N data and the probability of
the sea sate. The long term fatigue damage can be calculated by applying a linear damage
summation technique such as Miner’s Rule to the short term fatigue damages.
In order to determine the suitability of the structure for a design fatigue life, it is necessary to
calculate the total long term fatigue damage for the design fatigue life. One can use the total
fatigue damage and the design fatigue life to calculate the fatigue life of a detail.
FIGURE 1
Schematic Spectral-based Fatigue Analysis Procedure
A seakeeping analysis is to be performed at 75% of the vessel’s design speed for each wave frequency,
vessel heading angle and selected loading condition to be used in the spectral analysis. Normally, the
frequency range to be used is 0.1 rad/sec to 1.80 rad/sec in increments not larger than 0.1 rad/sec.
However, depending on the characteristics of the response, it may be necessary to consider a different
frequency range. The wave heading range is to be 0 to 360 degrees in increments not larger than 30
degrees.
Because of the variability of the loading conditions and their impact on fatigue strength predictions, it is
necessary to consider more than one loading condition in the fatigue analysis. At a minimum, two cases,
the most probable deepest and shallowest drafts that the vessel is expected to experience during its service
life, should be represented and used in the Spectral-based Fatigue Analysis process.
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.
● Displacement: ±1%
● Trim: ±0.5 degrees
● Draft:
– Forward ±1 cm
– Mean ±1 cm
– Aft ±1 cm
● LCB: ±0.1% of length
● SWBM: ±5%
Additionally, the longitudinal locations of the maximum and the minimum SWBMs and SWSFs and, if
appropriate, those of zero SWBM and zero SWSF should be checked to ensure proper distribution of the
SWBM along the vessel’s length.
For the load component types and structural responses of primary interest, analysis software formulations
derived from linear idealizations are deemed to be sufficient.
RF = 0 . 5 1 . 0 + tanh(0 . 35d)
where d is depth, in meters, of the field point below the still-water line.
In order to correctly account for the intermittent wetting effects, it is recommended that the vertical height
of the panels near the waterline be no greater than two times the longitudinal spacing. If the panels are too
large, the pressure reduction factor can be applied to the structural model after the wave pressures are
mapped from the hydrodynamic model.
Motion and load components are presented in ship fixed coordinate system in this section. The origin is
located at the intersection of the aft perpendicular, centerline and baseline. The positive sign convention is
longitudinal forward, transverse port and vertical upwards (See 4/1 FIGURE 1).
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.
R = distance vector from the vessel’s center of gravity to the calculation point
The liquid tank pressure for each of the tank boundary points can be calculated from the following
equation, which is expressed as a combination of the static and dynamic pressure components:
1/2
p = po + ρℎi (g + aV)2 + (gT + aT)2 + (gL + aL)2
where
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 ℎi is calculated in the direction of g.
FIGURE 2
Liquid Pressure on a Completely Filled Tank
FIGURE 3
Liquid Pressure on a Partially Filled Tank
The bulk cargo pressures acting on the internal surfaces of cargo holds are to be calculated and applied to
the structural model for analysis. Static and dynamic bulk cargo pressures should be included in the
analysis assuming there is no relative motion between the cargo hold and contained bulk cargo.
The following definitions are used in the formulas for the static and dynamic load calculations in this
Paragraph.
n = surface unit normal vector in the ship coordinate system, towards the inside of the cargo
hold
= nx i + ny j + nzk
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 α
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 = 0 otherwise
pt = ρℎTa sinα0 sinα cosα 0° ≤ α ≤ 90°
pt = 0 otherwise
where
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 ℎT is calculated in the direction of g.
FS = mg
where
The static load due to a stack of containers may be summed and applied to appropriate nodes on
the cargo hold bottom plate. Total vertical load due to the containers on the deck may be applied to
the appropriate nodes on the hatch coaming top plates.
FV = maV
where
The transverse component (see 4/3 FIGURE 5) of dynamic container load can be calculated using the
following equation which is expressed as a combination of the quasi-static and inertial load components:
FT = m(gT + aT)
where
The transverse load due to containers may be distributed to appropriate nodes on the bulkhead structure via
the container cell guide. The total transverse load due to the containers on deck may be applied to the
appropriate nodes on the hatch coaming top plates via the container lashing system.
The longitudinal component of dynamic container load can be calculated using the following equation
which is expressed as a combination of the quasi-static and inertial load components:
FL = m(gL + aL)
where
FIGURE 5
Vertical and Transverse Container Load Components
The static load acting on a node of the lightship structure and equipment can be expressed as:
FS = mg
where
The dynamic components of the load acting on the node in vertical, transverse and longitudinal directions
can be calculated from the following equations:
FV = maV
FT = m(gT + aT)
FL = m(gL + aL)
where
The formulas for FT and FL above represent a combination of the quasi-static and inertial load components.
The stress transfer function for each location of fatigue strength evaluation needs to be calculated for each
load case. It is likely that the number of analysis load cases will be very large. The recommended
parameters for a spectral-based fatigue analysis for each considered loading condition are 18 wave
frequencies (0.1 to 1.80 rad/s in increments of 0.10 rad/s), 12 wave headings (0 to 360 degrees in 30 degree
increments) and the real and imaginary wave components. In the case of a standard analysis, the total
number of load cases per loading condition will be (18 · 12 · 2) = 432. In addition, a static load case may
be needed. Provided that the shape of each RAO and its local maxima are adequately captured, larger
frequency steps may be used to reduce the number of analysis load cases.
Reference should be made to additional ABS guidance on modeling and analysis techniques for vessel
structure (e.g., the ABS Guidance Notes on SafeHull Finite Element Analysis for Hull Structures). While
the partial hull model used for a SafeHull analysis is different from the full ship model used for an SFA
analysis, the modeling practices and considerations presented are applicable. Detailed modeling guidance
relating to items such as element types, mesh sizes, dependence between local and global models, etc., are
also presented in the ABS Guidance Notes on SafeHull Finite Element Analysis for Hull Structures.
While creating the model, a sufficient number of nodes, elements and degrees of freedom are to be used to
adequately represent the stiffness and inertial properties of the vessel.
Three types of elements are typically used to approximate the geometry, configuration, and stiffness of the
vessel structure:
● The selection of an S-N data class (see Subsection 6/3 and Appendix A1) that applies for each
instance.
● Creating the fine mesh FEM adjacent to the detail and calculating the stress at the hot spot location via
extrapolation of the stress distribution.
4/7 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 A1) 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 Marine Vessel Rules.
FIGURE 6
Definition of Hot Spot Stress
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 A2.
The Bretschneider or two parameter Pierson-Moskowitz wave spectrum is recommended for the North
Atlantic, described by the following expression:
Hs2 2π 4 −5 2π 4 −4
S(ω) = 4π Tz ω exp − π1 Tz ω
where
Spectral-based fatigue analysis typically uses environmental data for ocean waves that are given in a wave
scatter diagram format. A wave scatter diagram consists of cells, or sea-states, which are short-term
descriptions of the sea in terms of joint probability of occurrence of a significant wave height, Hs, and a
characteristic period. Each cell contains three data items:
Section 5, Tables 1 and 2 present the wave scatter diagram data that is to be used in the spectral-based
fatigue analysis of a vessel classed for unrestricted service and for world-wide trading service, respectively.
It can be assumed that there is an equal probability of each vessel heading relative to the direction of the
waves.
Once the stress transfer functions are obtained via a structural analysis and the wave spectra are selected
for each of the cells in the wave scatter diagram, the stress spectra and spectral moments can be calculated.
With the available stress spectra, the Rayleigh probability density function (pdf) describing the short term
stress-range distribution, the zero up-crossing frequency of the stress response, and the stress bandwidth
parameter used in calculating rainflow correction can be calculated for each sea state in the wave scatter
diagram. The aforementioned terms are defined in Appendix A2 which presents a detailed example of a
closed form spectral-based fatigue analysis procedure. The stress range is considered to be twice the stress
amplitude.
TABLE 1
ABS Wave Scatter Diagram for Unrestricted Service Classification
(Representing the North Atlantic)
3.5 2 113 1332 4599 6488 4716 2092 642 149 28 20161
10.5 1 9 37 70 76 53 26 10 282
11.5 4 18 36 42 32 17 7 156
12.5 2 9 19 24 19 11 4 88
13.5 1 4 10 14 12 7 3 51
>14.5 1 5 13 19 19 13 7 77
Sum over 8 326 3127 12779 24880 26874 18442 8949 3335 1014 266 100000
All
Heights
Note:
* Wave heights taken as significant values, Hs
** Wave periods taken as zero crossing values, Tz
TABLE 2
ABS Wave Scatter Diagram for World-Wide Trading Service Classification (1
October 2017)
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
8.5 1 7 21 44 72 82 63 34 14 5 343
9.5 3 10 20 34 42 34 20 9 3 175
10.5 1 5 10 17 22 19 11 5 2 92
11.5 1 2 5 9 12 10 7 3 1 50
12.5 1 3 5 6 6 4 2 1 28
13.5 1 1 2 3 4 3 1 1 16
>14.5 1 2 3 5 5 4 2 1 23
Sum over 2500 10691 19704 22543 19628 13425 7127 2983 1024 297 78 100000
All
Heights
Note:
* Wave heights taken as significant values, Hs
** Wave periods taken as zero crossing values, Tz
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 A2 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 A1 of this Guide and Part 5C of the Marine 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 Marine 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.
Note: source United Kingdom’s Dept of Energy (HSE) Guidance Notes, 4th Edition.
It may be necessary to make adjustments to the S-N data to reduce the fatigue capacity of the detail in
some instances. Conversely, if the detail has characteristics beneficial to fatigue performance, the S-N data
may be adjusted to increase its fatigue capacity.
A reduction in fatigue capacity may be required to account for factors such as the absence of corrosion
protection (coating) for structural steel or relatively large plate thickness. The reduction is to be done in
accordance with standard ABS practice. The fatigue capacity of a detail may be increased to account for
conditions such as compressive mean stress effects – a high compressive portion of the acting variable
stress range – or the use of weld improvement techniques.
The use of a weld improvement technique such as weld toe grinding or peening to relieve ambient residual
stress can be effective in increasing fatigue life. However, such effects should not be considered in the
design of the structure. Typically, an increase in fatigue capacity through weld improvement techniques
will be considered only for situations arising during construction, operation or future reconditioning of the
structure. If the design fatigue life cannot be obtained using preferred design measures such as layout
refinement, geometric changes, scantling adjustment or weld profile modification to reduce fatigue
damage, an exception may be granted.
Weld improvement techniques such as weld toe grinding or ultrasonic peening may be used to improve
fatigue life only if the calculated fatigue life is sufficient when the weld improvement techniques are not
considered. The required calculated fatigue life without considering weld improvements for various design
fatigue lives are listed in the table below:
20 years 15 years
25 years 17 years
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 A2 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.
The fatigue life for each loading condition is to be calculated separately. If the exposure time ratios have
been considered in the fatigue life calculations for each loading condition, the combined fatigue life is
given by:
1
LC = 1 + 1 + · · · + 1
αs
L1 L2 Ln
where
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 p2 p
αs + +...+ n
L1 L2 Ln
where
∑ni = 1 pi = 1 . 0
FIGURE 1
S-N Curves
SmN = A
where
S = stress range
N = number of cycles to failure
A, m = parameters representing the intercept and inverse slope of the upper (left) portion of the S-
N Curve. These change at N = 107 cycles to C and r, respectively. Values of these
parameters are given in the following table.
TABLE 1
Parameters for Basic S-N Design Curves
A m C r
Class (For MPa units) (For MPa units)
B 1 . 013 × 1015 4 1 . 020 × 1019 6
Refer to Part 5C of the Marine Vessel Rules for guidance on the categorization of structural details into the
indicated 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.
Notes:
1 (1 October 2017) This Appendix is referred to in Section 6. It is provided to describe the formulations comprising
a Spectral-based Fatigue Analysis approach, which can be employed to satisfy the criteria to obtain the SFA
(years) or SFA (years, WWT) Classification notations. However, it is at this formulation level that there are
multiple valid methods that may be selected. For that reason, it is emphasized that the contents of this Appendix
are provided primarily to illustrate principle rather than to give mandatory steps for the Spectral-based Fatigue
method.
2 The procedure described below considers the use of a wave scatter diagram (e.g., 5/5 TABLE 1). Where a different
base period for the wave scatter diagram is employed, the procedure must be suitably modified.
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.
2
Sσ ωHs, Tz, θ = Hσ ωθ Sη ωHs, Tz 1
3) Calculate the spectral moments at the same ship forward speed, V, used in the seakeeping analysis.
The nth spectral moment, mn, is calculated as follows:
∞
mn = ∫ (ω − Vω2cosθ/g)nSσ(ωHs, 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 A2/3.3
FIGURE 1). Applying the wave spreading function modifies the spectral moment as follows:
θ + 90 ∞
2 n
mn = ∫ π cos2 α − θ · ∫ ω − Vω2cosα/g Sσ ωHs, Tz, α dω dα 3
θ − 90 0
FIGURE 1
Spreading Angles Definition
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σ
1 m2
f= 2π m0 5
Bandwidth Parameter:
2
m2
ε= 1− m0m4 6
where
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
D = ∑i = 1 di = ∑i = 1 Ni 7
i
where
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 m
Di = A ∫ ktkmskℎs f0ipig s ids 8
0
where
= t n
22 for t ≥ 22mm
= 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
σms = mean stress
s4 = long-term stress range corresponding to the representative probability level of 10-4
m, A = physical parameters describing the S-N curve
T = design fatigue life, in seconds
f0i = zero-up-crossing frequency of the stress response, Hz
pi = joint probability of Hs and Tz
gi = probability density function governing s in the i-th sea state
s = specific value of stress range
Summing Di over all of the sea-states in the wave scatter diagram leads to the total cumulative
damage, D. Therefore:
∞
m f0T
D = kℎktkms A ∫s
m ∑M
i = 1 f0ipigi /f0 ds 9
0
where
Introducing long-term probability density function, g(s) of the stress range as:
f0ipigi
i
g s = 10
f0ipi
i
and
NT ∞ m
D = (kℎktkms)m A ∫ s g(s)ds 11
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.
where
where
For bi-linear S-N curves (see Appendix A1) where the negative slope changes at point Q = NQ, SQ from
m to r = m + Δm Δm > 0 and the constant A changes to C, the expression for damage, as given in
equation 12, is as follows:
m m
D= T
A 2 2 Γ m/2 + 1 ∑iM= 1 λ m, εi μif0ipi kℎktkmsσi 14
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:
SQ SQ
∆m
smgids − A kℎktkms sm + ∆ mgids
C
0 0
μi = 1 − ∞ 15
∫ smgids
0
∆ m/2 ∆m
Γ0 m/2 + 1, vi − 1/vi kℎktkms Γ0 r/2 + 1, vi
μi = 1 − Γ m/2 + 1
16
where
SQ 2
vi =
2 2σi