Chapter 5 - Seakeeping Models

Chapter 5 – Seakeeping Models
5.1 Hydrodynamic Concepts and Potential Theory The study of ship dynamics has traditionally
been covered by two main theories:
5.2 Seakeeping and Maneuvering Kinematics maneuvering and seakeeping.
5.3 The Classical Frequency-Domain Model
5.4 Time-Domain Models including Fluid Memory Effects Maneuvering refers to the study of ship motion
5.5 Identification of Fluid Memory Effects in the absence of wave excitation (calm water).
The maneuvering equations of motion are
derived in Chapter 6 under the assumption that
the hydrodynamic potential coefficients and
radiation-induced forces are frequency
independent (constant).
Seakeeping refers to the study of motion of

marine craft on constant course and speed
when there is wave excitation. This includes the
trivial case of zero speed. In sea keeping
analysis, a dissipative force known as fluid
memory effects (Cummins 1962) is introduced
due to frequency-dependent potential
coefficients.
1
Lecture Notes TTK 4190 Guidance, Navigation and Control of Vehicles (T. I. Fossen)
Chapter Goals
• Understand seakeeping coordinates {s} and how they relate to BODY and NED
• Understand frequency-dependent hydrodynamic matrices and their application
• Application of Cummins equation to transform frequency-dependent matrices to the time domain
Computed using hydrodynamic seakeeping codes
Transformation from SEAKEEPING {s}

to BODY axes {b}
Linear model used for

control design and simulation
2
Chapter 5 - Seakeeping Models
Equations of Motion
Seakeeping theory is formulated in equilibrium (SEAKEEPING) axes {s} but it can be transformed to BODY axes {b}
by including fluid memory effects represented by impulse response functions.
The transformation is done within a linear framework such that additional nonlinear damping must be added in the
time-domain under the assumption of linear superposition.
μr is an additional term representing the fluid memory effects. -In the absence of wave excitation, {s} coincides with {b}.
- Under the action of the waves, the hull is disturbed
from its equilibrium and {s} oscillates, with respect to its
equilibrium position.
3
5.1 Hydrodynamic Concepts and Potential Theory
Strip Theory (2-D Potential Theory)
For slender bodies, the motion of the fluid can be formulated as a 2-D problem. An accurate estimate of the hydrodynamic
forces can be obtained by applying strip theory (Newman 1977; Faltinsen 1990; Journee and Massie 2001).
The 2-D theory takes into account that variation of the flow in the cross-directional plane is much larger than the variation in
the longitudinal direction of the ship.
The principle of strip theory involves dividing the submerged part of the craft into a finite number of strips. Hence, 2-D
hydrodynamic coefficients for added mass can be computed for each strip and then summed over the length of the body to
yield the 3-D coefficients.
There exists more than 20 commercial strip theory programs and several opens source university codes.
4
5 ShipX (VERES) - VEssel RESponse program is a Strip Theory Program which calculates wave-induced loads on and motions of mono-hulls and barges in deep to very shallow
water.Lecture
The program is based
Notes TTK 4190 on the famous
Guidance, paper and
Navigation by Salvesen,
Control ofTuck and(T.
Vehicles Faltinsen (1970). “Ship Motions and Sea Loads”. Transactions of SNAME.
I. Fossen)
Panel Methods (3-D Potential Theory)
For potential flows, the integrals over the fluid domain can be transformed to integrals over the boundaries of the fluid
3D Visualization of the Wamit file: supply.gdf
domain. This allows the application of panel or boundary element methods to solve the 3-D potential theory problem.
-2
Z-axis (m)
-4
3-D panelization of
a supply vessel
-6
-8
-10
-12
-40
-30
-20
30
-10
20
0 10
10 0
20 -10
30 -20
40 -30
X-axis (m) Y-axis (m)
Panel methods divide the surface of the ship and the surrounding water into discrete elements (panels). On each of these
elements, a distribution of sources and sinks is defined which fulfill the Laplace equation.
Commercial code: WAMIT (www.wamit.com)

