# A Practical Approach to Numerically Predicting a Maneuvering Vessel in Waves Oriented to Maritime Simulator

Authors: Shengke Ni; Zhengjiang Liu; Yao Cai; Teng Zhang
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

To improve the behavioral realism of maritime simulator, the mathematical model of ship heave and pitch motion in regular waves is proposed. To avoid the influence of irregular frequency by two-dimensional Frank source and sink distribution method, the multiparameter conformal mapping method was adopted to solve the hydrodynamic problem of ship’s transverse sections, and then the integration of the hydrodynamic coefficients and the wave exciting forces for the whole ship hull was obtained using the Salvesen–Tuck–Faltinsen (STF) strip method. The Abkowitz model was used for the horizontal maneuverability motion equations, and the ship heave and pitch motion mathematical model in regular waves can be built considering the maneuvering factors. To verify the credibility and applicability of the improved model, the Mariner vessel with a Froude number of 0.2 was taken as the simulation case; when Mariner vessel sails in head waves, the maximum relative error between the calculated results and the experimental results is 23.6% for amplitude response operators of ship heave and pitch motion, and the calculated results and the experimental results have the same trend. When turning maneuvers in regular waves are carried out, the time history of ship pitch and pitch motions in regular waves does not show the sinusoidal periodic change any more, but the time history of heave and pitch motion still possesses periodic change to some extent, at which the period is about the time required for the turning motion to finish a circle. When turning maneuvers are carried out for 3000 seconds, the time required for computer operation is about 487.6 s, and the real-time requirements of maritime simulator are satisfied. Thus, the mathematical model of ship heave and pitch motion established in this paper can be effectively applied to the maritime simulator.


---

## Body

## 1. Introduction

As an efficient and economical virtual simulation training platform, maritime simulator illustrated in Figure1 has been widely used in navigation training, maritime analysis, and port evaluation, and the ship motion mathematical model is the key to improving behavioral realism for maritime simulator by Jin and Yin [1]. At present, most of the ship motion mathematical models for maritime simulator are generated using commercial software, which is quite expensive and difficult to update, and greatly affects the sustainable development of maritime simulator. Therefore, it is urgent to establish ship motion mathematical model with high accuracy and applicability.Figure 1
Schematic diagram of maritime simulator.The hydrodynamic coefficients and the wave exciting force coefficients of the ship can be obtained by empirical formula, and the mathematical model of ship heave and pitch motion in regular waves was built by Li and Yang [2]. By supposing the ship as “box ship,” the wave forces calculation model was established for dynamic positioning simulator. The empirical formula and the “box ship” assumption reduce the accuracy of ship motion mathematical model, which limits the scope of application of ship type in Qian et al. [3].With the rapid improvement of computer computing and storage capacity, the numerical prediction of ship motion based on two-dimensional (2D) and three-dimensional (3D) theory has been greatly developed. In the field of two-dimensional theory, Hou [4] adopted the ordinary strip method to predict ship heave and pitch motions, and numerical results of the S175 container ship can be obtained and analyzed under different wavelengths and wave heights. The ordinary strip method is only applicable to the head sea condition, which greatly limits the application scope. The STF strip method was presented to predict ship motions in regular waves by Salvesen et al. [5], and the comparisons between numerical results and the experimental results show good agreement for most cases. Although the hydrodynamic coefficients of transverse sections were calculated by Frank source and sink distribution method, they were easy to be affected by irregular frequency, which leads to numerical divergence. Guevel and Bougis [6] solved the radiation-diffraction problem of a ship sailing in waves via the three-dimensional frequency domain Green function method, and the numerical analysis of the mixed distribution method and the source distribution method was carried out. The geometric characteristics of the hull surface and flow field could be taken into account by the three-dimensional frequency domain method. Newman [7] found that the calculation of the three-dimensional frequency Green function with forward speed is quite complex and time-consuming, and the accuracy of the numerical prediction for ship motions was quite poor especially when the ship speed was high. In the three-dimensional frequency domain. Nakos [8] carried out the numerical prediction of ship heave and pitch motion responses via high-order Rankine source method in the paper [9]. In the time domain, Chen et al. [10] adopted three-dimensional Rankine source function to obtain the hydrodynamic problems and motions of ship in regular waves. However, radiation condition is hard to achieve in short wavelengths for Rankine source method. Sun et al. [11] established the hydrodynamic analysis model of the ship in still water by the three-dimensional transient free-surface Green function. Clément [12] discovered that the three-dimensional transient free-surface Green function can be a solution to an ordinary differential equation. And Li et al. [13] solved three-dimensional transient free-surface Green function by modified precise integration method. Sun and Ren [14] further extended the three-dimensional transient free-surface Green function to diffraction problem and carried out numerical predictions of ship motions in regular waves. The three-dimensional transient free-surface Green function was not stable for the evaluation the hydrodynamic forces of the ship with flare, which reduced the scope of application of the ship type. Datta et al. [15] carried out modifications to the ship hull which is non-wall side, so that the three-dimensional transient free-surface Green function can be applied to solve the hydrodynamic problem without numerical divergence.However, the ship needs maneuvering in waves to achieve the scheduled navigation status; the influence of maneuver factors on ship motions in regular waves should be considered effectively. Yasukawa [16, 17] conducted the turning test, stopping test, and zigzag test in regular waves for SR108 container ship model. The free-running model maneuvering test could obtain reliable results, but experimental methods are not easy to be carried out. Based on the two-dimensional Frank source sink distribution method, the wave forces and hydrodynamic forces acting on the ship were calculated and directly added to the ship maneuvering motion equations in still water, and maneuvering motion and swaying motion of ship in regular waves can be predicted by Zhu et al. [18]. The waves forces acting on the ship were computed on the basis of “box ship” assumption and were added to the maneuvering mathematical group (MMG) equations, where the mathematical model of ship motion with six degrees of freedom in waves was established for the maritime simulator by Zhang and Yin [19]. The interaction among ship, propeller, and rudder was not well considered in the model of MMG. Based on the high-order Rankine source method, Zhang et al. [20] established the boundary integral equation of disturbance velocity potential considering the horizontal maneuvering motion, and the coupled equations of ship maneuvering motion and swaying motion were solved using the two-time-scale approach. Wang et al. [21] carried out the CFD (computational fluid dynamics) simulations of free-running ship under different operational conditions. Although the complex characteristics of flow field were considered by CFD-based method, the computational effort required could still be enormous. Several hundreds of CPU (Central Processing Unit) hours can be taken to simulate simply several seconds maneuver, which is difficult to meet the real-time simulation requirements of maritime simulator.To verify the proposed ship motion mathematical model for maritime simulator, within two-dimensional strip theory, the multiparameter conformal mapping method is employed to solve the hydrodynamic forces and wave exciting forces acting on the ship’s transverse sections, and the influence of wave direction angle and wave length on ship heave and pitch motion is studied. The Abkowitz model is employed to predict the ship heave and pitch motion in waves under maneuvering conditions, and the effects of rudder angle on heave and pitch motion are studied. The ship heave and pitch motion model can be integrated into maritime simulator by Visual Studio software, where the visualization of ship motion mathematical model can be achieved. Mariner vessel and YUKUN vessel are taken as study cases, and comparisons are presented between calculated results and experimental data in the paper [5]. Thus, the ship motion model proposed for maritime simulator in this paper is validated.


