InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Robotics » "Kinematics", book edited by Efren Gorrostieta Hurtado, ISBN 978-953-51-3688-0, Print ISBN 978-953-51-3687-3, Published: December 20, 2017 under CC BY 3.0 license. © The Author(s).

Chapter 2

The Inertia Value Transformation in Maritime Applications

By Holger Korte, Sven Stuppe, Jan-Hendrik Wesuls and Tsutomu Takagi
DOI: 10.5772/intechopen.71445

Article top


Relation of wind speed (Beaufort wind scale) and wave height (Douglas sea scale) for fully developed seas (modified to practice).
Figure 1. Relation of wind speed (Beaufort wind scale) and wave height (Douglas sea scale) for fully developed seas (modified to practice).
Structure of a wave spectrum as superposition of different wave components.
Figure 2. Structure of a wave spectrum as superposition of different wave components.
Comparison of Pierson-Moskowitz spectrum (top) and modified Pierson-Moskowitz spectrum (bottom) for several wind speeds.
Figure 3. Comparison of Pierson-Moskowitz spectrum (top) and modified Pierson-Moskowitz spectrum (bottom) for several wind speeds.
Comparison of modified Pierson-Moskowitz spectra and JONSWAP spectra for 350 km, 200 km and 50 km fetch and several wind speeds.
Figure 4. Comparison of modified Pierson-Moskowitz spectra and JONSWAP spectra for 350 km, 200 km and 50 km fetch and several wind speeds.
6DOF motions of a free-floating ship.
Figure 5. 6DOF motions of a free-floating ship.
CAD snapshot of the simulated crew transfer vessel.
Figure 6. CAD snapshot of the simulated crew transfer vessel.
Comparison of the ship position and orientation for head seas. Solid line for body-fixed reference frame. Dashed line for inertial reference frame.
Figure 7. Comparison of the ship position and orientation for head seas. Solid line for body-fixed reference frame. Dashed line for inertial reference frame.
Comparison of the ship position and orientation for beam seas. Solid line for body-fixed reference frame. Dashed line for inertial reference frame.
Figure 8. Comparison of the ship position and orientation for beam seas. Solid line for body-fixed reference frame. Dashed line for inertial reference frame.

The Inertia Value Transformation in Maritime Applications

Holger Korte1, Sven Stuppe1, Jan-Hendrik Wesuls1 and Tsutomu Takagi2
Show details


Due to recent developments in computer technology, computer-aided investigations of structural movements in a maritime environment have become more relevant during the last years. With regard to mechanically coupled multibody systems in fishery and offshore operations, the analysis of such systems is in the focus of research and development. To analyse multibody systems, forces and moments of all included bodies have to be defined within the same reference frame, which requires a transformation algorithm. Showing the correctness of the transformation algorithm, it must be also applicable for six degrees of freedom (6DOF) motions of a free floating single body in seaways. Therefore, the computation of irregular waves is discussed before the traditional motion description of a floating structure by using the Kirchhoff equations. With these basics, an approach to calculate the motion equations of single bodies within the earth-fixed reference frame is presented before the method of the inertia value transformation. To compare the body-fixed and earth-fixed calculation method, a free-floating crew transfer vessel in irregular waves is simulated and the results are discussed. Finally, the inertia value transformation will be proved by the energy conservation principle on the example of a pure rotating rigid body with none digital calculations.

Keywords: inertia value transformation, wave-disturbed ship motions, wave-structure interaction, six degrees of freedom (6DOF), hydromechanics, inertial kinematics, Euclidian room

1. Introduction

Modern simulation techniques enable a more profound analysis of various engineering problems in an early design stage. In maritime kinematics, the focus is on the behaviour of offshore structures under natural environment conditions like wind, waves and current. Detailed knowledge of loads and motions is required for proper dimensioning of efficient and safe systems. Furthermore, simulations are increasingly used for the design of controller-based automation and assistance systems. In addition to the accurate calculation of the structural movements, the calculation speed of a simulation is an important quality feature for these applications, e.g. real-time constraints. To determine the kinematics model of a free-floating structure in a seaway, different approaches are commonly used. The equations of motions can either be defined in the inertial reference frame or alternatively in the body-fixed reference frame. The difference between dry mechanics and maritime mechanics is based on additional hydrodynamic effects, namely, the hydrodynamic added mass, e.g. [13]. The hydrodynamic mass force is an inertia force. The relative acceleration between incompressible fluid and structure induces a pressure field, which results in a hydrodynamic force that is formulated as the product of the relative acceleration and a ‘virtual’ mass. The size of the hydrodynamic added mass depends primarily on the direction of movement and the geometry of the structure [4].

