Open access peer-reviewed chapter

# The Inertia Value Transformation in Maritime Applications

Written By

Holger Korte, Sven Stuppe, Jan-Hendrik Wesuls and Tsutomu Takagi

Submitted: April 21st, 2017 Reviewed: October 4th, 2017 Published: December 20th, 2017

DOI: 10.5772/intechopen.71445

From the Edited Volume

## Kinematics

Chapter metrics overview

View Full Metrics

## Abstract

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.

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.

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

TP=2πω0E1

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

mn=0ωnSωE2

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

Hm0=4m0E3

and the significant wave height is

HS=3.81m0E4

The period of the characteristic wave height Hm0is defined as

T1=m0m1E5

The mean zero upcrossing period is

TZ=m0m2E6

Variance σ2and standard deviation σare

σ2=m0E7
σ=m0=0.263HSE8

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

SPMω=Aω5eBω4E9

with parameters Aand B

A=0.0081g2E10
B=0.74gVwind4=2.814HS2E11

The relation between significant wave height and wind speed is:

HS=1.95g2Vwind2E12

### 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 Aand Bwere modified and can be determined from the observable weather parameters HSand TZ:

A=4π3HS2TZ4E13
B=16π3TZ4E14

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.

### 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]

SJONω=αg2ω5e54ω0ω4γrE15

with

r=eωω022σ2ω02E16
γ=3.3E17
α=0.076V102FWg0.22E18
ω0=22g2FWV1013E19
σ=0.07,ωω00.09,ω>ω0E20

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 ω0of 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.

## 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.

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(bfor body) in the traditionally used body-fixed reference system is defined in Eq. (22). In the matrix, mis the physical mass and mh,ijis 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:

mb,ij=m+mh,ijE21
m¯b=m+mh,11000m+mh,22000m+mh,33E22

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)):

dPdt+ω×P=X,Y,ZE23
dLdt+ω×L+v×P=K,M,NE24

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,zand 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=C¯rot,beωbE25

with

C¯rot,be=1sΦtΘcΦtΘ00=C¯rot,eb1E26

(efor earth or inertial, bfor body, cfor cosine, sfor sine and tfor 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

C¯be=cΘcΨsΦsΘcΨcΦsΨcΦsΘcΨ+sΦsΨcΘsΨsΦsΘsΨ+cΦcΨcΦsΘsΨsΦcΨsΦcΘcΦcΘE27

with

C¯beC¯be=100010001=C¯beC¯beE28

and

C¯eb=C¯beE29

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

J¯e=C¯beJ¯bC¯beE30

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

m¯e=C¯bem¯bC¯beE31

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

C¯̇beC¯be=ω~e=0Ψ̇Θ̇Ψ̇0Φ̇Θ̇Φ̇0=ωe×E32

### 5.2. Rotation

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

Le=J¯eωeE33

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

dLedt=Me=J¯eω̇e+J¯̇eωeE34

The temporal derivative of the inertia matrix is

J¯̇e=C¯̇beJ¯bC¯be+C¯beJ¯bC¯̇beE35

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

Me=J¯eω̇e+C¯̇beJ¯bC¯be+C¯beJ¯bC¯̇beωeE36

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

Me=J¯eω̇e+ω~eJ¯eωe=J¯eω̇e+ωe×J¯eωeE37

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

Mb=C¯ebMe=C¯ebJ¯eC¯ebC¯ebω̇e+C¯ebω~eC¯ebC¯ebJ¯eC¯ebC¯ebωeE38

It follows with

Mb=J¯bω̇b+ω~bJ¯bωb=J¯bω̇b+ωb×J¯bωbE39

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×Pis 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)):

ωe=C¯ebω̇bdtC¯ebω̇edtE40

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:

Pe=m¯eveE41

The force is the temporal derivative of the momentum:

dPedt=Fe=m¯ev̇e+m¯̇eveE42

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

m¯e=C¯bem¯bC¯beE43

The other direction is

m¯b=C¯bem¯eC¯beE44

The temporal derivative of the mass matrix is

m¯̇e=C¯̇bem¯bC¯be+C¯bem¯bC¯̇beE45

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

m¯̇e=C¯̇beC¯bem¯e+m¯eC¯beC¯̇beE46

and by using the cross product-operator relation

m¯̇e=ω~em¯e+m¯eω~eE47

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

Fe=m¯ev̇e+ω~em¯evem¯eω~eve=m¯ev̇e+ωe×m¯evem¯eωe×veE48

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

ve=C¯bevbE49

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

v̇e=C¯̇bevb+C¯bev̇bE50

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

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

v̇e=C¯̇beC¯beC¯bevb+C¯bev̇bE51

This is identical to

v̇e=ω~eve+C¯bev̇b=ωe×ve+C¯bev̇bE52

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

Fe=m¯eC¯bev̇b+ω~em¯eveE53

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

