Seakeeping: Complex numbers

Compute:
(3 4i ) ( 2 2i )
(3 4i ) ( 2 2i )
(3 4i ) (2 2i )
3 4i

2 2i
Euler's formula brings trigonometry and exponential function together:

eix cos x i sin x

where i 1 .
eix A sin( x ) , where A
is the complex amplitude with real part
Show for one example that Re A

) / Im( A
) .
Re() and imaginary part Im() and absolute value A = A and the phase shift is tan Re( A
= 3+4i.
Take as an example A

Solution
We use elementary mathematical operations and get:
(3 4i ) ( 2 2i ) 5 2i
(3 4i ) ( 2 2i ) 1 6i
(3 4i ) (2 2i ) 6 8i 6i 8 14 2i
3 4i 3 4i 2 2i 2 14i

0.25 1.75i
2 2i 2 2i 2 2i
8
Compute the phase shift:
)
Re( A
Re(3 4i )
3
arctan
arctan
arctan 36.87

Im(3 4i )
4
Im( A)
The (real valued) amplitude is:

A A

) 2 Im( A
)
Re( A

32 4 2 5

Re A eix Re (3 4i )(cos x i sin x) 3 cos x 4 sin x

and

A sin( x ) 5 sin( x 36.87) 5(sin x cos 36.87 sin 36.87 cos x) 4 sin x 3 cos x

The two expressions are equal.

Seakeeping: Celerity of wave

A surf board travels with Froude number Fn = 0.5 on a deep-water wave of 4 m length. How long must the
surf board be to ride on the wave, i.e. to have the same speed as the wave?
Solution
The speed of the surf board is: V Fn gL
g
2
Both velocities must be the same for surf riding:

The celerity of the wave is: c

V c Fn gL

g
2

2 .5 m
2
2 Fn 2 0.52

Seakeeping: Velocity and pressure in elementary wave

Consider a regular deep-water wave of wave length =100 m and wave height 5 m.
a) What is the orbital velocity in 10 m depth?
b) What is the amplitude of the pressure fluctuation (absolute and in percent of the hydrostatic
pressure)? The pressure fluctuation is (linearised): p f t
Solution
a) The amplitude of the wave is half the wave height: h = 2.5 m.
2
2

0.0628 m1
The wave number is: k

100
The frequency is: g k 9.81 0.0628 0.785 s1
The orbital velocity is the maximum of the x-velocity (for deep waves). The x-velocity is:

v x x Re h e kz e i (t kx )

Thus the amplitude is:

v x h e kz 0.785 2.5 e 0.062810 1.05 m/s
b) The pressure fluctuation (=pressure without hydrostatic part) is given by:

p f t g Re h e kz e i (t kx )

The amplitude of the pressure fluctuation is then:

p f gh e kz 1025 9.81 2.5 e 0.062810 13.4 kPa

The hydrostatic pressure is

p stat gz 1025 9.81 10 100.6 kPa
Thus the pressure fluctuation is 13% of the hydrostatic pressure.

Seakeeping: Velocity and acceleration in shallow-water wave

The potential of a regular wave on shallow water is given by:

Re ich sinh( kH ) cosh(k ( z H ))ei (t kx )

wave length
= 100 m
wave amplitude h = 3 m
water depth
H = 30 m
Determine velocity and acceleration field at a depth of z = 20 m below the water surface!
The following parameters are given:

Solution
The velocity is derived by differentiation of the shallow-water potential:

h
h
vx x Re
cosh(k ( z H ))ei (t kx )
cosh( k ( z H )) cos(t kx)
sinh( kH )
sinh( kH )

ih
h
vz z Re
sinh( k ( z H ))ei (t kx )
sinh( k ( z H )) sin(t kx)
sinh( kH )
sinh( kH )

This derivation used the relation c = /k. The individual values are:
2
2
k

0.06283 m1

100
gk tanh( kH ) 9.81 0.06283 tanh(0.06283 30) 0.767 s1
This yields:
vx 0.861 cos(0.767t 0.06283 x )

vz 0.480 sin(0.767t 0.06283 x)

The accelerations are obtained by differentiating the velocities with respect to time:
a x xt 0.661 sin(0.767t 0.06283 x )
a z zt 0.368 cos(0.767t 0.06283 x )
The velocities have a phase shift of 90 and different amplitudes. The particles will thus trace an ellipse.

Seakeeping: Wave breaking

Wave breaking occurs theoretically when particle velocity in x-direction is larger than the celerity c of the
wave. In practice, wave breaking occurs for wave steepness h/ exceeding 1/14. What is the theoretical and
practical limit for h concerning wave breaking in a deep-water wave of = 100 m?
Solution
Practically the limit is h = /14 = 100/14 = 7.143 m.
2
2

0.0628 m1
The wave number is: k

100
kh
The maximum particle velocity in x-direction is: vx , max h e
The celerity is: c