## 2. Mathematical Formulations

### 2.1. Coordinate Systems

As shown in Figure2, there are two coordinates systems used for solving the ship hydrodynamic and motion problems; the space-fixed coordinate system OXYZ is adopted to keep track of position and Euler angles of the ship, OXY plane describes the mean free surface with oz-axis pointing upward. The body-fixed coordinate system OXYZ is coincident with space-fixed coordinate system OXYZ at time t = 0. The ox-axis points toward the bow, the oy-axis points toward the port side, and the origin o is fixed at the midship. In the body-fixed coordinate system OXYZ, the u and v are translational velocities, and r is rotational rate.ψ is heading angle of the bow. δ is rudder angle. χ is wave angle of incident wave, and χ=0° represents the wave angle propagating positively along the ox-axis.Figure 2
Definition of coordinate systems.The relationship of the abovementioned two coordinate systems at timet should be(1)X=xcosψ−ysinψ+xo,Y=xsinψ+ycosψ+yo,Z=z,where xo,yo,0 are the coordinates of origin o in the space-fixed coordinate system.In the space-fixed coordinate system, the wave elevationζIt at origin o is given in the indefinite water depth as [22](2)ζ1t=ζ0coskxocosχ+yosinχ−ωt,where ζ0 is the amplitude of incident wave; ω is the absolute frequency of incident wave, k=ω2/g denotes wave number; and g denotes gravitational acceleration.In the body-fixed coordinate system, the incident wave potentialϕIx,y,z,t can be(3)ϕIx,y,z,t=ζ0gωekzsinkϖ−ωt,where ϖ=xcosχ−ψ+ysinχ−ψ+xocosχ+yosinχ.


### 2.2. Establishment of Abkowitz Model

The Abkowitz model describes horizontal maneuvering motion, which includes surge, sway, and yaw motions [23]. In the body-fixed coordinate, the dynamic equations of horizontal maneuvering motion can be given as(4)m−Xu˙u˙=f1u,v,r,δ,W,m−Yv˙v˙+mxG−Yr˙r˙=f2u,v,r,δ,W,mxG−Nv˙v˙+Izz−Nr˙r˙=f3u,v,r,δ,W,where m denotes mass of ship; xG denotes the x-coordinate value of the ship’s center of gravity; Izz denotes the moment of inertia about z-axis; Xu˙, Yv˙, and Yr˙ denote added mass on the ox-axis and oy-axis; Nv˙ and Nr˙ denote add moment of inertia about the oz-axis; f1, f2, and f3 are nonlinear force functions with respect to u, v, r, δ, and W; u denotes surge velocity; v denotes sway velocity; r denotes yaw velocity; δ denotes rudder angle; and W is corresponding to second-order wave force.Under several assumptions,f1, f2, and f3 can be given as shown by Bai [24]. In order to solve equation (4), hydrodynamic coefficients need to be nondimensionalized by Fossen [23], and the differential equation of horizontal maneuvering motion can be solved accurately by the fourth-order Runge–Kutta method.


### 2.3. Multiparameter Conformal Mapping Transformation of Ship’s Transverse Sections

