Computational Model for Simulation of Lifeboat Free-Fall during Its Launching from Ship in Rough Seas

: In order to improve the accuracy of the freefall of lifeboat motion simulation in a ship life-saving simulation training system, a mathematical model using the strip theory and Kane’s method is established for the freefall of the lifeboat into the water from a ship. With the boat moving on a skid, the model of the ship’s maneuvering mathematical group (MMG) is used to model the motion of the ship in the waves. Based on the formula of elasticity and friction theory, the forces of the skid acting on the boat are calculated. When the boat enters the water, according to the analytical solution theory of slamming, the slamming force of water entry is solved. The simulation experiments are carried out by the established model. The results of the numerical simulation are compared with the calculation results of the hydrodynamics software Star CCM + at water entry under initial condition A in the paper. The position and velocity of the center of gravity of the boat, the angle, and velocity and acceleration of pitch calculated by the two methods are in good agreement. There is a little di ﬀ erence between the values of translation acceleration calculated by the two methods, which is acceptable. This shows that our numerical algorithm has good accuracy. A qualitative analysis is performed to ﬁnd the safe point of water entry under the condition of di ﬀ erent wave heights and two situations of a ship encountering waves. Finally, the model is applied to the ship life-saving training system. The model can meet the system requirements and improve the accuracy of the simulation.


Motivation
The lifeboat is the main life-saving equipment onboard ships. When a shipwreck accident occurs, the crew onboard can quickly escape from the ship by the lifeboat. During the freefall of the lifeboat, the hull of the boat usually slams into the water at high speed, with huge instantaneous impact pressure. If the crews make a mistake when launching the lifeboat, it can cause serious damage to the hull structure and threaten the personal safety of the crew [1]. According to the International Convention on Standards of Training, Certification and Watchkeeping for Seafarers, 1978, as amended in 1995 (STCW 78/95), the crew must be trained and pass a lifeboat assessment before boarding [2]. Training is limited by factors such as time and costs. In recent years, virtual reality technology has developed rapidly, and has been used in training for marine life-saving [3,4]. In order to improve the immersion experience and the reality of the ship's life-saving training system, a mathematical model is established for the motion of the lifeboat's freefall during its launch from the ship. Section 2 describes the mathematical model of the motion of boat on the skid. Section 3 describes the mathematical model of boat's water entry. The setup, results, and analysis of the simulation experiments are presented in Section 4. Section 5 describes an application of the mathematical model and Section 6 gives conclusions.

Motion of the Boat on the Skid
After the lifeboat is unhooked and released, the lifeboat slides down the skid away from the mother ship by its own gravity. When the center of gravity of the boat slides out of the lowest point of the skid, the boat begins to rotate owing to the vector of the gravity and the force of the skid not acting on the boat in a straight line. Figure 1 shows the central longitudinal section of the ship. There are three Cartesian coordinate systems:

Coordinate System
is the inertial coordinate system, with unit vectors of the three axes 1 , 2 , 3 ; the coordinate system of 0 0 0 0 is fixed on the ship; 0 is located at the center of gravity of the ship, the direction of 0 0 points to the bow, the direction of 0 0 points to the keel, and 0 , 0 , 0 are the unit vectors of the axes. The coordinate system of 1 1 1 1 is fixed on the skid and 1 is located at the center of the upper end of the skid, with 1 1 pointing down along the skid and 1 1 pointing up the perpendicular to the slope of the skid; 1 , 1 , 1 are the unit vectors of the three axes; is the angle between the plane of the skid and the horizontal plane. As shown in Figure 2, is the center of gravity of the lifeboat, 2 2 2 2 is the coordinate system attached to the boat, 2 is located at , and 1 , 2 , 3 are the unit vectors for the three axes.   Section 2 describes the mathematical model of the motion of boat on the skid. Section 3 describes the mathematical model of boat's water entry. The setup, results, and analysis of the simulation experiments are presented in Section 4. Section 5 describes an application of the mathematical model and Section 6 gives conclusions.

Motion of the Boat on the Skid
After the lifeboat is unhooked and released, the lifeboat slides down the skid away from the mother ship by its own gravity. When the center of gravity of the boat slides out of the lowest point of the skid, the boat begins to rotate owing to the vector of the gravity and the force of the skid not acting on the boat in a straight line. Figure 1 shows the central longitudinal section of the ship. There are three Cartesian coordinate systems:

Coordinate System
is the inertial coordinate system, with unit vectors of the three axes 1 , 2 , 3 ; the coordinate system of 0 0 0 0 is fixed on the ship; 0 is located at the center of gravity of the ship, the direction of 0 0 points to the bow, the direction of 0 0 points to the keel, and 0 , 0 , 0 are the unit vectors of the axes. The coordinate system of 1 1 1 1 is fixed on the skid and 1 is located at the center of the upper end of the skid, with 1 1 pointing down along the skid and 1 1 pointing up the perpendicular to the slope of the skid; 1 , 1 , 1 are the unit vectors of the three axes; is the angle between the plane of the skid and the horizontal plane. As shown in Figure 2, is the center of gravity of the lifeboat, 2 2 2 2 is the coordinate system attached to the boat, 2 is located at , and 1 , 2 , 3 are the unit vectors for the three axes. Figure 1. The three coordinate systems are the inertial coordinate system , the coordinate system 0 0 0 0 fixed on the ship, and the coordinate system 1 1 1 1 fixed on the skid.  Generalized coordinates q r of the boat in the coordinate system fixed on the skid and the coordinate system attached to the boat with the unit vectors b r (r = 1, 2, 3).

Mathematical Model of Boat
The relative position of o 0 and o 1 remains unchanged in o 0 x 0 y 0 z 0 . The vector o 0 o 1 is represented by r and a is the angle of r with the horizontal line. Based on the theorem for composited velocity and acceleration, the velocity and acceleration of o 1 are where v o 0 , ω o 0 are the velocity and acceleration of the ship. When considering only the ship's longitudinal motion, the ship's pitch angle is θ ; u and w are components of the ship's velocity in the direction of axis ox 0 and oz 0 . The velocity and acceleration of The velocity and acceleration of G b are The angular velocity of the rigid body is irrelevant to the choice of the base point, so θ. The angular velocity and acceleration of the boat are The lifeboat is divided into n cross sections with equal thickness along the length of the boat. The coordinates of the center of each section S k , (k = 1, . . . ., n) are expressed as x s k , 0, 0 in the coordinate system of the boat, and l b is the length of boat. The thickness of each section is t s = l b /n ; the velocity at the center of each cross section is Generalized coordinates are q r , generalized velocities are u r = . q r (r = 1, 2, 3), partial velocities associated with points G b , S k and the angular velocity are written as follows [30]: The force analysis of the boat at the skid is shown in Figure 3. Forces acting on the boat are gravity G, supporting force F n , and friction F f .   The generalized force of gravity and inertia are, respectively, where is the mass of the boat, is the acceleration of gravity, and 2 is the moment of inertia around 2 . For the supporting force , this paper regards the contact roller as a spring with stiffness coefficient and damping coefficient . The frictional contact force is modeled using a coefficient of . The coordinates of the center of each section ( = 1, . . . . , ) are ( 1 + 3 , 0, 2 + 3 ) in the boat coordinate system. We define them as follows: The velocity of each section relative to its contacting point is / ′ . The positions of each contact roller are known; there are contact rollers on each side of the skid. The coordinates of the top of the contact roller on each side of the axes of 1 1 and 1 1 at the skid coordinate system are , ( = 1, . . . , ) as − 2 3 < < + 2 3 . The contacting force between and contact rollers [31,32] where The associated generalized force acting on the by the skid is The combined generalized force of the skid that acted on the boat is = ∑ ( = 1,2,3). The equation of the boat motion is where = + . Finally, the fourth-order Runge-Kutta method is used to solve the differential equation. The generalized force of gravity and inertia are, respectively, where m b is the mass of the boat, g is the acceleration of gravity, and I 2 is the moment of inertia around oy 2 .
For the supporting force F n , this paper regards the contact roller as a spring with stiffness coefficient k r and damping coefficient c r . The frictional contact force is modeled using a coefficient of µ r . The coordinates of the center of each section S k (k = 1, . . . ., n) are q 1 + x s k cosq 3 , 0, q 2 + x s k sinq 3 in the boat coordinate system. We define them as follows: The velocity of each section relative to its contacting point is v s k /s k . The positions of each contact roller are known; there are n p contact rollers on each side of the skid. The coordinates of the top of the contact roller on each side of the axes of o 1 x 1 and o 1 z 1 at the skid coordinate system are x p i , z p i = 1, . . . , n p as x p i − t s 2 cosq 3 < α s k < x p i + t s 2 cosq 3 . The contacting force between S k and contact rollers [31,32] is where The associated generalized force acting on the S k by the skid is The combined generalized force of the skid that acted on the boat is F r c = F r s k (r = 1, 2, 3). The equation of the boat motion is where F r = F r c + F r G . Finally, the fourth-order Runge-Kutta method is used to solve the differential equation.

