#### Abstract

The nonlinear hydroelastic response of very large floating structures (VLFSs) or an ice sheet floating on the surface of deep water, idealized as a semi-infinite thin elastic plate, is investigated analytically in the case of nonlinear incident waves. Assuming that the fluid is inviscid and incompressible and the motion is irrotational, we consider incident progressive waves with a given angular frequency within the framework of potential flow theory. With the aid of the homotopy analysis method (HAM), the convergent analytical series solutions are derived by solving the simultaneous equations in which we apply a convergence-control parameter to obtain convergent solutions with relatively few terms. The clear calculation results are represented to show nonlinear wave-plate interaction. The effects of different physical parameters, including incident wave amplitude, Young’s modulus, the thickness and density of the plate on the wave scattering, and the hydroelastic response of the floating plate, are considered. We find that the variations of the plate stiffness, thickness, and density greatly change amount of wave energy which is reflected into the open water region and is transmitted into the plate-covered region. Further, the hydroelastic response of the plate also can be affected by the amplitude of incident wave.

#### 1. Introduction

In recent decades, there have been intensive researches on hydroelastic interaction between water waves and elastic plates in the area of ocean engineering and polar engineering in view of their academic research and applications. The very large floating structures (VLFSs), which are expected to serve as floating airports, artificial floating islands, ultralarge ships, or even mobile offshore bases, are mathematically assumed as elastic plates floating on a fluid for their strong flexibility and huge horizontal scales compared to vertical scales and characteristic wavelength. Accordingly, the hydroelastic deformation in the response of the structure to water waves is a predominant factor to design and maintain a safe and economic VLFS. Early studies in this field were mainly focused on hydroelastic interaction between the water waves and the large ice sheets which were also idealized as floating elastic plates [1], and the elastic plates are usually assumed to be semi-infinitely long in comparison with the wavelength of the incident wave [2, 3]. A comprehensive summary of the early work can be found in some review papers such as [2–5]. Unfortunately, mathematical methods on solving nonlinear problems are not well developed, so most of the present analytical researches on the hydroelastic interaction between water waves and an elastic plate still are in the scope of linear theory, which only can describe small-amplitude waves. Taking a semi-infinite elastic plate by using the eigenfunction expansion method, for example, Fox and Squire [4] investigated the reflection and transmission of ocean waves at the margin of an ice sheet which is idealized as a semi-infinite elastic plate and then determined the expansion coefficients via matching the velocity potential and the pressure at the interface of the open water region and the plate-covered region and enforcing the zero shear and bending moment at the free end. However, Fox and Squire [4] also stated that the beat Lagrange’s multiplier is very difficult to choose and the mathematical calculations become impossible for the unwise choice of the multipliers. Fox and Squire [3] subsequently improved their original method by removing the free end boundary conditions from the error function and then reducing the number of the Lagrange’s multiplier from three to one and analyzed the hydroelastic response of a semi-infinite elastic with a free edge due to the obliquely incident waves. Following the framework of Fox and Squire [3], Sahoo et al. [6] developed a new inner product based on the Fourier analysis within which the original eigenfunctions for the plate-covered region were orthogonal. It is observed that the hydrodynamic behaviors of the floating semi-infinite elastic plate with different edge conditions, including a free edge, a simply supported edge, and built-in edge condition, depend on the wave conditions, the geometrical settings, and the edge conditions. Meanwhile, Teng et al. [7] optimized Fox and Squire’s method [3] by removing the Lagrange multipliers from the error function completely and improved the method of Sahoo et al. [6] through employing the eigenfunctions in the open water region and the plate-covered region, respectively. It is demonstrated with examples that the solver for the associated linear system is simplified. In Xu and Lu [8], the eigenfunction expansion method was optimized by using the orthogonality property of eigenfunctions in the open water region. And numerical analysis showed that this modified method was effective and had higher convergence than the previous results in Sahoo et al. [6] and Teng et al. [7]. With the aid of the methods of matched eigenfunction expansion and the inner product of the two-layer fluid, Lin and Lu [9] extended the study of Fox and Squire [3] to the case that a semi-infinite elastic plate floating on a two-layer fluid of finite depth is subjected to the obliquely incident waves and found a critical angle for the incident waves of the surface wave mode and three critical angles for the incident waves of the interfacial wave mode, which are related to the existence of the propagating waves.

