In this paper, we present a method for calculating reaction forces for a crane mounted on a ship moving in waves. The method is used to calculate the reaction forces between the crane base and the vessel deck. This includes the case where the crane is mounted on the platform that keeps the base of the crane horizontal when the vessel is moving in roll and pitch. The wave motion of the ship is modeled with force response amplitude operators (RAOs) based on the JONSWAP wave spectrum. The combined equations of motion for a vessel and a crane are derived using Kane’s equations of motion, where velocities and angular velocities are formulated in terms of twists, and the associated partial velocities and partial angular velocities are given as lines in Plücker coordinates. The unknown reaction forces are represented as wrenches and are determined using screw transformations. The method is used to study the effect of the roll and pitch compensation platform in numerical simulations. The efficiency of the platform is evaluated in terms of the magnitude of reaction forces and crane payload sway angles.

# Determination of Reaction Forces of a Deck Crane in Wave Motion Using Screw Theory OPEN ACCESS

**Andrej Cibicik**

**Geir O. Tysse**

**Olav Egeland**

^{1}Corresponding author.

Contributed by the Ocean, Offshore, and Arctic Engineering Division of ASME for publication in the Journal of Offshore Mechanics and Arctic Engineering. Manuscript received November 29, 2018; final manuscript received April 29, 2019; published online June 19, 2019. Assoc. Editor: Marcelo R. Martins.

*J. Offshore Mech. Arct. Eng*141(6), 061604 (Jun 19, 2019) (12 pages) Paper No: OMAE-18-1203; doi: 10.1115/1.4043701 History: Received November 29, 2018; Accepted April 29, 2019

Crane operations are important in the offshore industry as a vital element in the supply chain of offshore installations. Offshore cranes may be large and heavy structures, which apply significant dynamic reaction forces and moments to the structure of the vessel. For small vessels, the mass of a crane may be significant compared to the mass of the vessel, which means that the motion of the crane will influence the motion of the ship. In such cases, it will be useful to determine the reaction forces from a dynamic analysis, where the equations of motion for the combined crane and vessel system are used. The determined reaction forces will be needed in the design of the connection between the crane and the deck. Alternatively, the reaction forces can be used to determine the design specifications for a roll and pitch compensation platform, when it is used to keep the base of the crane horizontal, while the vessel is moving in waves. The dynamic analysis based on the combined equations of motion can be used to determine the reaction forces for the cases with and without the motion compensation platform, which may be useful for the evaluation of the platform efficiency.

The dynamic coupling between crane load and vessel motion was studied in Ref. [1] for a floating crane barge. The response of the barge was found when the barge and the crane are modeled as one rigid body. A simulator for a vessel with a crane was presented in Ref. [2], where the authors used an object-oriented modeling approach and co-simulation with the functional mock-up interface [3].

The combined dynamics of a ship and a deck crane were derived in Ref. [4], where Lie Groups and Lagrangian mechanics were used to derive the equations of motion for a general vehicle-manipulator system where the vehicle has six degrees-of-freedom (DOFs). Earlier work on the same topic was presented in Ref. [5], where the equations of motion of a spacecraft-manipulator system were derived using Kane’s equation of motion. Kane’s equations of motion [6] for multibody systems are based on the Newton–Euler formulation, where the forces of constraints are eliminated using the principle of virtual work. In Kane’s method, the principle of virtual work is formulated in terms of partial velocities and partial angular velocities, which serve as projection operators that project inertial and external forces onto the directions associated with the generalized speeds. In this way, a minimal set of ordinary differential equations is obtained. Kane’s method was originally formulated at the component level using coordinate-free vectors [6], and it was used in a very detailed derivation of the equations of motion for a Stanford manipulator in Ref. [7]. An alternative formulation of Kane’s equations of motion for serial-link mechanisms was achieved by Angeles and Ma [8], where a matrix formulation was introduced using the link twists to determine the projection operators in the form of a natural orthogonal complement. This was further developed in Ref. [9], where the partial velocities and the partial angular velocities used in the derivation of Kane’s equations of motion were represented by lines in Plücker coordinates [10], and this was demonstrated in the development of the equations of motion for the combined dynamics of a ship and a crane mounted on the ship.

