Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

BV - Ni461-1998 Guidance Notes

Download as doc, pdf, or txt
Download as doc, pdf, or txt
You are on page 1of 106

GUIDANCE NOTE

Quasi-dynamic
analysis of mooring
systems using
ARIANE software
FILE NAME:
/var/www/apps/conversion/tmp/scratch_2/355044578.doc

STATUS OF DOCUMENT on 18/02/1998:


First draft from LL (conversion required)
Converted by LR from LL's draft
Under examination by Working Group / Specialists
Editorial verification performed by LR
Under final examination by Technical Verificators (no return)
Never submitted to Ad Hoc Committee
No finalisation of comments in meetings on 12/02/98 & 17/02/98
Corrections from MF given on ________
Corrections from DL given on ________
Corrections integrated on ______________
Editorial verification performed by LR on _____________
Ready for printing
QUASI-DYNAMIC ANALYSIS
OF
MOORING SYSTEMS
USING ARIANE SOFTWARE

Guidance Note

NI 461 DTO R00 E


Tentative issue
May 1998
MARINE DIVISION GENERAL CONDITIONS
ARTICLE 1 4.3. The Society exercises due care and skill and acts in a
1.1. BUREAU VERITAS is a Society, the purpose of whose professional and workmanlike manner :
Marine Division is the classification of vessels, sea and river - in the selection of its surveyors,
units, offshore structures and craft of all kinds, collectively - in the performance of its services, considering the state
hereinafter referred to as "ships". of currently available technology at the time the services
The Society : are being performed.
- sets forth and develops Rules, Guidance Notes and
other documents; 4.4. Surveys conducted by the Society include, but are not
- publishes Registers; limited to, visual inspections, non-destructive testing, laboratory
- issues Certificates, Attestations and Reports following its testing and underwater inspection by divers. Survey practices
interventions. and procedures are selected by the Society in its sole and
absolute discretion based on its experience and knowledge and
1.2. BUREAU VERITAS also participates in the application according to generally accepted professional standards in the
of National or International Regulations or Standards, in Industry.
particular, by delegation from different Governments.
ARTICLE 5
1.3. BUREAU VERITAS can also carry out Technical 5.1. No report, statement, notation on a plan, review,
Assistance missions and provide special services exceeding the certificate of class or any document or information issued
scope of classification, to which Particular Conditions will then or given as part of the services provided by the Society
apply. shall be given any legal effect of implication other than a
representation that the ship, structure, item of material,
1.4. BUREAU VERITAS is neither, and may not be equipment or machinery or any other item covered by such
considered as, an underwriter, Consulting Engineer, Naval document or information has met one or more of the
Architect, Shipbuilder or a Shipowner, and cannot assume the Society's Rules. Any such representation is issued solely
obligations incumbent on such functions, even through its for the use of the Society, its committees and clients or
experience enables it to answer enquiries concerning matters not other duly authorized entities and for no other purpose.
covered by its Rules, Guidance Notes or other documents.
5.2. The classing of a ship, or the issuance of certificate with
ARTICLE 2 relation to, or in furtherance of the classification of a ship or the
2.1. Classification is the expression of confidence given by performance of services by the Society, is not a representation,
the Society to a ship, for a particular use or service and for a statement or warranty of seaworthiness, structural integrity,
certain period of time, by reference to its Rules, Guidance Notes quality or fitness for a particular purpose or service of any ship,
and other documents: this opinion expresses only the strict structure, material, equipment or machinery surveyed by the
personal view of the Society. It is represented by a class Society, beyond any representation set forth in the Rules of the
published in the Register. Society.

2.2. The documents issued pursuant to surveys carried out 5.3. Nothing contained herein, or in any information,
by the Society's surveyors according to the methods and report, certificate or like document issued in connection
modalities mentioned in Articles 3 and 4 reflect the condition of with or pursuant to the performance by the Society of its
the ship at the time of the survey. It is the Owner's or his Agent's services, shall be deemed to relieve any designer, naval
responsibility to maintain that condition until the next survey architect or engineer, builder, manufacturer, shipyard,
required by the Rules. seller, supplier, contractor or subcontractor, repairer, or
owner, operator or other entity from any warranty or other
ARTICLE 3 contractual obligations express or implied or from any fault
BUREAU VERITAS whatsoever, nor to create any right, claim or benefit in any
- appraises specifications and drawings for the third party.
construction of ships by reference to its Rules, In particular, the Society does not declare the
- inspects materials, machinery, equipment, fittings and acceptance or commissioning of a ship; that being the
other items used in the construction of ships, by exclusive responsibility of the owner.
reference to its Rules, The validity, application, meaning and interpretation of a
- surveys the construction of ships to ensure that they are Classification Certificate or any similar document or information
built in accordance with the Society's published Rules, issued by the Society in connection with, or in furtherance of the
- classes ships in accordance with its Rules and performance of its services, is governed by Rules of the Society
publishes the ship's names in its Register, which shall be the sole judge thereon, in its absolute discretion.
- carries out periodic surveys of ships in service in order to
ensure that the requirements for the specified class are ARTICLE 6
met. It is the owner's duty to report to the Society any 6.1. In providing services, information or advice, neither
incident or event affecting the condition of the ship or the the Society nor any of its servants or agents warrants the
need for repairs, whatever their nature or extent. accuracy of any information or advice supplied.
Furthermore, all express and implied warranties are
ARTICLE 4 specifically disclaimed. Except as set out herein, and also
4.1. The Rules, Guidance Notes and other documents in the case of the surveys mentioned in Article 1.2., neither
published by the Society endeavour to meet the state of currently the Society, nor any of its servants or agents, shall be liable
available technology at the time they are published. for any loss, damage or expense of whatever nature
Committees, consisting of personalities from the International sustained by any person, in tort or in contract, due to any
Maritime Industry, contribute to the development of these act or omission or error of whatever nature, whether or not
documents. negligent, and howsoever caused in any information or
The Society is not responsible for any inadequacy or failure advice given in any way whatsoever by or on behalf of the
of such Rules, Guidance Notes or other documents as a result of Society.
future development of techniques which could not reasonably
have been foreseen at the time of their publications. 6.2. Nevertheless, if any person uses the Society's
services, and suffers loss, damage or expense thereby,
4.2. BUREAU VERITAS is sole qualified to decide upon the which is proved to have been due to any negligent act,
meaning, the interpretation and application of the documents omission or error of the Society, its servants or agents, or
mentioned above. Any reference to them is of no value if it does any negligent inaccuracy in information or advice given by
not involve, accompany or follow a Society's intervention. or on behalf of the Society, then the Society will pay
compensation to such person for his proved loss, up to,
but not exceeding, five times the amount of the fee - if any -
charged by the Society for the particular service, ARTICLE 8
information or advice and, if no fee is charged, FRF 50,000. 8.1. Any intervention of BUREAU VERITAS, whatever its
Where the fees are related to a number of services, the nature, whether completed or interrupted for any cause
amount of fees shall be apportioned for the purpose of the whatsoever, shall involve the payment of fees upon receipt of the
calculation of the maximum compensation, by reference to invoice and the reimbursement of the expenses incurred.
the estimated time involved in the performance of each Interests may be demanded in case of late payment.
service.
Any liability for consequential damages is specifically 8.2. The class of a ship may be withheld or, if granted,
excluded. suspended or withdrawn, in the event of non-payment of fees.
In any case, whatever the amount of the fees, the
maximum damages payable by the Society shall be not ARTICLE 9
more than FRF 5,000,000. 9.1. All plans, drawings, specifications and information
provided to the Society, or of which the Society may become
6.3. Any claim for damages of whatever nature by virtue reasonably aware in connection with the performance of its
of the provisions set forth herein must be made in writing, services, shall be treated as confidential by the Society and shall
and notice must be provided to the Society within three not, without the prior consent of the party furnishing such
months of the date the services were first supplied or the document, be used for any purpose other than that for which they
damages first discovered. Failure to provide such notice are furnished.
within the time set forth herein will constitute an absolute
bar to the pursuit of such claim against the Society. 9.2. Records and other documents whatsoever relating to
ships classed with the Society shall not be disclosed or
ARTICLE 7 communicated by the Society to any third party but with the prior
Requests for interventions shall, in principle, be submitted in approval of the owner, unless it is in pursuance of an order of the
writing. They entail the acceptance without reservation, of the flag authorities or an injunction of a Court having jurisdiction.
present General Conditions and of Particular Conditions which
may be attached thereto. ARTICLE 10
Should a technical disagreement arise between the
requesting party and the BUREAU VERITAS surveyor, the
Society may, at the request of that party, designate another of its
surveyors.

ARTICLE 11
11.1. Any dispute over the interpretation, or the
enforcement of these General Conditions except as regards
article 8, or the performance of BUREAU VERITAS'
services, shall be submitted to arbitration in London.
Arbitration proceedings will be conducted according to the
English law on arbitration. The arbitrator shall be selected
by the parties and, in case they could not reach an
agreement, he would then be nominated by the President of
the Law Society of England and Wales upon request of any
party.

11.2 Disputes arising out of the enforcement of the


provisions of the article 8 shall be submitted to the Court of
Nanterre - France.
CONTENTS
SECTION 1. GENERAL.....................................................................................................1
1.1. Introduction...................................................................................................................................... 1
1.2. Scope of the present document..................................................................................................... 1
1.3. Limitations of the calculation methodology..................................................................................2

SECTION 2. ANALYSIS PROCEDURE.............................................................................3


2.1. Basic principle.................................................................................................................................. 3
2.2. Static mooring line response...................................................................................................... 5
2.3. Initial tensions................................................................................................................................ 13
2.4. Modelling of the environment....................................................................................................... 18
2.5. External loads contributing to the low frequency response..................................................21
2.6. Wave frequency response.......................................................................................................... 30
2.7. Simulation procedure.................................................................................................................... 33

SECTION 3. CRITERIA....................................................................................................39
3.1. Environmental criteria.................................................................................................................. 39
3.2. Design criteria................................................................................................................................ 40
3.3. Fatigue life...................................................................................................................................... 41

APPENDIX A SENSITIVITY TO INPUT PARAMETERS..................................................42


A.1 Account for mooring stiffness in RAO's calculations................................................................42
A.2 Linear damping matrix................................................................................................................... 42
A.3 Instant of line failure...................................................................................................................... 42

APPENDIX B TYPICAL VESSEL DATA..........................................................................43


B.1 Generic tankers.................................................................................................................................. 44
B.2 Semi-submersible units.................................................................................................................... 80
B.3 NKP barge........................................................................................................................................... 89

APPENDIX C CALCULATION EXAMPLES.....................................................................92


C.1 Fitting of the material elastic curve by a 5th order polynomial.................................................92
C.2 Characteristic equivalent to a CALM............................................................................................92
C.3 Characteristic equivalent to a fender...........................................................................................92

Bureau Veritas Guidance Note for the quasi-dynamic calculation of mooring systems using ARIANE sofware
Section 1 GENERAL

SECTION 1 GENERAL

2 Introduction 3 Scope of the present


document
The assessment of environmental loads on a
floating body and the derivation of design
mooring loads are generally based on an 4
engineering approach combining both model tests
and numerical simulations. Provided they are This document aims at providing:
properly carried out, model tests allow for the
a complete methodology for performing
accurate consideration of a large number of
numerical simulations (Section 12),
parameters such as the influence of the current,
the effect of the floating body velocity on the design requirements consistent with the
wave drift forces, the wave drift damping, the line above methodology for evaluating the
dynamics and, in some instances, the effects of fitness-for-purpose of catenary mooring
risers. lines and anchors in extreme and fatigue
conditions (Section 169),
However, model tests remain expensive and acceptance criteria as regards some other
the number of combinations of environmental and aspects such as minimum clearance to
vessel loading conditions that it is feasible to test prevent risks of interferences.
is inevitably small compared to the number of
actual possible combinations. Hence, numerical In other words, environmental and design
simulations appear to be the best way to criteria together with associated safety factors
complete the database derived from model tests. and fatigue criteria that are specified in Section
169 have been calibrated with respect to the
Whatever the method used, the evaluation of a calculation procedure defined in Section 12
station keeping system either permanent, such as
for a floating production system, or temporary, The attention is drawn on the fact that other
such as for a drilling rig, must address the approaches may therefore require the definition
following topics: of different safety factors.
definition of design loads in mooring
system components (e.g. mooring lines, 5
anchors, connecting hardware on the
moored vessel) with the aim of deriving Operational procedures for installation,
their minimum acceptable scantling anchoring line deployment, desinstallation,
definition of design motions with a view to inspection and maintenance, etc. are not covered
assessing risks of interference between herein.
mooring lines and moored vessel with
adjacent structures or to deriving relevant 6
data for the design of production or
operational equipment (e.g. connection The integrity of risers, if any, connected to the
bridges, risers, etc.) moored vessel, is not addressed. Neither is the
definition of design accelerations of the verification of bridges or any sub-structure or
moored vessel equipment the design of which may be governed
derivation of long term distribution of by the response of the mooring system, nor the
verification of the resistance of the connecting
aforesaid parameters in order to address
hardware on the vessel side, e.g. turret, fairleads,
risks of fatigue failure
winches, windlasses, brakes, etc.
In the above, the word design refers to both
extreme and service conditions in intact and
damage configuration of the station keeping
system.

Bureau Veritas Guidance Note for the quasi-dynamic calculation of mooring systems using ARIANE sofware
1
Section 1 GENERAL

Also is excluded from the scope of the present


