#### Abstract

A fully coupled nonlinear three-dimensional (3D) hydroelastic method is developed to investigate vibrational responses of a large ship with a pronounced bow flare subjected to high seas. This numerical model consists of a 3D boundary element method, 1D Euler-Bernoulli beam model, and a 2D generalized Wagner model. Green water loads were considered. Experimental study was carried out in a towing tank on a self-propelled segmented model with nonuniform steel backbones. The ship model was tested in regular incident waves of large amplitude. Impact pressure and nonlinear vertical bending moments were measured and compared with numerical predictions. The proposed nonlinear model produced similar results to the experimental model. Furthermore, the effects of elastic modes and nonlinearities on the numerical results were analyzed.

#### 1. Introduction

Detailed knowledge of wave-induced structural vibration of marine vessels sailing in seaways, more often than not in adverse weather, is of pivotal significance for the evaluation of seakeeping ability. Large vessels, their lowest eigenfrequencies of elastic vibration lying in the vicinity of the harmonics of ocean waves, are highly susceptible to severe transient vibrational responses (whipping) caused by wave impact on either the bow or stern region. This occurs when pitching motion and forward speed lead to emergence of the ship hull above the water surface and subsequent penetration into the wave surface. A pronounced bow flare will further aggravate the severity of whipping due to the large bow-wave contact area. Since the structural vibration of large vessels significantly affects the surrounding fluid flow and therefore the hydrodynamic forces, hydroelastic analysis should be used to analyze ship vibratory behaviors in waves.

Beginning with the pioneering work by Bishop and Price [1], significant developments have been achieved in theoretical and experimental methods of predicting the hydroelastic loads acting on ships. Wu et al. [2] extended the 2D linear hydroelastic theory to analyze the vibrational responses of slender ships in waves by coupling Euler-Bernoulli beam theory with linear 2D strip theory. Hermundstad et al. [3] proposed a 2D linear hydroelastic method to study the vibrational responses of high speed ships. Price and Wu [4] presented a 3D linear hydroelastic method to study structural responses of a small waterplane area twin hull (SWATH) in which 3D effects near the bow and stern cannot be omitted.

Recently the occurrence of extreme sea conditions is more frequently reported than before, thus highlighting the importance of nonlinear hydroelastic theories. Mikami and Kashiwagi [5] proposed a 2D nonlinear hydroelastic theory to study the characteristics of ship hull vibration. The calculations of incident wave forces (Froude-Krylov forces) and hydrostatic forces are performed for the instantaneous wet hull. In order to further analyze whipping responses, Rajendran et al. [6] presented a comparison of 2D nonlinear computational results and experimental results. The bow flare slamming loads were calculated based on von Kármán model.

Y. Kim and J.-H. Kim [7] conducted a benchmark study on motions and loads of a 6750-TEU containership to compare different numerical methods ranging from 2D theories to 3D theories. On average, good agreement was found between computational and experimental results with respect to linear responses. However, for larger wave amplitudes, the computational results were more scattered. Tuitman [8] developed a practical method to calculate the response of the ship structure due to the seakeeping loading including slamming loads. Kim et al. [9] presented a 3D nonlinear hydroelastic method to study whipping responses by coupling a 3D finite element method (FEM) with a 3D Rankine panel based boundary element method (BEM). The nonlinear Froude-Krylov forces and hydrostatic forces were computed from the instantaneous wetted surface. The bow flare impact forces were computed based on the 2D generalized Wagner method. The 2D generalized Wagner method is much more efficient than 3D CFD models and is well tested by drop tests of ship cross sections [10]. However, more experimental validation is needed for the whole elastic vessel model in waves.

