Ship Loading Influence on the Slamming Impact of Typical Sections of an S-175 Container Ship

This paper investigates the influence of ship-loading condition on slamming during water entry. Three typical sections of the S-175 container ship, namely the bow, parallel middle body and stern, under three different loading conditions are studied. Full-sized models are established and simulated by commercial software LS-DYNA based on the explicit finite element method (FEM) using the arbitrary Lagrangian–Eulerian (ALE) algorithm. At first, validation is carried out by simulating the Wave Induced Loads on Ships Joint Industry Project II (WILS JIP-II) ship section entering the water and by verifying that the response is in good agreement with published experimental data. Then, nine different cases with three typical sections of the container ship and three different loadings, including the no-load (lightship weight), half-load and full-load weights of the ship, are investigated. Finally, the influence of the ship loading and sectional shape on the water impacts is analyzed and discussed. The present study is useful for the analysis of loading effects on ship slamming at the early stage of ship design.


Instruction
Slamming consists of a powerful water impact between ship hulls and waves. The hydrodynamic force caused by slamming may lead to springing of the ship [1]. For container ships, the ongoing drive to increase container carrying capacity has resulted in rapid increase of ship size, which makes them susceptible to whipping and springing [2,3]. Prediction of slamming impact pressures is important in the structural design of all vessels, especially high-speed craft. The time interval of water entry is extremely short. Thus, large impulsive pressure loads may cause local deformation or failure of the ship structure. In rough sea, the slamming phenomenon cannot be avoided even at the stern or midship. For the S-175 container ship, its bow and stern suffer mainly from flare slamming, while for the midship, bottom slamming is a concern.
In the 1930s, Von Karman [4] and Wagner [5] developed theoretical formulas for the maximum pressures at the bottom of the float during landing, which agreed very well with the experimental results. Zhao et al. [6] developed two different theoretical methods for predicting slamming loads on two-dimensional sections, which were validated by comparison with experimental results from drop tests of ship cross sections. Tveitnes et al. [7] developed a test rig that drives a wedge section, with end plates, down guides to enter the water vertically at near constant velocity to measure the entry force. Wang et al. [8] experimentally and numerically investigated the water entry of a freefall wedge with a focus on the evolution of the pressure on the impact sides and the top side. A series of free-drop model tests were carried out by Wang [9], and the accelerations, slamming pressures and responses were measured. Based on the experimental results of Wang [9], an uncoupled method, combining Wagner theory and the finite element method, was presented to analyze the slamming problem for the 3D structure by Luo [10]. Mohammad et al. [11] addressed the water entry experiments of compliant slender bodies. Wang et al. [12] presented a new Computational Fluid Dynamics (CFD) method to simulate the physical process of a 2D elastic wedge entering water. Aquelet et al. [13] studied the water entry problem based on the Euler-Lagrange coupling algorithm. Yang et al. [14] studied the effect of an elastic structure subjected to water-entry slamming based on the double asymptotic approximation (DAA) method. Zhu et al. [15] studied the fully nonlinear free-surface deformations of initially calm water caused by the water entry and water exit of a horizontal circular cylinder, with both forced and free vertical motions using CIP (constrained interpolation profile) method. Watanabe et al. [16] conducted a model test in both regular and irregular waves to investigate the effects of bow flare on the deck wetness and asymmetry of the vertical wave bending moment. Wang et al. [17] carried out a freefall slamming experiment on a steel model of a hull section. The slamming pressure and structural elastic response stress on both sides of the steel model were measured. Hermundstad and Moan [18] used nonlinear strip theory to calculate the relative motions of ship and wave and found that water pileup around the bow due to the forward speed of the vessel significantly increases the slamming pressures. For the S-175 container ship, a series of model tests and numerical simulations have been carried out worldwide. Xu et al. [19] carried out a model test for an S-175 container ship model in regular waves and compared the experimental results with the results calculated by linear strip theory. Ramos et al. [20] conducted experiments for ship motions and sea loads on a segmented model of the S-175 container ship in ballast conditions and compared the vertical motions and bending moments between the measurements and the theoretical results. Fonseca et al. [21] conducted an experimental investigation of the nonlinear effects on the vertical motions and loads on the S-175 model advancing in irregular waves. Based on the explicit finite element software LS-DYNA, He et al. [22] simulated the water impact of three section models, and the water-entry velocity and the slamming pressures on the bottom were analyzed. A comprehensive review of research in loads on ship have been presented by Hirdaris et al [23]. As the further investigation of slamming, the study of hydroelasticity has been carried out in recent years. For example, Spyridon et al. [24] investigated the sensitivity to whipping and springing responses across a range of headings and sea states by two-dimensional hydroelasticity analysis. Recent advances and future trends of hydroelasticity has been presented by Hirdariset al. [25].
For transport vessels, the loading conditions often change due to different transportation tasks, so it is important to study the loading influence on the structural slamming. However, most studies on water slamming are focused on the material properties of the structure, water inflow rate, bottom angle, etc., and most work on S-175 container ships is focused on the effect of wave load on structure response. Few studies have been conducted on the loading influence on water slamming. Therefore, based on LS-DYNA, the loading influence on free-water slamming was studied. To analyze and reveal the loading influence on water slamming, systematic studies on the three S-175 container ship section models under different loading conditions were carried out. First, a two-dimensional WILS JIP-II hull section model was established. By analyzing and comparing the numerical simulation results and experimental data, the accuracy and reliability of LS-DYNA for simulating water slamming problems were verified. Then, according to the mass relationship between the lightship weight, half-load weight and full-load weight of a container ship [26], section models under different loading conditions were established. Finally, by calculating parameters such as speed and pressure during the water-entry process, the loading influence on water slamming was analyzed systematically.