It is supposed that horizontal low and wave
document any design recommendation or
frequency phenomena do not interfere.
performance consideration with regard to the
Compliance with this assumption is reasonably
components of the mooring line. In this respect,
satisfied if the natural period of the mooring
the safety factors that are given in 183 are
system in surge, sway and yaw is greater than
applicable to the relevant resistance
five times the zero-up crossing period of the
characteristics mentioned in the suppliers
wave.
catalogue, e.g. break strength, holding power,
etc., and to mooring line components in good
condition utilised in the way they are designed for. 10
The adequacy of the type of anchoring point The variation of the suspended line weight with
selected, e.g. drag anchor, driven pile, suction the motion of the moored vessel is supposed not
caisson, etc., with the soil capacity is not to significantly modify the average vessel
discussed. draught, trim or list angles.

11
7 Limitations of the
calculation methodology Assumptions set out in 8 to 10 are found to be
satisfied for conventional catenary mooring
systems operating in water depths not exceeding
8 450m.

The mooring system is assumed not to be


subject to resonance at the wave frequency. In
addition, out-of-horizontal-plane low frequency
motions are supposed to be negligible. Hence,
this methodology may not be appropriate to
Tension Leg Platforms or certain types of semi-
submersibles operating at very large draughts. 9

Bureau Veritas Guidance Note for the quasi-dynamic calculation of mooring systems using ARIANE sofware
2
Section 2 ANALYSIS PROCEDURE

SECTION 12 ANALYSIS PROCEDURE

13 Basic principle is the origin of the fixed axis system,


X is positive northward,
14 General Y is positive eastward,
Zis positive downward,
15
XY is in the plane of the still water
The methodology is based on time domain level.
simulations of the mooring system response in
specified environmental conditions. 20 Vessel axis system

Considering the limitations mentioned in 7, the The vessel axis system Oxyz is defined as
calculation process does not address the line follows:
dynamics. The safety factors given in 183 have O is the origin of the vessel axis
been calibrated accordingly. The mooring lines system.
are subsequently considered as non linear O is on the keel line of the vessel which
massless springs which can freely follow their is assumed to be longitudinally
fairleads without inducing to the vessel any other symmetric.
loads than their static spring reactions in the
Oxis positive forward,
instantaneous vertical plane of the mooring line.
Oyis positive rightward,
16 Oz is positive downward.
The calculation procedure consists of the
Throughout this document, Oxy remains
determination of the low frequency response of
horizontal in the low frequency motion of the
the moored vessel under the effect of waves,
vessel.
wind and current, followed by the superimposition
of the wave frequency motions. It is assumed that
The position of the vessel in the fixed axis
low and wave frequency components do not
system is identified by the three co-ordinates XO
significantly interfere with each other because of
YO and ZO of the origin O and by the angle that
very different time scales. As a consequence,
the projection of Ox onto XY makes with X.
they are assessed separately in the framework of
this approximation and added together at the end The angle is called the vessel heading. The co-
of each time step of the simulation. ordinates XG, YG and ZG of its centre of gravity
may also be used.
17
See Figure 20
The mooring line tension is derived at the end
of each time step from the static catenary The position of the vessel centre of gravity
response obtained for the instantaneous position within the vessel axis system is identified by its
of the fairlead. three co-ordinates x G, yG and zG. The vessel
centre of gravity is assumed to be in the plane
Oxz.
18 Axis systems

19 Fixed axis system


Throughout this document, the position of any
object is identified by its co-ordinates in a fixed
axis system XYZ which is defined as follows:
Section 2 ANALYSIS PROCEDURE

FIGURE 2.2.2.
AXIS SYSTEMS m 0 0
M 0 m 0

0 0 I
X

The values of m and I should be
derived from the lightship data and the
actual loading conditions of the vessel.
XG G x
{F(t)} is the three-component vector of
XO O
y those horizontal loads, applied to the
xG centre of gravity of the vessel at instant
t, which are contributing to the low
frequency response:
YO YG Y
FX (t)

21 Mooring line F(t) FY (t)
M
The mooring line position is identified by the / G (t)
co-ordinates of its anchoring point in the fixed
axis system and by the angle that the vertical 24
plane passing by both the anchor and the fairlead The breakdown and the way of computing F X(t),
makes with ZX. This angle is called the azimuth FY(t) and M /G(t) are detailed in 100
of the mooring line.

Within this vertical plane, any point of the 25 Wave frequency response
mooring line is identified by its horizontal and
vertical co-ordinates with respect to the anchor. 26
It is assumed that the wave frequency motions
22 Low frequency response of the vessel are not significantly disturbed by the
variation of the mooring stiffness with the low
23 frequency offset. An average mooring stiffness
can therefore be used for pre-determining the
The low frequency component is obtained by Response Amplitude Operators (RAO's) of the
solving at each time step a vectorial differential vessel.
equation of the following form:
27
M X
F(t)
At each time step, the six wave frequency
motions of the vessel centre of gravity are added
Where:
to its low frequency position. To do so, the
{X}is the three-component vector amplitude of each component of the wave signal
characterising the horizontal position of is multiplied by the RAO's of the centre of gravity
the vessel centre of gravity G, i.e. the of the vessel and the summation is carried out
vessel heading and the co-ordinates with due account for time and space phases as
XG and YG in the fixed axis system. detailed in 131
[M] is the horizontal mass matrix of
the vessel calculated at its centre of 28
gravity. It is a 3-by-3 diagonal matrix the Equations of rigid body motions lead to the
components of which are the actual instantaneous position of each fairlead. The
mass of the vessel for surge and sway tension in the mooring line is then computed by
motions (m), and the moment of inertia interpolation of the static catenary responses pre-
in yaw (I): calculated for different vertical positions of the
fairlead as explained in 29
Section 2 ANALYSIS PROCEDURE

29 Static mooring line modelled by a series of straight segments the


number of which depends on the complexity of
response the actual shape.

30 Main assumptions The same profile is used whatever the fairlead


position. This implies that the seabed shape is
31 locally approximated by a generalised cone
sector the symmetry axis of which is the vertical
The assumptions set out in 32 to 36, together passing by the anchoring point as shown on
with those specified in Section 1 and 13, make Figure 36
that the load induced by any mooring line to the
moored vessel in its low frequency motion This approach is rigorously exact if the seabed
depends only on the anchor-to-fairlead distance. is a horizontal plane and remains very accurate
for moderate variations of the seabed shape. It
The mooring line response can therefore be remains very accurate whatever the seabed
pre-computed using the equations detailed in 42 shape if the fairlead motions out of the initial
to 68, to build up the mooring line characteristic horizontal plane of the mooring line are small in
as defined in 37. This characteristic can next be comparison with the distance between the anchor
used throughout the simulation, whatever the and the fairlead. This latter condition is usually
position taken by the moored vessel. verified.

32
The upstream end of the mooring line is FIGURE 36
anchored to a fixed point which is located on or SEABED PROFILE
above the seabed. The downstream end is
connected to a fairlead of the moored vessel
which can be either emerging or immersed. The
mooring line cannot penetrate the seabed.

33
The mooring line is made up of a series of
homogeneous segments attached end-to-end, the
bending stiffness of which is negligible. A
homogeneous segment is characterised by
constant mechanical properties over its whole
length. A buoy or a sinker may be connected at
the downstream end of any segment.

34
At any time, the mooring line is assumed to be
in the vertical plane passing by its anchoring point
and its fairlead. This implies that wave, current,
wind and dynamic loads on any of the mooring
line components are neglected. It also implies
37 Mooring line characteristics
that friction effects, transverse to those parts of
the line laying on the seabed, are not taken into 38
account. The characteristic of a mooring line is the
curve which gives, for a given paid-out length and
35 a fixed altitude of the fairlead, the horizontal
component of the tension at the fairlead as a
Friction effects in the plane of the line tend to
function of the horizontal distance between the
decrease the tension towards the anchor. Friction
fairlead and the anchor.
loads are proportional to the perpendicular
component of the soil reaction.

36
The profile of the seabed is determined in the
initial vertical plane of the mooring line. It is
Section 2 ANALYSIS PROCEDURE

The tension-offset curve of a mooring line is length of the equations developed in 42 to 68 with
the curve which gives, for a given paid-out length due account for the applicable boundary
and a fixed altitude of the fairlead, the axial conditions.
tension at the fairlead as a function of the
horizontal distance between the fairlead and the Figure 38-II is an example of characteristic and
anchor. tension-offset curve obtained for the line
described in Figure 38-I.
The characteristic and the tension-offset curve
are derived from the integration over the full line

FIGURE -I
COMPOSITE MOORING LINE ON A COMPLEX SEABED PROFILE

-25,17662
-25,294 Chain 3" Gr3 Tension 910 kN
500
-25,41009 (1410 m) Distance 1536 m
-25,52492
400
-25,63849
-25,75083
-25,86195
300 Chain 3" Gr3
-25,97187 (450 m) Buoy (30t)
-26,0806
200
-26,18816
-26,29456
-26,39982
100
-26,50396
-26,60698
-26,796160
-26,79616
-47,88858
-100
-68,02418
Anchor Pipeline
-68,02418
-200
-68,25404
0 200 400 600 800 1000 1200 1400 1600 1800 2000
-68,47888
-68,69885
Section 2 ANALYSIS PROCEDURE

FIGURE 38-II
MOORING LINE CHARACTERISTIC AND TENSION-OFFSET CURVE

2550
2350
2150
Tension at fairlead kN)

1950 Horizontal component Axial tension at


of the tension at fairlead fairlead
1750
1550
1350
1150
950
750
1525 1575 1625 1675 1725 1775
Distance anchor-fairlead (m)

39 As leeward lines are liable to slacken, the


threshold distance under which the horizontal
For the determination of the low frequency
component of the tension at fairlead is equal to
response of the moored vessel, the characteristic
zero, should be determined in as accurate a
of each mooring line should be calculated for the
manner as possible. In this respect, catenary
altitude of the fairlead corresponding to the
equations cannot be used. One means is to
hydrostatic equilibrium of the vessel in still water.
compute the characteristic for two distances a bit
Calculations should be carried out for the full
bigger than the threshold and to extrapolate
range of anchor-to-fairlead distances that the
downwards as shown on Figure 39
fairlead is expected to meet during the simulated
sea state.
Section 2 ANALYSIS PROCEDURE

FIGURE
DETERMINATION OF THE ZERO OF THE CHARACTERISTIC BY LINEAR EXTRAPOLATION

40 Besides, environmental loads may not be


negligible on such systems. Suitable provisions
The tension-offset curve of each mooring line
are subsequently to be taken to adequately
should be calculated for at least three altitudes of
transfer such environmental loads to the vessel
the fairlead with a view to determining by
and thus ensure that the methodology is still
interpolation the tensions due to the combined
applicable. Possible aerodynamic and
low and wave frequency motions of the vessel.
hydrodynamic interferences may however
The three altitudes should be selected so as to
invalidate the whole calculation process.
appropriately cover the range of vertical positions
that the fairlead is expected to meet during the
simulated sea state. 42 Material elastic behaviour
41 43
Certain mooring systems such as SALM's or The elastic curve of the material should be
CALM's can be represented by an equivalent approximated by a polynomial function of the
characteristic for the purpose of the low frequency tension. If an elementary length dl of the material
response analysis. However, the restoring force of at rest is thus submitted to a tension T applied to
such a system may also vary with the direction of both ends, the length becomes dl+ (dl), is
the loads transmitted by the vessel attached to it. proportional to the initial length dl and increases
In such a case, several directional characteristic with T so that a relationship of the following form
curves should be determined in order to cover the can be written:
range of all load directions expected during the
simulated sea state. Interpolations on the load
direction are therefore to be carried out at each
time step of the simulation to derive the
instantaneous characteristic.
Section 2 ANALYSIS PROCEDURE

k n In practice, the function , which is worth 0

dl (dl) dl k T k1
when equals 0, should first be approximated by
a polynomial in the following way:

k 1 k n1
The coefficient 1 in the above summation is
equal to 1 since (dl) is equal to 0 if no tension is
k k
applied. k 1

The minimum degree of the polynomial and and the coefficients (k)1kn defined in 43
therefore the value of n depends on the nature of should then be derived from the following
the material as specified in 44 and 45 relations:

44 1 1
For mooring chains and wire ropes, n is equal
to 2 and 2 is given by the following equation: k 1
k k 1
for 2 k n
4 BL
2
E 2
An easy method to determine the coefficients
of a polynomial that approximates a given curve
where is the nominal diameter of the chain crossing 0 is presented in Appendix C.
or wire rope and E the equivalent Young modulus.
Without further information, the values of E given
by Table 44 can be used. 46 Suspended segments (catenary
equations)
TABLE 44 47
EQUIVALENT YOUNG MODULUS 48 to 50 give the relations between the main
OF WIRE ROPES AND MOORING CHAINS variables at both ends of a homogeneous
Mooring chain Six-strand wire Spiral-strand segment when fully suspended in a uniform fluid.
rope wire rope
A fully suspended homogeneous segment is
E (kN/mm 2) 112 98 152
characterised by the following data:
l its at-rest length,
p its apparent weight per unit length in
45 the fluid,
For synthetic materials such as nylon, (k) 1 k n the coefficients
polypropylene, kevlar, etc., the value of n should as defined in 42 of the polynomial used
not be taken less than 4. A higher value may to approximate the elastic curve of the
however be needed to fit the actual elastic curve material that it is made of.
over the complete range of tensions that the
mooring line is expected to encounter. At any point along the segment, the main
variables are:
It should be noted that the elastic curve of
t the tangent of the angle with the
synthetic ropes is usually presented by the
horizontal,
manufacturers as a curve giving the relative
elongation of the material as a function of the T the axial tension,
percentage of the breaking load (BL): x the horizontal co-ordinate,
y the vertical co-ordinate.
dl
where T B L
dl Subscript u is used to identify the variables at
the upstream end of the segment (e.g. Tu) while
subscript d is used for the downstream end (e.g.
Td).
Section 2 ANALYSIS PROCEDURE