6
WAMIT
WAMIT® is the most advanced set of tools available for analyzing wave interactions with offshore platforms and other
structures or vessels.
WAMIT® was developed by Professor Newman and coworkers at MIT in 1987, and it has gained widespread recognition
for its ability to analyze the complex structures with a high degree of accuracy and efficiency.
Panelization of semi-submersible using WAMIT user supplied tools
Over the past 30 years WAMIT has been licensed to more than 100 industrial and research organizations worldwide.
7
Potential theory programs typically compute:
• Frequency-dependent added mass, A(w)
• Potential damping coefficients, B(w)
• Restoring terms, C
• 1st- and 2nd-order wave-induced forces and motions
(amplitudes and phases) for given wave directions and frequencies
• … and much more
BUT THE CODE GIVES YOU NO VISCOUS DAMPING TEMRS We need to add this manually (Section 5.3)
One special feature of WAMIT is that the program solves a boundary value problem for zero and infinite added mass. These
boundary values are useful when computing the retardation functions describing the fluid memory effects.
Processing of hydrodynamic data:
Matlab MSS Toolbox: http://github.com/cybergalactic/MSS
The toolbox reads output data files generated by the hydrodynamic programs:
• ShipX (Veres) – www.sintef.com
• WAMIT – www.wamit.com
and processes the data for use in Matlab/Simulink.
8
MSS Hydro
ShipX Veres and WAMIT data sets in the toolbox
9
MSS Hydro
Load WAMIT data into Matlab vessel structure
10
MSS Hydro
Load ShipX Veres data into Matlab vessel structure
11
MSS Hydro
Viscous damping is added manually

(see Section 5.3)
12
MSS Hydro
13
Next Step: Transform the Hydrodynamic Model from
SEAKEEPING to BODY
Computed using hydrodynamic

seakeeping codes
Transformation from SEAKEEPING) axes {s} to BODY axes {b}

Using kinematic transformation (Sections 5.2-5.3)

14
5.2 Seakeeping and Maneuvering Kinematics
Seakeeping Theory (Perturbation Coordinates)
The SEAKEEPING reference frame {s} is not fixed to the craft; it is
fixed to the equilibrium state
! ! !! 1 , ! 2 , ! 3 , ! 4 , ! 5 , ! 6 " T #
! ! !"
!" ! !#" ! !$
-In the absence of wave excitation, {s} coincides with {b}.

! sb ! !! 4 , ! 5 , ! 6 " ! ! !"#, "#, "$" ! # - Under the action of the waves, the hull is disturbed
from its equilibrium and {s} oscillates, with respect to its
Transformation between {b} and {s} equilibrium position.
0 0 0 0 0 0 v nns ! !U cos !, U sin!, 0" ! # ! 0 $!

!! ! ! !U!L!" " e 1 " #
0 0 0 0 0 1 " ! 0 # $" #
!!# ! !# !UL! # ! nns ! !0, 0, 0" ! #
0 0 0 0 !1 0
# #
" $#
L :! " "!
" ns ! !0, 0, ! #
0 0 0 0 0 0
e 1 ! !1, 0,0,0, 0, 0" !
0 0 0 0 0 0
0 0 0 0 0 0
15
5.3 Cummins Equation
Seakeeping Analysis
The seakeeping equations of motion are assumed to be inertial.
Equations of Motion
! ! !" ! !!x, !y, !z, !", !#, !$" ! #
MRB !" ! #hyd " #hs " #exc #
Cummins (1962) showed that the radiation-induced hydrodynamic

forces in an ideal fluid can be related to frequency-dependent added t
mass A(ω) and potential damping B(ω).