Mathematical Model of Ship Motion
Based on the model of MMG [33], the coordinate system is as shown in Figure 1. The ship dynamic equation is where the variables with subscript H are the forces and moments on the hull; the variables with subscript wave are wave forces. m is the mass of the ship; m x , m y and m z are, respectively, added masses in the direction of the axes ox 0 , oy 0 and oz 0 ; I xx , I yy and I zz are the rotational moment of inertia around axes ox 0 , oy 0 and oz 0 ; J xx , J yy and J zz are the added moment of inertia around axes ox 0 , oy 0 , oz 0 ; u, v, w are velocities in the direction of axes ox 0 , oy 0 and oz 0 ; and p, q, r are the angular velocities around axes ox 0 , oy 0 , oz 0 . m z = m ; other added masses are calculated as follows: where c b is the block coefficient, d is the draft, B is the ship width, and L p is the length between perpendiculars. I xx , I yy , I zz , J xx , J yy and J zz are calculated as follows: where c p is the prismatic coefficient. The kinematics equation is where x, y, z and ϕ, θ, ψ are displacements and Euler angles relative to the inertial coordinate system, respectively.
The forces and moments on the hull are where N vr , and N vr are hydrodynamic derivatives, which are easily calculated according to the ship parameters [34]; ∆ is the ship displacement, GM is the metacentric height, z H is the coordinate of action point of Y H in the direction of oz 0 , and x c is the distance between the center of gravity and the center.
For the wave force, it is estimated, based on the Frude-Krenov assumption, that the hull is simplified to a box, and the six-degrees-of-freedom wave force and moment are as follows [34]: where a is the amplitude of the wave, k is the number of waves, ω e is the encounter frequency, χ is the encounter angle, ρ is the water density, and L is the ship waterline length.
The motion parameters of the ship can be obtained by solving the differential Equations (16) and (19) using the fourth-order Runge-Kutta method.

Water Entry
The geometric orientation relationship of the inertial and boat coordinate system is represented by β 1 , β 2 , β 3 . The transformation relationship between the two coordinate systems is as follows, where The generalized coordinates are The generalized velocities are where ω b , v b are the velocity and angular velocity of the lifeboat, respectively. The relationships between generalized coordinates and generalized velocities are The acceleration and angular acceleration of the boat are The partial angular velocities and velocities of the boat are .
The generalized inertial force is where I 1 , I 2 , I 3 is the moment of inertia around ox 2 , oy 2 , oz 2 . The generalized force caused by the gravity of the boat is The generalized force caused by air drag is where v b/W is the vector difference between the lifeboat velocity and the wind velocity v ρ w is the air density, A w is the projected area of the lifeboat in the plane perpendicular to v b/W , and C w is the coefficient of air drag. The generalized force caused by fluid drag is where F D s k is the fluid drag acting on S k and v s k r is the partial velocity of S k . The velocity of S k is as follows: where v s k . (41) The components of the fluid drag acting on S k are as follows [22]: where ρ a is the density of seawater and A s k i is the area of s k perpendicular to b i below the water surface. The section s m has maximum area A s m is the fluid drag coefficient, and v s k /R i is the component of the velocity of the section s k relative to the wave surface on the axis of b i . The generalized force caused by buoyancy is where V s k is the volume of S k in water and F B s k is the buoyancy acting on s k . The generalized force caused by the slamming force is where F p s k is the slamming force that acted on section s k . Previous scholars adopted the momentum theory to calculate the slamming force. When an object enters the water, part of its initial momentum will be transferred to the surrounding water. Assuming that the momentum conversion process is irreversible, the slamming force acting on the boat can be calculated by the rate of its momentum change: where v 0 is the velocity of the object before entering the water and v is the velocity of the object after entering the water; m 0 is the mass of the object and m a is the added mass of the object. Another method is based on energy theory, which was derived by Wu [35]. It has more physical significance and is widely used [36][37][38]. The relationship between slamming force and added mass is In order to consider the effect of the wave on the boat, the incident wave force F in is added into the slamming force: where V sub is the submerged volume and a f is the fluid acceleration. In this paper, the slamming force F p s k is calculated as the sum of energy theory and incident wave force for the first time. p s k i are the components of F p s k on the axis of b i : where a f /s k is the wave surface acceleration at section s k at the coordinate system of the boat, a s k have coordinates (x g s k , y g s k , z g s k ) in an inertial coordinate system: The depth of s k is where C ii are elements of the matrix of [ N C b ], the subscript is its location in the matrix, h 0 is the vertical distance from the center of gravity to the keel, h 0 is 1.15 m in Figure 4, and η(x, t) is a known function describing the wave surface.
where are elements of the matrix of [ ], the subscript is its location in the matrix, ℎ 0 is the vertical distance from the center of gravity to the keel, ℎ 0 is 1.15 m in Figure 4, and ( , ) is a known function describing the wave surface.
The added mass is 2 = 3 in this paper. 3 depends on the depth and the instantaneous half width of at the wave surface: where ( , ℎ ) is half the width of at depth ℎ . When the depth of is greater than 1 , the added mass is calculated according to the half-width at depth 1 , as shown in Figure 4. We take the derivative of 3 to time 3 = 3 ℎ ℎ with ℎ = max ( 3 / , 0).  To calculate 1 ℎ , define where ℎ = − ( , ). The length of water entry in the direction of 1 is 1 : where C(x s k , h s k ) is half the width of s k at depth h s k . When the depth of s k is greater than d 1 , the added mass is calculated according to the half-width at depth d 1 , as shown in Figure 4. We take the derivative of m where The length of water entry in the direction of b 1 is l 1 : where m ax is the added mass of the boat in the direction of b 1 when the boat is completely in the water.
where k m is a coefficient depending on , 0).