Arbitrary Lagrangian-Eulerian Algorithm
An explicit finite element method with the ALE (Arbitrary Lagrangian-Eulerian) algorithm was used, which can avoid numerical diffusion and instability. The ALE algorithm is similar to the Euler algorithm; the difference is that the elements in ALE can motion freely. Lagrange time-step is employed, when deformation occurs to the elements due to the motion of the substance. So the ALE time step is executed: (1) To keep the topological relation, the boundary of the substance is kept with the inner elements re-meshed; (2) the element parameters (e.g. density, energy) and the node speed are transformed to the re-meshed elements. Figure 1 shows the difference between the Lagrange, Euler and ALE algorithms. employed, when deformation occurs to the elements due to the motion of the substance. So the ALE time step is executed: (1) To keep the topological relation, the boundary of the substance is kept with the inner elements re-meshed; (2) the element parameters (e.g. density, energy) and the node speed are transformed to the re-meshed elements. Figure 1 shows the difference between the Lagrange, Euler and ALE algorithms. The governing equations can be given by the mass, momentum and energy conservation equations: where t indicates time; ρ indicates density; the relative velocity w is the difference between the material velocity v and the grid velocity u; xi,j represents the Euler coordinates; ij σ is the stress tensor; and E is energy.

Penalty Function Coupling Method
The penalty function coupling method [13] is employed to deal with the fluid-structure-coupling problem in the present study. Figure 2 shows the law of the penalty function coupling method. This coupling method calculates the pressure at each point on the contact surface by the relative displacement d and the derivative of time. Then, the pressure is applied to the nodes on the contact surface, which blocks the fluid from penetrating the structure. The coupling force F is described as: The governing equations can be given by the mass, momentum and energy conservation equations: where t indicates time; ρ indicates density; the relative velocity w is the difference between the material velocity v and the grid velocity u; x i,j represents the Euler coordinates; σ ij is the stress tensor; and E is energy.

Penalty Function Coupling Method
The penalty function coupling method [13] is employed to deal with the fluid-structure-coupling problem in the present study. Figure 2 shows the law of the penalty function coupling method. This coupling method calculates the pressure at each point on the contact surface by the relative displacement d and the derivative of time. Then, the pressure is applied to the nodes on the contact surface, which blocks the fluid from penetrating the structure. The coupling force F is described as: where k is the stiffness coefficient of the penalty function and d is the relative displacement between the slave substance and the master substance. Force F acts on the slave substance and the master substance at the same time in the opposite direction to ensure the balance of forces on the contact surface. F f at the fluid nodes is the fluid force, while F s is the force at the structure coupling node.
where f P is the scaling factor of penalty function stiffness; K indicates the bulk modulus of the fluid; and V is the volume of the fluid unit containing the master material node. A is the average area of the structural elements associated with the slave nodes on the coupling surface. By choosing an appropriate value of f P , the mutual penetration between the master substance and the slave substance can effectively be controlled, while the overall stiffness of the system is not significantly affected.
(a) (b) Figure 2. Sketch of the penalty function coupling algorithm: (a) at time t = t n ; (b) at time t = t n +Δt