48 52 Laid down segments (taut


The relations between the tangents and the thread equations)
axial tensions at both ends of the suspended
segment are: 53
54 to 56 give the relations between the main
pl 2 taut thread variables at both ends of a
t d tu 1 tu
Tu homogeneous segment which runs the length of a
rectilinear support that makes an angle with the
horizontal.
Td Tu

2 2
1 td 1 tu In addition to the data specified in 47, a taut
thread is characterised by the coefficient of
friction with its support. This coefficient is noted
49
in the following and considered as the tangent of
The relation between the horizontal co- an angle , so that tan. Without further
ordinates at both ends of the suspended segment information, can be given the value of 1 for
k mooring chains ( 45 deg) and 0.7 for wire
k n
k Tu ropes ( 35 deg).

d u
is: x x
k 1 p 1 t
2

Fxk t d Fxk tu where the Fx k
At any point along the segment, the main

u variables are the same as those defined for


catenary equations in 47
functions are defined as follows:
54
Fx1 t ln t 1 t2 The relations between the tangents and the
axial tensions at both ends of the laid-down
Fx 2 t t segment are:

1 k2 t d t u tan
t 1 t 2
(k 2) 2
Fx k t Fx t
k 1 k 1 k2
sin
for k greater than 2 Td Tu p l
cos
50
55
The relation between the vertical co-ordinates
at both ends of the suspended segment is: The relation between the horizontal co-
ordinates at both ends of the laid-down segment
k is:
k n
k Tu
y d yu
Fyk t d Fyk tu where cos cos kn k k k

the Fyk

k 1 p 1 t
2
x d xu Td Tu
u
functions are defined as follows: p sin k1 k
Fyk t
1

k 2
1 t2 56
k
The relation between the vertical co-ordinates
51 at both ends of the laid-down segment is:

A suspended segment that crosses the free


cos sin kn k k k

surface obeys the catenary equations on both
sides. All variables are continuous at the crossing
yd yu Td Tu
p sin k1 k
point, except p.
Section 2 ANALYSIS PROCEDURE

57
x x
If the slope of the rectilinear support balances
the friction coefficient (case where 0), the y y
formulae given in 55 and 56 become:
T T
k n

2 2

xd xu l cos k Tk1
1 t 1 t

T t T t
k 1 P

2 2
1 t 1 t

k n
y d yu l sin k T
Where:
k 1
P is the net action of the buoy or sinker on
k 1 the point of the line where it is attached
the minus sign as an exponent (e.g. T)
58 characterises the value of the variable just
before the point where the buoy or sinker is
A segment a part of which is laid down while attached,
the rest is suspended obeys the taut thread
equations for its laid down part and the catenary the plus sign as an exponent (e.g. T)
equations for its lifted part. All variables are characterises the value of the variable just
continuous at the junction. after the point where the buoy or sinker is
attached,

59 Buoyant segments (p0) 65


Buoys and sinker should be carefully modelled
60 so that the net action that they exert to the
A buoyant segment is a segment whose mooring line remains correct whatever the
apparent weight per unit length in water is tension in the line. Indeed, a subsea buoy initially
negative. If fully immersed, it takes the shape of immersed may emerge with the part of the line
a reverse catenary. If released, it floats as a where it is attached when the tension increases
straight line at the free surface This is the case of and then behave as a pure weight. In the same
certain synthetic hawsers. way, a sinker made up of heavy chains may fully
or partly lay on the seabed for low tensions and
61 therefore provide a lesser action on the line.

By using the appropriate value of p, equations 66


of 46 apply to such a segment, whether it is fully
suspended in air or fully immersed in water. As an illustration, Figures 66-I and 66-II
represent a schematic view of a buoy and its
62 associated curve which gives the net action
(weight minus buoyancy) as a function of the
By taking 0, equations of 54 and 57 vertical distance between the seabed and the
apply to the buoyant segment when floating at the point of the line to which the buoy is attached by
free surface. means of a pennant line.

63 Buoys and sinkers

64
A Buoy (respectively a sinker) induces an
upward (respectively downward) load to the point
of the line to which it is attached so that the
tangent is discontinuous at such a point. The
relations between the main catenary variables
just before and just after the buoy or sinker are as
follows:
Section 2 ANALYSIS PROCEDURE

FIGURE 66-I FIGURE 67-I


CYLINDRICAL BUOY WITH PENNANT LINE SINKER MADE UP OF HEAVY CHAINS

Cylindrical buoy
Total Buoyancy: B
Total Weight : W
B>W

Mooring line
Heavy chain sinker
Total Weight : W
Pennant Total Buoyancy : B
W>B

H h
Mooring
line H

Seabed

FIGURE 67-II

FIGURE 66-II NET ACTION OF THE SINKER VERSUS


WATER DEPTH
NET ACTION OF THE BUOY VERSUS
WATER DEPTH

P >H
P
rH/
H H H
0
rB/W
h

BW r(WB)

Hr H
0
H Hr h WrB
W

68 Seabed shape
67
In a similar manner, Figure 67-II represents the 69
net action of a sinker as a function of the vertical The seabed shape in the initial vertical plane of
distance between the seabed and the point of the the mooring line should be approximated by a
line to which it is attached. Figure 67-I is a sufficient number of straight segments. Using a
schematic view of the system. single segment obtained by average or least
mean square methods may not be accurate
enough. Figures 69-I and 69-II show both an
acceptable and an unacceptable modelling of the
same seabed shape. Figure 69-III clearly
demonstrates the sensitivity of the line
characteristic to the seabed modelling.
Section 2 ANALYSIS PROCEDURE

FIGURE 69-I FIGURE 69-III


Line elevation using a seabed shape SENSITIVITY OF THE CATENARY
approximated by the average plane RESPONSE TO THE SEABED MOORING
UNACCEPTABLE
Sensitivity of the catenary response to the seabed modelling

Line elevation using a seabed shape 1200


approximated by the average plane

Horizontal restoring force (kN)


UNACCEPTABLE 1000
0
Precise seabed profile
Fairlead 800
-20 Average seabed profile
Actual seabed shape
600
Water depth (m)

-40
400
-60

Anchor 200
-80
0
-100 575 580 585 590 595 600 605 610 615

Distance anchor to fairlead (m)


-120
0 100 200 300 400 500 600

Distance anchor to fairlead (m)


70
Four possible configurations may occur where
the mooring line crosses a seabed peak. They are
FIGURE 69- II all represented on Figure 70 In any case, the
Line elevation using a precise seabed conditions of continuity are as follows:
description ACCEPTABLE
x x

Line elevation using a precise seabed description y y


ACCEPTABLE Fairlead
0
T T
-20
Water depth (m)

-40
Where:
Anchor
-60
the minus sign as an exponent (e.g. T)
-80 characterises the value of the variable just
-100 before the seabed peak,
-120 the plus sign as an exponent (e.g. T)
0 100 200 300 400 500 600
characterises the value of the variable just
Distance anchor to fairlead (m)
after the seabed peak.

FIGURE 70
MOORING LINE CROSSING A SEABED
PEAK

CASE 1 CASE 2

CASE 3 CASE 4
Section 2 ANALYSIS PROCEDURE

71 L,D,T 0
Unless friction on both sides of a seabed
where:
trough perfectly balances the tension, the
mooring line flies over the trough according to L is the paid-out length of the line,
catenary equations developed in 46 Both ends of D is the horizontal distance between
the catenary are tangent to the seabed. the anchor and the fairlead,
T is the tension at fairlead.
FIGURE 71
In practice, an iterative process is needed in
MOORING LINE CROSSING A SEABED most cases to ensure that the three parameters
TROUGH are compatible together.

As an illustration, let us suppose that the


anchor and vessel positions are known. The
distance D can therefore be calculated. If the
initial tension T is also part of the data, the next
step is to determine the line length L to deploy.
This can be done by building up the tension-offset
curve as defined in 37 for various paid out
lengths. The right length is that for which the
72 Initial tensions curve passes simultaneously by D and T as
shown on Figure 75
73 Mooring pattern
Another situation may happen where the line
74 length L and the initial tension T are known
whereas the anchor position is not determined. If
The mooring pattern is the theoretical the seabed is flat, the tension-offset curve is the
description of the mooring system as installed on same wherever the anchor is. The curve can
site in the absence of any external load. The therefore be built up first and used afterwards to
mooring pattern thus includes the general layout locate the anchor. Otherwise, iterations should be
of the mooring system (plan view) and the performed because the initial plane of the line
elevation of each line in its initial vertical plane. and subsequently the seabed profile depend on
To achieve such a description, the following the anchor position.
information is needed:
Site bathymetry,
FIGURE 75
Horizontal co-ordinates of the centre of the
vessel axis system (XO,YO) or vessel DETERMINATION OF THE LINE LENGTH
centre of gravity (XG,YG) in the fixed axis WHEN THE TENSION AND THE ANCHOR-
system, TO-FAIRLEAD DISTANCE ARE KNOWN
Vessel heading (),
Vessel draught, KG,
Fairlead co-ordinates in the vessel axis
system,
Anchor co-ordinates in the fixed axis
system,
Mooring line composition,
Paid-out line lengths.

75
The knowledge of all items listed in 74
automatically settles the initial tensions since a
relation of the following form exists for each line
as soon as its vertical plane is known:
Section 2 ANALYSIS PROCEDURE

76 Optimisation of initial tensions mooring system in the prevailing storm direction.


The maximum restoring force of the mooring
77 system in a given direction is that external force
that the system can withstand in this direction
The methodology does not provide for without having any line tension exceeding the
modifications of the line lengths during the storm. maximum allowable, i.e. the breaking load
divided by the applicable safety factor as defined
In return, the line lengths can be adjusted, in 183
during the installation and deployment of the
mooring lines or later on during calm weather The adjustment can be obtained by plotting the
periods, in order to optimise their initial tensions restoring force curve in such direction (external
and improve the overall response of the mooring load versus vessel offset) for various sets of initial
system. 78 and 79 give basic optimisation tensions and by marking on each curve the point
criteria. where the most critical line has reached its
maximum allowable tension. The best set of initial
78 tensions is that which leads to the greatest
The initial tensions can be adjusted so as to restoring force at this point.
optimise the maximum restoring force of the
See Figure 78

FIGURE
OPTIMIZATION OF THE INITIAL TENSIONS TO OBTAIN
THE MAXIMUM RESTORING FORCE IN A GIVEN DIRECTION
Section 2 ANALYSIS PROCEDURE

FIGURE 79
OPTIMIZATION OF THE INITIAL TENSIONS TO MINIMIZE THE TENSION OF
THE MOST LOADED LINE IN SPECIFIC ENVIRONMENTAL CONDITIONS

79 Because a mooring system does not respond


linearly, their is no universal method for
In particular conditions, some mooring systems
optimising its initial tensions. In any case,
are very little sensitive to slow drift effects. The
whatever the criteria used, the benefit of the
response of such a system is therefore that of the
optimisation should be verified for all
wave frequency superimposed on the mean offset
environmental conditions that the mooring system
due to environmental loads. The initial tensions
is likely to encounter.
can then be adjusted for the prevailing storm
direction so as to minimise the tension in the
most critical line under the combination of the 81 Mooring stiffness
mean offset and maximum first order motion.
82
The adjustment can be obtained by plotting for
various sets of initial tensions, the tension of the The mooring stiffness is the 6-by-6 matrix
most critical line as a function of the vessel offset which links elementary external loads applied to
in the storm direction and by marking on each the vessel with its resulting elementary
curve the point corresponding to the summation displacements around and about a given position.
of mean and maximum wave frequency offsets.
The best set of initial tension is that which leads In practice, the stiffness induced by the
to the lowest tension in the most critical line. mooring system for out-of-horizontal plane
motions is normally negligible in comparison with
See Figure 79 the hydrostatic stiffness.

80 The mooring stiffness matrix can subsequently


be reduced to a 3-by-3 matrix as follows:
Examples of 78 and 79 show that different
criteria may lead to different optimisations. FX K XX K X
K XY X
Indeed, the same mooring system and the same
sets of initial tensions have been used in both FY K YX K YY K Y Y
cases for a given storm direction. However, the M K X K Y K
results are not similar at all: where the maximum
restoring force is searched, the initial tension set K
N4 is the best while N2 is the best for specific
environmental conditions. where:
Section 2 ANALYSIS PROCEDURE

[K] is the mooring stiffness matrix the is the elementary rotation of the
coefficients of which depend on the co- vessel about the centre of its axis
ordinates XO and YO of the centre of the system
vessel axis system and on the vessel
heading 83
FX is the component of the The matrix [K] is the summation of the stiffness
elementary load parallel to X applied matrices of all mooring lines:
to the centre O of the vessel axis
NL

K K
system
i
FY is the component of the
i 1
elementary load parallel to Y applied
to the centre of the vessel axis system
where:
M is the elementary yaw moment
[Ki]is the stiffness matrix of mooring line
applied to the centre of the vessel axis
number i
system
NL is the total number of mooring lines
X is the elementary displacement
parallel to X of the centre of the vessel The coefficients of the stiffness matrix of each
axis system mooring line are determined by the following
Y is the elementary displacement equations:
parallel to Y of the centre of the vessel
axis system