Experimental studies play a vital role not only in validating applicability of numerical methods but also in discovering new phenomena [11]. Zhu and Moan [12] conducted an experimental study on whipping responses of two large containerships and the experimental results were compared with empirical formulas. The discrepancy found implied that more data from real ships was needed to further improve the accuracy of empirical methods. Jialong et al. [13] presented an experimental and numerical research of hull girder vibration of a large ship. The impact wave forces were calculated based on momentum impact theory. The instantaneous wetted surface was approximately obtained by intersecting the original hydrodynamic mesh with the still water surface. Jialong et al. [14] conducted a large-scale model test in real coastal sea conditions. The bow flare slamming and the subsequent whipping phenomena were captured and analyzed.

The gradual march towards a complete understanding of fluid-structure interaction is ongoing. Although various formulations have been developed, some open questions remain on the nonlinear slam-induced responses of large vessels in severe seas, and the available experimental data are not very large. In the present research, a nonlinear time domain hydroelastic method was used to simulate whipping of a large ship in high seas. A 3D BEM, 2D generalized Wagner method, 1D beam model, and 1D dam breaking model were fully coupled in this numerical model. Experimental test of a large ship in severe seas was carried out to validate the numerical method. The numerically and experimentally acquired slamming pressures were compared. A comprehensive comparison and analysis of the measured and simulated structural response of the vessel were provided. In addition, the effects of nonlinearity and elasticity on the response were analyzed.

#### 2. Theory

Figure 1 demonstrates the general formulation of the problem of wave-induced structural vibration of vessels sailing in regular incident waves. The vessel initially travels at a constant velocity . The fluid domain is assumed to be infinite. Two inertial frames of reference are introduced to analyze this hydroelastic problem: the earth fixed Cartesian frame with its coordinate plane positioned at the calm waterplane and the translational Cartesian frame which translates invariantly at the velocity and remains parallel to . Then is appropriate to describe motion variations of the vessel. At the initial instant, the two frames are identical in space. The Galilean transformation relation between these two frames at time is shown as follows:

In (1), is the heading angle and in head seas is defined to be 0. In Figure 1, is the infinite fluid domain, is the infinite radius of , is the infinite depth of , is the boundary of , and is the free surface. denotes the structural domain and denotes the contact surface between and .

Figure 2 presents three representative states of the vessel: initial state, static hydroelastic state, and the transient hydroelastic state. The initial state corresponds to the initially undeformed vessel in calm water, which serves as the input for the computation of the following static hydroelastic state. The static hydroelastic state represents the deformed configuration of the vessel in calm water due to the uneven distribution of mass. The static hydroelastic state is the ship configuration based on which some hydrodynamic forces are computed. The transient hydroelastic state is the instantaneous deformed configuration of the vessel in waves. Some hydrodynamic forces are calculated based on the transient hydroelastic state, and these are explained later in this section. Figure 2(c) shows the transient hydroelastic state in severe sea conditions, where slamming events and green water phenomena occur.

**(a)**

**(b)**

**(c)**

##### 2.1. Structural Domain

The vessel structure as a whole is assumed to be homogeneous, isotropic, and elastic. Instead of adopting complicated and time-consuming 3D FE model, the elastic hull is modeled as a nonuniform Euler-Bernoulli beam. The beam model is more suitable, especially when 3D FE model is not available. In the inertia translational reference frame , the equation of motion of the elastic vessel can be written as follows [15]:where , , and are mass, damping, and stiffness matrices and is the external force matrix. , , and are the nodal displacement, velocity, and acceleration matrices, respectively.

Based on the principle of superposition, nodal variables , , and can be expressed as a sum of the distortion in the principal modes of the vessel:where denotes the principal mode matrix and denotes the principal coordinate matrix.

Substituting (3) into (2), multiplying by , the equation of motion expressed in terms of principal coordinates is obtained:where , , , and . , , , and are called generalized mass, damping, stiffness, and external force matrices, respectively.

##### 2.2. Fluid Domain

###### 2.2.1. 3D Components of Wave Forces