The acceleration of the section is
where a s k where L f is the distance from the center of gravity of the boat to the bow. In summary, the equation of motion is The motion parameters of the lifeboat can be solved by Equation (64) by the fourth-order Runge-Kutta method.

Experimental Setup and Implementation
The basic information of the boat, ship, and skid is shown in Tables 1-3. Air drag coefficient Figure 5 gives the structural dimension diagram of the ship and boat in the simulation experiment.      The simulation experiment of the whole algorithm was realized by MATLAB. The flow of the algorithm consists of two stages that are separated by the time point of the boat out of the skid. In the first stage, the motion of the skid is calculated according to parameters of ship motion solved by Equations (16) and (19). The generalized contact force can be solved by Equation (14), the generalized force caused by gravity and inertia can be solved by Equations (9) and (10), and the generalized velocities and generalized coordinates of the boat can be solved by Equation (15). In the second stage, the generalized force caused by inertial, gravity, air drag, fluid drag, buoyancy, and slamming force can be solved by Equations (32) (26) and (27). For the two stages, 20 lifeboat segments were used. For convenience of analysis and understanding, the data of the boat's translation and pitch are shown in the coordinates obtained by the inertial coordinate system, rotating 180° around the axis of . The results of a simulation experiment by the numerical algorithm in this paper were compared with the results of Star CCM+ at the initial condition A, as shown in Table 4. The version of Star CCM+ used was 13.02. These computations were performed using an overlapping grid, as shown in The simulation experiment of the whole algorithm was realized by MATLAB. The flow of the algorithm consists of two stages that are separated by the time point of the boat out of the skid. In the first stage, the motion of the skid is calculated according to parameters of ship motion solved by Equations (16) and (19). The generalized contact force can be solved by Equation (14), the generalized force caused by gravity and inertia can be solved by Equations (9) and (10), and the generalized velocities and generalized coordinates of the boat can be solved by Equation (15). In the second stage, the generalized force caused by inertial, gravity, air drag, fluid drag, buoyancy, and slamming force can be solved by Equations (32), (34), (36), (38), (44) and (47); the generalized velocities of the boat can be solved by Equation (64); and the generalized coordinates of the boat can be solved by Equations (26) and (27). For the two stages, 20 lifeboat segments were used. For convenience of analysis and understanding, the data of the boat's translation and pitch are shown in the coordinates obtained by the inertial coordinate system, rotating 180 • around the axis of oy.
The results of a simulation experiment by the numerical algorithm in this paper were compared with the results of Star CCM+ at the initial condition A, as shown in Table 4. The version of Star CCM+ used was 13.02. These computations were performed using an overlapping grid, as shown in Figure 6a, in order to provide accurate wave representation at any location irrespective of the lifeboat position. The behavior of two fluids (air and water) in the same continuum is modeled by using the model of volume of fluid (Figure 6b). Due to the existence of two fluids in different phases, the Euler multiphase model is activated and the gravity model is used to take into account the gravity effect of the two fluids. The initial condition of water entry is obtained by our numerical algorithm: the boat enters the water at the crest of the wave, the pitch angle of the boat when entering the water is 50.7 • with horizontal, vertical, and rotational velocity of 5.9 m/s, −11 m/s, and 0.46 rad/s. Table 4. Details of initial condition A.