TH i 2
K iXX sin A i k i cos2 A i
Di
TH i
K iXY k i sinA i cosA i
Di
TH i
K iX R i sinA i cos A i i k i cosA i sin A i i
Di
TH i
K iYY cos2 A i k i sin 2 A i
Di
TH i
K iY R i cosA i cos A i i k i sinA i sin A i i
Di
R
K i R i TH i 1 i cos A i a i cos A i i k i R i sin 2 A i i
Di

K iX K iX Di is the horizontal distance between


the anchor and the fairlead of mooring
K iY K iY line number i
K iYX K iXY Ai is the azimuth of mooring line
number i as defined in Section 1
where: i is the azimuth characterising the
THi is the horizontal component of position of the fairlead of line number i
the tension applied at the fairlead of i atan(yi/xi) where is the
mooring line number i vessel heading and x i and yi are the co-
ordinates in the vessel axis system of
the fairlead of line number i
Section 2 ANALYSIS PROCEDURE

Ri is the horizontal distance between cos sin 0


the centre O of the vessel axis system
P sin cos 0
and the fairlead of mooring line number O
i 0 0 1
where is the vessel heading.
R i xi 2 y i 2
Hence:
ki is the tangent of the characteristic
curve of line number i at point (THi, Di) cos sin 0
0
P
T
sin cos
dTH i O
ki 0 0 1
dD i
1) Stiffness matrix calculated in the vessel
See Figure 83 axis system at the vessel centre of gravity
G:

K P K P
T
FIGURE 83 G G G

TANGENT TO THE CHARACTERISTIC The matrix P is defined as


G

follows:
cos sin 0 1 0 0
P sin cos 0 0 1 x G
G
0 0 1 0 0 1

where x G is the algebraic distance of the


vessel centre of gravity G from the origin O of the
vessel axis system.

Hence:

1 0 0 cos sin 0
P 0 1 0 sin 0
T
cos
G
84 0 x G 1 0 0 1
If drift loads on the vessel and wave frequency
motions are derived from diffraction radiation
analysis, the mooring stiffness of the system in its
average position should be determined and used
85 Modelling of the
in such analysis as specified in 100 and 131 environment
However, most existing diffraction radiation 86 General
computer programs use the stiffness matrix
K O or K G which is calculated in the vessel Global characteristics of wave, wind and
axis system at either the origin O of the vessel current are assumed to remain constant during a
axis system or at the vessel centre of gravity G. storm as specified in Sections 87 to 96
Such a matrix can be derived from the above
matrix K as explained hereafter:
87 Wave
1) Stiffness matrix calculated in the vessel
axis system at the origin O: 88

K P K P The wave is characterised by its incidence and


T
O O O
an energy spectrum appropriate to the mooring
The matrix P is defined as
O
site. Spreading functions are not considered in
this approach.
follows:
The Pierson-Moskowitz and JONSWAP
formulations are generally used. The analytical
Section 2 ANALYSIS PROCEDURE

formulae of their energy density function are K is a factor such as


given in 89 and 90 HS 2
S J d
Other spectra can be used if more relevant. 0 16
Where the mooring analysis aims at corroborating
the results of model tests, it is preferable to use a The peakedness function J is defined as
discretized formulation corresponding to the follows:
spectrum actually generated in the tank test.


2
p
exp 1

2 p
89

The energy density function of the Pierson-


Moskowitz spectrum is defined as follows: where:
is the peakedness factor, typically
2
H Z
4
1 Z4
SPM S
taken equal to 3,3 in the North Sea;
exp
4 5
4
is the peakedness scale, typically
taken in the North Sea as 1 0,07 for
where: p and as 2 0,09 for p.

SPM is the energy density function of Because the energy density function of the
the spectrum; JONSWAP spectrum has no analytical primitive
is the wave circular frequency; function, the factor K should be computed by
numerical integration within the range [m , M]
HS is the significant wave height;
such as:
Z is the circular frequency
corresponding to the zero-up crossing m p M
period Tz Z 2 TZ .
and
Using the peak period, the energy density can
also be expressed as follows: SJ m SJ p / 100
5 HS
2
p
4
p
4
and
SPM exp 125
.
16 5
4
SJ M SJ p / 100
where:
See Figure 90
p is the circular frequency
corresponding to the peak period Tp FIGURE 90
p 2 Tp TRUNCATION OF THE JONSWAP
SPECTRUM FOR THE CALCULATION OF
90 THE K FACTOR
The energy density function of the JONSWAP
spectrum is defined as follows:

S J K SPM J

where:
SPM is the energy density function of
the Pierson-Moskowitz spectrum as
defined in 89;
J is the peakedness function of the
JONSWAP spectrum;
Section 2 ANALYSIS PROCEDURE

91 hundred Airy waves representative of the


spectrum:
The wave signal should be modelled by a
minimum of 100 elementary Airy wave 1) the energy density function of the spectrum
components representative of the whole is divided into 100 intervals of same length
spectrum. The instantaneous surface elevation at between m and M such as:
any point of co-ordinates (x,y) with respect to the
origin of the phases can thus be expressed as m p M
follows:

x, y, t
S m S M S p 100

M m 100
100
ai cos i t i k i x cos k i y sin
i 1
where p is the circular frequency
where: corresponding to the peak period
(x,y,t) is the surface elevation at ( p 2 Tp )
of the spectrum.
instant t at point of co-ordinates (x,y);
2) within each interval, a wave circular
ai i i are the amplitude, the frequency k is randomly selected.
circular frequency and the phase of the
elementary Airy wave number i; 3) to each k is attached an angle k
randomly drawn within [0,2].
ki is defined by the following formula in
which H is the water depth and g the 4) the amplitude ak of the wave component
gravity acceleration: number k is derived from the relation
between the wave amplitude and the
k i tanh(k iH) 2 g energy density function:
is the angle characterising the
ak 2 S k
direction where the wave comes from (it
is the angle of the direction of
propagation augmented by which
93 Wind
explains the positive signs in the above
formula).
94
See Figure 91 The wind is characterised by its incidence and
its 10-mn average speed at 10m above the mean
sea level.
FIGURE 91
DIRECTION OF WAVE PROPAGATION AND 95
SURFACE ELEVATION Where the wind speed is significantly
fluctuating and liable to generate a non negligible
low frequency response of the mooring system,
the wind speed should be characterised by its 1-
hour average value at 10m above the mean sea
level together with the reduced parameters of the
type of gust spectrum applicable to the site. For
tanker shaped units, the modelling of the wind by
means of a dynamic wind spectrum has been
shown not to significantly influence the results.
This conclusion may not be valid for large semi-
submersibles.

A method similar to that described for waves in


87 should be implemented to determine one
hundred or more representative wind
92
components. In return, it is not required to
The Montecarlos method can be used as account for spatial coherence so that the
follows to select a family (a k, k,k)1k100 of one
Section 2 ANALYSIS PROCEDURE

complete wind signal at any point can be written Oxy, all equations are projected on this latter.
as follows: Furthermore, all moments are computed at the
origin O of this system because its position within
100 the vessel does not depend on the loading
VW t V1h v k cos 2fk t k condition. The same axis system can thus be
used to record permanent data (e.g. fairlead co-
k 1 ordinates, lightship centre of gravity, etc.) as well
as other data which vary with the vessel or site
where: conditions (e.g. RAO's, load transfer functions,
etc.).
Vw(t) is the wind speed at instant t;
V1h is the wind speed at 10 metres Refer to 18 for the definition of the axis
above the still water level averaged systems.
over one hour;
(v k,f k,k)1k100 is the family of one
hundred wind components FIGURE 103
representative of the gust spectrum. ANGLES CHARACTERIZING THE
DIRECTIONS OF WAVE, WIND AND
96 Current CURRENT

97
The current is characterised by its incidence
and its average velocity over the draught of the
moored vessel.

98
As detailed in 114, the fluid velocity and
incidence relative to the vessel should be used to
derive current loads.

99
It is assumed that the influence of the current
on the wave kinematics is negligible and that
wave and current loads can be independently
calculated and simply added as detailed in 100
103
This assumption remains valid for current
velocities not exceeding 2.5 m/s. The wave, current and wind incidences are
respectively identified by the angles H , W and
C that the directions they come from make with
the fixed X-axis. Hence, a wave propagating
100 External loads contributing along the X-axis from the positive values towards
to the low frequency the negative values is characterised by an angle
response H equal to zero.

See Figure 103


101 General
104
102
The low frequency surge and sway of the
The vectorial differential equation presented in vessel are identified by the horizontal motions of
23 expresses the fundamental principle of the origin O of its axis system. The components
dynamics applied to the centre of gravity of the of the absolute speed of O projected on the
vessel in its low frequency motion relative to the vessel axis system are thus noted u and v
fixed axis system XY. (respectively surge and sway). The components
of the speed relative to the fluid at the same point
However, external load transfer functions being are noted u and v . These four parameters are
generally available in the vessel axis system
linked by the following relationships:
Section 2 ANALYSIS PROCEDURE

u u VC cos C I is as defined in 23;



v v VC sin C
xG is the algebraic distance of the
vessel centre of gravity G from the
origin O of the vessel axis system;
where: Fx is the projection on the vessel x-axis
VC is the absolute current velocity; (surge axis) of the external loads
applied to the vessel;
C is the current angle as defined in
Fy is the projection on the vessel y-axis
103;
(sway axis) of the external loads applied
is the vessel heading as defined in to the vessel;
18
M /O is the yaw moment of
105 external loads calculated at the origin of
the vessel axis system.
Using the variables of relative fluid velocity, the
equations of the slow drift motions projected on 106
the vessel axis system become:
External loads include hydrodynamic, mooring,
damping, wave drift, wind and current loads but
m u v 2 Fx
xG also other loads of various natures that are liable

to contribute to the low frequency response.
These are for instance the loads induced by
thrusters used to assist the mooring system or
those due to numerous risers.
I m xG
m x G v u
M /O
As previously said, all loads are projected on
where: the vessel axes, all moments being calculated at
the origin O, so that three equations of the
m is as defined in 23; following form can be written:

Fx FHx FM x FB x FD x FW x FC x FTx FRx

Fy FHy FM y FB y FD y FW y FC y FT y FRy

M / O M H / O M M / O M B / O M D / O M W / O M C / O M T / O M R / O

where the following subscripts are used to 107 Hydrodynamic loads


identify the origin of each term:
H for hydrodynamic loads, 108
M for mooring loads, The hydrodynamic loads are determined
according to the theory of manoeuvrability:
B for damping loads,
D for wave drift loads, FHx Ma xx u Ma yy v
Ma y
2
W for wind loads,
C for current loads, FHy Ma yy v Ma xxu
Ma y

T for thruster loads,
R for riser loads.
M H / O Ma x GMa y
Ma y x GMa yy v
The way of computing all above load Ma y
x GMa xx u
Ma yy Ma xx u v
components at each time step of the simulation is where:
detailed in the following sections.
Ma xx 0 0
Ma 0 Ma yy Ma y

0 Ma y Ma

Section 2 ANALYSIS PROCEDURE

is the asymptotic mass matrix calculated in the position of the vessel as close as possible to its
vessel axis system at the vessel centre of gravity average position during the storm.
G, so that:
112
Maxx is the asymptotic added mass in
surge; In practice, the diagonal terms of the QTF's are
Mayy is the asymptotic added mass in available for a limited number of incidences
sway; relative to the vessel heading. During the
simulations, the functions relevant to the
Ma is the asymptotic added mass in instantaneous incidence should be obtained by
yaw; interpolations between those of the closest two
May is the asymptotic added mass in incidences which are available. The step between
sway-yaw coupling. two consecutive incidences is therefore to be
limited to avoid invalid results.
109
Similarly, for a given incidence, the
The coefficients of the asymptotic added mass
discretization versus the wave circular frequency
matrix can be derived from model tests or 3-D
should be refined enough to correctly account for
diffraction-radiation analysis using a recognised
all peaks and troughs of the functions within an
computer program.
interval where the minimum value taken by the
energy density function of the wave spectrum is
The values of these coefficients are those
not greater than 1% of its maximum. Reference is
obtained when the wave circular frequency tends
made to 95
towards zero.
Discretization criteria are given in Table 112
110 Wave drift loads

111 TABLE 112