By assuming that the fluid flow is nonviscid, incompressible, and irrotational, the potential flow theory can be employed to model the fluid in the problem of fluid-structure interaction. In the reference frame , the fluid velocity potential of the fluid-structure system in steady state can be written as [16]where is the spatial part of , is the imaginary unit, and is the wave encounter frequency. is defined aswhere is the natural wave circular frequency and is the wave number.

can be further decomposed into three parts:where denotes the incident wave potential, denotes the diffraction wave potential, and denotes the radiation wave potential of the th degree of freedom of the vessel. is greater than 6 and indices represent the elastic degrees of freedom. , , and can be solved based on the boundary element method [16].

The generalized diffraction forces and radiation forces are calculated based on the wet contact surface of the fluid-vessel system in the static hydroelastic state as illustrated in Figure 2:where is the fluid density, is the incident wave amplitude, is the normal vector pointing into the fluid domain, is the th principal mode of the vessel, and is the infinitesimal area element of the wet contact surface.

The generalized Froude-Krylov forces and hydrostatic restoring forces are computed based on the wet contact surface of the fluid-vessel system in the transient hydroelastic state as demonstrated in Figure 2:where is the vertical displacement of the th degree of freedom and is the mass density of the vessel.

###### 2.2.2. Impact Wave Forces

Slamming loads acting on deformable large vessels with considerable bow flares are calculated based on the 2D generalized Wagner method [17]. Some selected cross sections are chosen to calculate the 2D slamming pressure [10]. Figure 3 displays the description of 2D water entry of cross sections. is the velocity of the section, and are the coordinates of the contact point at time , and is the contour of the cross section. In the present study, is the combination of rigid and elastic body motions of the vessel. Note that the effects of pitch on the motion of those cross sections are confined to the vertical components because the ship is slender.

The 2D potential flow theory is used to model this fluid-structure problem and the slamming pressure acting on the wet contour of each cross section can be solved by the conformal mapping method. Figure 4 demonstrates how conformal mapping is performed. The real free surface contour is linearized as the horizontal line starting from the contact point. The line segments AB and CD from the original real plane are mapped to intervals and in the complex plane , respectively. The wet contour BC is mapped to the interval . The slamming pressure can be obtained as follows:where is the vertical displacement of the cross section. The expressions of and can be found in [17]; here we do not explain these two expressions for brevity.

**(a)**

**(b)**

The solution in the original plane can be obtained based on the inverse of the conformal mapping. Then the generalized slamming forces can be calculated aswhere is the normal vector pointing out of the wet contour of the cross section and is the infinitesimal line segment of the wet contact line of the section.

###### 2.2.3. Green Water Loads

In extreme sea conditions, after the bow impact events, the incident wave may overtop the bow and run up onto the deck. In this paper, green water loads on deck are incorporated and simulated based on the dam-breaking Model. This method models the green water phenomenon as the process of sudden collapse of a dam. Before the sudden collapse, the water stored in the dam is static. Further we assume that during the process the ground exerts no resistant forces to the fluid. We introduce a 2D bow-fixed reference frame in which the origin is located at the tip of the bow. The simplified Saint-Venant equations [18] can be expressed aswhere and are the depth and velocity of the green water flow, respectively. is the gravitational acceleration. Based on the method of characteristics, the depth and velocity can be solved as follows:where is the initial depth of the water in the intact dam. Velocity is a linear function of , and green water pressure on deck can be easily calculated from . We assume that does not change in direction; then the generalized green water loads can be calculated as

##### 2.3. Hydroelastic Vibration Domain

Based on the analyses of the structural and fluid domain, now the hydroelastic equation of motion can be rewritten as follows:where the right-hand side (RHS) of (4) is decomposed into the six components shown in the RHS of (16). , , and are the generalized incident, diffraction, and radiation wave force matrices, respectively, is the generalized hydrostatic restoring force matrix, and and are the generalized slamming and green water force matrices, respectively.