The forces of constraint (i.e., the reaction forces) in a multibody system can be determined using Lagrange multipliers, which results in a differential–algebraic system of equations, which can be solved for the reaction forces [11,12]. This modeling approach might require stabilization [13], and the obtained reaction forces are related to the constraint equations and may need conversion to the defined coordinates. Another approach is to use a minimal set of ordinary differential equations as in Kane’s formulation, where the constraint forces are eliminated from the equations. The constraint forces can then be brought to evidence using auxiliary generalized speeds, which define fictitious velocities and partial velocities in the directions of the unknown reaction forces [6]. This was investigated for a knuckleboom crane in Refs. [14,15], where the method of auxiliary generalized speeds was formulated in terms of twists and screws, which led to the systematic derivation procedure of the auxiliary partial velocities used in the determination of the constraint forces.

In this work, we extend the results of Refs. [9,14], first by deriving the equations of motion for a crane and a vessel system including a roll and pitch compensation platform between the ship deck and the crane base. Next, a procedure is developed for the determination of the reaction forces between the crane and the platform and between the platform and the vessel. A new feature of the proposed method is that the auxiliary partial velocities are given by lines in Plücker coordinates, and the reaction forces are given in terms of wrenches that are transformed with screw transformations. This gives a formulation with a clear geometric interpretation that is helpful in the development of the dynamical model. The performance of the derived model is studied by numerical simulations where the ship motion in waves is simulated using the JONSWAP spectrum [16–18] in combination with force response amplitude operators (RAOs) [19,20]. In addition, we investigate the efficiency of the motion compensation platform. Efficiency is evaluated in terms of the magnitude of the reaction forces at the deck/platform and platform/crane interfaces and in terms of the payload sway angles.

The rest of this article is organized as follows: Section 2 presents the theoretical preliminaries. In Sec. 3, we show the detailed derivation of the dynamical model, while in Sec. 4, we provide the procedure for the determination of reaction forces. In Secs. 5 and 6, we present the simulation results and provide the discussion of the results. Section 7 presents conclusions.

In this section, we present the theoretical preliminaries of this work. First, we present the procedure for modeling a marine vessel in wave motion. Then, we introduce twists, which are screw representation of linear and angular velocities of rigid bodies, and wrenches, which are screw representation of forces and torques on rigid bodies. Then, we present a general method for modeling dynamics of open-chain manipulators, which is based on Kane’s method, where partial angular velocities and partial linear velocities are represented as lines or screws. The modeling procedures given in the preliminaries will later be used to model a coupled crane/vessel system.

In this section, the mathematical model of a ship moving in waves is presented. This is done by the well-established method where force RAOs are used to calculate the wave forces on the ship and the waves are described with the JONSWAP wave spectrum [9,17]. The JONSWAP wave spectrum describes the distribution of waves in the North Sea and is given by

*H*

_{s}is the significant wave height, which is an approximation of the mean peak-to-peak wave height of the highest one-third of the waves, ω

_{p}is the peak wave frequency, γ = 3.3 and

*S*(ω) is approximated with a discrete wave spectrum

*S*(ω

_{i}) where ω

_{i},

*i*= 1…

*N*, is selected randomly in the interval $[\omega \xafi\u2212\Delta \omega /2,\omega \xafi+\Delta \omega /2]$, where the frequency range is divided into

*N*equal intervals Δω with a mean value $\omega \xafi$ in each interval.

Short crested irregular waves are modeled by discretizing the direction χ of the waves into *M* equal intervals Δχ with a mean value $\chi \xafj$ of each interval. Then, the discrete angles χ_{j}, *j* = 1… *M*, are selected randomly in each interval $[\chi \xafj\u2212\Delta \chi /2,\chi \xafj+\Delta \chi /2]$, and the spreading function is set to

_{0}is the dominant wave direction. The resultant wave elevation for each frequency ω

_{i}is

*F*

_{k}(ω, χ) for degree-of-freedom

*k*is a transfer function, which is given in terms of its amplitude $|F(\omega ,\chi )|$ and phase $\u2220F(\omega ,\chi )$, which are calculated from the geometry of the hull. The wave forces in the degree-of-freedom