1
e kh
kh
The above expression is solved iteratively for kh, using kh e kh . The iteration starts with the practical
limit: k0h = 7.143 0.0628 = 0.4488, then yields:
The condition c = vx,max yields:

1
0.4488
11
0.5667

2
0.6384
12
0.5674

3
0.5281
13
0.5670

kh 0.5671

4
0.5897
14
0.5672

5
0.5545
15
0.5672

0.5671
9.026 m
0.0628

6
0.5744
16
0.5671

7
0.5631
17
0.5671

8
0.5695
18
...

9
0.5658

10
0.5679

Seakeeping: Wave with two buoys

The position of two buoys is given by their (x,y) coordinates in [m] as sketched below. The buoys are excited
by a regular wave of = 62.8 m, amplitude h = 1 m, and angle = 30 to the x-axis.
a) What is the maximum vertical relative motion between the two buoys if they follow exactly the
waves?
b) What is the largest wave length to achieve maximum vertical relative motion of twice the wave
amplitude?

Solution
2
2

0.1 m1

62.8
Transform coordinates in local -system. Points on = const. have same -values for regular waves.
The two buoys have then coordinates:
1 = 0
2 = x2 cos 30 y2 sin 30 = 30 cos 30 10 sin 30 = 21 m
The relative difference between the wave elevations is:

a) The wave number of this wave is: k

1 2 Re h ei (t k1 ) Re h ei (t k 2 ) Re h eit 1 e ik 2 )

The amplitude (= maximum) of this relative motion is given by:

h e it 1 e ik 2 ) h
1

1 cos(k 2 ) 2 sin( k 2 ) 2
1 cos(0.1 21) 2

sin(0.1 21)

= 1.73 m
b) The maximum difference between wave elevations is twice the wave amplitude if the buoys are
spaced by an odd multiple of half the wave length. The longest wave is obtained for a spacing of
half the wave length:

21 m

= 42 m
2
All other wave lengths fulfilling the criterion are even shorter.

Seakeeping: Encounter frequency

A ship travels with 28.28 kn in deep sea in regular sea waves. The ship travels east, the waves come from
southwest. The wave length is estimated to be between 100 m and 200 m. The encounter period Te is
measured at 31.42 s.
a) What is the length of the seaway?
b) Two days before, a storm started in an area 1500 km southwest of the ship's position.
Can the waves have their origin in this storm area?

Solution
The speed is V = 28.28 0.5144 m/s = 14.55 m/s.
The encounter period yields the encounter frequency:
The value of is:

2
2

0.2 s1
Te
31.42

eV 0.2 14.55

0.2966
g
9.81

There are three possible frequencies which could excite this encounter frequency:

g
9.81
1 1 4 cos
1 1 4 0.2966 cos 45 1.123 Hz
2V cos
2 14.55 cos 45
g
9.81
2
1 1 4 cos
1 1 4 0.2966 cos 45 0.6684 Hz
2V cos
2 14.55 cos 45
g
9.81
3
1 1 4 cos
1 1 4 0.2966 cos 45 0.2854 Hz
2V cos
2 14.55 cos 45

The corresponding wave lengths are:

2g 2 9.81
1 2
49 m
1
1.1232

2g 2 9.81

138 m
22
0.66842
2g 2 9.81
3 2
756 m
3
0.28542

Only 2 = 138 m fits the observed bandwidth of wave lengths.

Wave groups travel with group velocity. For deep ocean water, this is
1 g
1 9.81 138
cgr

7.34 m/s
2 2
2
2
At this speed, waves can travel within two days:
s = 48 h 60 min/h 60 s/min 7.34 m/s = 1268 km < 1500 km
The waves can thus not originate from the storm area.

Seakeeping: Undamped free heave oscillation

The differential equation of a free undamped heave oscillation is:

m m33 z cz 0

m is the displacement of the ship, m33 the added mass for heave motion. c is the restoring force coefficient.
The (circular) natural frequency is subsequently:
c
z
m m33
Consider the ship with the following particulars:
L = 125 m, B = 17 m, T = 7 m, CB = 0.75, Cwp = 0.8, m33 = 0.9m.
a) Determine the (circular) natural frequency and natural period of heave!
b) Initial conditions are: At t = 0 we have z = 0 and z =0.5 m/s. What is the maximal load a winch of
4000 kg mass exerts on its foundation?
Solution
a) The natural frequency follows from:
g Cwp L B
g Cwp
c
z

m m33
CB L B T (1 m33 / m)
CBT (1 m33 / m)

9.81 0.8
0.89 Hz
0.75 7 (1 0.9)

2
2

7.08 s
z
0.89
b) The winch exerts its maximal load when gravity g and amplitude of heave oscillation are
superposed. We write the heave motion as:
Tz

z Re K ei z t

K K1 iK 2 is determined by the initial conditions:

zt 0 Re K ei z 0 Re K K1 0

z t 0 Re i z K ei z 0 K 2 z 0.5

K2

Then the amplitude of the heave acceleration is given by:

z z2 K 2 0.89 2 0.56 0.45 m/s2
The maximal load of the winch is then:
mwinch ( g z) 4000 (9.81 0.45) 41 kN

0.5
0.5

0.56 m
z
0.89

Seakeeping: Power requirements for a wave maker

A wave maker is to be designed for a towing tank of width B = 4 m and depth H = 2.5 m. The wave maker
shall be designed for a wave of 5 m length and 0.2 m amplitude.
a) What is the power requirement for the motor for the wave maker if we assume 30% total efficiency
between motor and wave? (For the considered wave length, the depth can be regarded as deep.
The power requirement of the wave is (energy/meter [in direction of wave propagation]) * group
velocity, as the energy in a wave is transported with group velocity.)
b) After switching the wave maker off, for a long time there is still a wave motion with period 40 s
observed in the tank. How long is the tank if the motion is due to the lowest natural frequency of
the tank?
Solution
a) The group velocity is:
1
1 g
1 9.81 5
cgr c

1.397 m/s
2
2 2
2
2
The average energy per area is:
1
1
E gh 2 1000 9.81 0.2 2 196.2 N/m
2
2
Thus the power of the wave is:
P E B cgr 196.2 4 1.397 1096 W
The power of the motor then needs to be:
Pm P / 0.3 3.65 kW
b) The wave length is now long compared to the depth. Thus we have to use finite water depth
expressions:
2
k tanh( kH )

2
0.157 2
0.157 Hz, this yields:
k tanh( kH )
0.002515 m1
T
9.81
This equation has to be solved iteratively. As the convergence is slow, it is useful to start with a good
estimate. For small x, tanh(x) x. This would yield k 0.002515 / 2.5 0.03172 .
We solve the problem using Newton's iteration:
f (k )
k k
f ' (k )
With f k tanh( 2.5k ) 0.002515 and f ' tanh(2.5k ) k / cosh 2 (2.5k ) , we obtain within
7 iterations:
2
2
k 0.03175

198 m.
k
0.03175
With

The lowest natural frequency is at /2 = L, i.e. the tank has a length of 99 m.

Seakeeping: Wave maker for required wave length and amplitude

A wave maker consists of a cylinder which moves up and down at a wall. The wave maker has a cross
section of half a Lewis section with Cm = 0.8.
The draft of the section is T = 1 m, the half-width is B/2 = 0.6 m. Determine the amplitude and frequency of
motion to create waves of 5 m length and 0.2 m amplitude.
Note: The ratio between wave amplitude and body motion amplitude is Az .
Solution
The frequency corresponding to = 5 m is:
2g
2 9.81

3.51 Hz

5
Important parameters for Lewis section curves:
2 B 3.512 0.6
2g

9.81

0.754

B
0.6

0.6
2T
1
The curves for Lewis sections give Az = 0.52. Thus:
h
0.2
u3

0.38 m
Az
0.52
H

Seakeeping: Wave maker with Lewis section moves floating body

A wave maker consists of a cylinder which moves up and down at a wall. The wave maker has a cross
section of half a Lewis section with Cm = 0.8. The draft of the section is T = 1 m, its half-width B/2 = 0.6 m.
The wave maker performs heave motions with 0.6 m amplitude and = 3.51 Hz. The basin is 10 m wide.
The mass of the wave maker itself is negligibly small. Use linear theory in your solution approach.
a) What is the amplitude of the force needed to move the wave maker?
b) What is the average power supply if an efficiency of 70% is assumed?
c) In the middle of the basin is a cylinder of Cm = 0.8, B = 1.2 m, T = 1 m, L = 5 m, arranged with axis
parallel to the wave maker axis. The cylinder has all six degrees of freedom suppressed. What is the
exciting force on the cylinder?
d) If the cylinder in the middle of the basin is free to heave, what is its heave amplitude?
Assume deep water.
Solution
a) The basic equation for heave motion is:

m m
2

33

in33 c33 u3 f

(The force is negative as now the body excites the water and not the water the body.) All quantities
are in this equation per length, i.e. for the final computation they need to be multiplied by the tank
width. The individual quantities are:
2 B 3.512 0.6

2g

9.81

0.754

B
0.6

0.6
2T
1
The curves for Lewis sections give Cz = 0.6, Az 0.52 . Thus:
B2
1.22
m33 C z
0.6 1000
339 kg/m
8
8
H

2
g 2
2 1000 9.81
n33 A 3 0.52
601.8 kg/ms

3.513
2
z

c33 gB 1000 9.81 1.2 11770 N/m2

The mass of the wave maker is negligible: m = 0
Thus:

3.51 (0 339) i 3.51 601.8 11770 0.6 f

2

f (4556 1267i ) N/m

This force is for waves radiated in both directions, i.e. so far we assumed symmetry. The total force
amplitude is then given by:
1
1
| F | | f | BTank
45562 1267 2 10 23.6 kN
2
2
b) The power needed to drive the wave maker is given by:
E BTank cgr gh 2
c 1
P

BTank kN