For large elastic vessels with pronounced bow flare in severe sea conditions, it is of great importance to appropriately incorporate some nonlinearities in those generalized hydrodynamic force matrices of (16). This may not be purely mathematically consistent because (16) is basically linear. However, this treatment does allow evident improvements to the simulation of fierce fluid-structure interaction. As formulated in Section 2.2, among the six generalized hydrodynamic force matrices, only and are calculated based on the unchanged static hydroelastic state of the vessel. This is because of the difficulty of recalculating and using 3D BEM based on the updated wet configuration of the vessel. In contrast, the other four generalized forces are calculated based on the actual instantaneous wet contact surface because the computational burden is acceptable. In other words, except for and , all the other generalized external forces are nonlinear.

##### 2.4. Necessary Numerical Details

###### 2.4.1. Mesh Generation

Mesh generation is the prerequisite of the computation of the generalized external forces. The hull mesh for the initial state of the vessel is generated by employing the method of cubic spline function [19]. Then the transient elastic hull mesh can be generated according to the superposition of the distortion in the principal modes. Figure 5 displays the initial hull mesh and the transient elastic hull mesh.

**(a)**

**(b)**

###### 2.4.2. Nonlinear and

The nonlinear generalized force matrices and are calculated based on the instantaneous wet contact surface and therefore the acquisition of the exact wet hull plays a significant role. The method used in the present study is called the coarse-fine mesh strategy (see Figure 6), which is illustrated as follows: the initial hull mesh is classified as the coarse mesh. The coarse mesh panels are checked at every time step for their position relative to the incident wave profile. Those panels whose vertices are all below (or above) the incident wave contour are considered as the ordinary fully wet (or dry) panels. Special care is given for those partially wet ones. Each partial wet panel is further subdivided into smaller subpanels which are called the fine mesh. Then the reevaluation of the relative position of the fine mesh with respect to the incident wave profile is carried out. The incident wave profile within every fine mesh panel is approximated as the straight line. Based on this strategy, and are calculated by the formulations given in Section 2.2.

###### 2.4.3. Nonlinear

The extraction of the cross sections of the vessel can be carried out based on the information of the hull geometry as shown in Figure 7. For bow and stern areas, more cross sections should be extracted not only due to the fact that the slamming phenomena are more notable in these two areas but also for the compensation of the evident 3D effects associated with these regions.

As mentioned above, the generalized Wagner method assumes that there is no flow separation in the slamming phenomena. Therefore modifications should be made for those sections whose breadth of the outer contour does not increase monotonically with respect to the height of the contour. Figure 8 takes one of the cross sections of the bulb bow as an example to demonstrate those modifications.

For the cross sections where the outer contour consists of two disjoint curves as shown in Figure 8, we use the smooth curve called the connector to connect the two separated modified curves. The wet contour of the connector is excluded in the computation of .

###### 2.4.4. Nonlinear

As formulated in Section 2.2.3, the green water phenomenon is assumed to be one-dimensional. Figure 9 shows green water depth distribution at time along the deck. The mesh for the deck is generated as shown in Figure 5. in (14) is determined by the relative position of the deck and the incident wave profile at time . Then is computed on the wet deck mesh based on (15).

###### 2.4.5. Hydroelastic Response

In the present study, we adopt the well tested explicit fourth-order Runge-Kutta algorithm [20] to solve (16). Then the hydroelastic response can be acquired based on the principle of modal superposition:where , , and are the displacement, vertical bending moment, and shear force; , , and are the th mode of displacement, bending moment, and shear force, respectively.

#### 3. Experimental Setup