Fb=C¯ebFe=C¯ebm¯eC¯ebC¯ebC¯bev̇b+C¯ebω~eC¯ebC¯ebm¯eC¯ebC¯ebveE54

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

Fb=m¯bv̇b+ω~bm¯bvb=m¯bv̇b+ωb×m¯bvbE55

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.3mand 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.

### 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.1mand a peak period of ω0=0.7rad/s. Swell is neglected.

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 transformationbetween 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]:

A¯e=C¯beA¯bC¯be1E56

whereby C¯bedescribes 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,be1for 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=ΦΘΨTusing the well known transformation matrix, compare with [1] (Eq. (27)). The components Cijof rotation matrix are

C¯be=C11C12C13C21C22C23C31C32C33=cΘcΨsΦsΘcΨcΦsΨcΦsΘcΨ+sΦsΨcΘsΨsΦsΘsΨ+cΦcΨcΦsΘsΨsΦcΨsΦcΘcΦcΘE57

The rotatory inertia value J¯dmof the mass point calculates

J¯dmdm=dmRy2+Rz2RxRyRxRzRxRyRx2+Rz2RyRzRxRzRyRzRx2+Ry2E58

by angular speed

Π̇e=ωe=Φ̇Θ̇Ψ̇TE59

and lever arm Rm=RxRyRzTfrom 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]:

Erot=12dmωeTC¯beJ¯dmC¯ebωeE60

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

C¯ebωe=C11Φ̇+C21Θ̇+C31Ψ̇C12Φ̇+C22Θ̇+C32Ψ̇C13Φ̇+C23Θ̇+C33Ψ̇E61

Now, J¯dmC¯ebωecan be estimated:

Ry2+Rz2RxRyRxRzRxRyRx2+Rz2RyRzRxRzRyRzRx2+Ry2C11Φ̇+C21Θ̇+C31Ψ̇C12Φ̇+C22Θ̇+C32Ψ̇C13Φ̇+C23Θ̇+C33Ψ̇=C11Ry2Φ̇+C11Rz2Φ̇+C21Ry2Θ̇+C21Rz2Θ̇+C31Ry2Ψ̇+C31Rz2Ψ̇C12RxRyΦ̇C22RxRyΘ̇C32RxRyΨ̇C13RxRzΦ̇C23RxRzΘ̇C33RxRzΨ̇C11RxRyΦ̇C21RxRyΘ̇C31RxRyΨ̇+C12Rx2Φ̇+C12Rz2Φ̇+C22Rx2Θ̇+C22Rz2Θ̇+C32Rx2Ψ̇+C32Rz2Ψ̇C13RyRzΦ̇C23RyRzΘ̇C33RyRzΨ̇C11RxRzΦ̇C21RxRzΘ̇C31RxRzΨ̇C12RyRzΦ̇C22RyRzΘ̇C32RyRzΨ̇C13Rx2Φ̇+C13Ry2Φ̇+C23Rx2Θ̇+C23Ry2Θ̇+C33Rx2Ψ̇+C33Ry2Ψ̇E62

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

C¯ebωeTJ¯dmC¯ebωe=E1+E2+E3E63

getting three components of a sum E1to E3:

E1=C112Ry2Φ̇2+C112Rz2Φ̇2+C11C21Ry2Φ̇Θ̇+C11C21Rz2Φ̇Θ̇+C11C31Ry2Φ̇Ψ̇+C11C31Rz2Φ̇Ψ̇C11C12RxRyΦ̇2C11C22RxRyΦ̇Θ̇C11C32RxRyΦ̇Ψ̇C11C13RxRzΦ̇2C11C23RxRzΦ̇Θ̇C11C33RxRzΦ̇Ψ̇+C11C21Ry2Φ̇Θ̇+C212Rz2Θ̇2+C11C21Rz2Φ̇Θ̇+C212Ry2Θ̇2+C21C31Ry2Θ̇Ψ̇+C21C31Rz2Θ̇Ψ̇C12C21RxRyΦ̇Θ̇C21C22RxRyΘ̇2C21C32RxRyΘ̇Ψ̇C13C21RxRzΦ̇Θ̇C21C23RxRzΘ̇2C21C33RxRzΘ̇Ψ̇+C11C31Ry2Φ̇Ψ̇+C11C31Rz2Φ̇Ψ̇+C21C31Ry2Θ̇Ψ̇+C21C31Rz2Θ̇Ψ̇+C312Ry2Ψ̇2+C312Rz2Ψ̇2C12C31RxRyΦ̇Ψ̇C22C31RxRyΘ̇Ψ̇C32C31RxRyΨ̇2C13C31RxRzΦ̇Ψ̇C23C31RxRzΘ̇Ψ̇C31C33RxRzΨ̇2
E2=C11C12RxRyΦ̇2C12C21RxRyΦ̇Θ̇C12C31RxRyΦ̇Ψ̇+C122Rx2Φ̇2+C122Rz2Φ̇2+C12C22Rx2Φ̇Θ̇+C12C22Rz2Φ̇Θ̇+C12C32Rx2Φ̇Ψ̇+C12C32Rz2Φ̇Ψ̇C12C13RyRzΦ̇2C12C23RyRzΦ̇Θ̇C12C33RyRzΦ̇Ψ̇C11C22RxRyΦ̇Θ̇C21C22RxRyΘ̇2C22C31RxRyΘ̇Ψ̇+C12C22Rx2Φ̇Θ̇+C12C22Rz2Φ̇Θ̇+C222Rx2Θ̇2+C222Rz2Θ̇2+C22C32Rx2Θ̇Ψ̇+C22C32Rz2Θ̇Ψ̇C13C22RyRzΦ̇Θ̇C22C23RyRzΘ̇2C22C33RyRzΘ̇Ψ̇C11C32RxRyΦ̇Ψ̇C21C32RxRyΘ̇Ψ̇C31C32RxRyΨ̇2+C12C32Rx2Φ̇Ψ̇+C12C32Rz2Φ̇Ψ̇+C22C32Rx2Θ̇Ψ̇+C22C32Rz2Θ̇Ψ̇+C322Rx2Ψ̇2+C322Rz2Ψ̇2C13C32RyRzΦ̇Ψ̇C23C32RyRzΘ̇Ψ̇C32C33RyRzΨ̇2
E3=C11C13RxRzΦ̇2C13C21RxRzΦ̇Θ̇C13C31RxRzΦ̇Ψ̇C12C13RyRzΦ̇2C13C22RyRzΦ̇Θ̇C13C32RyRzΦ̇Ψ̇+C132Rx2Φ̇2+C132Ry2Φ̇2+C13C23Rx2Φ̇Θ̇+C13C23Ry2Φ̇Θ̇+C13C33Rx2Φ̇Ψ̇+C13C33Ry2Φ̇Ψ̇C11C23RxRzΦ̇Θ̇C21C23RxRzΘ̇2C23C31RxRzΘ̇Ψ̇C12C23RyRzΦ̇Θ̇C22C23RyRzΘ̇2C23C32RyRzΘ̇Ψ̇+C13C23Rx2Φ̇Θ̇+C13C23Ry2Φ̇Θ̇+C232Rx2Θ̇2+C232Ry2Θ̇2+C23C33Rx2Θ̇Ψ̇+C23C33Ry2Θ̇Ψ̇C11C33RxRzΦ̇Ψ̇C21C33RxRzΘ̇Ψ̇C31C33RxRzΨ̇2C12C33RyRzΦ̇Ψ̇C22C33RyRzΘ̇Ψ̇C32C33RyRzΨ̇2+C13C33Rx2Φ̇Ψ̇+C13C33Ry2Φ̇Ψ̇+C23C33Rx2Θ̇Ψ̇+C23C33Ry2Θ̇Ψ̇+C332Rx2Ψ̇2+C332Ry2Ψ̇2

Finally, the sum of these components E1, E2and E3has 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):

2Erotdm=C322+C332Rx2Ψ̇22C31C32RxRyΨ̇22C31C33RxRzΨ̇22C32C33RyRzΨ̇2+C312+C332Ry2Ψ̇2+C312+C322Rz2Ψ̇2+2C22C32+C23C33Rx2Θ̇Ψ̇2C22C31+C21C32RxRyΘ̇Ψ̇2C21C33+C23C31RxRzΘ̇Ψ̇2C22C33+C23C32RyRzΘ̇Ψ̇+2C23C33+C21C31Ry2Θ̇Ψ̇+2C22C32+C21C31Rz2Θ̇Ψ̇+2C13C33+C12C32Rx2Φ̇Ψ̇2C12C31+C11C32RxRyΦ̇Ψ̇2C11C33+C13C31RxRzΦ̇Ψ̇2C12C33+C13C32RyRzΦ̇Ψ̇+2C13C33+C11C31Ry2Φ̇Ψ̇+2C12C32+C11C31Rz2Φ̇Ψ̇+C222+C232Rx2Θ̇22C21C22RxRyΘ̇22C21C23RxRzΘ̇22C22C23RyRzΘ̇2+C212+C232Ry2Θ̇2+C212+C222Rz2Θ̇2+2C13C23+C12C22Rx2Φ̇Θ̇2C11C22+C12C21RxRyΦ̇Θ̇2C11C23+C13C21RxRzΦ̇Θ̇2C12C23+C13C22RyRzΦ̇Θ̇+2C11C21+C13C23Ry2Φ̇Θ̇+2C11C21+C12C22Rz2Φ̇Θ̇+C122+C132Rx2Φ̇22C11C12RxRyΦ̇2.2C11C13RxRzΦ̇22C12C13RyRzΦ̇2+C112+C132Ry2Φ̇2+C112+C122Rz2Φ̇2E64

### 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)):

Ekin=12dmωe×C¯beRmTωe×C¯beRmE65

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:

vdm=ω˜eC¯beRm=C21RxΨ̇C22RyΨ̇C23RzΨ̇+C31RxΘ̇+C32RyΘ̇+C33RzΘ̇C11RxΨ̇+C12RyΨ̇+C13RzΨ̇C31RxΦ̇C32RyΦ̇C33RzΦ̇C11RxΘ̇C12RyΘ̇C13RzΘ̇+C21RxΦ̇+C22RyΦ̇+C23RzΦ̇E66

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

2Ekindm=vdm2=vdmTvdm=E4+E5+E6E67

with.

E4=C212Rx2Ψ̇2+2C21C22RxRyΨ̇2+2C21C23RxRzΨ̇22C21C31Rx2Θ̇Ψ̇2C21C32RxRyΘ̇Ψ̇2C21C33RxRzΘ̇Ψ̇+C222Ry2Ψ̇2+2C22C23RyRzΨ̇22C22C31RxRyΘ̇Ψ̇2C22C32Ry2Θ̇Ψ̇2C22C33RyRzΘ̇Ψ̇+C232Rz2Ψ̇22C23C31RxRzΘ̇Ψ̇2C23C32RyRzΘ̇Ψ̇2C23C33Rz2Θ̇Ψ̇+C312Rx2Θ̇2+2C31C32RxRyΘ̇2+2C31C33RxRzΘ̇2+C322Ry2Θ̇2+2C32C33RyRzΘ̇2+C332Rz2Θ̇2
E5=C112Rx2Ψ̇2+2C11C12RxRyΨ̇2+2C11C13RxRzΨ̇22C11C31Rx2Φ̇Ψ̇2C11C32RxRyΦ̇Ψ̇2C11C33RxRzΦ̇Ψ̇+C122Ry2Ψ̇2+2C12C13RyRzΨ̇22C12C31RxRyΦ̇Ψ̇2C12C32Ry2Φ̇Ψ̇2C12C33RyRzΦ̇Ψ̇+C132Rz2Ψ̇22C13C31RxRzΦ̇Ψ̇2C13C32RyRzΦ̇Ψ̇2C13C33Rz2Φ̇Ψ̇+C312Rx2Φ̇2+2C31C32RxRyΦ̇2+2C31C33RxRzΦ̇2+C322Ry2Φ̇2+2C32C33RyRzΦ̇2+C332Rz2Φ̇2
E6=C112Rx2Θ̇2+2C11C12RxRyΘ̇2+2C11C13RxRzΘ̇22C11C21Rx2Φ̇Θ̇2C11C22RxRyΦ̇Θ̇2C11C23RxRzΦ̇Θ̇+C122Ry2Θ̇2+2C12C13RyRzΘ̇22C12C21RxRyΦ̇Θ̇2C12C22Ry2Φ̇Θ̇2C12C23RyRzΦ̇Θ̇+C132Rz2Θ̇22C13C21RxRzΦ̇Θ̇2C13C22RyRzΦ̇Θ̇2C13C23Rz2Φ̇Θ̇+C212Rx2Φ̇2+2C21C22RxRyΦ̇2+2C21C23RxRzΦ̇2+C222Ry2Φ̇2+2C22C23RyRzΦ̇2+C232Rz2Φ̇2

Summarising E4, E5and E6, follows Eq. (68).

2Ekindm=C112+C212Rx2Ψ̇2+2C11C12+C21C22RxRyΨ̇2+2C11C13+C21C23RxRzΨ̇2+2C12C13+C22C23RyRzΨ̇2+C122+C222Ry2Ψ̇2+C132+C232Rz2Ψ̇22C21C31Rx2Θ̇Ψ̇2C21C32+C22C31RxRyΘ̇Ψ̇2C21C33+C23C31RxRzΘ̇Ψ̇2C22C33+C23C32RyRzΘ̇Ψ̇2C22C32Ry2Θ̇Ψ̇2C23C33Rz2Θ̇Ψ̇2C11C31Rx2Φ̇Ψ̇2C11C32+C12C31RxRyΦ̇Ψ̇2C11C33+C13C31RxRzΦ̇Ψ̇2C12C33+C13C32RyRzΦ̇Ψ̇2C12C32Ry2Φ̇Ψ̇2C13C33Rz2Φ̇Ψ̇+C112+C312Rx2Θ̇2+2C11C12+C31C32RxRyΘ̇2+2C11C13+C31C33RxRzΘ̇2+2C12C13+C32C33RyRzΘ̇2+C122+C322Ry2Θ̇2+C132+C332Rz2Θ̇22C11C21Rx2Φ̇Θ̇2C11C22+C12C21RxRyΦ̇Θ̇2C11C23+C13C21RxRzΦ̇Θ̇2C12C23+C13C22RyRzΦ̇Θ̇2C12C22Ry2Φ̇Θ̇2C13C23Rz2Φ̇Θ̇+C212+C312Rx2Φ̇2+2C21C22+C31C32RxRyΦ̇2+2C21C23+C31C33RxRzΦ̇2+2C22C23+C32C33RyRzΦ̇2+C222+C322Ry2Φ̇2+C232+C332Rz2Φ̇2E68