The slow drift loads are derived from the DISCRETIZATION CRITERIA OF QTF's
diagonal terms of the Quadratic Transfer DIAGONAL TERMS
Functions (QTF's) of the unit according to the
(1)
formulae specified in 113 TYPE OF UNIT

DISCRETIZATION Semi- Tanker Prismatic


The QTF's are the functions which give in the CRITERIA submersible unit shaped unit barge
vessel axis system the mean loads applied to the
centre of gravity G of the vessel when this latter Maximum
incidence step 25 deg 30 deg 30 deg
is subjected to the action of a bichromatic wave
of unitary amplitude. These loads are proportional
Number of
to the square of the wave amplitude. They incidences on 8 7 7
obviously vary with the two circular frequencies 180 deg
which compose the bichromatic wave but also
with the water depth and the wave incidence Number of
relative to the vessel. circular 40 25 25
frequencies
The diagonal terms of the QTF's are those
obtained when the two components of the (1) Examples of QTF-diagonal-term functions are given in
bichromatic wave are identical. They Appendix B for these types of units
subsequently correspond to the mean loads
induced to the vessel when subjected to an Airy 113
wave. They can be obtained by model tests or by
The slow drift loads are computed by a method
a second order diffraction-radiation analysis using
based on the Newmans approximation. The
a recognised computer program. In any case,
formulae used, however, involve four summations
they should be determined with due account for
instead of two in the original formulation. This
the actual site water depth and the mooring
allows the problem of indetermination of the load
stiffness. If linear calculation methods are used,
sign to be solved:
the mooring stiffness should be evaluated for a
Section 2 ANALYSIS PROCEDURE

100
FDG t ak QTFG H,k ,k cos k t k
k1
100

ak QTFG H,k ,k cos k t k sign QTFG H,k ,k
k1
100
ak QTFG H,k ,k sin k t k
k1
100

ak QTFG H,k ,k sin k t k sign QTFG H,k ,k
k1
where: M

FDG(t) is one of the three


FD MEAN / G 2
m
QTF
G
H , , S d

components in the vessel axis system


of slow drift loads applied to the vessel where m and M are defined in 90
centre of gravity G at instant t, i.e. F Dx ,
FDy or M D/G such as M D/O M D/G 114 Current loads
xG FDy
H is the wave incidence relative to the 115
vessel heading at instant t,
Current loads are usually pre-determined by
i.e. H H model tests. Such tests consist of measuring on
QTFGHkk is the relevant the vessel which is maintained in a fixed position,
diagonal function at the vessel centre of the loads generated by a constant current
gravity G, interpolated for the attacking from a given direction. The loads thus
instantaneous wave incidence H. measured are projected on the vessel axes in
order to establish a tabulation of each component
sign(u) is equal to: 1 if u>0
as a function of the current incidence. It is
-1 if u<0 assumed that all load components are
0 if u=0 proportional to the square of the current velocity.
(ak,k,k)1k100 is the family of one
If the vessel is moving, such tabulations can be
hundred Airy waves representative of
used to interpolate the instantaneous loads,
the wave spectrum as explained in 91
provided that the equivalent incidence and
The average value of FD G(t) on the whole current velocity as defined in 116 are applied. In
duration of the simulation can be obtained by the this respect, the tabulations given by: Prediction
following equation: of Wind and Current Loads on VLCC's, Oil
Companies International Marine Forum, 1994
can be used for tanker shaped units under the
same condition.
Section 2 ANALYSIS PROCEDURE

However, particular attention should be paid, FIGURE 116-II


when using such tabulations, to the reference
EQUIVALENT CURRENTS APPLIED TO A
point where the moments are calculated (origin O
of the vessel axis system, vessel centre of gravity
FIXED VESSEL
G, midship section, etc.).

116
The equivalent incidence of the current should
be determined according to the following
equation:

v v VC sin C
tan C
u u VC cos C

where C is the angle of the equivalent current


with respect to the heading of the moving vessel.
The other variables intervening in this equation
are defined in 104 Note that C is not equal to C -
unless the vessel speed is equal to 0 (u v
0).

The equivalent current velocity UC is given by


the following formula:

2 117
UC u v
If the vessel is fixed during model tests, the
Note that UC is not equal to V C unless the tabulations cannot include any effect due to a
vessel speed is equal to 0 (u v 0). possible rotation of the vessel in the fluid. The
Molins yaw moment is therefore to be added to
See Figures 116-I and 116-II. that derived from model tests performed with a
fixed vessel:

FIGURE 116-I
EQUIVALENT CURRENTS APPLIED TO A
MOVING VESSEL
Section 2 ANALYSIS PROCEDURE

xfore

MCMolin/O 21CMolinTE V(, ) V(, ) V(,) 0 V(,) 0 d
x a ft
where:
M C Molin /O is the additional yaw moment
V (, ) is the transverse component of
to be applied to the origin O of the the relative fluid velocity at the
vessel axis system; algebraic distance from O:

is the sea water density (typically V (, ) (v )


1025 kg/m 3);
TE is the vessel draught; V( , ) is the total relative fluid
velocity at the algebraic distance from
xaft is the algebraic distance of the aft
end of the vessel from the origin O of its O: V( ,
) u ( v )
axis system; C Molin is a coefficient to be determined
xfore is the algebraic distance of the as per Table 117
fore end of the vessel from the origin O
of its axis system;

TABLE
VALUE OF THE C Molin COEFFICIENT

TYPE OF UNIT

MOLINs COEFFICIENT Semi-submersible unit Ship shaped tanker Prismatic barge

C
FCy 2

0.3 C Cy 2

FCy 2
Molin

0.75 (2) (1)


0.75 (2)

1 1
LTE VC 2 0.4 C Cy

LTE VC 2
2 2 2
In this Table:

(1) CC y 2 is the value of the current coefficient in sway for beam incidence. This coefficient can be derived from: Prediction of Wind

and Current Loads on VLCC's, Oil Companies International Marine Forum, 1994, with due consideration for the reference point where the
moments are calculated;

FCy 2 VC 2

(2) should be derived from model tests or relevant literature;

is the water density;

L is the vessel length;

TE is the vessel draught.

118 included in the equation of the hydrodynamic


moment M H /O that is given in 108, should be
If tabulated yaw moments of current result from
substracted from current loads:
model tests, the Munks moment which is already


M C Munk Ma yy Ma xx VC sin C cos C
2
Section 2 ANALYSIS PROCEDURE

119 Wind loads 121


If a dynamic wind spectrum is used to generate
120 the fluctuation of the wind speed around its 1-
Same as for current, wind loads are hour average (V1h) as specified in 93, it should be
proportional to the square of the wind speed and noted that the average load derived from such a
are normally pre-determined by model tests. signal is not equal to the load computed with a
However, the wind speed is averagely far much constant wind speed equal to V 1h. Indeed, wind
greater than the speed of the moored vessel loads are proportional to the square of the speed,
throughout a storm. The wind loads can therefore the average of which verifies:
be computed at each time step of the simulation
by using the absolute wind speed and the actual 100
VW t V1h v
wind incidence relative to the vessel heading, 2 2 1 2
excepted where the fluctuating part of the wind 2 k
speed is significant. k 1

For tanker shaped units, the tabulations given


by: Prediction of Wind and Current Loads on 122 Mooring loads
VLCC's, Oil Companies International Marine
Forum, 1994 can be used. 123
At each time step of the simulation, the anchor-
to-fairlead horizontal distance and azimuth are
calculated for each line. The azimuth of a
mooring line is defined by the angle that the
fairlead-to-anchor segment makes with X. See
Figure 123

Dk X O X k xk cos y k sin 2 YO Yk xk sin y k cos 2

YO Yk xk sin y k cos
tan A k
X O X k xk cos y k sin

where: FIGURE 123


Dk is the anchor-to-fairlead horizontal ANCHOR-TO-FAIRLEAD HORIZONTAL
distance of mooring line number k; DISTANCE
Ak is the fairlead-to-anchor azimuth of AND MOORING LINE AZIMUTH
mooring line number k;
(XO,YO) are the co-ordinates of O
in the fixed axis system;
(Xk,Yk)are the co-ordinates in the fixed
axis system of the anchor of mooring
line number k;
(xk,yk) are the co-ordinates in the vessel
axis system of the fairlead of mooring
line number k;
is the instantaneous vessel heading.
Section 2 ANALYSIS PROCEDURE

124
Many research works are being undertaken
The instantaneous horizontal component of the world-wide to identify and better understand the
tension exerted by each line on the vessel is main damping causes.
obtained by interpolations on the distances using
the discretized characteristics which have been For the time being, any damping model
pre-calculated as detailed in 29 requires to be calibrated and the boundaries of its
field of applicability to be clearly identified. This
Horizontal components are then projected on has been done for the model adopted herein and
the vessel axis system and added together as per reference is made to 7 where the limitations of
the following equations: the present approach are specified.

NL In this methodology, damping effects are partly

FMx THk cos A k


generated by current loads since they are
calculated on the basis of the relative fluid
k 1 velocity (see 114) and partly by additional terms
which are proportional to the absolute speed of
NL the origin O of the vessel axis system, according

FMy THk sin A k


to the following formulae:

k 1 FB x B xxu
FB y B yy v
M M /O M B / O B


NL where:
x k THk sin A k y k THk cos A k u is the absolute velocity in surge of

k1 the origin O of the vessel axis system;
v is the absolute velocity in sway of the
where: origin O of the vessel axis system;
THk is the horizontal component of the is the derivative function versus
tension at fairlead number k; time of the vessel heading;
NL is the total number of working lines. Bxx is the linear damping coefficient in
surge;
125 Damping loads Byy is the linear damping coefficient in
sway;
126 B is the linear damping coefficient
Different theories exist to explain and model in yaw.
the damping effects on a moored vessel.
127
However, all of them are based on either fully
empirical or semi-analytical formulations the The linear damping coefficients may govern
range of validity of which is necessarily limited. the results of the simulation. Great care must
The principal problem is that damping sources therefore be paid to their evaluation.
are multiple and of various natures:
The values proposed in Table 127 account for
Hull friction and wave making resistance (in
all sources of damping listed in 126 provided that
calm water or in current); loads induced by the current are computed
Hydrodynamic mooring line damping (in according to the method specified in 114 These
calm water or in current); values are consistent with the safety factors given
Bottom friction effects on mooring lines; in 183

Second order wave drift damping vessel;


Damping due to mooring line dynamics in
waves;
Wind damping.
Section 2 ANALYSIS PROCEDURE

TABLE
VALUES OF THE LINEAR DAMPING COEFFICIENTS FOR DIFFERENT TYPES OF MOORING
SYSTEMS

MOORING LINEAR DAMPING COEFFICIENT (1)

SYSTEM
B xx
(2)
B yy
(2)
B
(2)

Barge or


tanker in
spread 0.06 K Oxx m Ma xx
0.06 K Oyy m Ma yy
0.10 K O I Ma m Ma yy x G 2

mooring

Barge or
tanker on an g g 0.083 L2 B yy
0.01 m 0.02 m
SPM(3) B
L

Semi-
submersible 0.20 K Oxx m Ma xx
0.20 K Oyy m Ma yy
0.10 K O I Ma m Ma yy x G 2

unit at
operating
draught

(1) Valid for usual mooring systems.

(2)
K Oxx K
Oyy K
O
are the diagonal terms of the mooring stiffness matrix
K
O
as defined in 84 and evaluated at the average
position of the vessel during the storm;

Maxx Mayy Ma are the diagonal terms of the asymptotic added mass matrix of the vessel as defined in
108,

m is the mass of the vessel;


I is the moment of inertia in yaw calculated at the centre of gravity G of the vessel;

L is the length of the vessel;

B is the breadth of the vessel.

(3) These values assume that the origin O of the vessel axis system is in the midship section.

128 Other loads Recognised methods should be used for the


load calculations, with due account for possible
129 interferences between the thrusters themselves or
between the thrusters and the hull and for all
Thrusters are sometimes used to assist a other phenomena that may modify, alter or
passive mooring system. degrade the thruster performances.

The driving system of the thruster loads may 130


be very simple (e.g. constant load in a constant
direction relative to the vessel heading), more The presence of a large number of risers
complicated (e.g. constant load in an absolute connected to the moored vessel may generate
direction) or extremely sophisticated (e.g. net additional loads that should be evaluated and
thrusts derived from vessel position and taken into account if they are not negligible.
kinematics, line tensions, etc.).

In any case, the loads that the thrusters


actually produce on the moored vessel should be
computed at each time step as function of the
applicable parameters.
Section 2 ANALYSIS PROCEDURE

131 Wave frequency response The complete signal of the vessel motion thus
obtained allows the instantaneous anchor-to-
fairlead distance of each line and its fairlead
132 General altitude to be calculated. The instantaneous
tensions can next be derived by interpolations in
The low frequency response of the moored the pre-calculated tension-offset curves.
vessel is obtained by numerical resolution in the Reference is made to 29
time domain of the vectorial differential equation
discussed in 13 and 100 At the end of each time
step of this numerical integration, the wave 133 Response Amplitude Operators
frequency motions are added for output purpose
only. In other words, the wave frequency 134
response that is calculated at any step does not If an Airy wave is applied to a floating vessel,
impact at all on the low frequency response either moored or not, the fluctuating motion
calculated at further steps. observed in any direction, around or about its
centre of gravity, is also a harmonic function of
The wave frequency motions are obtained by time, with the same frequency as the wave but
linear summation of those due to each usually with a phase lag. The amplitude of each
component of the wave. It is assumed in this motion is proportional to that of the wave. It also
process that wave frequency motions are not varies, as well as the phase, with the water depth
significantly influenced by the variations of and the wave incidence relative to the vessel.
mooring stiffness with low frequency motions.
Wave frequency motions are therefore computed See Table 134 Note that "around" is used for
for the average mooring stiffness corresponding translation displacements while about is used for
to the mean vessel position during the storm. rotations.

TABLE
MOTIONS OF A VESSEL SUBMITTED TO AN AIRY WAVE

Wave at COG MOTIONS OF THE VESSEL AROUND AND ABOUT ITS COG

Incidence Surge Sway Heave


a R x , cos t x , a Ry , cos t y , a Rz , cos t z ,
Elevation Roll Pitch Yaw

a cos t
a R , cos t , a R , cos t , a R , cos t ,
The Response Amplitude Operators (RAO's) close as possible to its average position during
are the functions which give in the vessel axis the storm.
system the amplitude and phase of the vessel
motions around and about its centre of gravity 135
when subjected to the action of an Airy wave of
The RAO's are normally available for the same
unitary amplitude.
incidences as QTF's and interpolations should be
performed in the same way during the
In Table 134, the RAO's of surge for instance
simulations. Discretization criteria are also the
are the functions Rx and x which, for a given
same as for QTF's.
incidence relative to the vessel heading, depend
only on the wave circular frequency.
Reference is made to 112
The RAO's can be obtained by model tests or
The phase functions are generally not
by a first order diffraction-radiation analysis using
continuous and it is therefore preferable to
a recognised computer program. In any case,
perform the interpolations on the in-phase and
they should be determined with due account for * **
the actual site water depth and the mooring out-of-phase amplitudes R and R which are
stiffness evaluated for a position of the vessel as continuous, rather than on the actual amplitude
and phase R and . Keeping the formulations of
Section 2 ANALYSIS PROCEDURE

Table 134, the relations between these 136 Instantaneous centre-of-gravity


parameters are: motion and vessel rotation

R* R cos 137


At each time step, the wave signal should be
determined at the instantaneous low frequency
R** R sin position of the vessel centre of gravity by means
of the equation of 91

138
Each of the six wave frequency motions of the
vessel around or about the low-frequency position
of its centre of gravity are then given by an
equation of the following form:

100
G t ai R i , H cos i t i k i X G cos H k i YG sin H
*

i 1
100

a R , sin t k X cos H k i YG sin H


**
i i H i i i G
i 1

where: FIGURE 138


H is the absolute direction where the COMBINATION OF LOW AND WAVE
wave comes from; FREQUENCY HORIZONTAL MOTIONS
(ai , i , i , ki)1i100 characterise the
family of Airy waves representative of
the wave spectrum;
(XG , YG) are the low frequency co-
ordinates of the vessel centre of gravity
in the fixed axis system XY;
is the low frequency vessel heading;
G designates any of the six wave
frequency motions of the vessel around
or about the low frequency position of its
centre of gravity (these motions are
those projected on the low-frequency
vessel axis system Oxy);
R
*
i , H is the in-phase 139 Instantaneous fairlead motion
amplitude of motion G interpolated
140
in the RAO's for the wave frequency
i and the relative wave incidence H At each time step, the wave frequency motions
; make the vessel turn about its centre of gravity.
The three components of the rotation vector in
R
**
i , H same for out-of- the low frequency vessel axis system are the roll,
pitch and yaw angles calculated according to the
phase motion.
equation of 138 At the same time, the vessel
Figure 138 is an illustration of the moves by a translation whose vector components
superimposition of low and horizontal wave in the same axis system are the surge, sway and
frequency motions. heave displacements calculated in the same way.

Further to this rotation and translation, the


instantaneous co-ordinates of any fairlead of the
Section 2 ANALYSIS PROCEDURE

vessel with respect to the low frequency vessel


axis system are given by the following equation:

xF t x G x G t 11 t 21 t 31 t xF x G

y F t y G y G t 12 t 22 t 32 t yF y G
z t z z t t t t z z
F G G 13 23 33 F G

1 sin 2
where: sin
2

23 t x y z
2
( xF , y F , zF ) are the co-ordinates of 2

1 sin 2
2
the fairlead in the low frequency vessel sin
31 t y x z
axis system further to the wave 2 2
frequency rotation about G;
1 sin 2
2
sin
( xF , y F , zF ) are the permanent co- 32 t x y z
2 2
ordinates of the fairlead in the vessel
sin 2
2
axis system; 1
33 t cos

z2
(x G,yG,zG) are the permanent co- 2 2
ordinates of the vessel centre of gravity
in the vessel axis system. Note that in 141
principle yG is equal to 0;
The instantaneous fairlead position in the fixed
( x G , y G , z G ) are the three axis system is then given by the following
components, in the low frequency equation:
vessel axis system, of the wave
frequency motion of the vessel centre of X F X G cos sin 0 xF t x G
gravity (surge, sway and heave);
[ij(t)] is the matrix of the instantaneous YF YG sin cos 0 y F t y G
Z Z 0 0 1 zF t z G
rotation whose coefficients are specified F G
hereafter. The following notations are
used to simplify the expressions of the
where:
matrix coefficients:
x G t ( wave frequency roll) ( X F , YF Z F ) are the
instantaneous co-ordinates of the
y G t ( wave frequency pitch)
fairlead in the fixed axis system XY;
z G t ( wave frequency yaw )
( X G , YG Z G ) are the
2 2 2
x y z instantaneous co-ordinates of the vessel
centre of gravity in the fixed axis system
sin 2
2
1 XY;
11 t cos
2
x
2 2 is the low frequency heading of the
1 sin 2
2
sin vessel.
12 t z xy
2 2
142 Instantaneous line tension
1 sin 2
2
sin
13 t y xz
2 2 143
1 sin 2
2
sin The instantaneous fairlead-to-anchor horizontal
21 t z xy
2 2 distance is given by:

sin 2
2
1
22 t cos
2 2
y2 DF X F X A 2 YF YA 2

where:
Section 2 ANALYSIS PROCEDURE

DF is the anchor-to-fairlead horizontal these forces depend on the vessel heading, an


distance; iterative process should be implemented
(XA , YA) are the co-ordinates of
the anchor in the fixed axis system; If other loads such as riser or thruster loads are
existing, their mean values should be determined
(XF , YF) are the instantaneous co- for the calculation of the equilibrium position.
ordinates of the fairlead in the fixed axis Here again, an iterative process even more
system, as calculated in 141 complicated, may be necessary since such loads
may vary with the three equilibrium position
144
parameters.
The offset tension curves of the mooring line
have been pre-calculated for three altitudes of the 148
fairlead as explained in 29 Those two curves
If several equilibrium positions exist, the most
calculated for fairlead altitudes on each side of
stable should be selected.
the instantaneous one should be identified. The
axial tension at the fairlead is then to be
Figure 148 is an illustration of the four possible
interpolated between these two curves.
equilibrium positions of a tanker moored by a
single line under the action of the current.
See Figure 144
Obviously, the second position is totally unstable.
The first position is the most expected. However,
the last two positions, when they exist, are
FIGURE 144
generally more stable.
INTERPOLATION OF THE
INSTANTANEOUS TENSION BETWEEN
THE OFFSET-TENSION CURVES PRE- FIGURE 148
CALCULATED FOR 3 FAIRLEAD THEORITICAL EQUILIBRIUM POSITIONS
ALTITUDES OF A TANKER
MOORED BY A SINGLE LINE IN THE
CURRENT

145 Simulation procedure

146 Starting position

147
Numerical transitory responses are prone to
generate temporary instabilities that may bias the
results. In order to minimise such effects, the
simulation should be started from the equilibrium
position of the moored vessel under the effect of
mean wind, current and wave drift forces. As
Section 2 ANALYSIS PROCEDURE

149 The beginning of the response signal is


therefore to be truncated prior to any statistics
In order to avoid an undue numerical transitory
being carried out. The length of time to discard is
response at the beginning of the simulation, it is
something between half and twice the natural
recommended that the vessel not be suddenly
period of the mooring system.
submitted to the time varying components of
environmental loads. A possible solution is to use
See Figure 149
a ramp function to apply them progressively or to
start at an instant where all loads are equal to
their mean value.

FIGURE
TRUNCATION OF THE BEGINNING OF THE RESPONSE SIGNAL TO ELIMINATE NUMERICAL
TRANSIENT

Truncated period Low frequency Complete


Response Response
2900

2800

2700

2600
tension (kN)

2500

2400

2300

2200

2100

2000
0

100

400

700
200

300

500

600

time (s)

A more accurate approach may consist of the 152


computation of the first four moments of the
For each mooring line, the maximum tension at
response signal using increasing truncated
fairlead reached during each simulation should be
intervals until these moments remain stable.
identified and recorded. The mean and the
standard deviation of the five maxima thus
150 Intact condition obtained are then to be computed. The design
tension of the line in intact condition is
151 defined as this mean value plus two standard
deviations as per the following formulae:
For any sea state to be considered (see
Section 169), five simulations of at least three 5
1
hours each should be performed using different
sets of elementary waves representative of the
T T
5 k 1
k
whole spectrum. If a wind gust spectrum is used,
the same provision applies. 5

T
1 2
T k T
The response signals should be built up with a 5 k 1
time step equal to or less than one tenth of the
peak or zero-up crossing period of the wave TIC T 2 T
spectrum, whichever is the more appropriate.
where:
Section 2 ANALYSIS PROCEDURE

Tk is the maximum tension at fairlead Homogeneous segments which are partly


obtained during the k-th simulation; laying on the seabed and partly lifted are
therefore to be split up into sub-segments for
T is the mean of the maximum
searching the maximum tension. The same
tensions obtained by the five should be done for laying segments crossing
simulations; peaks or troughs of the seabed.
T is the standard deviation of the
maximum tensions obtained by the five 154 One-line failure
simulations;
TIC is the design tension of the line in 155
intact condition.
The procedure described in 156 to 158 should
153 be repeated for each line. Provisions of 153 are
also to be applied.
For a given tension at the fairlead, the
maximum tension over the whole line length can 156
be at another point. This is the case of composite
lines including sub-sea buoys or floating For each sea state to be considered and for
segments or when the seabed shape is not flat. each mooring line, the simulation to be selected
Moreover, segments submitted to lower tensions for the study of the one-line failure case should
than that at the fairlead may be more critical satisfy the following criteria:
because of a weaker strength than that of the the maximum fairlead tension is the closest
segment connected to the fairlead. to its design tension in intact condition as
defined in 152, and
The tension at any point of the line when its
fairlead is submitted to the design tension is this maximum occurs while the low
subsequently to be analysed. This can be done frequency component of the tension is
using the equations detailed in 29 However, the increasing.
following features of catenary and taut thread
If the second criterion is not met, the selection
equations allow the maximum tension in each
procedure should be resumed with the highest
segment to be easily found:
second maximum in the same simulation or the
the tension in a suspended homogeneous closest second simulation.
segment is maximum at one of its two
ends, whether the segment crosses the free The two instants around the maximum tension
surface or not, where the low frequency component of the
tension is respectively minimum and maximum
the tension in a taut thread homogeneous
should be identified.
segment is maximum at one of its two
ends.
See Figure 156
Section 2 ANALYSIS PROCEDURE

FIGURE
RANGE OF TIME TO BE SELECTED FOR THE ONE-LINE FAILURE CASE

Maximum tension in
2750
intact condition
2700

2650

2600

2550
tension (kN)

2500

2450

2400

2350

2300

2250
10000

10100

10400

10500

10600

10700
10200

10300

Half low frequency


period of interest time (s)

157 161
Using always the same sets of Airy waves and Two kinds of failure should be investigated:
wind components as those used in the simulation
1) the total loss of one thruster, the other
identified in 156, five simulations are repeated;
thrusters having two thirds of their
during these simulations, the line should be
maximum thrust capacity available,
broken at different times equally distributed
between the two instants identified in 156 2) loss of half the total thrust capacity.