To evaluate the capabilities of the methodology presented in this study, a model test of a large vessel with a pronounced bow flare was carried out in the ship model towing tank of Harbin Engineering University. The model was provided by a commercial company (China State Shipbuilding Corporation). The towing tank, as shown in Figure 10, is 108 m in length, 7 m in width, and 3.5 m in depth, with the wave absorption beach located at the front end and the electrically driven flap-type wave generator at the back end. The wave generator enables the generation of long-crest regular waves, the periods ranging from 0.4 s to 4 s, and irregular waves with the maximum significant wave height being 0.32 m. A wave probe is positioned close to the wave maker. A motion recorder apparatus, which is both installed in the model and attached to the towing carriage, is designed to obtain the time history of heave and pitch and fix the yaw angle. The maximum velocity of the towing carriage is 6 m/s. The schematic of the experimental setup is shown in Figure 11.

**(a)**

**(b)**

The elastically scaled model consists of seven segmented hulls made from epoxy, three rectangular nonuniform steel beams connecting those segments, a self-propelled system positioned in the stern, and the wired monitoring system used to acquire experimental data. Table 1 gives the principal particulars of the elastic model. The propelling system was designed to maintain the expected forward speed of the model in still water. Besides the aforementioned motion recorder, the tethered monitoring units in the model comprise (i) strain gauges glued on the surface of the beams and (ii) impact pressure sensors positioned in the bow area. The main particulars of those units are demonstrated in Table 2. Figure 12 presents necessary details of the arrangement of the beams and pressure sensors. Some complementary and schematic information of experimental setup is displayed in Figure 13.

**(a)**

**(b)**

#### 4. Comparison between Numerical and Experimental Results

The numerical model was established based on the experimental setup. For the nonuniform beam model, 20 elements were utilized and only the first three elastic modes of vertical bending were considered to model the elastic characteristics. The whole ship surface was divided into 2074 panels, 56 of which were the deck mesh for the computation of . One coarse panel was subdivided into 32 subpanels when necessary for the computation of and . 22 unevenly spaced cross sections were chosen to calculate . A series of hydroelastic experiments in regular waves were conducted. The present work focuses on those large wave amplitude conditions. Two representative experimental conditions were picked out for the following comparison and analysis (see Table 3). It should be mentioned here that, for the preservation of commercial benefits of the company that provided the vessel prototype, some parts of the pictorial details in the following experimental photographs were pixelated.

##### 4.1. Natural Vibrational Characteristics

The elastic modes were calculated based on an improved Euler-Bernoulli beam model. In this treatment, the shell of the vessel was incorporated and modeled as rigid elements assigned with mass. Then the effects of mass distribution and subsequent moment of inertia can be considered in the computation of the natural vibrational characteristics. In order to validate the result, a hammering test, which involves applying an arbitrary impact to the ship model in still water, was conducted before the model tests in waves. The hammering tests were repeated three times at three locations of the model and the averaged wet natural frequencies were compared with those acquired numerically (see Table 4). A satisfactory agreement appeared in the results of the first order, with increasingly larger discrepancies in higher orders. This was acceptable because the elastic backbones were designed based on the first-order vibrational characteristics. Moreover, the hydroelastic response components associated with the second and third orders were fairly small compared with those corresponding to the first order. The first three elastic natural modes are shown in Figure 14.

**(a)**

**(b)**

##### 4.2. Wave Calibration

Since the parameters of experimentally generated regular waves may slightly deviate from the designed ones, in the present study, a wave probe was installed near the wave maker to measure the uncertainties associated with the parameters of the generated waves. For the following comparison and analysis, numerical calculation will be based on the measured data from the wave probe to reduce as many preventable errors as possible.

Figure 15 shows the time history of the experimentally generated waves whose wave height and frequency are 181.25 mm and 1.66 s. The whole process can be divided into two regions: (i) region I describes the duration before the wave front reaches the wave probe and records the first five periods of the waves with ever-increasing wave heights; (ii) region II illustrates the relatively smooth periods of the waves. Among all the periods, 17 periods from 40 s to 70 s as displayed in Figures 15 and 16 are picked out, since those periods correspond to the effective experimental results. The fluctuations of the wave parameters are depicted in Figure 16, thus suggesting strongly the necessity of modifications of the corresponding numerical inputs in the following computation of slamming pressure and bending moments. In this paper, we replace the original wave parameters with the averaged value obtained from the data provided by the wave probe to perform numerical calculation. For brevity, the pictorial details of wave parameters of the second experimental condition are not presented here. The modified wave parameters of the two representative experimental conditions are presented in Table 5.