Due to the phenomenon of added or virtual masses, non-scalar and directed inertia values are necessary for describing translational and rotational motions. Probably, this results in an ordinary practice to describe floating structures within a body-fixed view, e. g. by the Kirchhoff motion equations of floating bodies [1, 3].

The investigation of mechanically coupled and multiple rigid body systems is in the focus of the authors and represents a special difficulty in this area [2, 5]. Within literature before millennium, only Paschen’s algorithm was found describing dynamics of real 3D systems for fishery and mine-hunting systems with flexible, non-elastic numerical elements and rigid bodies by partially neglect of added masses [6, 7]. Later, different marine multibody systems were described with algorithms of structural mechanics [8] or as pure planar motion descriptions, e.g. [913]. To describe the mutual interdependencies, the equations of motions of all involved bodies have to be set up in the same coordinate system, preferably the inertial or earth-fixed system. This includes the transformation of all vectorial entities including the inertia matrices (mass matrix, moments of inertia). Therefore, the use of a transformation algorithm is required, which describes the system in the inertial reference frame. The inertia value transformation here presented, also known as Kane’s method [14], was introduced by Korte and Takagi [2, 5] for fishery systems to analyse forces and motions of purse seines. Based on the principle of concentrated masses, all included nodes of the net structure are focused and connected with damped mass-spring elements. To connect the inertia values of the nodes, a transformation of the hydrodynamic added masses is required [15].

Korte et al. presented an application of the inertia value transformation for multibody systems in 2015 for a pure rotating body without elasticity [5]. The analysed scenario represents the transfer of offshore service staff from a crew transfer vessel (CTV) to a wind turbine, which is a real and environmentally affected gyro. CTV is constrained at the bow to the landing of the wind turbine in seas. The inertial reference system is located at the contact point of the structure. The ship can make an ideal rotation around all three axes of the inertial system. The motion equations of the system are defined within the inertial system to calculate the constraining forces. It was observed that temporal integration of the equations within the inertial reference frame leads to unpredictable motions of the ship. The simulation was unstable. During development of the multibody scenario, it was discovered that an additional and oppositely directed transformation of the rotational accelerations leads to a stable solution of the simulation.

To analyse a general applicability of the inertia value transformation for maritime applications, Korte et al. defined the equations of motions for a free-floating ship with six degrees of freedom (6DOF) in the inertial reference frame [17]. Motion equations were derived from momentum and angular momentum theorem. The results of the free-floating ship were compared with ordinary simulations in the body-fixed reference frame and kinematics (see [1]).

This contribution describes the applicability of the inertia value transformation in maritime applications while transferring the inertia values between a body-fixed and the inertial reference frame. The motion behaviour of the free-floating structure is in the focus of the discussion. In the following section, the modelling of the irregular seas as environmental reason of structure motions is presented. The six degrees of freedom of a floating body as well as the traditional motion equations in the body-fixed reference system are introduced afterwards. The derivation of the motion equations within the inertial reference frame by using the inertia value transformation is part of the main section. The possibility to transform the equations from inertial to body-fixed system is shown. For comparison of both methods, 6DOF simulation results of a free-floating CTV in waves are discussed. In the last section, a proof shows the energy conservation of the inertia value transformation for the case of a pure rotating rigid body.

2. Wave modelling

Depending on formation is natural irregular seas divided into wind seas and swell. Real occurred sea state phenomena are described in nautical practice as the superposition of a stochastic wind sea and two observable swells.

Wind seas are direct wind excited. The waves are rather short and steep. Due to friction, the wind transfers energy to the water surface. The resulting capillary waves with short wavelength increase the area of the free surface, whereby the effect of the energy transfer is amplified. This leads in higher and longer gravity waves, which influence the water column down to a water depth of a half wavelength. The waves propagate in the wind direction. Wind wave characteristics depend on wind speed, wind duration and fetch. The wave components of wind seas are located in the higher-frequency range of a wave spectrum. Due to containing higher frequency waves with small amplitudes, the water column oscillates inhomogeneously. Wind seas decrease by a disappearance of the wind excitation. The relation between wind classes (Beaufort wind scale) and sea-state classes (Douglas sea scale) can be seen in Figure 1.


Figure 1.

Relation of wind speed (Beaufort wind scale) and wave height (Douglas sea scale) for fully developed seas (modified to practice).