Figure3 is a diagram of conformal transformation between the plane x′o′y′ and the plane ξ′o′η′; the o′x′-axis points downward; o′y′-axis points starboard; yw denotes half-width of transverse section; D denotes the draft of transverse section; C0 denotes contour line of transverse section; n denotes the unit normal vector pointing into the transverse section; n⟶ denotes the unit normal vector of a quarter unit circle; γ denotes radian angle of a quarter unit circle; and R denotes the radius of unit circle.Figure 3
Conformal mapping transformation between half transverse section and a quarter unit circle.The outer part of the ship’s transverse section contourC0 in the z′ plane is transformed into the outer part of the unit circle in the ς′ plane by conformal mapping transformation:(5)z′=∑s=1M−1a2s−1ς′−2s+1,where z′=x′+iy′;ς′=ξ′+iη′=R⋅eiγ; 0≤γ≤π/2; a2s−1 is transformation parameters; and M is the number of conformal mapping transformation parameters.Let(6)w′=z′ς′=p′+iq′.Complex numberzi′ denotes the ith transformation point on the contour line C0, and complex number wi′ is ith complex number w′ and can be defined as(7)wi′=zi′ςi′=pi+iqi,i=1,2,3,…,M−1,M,where ςi′=eiγi, γi=i−1π/2M−1 is the ith radian angle in the first quadrant of the unit circle corresponding to the ith transformation point zi′, and pi and qi are the real part and the imaginary part of the complex number wi′, respectively.In combination with equations (5)–(7), the equation can be given as(8)pi=a−1+a1cos2γi+⋯+a2M−3cos2M−2γi.If equation (8) is solved, then M conformal transformation coefficients can be deduced as(9)a−1=12M−1p1+2∑i=2M−1pi+pM,a2l−1=1M−1p1+∑i=2M−1picosi−1M−1lπ+−1fpM,l=1,2,…,M−2,a2M−1=12M−1p1+2∑i=2M−1−1ipi+−1M−1pM.To obtain the conformal transformation coefficients, it is necessary to give the initial value of the transformation pointzi as(10)zi′=a−1ςi′+a1ςi′,where a−1=D+yw and a1=D−yw/2.Figure4 shows the calculation flow of M conformal mapping coefficients. Since the initial M transformation points are not on the data line of the ship’s transverse section, the conformal mapping coefficients can be solved iteratively by Westlake and Wilson [25].Figure 4
The flowchart of conformal mapping.In Figure4, σ is the designated error tolerance, which can be set according to the ship shape, scale, and other parameters.


### 2.4. Establishment of Ship Heave and Pitch Motion Equations in Regular Waves

The transverse section is assumed to carry out harmonic heave motion in the planex′o′y′ with an oscillation frequency ω and a small displacement amplitude h0:(11)x′=h0cosωt+ε0,where ε0 is the phase angle.Within the potential flow theory, radiation velocity potentialϕR in the ς′ plane can be transformed from z′ plane, and the following requirement should be satisfied with respect to the parameters R and γ [26]:(1)
The Laplace equation:(12)∂2∂R2+1R∂∂R+∂2∂γ2ϕR=0.(2)
Symmetrical equation abouto′x′-axis:(13)∂ϕR∂γ=0,γ=0,R>1.(3)
The linearized free-surface condition in infinite depth:(14)k∑i=0M−1−1i2i−1a2i−1R−2i−1⋅ϕR+∂ϕR∂γ=0,γ=π2,R>1.(4)
The boundary condition on the surface of transverse section:(15)∂ϕR∂R=x˙′∂x′∂R.(5)
The radiation condition at infinity:(16)ϕR⟶0,R⟶∞onz=0.In combination with equations (12)–(16), the radiation velocity potential ϕR can be obtained [27].From Bernoulli’s equation, the hydrodynamic forcefR on the transverse section of unit length in the o′x′-direction is obtained as(17)fR=−ρ∫C0∂ϕR∂tdy′.The hydrodynamic forcefR can also be expressed as(18)fR=−m′dx˙′dt−N′x˙′.In combination with equations (17) and (18), the hydrodynamic coefficients m′ and N′ can be solved.When ship advances with a constant forward speedU in regular waves, the linear differential equations of coupled heave and pitch motion can be given as(19)m+A33z¨+B33z˙+C33z+A35θ¨+B35θ˙+C35θ=Fzcosωet+εz,Jyy+A55θ¨+B55θ˙+C55θ+A53z¨+B53z˙+C53z=Fθcosωet+εθ,where θ is the pitch angle; Jyy is the moment of inertia about y-axis; ωe is the wave encounter frequency; εz is the phase angle of heave wave exciting force; εθ is the phase angle of pitch wave exciting force moment; A33, B33, A35, B35, A55, B55, A53, and B53 are the hydrodynamic coefficients; Fz and Fθ are the heave wave force amplitude and pitch wave moment amplitude, respectively, induced by incident wave; C53, C55, C33, and C35 are the hydrostatic restoring coefficients; for better consideration of speed and end of ship effects, based on STF strip method [5], the abovementioned coefficients can be obtained by integrating the hydrodynamic coefficients m′ and N′ of transverse section.To eliminate numerical divergence induced by the initial disturbance of the fluid, the right terms of equation (19) need to be multiplied by the ramp function fm [28], where fm can be set as(20)fm=121−cosπtTw,t<Tw,1,t≥Tw,where Tw can be set as two times of the incident wave period.When the initial timet = 0, the displacement and velocity value of the ship heave and pitch motion are set as z=0, z˙=0, θ=0 and θ˙; the time history of the ship heave and pitch motion can be obtained by solving equation (19) with the fourth-order Runge–Kutta method [4].When the ship advances in regular waves with a constant forward speed, the ship heave and pitch motion can be directly solved by equation (19).When the ship advances in regular waves under maneuvering conditions, equations (4) and (19) need to be coupled and solved to obtain the ship heave and pitch motion results.Figure5 shows the flowchart of numerical prediction of ship heave and pitch motion (Δt is time step).Figure 5
The flowchart for numerical prediction of ship heave and pitch motion.


## 2.1. Coordinate Systems