Item (Unit) Value
Wave height (m) 2 Period of wave (s) 8 Wind velocity (m/s) (0, 0, 0) Ship's initial position (m) (0, 0, 0) Ship's initial velocity (knot) (1, 0, 0) Situation of ship encountering waves (0: head the wave; 1: follow the wave) 1 using the model of volume of fluid (Figure 6b). Due to the existence of two fluids in different phases, the Euler multiphase model is activated and the gravity model is used to take into account the gravity effect of the two fluids. The initial condition of water entry is obtained by our numerical algorithm: the boat enters the water at the crest of the wave, the pitch angle of the boat when entering the water is 50.7° with horizontal, vertical, and rotational velocity of 5.9 m/s, -11 m/s, and 0.46 rad/s. Table 4. Details of initial condition A.

Comparison and Analysis of Experimental Results
The red line is the result of the whole motion process calculated by the model in this paper. The red blue is the result of the water entry process calculated by Star CCM+. Figure 7 shows the trajectory and pitch angle of the boat. The results show that the trajectory and pitch angle calculated by two methods are in good agreement. When the boat moves on the skid, the pitch angle begins to change with no obvious logic. The pitch angle begins to increase at the bottom of the skid. After the boat leaves the skid, the pitch angle of the boat gradually increases to about 53°. The pitch angle begin to decrease after increases in short time because the bow of the boat hits the water. The center of gravity of the boat is sinking, the bow of the boat is raised, resulting in the stern hitting the water. The buoyancy increases with the increase of the depth of the stern. Subsequently, the pitch angle of the boat begins to increase again, and the lifeboat begins to float. After it comes out of the water, the bow hits the water again.