Contrary to wind seas, which is described as evolving waves, is swell the full homogeneously oscillating water column. Swell is a decaying, nearly sinusoidal wave that has moved away from the formation area, e.g. a storm region. The frequencies are in the lower range of a wave spectrum. The energy density of waves is inversely proportional to the wave frequency. Long waves with a small frequency are more energetic than short waves.

Irregular waves are mathematically described as superposition of a finite number of wave components with different wave frequencies, different amplitudes and different phases. The representation is carried out with energy density spectra, which have a unique distribution of the wave components for each sea area. The spectrum contains the energy Snω of all included wave frequencies. Figure 2 shows the structure of a wave spectrum.


Figure 2.

Structure of a wave spectrum as superposition of different wave components.

Relevant wave parameters can be derived from the wave spectrum. The mean wave frequency ω0 is located at the maximum Sω of the spectrum. The peak period TP is

Further, wave properties can be determined with the so-called spectral moments mn, compared with [1]. Their definition is


The moments up to the second order, n=2, are of interest to characterise the sea state. The characteristic wave height of the spectrum is

and the significant wave height is

The period of the characteristic wave height Hm0 is defined as

The mean zero upcrossing period is

Variance σ2 and standard deviation σ are


There are several mathematical wave spectra for computer simulations of irregular waves available, namely, the Bretschneider spectrum, the Phillips spectrum, the Pierson-Moskowitz spectrum as well as the Joint North Sea Wave Project (JONSWAP) spectrum. The most important spectra are introduced in this section.

2.1. Pierson-Moskowitz spectrum

The Pierson-Moskowitz (PM) spectrum was introduced for fully developed seas and is commonly used since 1964 [18]. It is based on long-term observation data of a weather ship in the period from 1955 to 1960. It assumes the fully developed sea, including wind sea and swell, on the North Atlantic with unlimited water depth and fetch as well as steady wind for a long time. The spectrum is


with parameters A and B


The relation between significant wave height and wind speed is:


2.2. Modified Pierson-Moskowitz spectrum

After experiences with the Pierson-Moskowitz spectrum, the maritime community recommended a magnification of the spectrum [19, 20]. According to this recommendation, the parameters A and B were modified and can be determined from the observable weather parameters HS and TZ:

Figure 3 shows a comparison of the Pierson-Moskowitz spectrum and the Modified Pierson-Moskowitz (MPM) spectrum for several wind speeds. An increase of the energy density in the peak frequency of the MPM as well as a shift of the peak frequency to lower frequencies can be seen.


Figure 3.

Comparison of Pierson-Moskowitz spectrum (top) and modified Pierson-Moskowitz spectrum (bottom) for several wind speeds.

2.3. JONSWAP spectrum

The Pierson-Moskowitz spectra were developed for unlimited water depth. However, the most intensively used sea areas are often in regions of shelf seas with restrictions on the water depth and the fetch.

In 1968/1969, the characteristics of wave formations in sea areas with limited fetch and water depth were investigated exemplary for the North Sea in the international joint project ‘Joint North Sea Wave Project (JONSWAP)’. During a period of approximately 10 weeks, measurements were carried out and evaluated at points reaching up to 160 km seawards in the region of the island of Sylt. The JONSWAP spectrum is based on the PM spectrum and provided with a magnification factor γ for the peak distribution.

Hasselmann’s mathematical definition of the JONSWAP spectrum is [20]




Fetch is defined in metres and increases the wave energy linearly. A direct comparison of the JONSWAP spectrum with the MPM spectrum for a large-scale wind sea on the North Sea illustrates the lower total energy in Figure 4. The characteristic frequencies ω0 of the individual wind classes are higher than those of the PM spectrum. Due to limited water depth, the wave heights are reduced. Both, the maximum of the spectrum and the area under the curve are considerable smaller than those of the reference spectrum. A comparison of seaways with different fetch lengths confirms this trend (cf. Figure 4). A JONSWAP spectra are used for the here presented simulations.


Figure 4.

Comparison of modified Pierson-Moskowitz spectra and JONSWAP spectra for 350 km, 200 km and 50 km fetch and several wind speeds.

3. 6DOF motions of a free-floating offshore structure

In the following section, the six degrees of freedom of a ship are introduced. Figure 5 shows a sketch of a ship hull with a body-fixed coordinate system. Body-fixed coordinate system means that the coordinate system moves with the ship. Due to advantages in determining the moments of inertia, the origin of the system is preferably located in the ships centre of gravity. The motions of a ship are described in the body-fixed coordinate system. However, the ship’s position and the orientation are referenced in an earth-fixed coordinate system.