**(a)**

**(b)**

##### 4.3. Slamming Loads in Bow Area

As shown in Figure 12, 12 pressure sensors were positioned on the bow area to measure the impact pressure. Pressure sensor number 2 was picked out as an example to compare experimental and numerical results. Figures 17 and 18, as well as all the other experimental results in this study, reveal the inevitability of the presence of noise in the experimental pressure recordings. This is due to the effects of the surrounding metallic structures and other electromagnetic sources on the pressure sensor recording system. Furthermore, it can be deduced from Figures 17 and 18 that, for those experimental results in which the duration of the impact is comparable to the period of the noise, the common method of filtering the original experimental data to obtain the realistic impact peak value fails to be applicable. Therefore, the exact realistic peak impact value should be replaced by an estimated one, with the error being the peak noise value that is extracted by filtering operation (see the enlarged images in Figures 17 and 18). Table 6 gives the impact peak value intervals. Although the exact peak value cannot be obtained, other information, such as the local structure vibration near the sensors, can be deduced from the pressure oscillations through filtering operation as also shown in those enlarged images.

**(a)**

**(b)**

For both operational conditions (see Figure 20), the experimentally obtained impact peak values and periods vary from each other, resulting from the fluctuations of the parameters of experimentally generated waves and the variations of the speed of the test model. On average, the numerically obtained impact peak values are larger than those acquired experimentally. In Figure 20(b), the experimental result shows that, for the severe sea conditions, as exemplified in Figure 19, the real physical characteristics of impact events are actually rather complicated that the impact peak is followed by several smaller peaks and troughs. The current generalized Wagner method cannot fully capture these details. Nevertheless, Figure 20 in general shows reasonable agreement between experimental and numerical results.

**(a)**

**(b)**

**(a)**

**(b)**

##### 4.4. Hydroelastic Vibrational Vertical Bending Moments

Experimental and numerical nonlinear structural responses, extracted amidships for the two aforementioned operational conditions, were compared and analyzed in terms of the total vertical bending moments (VBM), low frequency VBM, and high frequency impact VBM (see Figures 21–23).

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

In general, there are satisfactory agreements between the experimentally and numerically determined VBM for these two operational conditions. In Figure 21(a), on average, the experimental results of the peaks slightly outweigh the numerical results, while the opposite phenomenon can be seen for the results of troughs. In contrast both the numerically obtained peaks and troughs outnumber experimental results in Figure 21(b). Figures 21(a) and 21(b) share a common feature that, for both experimental and numerical results, the values of the troughs are bigger than the values of the peaks which can be attributed to the characteristics of impact VBM demonstrated in Figure 23. The low frequency VBM in Figures 22(a) and 22(b) also shows similar asymmetry between peaks and troughs which can be mostly attributed to the nonlinearities associated with the vertical asymmetrical characteristics of the geometry of the hull. In Figure 22(a) the numerical method slightly underestimates peaks of the low frequency VBM, with the observation of relatively larger underestimation of troughs in Figure 22(b). From Figure 23, it can be seen that the numerical method underestimates the peak values while it overestimates the trough values compared with the experimental results. In addition, the impact VBM clearly manifests transient characteristics in these two sea conditions.

##### 4.5. Analysis on Effects of Nonlinearity

Linear numerical methods are applicable for the prediction of structural response of vessels in low-to-moderate seas. However, for large vessels in high seas, nonlinearities cannot be neglected. The purpose of this part is to determine the nonlinear effects on the numerical simulation.