### 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 1and 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 1can be found on the left-hand side and the coefficients of kinetic energy (cf. Eq. (68)) of step2on the right-hand side. As described previously, the condition of identical energy is only proved, if all coefficients ar equal.

Rx2Ψ̇2:C322+C332=C112+C212s2Φc2Θ+c2Φc2Θ=c2Θc2Ψ+c2Θs2Ψc2Θ=c2ΘE69
RxRyΨ̇2:2C31C32=2C11C12+C21C22sΘcΘsΦ=sΦsΘcΘc2ΨcΦcΘsΨcΨ+sΦsΘcΘs2Ψ+cΦcΘsΨcΨsΘcΘsΦ=sΦsΘcΘE70
RxRzΨ̇2:2C31C33=2C11C13+C21C23sΘcΘcΦ=cΨsΨcΘsΦ+c2ΨcΘsΘcΦsΨcΨcΘsΦ+s2ΨcΘsΘcΦsΘcΘcΦ=cΘsΘcΦE71
RyRzΨ̇2:2C32C33=2C12C13+C22C23c2ΘsΦcΦ=s2ΨcΦsΦsΨcΨc2ΦsΘ+sΨcΨs2ΦsΘ+c2Ψs2ΘsΦcΦc2ΨcΦsΦ+sΨcΨc2ΦsΘsΨcΨs2ΦsΘ+s2Ψs2ΘsΦcΦ=1cΦsΦ+s2ΘsΦcΦc2ΘsΦcΦ=c2ΘsΦcΦE72
Ry2Ψ̇2:C312+C332=C122+C222s2Θ+c2Θc2Φ=s2Ψc2Φ2sΨcΨsΘcΦsΦ+c2Ψs2Θs2Φ+c2Ψc2Φ+2sΨcΨsΘsΦcΦ+s2Φs2Θs2Ψ=c2Φ+s2Θs2Φ=c2Φ1+s2Θs2Φ=c2Φc2Θ+c2Φs2Θ+s2Θ1s2Θc2Φs2Θ+c2Θc2Φ=c2Φc2Θ+s2ΘE73
Rz2Ψ̇2:C312+C322=C132+C232s2Θ+c2Θs2Φ=s2Ψs2Φ+2sΨcΨsΘsΦcΦ+c2Ψs2Θc2Φ+c2Ψs2Φ2sΨcΨsΘsΦcΦ+s2Θs2Ψc2Φ=s2Φ+s2Θc2Φ=s2Φ1+s2Θc2Φ=s2Φc2Θ+s2Φs2Θ+s2Θ1s2Θs2Φs2Θ+c2Θs2Φ=s2Φc2Θ+s2ΘE74
Rx2Θ̇Ψ̇:2C22C32+C23C33=2C21C31cΨcΘsΦcΦ+sΨsΘcΘs2ΦcΨcΘsΦcΦ.+sΨsΘcΘc2Φ=sΨcΘsΘsΨsΘcΘ=sΨcΘsΘE75
RxRyΘ̇Ψ̇:2C22C31+C21C32=2C21C32+C22C31E76
RxRzΘ̇Ψ̇:2C21C33+C23C31=2C21C33+C23C31E77
RyRzΘ̇Ψ̇:2C22C33+C23C32=2C22C33+C23C32E78
Ry2Θ̇Ψ̇:2C23C33+C21C31=2C22C32cΨcΘsΦcΦ+sΨsΘcΘc2ΦsΨsΘcΘ=cΨcΘcΦsΦsΨcΘsΘs2ΦsΨsΘcΘ1s2Φ1=sΨcΘsΘs2ΦsΨcΘsΘs2Φ=sΨcΘsΘs2ΦE79
Rz2Θ̇Ψ̇:2C22C32+C21C31=2C23C33cΨcΘsΦcΦ+sΨsΘcΘs2ΦsΨsΘcΘ=cΨcΘsΦcΦsΨsΘcΘc2ΦsΨsΘcΘ1c2Φ1=sΨcΘsΘc2ΦsΨcΘsΘc2Φ=sΨcΘsΘc2ΦE80
Rx2Φ̇Ψ̇:2C13C33+C12C32=2C11C31sΨcΘsΦcΦ+cΨsΘcΘc2ΦsΨcΘsΦcΦ+cΨsΘcΘs2Φ=cΨcΘsΘcΨsΘcΘ=cΨcΘsΘE81
RxRyΦ̇Ψ̇:2C12C31+C11C32=2C11C32+C12C31E82
RxRzΦ̇Ψ̇:2C11C33+C13C31=2C11C33+C13C31E83
RyRzΦ̇Ψ̇:2C12C33+C13C32=2C12C33+C13C32E84
Ry2Φ̇Ψ̇:2C13C33+C11C31=2C12C32sΨcΘsΦcΦ+cΨsΘcΘc2ΦcΨcΘsΘ=sΨcΘsΦcΦcΨsΘcΘs2ΦcΨsΘcΘ1s2Φ1=cΨsΘcΘs2ΦcΨsΘcΘs2Φ=cΨsΘcΘs2ΦE85
Rz2Φ̇Ψ̇:2C12C32+C11C31=2C13C33sΨcΘsΦcΦ+cΨsΘcΘs2ΦcΨcΘsΘ=sΨcΘsΦcΦcΨsΘcΘc2ΦcΨsΘcΘ1c2Φ1=cΨsΘcΘc2ΦcΨsΘcΘc2Φ=cΨsΘcΘc2ΦE86
Rx2Θ̇2:C222+C232=C112+C312c2Ψc2Φ+2sΨcΨsΘsΦcΦ+s2Φs2Θs2Ψ+c2Ψs2Φ2sΨcΨsΘsΦcΦ+s2Θs2Ψc2Φ=c2Ψc2Θ+s2Θc2Ψ+s2Θs2Ψ=c2Ψ1+s2Θs2Ψ=c2Ψs2Θ+c2Ψc2Θ+s2Θ1s2Θc2Ψ=c2Ψc2Θ+s2Θ=c2Ψc2Θ+s2ΘE87
RxRyΘ̇2:2C21C22=2C11C12+C31C32sΨcΨcΘcΦs2ΨcΘsΘsΦ=cΨsΨcΘcΦ+c2ΨsΘcΘsΦsΘcΘsΦs2ΨcΘsΘsΦ=1s2Ψ1sΘcΘsΦs2ΨcΘsΘsΦ=s2ΨsΘcΘsΦE88
RxRzΘ̇2:2C21C23=2C11C13+C31C33sΨcΨcΘsΦs2ΨsΘcΘcΦ=sΨcΨcΘsΦ+c2ΨcΘsΘcΦsΘcΘcΦs2ΨsΘcΘcΦ=1s2Ψ1sΘcΘcΦs2ΨsΘcΘcΦ=s2ΨsΘcΘcΦE89
RyRzΘ̇2:2C22C23=2C12C13+C32C33c2ΨsΦcΦsΨcΨsΘc2Φ+sΨcΨsΘs2Φs2Ψs2ΘsΦcΦ=s2ΨsΦcΦsΨcΨsΘc2Φ+sΨcΨsΘs2Φ+c2Ψs2ΘsΦcΦ+c2ΘsΦcΦc2ΨsΦcΦs2Ψs2ΘsΦcΦ=s2ΨsΦcΦ+c2Ψs2ΘsΦcΦ+c2ΘsΦcΦsΦcΦc2Ψs2Ψs2Θ=sΦcΦs2Ψ+c2Ψs2Θ+c2Θc2Ψs2Ψs2Θ=s2Ψ+c2Ψs2Θ+c2Θ1s2Ψs2Θ+c2Ψs2Θ=s2Ψ+c2Θ+c2Ψs2Θ=s2Ψ+c2Ψs2Θ+c2ΘE90
Ry2Θ̇2:C212+C232=C122+C322s2Ψc2Θ+c2Ψs2Φ2sΨcΨsΘsΦcΦ+s2Ψs2Θc2Φ=s2Ψc2Φ2sΨcΨsΘsΦcΦ+c2Ψs2Θs2Φ+c2Θs2Φs2Ψc2Θ+c2Ψs2Φ+s2Ψs2Θc2Φ=s2Ψc2Φ+c2Ψs2Θs2Φ+c2Θs2Φs2Ψc2Θ1+c2Ψ1s2Φ+s2Ψs2Θc2Φ=s2Ψ1c2Φ+c2Ψs2Θs2Φ+1c2Θs2Φs2Ψc2Θs2Φ1+s2Ψc2Θc2Φ2+c2Ψs2Θs2Φ3+c2Ψc2Θs2Φ4+s2Ψs2Θc2Φ5=s2Ψs2Θc2Φ5+s2Ψc2Θc2Φ2+c2Ψs2Θs2Φ3+s2Ψc2Θs2Φ1+c2Ψc2Θs2Φ4E91
Rz2Θ̇2:C212+C222=C132+C332s2Ψc2Θ+c2Ψc2Φ+2sΨcΨsΘsΦcΦ+s2Ψs2Θs2Φ=s2Ψs2Φ+2sΨcΨsΘsΦcΦ+c2Ψs2Θc2Φ+c2Θc2Φs2Ψc2Θ+c2Ψc2Φ+s2Ψs2Θs2Φ=s2Ψs2Φ+c2Ψs2Θc2Φ+c2Θc2Φs2Ψc2Θ1+c2Ψ1c2Φ+s2Ψs2Θs2Φ=s2Ψ1s2Φ+c2Ψs2Θc2Φ+1c2Θc2Φs2Ψc2Θs2Φ1+s2Ψc2Θc2Φ2+c2Ψs2Θc2Φ3+c2Ψc2Θc2Φ4+s2Ψs2Θs2Φ5=s2Ψs2Θs2Φ5+s2Ψc2Θs2Φ1+c2Ψs2Θc2Φ3+s2Ψc2Θc2Φ2+c2Ψc2Θc2Φ4E92
Rx2Φ̇Θ̇:2C13C23+C12C22=2C11C21sΨcΨs2Φ+s2ΨsΘsΦcΦc2ΨsΘsΦcΦ+sΨcΨs2Θc2ΦsΨcΨc2Φs2ΨsΘsΦcΦ+c2ΨsΘsΦcΦ+sΨcΨs2Θs2Φ=cΨsΨc2Θ
sΨcΨs2Φ+sΨcΨs2Θc2ΦsΨcΨc2Φ+sΨcΨs2Θs2Φ=sΨcΨ+sΨcΨs2Θ=sΨcΨ1+1c2Θ=sΨcΨc2Θ=cΨsΨc2ΘE93
RxRyΦ̇Θ̇:2C11C22+C12C21=2C11C22+C12C21E94
RxRzΦ̇Θ̇:2C11C23+C13C21=2C11C23+C13C21E95
RyRzΦ̇Θ̇:2C12C23+C13C22=2C12C23+C13C22E96
Ry2Φ̇Θ̇:2C11C21+C13C23=2C12C22cΨsΨc2ΘsΨcΨs2Φ+s2ΨsΘsΦcΦc2ΨsΘsΦcΦ+sΨcΨs2Θc2Φ=sΨcΨc2Φ+s2ΨsΘsΦcΦc2ΨsΘsΦcΦsΨcΨs2Θs2ΦcΨsΨc2ΘsΨcΨs2Φ+sΨcΨs2Θc2Φ=sΨcΨc2ΦsΨcΨs2Θs2ΦsΨcΨc2Θs2Φ+s2Θc2Φ=sΨcΨc2Φs2Θs2Φc2Θs2Φ+s2Θc2Φ=c2Φs2Θs2Φ1s2Θ1+c2Φ+s2Θc2Φ=c2Φ+s2Θ1+c2Φ=c2Φ+s2Θ1+1s2Φ=c2Φs2Θs2Φ=c2Φs2Θs2ΦE97
Rz2Φ̇Θ̇:2C11C21+C12C22=2C13C23sΨcΨc2ΘsΨcΨc2Φs2ΨsΘsΦcΦ+c2ΨsΘsΦcΦ+sΨcΨs2Θs2Φ=sΨcΨs2Φs2ΨsΘsΦcΦ+c2ΨsΘsΦcΦsΨcΨs2Θc2ΦcΨsΨc2ΘsΨcΨc2Φ+sΨcΨs2Θs2Φ=sΨcΨs2ΦsΨcΨs2Θc2Φ
sΨcΨc2Θc2Φ+s2Θs2Φ=sΨcΨs2Φs2Θc2Φc2Θc2Φ+s2Θs2Φ=s2Φs2Θc2Φ1s2Θ1+s2Φ+s2Θs2Φ=s2Φ+s2Θ1+s2Φ=s2Φ+s2Θ1+1c2Φ=s2Φs2Θc2Φ=s2Φs2Θc2ΦE98
Rx2Φ̇2:C122+C132=C212+C312s2Ψc2Φ2sΨcΨsΘsΦcΦ+c2Ψs2Θs2Φ+s2Ψs2Φ+2sΨcΨsΘsΦcΦ+c2Ψs2Θc2Φ=s2Ψc2Θ+s2Θs2Ψc2Φ+c2Ψs2Θs2Φ+s2Ψs2Φ+c2Ψs2Θc2Φ=s2Ψ+c2Ψs2Θ=s2Ψ1+1s2Ψs2Θ=s2Ψs2Θ+s2Ψc2Θ+s2Θs2Ψs2Θ=s2Ψc2Θ+s2Θ=s2Ψc2Θ+s2ΘE99
RxRyΦ̇2:2C11C12=2C21C22+C31C32sΨcΨcΘcΦc2ΨcΘsΘsΦ=cΨsΨcΘcΦ+s2ΨsΘcΘsΦsΘcΘsΦc2ΨcΘsΘsΦ=1c2Ψ1sΘcΘsΦc2ΨcΘsΘsΦ=c2ΨsΘcΘsΦE100
RxRzΦ̇2:2C11C13=2C21C23+C31C33sΨcΨcΘsΦc2ΨsΘcΘcΦ=sΨcΨcΘsΦ+s2ΨsΘcΘcΦsΘcΘcΦc2ΨcΘsΘcΦ=1c2Ψ1sΘcΘcΦc2ΨcΘsΘcΦ=c2ΨsΘcΘcΦE101
RyRzΦ̇2:2C12C13=2C22C23+C32C33s2ΨsΦcΦ+sΨcΨsΘc2ΦsΨcΨsΘs2Φc2Ψs2ΘsΦcΦ=c2ΨsΦcΦ+sΨcΨsΘc2ΦsΨcΨsΘs2Φ+s2Ψs2ΘsΦcΦ+c2ΘsΦcΦsΦcΦs2Ψc2Ψs2Θ=sΦcΦc2Ψ+s2Ψs2Θ+c2Θs2Ψc2Ψs2Θ=c2Ψ+s2Ψs2Θ+c2Θ=1+s2Ψ+s2Ψs2Θ+1s2Θ=s2Ψ+s2Θ1c2Ψ1s2Ψc2Ψs2Θ=s2Ψs2Θc2ΨE102
Ry2Φ̇2:C112+C132=C222+C322c2Ψc2Θ+s2Ψs2Φ+2sΨcΨsΘsΦcΦ+c2Ψs2Θc2Φ=c2Ψc2Φ+2sΨcΨsΘsΦcΦ+s2Ψs2Θs2Φ+c2Θs2Φc2Ψc2Θ+s2Ψs2Φ+c2Ψs2Θc2Φ=c2Ψc2Φ+s2Ψs2Θs2Φ+c2Θs2Φc2Ψc2Θ1+s2Ψ1s2Φ+c2Ψs2Θc2Φ=c2Ψ1c2Φ+s2Ψs2Θs2Φ+1c2Θs2Φc2Ψc2Θs2Φ1+c2Ψc2Θc2Φ2+s2Ψs2Θs2Φ3+s2Ψc2Θs2Φ4+c2Ψs2Θc2Φ5=c2Ψs2Θc2Φ5+c2Ψc2Θc2Φ2+s2Ψs2Θs2Φ3+s2Ψc2Θs2Φ4+c2Ψc2Θs2Φ1E103
Rz2Φ̇2:C112+C122=C232+C332c2Ψc2Θ+s2Ψc2Φ2sΨcΨsΘsΦcΦ+c2Ψs2Θs2Φ=c2Ψs2Φ2sΨcΨsΘsΦcΦ+s2Ψs2Θc2Φ+c2Θc2Φc2Ψc2Θ+s2Ψc2Φ+c2Ψs2Θs2Φ=c2Ψs2Φ+s2Ψs2Θc2Φ+c2Θc2Φc2Ψc2Θ1+s2Ψ1c2Φ+c2Ψs2Θs2Φ=c2Ψ1s2Φ+s2Ψs2Θc2Φ+1c2Θc2ΦE104000
c2Ψc2Θs2Φ1+c2Ψc2Θc2Φ2+s2Ψs2Θc2Φ3+s2Ψc2Θc2Φ4+c2Ψs2Θs2Φ5=c2Ψs2Θs2Φ5+c2Ψc2Θs2Φ1+s2Ψs2Θc2Φ3+s2Ψc2Θc2Φ4+c2Ψc2Θc2Φ2E104

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.