Figure 5.

6DOF motions of a free-floating ship.

The translational motions are as follows:

  • Surge: Translation along the longitudinal x-axis

  • Sway: Translation along the transversal y-axis

  • Heave: Translation along the vertical z-axis

The rotational motions are as follows:

  • Roll: Rotation around the longitudinal x-axis

  • Pitch: Rotation around the transversal y-axis

  • Yaw: Rotation around the vertical z-axis

Another distinction is the classification into horizontal (surge, sway, yaw) and vertical (heave, roll, pitch) motions. The reason for this differentiation is the restoring forces caused by gravitation. These exist only for the vertical degrees of freedom. After a perturbation of the equilibrium, a ship always tends to return to it. This does not apply for horizontal motions. The result is the drift of a free-floating structure away from the original position as well as a change of the heading in case of wind, waves or current. For a driven ship or a ship in dynamic positioning mode, the ship’s actuators (propeller, rudder, thruster, etc.) control the horizontal motions.

4. Motion equations in the body-fixed reference system: Kirchhoff equations

As already mentioned, hydrodynamic inertial effects in the form of a hydrodynamic added mass have to be taken into account in marine applications. The mass matrix mb (b for body) in the traditionally used body-fixed reference system is defined in Eq. (22). In the matrix, m is the physical mass and mh,ij is the direction-dependent content of the hydrodynamic inertia. The mass matrix is constant in the body-fixed frame. This applies analogously for the matrix of moments of inertia:


In 1869, Kirchhoff published his work about the ‘Movement of a rotating Body in a Fluid’ [3]. In this work, he defined the motion equations of a floating body in a body-fixed reference system in analogy to Euler’s gyro equation. The Kirchhoff equations are a system of each three equations for translation (Eq. (23)) and rotation (Eq. (24)):


The so-called living forces and moments depict all external forces and moments including all hydrodynamic effects, weight, buoyancy and their effect and can be found at the right-hand side of the equations inclusively. The determination of these external forces and moments is especially difficult for the horizontal degrees of freedom with a lack of restoring forces. For the simulations in this work, a simplified model is used, which considers weight, buoyancy as well as potential damping. At present, viscous effects are neglected.

The accelerations have to be integrated twice to analyse the position and orientation of the free-floating ship. The first integration is executed in the body-fixed system.

To calculate the position x,y,z and orientation Φ,Θ,Ψ of a floating body, the velocities have to be transformed into the inertia system. The peculiarity is the differing transformation of the rotational vector from the body-fixed system into the inertia system and vice versa, e.g. using the transformation with Euler’s angles (Eq. (26), cf. [1]):




(e for earth or inertial, b for body, c for cosine, s for sine and t for tangent).

5. Motion equations in the inertial reference system

In the following section, the motion equations of a free-floating body in the earth-fixed reference system are shown. The equations are formally derived from momentum and angular momentum theorem. The transformation of the inertia values is discussed, and it is shown that a transformation of the motion equations leads to the Kirchhoff equations. It is known that the twice integration of the rotational equations results in an unstable solution. By using an additional and opposite directed transformation of the rotational accelerations, the simulations can be stabilised for longer time series.

5.1. Inertia value transformation

Kane introduced the transformation of inertia values in 1985 [14]. The used transformation matrix from the body-fixed into the earth-fixed reference system is





The transformation of the inertia matrix is [14, 21]


Analogously, it follows the transformation of the mass matrix including added mass [2]:


The cross product-operator relation for matrices can be found in the literature [22, 23]:


5.2. Rotation

The rotational equations follow from angular momentum theorem, which is defined in the inertial reference frame [21, 23]:


Taking the temporal change of the inertia matrix J¯e into account follows for the moment:


The temporal derivative of the inertia matrix is


Substituting Eq. (35) into Eq. (34) follows:


Using the cross product-operator relation follows the momentum equation in the inertial reference system:


The transformation of Eq. (37) into the body-fixed reference frame is


It follows with


Eq. (39) is identical to Euler’s gyroscope equation, which is defined in the body-fixed system and corresponds to the Kirchhoff rotation equation (Eq. (24)) if the term v×P is neglected.

As already mentioned, the integration of the equations in the inertial system results in unstable behaviour. The floating body makes unpredictable, chaotic movements. It is assumed that the instability is a consequence of the numerical inaccuracy as well as the twice integration of the equations. It was found within the project Mine Hunting 2000 that an additional and opposing transformation stabilises the simulation for longer time [6]. Due to the missing inertia value algorithm at that time, it could not be explained why. The authors assumed a numerical reason, which cannot be proved. That phenomenon of stabilisation was rebuilt for simulations of docked CTV at wind turbine tower, compare [24] with [16].