As shown in Figure2, there are two coordinates systems used for solving the ship hydrodynamic and motion problems; the space-fixed coordinate system OXYZ is adopted to keep track of position and Euler angles of the ship, OXY plane describes the mean free surface with oz-axis pointing upward. The body-fixed coordinate system OXYZ is coincident with space-fixed coordinate system OXYZ at time t = 0. The ox-axis points toward the bow, the oy-axis points toward the port side, and the origin o is fixed at the midship. In the body-fixed coordinate system OXYZ, the u and v are translational velocities, and r is rotational rate.ψ is heading angle of the bow. δ is rudder angle. χ is wave angle of incident wave, and χ=0° represents the wave angle propagating positively along the ox-axis.Figure 2
Definition of coordinate systems.The relationship of the abovementioned two coordinate systems at timet should be(1)X=xcosψ−ysinψ+xo,Y=xsinψ+ycosψ+yo,Z=z,where xo,yo,0 are the coordinates of origin o in the space-fixed coordinate system.In the space-fixed coordinate system, the wave elevationζIt at origin o is given in the indefinite water depth as [22](2)ζ1t=ζ0coskxocosχ+yosinχ−ωt,where ζ0 is the amplitude of incident wave; ω is the absolute frequency of incident wave, k=ω2/g denotes wave number; and g denotes gravitational acceleration.In the body-fixed coordinate system, the incident wave potentialϕIx,y,z,t can be(3)ϕIx,y,z,t=ζ0gωekzsinkϖ−ωt,where ϖ=xcosχ−ψ+ysinχ−ψ+xocosχ+yosinχ.


## 2.2. Establishment of Abkowitz Model

The Abkowitz model describes horizontal maneuvering motion, which includes surge, sway, and yaw motions [23]. In the body-fixed coordinate, the dynamic equations of horizontal maneuvering motion can be given as(4)m−Xu˙u˙=f1u,v,r,δ,W,m−Yv˙v˙+mxG−Yr˙r˙=f2u,v,r,δ,W,mxG−Nv˙v˙+Izz−Nr˙r˙=f3u,v,r,δ,W,where m denotes mass of ship; xG denotes the x-coordinate value of the ship’s center of gravity; Izz denotes the moment of inertia about z-axis; Xu˙, Yv˙, and Yr˙ denote added mass on the ox-axis and oy-axis; Nv˙ and Nr˙ denote add moment of inertia about the oz-axis; f1, f2, and f3 are nonlinear force functions with respect to u, v, r, δ, and W; u denotes surge velocity; v denotes sway velocity; r denotes yaw velocity; δ denotes rudder angle; and W is corresponding to second-order wave force.Under several assumptions,f1, f2, and f3 can be given as shown by Bai [24]. In order to solve equation (4), hydrodynamic coefficients need to be nondimensionalized by Fossen [23], and the differential equation of horizontal maneuvering motion can be solved accurately by the fourth-order Runge–Kutta method.


## 2.3. Multiparameter Conformal Mapping Transformation of Ship’s Transverse Sections

Figure3 is a diagram of conformal transformation between the plane x′o′y′ and the plane ξ′o′η′; the o′x′-axis points downward; o′y′-axis points starboard; yw denotes half-width of transverse section; D denotes the draft of transverse section; C0 denotes contour line of transverse section; n denotes the unit normal vector pointing into the transverse section; n⟶ denotes the unit normal vector of a quarter unit circle; γ denotes radian angle of a quarter unit circle; and R denotes the radius of unit circle.Figure 3
Conformal mapping transformation between half transverse section and a quarter unit circle.The outer part of the ship’s transverse section contourC0 in the z′ plane is transformed into the outer part of the unit circle in the ς′ plane by conformal mapping transformation:(5)z′=∑s=1M−1a2s−1ς′−2s+1,where z′=x′+iy′;ς′=ξ′+iη′=R⋅eiγ; 0≤γ≤π/2; a2s−1 is transformation parameters; and M is the number of conformal mapping transformation parameters.Let(6)w′=z′ς′=p′+iq′.Complex numberzi′ denotes the ith transformation point on the contour line C0, and complex number wi′ is ith complex number w′ and can be defined as(7)wi′=zi′ςi′=pi+iqi,i=1,2,3,…,M−1,M,where ςi′=eiγi, γi=i−1π/2M−1 is the ith radian angle in the first quadrant of the unit circle corresponding to the ith transformation point zi′, and pi and qi are the real part and the imaginary part of the complex number wi′, respectively.In combination with equations (5)–(7), the equation can be given as(8)pi=a−1+a1cos2γi+⋯+a2M−3cos2M−2γi.If equation (8) is solved, then M conformal transformation coefficients can be deduced as(9)a−1=12M−1p1+2∑i=2M−1pi+pM,a2l−1=1M−1p1+∑i=2M−1picosi−1M−1lπ+−1fpM,l=1,2,…,M−2,a2M−1=12M−1p1+2∑i=2M−1−1ipi+−1M−1pM.To obtain the conformal transformation coefficients, it is necessary to give the initial value of the transformation pointzi as(10)zi′=a−1ςi′+a1ςi′,where a−1=D+yw and a1=D−yw/2.Figure4 shows the calculation flow of M conformal mapping coefficients. Since the initial M transformation points are not on the data line of the ship’s transverse section, the conformal mapping coefficients can be solved iteratively by Westlake and Wilson [25].Figure 4
The flowchart of conformal mapping.In Figure4, σ is the designated error tolerance, which can be set according to the ship shape, scale, and other parameters.


## 2.4. Establishment of Ship Heave and Pitch Motion Equations in Regular Waves