In the linear case, the 2D impact pressure of (11) is integrated based on the unchanged wet contact line of the cross section of the vessel in the static hydroelastic state. The sectional vertical displacement of (10) is also calculated linearly. Figure 24 shows that this linear treatment largely underestimates the impact pressure compared with the proposed nonlinear method. In addition, after the impact peak, the following peak caused by the structural vibration feedback is also underestimated by the linear method.

**(a)**

**(b)**

Linear and nonlinear results of structural responses are compared (see Figures 25–27). In the linear model, all the generalized hydrodynamic forces are calculated based on the unchanged wet contact surface. As shown in these three figures, the linear results are generally smaller than the nonlinear and experimental results. The most evident difference between the linear case and the nonlinear case is that the peaks and troughs determined linearly are nearly the same. This is because the linear model cannot capture the change of the geometry of the wet contact surface caused by the motion of the vessel.

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

##### 4.6. Analysis on Effects of Elasticity

A conventional approach takes two steps to calculate whipping responses. The first step assumes that the vessel is a rigid body and the motion and VBM are computed based on this assumption. The second step considers the vessel as an elastic body and the results of motion in the first step serve as the input for the computation of impact VBM. Then the impact VBM is superimposed on the VBM in the first step for the total VBM. This part compares the results determined by this conventional method and the proposed nonlinear method to investigate the effects of elasticity on the impact pressure and structural responses. The comparisons are demonstrated from Figures 28–31. The results computed by the conventional model are tagged “rigid result” in Figures 28, 29, 30, and 31.

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

Figure 28 shows that the impact pressure calculated based on the rigid body outweighs the elastic body based results. Also the rigid body impact pressure decreases smoothly and monotonically after the impact peak. This is because the motion input for the computation of impact pressure is totally determined by the rigid body model and no structural vibration feedback is available. From Figures 29–31, it can be found that, in general, the rigid body method overestimates the structural responses. This is also due to the absence of elastic vibration feedback in the rigid model because the elastic vibration determined by elastic modes will mitigate the uneven external hydrodynamic pressure distribution and the mitigated external loads will reduce the structural responses.

#### 5. Conclusions

In the present research, impact pressure and vibrational responses of a large vessel with a pronounced bow flare were studied by a nonlinear hydroelastic method, and experimental tests were conducted to validate the proposed numerical model. Main findings are concluded as follows:(1)The fluctuations in the parameters of the experimentally generated regular waves seemed inevitable. Therefore, for the following comparison between experimental and numerical results, the corresponding numerical wave parameters needed to be modified based on the measuring data of the wave probe.(2)The existence of noise in the experimental impact pressure makes the filtering operation necessary. In some cases where the period of noise was very close to the duration of the impact, the exact peak impact pressures were difficult to obtain and should be treated as estimated values.(3)The impact pressures in the two representative operational conditions were overestimated by the fully coupled nonlinear method. As for the vibrational responses, the numerical and experimental results were in reasonable agreement.(4)The linear method underestimated the impact pressure and structural responses compared with the proposed nonlinear method and experimental model and produced similar peaks and troughs in VBM. This indicates that nonlinear effects should be considered in severe sea conditions.(5)The conventional rigid body method overestimated the impact pressure and the structural responses in comparison with the results determined by the proposed hydroelastic model and experiments. This is because of the inability of the rigid model to consider the elastic vibration feedback of the vessel. Therefore the hydroelastic model is much more suitable than the rigid body model in terms of computing structural response for large vessels sailing in high seas.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This study was in part supported by the National Natural Science Foundation of China (no. 51079034) and partly funded by China Scholarship Council (no. 201606680003). The support is highly appreciated. The warm help and decisive contributions from all the researchers who participated in the experimental tests and gave valuable advice to the authors are credited. The authors express sincere gratitude to China State Shipbuilding Corporation (CSSC) for the necessary support of details of vessel prototype. In addition, the administrative support from College of Shipbuilding Engineering of Harbin Engineering University is gratefully acknowledged.