Comparison and Analysis of Experimental Results
The red line is the result of the whole motion process calculated by the model in this paper. The red blue is the result of the water entry process calculated by Star CCM+. Figure 7 shows the trajectory and pitch angle of the boat. The results show that the trajectory and pitch angle calculated by two methods are in good agreement. When the boat moves on the skid, the pitch angle begins to change with no obvious logic. The pitch angle begins to increase at the bottom of the skid. After the boat leaves the skid, the pitch angle of the boat gradually increases to about 53 • . The pitch angle begin to decrease after increases in short time because the bow of the boat hits the water. The center of gravity of the boat is sinking, the bow of the boat is raised, resulting in the stern hitting the water. The buoyancy increases with the increase of the depth of the stern. Subsequently, the pitch angle of the boat begins to increase again, and the lifeboat begins to float. After it comes out of the water, the bow hits the water again.  Figure 8 shows the ratio of the acceleration of the center of gravity to . When the boat is on the skid, the acceleration value of the center of gravity of the boat begins to vibrate. When the boat is at the low end of the skid, the vibration frequency of the acceleration value increases and the amplitude decreases. When the boat leaves the skid, the horizontal acceleration of the boat is zero and the vertical acceleration is before entering the water. The horizontal acceleration of the boat is negative in the initial stage of water entry because of the slamming force and the fluid drag force,  Figure 8 shows the ratio of the acceleration of the center of gravity to g. When the boat is on the skid, the acceleration value of the center of gravity of the boat begins to vibrate. When the boat is at the low end of the skid, the vibration frequency of the acceleration value increases and the amplitude decreases. When the boat leaves the skid, the horizontal acceleration of the boat is zero and the vertical acceleration is g before entering the water. The horizontal acceleration of the boat is negative in the initial stage of water entry because of the slamming force and the fluid drag force, and the acceleration becomes positive when the boat floats up because of its large buoyancy. The vertical acceleration of the boat is almost positive because of the upward force, and its negative value caused by gravity and fluid drag when the boat floats up. There is a little difference between the acceleration values calculated by the two methods, especially at the initial stage of water entry. The difference is acceptable. This phenomenon may be due to the fact that the deformation of the water caused by boat is not considered in this algorithm. . Trajectory (a) and pitch (b) of the lifeboat freefall during its launch from the ship. Figure 8 shows the ratio of the acceleration of the center of gravity to . When the boat is on the skid, the acceleration value of the center of gravity of the boat begins to vibrate. When the boat is at the low end of the skid, the vibration frequency of the acceleration value increases and the amplitude decreases. When the boat leaves the skid, the horizontal acceleration of the boat is zero and the vertical acceleration is before entering the water. The horizontal acceleration of the boat is negative in the initial stage of water entry because of the slamming force and the fluid drag force, and the acceleration becomes positive when the boat floats up because of its large buoyancy. The vertical acceleration of the boat is almost positive because of the upward force, and its negative value caused by gravity and fluid drag when the boat floats up. There is a little difference between the acceleration values calculated by the two methods, especially at the initial stage of water entry. The difference is acceptable. This phenomenon may be due to the fact that the deformation of the water caused by boat is not considered in this algorithm.  Figure 9 shows the velocity of the center of gravity of the boat. The horizontal velocity of the boat gradually begins to increase at the skid, and remains the same in the air before entering the water; it begins to decrease due to the negative acceleration at the beginning of the water entry phase, and begins to increase when the boat floats up. The value of vertical acceleration is decreasing because the boat moves down the skid. The value of vertical acceleration decreases at the same rate because of g in the air before entering the water. Its value gradually increases until the sign is positive because of the slamming force and the buoyancy at the water entry phase; its sign becomes negative when the boat falls into the water again. The velocity values calculated by the two methods are in good agreement.  Figure 9 shows the velocity of the center of gravity of the boat. The horizontal velocity of the boat gradually begins to increase at the skid, and remains the same in the air before entering the water; it begins to decrease due to the negative acceleration at the beginning of the water entry phase, and begins to increase when the boat floats up. The value of vertical acceleration is decreasing because the boat moves down the skid. The value of vertical acceleration decreases at the same rate because of g in the air before entering the water. Its value gradually increases until the sign is positive because of the slamming force and the buoyancy at the water entry phase; its sign becomes negative when the boat falls into the water again. The velocity values calculated by the two methods are in good agreement.  Figure 10 shows the velocity and acceleration of pitch. When the boat is on the skid, the value of the velocity of the pitch begins to vibrate. When the boat leaves the skid, the velocity of the pitch remains the same before entering the water. The velocity of the pitch decreases to a negative value when the boat enters the water. Its sign is positive when the boat's bow begins to fall into the water again. The acceleration of the pitch has the same rule of change as the horizontal acceleration. The velocity and acceleration of pitch calculated by the two methods are in good agreement.
In summary, the position and velocity of the center of gravity of the boat, the angle, and  Figure 10 shows the velocity and acceleration of pitch. When the boat is on the skid, the value of the velocity of the pitch begins to vibrate. When the boat leaves the skid, the velocity of the pitch remains the same before entering the water. The velocity of the pitch decreases to a negative value when the boat enters the water. Its sign is positive when the boat's bow begins to fall into the water again. The acceleration of the pitch has the same rule of change as the horizontal acceleration. The velocity and acceleration of pitch calculated by the two methods are in good agreement. Figure 10 shows the velocity and acceleration of pitch. When the boat is on the skid, the value of the velocity of the pitch begins to vibrate. When the boat leaves the skid, the velocity of the pitch remains the same before entering the water. The velocity of the pitch decreases to a negative value when the boat enters the water. Its sign is positive when the boat's bow begins to fall into the water again. The acceleration of the pitch has the same rule of change as the horizontal acceleration. The velocity and acceleration of pitch calculated by the two methods are in good agreement.
In summary, the position and velocity of the center of gravity of the boat, the angle, and velocity and acceleration of pitch calculated by the two methods are in good agreement. There is a little difference between the values of translation acceleration calculated by the two methods. The difference is acceptable. This shows that the numerical algorithm in this paper has good accuracy. The boat experiences three water impacts (bow impact, stern impact, bow impact) in the previous analysis. In order to analyze the local acceleration of different positions in the boat, the local acceleration of the bow, midship, and stern is calculated by Equation (1). Figure 11 shows the acceleration of the bow (a,b), midship (c,d), and stern (e,f) in the boat coordinate system. The acceleration values in all figures are the ratio of actual values to the acceleration of gravity .
During the whole water entry process, the acceleration curves have three peaks because of three water impacts (Figure 11a,c) The first two peaks are close; there is a positive maximum acceleration in the direction of 3 at the bow at the first impact as the first peak of the curve (Figure 11a). There is a positive maximum acceleration in the direction of 3 at the midship at the second impact as the second peak of the curve (Figure 11c). There are two contrary extreme values of the acceleration at the stern in the direction of 3 (Figure 11e): a negative extremum due to the first impact, and a positive extreme value due to the second impact. The negative maximum acceleration in the direction of 1 at the bow, midship, and stern is because of the first impact (Figure 11b,d,f). The maximum acceleration is about 3 in the direction of 3 , derived from the bow's first impact. In summary, the position and velocity of the center of gravity of the boat, the angle, and velocity and acceleration of pitch calculated by the two methods are in good agreement. There is a little difference between the values of translation acceleration calculated by the two methods. The difference is acceptable. This shows that the numerical algorithm in this paper has good accuracy.
The boat experiences three water impacts (bow impact, stern impact, bow impact) in the previous analysis. In order to analyze the local acceleration of different positions in the boat, the local acceleration of the bow, midship, and stern is calculated by Equation (1). Figure 11 shows the acceleration of the bow (a,b), midship (c,d), and stern (e,f) in the boat coordinate system. The acceleration values in all figures are the ratio of actual values to the acceleration of gravity g.
During the whole water entry process, the acceleration curves have three peaks because of three water impacts (Figure 11a,c) The first two peaks are close; there is a positive maximum acceleration in the direction of b 3 at the bow at the first impact as the first peak of the curve (Figure 11a). There is a positive maximum acceleration in the direction of b 3 at the midship at the second impact as the second peak of the curve (Figure 11c). There are two contrary extreme values of the acceleration at the stern in the direction of b 3 (Figure 11e): a negative extremum due to the first impact, and a positive extreme value due to the second impact. The negative maximum acceleration in the direction of b 1 at the bow, midship, and stern is because of the first impact (Figure 11b,d,f). The maximum acceleration is about 3 g in the direction of b 3 , derived from the bow's first impact.