All five simulations should be run from the These two cases lead to the same remaining
beginning in order to ensure the same initial thrust capacity for a unit equipped with four
numerical transient. They can be terminated after identical thrusters.
two low frequency cycles following the line failure.
162
The maximum tension obtained at the fairlead
The design tensions in the mooring lines
of each remaining line after the line failure should
should be determined by a method similar to that
be identified for all five simulations. For this
described in Section 154 for the one-line failure
particular one-line-failure case, the tension in
case. The instants of failure, however, can be
damage condition of each remaining line is the
randomly selected during the simulation identified
average of its five maxima thus obtained.
in 156 in a sufficient number to ensure that the
statistics derived from the response samples are
158
reasonably representative.
For a given line, the design tension in
damage condition is the maximum of all
possible one-line failure cases. 163 Fatigue damage

164
159 Thruster failure
A fatigue analysis should be performed for
installations intended to stay moored on site for a
160
period longer than two years. Change of vessel
The thruster failure is assumed not to be draught, if any during the period of exposure,
concomitant with a mooring line failure. should be accounted for.
Section 2 ANALYSIS PROCEDURE

165
167
For each environmental condition selected as
The method used for counting and classing the
specified in 173, the fatigue damage in any
tension ranges should account for both low and
component of the mooring line is obtained by
wave frequency cycles.
means of the Miners ratio calculated for one year
(31 557 600 seconds):
168
3.15576 10
7
n jk The T-N curve used should be appropriate to
Dj p j
dj
N k
the mooring line component and to the type of
k
excitation encountered. Without further
information, the T-N curves specified in Table 168
where: can be used for mooring line components
Dj is the fatigue damage accumulated submitted to pure tension cycles.
over one year by the component under
the environmental condition number j;
TABLE 168
pj is the probability of occurrence of the
environmental condition number j (the T-N CURVES OF MOORING LINE
sum of the probabilities of all selected COMPONENTS
environmental conditions should be
equal to 1); T-N CURVE EQUATION

dj is the duration of the simulation of


BL
the environmental condition number j N
where:
(normally 10 800 seconds as specified T
in 150); BL is the component breaking load
T is the tension range (double amplitude)
njk is the number of cycles within the
tension range interval number k
Component
encountered by the component under
the environmental condition number j;
Common chain link 3,36 370
Nk is the number of cycles to failure at
tension range k as given by the Kenter link 3,36 90
appropriate T-N curve.
Six-strand wire rope 4,09 231
The total damage D accumulated over one
year by the component is then given by the Spiral-strand wire rope 5,05 166
following formula:

3.15576 10
7
n jk
D Dj p
i j k
j
dj Nk

166
The maximum width of each tension range
interval to be considered should be less than or
equal to one thousandth of the mooring line
component breaking load, or 10 kN, whichever is
less.
Section 3 CRITERIA

SECTION 169 CRITERIA

170 Environmental criteria 177


For temporary installations in the vicinity of
171 General other offshore structures, N should not be taken
less than 10. Special considerations as regards
172 the risk due to collision may however require that
a bigger figure be used.
The mooring system in design condition should
be verified for the maximum environmental 178
conditions which are defined as the combination
of the N-year return period wave with its If a unit working alongside another platform
associated wind and current. The value of N can move away and reach a pre-defined stand-off
depends on the type of installation as specified in position within a time reasonably compatible with
174 However, other combinations having the the local meteo forecast, the provisions of 176
same or a greater probability of occurrence can be applied to the stand-off position. The
should be used instead, if they lead to a more actual environmental conditions beyond which the
severe response of the mooring system. For unit is to move away should be specified in the
instance, overrunning of a vessel moored to an operation manual and the mooring system in
SPM system generally occurs in relatively mild working position should be checked for such
sea state conditions. conditions.

173 179

When and where a fatigue analysis is required Above return periods apply to manned
as specified in 163, a series of environmental installations. For unmanned installations, criteria
conditions representative of the site should be should be defined on a case by case basis.
defined. To do so, directional wave scatter
diagrams applicable to the site can be used to 180 Directionality
determine the wave statistical parameters and
their associated probabilities of occurrence. The 181
directions of wind and current associated to each
wave condition should be defined in the same The mooring system in either design or fatigue
way as explained in 180. The complete scatter condition should be checked for at least eight
diagram of each direction can be reduced by wave incidences covering 360. Directions of
appropriate means in order to limit the total wind and current should be those associated with
number of environmental conditions. As a general each wave incidence. Without relevant
rule, 5 to 10 environmental conditions per information about the relative direction of wind
direction is satisfactory. and current, the following combinations are at
least to be verified for each wave incidence:

174 Return period for design wind and current acting in the same
condition direction as wave;
wind acting in the same direction as wave,
175 current crossing wave with an angle of
22.5;
For permanent installations, N should not be
taken less than 100. current acting in the same direction as
wave, wind crossing wave with an angle of
176 22.5;
For temporary installations far away from any parallel wind and current crossing wave
other offshore structure, N can be reduced to five with an angle of 22.5;
times the duration in years of the operation, wind and current crossing wave, each on
without nevertheless being less than 5.
one side with an angle of 22.5.
Section 3 CRITERIA

See Figure 181 mooring system response to slight variations of


the wave incidence about the main directions
selected.
FIGURE 181
COMBINATIONS OF WAVE, WIND AND
CURRENT DIRECTIONS 183 Design criteria

Wave 184 Mooring lines


Wind
185
Current
Two types of lines are defined:
Parallel wind wave current Wave incidences 1) Lines of TYPE I are those lines the failure
of which leads to a transient response that
moves the moored vessel towards an
installation in close proximity, if any.
-22.5 -22.5
2) Other lines are of TYPE II
+22.5 +22.5
186
The safety factor is defined as the ratio
between the catalogue breaking load of the
Crossing current Crossing wind mooring line component and the maximum
tension occurring over its length when the design
tension as determined in 145 is applied to the
+22.5 fairlead.
-22.5 -22.5
187
+22.5 +22.5
At any point of the line, the safety factor should
-22.5 be greater than or equal to those specified in
Table 187

Parallel wind and current Symmetric wind and current


crossing wave crossing wave
TABLE 187
MINIMUM SAFETY FACTORS APPLICABLE
182 TO MOORING LINES
If symmetry considerations are applicable or
prevailing directions of one or more parameters TYPE OF LINE MOORING SYSTEM CONDITION
are known, the number of incidences or
combinations may consistently be reduced. Intact Damage (1)

However, symmetry considerations should not


lead to irrelevant simplifications. TYPE I 2,20 1,75

For instance, if maximum environment TYPE II 1,75 1,25


conditions are nearly the same whatever the
wave incidence, the response of a fully symmetric (1) Damage condition applies to both one-line failure and
turret mooring system in such conditions should thruster failure cases
be verified for at least wave incidences parallel to
one mooring line and wave incidences parallel to
the secant of two mooring lines. 188 Anchors

Besides, it should be noted that a precise 189


knowledge of the environment directionality is For drag anchors, the safety factor is the ratio
very rare. In addition, actual incidences of wave, between the maximum holding power applicable
current and wind may vary during a storm. It is to the mooring site, and the tangent-to-the-
therefore needed to check the sensitivity of the seabed component of the tension in the line at the
Section 3 CRITERIA

anchoring point when the design tension is


193
applied to the fairlead.
In any condition, the distance between any
The perpendicular-to-the-seabed component of point of the mooring line and any external object,
the load applied to a drag anchor when the line is either suspended or laying on the seabed, the
submitted to its design tension at fairlead, should contact with which could be damageable for
remain less than 20% of its wet weight projected either the mooring line or the aforesaid object,
onto the same direction. should remain greater than or equal to five
metres.
For pile driven anchors, the safety factor is the
ratio between the maximum load that the pile can
withstand in compliance with the requirements of
applicable structural and soil mechanics codes 194 Fatigue life
and standards, and the tension in the line at the
anchoring point when the design tension is 195 General
applied to the fairlead.
196
190
The fatigue life, in years, of any component of
The minimum safety factors applicable to
the mooring line is defined as the invert of its
anchors are specified in Table 190
fatigue damage as calculated in 163

197
TABLE 190
The service life is the duration, in years, that
MINIMUM SAFETY FACTORS APPLICABLE the moored vessel is intended to stay on site.
TO ANCHORS
198
TYPE OF LINE MOORING SYSTEM CONDITION
The fatigue life of the mooring lines should be
determined and compared to the service life of
Intact Damage (1)

the installation if this latter is greater than 2 years.


TYPE I 2,00 1,60
199 Minimum fatigue life
TYPE II 1,60 1,15
The fatigue life of any component of the
(1) Damage condition applies to both one-line failure and mooring line (i.e. chains, wire or synthetic ropes,
thruster failure cases shackles, Kenter links, etc.) should be at least 5
times greater than the service life of the
191 Clearance installation. This factor accounts for both the line
dynamics which is not included in the calculation
methodology and for the uncertainties in Miner
192
theory, load prediction, T-N curves and other data
In any condition, the distance between any used in the fatigue assessment process.
external object which the moored vessel is
normally not in contact with, and any part of the
vessel should remain greater than or equal to one
metre.
Section 3 CRITERIA

APPENDIX A SENSITIVITY TO INPUT


PARAMETERS

A.1 Account for mooring stiffness in RAO's calculations

Effect of mooring stiffness on RAO


2,0

1,8
0
1,6 0,1
Increasing mooring stiffness 0,2
1,4
0,3
Motion (m/m)

1,2
0,4
1,0 0,5
0,6
0,8
0,7
0,6
0,8
0,4 0,9
1
0,2

0,0
0,2 0,4 0,6 0,8 1,0 1,2 1,4
Wave circular frequency (rad/s)

A.2 Linear damping matrix A.3 Instant of line failure

428 t / s 0 0

0 1009 t / s 0 Effect of the instant of failure on
0 0 1,123.109 tm2 / s the peak tension in adjacent line
Maximum tension
in adjacent line

4000
3500
(kN)

3000
2500
2000
1500
5950 6000 6050 6100 6150

Instant of failure (s)


Appendix B TYPICAL VESSEL DATA

APPENDIX B TYPICAL VESSEL DATA


Appendix B TYPICAL VESSEL DATA
B.1.1 Generic 140,000 tdw tanker ballast draught 7.40 m waterdepth 25 m

B.1 Generic tankers B.1.1 Generic 140,000 tdw tanker


Longitudinal drift force Transverse drift force Drift yaw moment
0 0 4000
-100 2000
-20
0

QTF (kN.m/m)
-200
QTFx (kN/m)

QTFy (kN/m)
-40 -2000
-300 -4000
-60
-400 -6000
-80 -8000
-500
-10000
-100 -600 -12000
-120 -700 -14000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45

Longitudinal drift force Transverse drift force Drift yaw moment


200 0 40000
150 -200 30000
100 20000

QTF (kN.m/m)
-400
QTFx (kN/m)

QTFy (kN/m)
50 10000
-600
0 0
-800
-50 -10000
-100 -1000 -20000
-150 -1200 -30000
-200 -1400 -40000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5

Longitudinal drift force Transverse drift force Drift yaw moment


120 0 20000

100 -100
15000

QTF (kN.m/m)
-200
QTFx (kN/m)

QTFy (kN/m)

80
-300 10000
60
-400 5000
40
-500
20 0
-600
0 -700 -5000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180

Surge RAO Sway RAO Yaw RAO


2 2 1.2
1.8 1.8
1.6 1.6 1
1.4 1.4
)
Appendix B TYPICAL VESSEL DATA
B.1.1 Generic 140,000 tdw tanker ballast draught 7.40 m waterdepth 25 m

B.1 Generic tankers B.1.1 Generic 140,000 tdw tanker


Longitudinal drift force Transverse drift force Drift yaw moment
0 0 4000
-100 2000
-20
0

QTF (kN.m/m)
-200
QTFx (kN/m)

QTFy (kN/m)
-40 -2000
-300 -4000
-60
-400 -6000
-80 -8000
-500
-10000
-100 -600 -12000
-120 -700 -14000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45

Longitudinal drift force Transverse drift force Drift yaw moment


200 0 40000
150 -200 30000
100 20000

QTF (kN.m/m)
-400
QTFx (kN/m)

QTFy (kN/m)
50 10000
-600
0 0
-800
-50 -10000
-100 -1000 -20000
-150 -1200 -30000
-200 -1400 -40000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5

Longitudinal drift force Transverse drift force Drift yaw moment


120 0 20000

100 -100
15000

QTF (kN.m/m)
-200
QTFx (kN/m)

QTFy (kN/m)

80
-300 10000
60
-400 5000
40
-500
20 0
-600
0 -700 -5000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180

Surge RAO Sway RAO Yaw RAO


2 2 1.2
1.8 1.8
1.6 1.6 1
1.4 1.4
)
Appendix B TYPICAL VESSEL DATA
B.1.1 Generic 140,000 tdw tanker ballast draught 7.40 m waterdepth 25 m

B.1 Generic tankers B.1.1 Generic 140,000 tdw tanker


Longitudinal drift force Transverse drift force Drift yaw moment
0 0 4000
-100 2000
-20
0

QTF (kN.m/m)
-200
QTFx (kN/m)

QTFy (kN/m)
-40 -2000
-300 -4000
-60
-400 -6000
-80 -8000
-500
-10000
-100 -600 -12000
-120 -700 -14000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45

Longitudinal drift force Transverse drift force Drift yaw moment


200 0 40000
150 -200 30000
100 20000

QTF (kN.m/m)
-400
QTFx (kN/m)

QTFy (kN/m)
50 10000
-600
0 0
-800
-50 -10000
-100 -1000 -20000
-150 -1200 -30000
-200 -1400 -40000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5

Longitudinal drift force Transverse drift force Drift yaw moment


120 0 20000

100 -100
15000

QTF (kN.m/m)
-200
QTFx (kN/m)

QTFy (kN/m)

80
-300 10000
60
-400 5000
40
-500
20 0
-600
0 -700 -5000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180

Surge RAO Sway RAO Yaw RAO


2 2 1.2
1.8 1.8
1.6 1.6 1
1.4 1.4
)
Appendix B TYPICAL VESSEL DATA
B.1.1 Generic 140,000 tdw tanker ballast draught 7.40 m waterdepth 25 m

B.1 Generic tankers B.1.1 Generic 140,000 tdw tanker


Longitudinal drift force Transverse drift force Drift yaw moment
0 0 4000
-100 2000
-20
0

QTF (kN.m/m)
-200
QTFx (kN/m)

QTFy (kN/m)
-40 -2000
-300 -4000
-60
-400 -6000
-80 -8000
-500
-10000
-100 -600 -12000
-120 -700 -14000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45

Longitudinal drift force Transverse drift force Drift yaw moment


200 0 40000
150 -200 30000
100 20000

QTF (kN.m/m)
-400
QTFx (kN/m)

QTFy (kN/m)
50 10000
-600
0 0
-800
-50 -10000
-100 -1000 -20000
-150 -1200 -30000
-200 -1400 -40000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5

Longitudinal drift force Transverse drift force Drift yaw moment


120 0 20000

100 -100
15000

QTF (kN.m/m)
-200
QTFx (kN/m)

QTFy (kN/m)

80
-300 10000
60
-400 5000
40
-500
20 0
-600
0 -700 -5000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180

Surge RAO Sway RAO Yaw RAO


2 2 1.2
1.8 1.8
1.6 1.6 1
1.4 1.4
)
Appendix B TYPICAL VESSEL DATA
B.1.1 Generic 140,000 tdw tanker ballast draught 7.40 m waterdepth

Longitudinal drift force Transverse drift force Drift yaw moment


0 0 4000
-10 2000
-100
-20 0

QTF (kN.m/m)
-30 -200
QTFx (kN/m)

QTFy (kN/m)
-2000
-40 -300 -4000
-50
-400 -6000
-60
-70 -8000
-500
-80 -10000
-600 -12000
-90
-100 -700 -14000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45

Longitudinal drift force Transverse drift force Drift yaw moment


200 0 25000
20000
150 -200
15000

QTF (kN.m/m)
100 -400 10000
QTFx (kN/m)

QTFy (kN/m)
50 -600 5000
0
0 -800 -5000
-50 -1000 -10000
-15000
-100 -1200
-20000
-150 -1400 -25000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5

Longitudinal drift force Transverse drift force Drift yaw moment


100 0 20000
90 18000
-100
80 16000
14000

QTF (kN.m/m)
70 -200
QTFx (kN/m)

QTFy (kN/m)

60 12000
-300 10000
50
-400 8000
40 6000
30 -500 4000
20 2000
-600
10 0
0 -700 -2000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180

Surge RAO Sway RAO Yaw RAO


1 1 0.4
0.9 0.9 0.35
0.8 0.8
0.3
0.7 0.7
/m)

0.6 0.25
m)

m)