More importantly, it is inevitable that large-amplitude waves will occur more often in the future with the growing frequency of extreme weather events on a warmer and warmer Earth. Thus, some scholars investigated nonlinear hydroelastic interaction between semi-infinite elastic plates and water waves by using the well-known perturbation method. Forbes [10] studied the nonlinear interaction between two-dimensional periodic waves with a constant speed and an elastic ice sheet floating on a fluid of infinite depth by using the perturbation expansion in the half-wave height and found that approximate solutions for the periodic waves had certain features in common with capillary-gravity waves. Forbes [11] subsequently improved their original perturbation method by introducing Newton-Raphson techniques to approximate the Fourier coefficients and confirmed the existence of multiple solutions for the nonlinear equations describing the shape of nonlinear periodic waves, and extremely large-amplitude waves were also found to exist. With the perturbation method in a similar way, Vanden-Broeck and Părău [12] extended the study of Forbes [11] to the case of large-amplitude periodic waves and unfortunately found that the perturbation method cannot give any information for very steep waves because this method has to depend on small physical parameters. Further, Milewski et al. [13] applied the asymptotic and numerical methods to analyze the hydroelastic solitary waves propagating under a semi-infinite elastic plate. It is found that wave packet solitary waves bifurcate from nonlinear periodic waves of minimum speed for the unforced problem, and when the problem is forced by a moving load, steady responses are possible at all subcritical speeds for small-amplitude forcing, while there was a transcritical range of forcing speeds for which there are no steady solutions for larger loads. Finally, it is also noted that these elastic models do not have a clear conservation form for the elastic potential energy.

It is well known that the traditional perturbation and asymptotic techniques depend on the small physical parameters and approximations of nonlinear problems usually break down when the nonlinearity becomes strong. Therefore, they are only valid for weakly nonlinear problems. In this paper, we apply the homotopy analysis method (HAM) developed by Liao [14], a powerful analytic method for highly nonlinear problems, to consider the nonlinear hydroelastic response of a semi-infinite elastic plate floating on a fluid due to incident progressive waves. Further, we investigate the dynamic influences of some important physical parameters, including incident wave amplitude, Young’s modulus, and the density and the thickness of the plate, on nonlinear hydroelastic response of the floating plate and the wave scattering.

#### 2. Mathematical Formulation

We consider the nonlinear hydroelastic interaction between incident progressive waves and a semi-infinite elastic plate floating on infinitely deep water for the two-dimensional case, as shown in Figure 1. Cartesian coordinates are chosen such that represents the undisturbed water surface. The -axis points horizontally rightward and the -axis points vertically upward. The semi-infinite elastic plate, which floating on a fluid, expends rightward from to infinity along the -axis without draft. Then the whole fluid domain is divided into two regions: an open water region () and a plate-covered region (). Under the assumptions that the fluid is inviscid and incompressible and the motion is irrotational, the velocity potential satisfies the Laplace equation based on the potential flow theory:

We consider the nonlinear hydroelastic interaction between incident progressive waves and a semi-infinite elastic plate floating on infinitely deep water for the two-dimensional case, as shown in Figure 1. Cartesian coordinates are chosen such that represents the undisturbed water surface. The -axis points horizontally rightward and the -axis points vertically upward. The semi-infinite elastic plate floating on a fluid expends rightward from to infinity along the -axis without draft. Then the whole fluid domain is divided into two regions: an open water region () and a plate-covered region (). Under the assumptions that the fluid is inviscid and incompressible and the motion is irrotational, the velocity potential satisfies the Laplace equation based on the potential flow theory:where is the wave surface elevation. The bottom boundary condition at deep flat rigid seabed readsThe nonlinear kinematic boundary condition on the unknown surface can be written asIn the open water region, the nonlinear dynamic boundary condition on the unknown surface readswhere is the gravitational acceleration. Combining (4) and (5) gives the boundary conditionwhere

In the plate-covered region, under the assumption that any fluid particle which is once between the thin elastic plate and the water surface will remain on it, the dynamic boundary condition on the fluid-plate interface can be written aswhere is the fluid density. For a thin homogeneous elastic plate with uniform mass density and constant thickness , is the mass per unit area of the plate. is the flexural rigidity of the plate which is expressed as by Young’s modulus , Poisson’s ratio , and the thickness of the plate. We partially combine (4) and (8) to gain the boundary conditions for the plate-covered region as follows:

Further, considering the continuities of the pressure and the velocity on the boundary between the plate-covered region and the open water region, we express the matching conditions as