*k*are then determined by

**M**

_{0,A}=

**M**

_{0}+

**A**(∞), and

**D**=

**B**(∞). The matrices $A\xafr$, $B\xafr$, $C\xafr$ and the vector

**x**are the terms from the state-space model of the fluid memory effect [21]. The term

**M**

_{0}is the associated mass matrix,

**A**(ω) is the frequency-dependent added mass matrix,

**B**(ω) is the frequency-dependent damping matrix, and

**G**is the restoring force matrix [22,23]. The vector $\nu 00$ is a six-dimensional vector including the linear velocity of the origin of frame 0 and the angular velocity of the vessel expressed in the vessel-fixed frame

In this section, the necessary background in screw theory is presented [10,24,25]. Screws are used in the following sections to formulate the equations of motion with Kane’s method and to calculate the reaction forces and moments. The use of screws gives a clear geometric interpretation and transformation rules that are useful in the derivation.

Consider two coordinate frames *i* and *j*, and let $Rji$ be the rotation matrix from frame *i* to frame *j*, which satisfies $(Rji)\u22121=(Rji)T=Rij$. Let **u**^{j} be a vector given in the coordinates of *j*. Then, the same vector given in the coordinates of frame *i* will be $ui=Rjiuj$. Let $a^$ denote the skew-symmetric matrix form of a vector **a** = [*a*_{1} *a*_{2} *a*_{3}]^{T}, so that

**b**. It follows that $u^i=Rjiu^jRij$ [25]. Consider the screw

**u**

^{j}and

**w**

^{j}given in the coordinates of

*j*and referenced to the point

*B*. The screw will satisfy the screw transformation

*i*and is referenced to the point

*A*. The screw transformation $VABij$ is given in terms of screw reference transformation matrix $UABi$ from

*B*to

*A*and the screw coordinate transformation matrix $R\xafji$, where

*A*to

*B*. It is noted that $UABiR\xafji=R\xafjiUABj$, which gives the alternative expressions

A line can be described by the screw

*i*and referenced to the point

*A*. Here,

**a**

^{i}is the unit direction vector of the line in the coordinates of

*i*and $m/Ai=p^Aiai$ is the moment of the direction vector

**a**

^{i}about the reference point

*A*, where $pAi$ is the vector from

*A*to a point on the line. The elements of the line $L/Ai$ are the Plücker coordinates of the line [10]. It is noted that the line $L/Ai$ is a screw, and it satisfies the screw transformations (11). Moreover, if the line is through the reference point

*A*, then $pAi=0$ and $L/Ai=[(ai)T0T]T$.

Consider the screw

*h*is the pitch. The line $L/Ai$ is said to be the screw axis of $s/Ai$. It is seen that if σ = 1 and

*h*= 0, then $s/Ai=L/Ai$. A line is therefore said to be a screw with zero pitch. Moreover, if

*h*→ ∞ and σ = 1/

*h*, then $s/Ai\u2192N/Ai$, where

**a**

^{i}.

Let the displacement from frame *i* to frame *j* be described by a translation of the origin given by a vector $piji$ followed by a rotation $Rji$ about the translated origin. This means that the displacement is given by the homogeneous transformation matrix

*j*relative to frame

*i*in the coordinates of

*i*and $\omega iji$ is the angular velocity of

*j*relative to

*i*in the coordinates of

*i*. The time derivative (18) can be written as follows [25]:

*j*relative to frame

*i*referenced to the origin of frame

*i*, and given in the coordinates of

*i*, while $tij/jj$ is the same twist referenced to the origin of frame

*j*and given in the coordinates of

*j*. It is noted that a twist is a screw, and it satisfies the screw transformation (11). From Eq. (19), it is seen that twists in the matrix form are transformed according to

A wrench is a screw representation of forces and torques. The wrench with a line of action through the origin of frame *i* referenced to *i* and given in the coordinates of *i* is given as follows:

*i*and given in the coordinates of

*i*. The lines are through the origin of frame

*i*, $Lxi/ii$ is along the

*x*

_{i}axis, $Lyi/ii$ is along the

*y*

_{i}axis, and $Lzi/ii$ is along the