The transverse section is assumed to carry out harmonic heave motion in the planex′o′y′ with an oscillation frequency ω and a small displacement amplitude h0:(11)x′=h0cosωt+ε0,where ε0 is the phase angle.Within the potential flow theory, radiation velocity potentialϕR in the ς′ plane can be transformed from z′ plane, and the following requirement should be satisfied with respect to the parameters R and γ [26]:(1)
The Laplace equation:(12)∂2∂R2+1R∂∂R+∂2∂γ2ϕR=0.(2)
Symmetrical equation abouto′x′-axis:(13)∂ϕR∂γ=0,γ=0,R>1.(3)
The linearized free-surface condition in infinite depth:(14)k∑i=0M−1−1i2i−1a2i−1R−2i−1⋅ϕR+∂ϕR∂γ=0,γ=π2,R>1.(4)
The boundary condition on the surface of transverse section:(15)∂ϕR∂R=x˙′∂x′∂R.(5)
The radiation condition at infinity:(16)ϕR⟶0,R⟶∞onz=0.In combination with equations (12)–(16), the radiation velocity potential ϕR can be obtained [27].From Bernoulli’s equation, the hydrodynamic forcefR on the transverse section of unit length in the o′x′-direction is obtained as(17)fR=−ρ∫C0∂ϕR∂tdy′.The hydrodynamic forcefR can also be expressed as(18)fR=−m′dx˙′dt−N′x˙′.In combination with equations (17) and (18), the hydrodynamic coefficients m′ and N′ can be solved.When ship advances with a constant forward speedU in regular waves, the linear differential equations of coupled heave and pitch motion can be given as(19)m+A33z¨+B33z˙+C33z+A35θ¨+B35θ˙+C35θ=Fzcosωet+εz,Jyy+A55θ¨+B55θ˙+C55θ+A53z¨+B53z˙+C53z=Fθcosωet+εθ,where θ is the pitch angle; Jyy is the moment of inertia about y-axis; ωe is the wave encounter frequency; εz is the phase angle of heave wave exciting force; εθ is the phase angle of pitch wave exciting force moment; A33, B33, A35, B35, A55, B55, A53, and B53 are the hydrodynamic coefficients; Fz and Fθ are the heave wave force amplitude and pitch wave moment amplitude, respectively, induced by incident wave; C53, C55, C33, and C35 are the hydrostatic restoring coefficients; for better consideration of speed and end of ship effects, based on STF strip method [5], the abovementioned coefficients can be obtained by integrating the hydrodynamic coefficients m′ and N′ of transverse section.To eliminate numerical divergence induced by the initial disturbance of the fluid, the right terms of equation (19) need to be multiplied by the ramp function fm [28], where fm can be set as(20)fm=121−cosπtTw,t<Tw,1,t≥Tw,where Tw can be set as two times of the incident wave period.When the initial timet = 0, the displacement and velocity value of the ship heave and pitch motion are set as z=0, z˙=0, θ=0 and θ˙; the time history of the ship heave and pitch motion can be obtained by solving equation (19) with the fourth-order Runge–Kutta method [4].When the ship advances in regular waves with a constant forward speed, the ship heave and pitch motion can be directly solved by equation (19).When the ship advances in regular waves under maneuvering conditions, equations (4) and (19) need to be coupled and solved to obtain the ship heave and pitch motion results.Figure5 shows the flowchart of numerical prediction of ship heave and pitch motion (Δt is time step).Figure 5
The flowchart for numerical prediction of ship heave and pitch motion.


## 3. Numerical Simulation and Result Analysis

### 3.1. Study Cases

Mariner vessel and YUKUN vessel are adopted as the study cases in this paper, and parameters of the ships are shown in Table1. L, B, and D denote length of ship, breadth of ship, and draft of ship, respectively; kyy denotes radius of pitch inertia moment; ∇ denotes volume of displacement; Cb denotes block coefficient; Ar denotes area of rudder; and D0 denotes diameter of propeller.Table 1
Main particulars of the Mariner vessel and YUKUN vessel.
ShipL (m)B (m)D (m)kyy∇ (m3)CbAr (m2)D0 (m)Mariner160.9323.178.230.25L185410.6130.016.71YUKUN105185.40.25L5710.20.5611.83.8


### 3.2. Conformal Mapping of the Transverse Section for Mariner Vessel

The Mariner vessel is divided into twenty-one transverse sections. The transverse sections are numbered from the stern to bow. No. 0 is for the first section numbered, No. 1 is the second transverse section, and so on. As shown in Figure6, “Lewis mapping” denotes mapping of transverse section by two conformal mapping coefficients (a−1 and a1), “Current mapping” denotes mapping of transverse section by multiparameter conformal mapping coefficients (M>2), and “Actual transverse section” denotes data points obtained from the actual transverse section contour line. The designated error tolerance σ is set as B/40 for “Current mapping.” The number of mapping coefficients for No. 4 and No. 5 transverse sections are 17 and 19, respectively, by “Current mapping.”Figure 6
Results of conformal mapping for transverse section No. 4 and No. 5.As seen from Figure6, for No. 4 and No. 5 of the ship’s transverse section, the data points of actual ship’s transverse section are all located on the lines generated by “Current mapping,” while only a few of the data points on of actual ship’s transverse section are located on the lines generated by the “Lewis mapping” [29]. Therefore, the actual ship’s transverse sections can be better fit with multiparameter coefficient conformal method, and an effective solver can be provided for ship hydrodynamic solution.


### 3.3. Numerical Simulation and Analysis of the Motions of Mariner Vessel in Regular Waves