!hyd ! !Ā"# ! " K
$ !t ! !""%!!"d! #
0
Ā ! A!!"
"
2
K!t" ! ! ! 0 B!"" cos!"t"d" #
Cummins, W. E. (1962). The Impulse Response Function and Ship Motions. Technical Report 1661. David Taylor Model Basin. Hydrodynamics Laboratory, USA.
16
5.3 Frequency-Dependent Potential Coefficients
Frequency-dependent added mass A22(ω) and potential damping B22(ω) in sway
Ā ! A!!"
17
5.3 Retardation Functions and Fluid Memory Effects
t
!hyd ! !Ā"# ! " K
$ !t ! !""#!!"d! #
0
Matrix of retardation functions

"
2
K!t" ! ! ! 0 B!"" cos!"t"d" #
Cummins Equation
If linear restoring forces τhs = -Cξ are included in the model, this
results in the time-domain model:
t
!MRB ! A!!""!" ! " K
# !t # !"!$!!"d! ! C! " %exc #
0
The fluid memory effects can be replaced by a state-space model

to avoid the integral, see Section 5.5
In mathematics, an integro-differential equation is an equation that involves both integrals and derivatives of a function
18
5.3 Forced Oscillations
Wave excitation
In an experimental setup with a restrained scale model, it is possible to vary the wave excitation frequency ω and the
amplitudes fi for i = 1, .., 6 of the excitation force. Hence, by measuring the position and attitude vector ξi the response of the
second-order system can be fitted to a linear model (e.g. by least squares) for each frequency ω
This gives two frequency-dependent model matrices A(ω) and Btotal(ω)
19
5.3 Frequency-Dependent Hydrodynamic Coefficients
The matrices A(ω), Btotal(ω) and C represent a hydrodynamic mass– damper–spring system which varies with the
frequency of the forced oscillation. The added mass matrix should not be understood as additional mass due to a finite
amount of water that is dragged with the vessel
20
5.3 Incorrect Mixture of Time and Frequency in
Hydrodynamic Models
Hydrodynamic mass—damper—spring system
This model is rooted deeply in the literature of hydrodynamics and this false time-domain model has been discussed
eloquently in the literature since the model mixes time and frequency incorrectly (abuse of notation). Hence, it is not an ODE.
An Ordinary Differential Equation (ODE) is an equation with a function and one or more of its derivatives. It is a function of
only one independent variable, typically the time.
Consequently, the correct approach is to integrate Cummins time-domain model to find the solution for all frequencies
NB! If the matrices A(ω) and Btotal(ω) are computed at one and only one constant frequency (e.g. the zero frequency which is
used in maneuvering theory) the hydrodynamic mass damper-spring-system will be an ODE, which can be solved by numerical
integration.
21
5.3 Viscous Damping
When running seakeeping codes it is important to include an external viscous damping matrix BV(ω)
in order to obtain accurate estimates of the vessel responses.
BV(ω) can be chosen as a diagonal matrix
where
α > 0 is the exponential rate

βi > 0 (i = 1, 2, 6) are linear viscous skin friction coefficients in surge, sway and yaw
NITTC(A1) is equivalent linear surge resistance as a function of the surge velocity amplitude A1
βIKEDA(ω) is frequency-dependent roll damping based on the theory of Ikeda et al. (1976)
(other models for viscous roll damping can also be used)
One useful property of the exponential function is that linear skin friction only affects low-frequency motions.
Ikeda, Y., K. Komatsu, Y. Himeno and N. Tanaka (1976). On Roll Damping Force of Ship: Effects of Friction of Hull and Normal Force of Bilge Keels.
Journal of the Kansai Society of Naval Architects 142, 54–66.
22
5.3 Viscous Damping (cont.)
Surge, sway and yaw

Specify three time constant T1, T2, and T6 satisfying
This gives
Heave and pitch

Specify a percentage increase in damping e.g. β3 = p B33(ω3) where p > 0.
Then solve for β3 and β5 in
23
5.3 Quadratic Surge Resistance
Quadratic surge resistance can be approximated using a describing function
(similar to the equivalent linearization method)
y ! N!A"u # Sinusoidal input u ! Asin!!t" #
For static linearities, displaying no dependence upon the derivatives, the describing function for the
particular odd polynomial nonlinearity
Quadratic damping approximation
y ! c 1 x " c 2 x|x|"c 3 x 33 # X ! !X |u|u |u|u
" N ITTC !A1 "u #
is (Gelb and Vander Velde 1968)
2
N!A" ! c 1 " 8A c 2 " 3A c 3 #
3! 4
This assumes that the surge velocity