Equation of State
A linear polynomial equation of state (EOS) is used for the internal energy. The pressure is given by: where P is pressure; C0 to C6 indicate the polynomial equation coefficients; E is internal energy; μ is given by the following equation: where 0 α and α indicate the initial specific volume and specific volume, respectively. The value of C0 to C6 is given based on Li et al. [27]. The EOS parameters of each material are shown in When the stiffness coefficient of the penalty function k tends to infinity, the relative displacement d tends to zero. Though the contact surface cannot be penetrated, it affects the overall stiffness of the system. The stiffness coefficient of the penalty function is given by the following equation: where P f is the scaling factor of penalty function stiffness; K indicates the bulk modulus of the fluid; and V is the volume of the fluid unit containing the master material node. A is the average area of the structural elements associated with the slave nodes on the coupling surface. By choosing an appropriate value of P f , the mutual penetration between the master substance and the slave substance can effectively be controlled, while the overall stiffness of the system is not significantly affected.

Equation of State
A linear polynomial equation of state (EOS) is used for the internal energy. The pressure is given by: where P is pressure; C 0 to C 6 indicate the polynomial equation coefficients; E is internal energy; µ is given by the following equation: where α 0 and α indicate the initial specific volume and specific volume, respectively. The value of C 0 to C 6 is given based on Li et al. [27]. The EOS parameters of each material are shown in Table 1. E 0 and V 0 indicate the initial internal energy per unit reference volume and initial relative volume, respectively. Table 1. Parameters of the test model.

WILS JIP-II Hull Section Model
To validate the numerical model, the water slamming problem of the Wave Induced Loads on Ships Joint Industry Project II (WILS JIP-II) hull section model was simulated based on the model experiment [28]. The model is free to release at a distance of 0.3 m from the center of gravity to the free surface of the water. Test parameters and pressure gauges of the model are shown in Table 2 and Figure 3, respectively.  Table 1. E0 and V0 indicate the initial internal energy per unit reference volume and initial relative volume, respectively.

WILS JIP-II Hull Section Model
To validate the numerical model, the water slamming problem of the Wave Induced Loads on Ships Joint Industry Project II (WILS JIP-II) hull section model was simulated based on the model experiment [28]. The model is free to release at a distance of 0.3 m from the center of gravity to the free surface of the water. Test parameters and pressure gauges of the model are shown in Table 2 and Figure 3, respectively.

Validation of the Numerical Model
The comparison of the pressure on the gauges between the present results and available experimental data is shown in Figure 5. The solid line is the numerical result, while the dot is the experimental data from Hong et al. [28]. From Figure 5, it can be seen that the numerical results generally agree well with the experimental data, while an overestimated peak of the pressure by the numerical model is also found, especially for P1 gauge. The reason of this overestimated pressure peak may be due to the three-dimensional effect in the experiment. Another reason can be considered as the numerical inaccuracy in capturing the peak of pressure, since the peak of the pressure is sensitive to the numerical schemes, numerical grids, time steps and so forth. Although the peak of the pressure is overestimated, the above comparison generally validates the present LS-DYNA numerical model.

Ship Model and Initial Condition Setting
Three typical cross sections of an S-175 model, including bow, parallel middle body and stern,

Validation of the Numerical Model
The comparison of the pressure on the gauges between the present results and available experimental data is shown in Figure 5. The solid line is the numerical result, while the dot is the experimental data from Hong et al. [28]. From Figure 5, it can be seen that the numerical results generally agree well with the experimental data, while an overestimated peak of the pressure by the numerical model is also found, especially for P1 gauge. The reason of this overestimated pressure peak may be due to the three-dimensional effect in the experiment. Another reason can be considered as the numerical inaccuracy in capturing the peak of pressure, since the peak of the pressure is sensitive to the numerical schemes, numerical grids, time steps and so forth. Although the peak of the pressure is overestimated, the above comparison generally validates the present LS-DYNA numerical model.