To verify the reliability and accuracy of the mathematical model of ship motion in waves established in this paper, the numerical simulation of the motion for Mariner vessel in regular waves is carried out. Mariner vessel advances at constant speed Froude number Fn = 0.2 (Fn =U/gL).As shown in Figures7 and 8, z0 denotes amplitude of heave motion, z0/ζ0 denotes response amplitude operator of heave motion, θ0 denotes amplitude of heave motion, θ0L/2ζ0 denotes response amplitude operator of heave motion, and λ/L denotes wavelength to ship length ratio. “Current method” denotes numerical results obtained by the presented method in this paper; “OSM method” denotes numerical results obtained by OSM (ordinary strip method); “Salvesen method” denotes numerical results obtained by Salvesen et al. [5], where the hydrodynamic coefficients are solved by Frank source and sink method; and “Experiment” denotes experimental results given by Salvesen et al. [5].Figure 7
Heave response amplitude operators of Mariner vessel in head waves.Figure 8
Pitch response amplitude operators of Mariner vessel in head waves.In Figures7 and 8, for the response amplitude operator (RAO) of the heave motion of Mariner vessel, when the length to wavelength ratio λ/L is in the range of 0.8∼2.1, the numerical results of “Current method”; both “Salvesen method” and “OSM method” are in good agreement with the results of “Experiment.” Relative error between “Current method” and “Experiment” is within 8.8%, the relative error between the “Salvesen method” and “Experiment” is within 9.1%, and the relative error between “OSM method” and “Experiment” is within 9.7%. For the response amplitude operator of the pitch motion of the Mariner vessel, when the length to wavelength ratio λ/L is in the range of 0.8∼2.1, the relative error between “Current method” and “Experiment” is within 12.9%, relative error between the “Salvesen method” and “Experiment” is within 3.5%, and relative error between “OSM method” and “Experiment” is within 14.6%. The numerical results of “Current method” are in better agreement with the “Experiment” results than “OSM method,” especially in the peak area of pitch response amplitude operator curve; “OSM method” results are lower than other methods.As the wave angleχ is in the range of 0°∼360° and the amplitude of wave ζ0 is 1 m, Figures 9 and 10 show the amplitude of heave and pitch motion under different λ/L.Figure 9
Amplitude of heave motion for Mariner vessel under differentλ/L.Figure 10
Amplitude of pitch motion for Mariner vessel under differentλ/L.In Figure9, the amplitude of heave motion is symmetrical with respect to the wave angle χ=180°. With the increase of λ/L, the fluctuation range of the amplitude of heave motion varying with the wave angle is narrowed, and the amplitude of heave motion gradually tends to 1 m.In Figure10, the amplitude of the ship’s pitch motion is symmetrical with respect to the wave angle χ=180°, and the wave angle χ=90° is located at the trough of the pitch motion amplitude curve; due to the asymmetry of the hull geometry around midship, the amplitude of pitch motion is not 0, but the value is about 0.02°.For Figures11–16, the parameters of incident wave are set as follows: when λ/L = 1.0, the amplitude of incident wave is ζ0 = 1 m, the wave angle is χ=135° at t = 0; the rudder angles of Mariner vessel are δ=0°, δ=15°, and δ=25°, respectively.Figure 11
The time history of heave motion for Mariner vessel (δ=0°).Figure 12
The time history of pitch motion for Mariner vessel (δ=0°).Figure 13
The time history of heave motion for Mariner vessel (δ=15°).Figure 14
The time history of pitch motion for Mariner vessel (δ=15°).Figure 15
The time history of heave motion for Mariner vessel (δ=25°).Figure 16
The time history of pitch motion for Mariner vessel (δ=25°).As shown in Figures11–16, when the rudder angle δ=0° is applied to the Mariner vessel, the ship advances forward straightly with constant speed in regular waves, and the time histories of the heave and pitch motion change sinusoidal periodically. When a nonzero rudder angle is applied to Mariner vessel, the ship carries out turning motions in regular waves, the wave encounter frequency for Mariner vessel changes all the time, time history of the ship’s heave and pitch motion is not sinusoidal and still possesses periodicity to some extent, and its period is about the time needed for the Mariner vessel to finish a turning circle.HpZ440 workstation is used as the simulation platform in this paper. The hardware parameters are Intel® Xeon™ e5-1650v3 CPU @ 3.5 GHz, memory RAM is 32.00 GB. To verify the real-time performance of the proposed mathematical model of ship motion, the ship pitch and heave motion under a certain rudder angle for Mariner vessel are taken as an example; when the Mariner ship carried out motion for 3000 s, the time consuming of CPU is about 487.6 s, which is much lower than the actual time of the ship motion, thus providing a solid guarantee for the real-time performance of the maritime simulator.


### 3.4. Visualization of the Motions for the YUKUN Vessel in Regular Waves

Based on the Visual Studio software, the numerical simulation program of the mathematical model of the ship heave and pitch motion in regular waves is developed, which is connected to the simulation platform of the maritime simulator to realize the real-time solution and visualization of the ship heave and pitch motion model.YUKUN vessel advances straightly at the design speed of 16.7 kn (1 kn = 0.514 m/s) in regular waves. The parameters of incident waves are set as follows: the amplitude of wave isζ0 = 2 m; λ/L is 1; the wave angles are χ=0°, χ=90°, and χ=180°, respectively.From Figures17 and 18, when YUKUN ship advances straightly with constant speed in regular waves, under the action of different wave angles, the ship heave and pitch motion change sinusoidal period, and the numerical value does not diverge. It proves the robustness and applicability of the program developed in this paper, which can further model the motion of other ship types for the maritime simulator.Figure 17
The time history of heave motion for YUKUN vessel (χ=0°, χ=90°, and χ=180°).Figure 18
The time history of pitch motion for YUKUN vessel (χ=0°, χ=90°, and χ=180°).Under the above wave parameters and working conditions, when the wave angle isχ=180°, Figure 19 is the visual schematic diagram of the motion attitude of the YUKUN vessel on the simulation platform of the maritime simulator.Figure 19
The heave displacement of ship is 1.01 m and pitch displacement of ship is−1.96°, respectively, when t is 44.7 s.From Figure19, the ship motion mathematical model built in this paper could be applied to the simulation platform of maritime simulator.