is harmonic with amplitude A1, i.e.
For a ship moving at high speed, the amplitudes A2 and A6 will be much smaller than A1.
Hence, it is common to neglect these terms in seakeeping analysis.
24
MSS Hydro
Plots the hydrodynamic matrices A, B, C

and Bv as a function of frequency
25
5.3 Potential Coefficients
Longitudinal added
mass coefficients as
a function of
frequency.
26
5.3 Potential Coefficients (cont.)
Lateral added mass

coefficients as a
function of frequency.
27
Longitudinal potential
damping coefficients as a
function of frequency.
Exponential decaying viscous

damping is included for B11.
28
Lateral potential damping

coefficients as a function
of frequency.
Exponential decaying
viscous damping is
included for B22 and B66
while viscous IKEDA
damping is included in B44
29
high speed,
at high the
For
speed, amplitudes
athe
ship
amplitudes 2 and
moving atAhigh willamplitudes
A6 the
speed, be muchAsmaller than
2 and A6than
will be much smaller than
A 2 and A6 will be much smaller
mon to neglect these terms in seakeeping analysis.
A1 . Hence,
ommon to neglect it is
these 5.3 Response Amplitude Operator (RAO)
common
terms to neglect these
in seakeeping terms in seakeeping analysis.
analysis.
amplitude 5.3.3 Response amplitude operators