2
2
The wave amplitude is:
h | u3 | Az 0.6 0.52 0.312 m
The celerity is:
g 9.81
c

2.79 m/s
3.51
Thus:

1000 9.81 0.3122

2.79 1
10

9.5 kW
2
2 0.7
Now we check the result, using a different approach. The power needed is due to the damping force
times the velocity. The mass force and the restoring force are phase shifted by 90 to the velocity
and thus do not contribute to the power. Again we take into account that we have only half a cross
section:
P

1 Im( f )

BTank
2
2

( z )

1 1 1267
1

10 (3.51 0.6)
9.5 kW
2 2
0. 7

c) We read from the diagrams for Lewis sections:

Re fe
0.25 f e , r 0.25 gBh 0.25 1000 9.81 1.2 0.312 918.2 N/m
gBh
Im fe
0.25 f e, i 918.2 N/m
gBh

This then gives the absolute value for the three-dimensional force:
Fe f e2,r f e2,i L 918.2 2 918.2 2 5 6.5 kN
d) The basic equation for heave motion is:

m m
2

33

in33 c33 u3 fe

The mass (per length) follows from Archimedes:

m Cm B T 1000 0.8 1.2 1.0 960 kg/m
As before, we have then:

3.51 960 339 i 3.51 601.8 11770 u

2

4233.8 2112 .3i u3 (918.2 918.2i)

u3

0.087 2 0.260 2 ) 0.274 m

(918.2 918.2i )

u3 ( 0.087 0.260i ) m

Seakeeping: Cylinder with Lewis cross section in regular waves

A cylinder of Lewis cross section (L = 10 m, B = 1 m, T = 0.4 m, Cm = 0.8) floats parallel to the wave crests
of regular waves = 5 m, h = 0.25 m) which excite heave motions. Assume that the form of the free surface
is not changed by the cylinder. What is the amplitude of relative motion between cylinder and free surface
predicted by linear theory?
Solution
Fundamental equation for heave motion:

m m
2

33

in33 c33 u3 fe

The individual quantities needed are:

2g
2 9.81

3.51 Hz

5
2 B 3.512 1
2g

2 9.81

0.628

B
1

1.25
2T 2 0.4
*
*
Then the curves for Lewis sections yield: Cz = 0.62, Az 0.68 , f er 0.39 , f ei 0.38 . This in turn
gives:
B2
12
m33 C z
0.62 1000 243 kg/m
8
8
m Cm B T 1000 0.8 1.0 0.4 320 kg/m
H

n33 Az2

2
g 2
2 1000 9.81

0
.
68

1029 kg/ms
3
3.513

c33 gB 1000 9.81 1.0 9810 N/m2

f e gB f er* f ei*i h 1000 9.81 1.0 (0.39 0.38i ) 0.25 956.5 932.0i N/m
We assume = 1000 kg/m3 here for convenience. As appears in all terms, the final result is not influenced
by the choice of .
This yields the heave motion amplitude:

3.51 243 320 i 3.51 1029 9810 u

2

956.5 932.0i

956.5 932.0i
0.287 0.0366i m
2874 3612i
The amplitude of relative motion is determined from:
u3

r u3 h 0.287 0.0366i 0.25 0.037 0.0366i

0.037 2 0.0366 2 = 0.052 m

Seakeeping: Lewis section in forced heave

An infinite cylinder of width B = 2 m, draft T = 1 m and with Lewis cross section of coefficient Cm = 0.8
floats in equilibrium in fresh water. Then a harmonic force per length with period Te = 3.14 s and amplitude
fe=1000 N/m is applied. What motion results after a long time (when the initial start-up has decayed)?
What wave length and what wave amplitude are generated, assuming a deep basin?
Solution
The motion has only one degree of freedom, namely heave. The basic equation is then:

m m
2

33

in33 c33 u3 fe

u3 is the complex amplitude of heave motion. The individual quantities needed are:
2
2

2 Hz
Te
3.14
2 B
22 2

0.408
2g
2 9.81
B
2
H

1
2T 2 1
Then the curves for Lewis sections yield: Cz = 0.7, Az 0.51 . This in turn gives:

B2
22
0.7 1000
1100 kg/m
8
8
m Cm B T 1000 0.8 2.0 1.0 1600 kg/m

m33 C z

n33 Az2

2
g 2
2 1000 9.81

0
.
51

3129 kg/ms
3
23

c33 gB 1000 9.81 2.0 19620 N/m2

This yields the heave motion amplitude:

2 1600 1100 i 2 3129 19620 u

2

1000

1000
0.0754 0.0535i m
8820 6258i
The amplitude is:
u3 0.07542 0.05352 0.0925 m
The wave length follows from:
2 2
2g

2 15.4 m
g

The wave amplitude follows from Az :

h Az u3 0.51 0.0925 0.047 m
u3

Seakeeping: Power generator