During the presented motion of the free-floating vessel, the changed rotational accelerations stabilise the system for longer periods (cf. Eq. (40)):


By using the additional transformation in the multibody application described in Ref. [16], the motions of the vessel on the wind turbine tower were stable at all times.

5.3. Translation

Consequently, the derivation of the translational equations follows the derivation of the rotational equations from the last section. The hydrodynamic added mass force is generally an external force, which defined the right-hand side of the motion equations. For the derivation of the motion equations, the hydrodynamic mass is considered as intrinsic property of the free-floating body, which has to be taken into account for every accelerated marine system. Analogously, as for the rotational inertia matrix, the mass matrix of a rotating body changes. This temporal change is equal to zero in case of pure translation.

The momentum theorem is defined in the inertial reference system:


The force is the temporal derivative of the momentum:


With the transformation matrix C¯be, the transformation of the mass matrix into the inertial system follows:


The other direction is


The temporal derivative of the mass matrix is


By substituting Eq. (43) into Eq. (45) follows


and by using the cross product-operator relation


Substituting Eq. (47) into Eq. (42) results in the translational equation for a floating body in the inertial reference frame:


Eq. (48) contains only values of the inertial reference system. Due to change of the orientation and the non-scalar characteristic of the mass matrix, Eq. (48) contains two additional terms. These terms disappear for each application where the hydrodynamic added mass is not considered and the mass is scalar, e.g. in aerospace industries or robotics. In case of pure translation, these terms are equal to zero too, and the equation results in Newton’s second law.

In the following equation, the transformation of Eq. (48) into the body-fixed reference system is shown. The relation of the velocity in both systems is


In case of simultaneous rotation and translation, the rotation matrix changes. The acceleration in inertial frame is


(Remark: In case of pure translation, v̇e=Cbev̇b).

By inserting identity matrix, Eq. (28) follows:


This is identical to


Substituting Eq. (52) into Eq. (48), the term m¯eω~eve disappears, and Eq. (48) becomes


The transformation of Eq. (53) into the body-fixed reference frame is


Replacing the values from inertial system with values from body-fixed system follows:


Eq. (55) corresponds to the translational Kirchhoff equation (Eq. (23)) in the body-fixed system.

6. Simulations in body-fixed and inertial reference frame

The following section shows simulation comparisons of a free-floating vessel in both reference systems. 6DOF motion simulations are performed for various wave conditions. The hydrodynamic model of the simulation is implemented in MATLAB/Simulink.

6.1. Parameterisation of the ship

The simulated ship is a crew transfer vessel, which is used for the transfer of offshore service staff in the German Bight. Figure 6 shows a snapshot of the ship’s CAD model. It is a catamaran hull with a length of LoA=22.0m, a breadth of B=8.3m and a mass of m=60t. For the simulation, the hydrodynamic parameters like hydrodynamic added masses and moments of inertia, as well as the potential damping coefficients of the CTV, are required. They were determined within the project ‘Safe Offshore Operations (SOOP)’ using the potential radiation and diffraction programme WAMIT (cf. [24]). The ship is discretised station-wise to calculate wave-induced forces and moments. The calculated forces are buoyancy and weight as well as potential damping force.


Figure 6.

CAD snapshot of the simulated crew transfer vessel.

6.2. 6DOF simulations

The following figures show the comparison of the CTV 6DOF body motions in irregular seas. The ship is free floating with no initial velocity. Figure 7 shows the simulation results for the case of head seas and Figure 8 for beam seas. The wind wave parameters are a significant wave height of HS=2.1m and a peak period of ω0=0.7rad/s. Swell is neglected.


Figure 7.

Comparison of the ship position and orientation for head seas. Solid line for body-fixed reference frame. Dashed line for inertial reference frame.


Figure 8.

Comparison of the ship position and orientation for beam seas. Solid line for body-fixed reference frame. Dashed line for inertial reference frame.

The figures show from top to bottom the registered wave elevation at the ships centre of gravity, the x-, y- and z-position as well as the orientation angles Φ, Θ and Ψ. Blue curves are for the body-fixed reference frame and red curves for simulations in the inertial reference frame.