Validation of the Numerical Model
The comparison of the pressure on the gauges between the present results and available experimental data is shown in Figure 5. The solid line is the numerical result, while the dot is the experimental data from Hong et al. [28]. From Figure 5, it can be seen that the numerical results generally agree well with the experimental data, while an overestimated peak of the pressure by the numerical model is also found, especially for P1 gauge. The reason of this overestimated pressure peak may be due to the three-dimensional effect in the experiment. Another reason can be considered as the numerical inaccuracy in capturing the peak of pressure, since the peak of the pressure is sensitive to the numerical schemes, numerical grids, time steps and so forth. Although the peak of the pressure is overestimated, the above comparison generally validates the present LS-DYNA numerical model.

Ship Model and Initial Condition Setting
Three typical cross sections of an S-175 model, including bow, parallel middle body and stern,

Ship Model and Initial Condition Setting
Three typical cross sections of an S-175 model, including bow, parallel middle body and stern, are selected, which correspond to station numbers 19, 10 and 1, respectively. The three sectional models are established and simulated in full-scale size. Tables 3 and 4, respectively, list the main dimensions of the S-175 container ship and the water-entry velocity of each section. The velocity is given by He et al. [29]. The parameters of materials are shown in Table 5. The sectional models are set to be rigid parts to avoid the influence of elastic response. Figure 6 shows the body plan of S-175, and Figure 7 shows the cross section of S-175. The positions of the pressure gauges P1 (x = 1.2 m), P2 (x = 3 m) and P3 (x = 6 m) are shown in Figure 7. Transverse center of gravity (m) 0 (from midcourt line) Vertical center of gravity (m) 8.5 (from keel strake) Table 4. Velocities of water entry for each section.

Model
Velocity (m/s) are selected, which correspond to station numbers 19, 10 and 1, respectively. The three sectional models are established and simulated in full-scale size. Tables 3 and 4, respectively, list the main dimensions of the S-175 container ship and the water-entry velocity of each section. The velocity is given by He et al. [29]. The parameters of materials are shown in Table 5. The sectional models are set to be rigid parts to avoid the influence of elastic response. Figure 6 shows the body plan of S-175, and Figure 7 shows the cross section of S-175. The positions of the pressure gauges P1 (x = 1.2 m), P2 (x = 3 m) and P3 (x = 6 m) are shown in Figure 7.    The grids for the bow section are shown in Figure 8.  Table 6.  Benefiting from the symmetry of the hull structure and the simulation problem, the present study uses a half computational domain. The left boundary of the calculation domain (x = 0) is set to a symmetric boundary. Other boundaries are set to nonreflective boundaries. The water depth, height above the water (air area) and half width of the computational domain are set to 30, 40 and 100 m, respectively. A nonuniform grid is used except for the grids near the water surface, which are evenly arranged. Thin shell elements are used for the ship structure, while the fluid domain is established by 1 point ALE multi-material elements [30]. The minimum size of the mesh is 0.05 m. The grids for the bow section are shown in Figure 8. The grids for the bow section are shown in Figure 8.  Table 6. Table 6. Loading set of each section.