*z*

_{i}axis. Next, consider the screws with infinite pitch

*i*and given in the coordinates of

*i*. Then, the wrench (28) associated with the three orthogonal forces and the three orthogonal torques can be written as a sum as follows:

*a*

_{j}for

*j*= 1, 2, 3 stands for the axes

*x*

_{i},

*y*

_{i},

*z*

_{i}, $\rho fj$ are the magnitudes of the forces and $\rho mj$ are the magnitudes of the torques in

*x*

_{i},

*y*

_{i},

*z*

_{i}directions. The expression (31) can alternatively be written in the matrix form as follows:

**ρ**is the vector of the magnitudes of forces and torques

*j*and expressed in the coordinates of

*j*by

*j*and given in the coordinates of

*j*can then be found by

Suppose that the wrench $wi/mjj$ is acting on a rigid body *i* that is moving with a twist $tni/mjj$. Then, the power of the wrench on the twist is

**Π**is the interchange operator [10] defined by

Dynamics of a crane can be modeled in the same manner as dynamics of a robotic manipulator arm. Therefore, in this section, we present the detailed procedure for deriving the equations of motion for robotic manipulator arms, which is based on Ref. [14]. The equation of motion of a rigid link *i* is derived using the Newton–Euler approach. We describe dynamics about the COG and formulate equations in a convenient matrix form

*i*relative to the inertial frame, $\omega nii$ is the angular velocity of link

*i*relative to the inertial frame, and $Mmii$ is the inertia tensor of link

*i*. The external force $fmii$ and external moment $ni/mii$ are applied at the COG of link

*i*, $fmii(c)$, and $ni/mii(c)$ are the forces and moments of constraint. All terms are given in the coordinates of frame

*i*, which is fixed in link

*i*.

The twist of link *i* relative to the inertial frame *n* and expressed in the coordinates of *i* is given by

*i*and is given in the coordinates of

*i*. This means that $vni/ii$ is the velocity of the origin of frame

*i*. This twist can be transformed, so that it is referenced to the COG point

*m*

_{i}of link

*i*with the transformation

*i*. The matrix $Umi,ii$ is the screw reference transformation matrix (12)

*i*to the origin of frame

*i*in the coordinates of

*i*. In the same way, the twist (40) can be referenced to the origin of frame

*i*+ 1, which is fixed in link

*i*+ 1, through the transformation

*i*+ 1 by the screw coordinate transformation

*i*+ 1 to

*i*. The screw transformation for the twist (40) to be referenced to the COG of link

*i*+ 1 and to be given in the coordinates of

*i*+ 1 can then be written as follows:

*i*+ 1 relative to frame

*n*can be calculated recursively from the twist of link

*i*according to

The links are connected with joints of a single rotational degree-of-freedom. The axis of rotation of the joint between link *i* − 1 and link *i* passes through the origin of frame *i*. This means that the twist of link *i* relative to link *i* − 1 will be given as a rotation about a line that passes through the origin of frame *i*. Then, the twist of frame *i* relative to *i* − 1 with the reference to *i* and given in the coordinates of *i* is

*i*and

*u*

_{i}is the generalized speed of joint

*i*. The line of the joint

*j*$Lj/mii$ referenced to the COG of link

*i*and expressed in the coordinates of

*i*is

*n*

_{q}is a number of links. The twist of link

*i*relative to the inertial frame, referenced to the COG of link

*i*and expressed in the coordinates of

*i*, can be given as a sum as follows:

**u**is a vector of generalized speeds, where the elements

*u*

_{i}are given in Eq. (49). The derivative of (53) is given by

*j*line referenced to the COG of link

*i*and expressed in the coordinates of

*i*is defined as

*m*

_{i}to the origin of frame

*i*, $pi\u2212k+1,i\u2212ki$ is a distance vector from frame

*i*−

*k*+ 1 to frame

*i*−

*k*, and $\omega ^iji$ is a skew-symmetric form of the angular velocity of frame

*j*with respect to frame

*i*, all expressed in the coordinates of

*i*. The equations of motion for link

*i*(39) can be written as follows:

*i*and the auxiliary skew-symmetric matrix

**W**