At the beginning of the presented simulations, the motion behaviour in both reference systems shows an identical behaviour. In further process, however, an increasing deviation of the ship’s heading in the inertial reference is registered. The result is a changing encounter angle of ship and waves. For the case of head seas, the ship begins to move in transverse direction and to roll. In the beam seas, simulation leads the changed encounter angle in a pitch motion. In total, it can be seen that the comparison of vertical motions shows an identical behaviour. As already mentioned, the yaw motion is controlled in case of a driven ship. It is assumed that the motions in both reference systems are identical in this case. Analysis of yaw motion is the subject of the present work. It has to be pointed out that the simulations in the inertial reference system were unstable after longer simulation time. Instability occurs abruptly. Intensive analysis of the instability reason is also the subject of the present work.

7. Proof of energy conservation for rotatory inertia value transformation

Problem: Common practice in calculating rotational motions of rigid bodies applies the so-called angular velocity transformation between a first integration over time in body-fixed frame and a second within inertial frame (see [1]). This is contrary to the so-called linear velocity transformation, which uses the ordinary vector transformation algorithm as applied in our solution [2], e.g. Euler’s angle transformation. The derivation of the rotation matrix C¯rot,be (Eq. (26)) is carried out by observations and describes the rotation around the three body-fixed coordinate axes sequentially in time.

In 1995, Blass and Gurevich postulates ‘Matrix Transformation Is Complete for the Average Case’ [25]:


whereby C¯be describes the ordinary transformation matrix for vector values from coordinate system b, here body fixed, to the inertial or earth-fixed ones e. Which transformation is correct, way I in comparison to Eq. (25) or way II regarding Eqs. (30) and (49)?

I:            if ωe=C¯rot,beωb                                         is true; also,

               Me=C¯rot,beJ¯bC¯rot,be1C¯rot,beωb       might be possible being true, but wrong in               our opinion, because C¯rot,be1 for each Θ0

or II:        if ωe=C¯beωb                                               is true, contrary to the above

                Me=C¯beJ¯bC¯beC¯beωb                          have to be true!

Matrix mathematics cannot give answer, because both assumptions are valid.

Authors expect calculus II as valid due to the correct transformation of the moment vector.

Its correctness may prove only by the general use of one physical principal of conservation. For that purpose, the energy conservation shall be applied. Therefore, the pure rotating energy from one mass point of the rigid body (see step 1) has to be compared with the kinetic energy of same mass point and motion applying continuum physics within FEM methods (see step 2). Furthermore and in accordance with authors who claim of wrong description of rotatory kinetic energy within the body-fixed frame, all used rotation speed components are time derivatives of inertial-fixed Euler’s angles. Within the calculation of step 3, it has to be shown that all coefficients of both energy equations are equal, otherwise the proof fails.

7.1. Calculation of rotatory kinetic energy of a mass point from a rotating body (step 1)

Vector value transformation based on Eulerian angles Πe=ΦΘΨT using the well known transformation matrix, compare with [1] (Eq. (27)). The components Cij of rotation matrix are


The rotatory inertia value J¯dm of the mass point calculates


by angular speed


and lever arm Rm=RxRyRzT from the bodies’ centre of gravity to the mass point in body-fixed frame b. By vectorial depiction, the rotatory kinetic energy can be formulated as follows, e.g. Ginsberg [21]:


The first calculation step can be C¯ebωe (cf. Eq. (61)):


Now, J¯dmC¯ebωe can be estimated:


In the following equation, the transposed vector from Eq. (61) is multiplied with the vector of Eq. (62):


getting three components of a sum E1 to E3:


Finally, the sum of these components E1, E2 and E3 has to be multiplied by the factor ½dmto get the rotatory kinetic energy. Because that factor is implemented also within the second energy equation (step 2), the energy components can be summarised and compared directly with Eq. (68). By sorting all components to the lever arm and rotating speed combinations, we get Eq. (64):


7.2. Calculation of kinetic energy from a mass point of a pure rotating body (step 2)

The energy can also be calculated by the tangential speed vector, preferably within the inertial frame (Eq. (65)):


The cross product can be calculated by a so-called skew operator ω~e (Eq. (32)). So, it is possible to calculate the tangential speed vector of the mass point within the inertial frame:


From Eq. (66) the square of speed can be calculated (see Eq. (67)):




Summarising E4, E5 and E6, follows Eq. (68).


7.3. Comparison of the coefficients from both energy calculations for the mass point of pure rotating body (step 3)