A cylinder for power generation is moored in regular waves such that the cylinder axis is parallel to the wave
crests. The cylinder dimensions are L = 50 m, B = 10 m, T = 5 m, the cross section is a Lewis section with
Cm = 0.8. The waves have = 50 m and h = 1 m. The mooring adds a restoring force coefficient c = 2.5106
N/m and a damping (from the generator) of d = 1.5106 kg/s. The mooring remains at all times under tension
as the buoyancy of the cylinder is much larger than its weight due to its mass of 10 6 kg. The sea water has a
density of 1025 kg/m3.
What is the amplitude of heave motions for the cylinder?
What would be the amplitude without generator?
Solution
The basic equation of motion is:

m m
2

33

L i (d n33 L) (c c33 L) u3 Lf e

The individual quantities needed are:

2g
2 9.81

1.11 Hz

50
2 B 1.112 10
2g

2 9.81

0.628

B
10

1
2T 2 5
*
*
Then the curves for Lewis sections yield: Cz = 0.61, Az 0.65 , f er 0.39 , f ei 0.36 . This in turn
gives:
B2
102
L m33 L C z
50 0.61 1025
1.2275 106 kg
8
8
H

L n33 L Az2

2
g 2
2 1025 9.81

50

0
.
65

1.5235 106 kg/s

3
1.113

L c33 L gB 50 1025 9.81 10 5.0275 106 N/m

L f e L gB f er* f ei*i h 50 1025 9.81 10 (0.39 0.38i ) 1 (1.961 1.81i ) 106 N

Then the basic equation yields:

1.11 1 1.2275 i1.11 (1.5 1.5235) (2.5 5.0275) u

2

1.961 1.81i

1.961 1.81i
0.453 0.061i m
4.783 3.356i
The amplitude is:
u3 0.4532 0.0612 0.46 m
The basic equation modifies for the case of no generator (i.e. d = 0):
u3

m m
2

33

L i (n33 L) (c c33 L) u3 Lf e

1.11 1 1.2275 i1.11 1.5235 (2.5 5.0275) u

2

1.961 1.81i
0.483 0.208i m
4.783 1.691i
u3 0.4832 0.2082 0.53 m

u3

1.961 1.81i

Seakeeping: Catamaran in waves

A raft consists of two cylinders with D = 1 m diameter and L = 10 m length. The two cylinders have a
distance of 2e = 3 m from centre to centre. The raft has no speed and is located in regular waves coming
directly from abeam. The waves have = 3 m and wave amplitude h = 0.5 m. Assume the two cylinders to
be hydrodynamically independent, i.e. waves created by one cylinder are not reflected at the other cylinder.
The raft has a draft of 0.5 m, centre of gravity in the centre of the connecting plate, radius of moment of
inertia for rolling is kx = 1 m. The density of the water is 1000 kg/m3. The Cm = /4 is close enough to 0.8 to
use the Lewis section curves for this Cm. Neglect sway motion. What is the maximum roll angle?
Hint: The centre of gravity of a semi-circle is 4r/(3) from the flat baseline, where r is the radius.

Solution

D 2
12
L 2 1000
10 2 7854 kg
8
8
2
2
The mass moment of inertia is: m k x 7854 1 7854 kgm2
2 2
k

0.6667 m1
Wave number:

3
2g
2 9.81
Wave frequency:

2.56 s1

3
The metacentric height is:
4
D 3 L / 12 DLe 2 2

GM KB BM KG 1
D
T
3
D 2 L / 4

The mass of the raft is:

4
13 10 /12 1 10 1.52 2 1 5.23 m

1
0 .5
3
12 10 / 4

B
1. 0
2 B 2.56 2 1.0

0.333

1
Parameters for Lewis curves:
and
2g
2 9.81
2T
2 0.5
*
*
Then the curves for Lewis sections yield: Cz = 0.8, Az 0.4 , f er 0.6 , f ei 0.25

2 ( m44 ) in44 c44 u4 f e 4

The basic equation of roll motion is:
The hydrodynamic added mass m44 and damping n44 (moments) are derived from the added mass for heave
(forces) times lever e:
m44u4 2 m33 u3 e

The factor 2 is due to the two cylinders. The term in parentheses is the vertical force. With u3 e u4 , we
get for the whole 3-d raft:
B 2

1.02
m44 2
C z L e 2 2 1000
0.8 10 1.52 14137 kgm2
8
8

g 2
1000 9.812
10 1.52 41300 kgm2/s
n44 2 Az2 3 L e 2 2 0.4 2
3

2
.
56

The restoring force constant for roll motion is:

c44 mg GM 7854 9.81 5.23 403 kNm

The exciting moment has to consider the phase shift in the wave between the two hulls. The force on the left
floater is then:
fl f er* if ei* gBh L e ike
The force on the right floater is:
fr f er* if ei* gBh L eike
The phase shift is thus in the last term. Combine the two forces with lever e to get:

f e 4 e f er* if ei* gBh L e ike eike

1.5 0.6 0.25i 1000 9.81 1.0 0.5 10 2i sin(0.6667 1.5)