_{i}are given by

*i*,

*m*

_{i}is the mass of link

*i*, and $\omega ^nii$ is the skew-symmetric form of the angular velocity vector taken from the first three rows of $tni/mii$. According to the D’Alembert’s principle and the principle of virtual work, constraint forces are cancelled from the equation of motion of the whole multibody system

We consider a mechanical system given in Fig. 1. The system is a knuckleboom crane with a pitch and roll compensation platform mounted on a vessel in wave motion. The system is described as a serial-link mechanism with eight links, representing the rigid bodies of the system. Body 0 is the vessel, bodies 1 and 2 describe the motion compensation platform, body 3 is the crane king, body 4 is the first boom, body 5 is the second boom, and bodies 6 and 7 describe the payload. Relative motion between bodies *i* and *i* + 1 is modeled with single degree-of-freedom joints. Each body *i* has a body-fixed frame *i* and a mass *m*_{i}. The COG of body *i* for *i* = 1…7 is located at a distance *d*_{i} from the origin of frame *i*.

The position and orientation of frame 0 relative to frame *n* (the inertial frame) is given by the vector

*x*

_{0}

*y*

_{0}

*z*

_{0}]

^{T}is the position of the origin of frame 0, ϕ is the roll angle, θ is the pitch angle, and ψ is the yaw angle. The rotation matrix from frame

*n*to frame 0 is

**R**

_{z}(ψ) is the rotation matrix for rotation ψ about the

*z*axis,

**R**

_{y}(θ) is the rotation matrix for rotation θ about the

*y*axis, and

**R**

_{x}(ϕ) is the rotation matrix for rotation ϕ about the

*x*axis. The rotation matrix from frame 0 to frame 1 is $R10=Rx(\pi )Rx(q1)$, where

*q*

_{1}is the rotation angle about the

*x*

_{1}axis. The rotation matrix from frame 1 to frame 2 is $R21=Ry(q2)$, where

*q*

_{2}is the rotation angle about the

*y*

_{2}axis. The rotation matrix from frame 2 to frame 3 is $R32=Rz(q3)$, where

*q*

_{3}is the rotation angle about the

*z*

_{3}axis. The rotation matrix from frame 3 to frame 4 is $R43=Rx(\pi 2)Rz(\pi 2)Rz(q4)$, where

*q*

_{4}is the rotation angle about the

*z*

_{4}axis. The rotation matrix from frame 4 to frame 5 is $R54=Rz(q5)$, where

*q*

_{5}is the rotation angle about the

*z*

_{5}axis. The rotation matrix from frame 5 to frame 6 is $R65=Rz(q6)$, where

*q*

_{6}is the rotation angle about the

*z*

_{6}axis. The rotation matrix from frame 6 to frame 7 is $R76=Ry(q7)$, where

*q*

_{7}is the rotation angle about the

*y*

_{7}axis.

The configuration of the system is given by the vector of generalized coordinates

*n*referenced to the origin of frame 0 and expressed in the coordinates of frame 0 and $ui=q\u02d9i$. The lines of the crane joints in the system are given as

*n*, and expressed in the coordinates of frame 0 is

*a*

_{j}for

*j*= 1, 2, 3 stands for the axes

*x*

_{0},

*y*

_{0}, and

*z*

_{0}, the terms $(\omega n00)aj$ and $(vn0/00)aj$ are the projections of $\omega n00$ and $vn0/00$, respectively, on

*a*

_{j}axis, while $Laj/00$ and $Naj/00$ are defined in Eqs. (29) and (30). The sum (72) can be written in a matrix form as follows:

*i*and expressed in the coordinates of

*i*by the screw transformation

*j*$Lj/mii$ referenced to the COG of body

*i*and expressed in the coordinates of

*i*can be obtained by $Lj/mii=Vmi,jijLj/jj$ as in Eq. (50). The projection matrix of body 0 is defined as

*i*for

*i*= 1…7 is defined as

*n*

_{q}= 7 is the number of the crane and platform DOFs. The twist of body

*i*referenced to the COG of body

*i*and expressed in the coordinates of

*i*can be given as a sum

*i*= 1…7 is defined as