Qualitative Analysis
According to the numerical calculation method used in this paper, the local acceleration extremum is calculated under different initial conditions. The wave heights are 1, 2, 3, 4, and 5 m. By adjusting the initial phase of the wave, the points of water entry of the boat are shown in Figure 12. Figure 12a shows the ship heading the wave and Figure 12b shows the ship following the wave. When the ship is heading or following the wave, only the surge, heave, and pitch of the ship are considered. Other initial conditions (period of wave, wind velocity, ship's initial position, ship's initial velocity) are the same as in condition A. Figures 13-17 show the trajectory and pitch angle of the lifeboat when the ship is heading the wave. Figure 18 shows the maximum value of the acceleration in the directions of b 1 and b 3 at the bow (Figure 18a,b), midship (Figure 18c,d), and stern (Figure 18e,f) of the boat at eight different points of water entry when the ship is heading waves. Figures 19-23 show the trajectory and pitch angle of the lifeboat when the ship is following waves. Figure 24 shows the maximum value of the acceleration in the directions of b 1 and b 3 at the bow (Figure 24a,b), midship (Figure 24c,d), and stern (Figure 24e,f) of the boat at eight different points of water entry when the ship is following waves. As there are no experimental data for comparison, the numerical calculation results may not be very accurate, but can be used in a qualitative analysis.