operators
nse amplitude operators
Equation (5.56) can be transformed to the frequency domain by using the Laplace trans-
en transformed
be transformed to to
formation. Laplace
the frequency
the frequency
Moreover, domain
domainby
application ofusing
by using ¨ the=
L{⇠(t)} Laplace
Laplace trans-
⇠(s) trans-
andFluid
s2including L{⇠}˙ = s⇠(s) together
114
114 ¨ ¨ Time-Domain
Time-Domain˙˙ Models
Models including Memory
Fluid Memory Effects
Effects
, application
over, withofs of
application
114 =L{L{
j!⇠(t)}
2 2
⇠(s)and
and
⇠(t)}==sTime-Domain
gives s⇠(s) L{
L{⇠}
Models ⇠} = s⇠(s)
including s⇠(s) together
together
114 Time-Domain Models including Fluid Memory Effects
AssumeAssume
that
exc⌧ (that
Assume that
2⌧
!j!) F((⇣j!)
=⌧exc
[M exc j!)
+=
for i=A(!)]
i⇣ for
⇣2,for
=F1, ...,i6ij!B
=1,1,
=are 2,2,...,
...,(!)
harmonic are
6 6are harmonic
harmonic
excitation
+ C ⇠( j!)excitation
excitation
forces propor- ( forces
= ⌧ forces propor-
j!)propor- (5.73)
exc
114tional totional
tional
an totoananincoming
incoming regular
i
RB
incoming wave ⇣ = ⇣Time-Domain
regular
regular wave
awave
e j!t ⇣where
total j!t
⇣==⇣a⇣⇣eaaej!t
is the Models
where
where
wave⇣aamplitude including
⇣isa the
is wave
the andwave Fluidand
amplitude
Fi amplitude Memory
Fand
i Fi Effects
MRB + A(!)]
+ A(!)]denotes j!B j!B
denotes
denotes
the the (!)
the
proportional
total (!) + +C
proportional
proportional
gain. C ⇠(
Linear ⇠(
gain.
gain. j!)
j!)
theory =
Linear
Linear=theory
implies ⌧⌧
theory
that
exc j!)
j!)
exc((⇠implies
i implies
j!t ⇠ =(5.73)
= ⇠¯i e thatthat i ⇠i (5.73)
where ⇠⇠¯¯
ii e
= ¯j!t
denotes
⇠ e where
j!t ¯i denotes
⇠
where ⇠¯ denotes
total
Harmonic Assume
excitations that ⌧exc ( j!) = Fi ⇣ for i = Regular
i 1, 2, ..., wave 6i are harmonic excitation forces propor-
thetheamplitudes.
the amplitudes. amplitudes.
Consequently, Consequently,
Consequently, j!t
Assume that C⌧exc!(2 [M j!) tional toforan iincoming j!tregularare wave
¯eharmonic ⇣⇣ e= j!t⇣a e where a is the wave amplitude and Fi
forces ⇣propor-
=+!22A(!)]
F[M i ⇣RB = 1, (!)2, a excitation
j!B j!B¯...,
⇠e 6=(!)Fi ⇣a⇠e j!t
j!t
= Fi(5.74) j!t (5.74)
CCRB + A(!)]
total
j!B total ¯ j!t (5.74) ⇠i F=i ⇠¯i e j!t where ⇠¯i denotes
tional to an incoming denotes the proportional gain. Linear ⇣a is thetheory wave implies that
! [M RB + A(!)] j!t
total (!) ⇠e = F i ⇣ a e
regular wave ⇣ = ⇣a e where amplitude and
denotesLet theRAOproportional
i (!)RAO
Let denote theresponse
the
i (!) denote
amplitudes.
gain. Linear
responsetheory
theamplitude Consequently,
operator
amplitude implies
between ⇣a that
operator and ⇠¯i⇠for
between i = i⇣a=⇠¯ j!t
i e ⇠¯i for
and i = ⇠¯i denotes
where
Let
1, 2, ..., 6.1,Hence, RAOtheHence,
decoupled
(!) denote the response
transfer functions amplitude
become operator between ⇣ and ⇠¯i for i =
the amplitudes. 2, ...,
RAO Consequently,
6. i
from
the
wave
decoupled
amplitude
transfer functions
z functions
to 2response
become
x
a
¯ j!t = Fi ⇣a e j!t
1, 2, ..., 6. Hence,
⇠¯i
the decoupled transfer
CFi a ! [MRB +i A(!)]
become j!Btotal (!) ⇠e (5.74)
RAOi (!) = =2 ⇠¯i 2 Fi (5.75)
j!t j!t
CRAO ⇣a! i (!)[M
C=ii ⇠¯i! =[M + A(!)]+2Aii (!)]j!B + jBtotal,ii (!)
Fi (!)! ¯
⇠e = Fi ⇣a(5.75)
e (5.74)
RAOi (!) = ⇣RB
a = CRB,ii
ii ! [MRB,ii + A total
ii (!)] + jBtotal,ii (!)!(5.75)
⇣a Cii ! and
Notice that RAOi (!) is a frequency-dependent
2[MRB,ii +A
complex ii (!)] +
function. jBcommon
It is total,ii (!)!
Notice that RAO (!)Let
is a RAO (!) denote
frequency-dependent
i and
to only consider the absolute value of the response amplitude operator
i the response
complex function. amplitude
It is common operator between ⇣a and ⇠¯i for i =
D RAFT M AN
Let RAO
Notice Amplitude
to only consider
that
(!) denote
RAO the
1, and
(!) thephase
absolute value
response of the response
amplitude
frequency-dependent amplitude
and operator
i 2, ..., 6. Hence, the decoupled transfer functions
is a complex operator
between
function. It is ⇠¯i for i =
andbecome
⇣acommon
D RAFT M A
i Fi
to
|RAO only
1, 2, ..., 6. Hence, consider
= q the absolute value of the response amplitude operator
(5.76)
|RAOthe decoupled transfer2 functions become
i (!)| Fi
D RAFT M
(Cii = !q
i (!)| 2 [M
RB,ii + Aii (!)]) + (!B¯ total,ii (!))2
⇠ (!)]) (5.76)
Fi Note the similarity to
Bode magnitude and phase
+ AiiiF
2 2
|RAOthe
The phase between q (CiiRAO
⇠¯i=excitation
i (!)|
wave
! 2 [MRB,ii=
and thei (!)
i
ship motions is =Fi 2
+ (!Btotal,ii (!))
plots, for which magnitude
(5.76) (5.75)
is logarithmic and
RAO (!) between
The iphase = =
the (C
wave ! 2 [MRB,ii
ii✓ excitation
⇣ ship
theAa
and + ii (!)])
C
◆ ii ! 2 [M(!))2
is total,iiRB,ii
+ (!B
motions
+ Aii (!)](5.75)
+ jB
given in
total,ii (!)!
decibels while phase is plotted in degrees
⇣a= atan Cii ! 2!B [M total,ii (!)+ Aii (!)] + jBtotal,ii (!)!
\RAOi (!) ✓ RB,ii (5.77)
◆ using a common logarithmic frequency axis.
The phase between the wave Cexcitation
ii and the
! 2 [MRB,ii Aship
+ !B ii (!)]motions
(!) is
Notice i (!)that is a frequency-dependent (5.77)and complex function. It is common
total,ii
\RAO = atan RAO i (!)
30 Note the similarity to Bode magnitude and phase ✓Ciiplots, 2
! for
[MRB,ii Aii (!)] is ◆
which+magnitude
Notice that
logarithmic
RAO
and Notesi (!)
given
isonly
to
in \RAO
aGuidance,
decibels
frequency-dependent
consider
while phase is theandabsolute
plotted in
!Bof
degrees
and complex
value
total,ii (!)(T. of function.
the responseIt is common operator
amplitude
plots, afor
using common (5.77)
Lecture TTK 4190 Navigation Control Vehicles I. Fossen)
D
(!) = atan
Note the similarityi to Bode magnitude and phase 2 which magnitude is
D
5.4 Time-Domain Models including Fluid Memory Effects
Linear SEAKEEPING model based on Cummins equation

(linear theory where frequency-dependent matrices are replaced by
fluid memory effects)
Linear equations of motion expressed in BODY

Transform from SEAKEEPING to BODY coordinates using the linear
kinematic transformation
Unified maneuvering and seakeeping model expressed in BODY

(Replace the linear kinematic and Coriolis terms with their nonlinear
counterparts. Add nonlinear viscous damping and maneuvering coefficients)
31
5.4 Cummins Equation in SEAKEEPING Coordinates
Cummins (1962) Equation
t
!MRB ! Ā"!" ! ! # !t " !"!$!!"d! ! C
K # ! " %wind ! %wave ! "% #
"#
The Ogilvie (1964) Transformation gives
Ā ! A!!" #
"
! !t" ! 2
K ! ! 0 B total!"" cos!"t"d" #
From a numerical point of view is it better to integrate the difference

"
2
K!t" ! ! ! 0 #B total!"" # B total!""$ cos!"t"d" #
This can be don by rewriting Cummins equation as

t
!MRB ! A"!#$!" ! B total"!#!# ! " K"t # !#!#"!#d! ! C! " $wind ! $wave ! "$ #
0
32
5.4 Linear Time-Domain Seakeeping Equations in
BODY Coordinates
t
!MRB ! A"!#$!" ! B total"!#!# ! " K"t # !#!#"!#d! ! C! " $wind ! $wave ! "$ #
0
It is possible to transform the time-domain representation of Cummins equation from {s} to {b} using the kinematic
relationships
!! ! ! !U!L!" " e 1 " # ! ! !"

!!# ! !# !UL! # !" ! !#
This gives:
t
!MRB ! A"!#$!!" !UL!$ ! B total"!#!! !U"L!# " e 1 #$ ! # K"t " "#!!""#d" ! C!# " $wind ! $wave ! "$ " $%# #
0
The steady-state control force τ needed to obtain the forward speed U when τwind = τwave= 0 and δη = 0 is
!" ! B total!!"Ue 1 #
Hence,
t
!MRB ! A"!#$!!" !UL!$ ! B total"!#!! ! UL!#$ ! " K"t # "#!!""#d" ! C!# " $wind ! $wave ! $ #
0
33
BODY Coordinates
t
!MRB ! A"!#$!!" !UL!$ ! B total"!#!! ! UL!#$ ! " K"t # "#!!""#d" ! C!# " $wind ! $wave ! $ #
0
When computing the damping and retardation functions, it is common to neglect the influence of δη on the forward
speed such that
!! ! v !U!L!" " e 1 " ! v "Ue 1 #
Finally, let use replace ν by the relative velocity νr to include ocean currents and define M = MRB + MA such that
t
M!" ! C !RB ! ! C !A !r ! D!r ! " K!t # !"#!!!"#Ue 1 $d! ! G# " $wind ! $wave ! $ #
0
where
MA ! A!!"
C "A ! UA!!"L Linear Coriolis and
C "RB ! UMRB L centripetal forces due
D ! B total!!" to a rotation of {b} about {s}
G!C
34
BODY Coordinates
The integral in the following equation represents the fluid memory effects
7 K (t)
x 10 22
2.5
2
t
! :! ! 0 K!t " !"#"!!""Ue 1 $ d! # 1.5
""
1
0.5
Approximated by a state-space model

0
-0.5
! ! H!s"#" !Ue 1 $ # -1
0 5 10 15 20 25
tim e (s)
Impulse response function

x! " A rx # B r!!
# "
" " C rx 2
K!t" ! ! ! 0 #B!"" # B!""$ cos!"t"d" #
35
5.4 Nonlinear Unified Seakeeping and Maneuvering
Copyright © Bjarne Stenberg
Model with Fluid Memory Effects
Linear Seakeeping Equations (BODY coordinates)
Copyright © The US Navy
Unified Nonlinear Seakeeping and Maneuvering Model
• Add ocean currents (relative velocity)

• Use nonlinear kinematics
• Replace linear Coriolis and centripetal forces with their nonlinear counterparts
• Include maneuvering coefficients in a nonlinear damping matrix (linear superposition)
36
MSS Simulink Library, Models and Templates
Simulink demo files for marine craft
simulation, wave generation, state
estimation, etc.
User editable Simulink templates for DP and

maneuvering.
Motion and force RAOs are used in DP, while

the unified maneuvering model can only be
used with force RAOs
All models have fluid memory effects
The MSS Simulink library is loaded by

double-clicking mssSimulink.slx, which
contains many Simulink blocks that can be
37 used in own applications.
DP Motion RAO Template
Double click to load

vessel and
vesselABC data
structures
38
Simulink Demo: DP System for 6-DOF Semisubmersible
39
Simulink Demo: Wave Spectrum to Wave Elevation
40
5.5 Case Study: Identification of Fluid Memory Effects
The fluid memory effects can be approximated using frequency-domain identification. The main tool for this is the MSS
FDI toolbox (Perez and Fossen 2009)
When using the frequency-domain approach, the property that the mapping: has relative degree one is
exploited. Hence, the fluid memory effects μ can be approximated by a matrix H(s) containing relative degree one
transfer functions:
p r sr "p r!1 sr!1 "..."p 0

h ij !s" !
sn "q n!1 sn!1 "..."q 0
r ! n ! 1, n " 2
State-space model:
H!s" ! C r!sI ! A r" !1 B r #
41
5.5 Frequency-Domain
Identification using the
MSS FDI Toolbox
42
5.5 Frequency-Domain Identification using the MSS FDI
Toolbox 3.2
3
10
8 DoF 33
Ac(jw)
3.2
3
10
8 DoF 33
A
Achat(jw) order 4 Aest
2.8 2.8
FPSO identification results for h₃₃(s) 2.6 2.6
without using the infinite added mass
Added Mass
2.4 2.4
|Ac(jw)|
A₃₃(∞). The left-hand-side plots show 2.2 2.2
the complex coefficient 2 2
1.8 1.8
1.6 1.6
1.4 1.4
1.2 1.2
-2 -1 0 1 -2 -1 0 1
10 10 10 10 10 10 10 10
and its estimate while added mass A(jw) Freq. [rad/s] Frequency [rad/s]
and potential damping B(jw) are plotted

on the right-hand-side. 5 5
10
7
B
0 Best
4
-5
Phase Ac(jw) [deg]

-10 3
Damping
-15
2
-20
-25 1
-30
0
-35
-40 -1
-2 -1 0 1 -2 -1 0 1
10 10 10 10 10 10 10 10
Frequency [rad/s] Frequency [rad/s]
43
Chapter Goals - Revisited
• Understand seakeeping coordinates {s} and how they relate to BODY and NED
• Understand frequency-dependent hydrodynamic matrices and their application
• Application of Cummins equation to transform frequency-dependent matrices to the time domain
Computed using hydrodynamic seakeeping codes
Transformation from SEAKEEPING {s}

to BODY axes {b}

44

Chapter 5 - Seakeeping Models

Uploaded by

Document Informationclick to expand document information

Document Informationclick to expand document information

Copyright:

Available Formats

Chapter 5 - Seakeeping Models

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Chapter 5 - Seakeeping Models

Uploaded by

Copyright:

Available Formats

Chapter 5 – Seakeeping Models

Seakeeping refers to the study of motion of

Computed using hydrodynamic seakeeping codes

Transformation from SEAKEEPING {s}

Linear model used for

Commercial code: WAMIT (www.wamit.com)

Panelization of semi-submersible using WAMIT user supplied tools

Processing of hydrodynamic data:

Matlab MSS Toolbox: http://github.com/cybergalactic/MSS

ShipX Veres and WAMIT data sets in the toolbox

Viscous damping is added manually

Computed using hydrodynamic

Transformation from SEAKEEPING) axes {s} to BODY axes {b}

Linear model used for

-In the absence of wave excitation, {s} coincides with {b}.

0 0 0 0 0 0 v nns ! !U cos !, U sin!, 0" ! # ! 0 $!

! ! !" ! !!x, !y, !z, !", !#, !$" ! #

MRB !" ! #hyd " #hs " #exc #

Cummins (1962) showed that the radiation-induced hydrodynamic

mass A(ω) and potential damping B(ω).

Matrix of retardation functions

The fluid memory effects can be replaced by a state-space model

This gives two frequency-dependent model matrices A(ω) and Btotal(ω)

BV(ω) can be chosen as a diagonal matrix

α > 0 is the exponential rate

Surge, sway and yaw

Heave and pitch

y ! N!A"u # Sinusoidal input u ! Asin!!t" #

This assumes that the surge velocity

Plots the hydrodynamic matrices A, B, C

Lateral added mass

Exponential decaying viscous

Lateral potential damping

amplitude 5.3.3 Response amplitude operators

Linear SEAKEEPING model based on Cummins equation

Linear equations of motion expressed in BODY

Unified maneuvering and seakeeping model expressed in BODY

The Ogilvie (1964) Transformation gives

From a numerical point of view is it better to integrate the difference

This can be don by rewriting Cummins equation as

!! ! ! !U!L!" " e 1 " # ! ! !"

!! ! v !U!L!" " e 1 " ! v "Ue 1 #

Approximated by a state-space model

Impulse response function

Copyright © The US Navy

Unified Nonlinear Seakeeping and Maneuvering Model

• Add ocean currents (relative velocity)

User editable Simulink templates for DP and

Motion and force RAOs are used in DP, while

All models have fluid memory effects

The MSS Simulink library is loaded by

Double click to load

p r sr "p r!1 sr!1 "..."p 0

H!s" ! C r!sI ! A r" !1 B r #

FPSO identification results for h₃₃(s) 2.6 2.6

without using the infinite added mass

the complex coefficient 2 2