For the elastic plate freely floating on the water surface, it is reasonable that the motion of the plate satisfies the free edge conditions as follows:where is a given time. In this paper, we choose for matching conditions (10) and the free edge conditions (11). Now the unknown velocity potential and the plate deflection for two different water regions are governed by (2)–(4), (6), (9), (10), and (11). For the given angular frequency , the convergent analytical series solutions for and will be obtained based on the HAM in the subsequent section.

#### 3. Analytical Approach Based on the HAM

The homotopy analysis method (HAM) is an analysis approximation method developed by Liao [14] based on the concept of homotopy in algebraic topology. Different from the traditional perturbation method, the HAM is completely independent of any small or large physical parameters. Thus it is valid to highly nonlinear problems without any small or large parameter. More importantly, the HAM provides us with a convenient way to control and adjust the convergence of the approximate series solutions by means of introducing convergence-control parameters. Further, the HAM gives us great freedom in the choice of the initial guess, the basis solution functions, and the equation type of linear subproblems in order to simplify the solution procedure. Recently, the HAM has been successfully applied to solve many highly nonlinear problems in many fields. In hydrodynamics, for example, Liao [15] utilized the homotopy-based multiple-variable method to analyze the nonlinear interactions of periodic traveling waves in deep water. Xu et al. [16] further extended the study of Liao [15] to the case of steady-state system for the nonlinear interaction of two trains of waves propagating in water of definite depth by using the HAM in a similar way and found six different steady-state resonant waves. It is noted that Wang and Lu [17] employed the HAM to consider nonlinear hydroelastic progressive waves traveling in an elastic plate in deep water. And graphical representations of the series solutions show that several physical parameters have important influences on the hydroelastic response of the elastic plate floating on the surface of deep water. For details on the HAM, please refer to the two books by Liao [18, 19].

##### 3.1. Continuous Variations

Let denote the embedding parameter and let and denote the initial approximations of the velocity potential and wave elevation , respectively. Instead of solving these nonlinear equations (2)–(4), (6), (9), (10), and (11) directly, in the frame of the HAM we construct two continuous variations and which are called homotopies in topology, such that as increases from 0 to 1, deforms continuously from its initial approximation to the exact solution of the original problem, so does from the approximation to . Those two continuous variations are governed by a new family of nonlinear partial differential equations, namely, the so-called zeroth-order deformation equations:subject to the bottom conditionIn the open water region, the two boundary conditions on the unknown free surface ()where denotes a nonzero convergence-control parameter and the superscript identifies the open water region. is an auxiliary linear operator with the property . and are nonlinear differential operators defined bywith the definition

Similarly, in the plate-covered region, the so-called zeroth-order deformation equations are the governing equation (12) which is subject to the bottom condition (13) and the following two boundary conditions on the unknown plate-covered surface ():where the superscript identifies the plate-covered region. and are auxiliary linear operators with , (). and are nonlinear differential operators defined byNote that the definitions for , , , and are based on the four nonlinear boundary conditions (6) and (5) and (9) and (8), respectively. It should be emphasized that the velocity potential and wave elevation are defined as follows:And their initial approximations are defined byIn addition, the matching conditions of the velocity and pressure at the boundary between the plate-covered region and the open water region are deformed intoAnd the free edge conditions are

In order to obtain the series solutions for the velocity potential and wave elevation , by Taylor’s theorem, we expand and into the Taylor series about the embedding parameter at aswhere

Assuming that could be so properly chosen that the above two series in (24) converge at , then we have the so-called homotopy-series solutions aswith the definitions

As shown later in Section 3.3, the unknown terms and can be obtained by solving linear PDEs (38)–(41) for the open water region, (38), (39), (42), and (43) for the plate-covered region, and (44)–(45) at the boundary.

##### 3.2. Solution Expressions

Considering the progressive gravity wave propagating periodically in with a given angular frequency , the progressive waves can be expressed as the combinations of the wave traveling to the right and the wave traveling to the left in the open water region. From viewpoints of the physical background, we represent the progressive wave elevation in the open water region as follows:where and are unknown coefficients to be determined. is the wave number for the open water region which satisfies [16] , where the value of is not equal to 1. Since the dispersion relation is valid only for the linear waves, in this paper we consider the case of .

According to the linear wave theory, the governing equation (12) and the bottom condition (13), in the open water region the potential function , should be in the formwhere and are unknown coefficients to be determined.