Qualitative Analysis
According to the numerical calculation method used in this paper, the local acceleration extremum is calculated under different initial conditions. The wave heights are 1, 2, 3, 4, and 5 m. By adjusting the initial phase of the wave, the points of water entry of the boat are shown in Figure 12. Figure 12a shows the ship heading the wave and Figure 12b shows the ship following the wave. When the ship is heading or following the wave, only the surge, heave, and pitch of the ship are considered. Other initial conditions (period of wave, wind velocity, ship's initial position, ship's initial velocity) are the same as in condition A.  Figures 13-17 show the trajectory and pitch angle of the lifeboat when the ship is heading the wave. Figure 18 shows the maximum value of the acceleration in the directions of 1 and 3 at the bow (Figure 18a,b), midship (Figure 18c,d), and stern (Figure 18e,f) of the boat at eight different points of water entry when the ship is heading waves. Figures 19-23 show the trajectory and pitch stern of the boat appear at point P4; the minimum value generally appears at P1 and P8. The maximum values at the bow and midship are almost the same, slightly lower than the stern. In the direction of 3 , the maximum absolute values of acceleration at the bow and midship of the boat appear at points P3 and P2; the minimum value generally appears at points P7 and P8. There is no obvious rule about the values of acceleration at the stern; the acceleration is the largest at the bow and the smallest at the midship.       The ship follows the wave. As shown in Figures 21-23a, when the height of the wave exceeds 3 m, the boat will move to the side of the ship after water entry at points P1 and P8. As shown in Figures 22-23a, when the height of the wave exceeds 4 m, the same phenomenon appears at point P7. As shown in Figures 19-23b, the pitch angle of the boat has a large range of change at points P5 and P6. As shown in Figure 24, in the direction of 1 , the maximum values of acceleration at the bow and midship appear at points P6 and P7, and the minimum values generally appear at points P2 and       (e) (f) Figure 24. The maximum acceleration of the bow, midship, and stern in the coordinate system of the lifeboat in the direction of (a,c,e), and in the direction of (b,d,f), at different wave heights.
Therefore, when the ship is heading the wave, P7 and P8 can be selected to enter the water at a low wave. In order to prevent the boat from moving to the side, one can choose point P1 at a high wave, which will withstand a large change in pitch. When the ship is following the wave, P1 and P8 can be selected at a low wave. In order to prevent the boat moving to the side, one can choose P2 and The ship heads the wave. As shown in Figures 15-17a, when the wave height exceeds 3 m, the boat will move to the side of the ship after water entry at P7 and P8. As shown in Figures 16a and  17a, when the wave height exceeds 4 m, the same phenomenon appears at point P6. As shown in Figures 13-17b, at points P1 and P8, the pitch angle of the boat has a large range of change. As shown in Figure 18, in the direction of b 1 , the maximum values of acceleration at the bow, midship, and stern of the boat appear at point P4; the minimum value generally appears at P1 and P8. The maximum values at the bow and midship are almost the same, slightly lower than the stern. In the direction of b 3 , the maximum absolute values of acceleration at the bow and midship of the boat appear at points P3 and P2; the minimum value generally appears at points P7 and P8. There is no obvious rule about the values of acceleration at the stern; the acceleration is the largest at the bow and the smallest at the midship.
The ship follows the wave. As shown in Figures 21-23a, when the height of the wave exceeds 3 m, the boat will move to the side of the ship after water entry at points P1 and P8. As shown in Figures 22 and 23a, when the height of the wave exceeds 4 m, the same phenomenon appears at point P7. As shown in Figures 19-23b, the pitch angle of the boat has a large range of change at points P5 and P6. As shown in Figure 24, in the direction of b 1 , the maximum values of acceleration at the bow and midship appear at points P6 and P7, and the minimum values generally appear at points P2 and P3. The maximum values of acceleration at the stern appear at points P4 and P5, and the minimum values generally appear at points P1 and P2. The maximum values at the bow and midship of the boat are almost the same, and slightly lower than for the stern. In the direction of b 3 , the maximum values of acceleration at the bow and midship of the boat appear at points P4 and P5; the minimum value generally appears at points P1 and P8. The maximum absolute values of the two extremal values at the stern appear at P4 and P5. The minimum values generally appear at P1 and P8; the acceleration is the largest at the bow and the smallest at the midship.
Therefore, when the ship is heading the wave, P7 and P8 can be selected to enter the water at a low wave. In order to prevent the boat from moving to the side, one can choose point P1 at a high wave, which will withstand a large change in pitch. When the ship is following the wave, P1 and P8 can be selected at a low wave. In order to prevent the boat moving to the side, one can choose P2 and P3 at a high wave. In summary, it is safer to select points at the crest and behind the crest when entering the water.

Application
This paper applies an established mathematical model to a ship's lifesaving training system, improving the immersion and reality of the system. The system takes the Panama bulk carrier as a physical prototype, and uses 3D Studio Max (3ds Max) to build 3D models of the ship, lifeboat, and skid. The system constructs a virtual scene of a ship life-saving drill, and releases the lifeboat through a three-dimensional virtual operation, as shown in Figure 25.

Conclusions
This paper presents a computational model for the simulation of lifeboat freefall during its launching from a ship in rough seas. The model is applied to the ship's lifesaving training system. We can draw the following conclusions: