Research Article  Open Access
The Nonlinear Hydroelastic Response of a SemiInfinite Elastic Plate Floating on a Fluid due to Incident Progressive Waves
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 semiinfinite 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 convergencecontrol parameter to obtain convergent solutions with relatively few terms. The clear calculation results are represented to show nonlinear waveplate 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 platecovered 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 semiinfinitely 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 smallamplitude waves. Taking a semiinfinite 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 semiinfinite 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 platecovered 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 semiinfinite 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 platecovered region were orthogonal. It is observed that the hydrodynamic behaviors of the floating semiinfinite elastic plate with different edge conditions, including a free edge, a simply supported edge, and builtin 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 platecovered 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 twolayer fluid, Lin and Lu [9] extended the study of Fox and Squire [3] to the case that a semiinfinite elastic plate floating on a twolayer 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 largeamplitude 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 semiinfinite elastic plates and water waves by using the wellknown perturbation method. Forbes [10] studied the nonlinear interaction between twodimensional 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 halfwave height and found that approximate solutions for the periodic waves had certain features in common with capillarygravity waves. Forbes [11] subsequently improved their original perturbation method by introducing NewtonRaphson 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 largeamplitude waves were also found to exist. With the perturbation method in a similar way, VandenBroeck and Părău [12] extended the study of Forbes [11] to the case of largeamplitude 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 semiinfinite 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 smallamplitude 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 semiinfinite 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 semiinfinite elastic plate floating on infinitely deep water for the twodimensional 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 semiinfinite 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 platecovered 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 semiinfinite elastic plate floating on infinitely deep water for the twodimensional 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 semiinfinite 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 platecovered 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 platecovered 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 fluidplate 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 platecovered region as follows:
Further, considering the continuities of the pressure and the velocity on the boundary between the platecovered 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 convergencecontrol 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 homotopybased multiplevariable 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 steadystate 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 steadystate 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 socalled zerothorder deformation equations:subject to the bottom conditionIn the open water region, the two boundary conditions on the unknown free surface ()where denotes a nonzero convergencecontrol 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 platecovered region, the socalled zerothorder deformation equations are the governing equation (12) which is subject to the bottom condition (13) and the following two boundary conditions on the unknown platecovered surface ():where the superscript identifies the platecovered 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 platecovered 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 socalled homotopyseries 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 platecovered 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 platecovered 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 platecovered 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 platecovered 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 platecovered 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. HighOrder Deformation Equations
The linear PDEs for the unknown and can be derived directly from the zerothorder deformation equations. Namely, differentiating the zeroorder deformation equations (12)–(15) for the open water region, (12), (13), and (18) for the platecovered region, and (22) and (23) at the boundary times with regard to , then dividing them by and setting , we have the socalled thorder deformation equations:subject to the bottom conditionand in the open water region, the two boundary conditions at (, )In the platecovered 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 platecovered 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 ConvergenceControl Parameter
If we fix all model physical parameters in our nonlinear boundaryvalue problem, there will still be an unknown convergencecontrol 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 convergencecontrol 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 convergencecontrol 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 convergencecontrol parameter , the square error can decrease quickly to at the 5thorder approximation, as listed in Table 1. These results indicate the validity and convergence of the corresponding HAMbased 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 zerothorder 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 HAMbased nonlinear solutions of hydroelastic response of a semiinfinite 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 homotopyseries 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 platecovered 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 platecovered 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 platecovered 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 platecovered 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 platecovered region.
5. Conclusions
By the homotopy analysis method, convergent series solutions are derived for the free wave surface elevation and the deflection of a semiinfinite 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 platecovered region, we obtain solutions with fairly low residual error after computing only a small number of terms by introducing a convergencecontrol parameter. Further, in the platecovered 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 homotopyseries 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 platecovered region also becomes larger. With the increase of Young’s modulus of the plate, the semiinfinite 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 platecovered 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 platecovered region.
Appendix
Detailed Derivations of , , , and in (40) and (41) for the Open Water Region and , , , and in (42) and (43) for the PlateCovered 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 likepower 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 platecovered 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.
References
 A. G. Greenhill, “Wave motion in hydrodynamics,” The American Journal of Mathematics, vol. 9, no. 1, pp. 62–96, 1886. View at: Publisher Site  Google Scholar  MathSciNet
 D. V. Evans and T. V. Davies, “Waveice interaction,” Tech. Rep., DTIC Document, 1968. View at: Google Scholar
 C. Fox and V. A. Squire, “On the oblique reflexion and transmission of ocean waves at shore fast sea ice,” Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences, vol. 347, no. 1682, pp. 185–218, 1994. View at: Publisher Site  Google Scholar
 C. Fox and V. A. Squire, “Reflection and transmission characteristics at the edge of shore fast sea ice,” Journal of Geophysical Research: Oceans, vol. 95, no. C7, pp. 11629–11639, 1990. View at: Publisher Site  Google Scholar
 V. A. Squire, J. P. Dugan, P. Wadhams, P. J. Rottier, and A. K. Liu, “Of ocean waves and sea ice,” Annual Review of Fluid Mechanics, vol. 27, no. 1, pp. 115–168, 1995. View at: Publisher Site  Google Scholar
 T. Sahoo, T. Yip, A. T. Chwang et al., “On the interaction of surface waves with a semiinfinite elastic plate,” in Proceedings of the 10th International Offshore and Polar Engineering Conference, International Society of Offshore and Polar Engineers, 2000. View at: Google Scholar
 B. Teng, L. Cheng, S. X. Liu, and F. J. Li, “Modified eigenfunction expansion methods for interaction of water waves with a semiinfinite elastic plate,” Applied Ocean Research, vol. 23, no. 6, pp. 357–368, 2001. View at: Publisher Site  Google Scholar
 F. Xu and D. Q. Lu, “An optimization of eigenfunction expansion method for the interaction of water waves with an elastic plate,” Journal of Hydrodynamics, vol. 21, no. 4, pp. 526–530, 2009. View at: Publisher Site  Google Scholar
 Q. Lin and D. Q. Lu, “Hydroelastic interaction between obliquely incident waves and a semiinfinite elastic plate on a twolayer fluid,” Applied Ocean Research, vol. 43, pp. 71–79, 2013. View at: Publisher Site  Google Scholar
 L. K. Forbes, “Surface waves of large amplitude beneath an elastic sheet. Part 1. Highorder series solution,” Journal of Fluid Mechanics, vol. 169, pp. 409–428, 1986. View at: Publisher Site  Google Scholar  MathSciNet
 L. K. Forbes, “Surface waves of large amplitude beneath an elastic sheet. Part 2. Galerkin solution,” Journal of Fluid Mechanics, vol. 188, pp. 491–508, 1988. View at: Publisher Site  Google Scholar  MathSciNet
 J.M. VandenBroeck and E. I. Părău, “Twodimensional generalized solitary waves and periodic waves under an ice sheet,” Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 369, no. 1947, pp. 2957–2972, 2011. View at: Publisher Site  Google Scholar
 P. A. Milewski, J.M. VandenBroeck, and Z. Wang, “Hydroelastic solitary waves in deep water,” Journal of Fluid Mechanics, vol. 679, pp. 628–640, 2011. View at: Publisher Site  Google Scholar
 S. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems [Ph.D. thesis], Shanghai Jiao Tong University, Shanghai, China, 1992.
 S.J. Liao, “On the homotopy multiplevariable method and its applications in the interactions of nonlinear gravity waves,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 3, pp. 1274–1303, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 D. Xu, Z. Lin, S. Liao, and M. Stiassnie, “On the steadystate fully resonant progressive waves in water of finite depth,” Journal of Fluid Mechanics, vol. 710, pp. 379–418, 2012. View at: Publisher Site  Google Scholar
 P. Wang and D. Lu, “Analytic approximation to nonlinear hydroelastic waves traveling in a thin elastic plate floating on a fluid,” Science China Physics, Mechanics and Astronomy, vol. 56, no. 11, pp. 2170–2177, 2013. View at: Publisher Site  Google Scholar
 S. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, vol. 2, CRC Press, Boca Raton, Fla, USA, 2004. View at: MathSciNet
 S. Liao, Homotopy Analysis Method in Nonlinear Differential Equations, 2012.
Copyright
Copyright © 2015 Ping Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.