In the plate-covered region, if there is no gap between the bottom surface of the thin elastic plate and the top surface of the fluid layer, it is reasonable that the plate or the surface fluid still travels periodically in the . Then the plate deflection is mathematically represented as the form with (28):where is an unknown coefficient. is the wave number for the plate-covered water region which satisfies [16], where the linear dispersion relations[8]:With loss of generality, we consider the case of too.

Subsequently, we consider the expression of the potential function aswhere is an unknown coefficient to be determined. Note that the potential function defined by (32) also automatically satisfies the governing equation (12) and the bottom condition (13).

Equations (28) and (29) are called the solution expressions of and in the open water region, and in the plate-covered region (30) and (32) are the solution expressions, respectively, which have important roles in the frame of the HAM.

As mentioned by Liao [19], according to the solution expressions (29) and (32), we have extremely large freedom to choose the initial approximations and the auxiliary linear operator which greatly simplify the process of solving nonlinear problems. Then in the two regions we construct the initial approximations of the potential functionrespectively, where , , and are unknown coefficients. And we chooseas the initial guess of for the open water region and the plate-covered region, respectively, where is a given amplitude of the initial incident progressive waves. and are unknown coefficients.

It should be noted that the HAM can provide extremely large freedom in the choice of auxiliary linear operators; then for simplifying the computation procedure we can only choose linear operators of in and as the auxiliary linear operators and by means of the solution expressions (29) and (32) as follows:And we choose linear operators of in as the auxiliary linear operators by means of the solution expression (30) as

##### 3.3. High-Order Deformation Equations

The linear PDEs for the unknown and can be derived directly from the zeroth-order deformation equations. Namely, differentiating the zero-order deformation equations (12)–(15) for the open water region, (12), (13), and (18) for the plate-covered region, and (22) and (23) at the boundary times with regard to , then dividing them by and setting , we have the so-called th-order deformation equations:subject to the bottom conditionand in the open water region, the two boundary conditions at (, )In the plate-covered region, the governing equation (38) which is subject to the bottom condition (39) and the following two boundary conditions at (, )where and is the Heaviside step function.

The elaborated expressions for , , , and , , , are given in the Appendix.

At the boundary between the plate-covered region and the open water region, the matching conditions of the velocity and pressure of the free edge conditions areAnd the free edge conditions are

Now we see that the unknown terms and are governed by linear PDEs (38)–(41) for the open water region, (38), (39), (42), and (43) and (44) and (45) at the boundary. More importantly, these subproblems for and are not only linear but also decoupled and then can easily be solved.

##### 3.4. Optimal Convergence-Control Parameter

If we fix all model physical parameters in our nonlinear boundary-value problem, there will still be an unknown convergence-control parameter which is used to guarantee the convergence of the series solutions. The optimal value of corresponds to the fastest decrease of the total squared residuals of our nonlinear problem. The total squared residual is defined bywith the definitionswhere , , , and are given in Appendix, is the number of the discrete points, and . In this paper, we choose , and then the optimal convergence-control parameter is obtained by the minimum of .

#### 4. Results and Discussion

To verify the convergence of our approximate analytic solutions to the nonlinear hydroelasticity problem considered here by means of the HAM, we employ the dimensional case of m^{−1}, m, kgm^{−3}, , Pa, kgm^{−3}, s, and m and take these data hereinafter for calculation unless otherwise stated. The total residual square error at different orders of approximation versus the convergence-control parameter is shown in Figure 2. It indicates that as the order of approximation increases gradually, the total residual square error decreases in the interval , and the optimal value of is close to . Further, when the optimal convergence-control parameter , the square error can decrease quickly to at the 5th-order approximation, as listed in Table 1. These results indicate the validity and convergence of the corresponding HAM-based analytic solutions for velocity potential and the water deflection governed by (2)–(4), (6), (9), (10), and (11). It should be emphasized that the optimal in zeroth-order deformation equations such as (14) and (15) and (18) can help us to obtain relatively few terms of approximate solutions which have sufficiently small residual error.