## 3.1. Study Cases

Mariner vessel and YUKUN vessel are adopted as the study cases in this paper, and parameters of the ships are shown in Table1. L, B, and D denote length of ship, breadth of ship, and draft of ship, respectively; kyy denotes radius of pitch inertia moment; ∇ denotes volume of displacement; Cb denotes block coefficient; Ar denotes area of rudder; and D0 denotes diameter of propeller.Table 1
Main particulars of the Mariner vessel and YUKUN vessel.
ShipL (m)B (m)D (m)kyy∇ (m3)CbAr (m2)D0 (m)Mariner160.9323.178.230.25L185410.6130.016.71YUKUN105185.40.25L5710.20.5611.83.8


## 3.2. Conformal Mapping of the Transverse Section for Mariner Vessel

The Mariner vessel is divided into twenty-one transverse sections. The transverse sections are numbered from the stern to bow. No. 0 is for the first section numbered, No. 1 is the second transverse section, and so on. As shown in Figure6, “Lewis mapping” denotes mapping of transverse section by two conformal mapping coefficients (a−1 and a1), “Current mapping” denotes mapping of transverse section by multiparameter conformal mapping coefficients (M>2), and “Actual transverse section” denotes data points obtained from the actual transverse section contour line. The designated error tolerance σ is set as B/40 for “Current mapping.” The number of mapping coefficients for No. 4 and No. 5 transverse sections are 17 and 19, respectively, by “Current mapping.”Figure 6
Results of conformal mapping for transverse section No. 4 and No. 5.As seen from Figure6, for No. 4 and No. 5 of the ship’s transverse section, the data points of actual ship’s transverse section are all located on the lines generated by “Current mapping,” while only a few of the data points on of actual ship’s transverse section are located on the lines generated by the “Lewis mapping” [29]. Therefore, the actual ship’s transverse sections can be better fit with multiparameter coefficient conformal method, and an effective solver can be provided for ship hydrodynamic solution.


## 3.3. Numerical Simulation and Analysis of the Motions of Mariner Vessel in Regular Waves

To verify the reliability and accuracy of the mathematical model of ship motion in waves established in this paper, the numerical simulation of the motion for Mariner vessel in regular waves is carried out. Mariner vessel advances at constant speed Froude number Fn = 0.2 (Fn =U/gL).As shown in Figures7 and 8, z0 denotes amplitude of heave motion, z0/ζ0 denotes response amplitude operator of heave motion, θ0 denotes amplitude of heave motion, θ0L/2ζ0 denotes response amplitude operator of heave motion, and λ/L denotes wavelength to ship length ratio. “Current method” denotes numerical results obtained by the presented method in this paper; “OSM method” denotes numerical results obtained by OSM (ordinary strip method); “Salvesen method” denotes numerical results obtained by Salvesen et al. [5], where the hydrodynamic coefficients are solved by Frank source and sink method; and “Experiment” denotes experimental results given by Salvesen et al. [5].Figure 7
Heave response amplitude operators of Mariner vessel in head waves.Figure 8
Pitch response amplitude operators of Mariner vessel in head waves.In Figures7 and 8, for the response amplitude operator (RAO) of the heave motion of Mariner vessel, when the length to wavelength ratio λ/L is in the range of 0.8∼2.1, the numerical results of “Current method”; both “Salvesen method” and “OSM method” are in good agreement with the results of “Experiment.” Relative error between “Current method” and “Experiment” is within 8.8%, the relative error between the “Salvesen method” and “Experiment” is within 9.1%, and the relative error between “OSM method” and “Experiment” is within 9.7%. For the response amplitude operator of the pitch motion of the Mariner vessel, when the length to wavelength ratio λ/L is in the range of 0.8∼2.1, the relative error between “Current method” and “Experiment” is within 12.9%, relative error between the “Salvesen method” and “Experiment” is within 3.5%, and relative error between “OSM method” and “Experiment” is within 14.6%. The numerical results of “Current method” are in better agreement with the “Experiment” results than “OSM method,” especially in the peak area of pitch response amplitude operator curve; “OSM method” results are lower than other methods.As the wave angleχ is in the range of 0°∼360° and the amplitude of wave ζ0 is 1 m, Figures 9 and 10 show the amplitude of heave and pitch motion under different λ/L.Figure 9
Amplitude of heave motion for Mariner vessel under differentλ/L.Figure 10
Amplitude of pitch motion for Mariner vessel under differentλ/L.In Figure9, the amplitude of heave motion is symmetrical with respect to the wave angle χ=180°. With the increase of λ/L, the fluctuation range of the amplitude of heave motion varying with the wave angle is narrowed, and the amplitude of heave motion gradually tends to 1 m.In Figure10, the amplitude of the ship’s pitch motion is symmetrical with respect to the wave angle χ=180°, and the wave angle χ=90° is located at the trough of the pitch motion amplitude curve; due to the asymmetry of the hull geometry around midship, the amplitude of pitch motion is not 0, but the value is about 0.02°.For Figures11–16, the parameters of incident wave are set as follows: when λ/L = 1.0, the amplitude of incident wave is ζ0 = 1 m, the wave angle is χ=135° at t = 0; the rudder angles of Mariner vessel are δ=0°, δ=15°, and δ=25°, respectively.Figure 11
The time history of heave motion for Mariner vessel (δ=0°).Figure 12
The time history of pitch motion for Mariner vessel (δ=0°).Figure 13
The time history of heave motion for Mariner vessel (δ=15°).Figure 14
The time history of pitch motion for Mariner vessel (δ=15°).Figure 15
The time history of heave motion for Mariner vessel (δ=25°).Figure 16
The time history of pitch motion for Mariner vessel (δ=25°).As shown in Figures11–16, when the rudder angle δ=0° is applied to the Mariner vessel, the ship advances forward straightly with constant speed in regular waves, and the time histories of the heave and pitch motion change sinusoidal periodically. When a nonzero rudder angle is applied to Mariner vessel, the ship carries out turning motions in regular waves, the wave encounter frequency for Mariner vessel changes all the time, time history of the ship’s heave and pitch motion is not sinusoidal and still possesses periodicity to some extent, and its period is about the time needed for the Mariner vessel to finish a turning circle.HpZ440 workstation is used as the simulation platform in this paper. The hardware parameters are Intel® Xeon™ e5-1650v3 CPU @ 3.5 GHz, memory RAM is 32.00 GB. To verify the real-time performance of the proposed mathematical model of ship motion, the ship pitch and heave motion under a certain rudder angle for Mariner vessel are taken as an example; when the Mariner ship carried out motion for 3000 s, the time consuming of CPU is about 487.6 s, which is much lower than the actual time of the ship motion, thus providing a solid guarantee for the real-time performance of the maritime simulator.