To show that energy is identical in both approaches (step 1 and step 2), the coefficients of all 36 combinations of lever arms and angular velocities are compared individually with each other in step 3. The coefficients of the rotation energy (cf. Eq. (64)) of step 1 can be found on the left-hand side and the coefficients of kinetic energy (cf. Eq. (68)) of step2 on the right-hand side. As described previously, the condition of identical energy is only proved, if all coefficients ar equal.


The previous calculation demonstrated the validity of Eq. (35).

8. Conclusion

This work presents the method of inertial value transformation for maritime applications. Firstly, an introduction in irregular seas and an overview of ships degrees of freedom are given. In the following section, the traditional Kirchhoff motion equations in the body-fixed reference frame are introduced, which represent a hydrodynamic affected Eulerian gyro tied up a Newtonian body within the body-fixed view. The formal derivation of motion equations of a free-floating body in inertial coordinate system is presented in the main part. It is shown that the transformation of the equations into the body-fixed system leads to the well-known Kirchhoff motion equations. 6DOF simulations for a crew transfer vessel in head seas and beam seas illustrate the comparison of motions in both reference systems. It is mentioned that temporal integration of the motion equations in inertial system leads to unstable and chaotic motions of the ship. Rebuilding the Eulerian gyro at first within the project SOOP [6], Korte et al. [24] by introducing an additional and opposite directed transformation of the rotational accelerations, the present work shows their general applicability and necessity for free-moving bodies (6DOF) within inertial frame. By its consequent use, the motion behaviour of the ship can be stabilised over longer periods. A proof presents the energy conservation of inertia value transformation for a rotating body. Finally, the contribution has shown a failure in common motion calculation practice for vehicles.

The intention of motion equations in inertial reference system is the simulation of mechanically coupled multibody systems in seas. To analyse the interaction effects, the forces and moments of all included bodies have to be defined in the same reference system. This requires a transformation of the motion equations that can be realised with the presented method of inertia value transformation. For a multibody system, a CTV is fixed at the bow to an offshore wind turbine tower and can make ideal rotations; the method including the additional transformation leads to a stable system. The scenario describes the interaction of a fixed and a floating body. Further investigations in the field of multibody dynamics are planned for the future. A scenario of two ships, which are mechanically coupled in tandem and rotate ideally, is developed currently. The challenge in comparison to the first multibody system is the interaction of two floating bodies. Other applications of the method, which are planned in further work of the authors, are simulations of offshore crane processes or 3D simulations of a ROV, which is coupled to a mother ship via umbilical. For the parameterisation of controllers, the question of the real-time application is still in focus of research.


The authors wish to thank the Japanese Ministry of Education, Culture, Sports, Science, and Technology for supporting the 21st Century COE Program Center of Aquaculture Science and Technology for Bluefin Tuna and Other Cultivated Fish, the Lower Saxony Ministry for Science and Culture for supporting research programme Safe Offshore Operations (SOOP) and the German Federal Ministry for Economic Affairs and Energy for supporting the ZIM project Information System for Near Real-Time Logistics (IeK).