= (30956 74293 i) kgm2/s2
Then our basic equation becomes:

2.56 (7854 14137) i 2.56 41300 40300 u

(30956 74293i )
258.9 105.7i u4 (30.96 74.29i) u4 0.002 0.288i
2

The (real) roll amplitude is then

u4

0.002 2 0.2882 0.288 16.5

The assumption of linearity is thus no longer justifiable.

Seakeeping: Strip method roll

A circular cylinder (length L = 20 m, D = 1 m, radius of moment of inertia kxx = 0.4 m, mass m = 8000 kg)
floats in water of density = 1019 kg/m3. The centre of gravity is 0.1 m under the centre of the circle. The
cylinder is excited by regular waves with wave crests parallel to the cylinder axis. The waves have amplitude
h = 0.15 and length = 3.14 m. The following hydrodynamic properties are known:
a) hydrodynamic mass for sway: 0.4 m
b) hydrodynamic damping for sway given by ratio of amplitudes Ay 0.3
c) real part of horizontal exciting force: 0.8 orbital velocity of wave times damping coefficient as
given in b)
d) imaginary part of horizontal exciting force: 0.6 wave steepness times cylinder weight
The coordinate system should have its origin in the centre of the circle. Then the resultant force of the
hydrodynamic forces passes through the origin.
Determine the complex amplitude of roll motion!
Solution
We use a coordinate system with z pointing down, y pointing right in direction of wave propagation. x is then
the cylinder axis.
The displacement of the fully submerged cylinder is:
D 2
12
L
20 1019 16000 kg
4
4
As the cylinder mass is 8000 kg, the draft of the cylinder is T = D/2 = 0.5 m. The frequency is:
2g
2 9.81

4.43 Hz

3.14
The orbital velocity involves differentiation with respect to time (factor ), the steepness with respect to y
(factor k). The hydrodynamic quantities hidden in a) to d) are then:
m22 0.4 m 0.4 8000 3200 kg

n22 Ay2

2
g 2
2 1019 9.81

0
.
3

20 2040 kg/s
3
4.433

Re Fe 0.8 h n22 0.8 4.43 0.15 2040 1084 N

Im Fe 0.6 k h g m 0.6 2 h m 0.6 4.432 0.15 8000 14130 N

We can use the formulae for strip method and simplify them. As the waves come just from abeam, no
symmetric motions (heave, surge, pitch) are excited (with linear theory). Also, as the cylinder has always
the same cross section, no yaw motion is excited. The fundamental equation of motions is:
2 ( M A) iN S u Eh
Sway motion does not excite roll motion, neither does roll motion excite sway for a circular cylinder with
origin as chosen. There is no exciting roll moment either as in potential theory only normal pressures exist
and for the circular cylinder all normal vectors pass through the origin, i.e. excite no moment. Similarly,
there is no radiation moment, i.e. no damping and added mass for roll. Sway has no restoring term. The
metacentre of a circular cylinder is the centre of the circle. Thus GM = 0.1 m. Thus we can write for the two
degrees of freedom:

m
mz g

m22
0

mz g
1
m

z
zg

0
3200 0
m

0
0
0

z g 8000 800

k xx2 800 1280

 n
N 22
0
0
S
0

0
2040 0
m

0
0
0
0
0
m

mgGM
0

0
7848

Thus:

219800 15700 9037i

15700
25120 0

219800 9037i 15700 u2

15700
17272 u4

0 u2 1080 14130i
0

0 0 7848 u4
0

1080 14130i

The solution of this system of equations is:

u2
0.002 0.069i

0.002 0.063i
u4
Thus the roll amplitude is u4

0.0022 0.0632 = 0.063 = 3.6.

Seakeeping: Pontoon with crane in waves

Consider a pontoon with a heavy-lift derrick as sketched. The pontoon has L = 100 m, B = 20 m, D = 10 m,
mp = 107 kg. The load at the derrick has mass ml = 106 kg. The height of the derrick over deck is 10 m. This is
where the load can be considered to be concentrated in one point. The longitudinal position of the derrick is
20 m before amidships. A force F = 10 6 N acts on the forward corner. We assume homogeneous mass
distribution in the pontoon.
a) Consider the pontoon in air (without hydrodynamic masses) and determine the acceleration vector
!
u
b) Consider the pontoon statically in water and determine u!

Solution
All numbers are given in standard units. We use the coordinate system as in the book with z pointing down.
Origin is at K.
0
0

0
0

106
F

a) The general 6-component force vector is: F

7
F 10m
10
F 50m
5 107

0
0

The mass matrix needs the following expressions:

m m p ml 107 106 1.1 107
m z g m p z g , p ml z g , l 107 ( 5) 106 ( 20) 7 107

m xg m p xg , p ml z g ,l 107 0 106 20 2 107

The moments of inertia are computed for pontoon and load:

L
m
B3
B 2 107
20 2
B D 3 D
p D 2

102