## 3.4. Visualization of the Motions for the YUKUN Vessel in Regular Waves

Based on the Visual Studio software, the numerical simulation program of the mathematical model of the ship heave and pitch motion in regular waves is developed, which is connected to the simulation platform of the maritime simulator to realize the real-time solution and visualization of the ship heave and pitch motion model.YUKUN vessel advances straightly at the design speed of 16.7 kn (1 kn = 0.514 m/s) in regular waves. The parameters of incident waves are set as follows: the amplitude of wave isζ0 = 2 m; λ/L is 1; the wave angles are χ=0°, χ=90°, and χ=180°, respectively.From Figures17 and 18, when YUKUN ship advances straightly with constant speed in regular waves, under the action of different wave angles, the ship heave and pitch motion change sinusoidal period, and the numerical value does not diverge. It proves the robustness and applicability of the program developed in this paper, which can further model the motion of other ship types for the maritime simulator.Figure 17
The time history of heave motion for YUKUN vessel (χ=0°, χ=90°, and χ=180°).Figure 18
The time history of pitch motion for YUKUN vessel (χ=0°, χ=90°, and χ=180°).Under the above wave parameters and working conditions, when the wave angle isχ=180°, Figure 19 is the visual schematic diagram of the motion attitude of the YUKUN vessel on the simulation platform of the maritime simulator.Figure 19
The heave displacement of ship is 1.01 m and pitch displacement of ship is−1.96°, respectively, when t is 44.7 s.From Figure19, the ship motion mathematical model built in this paper could be applied to the simulation platform of maritime simulator.


## 4. Conclusions

(1)
In this paper, aiming at the problem that the empirical formula and three-dimensional theory are not suitable for ship motion modeling of maritime simulator, within two-dimensional strip theory, a new mathematical model of the ship heave and pitch motion in regular waves via multiparameter conformal mapping method has been established.(2)
To avoid the influence of the irregular frequency inherent in the two-dimensional Frank source and sink distribution method, this paper makes a high-precision fitting of the ship’s transverse section based on the multiparameter conformal mapping method, so as to provide a reliable solver for the study of ship’s hydrodynamic force and motion in waves within the two-dimensional strip theory.(3)
Based on the numerical calculation of Mariner vessel’s motion, the effects of different wavelength, wave angle, and rudder angle on the ship heave and pitch motion are analyzed. When Froude number of ship is set as 0.2, the relative error between the response amplitude operator of heave motion and the experimental value is 8.8%, which shows good agreement with the experimental results when Mariner vessel advances straightly in regular waves at a constant speed; under a certain rudder angle, the period of time history of the heave and pitch motion is about the time required for the ship to finish a circle in regular waves, and the consuming time of CPU is far less than the actual operation time of the ship, which effectively guarantees the real-time performance of the maritime simulator.(4)
The time history of heave and pitch motion for YUKUN vessel under various wave angles is calculated, and the stable and convergent numerical results are obtained. Based on Visual Studio, the simulation platform of maritime simulator is connected to realize the visualization of ship’s motion attitude. Thus, the mathematical model of ship heave and pitch motion proposed in present paper has achieved engineering application.(5)
Based on the theory of linear potential flow, this paper studies the ship heave and pitch motion in regular waves. In fact, when a ship sails in the ocean, the rolling motion of the ship is also present due to either direct wave forcing or interaction with other ship motions. Furthermore, during a turning circle due to the development of a rolling angle during this maneuver, especially if the ship’s speed during the maneuver is high and for special ship types (e.g., RO-RO ferries), this angle might be large enough to have a significant impact on ship motion. Because the viscous forces account for a large part, the theory of linear potential flow is difficult to meet the modeling requirements. It is necessary to incorporate more freedom for establishing mathematical model of ship motion under complex sea conditions and operational conditions.



---

*Source: 8361951-2020-10-26.xml*