1 - Fossen TI. Handbook of Marine Craft Hydrodynamics and Motion Control. 1st ed. Chichester UK: J. Wiley & Sons; 1994. DOI: ISBN 0 471 94113 1
2 - Korte H, Takagi T. Dynamic motion calculation of a flexible structure using inertia transformation algorithm. In: Workshop on Fishing and Marine Production Technology; 9-10. August; Trondheim. NTNU. 2004. pp. 1-11. DOI: 10.13140/RG.2.1.4006.4724
3 - Kirchhoff G. Über die Bewegung eines Rotationskörpers in einer Flüssigkeit. Crelles Journal. 1869;71:237-273
4 - Wendel K. Hydrodynamische Massen und hydrodynamische Massenträgheitsmomente. Jahrbuch der Schiffbautechnischen Gesellschaft. 1950;207-255
5 - Korte H, Takagi T. Transformation von Trägheitsgrößen zur Lösung der Mehrkörperdynamik in der Meerestechnik. In: 101. Jahrbuch der Schiffbautechnischen Gesellschaft. Hamburg: Schiffahrts-Verlag “Hansa”; 2007. pp. 313-325
6 - Korte H. Track control of a towed underwater sensor carrier. In: IFAC Control in Transportation Systems. Braunschweig: IFAC; 2000. pp. 89-94
7 - Paschen M. Contribution to Predict Trajectories of Pelagic Fishing Gears Affected by Ship Manoeuvres (in German) [Dissertation]. Rostock: University of Rostock; 1982
8 - Suzuki K, Takagi T, Shimizu T, Hiraishi T, Yamamoto K, Nashimoto K. Validity and visualization of a numerical model used to determine dynamic configuration of fishing net. Fisheries Science. 2003;69(4):1-18. DOI: 10.1046/j.1444-2906.2003.00676.x · Source: OAI
9 - Aamo O, Fossen TI. Finite element modelling of moored vessels. Mathematical and Computer Modelling of Dynamical Systems. 2001;7(1):47-75
10 - Abreu P, Morishita H, Pascoal A, Ribero J, Silva H. Marine vehicles with streamers for geotechnical surveys: Modeling, positioning and control. In: 10th IFAC CAMS; September; Trondheim. IFAC-PapersOnLine 49-23; 2016. p. 458-464
11 - Berntsen P, Aamo O, Leira B, Sørensen A. Structural reliability - based control of moored interconnected structures. Control Engineering Practice. 2008;16(4):495-504
12 - Korte H. Motion modelling of a manoeuvring ship at anchor in flowing waters. In: Yamane T, editor. Contributions on the Theory of Fishing Gears and Related Marine Systems; 5.-7. November; University of Kinki, Nara (Jp). Nara:2009. pp. 207-216. DOI: ISBN 4-946-421-13-6
13 - Yang K, Wang X, Wu C. Simulation platform for the underwater snakelike robot swimming based on Kane's dynamic model and central pattern generator. Journal Shanghai Jiaotong University. 2014;19(3):294-301
14 - Kane TR, Levison DA. Dynamics: Theory and Applications. 1st ed. New York: McGraw-Hill; 1985
15 - Takagi T, Miyata S, Fusejima I, Oshima T, Uehara T, Suzuki K, Nomura Y, Kanechiku M, Torisawa S. Potential of computer simulation for buoy-line type of purse seine fishing. In: Paschen M, editors. Contributions on the Theory of Fishing Gears and Related Marine Systems; 9-12. Ocktober; Rostock. Aachen: Shaker Verlag; 2013. p. 55-62. DOI: ISBN 978-3-8440-2251-3
16 - Korte H, Stuppe S, Wesuls J-H, Takagi T. Application of the inertia value transformation to simulate the transfer scenario of service staff from a crew transfer vessel to offshore wind turbines. In: Paschen M, O’Neil FG, Ivanovic A, editors. Contributions on the theory of fishing gears and related marine systems; 27-29. October; Aberdeen. Aachen: Shaker Verlag; 2015. p. 227-237. DOI: ISBN 978-3-8440-3955-9
17 - Korte H, Wesuls J-H, Stuppe S, Takagi T. Towards an approach to 6 DOF inertial kinematics. In: 10th IFAC CAMS; September; Trondheim. IFAC-PapersOnLine 49; 2016. p. 440-445
18 - Piersson WJ, Moskowitz L. A proposed spectral form for fully developed wind seas based on the similarity of S.A. Kitagorodskii. Journal of Geophysical Research. 1964;96(2)
19 - ITTC. Report of the Committee on Waves. Final Report and Recommendations to the 23rd International Towing Tank Conference. 2002;505-551
20 - Hasselmann K, Barnet TP, Bouws E, Carlson H, Cartwright DE, Enke K, Ewing JA, Gienapp H, Hasselmann DE, Krusemann P, Meerburg A, Müller P, Olbers DJ, Richter K, Seil W, Walden H. Measurements of wind-wave growth and swell decay during the joint North Sea wave project (JONSWAP). Ergänzungsheft zur Deutschen Hydrographischen Zeitschrift. 1973;Reihe A(8)(12):95
21 - Ginsberg JH. Advanced Engineering Dynamics. 2nd ed. New York: Cambridge University Press; 2004
22 - Woernle C. Dynamic of Multi-Body Systems. Rostock: University of Rostock, Lesson script; 2001
23 - Schiehlen W, Eberhard P. Technische Dynamik–Rechnergestützte Modellierung mechanischer Systeme im Maschinen- und Fahrzeugbau. 4th ed. Springer Vieweg; 2014. pp. 24-28. DOI: ISBN 978-3-658-06184-5
24 - Korte H, Ihmels I, Richter J, Zerhusen B, Hahn A. Offshore training simulations. In: 9th IFAC Manoeuvring and Control of Marine Crafts; 19–21 September 2012; Arenzano (IT). IFAC; 2012. pp. 37-42. DOI: ISSN 1474-6670
25 - Blass A, Gurevich Y. Matrix transformation is complete for the average case. SIAM Journal on Computing. 1995;24(1):3-29