Description

A Matlab Toolbox for Parametric Identification of Radiation-Force Models of Ships and Offshore Structures

Categories

Published

All materials on our website are shared by users. If you have any questions about copyright issues, please report us to resolve them. We are always happy to assist you.

Share

Transcript

ARC Centre of Excellence for Complex Dynamic Systems and Control, pp. 1–22
A Matlab Tool for Frequency-DomainIdentiﬁcation of Radiation-Force Models of Shipsand Oﬀshore Structures
Technical Report: 2009-02.0-Marine Systems SimulatorSeptember 2009
Tristan Perez
1
,
3
Thor I. Fossen
2
,
3
1
Centre for Complex Dynamic Systems and Control (CDSC), The University of Newcastle, Callaghan,NSW-2308, AUSTRALIA. E-mail:
Tristan.Perez@newcastle.edu.au
2
Department of Engineering Cybernetics, Norwegian University of Science and Technology, N-7491 Trond-heim, Norway. E-mail:
Fossen@ieee.org
3
Centre for Ships and Ocean Structures (CeSOS) Norwegian University of Science and Technology, N-7491Trondheim, Norway.
The University of Newcastle, AUSTRALIA MSS/09/02
ARC Centre of Excellence for Complex Dynamic Systems and Control–CDSC
Contents
1 Introduction 42 Dynamics of Ships and Oﬀshore Strucutres 53 Frequency-domain Models 54 Identiﬁcation of Radiation-force Models 6
4.1 Identiﬁcation when
A
∞
is Avalaible . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74.2 Order Selection . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84.3 Stability . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84.4 Passivity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94.5 Identiﬁcation when
A
∞
is not Availaible . . . . . . . . . . . . . . . . . . . . . . . . . . 9
5 Toolbox Description 10
5.1 FDIRadMod.m . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105.2 EditAB.m . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 125.3 Ident retardation FD.m . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 135.4 Ident retardation FDna.m . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 135.5 Fit siso fresp.m . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
6 Demos 137 Software Repository 148 Conclusion 14A Parameter Estimation Algorithm 16
2
Technical Report
Replacements
This is a new technical report.
Executive summary
This article describes a Matlab tool for parametric identiﬁcation of radiation-force models of marinestructures. These models are a key component of force-to-motion models used in simulators, motioncontrol designs, and also for initial performance evaluation of wave-energy converters. The softwaredescribed provides tools for preparing the non-parmatric data generated hydrodynamic codes andidentiﬁcation with automatic model-order detection. The identiﬁcation is considered in the frequencydomain.3
ARC Centre of Excellence for Complex Dynamic Systems and Control–CDSC
1 Introduction
One approach to develop linear time-domain models of marine structures consist of using potential-theory hydrodynamic codes to compute frequency-dependent coeﬃcients and frequency responses,and then use these data for system identiﬁcation in order to implement the Cummins equation, whichis a linearised vector equation of motion. If physical-model or full scale experiments are available,then mmathematical model based on the Cummins equation can be corrected for viscous eﬀects. Thisprocedure us summarised in Figure 1.A great deal of work has been reported in the literature proposing the use of diﬀerent identiﬁ-cation methods to obtain approximating fuild-memory models. Taghipour et al. (2008) and Perezand Fossen (2008b) provide an up-to-date review of the diﬀerent methods. In particular, the latterreference discusses the advantages of using frequency-domain methods for the identiﬁcation of ﬂuidmemory models. Since the data provided by hydrodynamic codes is in the frequency domain,identiﬁcation in the frequency domain is a natural approach, which does not require transformationof the data to the time domain. If not handled appropriately, the latter transformation can resultin errors due to the ﬁnite amount of frequency-domain data. More importantly, when performingfrequency-domain identiﬁcation, one can enforce model structure and parameter constraints; andthus, the class of models over which the search is done is reduced, and the models obtained presentcharacteristics in agreement with the hydrodynamic modelling hypothesis.In this article, we present a set of Matlab functions to perform identiﬁcation of radiationforces. We consider two cases. In ﬁrst case, information related to the inﬁnite-frequency added masscoeﬃcients is considered available. In the second case, these coeﬃcients are estimated jointly with theﬂuid-memory model (Perez and Fossen, 2008a). The second case is relevant for hydrodynamic codesbased on 2D-potential theory, which do not normally solve the boundary-value problem associatedwith the innite frequency.
!"#$%#"&'()* ,%#- .#-&/0*'/%& ,1(()&2 341'/%& 356-$)(-&72 8%#-9 :)7; <)2*%12 ,%$$-*/%&
Hull geometry and loading condition Non-parametric models: frequency response functions Parametric fluid memory model
Figure 1: Hydrodynamic modelling procedure.4
Technical Report
2 Dynamics of Ships and Oﬀshore Strucutres
The linearised equation of motion of marine structure can be formulated as
M
RB
¨
ξ
=
τ
.
(1)The matrix
M
RB
is the rigid-body generalised mass. The generalised-displacement vector
ξ
[
x,y,z,φ,θ,ψ
]
T
gives the position of the body-ﬁxed frame with respect to an equilibrium frame (
x
-surge,
y
-sway, and
z
-heave) and the orientation in terms of Euler angles (
φ
-roll,
θ
-pitch, and
ψ
-yaw).The generalised force vector and
τ
[
X,Y,Z,K,M,N
]
T
gives the respective forces and moments inthe six degrees of freedom. This force vector can be separated into three components:
τ
=
τ
rad
+
τ
visc
+
τ
res
+
τ
exc
,
(2)where the ﬁrst term corresponds to the radiation forces arising from the change in momentum of theﬂuid due to the motion of the structure and the waves generated as the result of this motion, thesecond term corresponds to forces due to ﬂuid viscous eﬀects, the third term corresponds to restoringforces due to gravity and buoyancy, and the fourth component represents the pressure forces due tothe incoming waves other forces used to control the motion of the marine structure.Cummins (1962) used potential theory to study the radiation hydrodynamic problem in thetime-domain for an ideal ﬂuid (no viscous eﬀects) and found the following representation:
τ
rad
=
−
A
∞
¨
ξ
−
t
0
K
(
t
−
t
′
)˙
ξ
(
t
′
)
dt
′
.
(3)The ﬁrst term in (3) represents pressure forces due the accelerations of the structure, and
A
∞
is aconstant positive-deﬁnite matrix called
inﬁnite-frequency added mass
. The second term representsﬂuid-memory eﬀects that capture the energy transfer from the motion of the structure to the radiatedwaves. The convolution term is known as a
ﬂuid-memory model
. The kernel of the convolution term,
K
(
t
), is the matrix of
retardation
or
memory functions
(impulse responses).By combining terms and adding the linearised restoring forces
τ
res
=
−
G
ξ
, the
Cummins Equation
(Cummins, 1962) is obtained:(
M
RB
+
A
∞
)¨
ξ
+
t
0
K
(
t
−
t
′
)˙
ξ
(
t
′
)
dt
′
+
G
ξ
=
τ
exc
,
(4)Equation (4) describes the motion of ships and oﬀshore structures in an ideal ﬂuid provided thelinearity assumption is satisﬁed. This model can then be embellished with non-linear componentstaking into account, for example, viscous eﬀects and mooring lines–see Figure 1.
3 Frequency-domain Models
When the radiation forces (3) are considered in the frequency domain, they can be expressed as follows(Newman, 1977; Faltinsen, 1990):
τ
rad
(
jω
) =
−
A
(
ω
)¨
ξ
(
jω
)
−
B
(
ω
)˙
ξ
(
jω
)
.
(5)5

Search

Similar documents

Related Search

a different reason for the building of SilburA Practical Method for the Analysis of Genetia different reason for the building of SilburA Practical Method for the Analysis of GenetiEuropean Group for the Study of Deviance and Identification of kinetic and transport phenoIdentification of Maintainers and restorers iCentre for the Study of Medicine and the BodyInitiative for Critical Studies of Men and MaIdentification of Biotypes and Secondary Endo

We Need Your Support

Thank you for visiting our website and your interest in our free products and services. We are nonprofit website to share and download documents. To the running of this website, we need your help to support us.

Thanks to everyone for your continued support.

No, Thanks

SAVE OUR EARTH

We need your sign to support Project to invent "SMART AND CONTROLLABLE REFLECTIVE BALLOONS" to cover the Sun and Save Our Earth.

More details...Sign Now!

We are very appreciated for your Prompt Action!

x