xx , p ( z 2 y 2 ) d m p
3
4
3
4
3
4
= 0.667 109
xx, l ml z g2 , l 106 20 4 0.4 109
Thus

xx xx , p xx , l 1.067 109 . Correspondingly we get:

xz , p

xz d m 0

xz , l ml x g , l z g , l 106 ( 20) 20 0.4 109

Thus

xz xz , p xz , l 0.4 109 .
yy , p ( x 2 z 2 ) d m

= 8.667 109

p B
m
L3
L2
107
1002
L D 3 D p D 2
10 2

3
4
3
4
3
4

yy ,l ml xg2 ,l z g2 ,l 106 (202 202 ) 0.8 109

Thus

yy yy , p yy , l 9.467 109

zz , p ( x 2 y 2 ) d m

p D
m
107

202 1002
LB 3 BL3 p B 2 L2
12
12
12

= 8.667 109

yy ,l ml xg2 , l y g2 ,l 106 (202 02 ) 0.4 109

Thus

zz zz , p zz , l 9.067 109

0
0
0
70
0
11
0
11
0
70
0
20

0
0
11
0
20
0
6
This yields the mass matrix: M 10

70
0
1067
0
400
0
70 0 20
0
9467
0

20
0
400
0
9067
0
Now we can use the fundamental equation F M u , where F and u are generalised 6-component
vectors. Since we solve manually, it is advisable to decouple the 66 system of equations into 33
systems for symmetric and anti-symmetric degrees of freedom:

1
11u
1
70u

3
11u
3
20u

2
11u
2
70u

4
70u
4
1067u

2
20u

4
400u

b) Basic equation of motion is:

5 0
70u
5 1
20u
5 50
9467u

6 0
20u
6 10
400u
6 0
9067u

1
u
3
u
5
u

0.036

0.101
0.005

2
u
4
u
6
u

(M A) i N S u F
2
e

Now we consider the quasi-static case, i.e. e = 0. F is then the exciting force. Then: S u F
m
1.1 107
3
5.5 m,
The hydrostatic data are: Aw L B 100 20 2000 m2, xw = 0, T
Aw 10 2000
GM = KB + BM KG =

T
B2

z g 2.75 + 3.06 6.36 = 2.45 m,

2 12T

T
L2

z g 2.75 151.51 6.36 147.9 m

2 12T
0
0
0
0 0
0 0
0
0
0

0 0 19.62
0
0
S 106
This gives:
0
264.4
0
0 0
0 0
0
0
15960

0
0
0
0 0
It is impossible to make any statement on u1, u2, and u6. The
simple system of equations:
19.62u3 1 u3 0.05 m
264.4u4 10 u4 0.038
15960u5 50 u5 0.0031
GM L

0
0

0
0

0
remaining degrees of freedom form a

Seakeeping: Probability of exceeding operational limit

A ship sails in a natural seaway of approximately Te = 6 s encounter period between ship and waves. The
bridge of the ship is located forward near the bow. During one 1 hour, 6 times a downward acceleration
occurred exceeding gravity acceleration g. For downward accelerations exceeding 1.5g, severe injuries for
the crew have to be expected. What is the probability for this happening if the ship continues sailing with the
same speed in the same seaway for 12 hours?
Solution
The number of amplitudes encountered per hour (=3600 s) is:
3600 3600
n

600
Te
6
The actual number of amplitudes for acceleration exceeding g in one hour was 6. Thus the probability that
one amplitude exceeding g occurs is W(1g) = 600/6 = 0.01. On the other hand, we can write this probability
formally:

W (1g ) e (1g )

/( 2 m0 )

0.01 m0 0.109 g 2

The probability for an amplitude exceeding 1.5g is correspondingly:

W (1.5 g ) e (1.5 g )

/( 2 m 0 )

e 1.5

/( 2 0.109 )

3.29 105

The number of amplitudes encountered in 12 h is 12 n 12 600 7200. Thus the probability for one
amplitude exceeding 1.5g in 12 h is:
W 7200 3.29 10 5 = 0.24 = 24 %

Seakeeping: Roll motion of cargo ship in waves

A new cargo ship type features relatively large beam and extreme sectional flare at the ship ends to allow
many container on deck. The ship has the following main data: LWL = 94m, B = 16 m, D = 8.40 m, T = 5 m,
radius of inertia (including added mass) k' = 0.4 B = 6.4 m. The service speed is V = 11.8 kn.
Investigate whether (respectively when) there is a danger of capsizing due to resonance in wave induced
rolling for possible metacentric heights (0.30 m < GM < 1.00 m). Consider the following scenarios:
a) waves from abeam ( < 150 m)
b) following waves ( < 150 m)
If there is a danger of capsizing, specify the conditions which are most critical. Assume deep water.
Solution
a) waves from abeam
Resonance appears for e = n. For elementary waves on deep water this yields:

2g
g GM
2 k ' 2 2 6.4 2 257.4