Section
Loading case Mass (kg) Bow Lightship weight 2809 The loading in the ship can affect the response of ship motion, correspondingly affecting the heave velocity (water-entry velocity) and slamming pressure. The mass of the bow, parallel middle body and stern sections at no-load (lightship weight), half-load and full-load conditions are shown in Table 6.  Figure 9 shows the time history of the velocity of the bow section under different loading conditions. The change of the velocity of the section under the three loading conditions shows a similar trend [17]. Because the mass of the section is light under no-load (lightship weight) conditions, the gravity of the bow is quickly balanced with the fluid force it receives from the water. It can be confirmed from Figure 9 that the turning point of the velocity of the section under the no-load (lightship weight) condition is earlier than for the half-load and full-load conditions. The time instant of the turning point in the velocity curve under full-load conditions is obviously earlier than for the no-load (lightship weight) condition. This is because the wave splashes from the bottom of the bow reaches the side of the bow and creates a slamming phenomenon. Figure 10 shows the scenarios of the flare slamming under no-load (lightship weight) and full-load conditions for a 14.6 m falling height of the bow. As shown in Figure 10, the larger the loading is, the more intense the splashing waves are, which means a greater slamming force is induced.  m falling height of the bow. As shown in Figure 10, the larger the loading is, the more intense the splashing waves are, which means a greater slamming force is induced.   m falling height of the bow. As shown in Figure 10, the larger the loading is, the more intense the splashing waves are, which means a greater slamming force is induced. The slamming of the bow is mainly divided into two phenomena: flare slamming and slamming when entering the water. Therefore, the pressure gauge P1 is set at the head of the bow, and P2 and P3 are set at the side surface of the bow. The slamming forces of the bow section under the three loading conditions are shown in Figure 11, and Figures 12-14 show the pressure at the three gauges under no-load (lightship weight), half-load and full-load conditions, respectively. The slamming of the bow is mainly divided into two phenomena: flare slamming and slamming when entering the water. Therefore, the pressure gauge P1 is set at the head of the bow, and P2 and P3 are set at the side surface of the bow. The slamming forces of the bow section under the three loading conditions are shown in Figure 11, and Figures 12-14 show the pressure at the three gauges under no-load (lightship weight), half-load and full-load conditions, respectively.   The slamming of the bow is mainly divided into two phenomena: flare slamming and slamming when entering the water. Therefore, the pressure gauge P1 is set at the head of the bow, and P2 and P3 are set at the side surface of the bow. The slamming forces of the bow section under the three loading conditions are shown in Figure 11, and Figures 12-14 show the pressure at the three gauges under no-load (lightship weight), half-load and full-load conditions, respectively.   The slamming of the bow is mainly divided into two phenomena: flare slamming and slamming when entering the water. Therefore, the pressure gauge P1 is set at the head of the bow, and P2 and P3 are set at the side surface of the bow. The slamming forces of the bow section under the three loading conditions are shown in Figure 11, and Figures 12-14 show the pressure at the three gauges under no-load (lightship weight), half-load and full-load conditions, respectively.     As shown in Figure 11, three stages in the slamming force of the bow are clearly observed even under different loading conditions. In the first stage, the bow begins to fall into the water, and the bow section is subjected to a small slamming pressure. In the second stage, flare slamming occurs at the side of the bow, which causes a large slamming pressure. In the third stage, the cavity between the hull and fluid disappears, and the backward flow around the hull side causes another high peak impact. From Figure 11, it can be seen that the slamming force of the bow section increases as the loading increases. In particular, the flare slamming force at the second stage mentioned above increased greatly, which confirms the idea that a larger loading induces a greater flare slamming force on the bow section.

Bow Section
It can be seen from Figures 12(b), 13(b), and 14(b) that two large pressure peaks occur at P2. One peak occurs at the moment when the flare slamming occurs, and the other peak occurs at the moment when the cavity disappears. The increase in the loading leads to a sharp increase in the first peak; the influence on the second peak is smaller. In other words, the loading has a strong influence on the flare slamming of the bow. As shown in Figure 11, three stages in the slamming force of the bow are clearly observed even under different loading conditions. In the first stage, the bow begins to fall into the water, and the bow section is subjected to a small slamming pressure. In the second stage, flare slamming occurs at the side of the bow, which causes a large slamming pressure. In the third stage, the cavity between the hull and fluid disappears, and the backward flow around the hull side causes another high peak impact. From Figure 11, it can be seen that the slamming force of the bow section increases as the loading increases. In particular, the flare slamming force at the second stage mentioned above increased greatly, which confirms the idea that a larger loading induces a greater flare slamming force on the bow section.

Parallel Middle-Body Section
It can be seen from Figure 12, Figure 13, and Figure 14 that two large pressure peaks occur at P2. One peak occurs at the moment when the flare slamming occurs, and the other peak occurs at the moment when the cavity disappears. The increase in the loading leads to a sharp increase in the first peak; the influence on the second peak is smaller. In other words, the loading has a strong influence on the flare slamming of the bow. The section shape of the parallel middle body is similar to a rectangle shape, which suffers mainly from bottom slamming. The slamming forces of the parallel middle body section under the three loading conditions are shown in Figure 16. Figures 17-19 show the pressure at each gauge under no-load (lightship weight), half-load and full-load conditions, respectively. The section shape of the parallel middle body is similar to a rectangle shape, which suffers mainly from bottom slamming. The slamming forces of the parallel middle body section under the three loading conditions are shown in Figure 16. Figures 17-19 show the pressure at each gauge under no-load (lightship weight), half-load and full-load conditions, respectively. The section shape of the parallel middle body is similar to a rectangle shape, which suffers mainly from bottom slamming. The slamming forces of the parallel middle body section under the three loading conditions are shown in Figure 16. Figures 17-19 show the pressure at each gauge under no-load (lightship weight), half-load and full-load conditions, respectively.  The section shape of the parallel middle body is similar to a rectangle shape, which suffers mainly from bottom slamming. The slamming forces of the parallel middle body section under the three loading conditions are shown in Figure 16. Figures 17-19 show the pressure at each gauge under no-load (lightship weight), half-load and full-load conditions, respectively. It can be seen from Figure 16 that the slamming force of the parallel middle body occurs mainly at the moment when it touches the water, and the influence of the loading on the slamming force is small. However, after the process of water entry becomes stable, the force on the parallel middle It can be seen from Figure 16 that the slamming force of the parallel middle body occurs mainly at the moment when it touches the water, and the influence of the loading on the slamming force is small. However, after the process of water entry becomes stable, the force on the parallel middle It can be seen from Figure 16 that the slamming force of the parallel middle body occurs mainly at the moment when it touches the water, and the influence of the loading on the slamming force is small. However, after the process of water entry becomes stable, the force on the parallel middle body gradually increases, and the larger the loading is, the greater the increase is. According to Wagner's theory, the slamming pressure at any position of a wedge can be expressed by:

Parallel Middle-Body Section
where ρ indicates the density of the coupled fluid; K w is the dimensionless pressure coefficient, which is determined by the deadrise angle of the wedge; and V n is the impact velocity of the wedge. At the stage where the water entry becomes stable, we can see from Figure 15 that, as the loading increases, the velocity of the parallel middle body becomes significantly larger, so the phenomenon mentioned above occurs.
By comparing the pressures at each gauge in Figures 17-19, it can be found that the increase in loading has little influence on P1 and P2, which are placed near the centerline of ship section, while the pressure at P3, near the side of the ship, is strongly affected. This is because the change in loading has little influence on the velocity of the section at the initial stage during water entry. According to Equation (6), the slamming pressures at P1 and P2 at the moment of contact with water are not sensitive to the different loading conditions. However, when the P3 gauge contacts the water, the difference in velocity of the section among the three loading conditions is larger. Therefore, the pressure under the full load condition is much greater than that under the no-load (lightship weight) condition.

Stern Section
The shape of the stern section is similar to that of the bow section. However, since the body-plan of the stern section is smoother than that of the bow, the water-entry process of the stern is much gentler than that of the bow. Figure 20  The shape of the stern section is similar to that of the bow section. However, since the body-plan of the stern section is smoother than that of the bow, the water-entry process of the stern is much gentler than that of the bow. Figure 20 shows  As shown in Figure 21, no significant change occurs in the velocity of the body at the moment when the bottom of the stern starts to enter the water. Because there is no flare slamming, the order of velocity turning points is no load (lightship weight), half load and full load. Compared with the case of the parallel middle body, the difference in velocity of the stern under the three loading conditions is smaller after the water-entry process is stable.  As shown in Figure 21, no significant change occurs in the velocity of the body at the moment when the bottom of the stern starts to enter the water. Because there is no flare slamming, the order of velocity turning points is no load (lightship weight), half load and full load. Compared with the case of the parallel middle body, the difference in velocity of the stern under the three loading conditions is smaller after the water-entry process is stable. Figure 21, no significant change occurs in the velocity of the body at the moment when the bottom of the stern starts to enter the water. Because there is no flare slamming, the order of velocity turning points is no load (lightship weight), half load and full load. Compared with the case of the parallel middle body, the difference in velocity of the stern under the three loading conditions is smaller after the water-entry process is stable. It can be seen from Figure 22 that since there is no flare slamming, the influence of the loading on the stern's slamming force is relatively small. Similar to the case of the parallel middle body, when the process of water entry becomes stable, the pressure on the stern gradually increases, and the larger the loading is, the greater the increase is. It can be seen from Figure 22 that since there is no flare slamming, the influence of the loading on the stern's slamming force is relatively small. Similar to the case of the parallel middle body, when the process of water entry becomes stable, the pressure on the stern gradually increases, and the larger the loading is, the greater the increase is. Comparing the pressures at each gauge in Figures 23-25, it can be seen that the peak value at each point increases significantly with increasing loading; as the loading becomes larger, the pressure after the peak is also larger. This is consistent with the trend of the slamming force.  Comparing the pressures at each gauge in Figures 23-25, it can be seen that the peak value at each point increases significantly with increasing loading; as the loading becomes larger, the pressure after the peak is also larger. This is consistent with the trend of the slamming force. Comparing the pressures at each gauge in Figures 23-25, it can be seen that the peak value at each point increases significantly with increasing loading; as the loading becomes larger, the pressure after the peak is also larger. This is consistent with the trend of the slamming force.