We formulate the equation of motion for the whole system by summing up the equations of motion of separate bodies premultiplied with the transpose of projection matrices. The forces and moments of constraint are cancelled from the equation of motion for the whole system by D’Alambert’s principle and the principle of virtual work as shown in Eq. (61). The equations of motion of the vessel (7) are defined in terms of $\nu 00$(8); however, in this work, we use the twist $tn0/00$(71) to represent velocities and angular velocities of the vessel. Therefore, some of the matrices in Eq. (7) need to be rearranged before the vessel and crane models can be coupled into one system. The mass, Coriolis, and gravity force matrices are transformed as

**M**

_{0,A},

**D**, and

**G**are defined in Eq. (7) and

**Π**is the interchange operator (38). The terms from the state-space model of the fluid memory effect are transformed as

The mass matrix of the whole system can be found by

**D**

_{0}is given in (81) and

**D**

_{i}for

*i*= 1…7 is defined in Eq. (60). The matrix of centrifugal and Coriolis forces is

**W**

_{0}is given in Eq. (81) and

**W**

_{i}for

*i*= 1…7 is defined in Eq. (60). The vector of generalized gravitational and buoyancy forces is given by

**G**

_{0}is the linearized gravitational and buoyancy force matrix given in Eq. (81), the wrench $wi/mii(g)$ for

*i*= 1…7 is generally defined as

*i*and referenced to the origin of frame

*i*for

*i*= 0 and referenced to the COG of body

*i*for

*i*= 1…5. The equations of motion for the whole system can now be written as follows:

**A**

_{r},

**B**

_{r}, and

**C**

_{r}are the auxiliary matrices from the state-space model of the fluid memory effect (82),

**τ**

_{thr}is the wrench of thruster control inputs, and

**τ**

_{w}is the wrench of wave forces estimated by Eq. (6). The term

**J**(

**η**) is the velocity transformation matrix [17].

In this section, we present the procedure for determination of reaction forces (we will, in general, refer to both forces and moments by just writing *forces*) between the vessel and the motion compensation platform (i.e., in joint 1), as well as between the platform and the crane king (i.e., in joint 3).

Define a vector of unknown magnitudes of the reaction forces

**ρ**

_{1}= [ρ

_{11}…ρ

_{16}]

^{T}is a vector of magnitudes of the reaction forces in joint 1,

**ρ**

_{3}= [ρ

_{31}…ρ

_{36}]

^{T}is a vector of magnitudes of the reaction forces in joint 3, and

**Π**is given in (38). In both cases, first three magnitudes are forces and last three magnitudes are moments. Auxiliary velocities and angular velocities can be expressed in terms of auxiliary twists, and then, partial auxiliary velocities and partial auxiliary angular velocities can be expressed as lines or screws. In the original formulation [6,14], the auxiliary velocities are first formulated and then partial differentiation is performed with respect to auxiliary generalized speeds. In the procedure presented in this paper, we do the derivation using lines and screws directly.

The wrench at the origin of frame 1 associated with the reaction forces and moments in joint 1 is given as

*i*and expressed in the coordinates of frame

*i*

*i*and expressed in the coordinates of frame

*i*

*i*in the system is defined as

**G**

_{0}, $\eta 0n$, and $wi/mii(g)$ are the same as in Eq. (85). The auxiliary vector of generalized crane and motion compensation platform control forces is given by

**M**and

_{f}**C**are defined as

_{f}**D**

_{0},

**D**

_{i},

**W**

_{0}, and

**W**

_{i}are the same as in Eqs. (83) and (84). Provided that

**P**

_{c}_{0}=

**0**, then Eq. (99) is simplified to

**u**are obtained from the simulation of Eq. (88).

The dynamical model of the marine vessel with the deck crane and the procedure for the determination of reaction forces are implemented numerically in this section. The lengths of the crane bodies are given in Table 1. The COG distances for the crane bodies are defines as *d*_{i} = *l*_{i}/2 for *i* = 1…6 and *d*_{7} = *l*_{7}. The masses of the crane bodies are given in Table 2.

The implemented vessel model is the same supply vessel model as used in Ref. [9]. The main vessel dimensions are as follows: the length between perpendiculars is 82.5 m, the breadth is 8.0 m, and the draught is 6.0 m. The mass of the vessel is 6362 mT. The hydrodynamic coefficients, force RAOs, parameters of the radiation force model, and rigid body mass matrix of the vessel were taken from the marine systems simulator [19]. The wave parameters considered in the simulations are significant wave height *H*_{s} = 5 m and peak frequency ω_{p} = 1.26 rad/s.

The vessel is initialized at $\eta 0n=[000000]T$, and a proportional-derivative (PD) controller is implemented for surge, sway, and yaw control, where the control task is to keep the controlled DOFs at zero. The crane is initialized at *q*3 = −90 deg, *q*_{4} = −45 deg, *q*_{5} = −90 deg, *q*_{6} = −45 deg, and *q*_{7} = 0 deg. We have implemented a PD controller with gravity compensation to keep *q*_{i} for *i* = 3, 4, 5 at its initial values, and the pendulum DOFs *q*_{6}, *q*_{7} are unactuated.

The roll and pitch compensation platform is initialized at *q*_{1} = 0 deg and *q*_{2} = 0 deg. Two control strategies are assumed in this paper, where in both, a PD controller is implemented. The first strategy is when the motion platform is controlled to stay horizontal (mode “MC on”). The second strategy is when the motion platform is controlled to stay parallel to the vessel deck (mode “MC off”). The second control strategy is used to simulated the absence of motion compensation for the comparison purposes.

Two simulations are run with two different crane locations on the vessel deck; see Fig. 2. In location I, the crane is placed in the origin of frame 0, which is in the middle of the deck seen in the transverse section. In location II, the crane is moved along the *y*_{0} axis by a distance of *e*_{c} = 3.0 m.

The results of the numerical simulation of the system with the crane in location I (see Fig. 2) are presented in this section. The results in time series with an active roll and pitch compensation platform (mode “MC on”) are shown as solid lines, while the results with a locked roll and pitch compensation platform (mode “MC off”) are shown as dashed lines. Roll, pitch, and heave time histories are shown in Figs. 3 and 4, while surge, sway, and yaw are close to zero throughout the simulation, and those graphs are omitted. The time histories for DOFs of the roll and pitch compensation platform are shown in Fig. 5. The time histories for the payload angles θ_{1} and θ_{2} (see Fig. 6) are shown in Figs. 7 and 8. The reaction moments and reaction forces on platform from the deck given in the coordinates of frame 1 are shown in Figs. 9–11.

The results of the numerical simulation of the system with the crane in location II (see Fig. 2) are presented in this section. The results in time series with an active roll and pitch compensation platform (mode “MC on”) are shown as solid lines, while the results with a locked roll and pitch compensation platform (mode “MC off”) are shown as dashed lines. Roll, pitch, and heave time histories are shown in Figs. 15 and 16, while surge, sway, and yaw are close to zero throughout the simulation, and those graphs are omitted. The time histories for DOFs of the roll and pitch compensation platform are shown in Fig. 17. The time histories for the payload angles θ_{1} and θ_{2} (see Fig. 6) are shown in Figs. 18 and 19. The reaction moments and reaction forces on platform from the deck are shown in Figs. 20–22.

The comparison of the maximum values of payload angles throughout the simulation for crane locations I and II are given in Figs. 26 and 27. The active roll and pitch motion compensation platform resulted in reduction of the pendulum sway of 66.5% for θ_{1} and 99.7% for θ_{2} in crane location I. Reduction of the crane pendulum sway in location II for θ_{1} and θ_{2} is accordingly 68.5% and 97.7%.

The relative comparison of the normalized maximum values of the reaction forces and moments on the platform from the vessel deck for crane locations I and II are given in Figs. 28 and 29. The relative comparison of the normalized maximum values of the reaction forces and moments on the crane king from the platform for crane locations I and II are given in Figs. 30 and 31. Reduction of the reaction moment and force components due to the active roll and pitch compensation platform is summed up in Tables 3 and 4.

The active roll and pitch compensation platform provided significant reduction in the maximum magnitudes of the payload sway angles, reaction forces, and reaction moments. The reduction of the payload sway angles improves the operational weather window of the cranes, while the reduction of the reaction forces and moments leads to the benefits in structural design and fatigue lifetime.