GM
GM
GM
k '2
Thus for 0.30 m < GM < 1.00 m only wave lengths 257 m < < 858 m lead to resonance. For the
given range below 150 m, there is no danger of capsizing.
b) following waves
Resonance appears for

n
1
n with n = 1, 2, 3, ...
e
2

For following waves (angle of encounter = 0) we have: e

With

2g /

we get: e

2
V
g

2g 2V

The first resonance lies at n = e . The natural frequency of roll is n

Thus we get for the first resonance:

g GM
1

k'
2

g GM
k'

2g V

k ' 2
2 3 k ' 2 V 2 k ' 2 V 2

2
g
g 2
3
Take all quantities in standard units (m and s), to get with V = 11.8 kn = 6.07 m/s:
6 .4 2
2 3 6.4 2 6.07 2 6.4 2 6.07 2 64.34 625.11 1518.3
GM

2
9.81

9.81 2
2
3
3
This is evaluated for various wave lengths in the given interval:
Solving for GM yields: GM

30
40
50
60
80 100 120
[m]
GM [m] 0.03 0.09 0.13 0.15 0.17 0.17 0.17

150
0.16

The case GM > 0.30 m is never reached. Thus there is no danger of capsizing for n = 1.
The second resonance lies at n = e. Due to the quadratic relation, we have thus four times the
values for GM:

 [m]
GM [m]

30
40
50
60
80 100 120
0.11 0.35 0.50 0.60 0.67 0.68 0.66

150
0.62

Thus for most wave lengths, there is a danger of resonance for smaller GM values.
The third resonance lies at n = 1.5 e.
[m]
GM [m]

30
40
50
60
80 100 120
0.11 0.35 0.50 0.60 0.67 0.68 0.66

150
0.62

The third resonance is also reached, but only for wave lengths much shorter than the ship length.
Thus the expected fluctuation in righting levers will be smaller as for second resonance.
Correspondingly the probability of capsizing will be lower.
Higher resonances will shift resonance to even shorter waves which will be uncritical.

Seakeeping: Roll motion of submarine in waves

Consider a surfaced submarine with speed V = 6 kn. Neglect the sail and model the submarine hull as a
circular cross section. Neglect damping.
a) What is the expected roll amplitude in elementary waves from abeam?
b) Evaluate the danger of parametric roll excitation in elementary following waves!
Ship data: L = 70.00 m, D = 8.00 m (diameter), T = 7.00 m, KG = 3.60 m, k' = 3.00 m (radius of the
transverse mass moment of inertia incl. added mass)
Wave data (deep water): = 70 m (wave length); h = 1.75 m (wave amplitude)
Solution:
a) waves from abeam
2g
2 9.81

0.93835 s1

70
For a circular cylinder, the metacentre is always in the centre of the circle:
D
8.00
GM KM KG
KG
3.60 0.40 m
2
2

Wave circular frequency:

The natural frequency for roll is: n

For

g GM

k '2

9.81 0.40
0.6603 s1
3.00 2

e 0.93835

1.4211 the roll excitation lies in the super-critical region. The roll amplitude
n
0.6603

for undamped roll is thus (with k = 2/):

1
1
2
| || u4 |
k h

1.75 0.15406 8.83

2
2
1 (e / n )
1 1.4211 70
b) following waves
Encounter frequency e is then:
2
2g 2V
2 9.81 2 6 0.5144
e
V

0.6613 s1
g

70
70
Although the roll natural frequency n is close to encounter frequency e, there is no parametric
excitation, because for a circular cylinder there is no fluctuation of righting levers. The righting lever
is always h = GM sin , with GM = 0.40 m = const.

Seakeeping: Cruise ship in waves from abeam

A small cruise ship is designed. As the destination island has no deep-water port, the passengers will
disembark via small boats. To allow normal disembarkation, the metacentric height should be always
sufficiently high to have maximum roll amplitude = 6 in regular waves from abeam up to wave steepness
max = 4.5 up to wave length of max = 60 m. Assume a linear restoring moment. For the small roll angles
considered, damping can be neglected.
The radius of inertia including added mass is k' = 4 m.
What is the resulting requirement for the metacentric height?
Solution
The requirement 6 for max = 4.5 means for the amplification function: Vratio

1.33
max 4.5

The equation for V yields then limiting ratios of frequency:

(e / n ) 2 0.25
1
Vratio

1
.
33

(e / n ) 2 1.75
1 (e / n ) 2
2
Thus the frequency interval to be avoided is 0.25 (e / n ) 1.75 . With the relations for a deep water
regular wave:
2g
e

and for the natural frequency of roll

g GM
n
k '2
we get for the interval to be avoided:
2 k '2
0.25
1.75
GM
For various wave lengths, we can then determine the interval to be avoided:

20 m
40 m
60 m

interval to avoid
2.87 m GM 20.11 m
1.44 m GM 10.05 m
0.96 m GM 6.70 m

One can quickly see that short wave lengths pose no problem. The limiting value is GM 0.96 m. This
requirement will hardly be possible for a passenger ship due to damaged stability requirements.