Loading Sensitivity Analysis
The sectional shape has a great influence on the velocity, force and pressure of the sections. The

Loading Sensitivity Analysis
The sectional shape has a great influence on the velocity, force and pressure of the sections. The sensitivity of different sections to loading conditions is also different. In this section, the data of the three sections under the three loading conditions are compared and analyzed to investigate the loading sensitivity on the water entry of different sections. Figure 26 shows the velocity of each section under each loading condition. It can be seen from Figure 26a that the velocity during water entry of the parallel middle body changes less than for the cases of the bow and the stern. By comparing Figure 26a-c, it can be found that as the loading increases, the peak value of the velocity of each section increases, and the velocity of the parallel middle body is most sensitive to the loading.

Loading Sensitivity Analysis
The sectional shape has a great influence on the velocity, force and pressure of the sections. The    Table 7. By comparing Figure 27a-c and Table 7, it can be found that the peak of the slamming force of the parallel middle body is not very sensitive to the loading, and the first peak of the bow is more affected by the loading due to the flare slamming.
After the process of water entry becomes stable, the slamming forces of the three sections increase with increasing loading. The parallel middle body is most affected by the loading.   Table 7. By comparing Figure 27a-c and Table 7, it can be found that the peak of the slamming force of the parallel middle body is not very sensitive to the loading, and the first peak of the bow is more affected by the loading due to the flare slamming. After the process of water entry becomes stable, the slamming forces of the three sections increase with increasing loading. The parallel middle body is most affected by the loading.   Table 7. By comparing Figure 27a-c and Table 7, it can be found that the peak of the slamming force of the parallel middle body is not very sensitive to the loading, and the first peak of the bow is more affected by the loading due to the flare slamming.
After the process of water entry becomes stable, the slamming forces of the three sections increase with increasing loading. The parallel middle body is most affected by the loading.   Figure 28 shows the variation of the slamming pressure at P2 under different loading conditions, and the peak value and their growth rates are shown in Table 8. Under each loading condition, the peak pressure of the parallel middle body is much larger than that of the bow and the stern, but the loading has little influence on the peak value. The peak of the slamming pressure of the bow is greatly influenced by the loading. The difference in the pressure peak between the no-load (lightship weight) and full-load conditions is nearly a factor of two. Therefore, in the process of water entry, the slamming force of the bow is most influenced by the loading due to flare slamming; because of the large difference in the velocity under different loading conditions, the slamming force of the parallel middle body is influenced by the loading more than that of the other two sections when the water entry becomes stable.   Figure 28 shows the variation of the slamming pressure at P2 under different loading conditions, and the peak value and their growth rates are shown in Table 8. Under each loading condition, the peak pressure of the parallel middle body is much larger than that of the bow and the stern, but the loading has little influence on the peak value. The peak of the slamming pressure of the bow is greatly influenced by the loading. The difference in the pressure peak between the no-load (lightship weight) and full-load conditions is nearly a factor of two. Therefore, in the process of water entry, the slamming force of the bow is most influenced by the loading due to flare slamming; because of the large difference in the velocity under different loading conditions, the slamming force of the parallel middle body is influenced by the loading more than that of the other two sections when the water entry becomes stable.

Conclusion
The influence of the loading and shape of sections on the slamming pressure was studied based

Conclusions
The influence of the loading and shape of sections on the slamming pressure was studied based on LS-DYNA. The process of the water entry of three typical sections, including the bow, parallel middle body and stern, under three loading conditions (no load (lightship weight), half load and full load) was simulated and discussed. By analyzing the velocity, force and pressure of the S-175 container ship sections under each loading condition, it can be concluded that: (1) For the bow section, the first peak of slamming force induced by the flare slamming is greatly influenced by the loading condition.
(2) For the parallel middle body, its slamming force is much larger than that of the bow and the stern section, but the peak is less sensitive to the loading. However, the loading has a greater influence on the slamming force after the water entry becomes stable. (3) For the stern section, it has a smaller slamming force than that of the bow and the parallel middle body sections, and the force and pressure change smoothly in the water-entry process. According to the above investigation, the ship slamming force is sensitive to the loading and the section shape, so both of them need to be considered carefully in the early stage of ship design.