Developing A Formula To Represent The Second-Order Wave Effects On Moored Floating Structures
Developing A Formula To Represent The Second-Order Wave Effects On Moored Floating Structures
Developing A Formula To Represent The Second-Order Wave Effects On Moored Floating Structures
of
SI
The effects of wave direction, frequency and the waterline shape of floating structures on the
wave mean drift force formula have been considered separately by several authors but there isn't
a general formula to take into account all of these effects. In this regard, Faltinsens wave drift
force formula has been modified by adding finite draft coefficient. The results obtained from
this formula which is dependent on wave frequencies, has been compared with Helvacioglus
experiments favorably in high wave frequencies. Moreover, the influence of the current on the
wave mean drift force has been taken into account by considering the current coefficient derived
from the ship added resistance formula. In addition, the formula for the calculation of the wave
drift damping has been extended to cover high wave frequencies as well as low wave
frequencies. The results compared with asymptotic formula showed good agreement in the high
frequency band.
Keywords: Hydrodynamics, Second-order wave effects, Floating structures
Ar
ch
ive
1. Introduction
In determining the total mooring
loads and motions of moored floating
structures such as FPSOs (Floating
Production, Storage and Offloading
systems) and SPMs (Single Point
Mooring systems), in addition to the first
order wave effects, the second-order
effects are also important. The secondorder wave effects yield additional
forces and responses which are
proportional to the square of wave
amplitudes. The presence of the secondorder effects can simply be illustrated
from
Bernoullis
equation
by
considering its quadratic velocity term
[1]. Helvacioglu[2] has done some
experiments in order to obtain drift force
coefficients. The second-order wave
forces can be derived from pressure
integration (near field) approach [3].
More details about the derivation can be
found in Pinkster [4] , Standing et al. [5],
13/E
1/E
www.SID.ir
Iranian Association of
Naval Architecture & Marine Engineering
1
2
=
2
1
2
2
2
V1 + V2 + V3
2
(3)
Where
,are velocity terms related to the
directions of axes of coordinate system.
Extending the first velocity term of Eq.
(3) for two wave components with
different wave amplitudes A1 and A2 ,
and different circular frequencies 1 and
2 propagating in an idealized sea state
leads to [1]:
2
V1 =
2
(4)
ive
of
SI
= ii
(1)
i =1
Ar
ch
Where
i is a small parameter (perturbation)
i is the ith order velocity potential so
for example; 2 denotes the second order
velocity potential
Now by considering a fixed coordinate
system the pressure in a point on the hull
of a floating structure (e.g. FPSO) can be
determined by writing down the
Bernoullis equation as:
p = P0 gz
1
t 2
(2)
+
+
cos (2 1t + 2 1 )
2
2
2
2
A2
cos (2 2 t + 2 2 )
2 2
A1 A2 cos (( 1 2 )t + 1 2 )
2
A1 A2 cos (( 1 + 2 )t + 1 + 2 )
2
A1
A
+ 2
2
2
):
cos (2 1t + 2 1 ) +
2
A2
cos (2 2 t + 2 2 ) +
2 2
2
Where
P0 is the static pressure (atmospheric
pressure),
is the gradient.
2/E
14/E
A1 A 2 cos ([ 1 + 2
]t + 1 + 2 )
www.SID.ir
( A1 A2 cos[(1 2 )t + 1 2 ]) :
ive
Ar
ch
sin
H dy = 1 2 g a2 Bsin 2 H
(5)
Where
H : is the tangential slope of the water
plane curve with respect to longitudinal
centre line.
B : is the diameter of the cylinder
a : is the wave amplitude
Also, the bar denotes the mean value.
(6)
Lr
Where
y
: is the same as sin H in the
l
Havelock's formula
L r : is the exposed water line part of a
floating structure as shown in
Fig. 1.
Incident
of
y
Fx = 1 g a2 dl
2
l
SI
Lr
Shadow
region
Ls
Wave
Pinkster
and
Oortmerssen
[11],
presented the direct pressure integration
method to obtain the mean wave forces
and moments. By considering a twodimensional piercing body subjected to a
very small wave length it can be
assumed that the hydrodynamic behavior
of the sea wave around the body is the
same as the hydrodynamic pattern
around a vertical plane wall (Fig. 2).
z
Wave
15/E
13/E
3/E
www.SID.ir
Iranian Association of
Naval Architecture & Marine Engineering
2 g a
e kz cos t cos kx
(7)
p = gz 1 1 + 1
t
2 y
z
2 g zdz = 2 g 1 z 2
2
0
a
0
= g a
Fi =
g a2
2
sin ( + )n dl
2
(8)
L1
Where;
F1 : is the wave surge drift force
component
F2 : is the wave sway drift force
component
F6 : is the wave yaw drift moment
component
n1 = sin
n 2 = cos
n6 = x cos y sin
L1 : is the integration domain, the non
shadow part of the water plane curve.
: is the wave propagation angle with
respect to x axis
: is tangential angle of water plane
curve
of
1 =
to
SI
ive
2
1 a = [ 2 g a ] a = 2 g a
t z = 0
2 2
1
+ 1 dz =
2 y z
0 1 4 g 2 a 2 2 2 kz 2
k e sin ky + cos2 ky dz
2
1442443
2 2
=1
ch
Ar
2 kz
16/E
4/E
14/E
1
= g a 2ke dz = g a 2
2
2
Shadow Zone
Wave
www.SID.ir
sin ( + 0)n dl
2
L1
ch
g a2
ive
Where = 90 , n1 = 1 , dl = dy and L1 is
the integration domain which for this
case is the floating structures beam.
Therefore the above formula would
become:
F1 =
Ar
F1 = CT
g a2
2
(9)
R1 ( ) = CT
(11)
of
F1 =
g a2
SI
F1
1
g a2 B
2
(10)
C cu ,i = 1 +
2U i
g
(12)
where
U 1 = U cos
U 2 = U sin = U 6
U : is the current speed
: is the current angle in respect to the x
axis.
13/E
17/E
5/E
www.SID.ir
Iranian Association of
Naval Architecture & Marine Engineering
D=10m
D=15m
D=20m
0.9
0.8
CT coefficient
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0.1
0.2
0.3
0.4
0.5
0.6
Frequency rad/sec.
0.7
0.8
0.9
ch
ive
of
SI
Ar
C cu ,i = 1 + 2Fni
(13)
where
Fni =
Ui
gDi
6/E
14/E
18/E
3.2.DERIVING
STEADY
DRIFT
FORCES AND MOMENTS FOR AN
FPSO
g a2
2
sin ( + )n dl
2
(14)
L1
2 U cos
F1 = C T 1 +
g
g a2
(15)
2
sin ( + )sin dl
L1
www.SID.ir
sin (0 + )sin(0)dx =0
Sway:
2 U sin g
F 2 = C T 1 +
g
2
sin (
2
2
a
(16)
+ ) cos dl
II
Part 3:
dl III = rd
L1
Yaw:
)( x cos
(17)
F1 = CT 1 +
g
2
1 2
2
2 2
1 2
r cos sin + + 0 + r sin + cos
344443443443
344443444434 4244
1
4
r cos
3
of
2 2
1
r sin 2 + cos
3 3
3
y sin )dl
L1
2U cos
F1 = CT 1 +
g
sin ( + )sin rd =
sin (
2
2
a
SI
2 U sin g
F 6 = C T 1 +
g
2
g a2
2
sin ( + )sin dl
2 L1
2 2U cos
CT 1 +
3
g
g a2 r cos
(18)
ch
ive
2U cos g a2
F1 = CT 1 +
g
2
I
3
14442444
part1
42444
3
144
II
III
part
2
14442444
3
part 3
F1 =
Ar
Part 1:
dl I = rd
sin ( + )sin rd =
2
1
2
2
r cos sin 2 +
3
3
3
where r =
B
2
Part 2:
dl II = dx ; = 0
2U sin
F2 = CT 1 +
g
g a2
sin 2 ( + ) cos dl
L1
19/E
13/E
7/E
www.SID.ir
Iranian Association of
Naval Architecture & Marine Engineering
2U sin g a2
F2 = CT 1 +
g
2
sin 2 ( + ) cos dl I +
144424443
part 4
14442444
III
part 5
1
44424443
II
part 6
2U sin g a2
F2 = CT 1 +
g
2
4
(19)
g a2
Part 4:
dl I = rd
L1
sin ( + )cosrd =
2
2U sin g a2
F6 = CT 1 +
g
Part 5:
dl II = dx ; = 0
sin (0 + )cos(0)dx = L
2
sin sin
of
II
I
44444
42444444
3
1
part 7
1444442444443
part 8
II
SI
2
2
ive
sin ( + )cos rd =
ch
Ar
2
2
2
a
2
2
4
r sin + L2 sin sin
3
or
20/E
8/E
14/E
Part 7:
dl I = rd ; x = L1 r sin ;
y = r cos
0
2
= rL1 sin sin cos
3
3
Part 8:
dl II = dx ; = 0
II
2
2
L2 L1
=
sin sin
2
Part 9:
dl I = rd ; x = L2 + r sin ;
y = r cos
www.SID.ir
2
= rL2 sin + sin cos
3
3
2U sin g a2
F6 = CT 1 +
g
5
4
S ( ) 0.11 T1
T1
=
exp 0.44
2 2
H s2T1
2
where:
S ( )d = m
ive
3.3.DRIFT
WAVE
IRREGULAR WAVES
LOADS
IN
ch
(21)
Ar
F ( , , )
d
Fi s = 2 s( ) i 2
i = 1,...,6
where
S ( ) is the sea spectrum
Fi ( , , ) is the ith mean force
component in regular waves with the
circular frequency of and arbitrary
direction of in the presence of
arbitrary current with angle of .
By substituting
from
Fi ( , , )
Eq.(18), Eq.(19) and Eq.(20) into
Eq.(21), the mean force components in
S ( )d = m
of
2U sin g a2
F6 = CT 1 +
g
2
2
2
1
L L + L L sin
sin
r
1
2
2
1
3
(20)
2 2U cos
gr cos d
F1s = 2 S ( ) CT 1 +
3
g
SI
2
rL1 3 sin 3 sin cos +
2
L2 L1 sin sin + rL 2 sin + 2 sin cos
2
2
3
3
m1 2
=
m0
T1
H s = 4 m0
H s2
4
C T 1 + Fn1
16
T1
L 4
r cos
g 3
(22)
H s2
4
CT 1 + Fn 2
T1
16
(23)
13/E
9/E
21/E
www.SID.ir
Iranian Association of
Naval Architecture & Marine Engineering
H s2
4
C T 1 + F n 2
16
T1
r sin L1 L 2 +
2
2
L 2 L1 sin sin +
r L1 + L 2 sin 2
F 6s = g
(24)
L1
F1s
Kx
(25)
F2s
=
Ky
(26)
F6s
K
(27)
6s =
Ar
ch
ive
where:
1s , 2s and 6s are steady drift
displacement in surge, sway and yaw
modes respectively.
F1s , F2s and F6s are steady drift loads in
irregular waves in surge, sway and yaw
modes respectively.
K x and K y are mooring stiffness
coefficients in x and y directions
respectively.
K is the rotational mooring stiffness
which for a turret moored FPSO can be
defined as:
2
K = Lm K x
where
L m is the distance between the turret
mooring point and FPSOs centre of
gravity.
3.5.WAVE DRIFT DAMPING
4
CT a2 r cos
3
(29)
(30)
C66 = CT a2
3 r sin (L1 L2 ) +
2
2
L2 L1 sin sin +
1 r (L1 + L2 )sin 2
of
s
2
C11 =
(28)
SI
1s =
C ii = C T a2 sin 2 ( + )ni dl
(31)
www.SID.ir
10
B66*sqrt(L/g)/[rho*B*(0.5H)2]
B11*sqrt(L/g)/[rho*B*(0.5H)2]
8
7
6
5
4
3
1000
500
2
1
0
5
6
w*sqrt(L/g)
10
10
8
6
4
10
ive
15
20
25
w*sqrt(L/g)
Ar
ch
120
100
B22*sqrt(L /g)/[rho*B*(0.5H)2]
20
4. CONCLUDING REMARKS
The wave direction, wave frequency
and waterline shape of a floating
structures were considered to obtain the
second-order mean wave loads on
moored floating. In this regards,
modification was made to Faltinsens
wave drift formula by adding finite draft
coefficient. The modified formula can
take into account the effects of wave
direction, frequency and the waterline
shape of floating structures. Meanwhile,
the obtained formula was applied to an
FPSO.. The results obtained from this
formula which is dependent on wave
frequencies, were compared with
Helvacioglus experiments [2] favorably
in high wave frequencies. Moreover, the
influence of the current on the wave
mean drift force was taken into account
by considering the current coefficient
derived from the ship added resistance
formula. It is found that, the presence of
current can increase the mean drift force
by up to 50 percent for a particular range
of wave frequencies and, the amount of
increase for a floating structure is
independent of its size.
For a 200,000 tdw FPSO, the mean drift
loads in irregular waves with H s =5 m
and T1 =16 sec are 5-15 percent of the
mean drift loads of the mentioned FPSO
when subjected to a regular wave with
of
B11*sqrt(L/g)/[rho*B*(0.5H)2]
12
80
60
40
20
15
w*sqrt(L/g)
SI
10
14
10
15
20
25
w*sqrt(L/g)
13/E
11/E
23/E
www.SID.ir
Iranian Association of
Naval Architecture & Marine Engineering
condition
Irregular waves
H s =5 m
ch
T1 =16 sec
U =1 m/s
Fn1 =0
Fn1 =0.02
Fn1 =0.013
Fn1 =0
FN 2 =0
FN 2 =0
FN 2 =0.033
FN 2 =0.047
1067
1044
988.8
123.6
133.42
130.54
123.6
=0
988.8
=0
U =1 m/s
ive
Regular waves
H =5 m
T =16 sec
Surge
(KN)
Sway
(KN)
Yaw
(KN.m)
Surge
(KN)
Sway
(KN)
Yaw
(KN.m)
U =0 m/s
=45
Ar
Current
Condition
Sea
condition
Regular waves
H =5 m
T =16 sec
Irregular waves
H s =5 m
T1 =16 sec
12/E
24/E
14/E
Surge
(KN)
Sway
(KN)
Yaw
(KN.m)
Surge
(KN)
Sway
(KN)
Yaw
(KN.m)
=90
SI
Sea
U =0 m/s
of
Current
Condition
= 0)
= 45 )
U =0 m/s
U =0 m/s
U =1 m/s
U =1 m/s
Fn1 =0
FN 2 =0
Fn1 =0.02
FN 2 =0
Fn1 =0.013
FN 2 =0.033
Fn1 =0
FN 2 =0.047
699.19
754.71
738.44
699.16
3865
3865
4082
4172
192340
192340
203150
207620
87.4
94.34
92.3
87.4
483
483
510.3
521.5
11683
11683
12339
12611
=0
=0
=45
=90
www.SID.ir
Regular
waves
H =5 m
T =16 sec
Irregular
waves
H s =5 m
T1 =16 sec
Surge
(KN)
Sway
(KN)
Yaw
(KN.m)
Surge
(KN)
Sway
(KN)
Yaw
(KN.m)
U =0 m/s
U =1 m/s
U =1 m/s
Fn1 =0
FN 2 =0
Fn1 =0.02
FN 2 =0
Fn1 =0.013
FN 2 =0.033
Fn1 =0
FN 2 =0.047
7320.4
7320.4
7731
7902
188620
188620
199220
203600
915.05
915.05
966.4
987
=0
11310
3. REFERNCES
=45
11310
11945
=90
12208
Ar
ch
ive
of
=0
Sea
condition
U =0 m/s
SI
Current
Condition
= 90 )
13/E
1/E
25/E
www.SID.ir
Iranian Association of
Naval Architecture & Marine Engineering
Ar
ch
ive
of
SI
26/E
2/E
14/E
www.SID.ir