Here we only change one physical parameter in the above case to exclude the effects of other ones on the nonlinear hydroelastic interaction between the incident progressive waves and the elastic plate. Firstly, we compare the HAM-based nonlinear solutions of hydroelastic response of a semi-infinite elastic plate floating on a fluid due to incident progressive waves with corresponding linear results also obtained by HAM. As shown in Figure 3, it can be seen clearly that our nonlinear homotopy-series approximation of the surface elevation varies more greatly than the linear approximation, which further indicates that hydroelastic response of the floating plate will be underestimated if only linear theory is considered. In Figures 4–6, we give the incident progressive wave profiles and hydroelastic deflections of the plate for different physical parameters. From the comparison between the open water and the plate-covered regions, it is easy to find that hydroelastic deflections of the plate are much smaller than the incident progressive wave profiles. This is due to the fact that, for the existence of a floating elastic plate, most of the wave energy is concentrated near the free surface. In other words, a large proportion of the wave energy is reflected back by the plate to the open water region and a small amount of wave energy is transmitted to the plate-covered region.

Practically speaking, in Figure 4 we give hydroelastic deflections of the plate and the incident wave profile by increasing incident wave amplitude from 0.05 to 0.5. It is easy to find that the amplitudes of the plate defection increase with the increment of the amplitude . The effect of Young’s modulus of the plate is shown in Figure 5. And we can see that the incident progressive wave elevation in the open region increases while the hydroelastic deflection of the plate decreases as Young’s modulus increases. This means the elastic plate becomes stiffer as Young’s modulus increases, and then most of the incident progressive waves which concentrate near the interface between the plate-covered region and the open water region are reflected back, while less wave energy is transmitted below the plate. Figure 6 shows the differences of the wave deflections in the open water region and the plate-covered region for the plate thicknesses m, m, and m. It is found that a large thickness decreases the plate deflection while increasing the interfacial wave elevation. And as shown in Figure 7, the plate density has analogous influence on the incident progressive wave elevation in the open region and the hydroelastic deflection of the plate in the plate-covered region.

#### 5. Conclusions

By the homotopy analysis method, convergent series solutions are derived for the free wave surface elevation and the deflection of a semi-infinite elastic plate floating on a fluid of infinite depth, respectively, to analytically study the nonlinear hydroelastic response of the plate with a free edge due to incident progressive waves. For the whole flow domain mathematically including an open water region and a plate-covered region, we obtain solutions with fairly low residual error after computing only a small number of terms by introducing a convergence-control parameter. Further, in the plate-covered region, it is noted that all of two nonlinear operators and contain the linear operators for both and mathematically. Fortunately, under the guarantee of the convergence of homotopy-series solutions, we may choose the auxiliary linear operator containing the derivatives of only and the auxiliary linear operator containing the derivatives of only, which greatly simplify the calculation of nonlinear hydroelastic interaction. In addition, the relation of series solutions for these two regions is obtained by matching the velocity and pressure at boundary (11). Numerical results show that all our series solutions are highly convergent and accurate.

The influences of several physical parameters on the hydroelastic responses of the plate are considered in detail. The amplitude of incident wave, Young’s modulus , the thickness , and density of the floating elastic plate have effects on the free wave surface elevation and the plate deflection. As the amplitude of the incident wave increases, the plate deflection in the plate-covered region also becomes larger. With the increase of Young’s modulus of the plate, the semi-infinite elastic plate deflection steps down, while the amplitude of free surface water elevation in the open water region becomes larger correspondingly. It is also found that the larger value of the thickness or density of the plate can sharpen the crest and flatten the trough of the water elevation for the open water region significantly, and in the plate-covered region the plate deflection decreases. All of those results indicate that when incident progressive waves in the open water region interact on a stiffer or thicker plate, more wave energy is reflected into the open water region, while less wave energy is transmitted into the plate-covered region.

#### Appendix

#### Detailed Derivations of , , , and in (40) and (41) for the Open Water Region and , , , and in (42) and (43) for the Plate-Covered Region

LetWe have a Maclaurin series for any :For , it follows from (A.1) and (A.2) thatwhereThus we have, for ,whereSubstituting expansions (A.1) and (A.5) into the boundary conditions (14) and (15), then equating the like-power of , we have two linear boundary conditions (40) and (41) for the open water region, respectively. Similarly, we can have two linear boundary conditions (42) and (43) for the plate-covered region, respectively. And the explicit expressions for , , , and , , , in these conditions are given bywhere

#### Conflict of Interests

There is no conflict of interests in this paper. The author obtained the approximate analytical solutions of the PDEs considered here by means of computer algebra system like Mathematica.

#### Acknowledgments

This research was sponsored by the Natural Science Foundation of Shanghai under Grant no. 14ZR1416200, the National Natural Science Foundation of China under Grant no. 11472166, and the Natural Science Foundation of Shandong Province of China under Grant no. ZR2013AL012.