We have presented a procedure for dynamic modeling of a coupled crane and vessel system when the vessel is moving in waves. We have included the case when a roll and pitch compensation platform is installed between the vessel and the crane. The kinematics of the model were derived by representing velocities and angular velocities of the bodies as twists and by representing partial velocities and partial angular velocities as lines in Plücker coordinates. In addition, we have presented a procedure for the determination of reaction forces in the deck/platform and platform/crane interfaces. The reaction forces were determined from algebraic relations, which were conveniently derived by representing the unknown reaction forces as wrenches. Since wrenches are screws, screw transformations were also used in the derivations.

The presented model was implemented, and the results of numerical simulations were provided. The analysis was carried out for two crane locations: in the middle of the deck and closer to the starboard. The simulation results were used to study the efficiency of a roll and pitch compensation platform installed between a crane and a vessel. The efficiency was evaluated in terms of the magnitude of the determined reaction forces and in terms of the magnitude of the payload sway angles. It was shown that the compensation of roll and pitch angles, such that the base of the crane stays horizontal, leads to significant reduction of both reaction forces and payload sway angles. The simulation results demonstrated that the proposed model can be used to determine dynamical forces for structural and fatigue analyses of the crane/vessel interface.

Further development of this work could be implementation of a crane/platform control system for damping out the payload oscillations. In addition, the crane model can be extended with a payload hoisting degree-of-freedom.

## References

*n*-Axis Serial Robotic Manipulators Using a Natural Orthogonal Complement,” Int. J. Rob. Res., 7(5), pp. 32–47. [CrossRef]

## References

*n*-Axis Serial Robotic Manipulators Using a Natural Orthogonal Complement,” Int. J. Rob. Res., 7(5), pp. 32–47. [CrossRef]

## Figures

Two locations of the crane considered in simulations, the offset *e*_{c} = 3.0 m along the *y*_{0} axis

Roll and pitch angles for crane location I. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode. The pitch data are superposed.

Heave of the vessel for crane location I. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode. The lines are superposed.

DOFs of the motion compensation platform for crane location I. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Orientation of the payload is described by the rotation R_{xn}(*θ*_{1})R_{y}′(*θ*_{2}) from the inertial frame

Payload orientation angle θ_{1} for crane location I. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Payload orientation angle θ_{2} for crane location I. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Reaction moments on platform for crane location I. Solid lines show results in “MC on” mode dashed lines in “MC off” mode.

Reaction forces on platform for crane location I. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Reaction force on platform for crane location I. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Reaction moments on crane for crane location I. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Reaction forces on crane for crane location I. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Reaction force on crane for crane location I. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Roll and pitch angles for crane location II. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode. The pitch data are superposed.

Heave of the vessel for crane location II. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode. The lines are superposed.

DOFs of the motion compensation platform for crane location II. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Payload orientation angle θ_{1} for crane location II. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Payload orientation angle θ_{2} for crane location II. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Reaction moments on platform for crane location II. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Reaction forces on platform for crane location II. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Reaction force on platform for crane location II. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Reaction moments on crane for crane location II. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Reaction forces on crane for crane location II. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

Reaction force on crane for crane location II. Solid lines show results in “MC on” mode and dashed lines in “MC off” mode.

## Tables

Reaction | τ_{7} | $ny1$ | $nz1$ | $fx1$ | $fy1$ | $fz1$ |
---|---|---|---|---|---|---|

Location I | 29.3 | 99.7 | 88.6 | 76.5 | 77.4 | 0.0 |

Location II | 42.2 | 97.7 | 90.0 | 81.5 | 84.2 | 0.0 |

Reaction | $nx3$ | $ny3$ | τ_{9} | $fx3$ | $fy3$ | $fz3$ |
---|---|---|---|---|---|---|

Location I | 99.6 | 24.5 | 99.8 | 76.3 | 99.2 | 0.0 |

Location II | 97.8 | 36.3 | 97.4 | 83.3 | 97.8 | 0.0 |

## Errata

Some tools below are only available to our subscribers or users with an online account.