## Acknowledgments

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).

## References

1. 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. 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. 3. Kirchhoff G. Über die Bewegung eines Rotationskörpers in einer Flüssigkeit. Crelles Journal. 1869;71:237-273
4. 4. Wendel K. Hydrodynamische Massen und hydrodynamische Massenträgheitsmomente. Jahrbuch der Schiffbautechnischen Gesellschaft. 1950;207-255
5. 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. 6. Korte H. Track control of a towed underwater sensor carrier. In: IFAC Control in Transportation Systems. Braunschweig: IFAC; 2000. pp. 89-94
7. 7. Paschen M. Contribution to Predict Trajectories of Pelagic Fishing Gears Affected by Ship Manoeuvres (in German) [Dissertation]. Rostock: University of Rostock; 1982
8. 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. 9. Aamo O, Fossen TI. Finite element modelling of moored vessels. Mathematical and Computer Modelling of Dynamical Systems. 2001;7(1):47-75
10. 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. 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. 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. 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. 14. Kane TR, Levison DA. Dynamics: Theory and Applications. 1st ed. New York: McGraw-Hill; 1985
15. 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. 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. 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. 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. 19. ITTC. Report of the Committee on Waves. Final Report and Recommendations to the 23rd International Towing Tank Conference. 2002;505-551
20. 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. 21. Ginsberg JH. Advanced Engineering Dynamics. 2nd ed. New York: Cambridge University Press; 2004
22. 22. Woernle C. Dynamic of Multi-Body Systems. Rostock: University of Rostock, Lesson script; 2001
23. 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. 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. 25. Blass A, Gurevich Y. Matrix transformation is complete for the average case. SIAM Journal on Computing. 1995;24(1):3-29

Written By

Holger Korte, Sven Stuppe, Jan-Hendrik Wesuls and Tsutomu Takagi

Submitted: April 21st, 2017 Reviewed: October 4th, 2017 Published: December 20th, 2017