0.6
Appendix B TYPICAL VESSEL DATA
B.1.1 Generic 140,000 tdw tanker ballast draught 7.40 m waterdepth

Longitudinal drift force Transverse drift force Drift yaw moment


0 0 4000
-10 2000
-100
-20 0

QTF (kN.m/m)
-30 -200
QTFx (kN/m)

QTFy (kN/m)
-2000
-40 -300 -4000
-50
-400 -6000
-60
-70 -8000
-500
-80 -10000
-600 -12000
-90
-100 -700 -14000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45

Longitudinal drift force Transverse drift force Drift yaw moment


200 0 25000
20000
150 -200
15000

QTF (kN.m/m)
100 -400 10000
QTFx (kN/m)

QTFy (kN/m)
50 -600 5000
0
0 -800 -5000
-50 -1000 -10000
-15000
-100 -1200
-20000
-150 -1400 -25000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5

Longitudinal drift force Transverse drift force Drift yaw moment


100 0 20000
90 18000
-100
80 16000
14000

QTF (kN.m/m)
70 -200
QTFx (kN/m)

QTFy (kN/m)

60 12000
-300 10000
50
-400 8000
40 6000
30 -500 4000
20 2000
-600
10 0
0 -700 -2000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180

Surge RAO Sway RAO Yaw RAO


1 1 0.4
0.9 0.9 0.35
0.8 0.8
0.3
0.7 0.7
/m)

0.6 0.25
m)

m)

0.6
Appendix B TYPICAL VESSEL DATA
B.1.1 Generic 140,000 tdw tanker ballast draught 7.40 m waterdepth

Longitudinal drift force Transverse drift force Drift yaw moment


0 0 4000
-10 2000
-100
-20 0

QTF (kN.m/m)
-30 -200
QTFx (kN/m)

QTFy (kN/m)
-2000
-40 -300 -4000
-50
-400 -6000
-60
-70 -8000
-500
-80 -10000
-600 -12000
-90
-100 -700 -14000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45

Longitudinal drift force Transverse drift force Drift yaw moment


200 0 25000
20000
150 -200
15000

QTF (kN.m/m)
100 -400 10000
QTFx (kN/m)

QTFy (kN/m)
50 -600 5000
0
0 -800 -5000
-50 -1000 -10000
-15000
-100 -1200
-20000
-150 -1400 -25000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5

Longitudinal drift force Transverse drift force Drift yaw moment


100 0 20000
90 18000
-100
80 16000
14000

QTF (kN.m/m)
70 -200
QTFx (kN/m)

QTFy (kN/m)

60 12000
-300 10000
50
-400 8000
40 6000
30 -500 4000
20 2000
-600
10 0
0 -700 -2000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180

Surge RAO Sway RAO Yaw RAO


1 1 0.4
0.9 0.9 0.35
0.8 0.8
0.3
0.7 0.7
/m)

0.6 0.25
m)

m)

0.6
Appendix B TYPICAL VESSEL DATA
B.1.1 Generic 140,000 tdw tanker fully loaded draught 16.37 m waterdepth 25 m

Longitudinal drift force Transverse drift force Drift yaw moment


0 0 10000
-100 5000
-50 -200

QTF (kN.m/m)
0
QTFx (kN/m)

QTFy (kN/m)
-300
-100 -5000
-400
-500 -10000
-150
-600
-15000
-200 -700
-800 -20000
-250 -900 -25000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45

Longitudinal drift force Transverse drift force Drift yaw moment


400 0 50000

300 -200 40000


-400

QTF (kN.m/m)
200 30000
QTFx (kN/m)

QTFy (kN/m)
-600
100 -800 20000
0 -1000 10000
-1200
-100 0
-1400
-200 -1600 -10000
-300 -1800 -20000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5

Longitudinal drift force Transverse drift force Drift yaw moment


350 0 20000
300 -100
15000
-200

QTF (kN.m/m)
250
QTFx (kN/m)

QTFy (kN/m)

-300 10000
200 -400
5000
150 -500
-600 0
100
-700
50 -5000
-800
0 -900 -10000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180

Surge RAO Sway RAO Yaw RAO


2 2 1.2
1.8 1.8
1.6 1.6 1
1.4 1.4 0.8
g/m)
/m)

/m)

1.2 1.2
Appendix B TYPICAL VESSEL DATA
B.1.1 Generic 140,000 tdw tanker fully loaded draught 16.37 m waterdepth 25 m

Longitudinal drift force Transverse drift force Drift yaw moment


0 0 10000
-100 5000
-50 -200

QTF (kN.m/m)
0
QTFx (kN/m)

QTFy (kN/m)
-300
-100 -5000
-400
-500 -10000
-150
-600
-15000
-200 -700
-800 -20000
-250 -900 -25000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45

Longitudinal drift force Transverse drift force Drift yaw moment


400 0 50000

300 -200 40000


-400

QTF (kN.m/m)
200 30000
QTFx (kN/m)

QTFy (kN/m)
-600
100 -800 20000
0 -1000 10000
-1200
-100 0
-1400
-200 -1600 -10000
-300 -1800 -20000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5

Longitudinal drift force Transverse drift force Drift yaw moment


350 0 20000
300 -100
15000
-200

QTF (kN.m/m)
250
QTFx (kN/m)

QTFy (kN/m)

-300 10000
200 -400
5000
150 -500
-600 0
100
-700
50 -5000
-800
0 -900 -10000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180

Surge RAO Sway RAO Yaw RAO


2 2 1.2
1.8 1.8
1.6 1.6 1
1.4 1.4 0.8
g/m)
/m)

/m)

1.2 1.2
Appendix B TYPICAL VESSEL DATA
B.1.1 Generic 140,000 tdw tanker fully loaded draught 16.37 m waterdepth 25 m

Longitudinal drift force Transverse drift force Drift yaw moment


0 0 10000
-100 5000
-50 -200

QTF (kN.m/m)
0
QTFx (kN/m)

QTFy (kN/m)
-300
-100 -5000
-400
-500 -10000
-150
-600
-15000
-200 -700
-800 -20000
-250 -900 -25000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45 Inc. = 0 Inc. = 22.5 Inc. = 45

Longitudinal drift force Transverse drift force Drift yaw moment


400 0 50000

300 -200 40000


-400

QTF (kN.m/m)
200 30000
QTFx (kN/m)

QTFy (kN/m)
-600
100 -800 20000
0 -1000 10000
-1200
-100 0
-1400
-200 -1600 -10000
-300 -1800 -20000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5 Inc. = 67.5 Inc. = 90 Inc. = 112.5

Longitudinal drift force Transverse drift force Drift yaw moment


350 0 20000
300 -100
15000
-200

QTF (kN.m/m)
250
QTFx (kN/m)

QTFy (kN/m)

-300 10000
200 -400
5000
150 -500
-600 0
100
-700
50 -5000
-800
0 -900 -10000
0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 0 0.2 0.4 0.6 0.8 1 1.2 1.4
(rad/s) (rad/s) (rad/s)

Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180 Inc. = 135 Inc. = 157.5 Inc. = 180

Surge RAO Sway RAO Yaw RAO


2 2 1.2
1.8 1.8
1.6 1.6 1
1.4 1.4 0.8
g/m)
/m)

/m)

1.2 1.2
Appendix B TYPICAL VESSEL DATA
B.1.1 Generic 140,000 tdw tanker fully loaded draught 16.37 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.1 Generic 140,000 tdw tanker fully loaded draught 16.37 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.1 Generic 140,000 tdw tanker fully loaded draught 16.37 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.2 Generic 210,000 tdw tanker ballast draught 8.86 m waterdepth 25 m

B.1.2 Generic 210,000 tdw tanker


Appendix B TYPICAL VESSEL DATA
B.1.2 Generic 210,000 tdw tanker ballast draught 8.86 m waterdepth 25 m
Appendix B TYPICAL VESSEL DATA
B.1.2 Generic 210,000 tdw tanker ballast draught 8.86 m waterdepth 25 m
Appendix B TYPICAL VESSEL DATA
B.1.2 Generic 210,000 tdw tanker ballast draught 8.86 m waterdepth 25 m
Appendix B TYPICAL VESSEL DATA
B.1.2 Generic 210,000 tdw tanker ballast draught 8.86 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.2 Generic 210,000 tdw tanker ballast draught 8.86 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.2 Generic 210,000 tdw tanker ballast draught 8.86 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.2 Generic 210,000 tdw tanker loaded draught 16.30 m waterdepth 25 m
Appendix B TYPICAL VESSEL DATA
B.1.2 Generic 210,000 tdw tanker loaded draught 16.30 m waterdepth 25 m
Appendix B TYPICAL VESSEL DATA
B.1.2 Generic 210,000 tdw tanker loaded draught 16.30 m waterdepth 25 m
Appendix B TYPICAL VESSEL DATA
B.1.2 Generic 210,000 tdw tanker loaded draught 16.30 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.2 Generic 210,000 tdw tanker loaded draught 16.30 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.2 Generic 210,000 tdw tanker loaded draught 16.30 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.3 Generic 280,000 tdw tanker ballast draught 10.02 m waterdepth 65 m

B.1.3 Generic 280,000 tdw tanker


Appendix B TYPICAL VESSEL DATA
B.1.3 Generic 280,000 tdw tanker ballast draught 10.02 m waterdepth 65 m
Appendix B TYPICAL VESSEL DATA
B.1.3 Generic 280,000 tdw tanker ballast draught 10.02 m waterdepth 65 m
Appendix B TYPICAL VESSEL DATA
B.1.3 Generic 280,000 tdw tanker ballast draught 10.02 m waterdepth 65 m
Appendix B TYPICAL VESSEL DATA
B.1.3 Generic 280,000 tdw tanker ballast draught 10.02 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.3 Generic 280,000 tdw tanker ballast draught 10.02 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.3 Generic 280,000 tdw tanker ballast draught 10.02 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.3 Generic 280,000 tdw tanker loaded draught 20.80 m waterdepth 65 m
Appendix B TYPICAL VESSEL DATA
B.1.3 Generic 280,000 tdw tanker loaded draught 20.80 m waterdepth 65 m
Appendix B TYPICAL VESSEL DATA
B.1.3 Generic 280,000 tdw tanker loaded draught 20.80 m waterdepth 65 m
Appendix B TYPICAL VESSEL DATA
B.1.3 Generic 280,000 tdw tanker loaded draught 20.80 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.3 Generic 280,000 tdw tanker loaded draught 20.80 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.1.3 Generic 280,000 tdw tanker loaded draught 20.80 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.2.1 Semi-submersible AKER-H3 operating draught 21.34 m waterdepth 200 m

B.2 Semi-submersible units


B.2.1 Semi-submersible AKER-H3
Appendix B TYPICAL VESSEL DATA
B.2.1 Semi-submersible AKER-H3 operating draught 21.34 m waterdepth 200 m

B.2 Semi-submersible units


B.2.1 Semi-submersible AKER-H3
Appendix B TYPICAL VESSEL DATA
B.2.1 Semi-submersible AKER-H3 operating draught 21.34 m waterdepth 200 m

B.2 Semi-submersible units


B.2.1 Semi-submersible AKER-H3
Appendix B TYPICAL VESSEL DATA
B.2.1 Semi-submersible AKER-H3 operating draught 21.34 m waterdepth 200 m

B.2 Semi-submersible units


B.2.1 Semi-submersible AKER-H3
Appendix B TYPICAL VESSEL DATA
B.2.2 Semi-submersible PACE SETTER operating draught 21 m waterdepth 800 m

B.2.2 Semi-submersible PACE SETTER


Appendix B TYPICAL VESSEL DATA
B.2.2 Semi-submersible PACE SETTER operating draught 21 m waterdepth 800 m
Appendix B TYPICAL VESSEL DATA
B.2.2 Semi-submersible PACE SETTER operating draught 21 m waterdepth 800 m
Appendix B TYPICAL VESSEL DATA
B.2.2 Semi-submersible PACE SETTER operating draught 21 m waterdepth 800 m
Appendix B TYPICAL VESSEL DATA
B.2.3 Semi-submersible SEDCO 700 operating draught 25.91 m waterdepth

B.2.3 Semi-submersible SEDCO 700


Appendix B TYPICAL VESSEL DATA
B.2.3 Semi-submersible SEDCO 700 operating draught 25.91 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.2.3 Semi-submersible SEDCO 700 operating draught 25.91 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.2.3 Semi-submersible SEDCO 700 operating draught 25.91 m waterdepth
Appendix B TYPICAL VESSEL DATA
B.3 NKP barge (220 m x 46 m x 10.5 m) draught 10.50 m waterdepth 170 m

B.3 NKP barge


Appendix B TYPICAL VESSEL DATA
B.3 NKP barge (220 m x 46 m x 10.5 m) draught 10.50 m waterdepth 170 m
Appendix B TYPICAL VESSEL DATA
B.3 NKP barge (220 m x 46 m x 10.5 m) draught 10.50 m waterdepth 170 m
Appendix B TYPICAL VESSEL DATA
B.3 NKP barge (220 m x 46 m x 10.5 m) draught 10.50 m waterdepth 170 m
Appendix B TYPICAL VESSEL DATA
B.3 NKP barge (220 m x 46 m x 10.5 m) draught 10.50 m waterdepth 170 m
Appendix C CALCULATION EXAMPLES

APPENDIX C CALCULATION EXAMPLES

C.1 Fitting of the material C.3 Characteristic equivalent


elastic curve by a 5th order to a fender
polynomial

C.2 Characteristic equivalent


to a CALM

Characteristic equivalent to a CALM

4000
tension (kN)

3000
2000
1000
0
0 20 40 60 80
Distance (m )

You might also like