Effect of Mooring Lines on the Hydroelastic Response of a Floating Flexible Plate Using the BIEM Approach

: A boundary integral equation method (BIEM) model for the problem of surface wave interaction with a moored ﬁnite ﬂoating ﬂexible plate is presented. The BIEM solution is obtained by employing the free surface Greens function and Green’s theorem, and the expressions for the plate deﬂection, reﬂection, and transmission coefﬁcients are derived from the integro-differential equation. Furthermore, the shallow water approximation model and its solution is obtained based on the matching technique in a direct manner. The accuracy of the present BIEM code is checked by comparing the results of deﬂection amplitude, reﬂection, and transmission coefﬁcients with existing published results and experimental datasets as well as the shallow water approximation model. The hydroelastic response of the moored ﬂoating ﬂexible plate is studied by analyzing the effects of the mooring stiffness, incidence angle, and ﬂexural rigidity on the deﬂection amplitude, plate deformations, reﬂection, and transmission coefﬁcients. The present analysis may be helpful in understanding the different physical parameters to model a wave energy conversion device with mooring systems over BIEM formulations.


Introduction
In recent years, researchers have been increasingly interested in studying the hydroelastic response of floating or submerged structures connected with mooring lines application to coastal or marine engineering problems in order to model floating breakwaters and wave energy converters (see [1][2][3][4]). These floating flexible structures need to survive the failure or any kind of damage in waves, wind, and currents during storm events. To model a floating flexible structure, the hydroelastic response and its analysis in different design parameters ensuring the reliability and safety of the flexible structures are of great importance as the floating structures can be located in unprotected seas and exposed to rigorous environmental conditions. Due to the flexible nature of the horizontal structure, the computational burden becomes too bulky. To overcome these difficulties, researchers often use numerical methods and solutions. The wave response analysis of a flexible floating structure by combining the two numerical methods boundary element method and finite element method (BEM-FEM) is shown in [5]. A B-spline Galerkin scheme was applied to calculate the hydroelastic response of a very large floating structure in waves in [6]. The hydroelastic response of a flexible floating structure was analyzed by applying a numerical method FEM in [7]. The numerical methods BEM and FEM have been applied to the hydrodynamic analysis of floating elastic surfaces in [8] using FEM; in [9] using BEM and the direct method in [10]. The hydroelastic response of a compliant floating raft modelled as sea ice flow using the boundary-integral method in [11]. The hydroelastic waves in a channel covered by cracked ice with arbitrary edge conditions at the walls expanded the solution into a series of cosine functions with unknown coefficients was studied in [12].
It may be noted that the numerical computations of the present solution based on integro-differential equation were carried out in a desktop machine with Intel ® core i7-4790 CPU with a 3.60 GHz processor and 16 GB of RAM. On average, each case was run roughly 8-10 min (simulation time). The present BIEM would be one of the alternative methods to compare and analyze the hydroelastic response of this problem by incorporating the structural deformations and flexible modes into the formulation. This method is computationally efficient and it also involves less computational cost than the different numerical methods used in [13][14][15].
On the other hand, the interaction of waves with floating or submerged flexible structures of finite depth connected with mooring lines can hardly be analyzed analytically. Some of the most relevant previous work relating to moored finite floating and or submerged flexible horizontal plate problems are discussed below.
There has been considerable progress in the study of real physical boundary value problems associated with moored finite horizontal floating structures in different water depths based on an analytical approach. Surface wave interaction with a moored finite floating elastic plate was formulated using linear water wave theory and the effect of different design parameters on the floating elastic plate were analyzed using the eigenfunction expansion method in [16]. The effect of mooring lines on the floating elastic plate in the presence of submerged flexible horizontal membrane under eigenfunction expansion method and application of orthogonal mode coupling relation was studied in [17]. Under the velocity decomposition method, the effect of a fixed submerged flexible membrane on a moored floating elastic plate in three dimensions was analyzed in [18]. Furthermore, the effect of a moored submerged flexible porous plate on the moored finite floating flexible plate was analyzed in two dimensions by studying different design parameters in [19]. Recently, a review on the numerical approaches in the hydroelastic response of horizontal elastic structures was presented in [20].
From the above literature, it is confirmed that until now there is no study reported to the public related to wave interaction with a moored finite floating flexible plate using BIEM and its hydroelastic response via BIEM. Therefore, this paper presents a BIEM model for linear water wave interaction with a moored floating flexible plate in finite water depth to analyze the effect of mooring lines on the floating flexible structure on various design parameters.
The derivation of the boundary integral equation along with the integrodifferential equation in terms of plate deflection for the moored finite floating flexible plates are presented by using three-dimensional free surface Green's function and applying Green's theorem. A BIEM code is developed to present numerical results of deflection amplitude, plate displacements, reflection, and transmission coefficients to study the hydroelastic response of a moored floating flexible plate. To show the level of accuracy of the BIEM code, the present results of deflection amplitude, reflection, and transmission coefficients are compared with existing published results and experimental datasets as well as the shallow water approximation (SWA) model. Furthermore, the hydroelastic response of moored finite floating flexible plate is studied by analyzing the effects of mooring stiffness, incident angle, and flexural rigidity of the moored floating flexible plate on the deflection amplitude, plate displacements, reflection, and transmission coefficients.

Model Formulation
Under linear wave theory, the mathematical formulation is considered in a threedimensions Cartesian Coordinate system and BVP is solved by applying BIEM in finite water depth. It is assumed that a floating flexible plate is placed at y = 0, 0 < x < b connected with mooring lines with stiffness s j , and z < ∞ over the impermeable bottom bed y = h. Hence, the whole fluid domain is divided into three fluid domains defined by open water domain, plate covered domain, and interface domain are referred as D O , D P , and ∂D, respectively (see Figure 1).
where g is the gravitational constant. The linearized kinematic condition at 0 y = in plate covered domain is given by where ( , ; ) x z t η is the deflection of the floating flexible plate. The floating plate condition is obtained by satisfying plate deflection, equating the hydrodynamic and dynamic pressure along with kinematic condition (4) at 0 y = as: where  and  are the flexural rigidity and the uniform compressive force of the plate, respectively.
At the moored edges 0, ; 0 x b y = = , the bending moment and shear force yield the conditions as ( ) It is considered that surface gravity wave interacts with floating flexible plate making an oblique angle β with the positive x-axis and the fluid is inviscid, incompressible, motion is irrotational and simple harmonic in time with angular frequency ω. Therefore, there exists velocity potential Φ(x, y, z; t) that can be defined as Φ(x, y, z, t) = Re φ(x, y, z)e −iωt where φ(x, y, z) is the complex wave potential with φ(x, y, z) = φ O (x, y, z) for D O and φ(x, y, z) = φ P (x, y, z) for D P which satisfies the Laplace equation in the whole fluid domain as and the impermeable sea bed condition, The linearized free surface condition in D 1 is obtained by the combination of kinematic and dynamic boundary conditions yields where g is the gravitational constant. The linearized kinematic condition at y = 0 in plate covered domain is given by where η(x, z; t) is the deflection of the floating flexible plate. The floating plate condition is obtained by satisfying plate deflection, equating the hydrodynamic and dynamic pressure along with kinematic condition (4) at y = 0 as: where X and are the flexural rigidity and the uniform compressive force of the plate, respectively. At the moored edges x = 0, b; y = 0, the bending moment and shear force yield the conditions as The continuity of pressure and velocity at the interface of the plate covered and the water surface are given by Finally, the complex wave potential φ(x, y, z) also satisfies the Sommerfeld radiation condition in the far-field.

Integro-Differential Equation Based on BIEM
In this Section, the solution of wave interaction with a moored floating flexible plate via integral equation and plate deflection in terms of integro-differential equation will be presented. The complex velocity potential φ O (x, y, z) is decomposed into the incident wave potential and the diffracted wave potential, which are denoted by φ I and φ S , respectively. So, the total potential in the domain D O can be expressed as where with ε(y) = cosh k(h − y)/cosh kh, I 0 is the incident wave amplitude, ω is the frequency and k 0 satisfies the gravity wave dispersion relation. Let G(x, y, z; x 0 , y 0 , z 0 ) refer to the three-dimensional Green's function with (x 0 , y 0 , z 0 ) and (x, y, z) are the source point and any point in the fluid domain, respectively. Therefore, the Green's function G(x, y, z; x 0 , y 0 , z 0 ) at the free surface y = y 0 = 0 is obtained as where C is the contour of integration, J 0 (αϕ) is the Bessel's function with and satisfies the governing equation with Dirac delta function δ along with free surface boundary condition (3) and bottom boundary condition (2). Now, applying Green's theorem to the velocity potentials φ O and φ P , for the domain (x, z) ∈ D O , one can get (14) and for (x, z) ∈ D P , result in where n denotes the outward normal to Ω (Ω is a smooth (regular) surface). By adding the two Equations (15) and using boundary condition (5) along with Equation (3) for the Green's function, one can derive as: The integro-differential equation can be derived by applying Green's formula, free surface condition for Green's function and Laplace equation to incident potential as where ϑ = ω 2 /g. Using the kinematic and free surface boundary conditions as in Equations (2) and (3) and substituting incident potential from Equation (11) into Equation (17), an integro-differential equation in terms of plate deflection η(x, z) is obtained as: The Greens function (12) satisfy the free surface boundary condition at y = y 0 = 0, bottom boundary condition, and radiation condition where the contour C defined over the integral is the complex α-plane from 0 to +∞ beneath the singularity α = α 0 (see Figure 2) chosen for satisfying the radiation condition and Bessel's function J 0 (αϕ). To overcome the singularity under the integrals of the integro-differential Equation (18) associated with the Green's function by using residue theorem based on the complex function theory. Furthermore, we applied the Sonine-Gegenbauer expression for J 0 (αϕ) to evaluate the integration with respect to z 0 (see [21]).
Using the kinematic and free surface boundary conditions as in Equations (2) and (3) and substituting incident potential from Equation (11) into Equation (17), an integro-differential equation in terms of plate deflection ( , ) x z η is obtained as: The Greens function (12) satisfy the free surface boundary condition at 0 0 y y = = , bottom boundary condition, and radiation condition where the contour C defined over the integral is the complexα plane from 0 to +∞ beneath the singularity 0 α α = (see Figure 2) chosen for satisfying the radiation condition and Bessel's function 0 ( ) J αϕ . To overcome the singularity under the integrals of the integro-differential equation (18) associated with the Green's function by using residue theorem based on the complex function theory. Furthermore, we applied the Sonine-Gegenbauer expression for 0 ( ) J αϕ to evaluate the integration with respect to 0 z (see [21]). In the next subsection, the displacement of the moored floating flexible plate, reflection and transmission coefficients will be determined in terms of series form.

Plate Displacement, Reflection and Transmission Coefficients
In this Section, the plate deflection can be expressed as a linear superposition of the eigenfunctions associated with the flexural gravity waves in the following form 0, , , , where n A and n B are the unknown amplitudes are to be determined with 0 sin κ α β = .
The wavenumber n γ satisfies the relation  In the next subsection, the displacement of the moored floating flexible plate, reflection and transmission coefficients will be determined in terms of series form.

Plate Displacement, Reflection and Transmission Coefficients
In this Section, the plate deflection can be expressed as a linear superposition of the eigenfunctions associated with the flexural gravity waves in the following form where A n and B n are the unknown amplitudes are to be determined with κ = α 0 sin β. The wavenumber γ n satisfies the relation γ n = λ 2 n − κ 2 where λ n satisfies the flexural gravity dispersion relation.
In Equation (20), S(λ n ) ≡ 0 is the dispersion relation associated with a floating horizontal flexible plate in finite water depth.
Substituting the Green's function as in Equation (12) and plate deflection as in Equation (19) into Equation (18), a system of two linear equations is obtained and another four linear equations are obtained using the moored edge conditions (6) and (7), the bending moment and shear force at Thus, there are six equations for the determination of the unknowns A n and B n . Using the relations in Equation (14), the amplitudes of reflected and transmitted waves R and T in terms of A n and B n are obtained as: where In the next section, the referred problem will be modelled in SWA in a direct manner.

Shallow Water Approximation Model
Proceeding similarly in the case of SWA as in [22,23], using the kinematic condition (4) and plate covered condition (5) Solving the Laplace equation using separation of variables along with free surface boundary condition, the velocity potential in D 0 associated with the reflection and transmission wave amplitudes R and T can be derived as: Using the conditions of pressure and velocity (8)(9)) at x = 0, b into the Equations (27) and (28), we obtain four linear equations and, also, we have another four linear equations from moored edge conditions (21)- (24). A MATLAB code is developed by using the system of eight linear equations to obtain the unknowns A n , B n , R, and T associated with Equations (27) and (28). Once the reflection and transmission wave amplitudes are determined, then the reflection and transmission coefficients can be computed by the following formulae C r = |R| and C t = |T|.

Numerical Results and Discussions
In numerical computation, the roots of the dispersion relation (20) in D P are in the upper complex plane. The dispersion relation (20) has two real complex roots, four complex roots, and many infinitely imaginary roots. In this work, for developing a numerical code, 10 different roots are utilized which give realistic physical solutions for λ and are positioned in the upper complex half-plane: one real positive root λ 0 , two complex roots λ 1 and λ 2 which are of the form ±a + ib, and seven imaginary roots of the form +iλ n for n = 3, 4, . . . , 9. Initially, to verify the accuracy of the present BIEM code, the obtained results of deflection amplitude, reflection and transmission coefficients are compared with existing published results of other calculations, experimental datasets, and SWA models. Then, the hydroelastic response of moored finite floating plate is analyzed by studying the effects of mooring stiffness, incidence angle, and flexural rigidity on the deflection amplitudes and plate displacements as well as reflection and transmission coefficients in detail. Figure 3 compares the present deflection amplitude with existing published results in [24] and experimental dataset in [25] with EI = 10 5 , = 1.5 EIρg, and b/h = 20 versus plate length. The comparison results show that the trend of the deflection amplitude between the present and [24] are similar. Slightly higher values of mooring stiffness s = 10 3.5 Nm −1 and lower compressive force = 0.02 EIρg were chosen which led to a similar trend between the results from the model [24] and experimental data [25] and the present BIEM model.  Figure 3. Comparison of the deflection amplitude with the model [24] and experimental data in [25]. , and / 20 b h = . In Figure 4a, the comparison shows that the trend of t C is in a similar pattern however, there are some deviations between the two models. This is because of the nonlinearity in [26] and the linearity in the present BIEM model. From Figure 4b,c, it is found that the trend of the r C and t C are similar and also the data between the BIEM and SWA are very close to each other. It is also seen that the SWA results are over or under predicted which is because of the neglection of evanescent modes in the case of SWA.

Comparison of Results with Existing Published Results, Experimental Datasets, and SWA Model
(a) Comparison of the t C : present (BIEM) and Green-Naghdi model [26]. In Figure 3, comparison results indicated that few of the experimental data points with BIEM are in agreement and not all data points but the trends are almost similar. These differences can be explained by the present BIEM model associated with the mooring lines connected at the edges of the flexible plate under compressive force, while these parameters were not included in the model [25]. Therefore, the present BIEM model reproduces a lower amplitude crest of the floating flexible plate that the model [25] cannot reproduce. Furthermore, in Figure 3, the deflection amplitude of the horizontal floating flexible plate from the [24] model is 0.09 m, while the BIEM model is 0.08 m, which is 11% larger. On the other hand, the deflection amplitude from the Experimental data [25] model is 0.2 m, while the present BIEM model is 0.8 m, which is 66% larger. However, the other two points from the model [25] are in almost agreement with the present BIEM model. Figure 4 shows the comparisons of C t between the present BIEM with (a) Green-Naghdi model ( [26], Figure 13), SWA of the (b) C r and (c) C t versus non-dimensional wavelength with EI = 10 5 , = 1.5 EIρg, and b/h = 20. In Figure 4a, the comparison shows that the trend of C t is in a similar pattern however, there are some deviations between the two models. This is because of the nonlinearity in [26] and the linearity in the present BIEM model. From Figure 4b,c, it is found that the trend of the C r and C t are similar and also the data between the BIEM and SWA are very close to each other. It is also seen that the SWA results are over or under predicted which is because of the neglection of evanescent modes in the case of SWA. , and / 20 b h = . In Figure 4a, the comparison shows that the trend of t C is in a similar pattern however, there are some deviations between the two models. This is because of the nonlinearity in [26] and the linearity in the present BIEM model. From Figure 4b,c, it is found that the trend of the r C and t C are similar and also the data between the BIEM and SWA are very close to each other. It is also seen that the SWA results are over or under predicted which is because of the neglection of evanescent modes in the case of SWA.
(a) Comparison of the t C : present (BIEM) and Green-Naghdi model [26].  Furthermore, in Figure 4a, the transmission coefficient from the Green-Naghdi theory model is 0.956 m, while the present BIEM model is 1 m, which is 4.37% larger. In Figure 3b: The reflection coefficient from the shallow water approximation model is 1 m, while the BIEM model is 0.9 m, which is 9% larger. In Figure 4c: The transmission coefficient from the Shallow water approximation model is 0.943 m, while the BIEM model is 0.84 m, which is 10.25% larger.
As a result of Figures 3 and 4, it is ensured that the present solution based on BIEM, is supported by the existing published results, experimental datasets and as well as SWA model.

Hydroelastic Response of the Moored Floating Flexible Plate Based on BIEM
In Figure 5 , and / 20 b h = are plotted. From Figure 5a,b, we can observe that the deflection amplitude decreases with an increase in the values of mooring stiffness whilst, the deflection amplitude becomes higher for higher incidence angle (shown in Figure 5c) which is because,when the incident Figure 4. Comparison of present BIEM with Green-Naghdi model [26] and SWA of the reflection and the transmission coefficients C r and C t . Furthermore, in Figure 4a, the transmission coefficient from the Green-Naghdi theory model is 0.956 m, while the present BIEM model is 1 m, which is 4.37% larger. In Figure 4b: The reflection coefficient from the shallow water approximation model is 1 m, while the BIEM model is 0.9 m, which is 9% larger. In Figure 4c: The transmission coefficient from the Shallow water approximation model is 0.943 m, while the BIEM model is 0.84 m, which is 10.25% larger.
As a result of Figures 3 and 4, it is ensured that the present solution based on BIEM, is supported by the existing published results, experimental datasets and as well as SWA model.

Hydroelastic Response of the Moored Floating Flexible Plate Based on BIEM
In Figure 5, the variations of the deflection amplitude of the floating flexible plate for different (a) mooring stiffness versus plate length, (b) mooring stiffness versus nondimensional wavenumber, (c) incidence angle with a certain value of mooring stiffness s = 10 3.5 Nm −1 , and (d) flexural rigidity for mooring stiffness s = 10 3.5 Nm −1 versus non-dimensional wavenumber kb with EI = 10 5 , = 1.5 EIρg, and b/h = 20 are plotted. From Figure 5a,b, we can observe that the deflection amplitude decreases with an increase in the values of mooring stiffness whilst, the deflection amplitude becomes higher for higher incidence angle (shown in Figure 5c) which is because, when the incident wave becoming perpendicular to the floating plate, it propagates along the floating horizontal plate in the z-direction. This is because, as the rigidity parameter of the floating increases, the plate behaves as a rigid structure which leads to less deflection or no deflection. From Figure 5d, it is observed that with a certain value of mooring stiffness s = 10 3.5 Nm −1 , the deflection amplitude of the floating flexible plate decreases with an increase in flexural rigidity.  wavelengths, the resonating pattern of the shear force vanishes. In Figure 7a,b, it is found that the observations are similar as in Figure 6a,b. Furthermore, it was observed that the bending moment attends higher values than those of shear force whilst, this trend is the reverse in case of the effect of incident angle as seen from Figures 6b and 7b with a certain 3.5 1 − Figure 5. Effects of (a,b) mooring stiffness, (c) incidence angle β with s = 10 3.5 Nm −1 , and (d) flexural rigidity with s = 10 3.5 Nm −1 on the deflection amplitude with EI = 10 5 , = 1.5 EIρg, and b/h = 20. Figures 6 and 7 show the effects of different mooring stiffness and incident angle on the bending moment and shear force on the moored floating flexible plate versus nondimensional wavelength L/b with EI = 10 5 , = 1.5 EIρg, and b/h = 20. From Figure 6a, it is observed that the bending moment acting on the floating flexible plate decreases with an increase in the values of mooring stiffness. However, for larger wavelengths, the resonating pattern of the shear force vanishes. In Figure 7a,b, it is found that the observations are similar as in Figure 6a,b. Furthermore, it was observed that the bending moment attends higher values than those of shear force whilst, this trend is the reverse in case of the effect of incident angle as seen from Figures 6b and 7b with a certain value of mooring stiffness s = 10 3.5 Nm −1 .    . It is observed that the plate deformation decreases with an increase in the values of mooring stiffness. This is because as the mooring stiffness of the plate increases, the plate becomes stiffer which leads to less deformation of the plate.   . It is observed that the plate deformation decreases with an increase in the values of mooring stiffness. This is because as the mooring stiffness of the plate increases, the plate becomes stiffer which leads to less deformation of the plate.   , and / 20 b h = . In Figure 9a, it is seen that for higher values of mooring stiffness, the reflection coefficient r C increases for / 8 L b > which is because as the stiffness becomes higher the flexible plate becomes stiffer which results in higher reflection (see [18]). On the other hand, in the case of t C the observations are similar to r C throughout the non-dimensional wavelength / L b . In Figure 9b, it is found that the values of the r C becomes lower as the values of incident angle increase whilst, the transmission coefficients decrease for higher values of incident angle. This is because as the stiffness becomes higher, the flexible plate becomes stiffer which results in a higher reflection (see [18]). Furthermore, the observations in Figure 9c are similar to Figure 9a. From Figure 9d, it is observed that the reflection coefficient increases with an increase in flexural rigidity whilst, the trend of transmission coefficient becomes the reverse to that of reflection coefficients.  Figure 9 demonstrates the effect of the mooring stiffness and incident angle on the reflection coefficient C r and the transmission coefficient C t versus (a) non-dimensional wavelength L/b and (b) wave period T(s) with EI = 10 5 , = 1.5 EIρg, and b/h = 20. In Figure 9a, it is seen that for higher values of mooring stiffness, the reflection coefficient C r increases for L/b > 8 which is because as the stiffness becomes higher the flexible plate becomes stiffer which results in higher reflection (see [18]). On the other hand, in the case of C t the observations are similar to C r throughout the non-dimensional wavelength L/b. In Figure 9b, it is found that the values of the C r becomes lower as the values of incident angle increase whilst, the transmission coefficients decrease for higher values of incident angle. This is because as the stiffness becomes higher, the flexible plate becomes stiffer which results in a higher reflection (see [18]). Furthermore, the observations in Figure 9c are similar to Figure 9a. From Figure 9d, it is observed that the reflection coefficient increases with an increase in flexural rigidity whilst, the trend of transmission coefficient becomes the reverse to that of reflection coefficients.

Conclusions
In the present paper, the new contribution is the addition of mooring edge boundary conditions to the BIEM formulation and the effect of mooring stiffness on the hydroelastic response of the floating plate. The BIEM code is compared with existing published results and experimental datasets from the other calculations as well as the SWA model. Furthermore, the hydroelastic response of moored finite floating flexible plate is studied by analyzing the effect of mooring stiffness, incident angle, and flexural rigidity versus different physical parameters on the deflection amplitude, the plate displacement, and the reflection and the transmission coefficients. It has been concluded that:

•
The mathematical model of the moored floating flexible plate using the BIEM technique and the analysis of its hydroelastic response is reported for the first time here.

•
Comparison of the results reveal that the present BIEM model results are supported by existing published results from other calculations and experimental datasets as well as the SWA model.

•
The deflection amplitudes of the moored floating flexible plate decrease for higher values of mooring stiffness and flexural rigidity which is because as these two parameters go on increasing the flexible plate becomes stiffer that leads to less deflection. Figure 9. Effects of (a) mooring stiffness versus non-dimensional wavelength, (b) incidence angle, (c) mooring stiffness, and (d) flexural rigidity versus wave period on the C r and C t for EI = 10 5 , = 1.5 EIρg, and b/h = 20.

Conclusions
In the present paper, the new contribution is the addition of mooring edge boundary conditions to the BIEM formulation and the effect of mooring stiffness on the hydroelastic response of the floating plate. The BIEM code is compared with existing published results and experimental datasets from the other calculations as well as the SWA model. Furthermore, the hydroelastic response of moored finite floating flexible plate is studied by analyzing the effect of mooring stiffness, incident angle, and flexural rigidity versus different physical parameters on the deflection amplitude, the plate displacement, and the reflection and the transmission coefficients. It has been concluded that:

•
The mathematical model of the moored floating flexible plate using the BIEM technique and the analysis of its hydroelastic response is reported for the first time here. • Comparison of the results reveal that the present BIEM model results are supported by existing published results from other calculations and experimental datasets as well as the SWA model.

•
The deflection amplitudes of the moored floating flexible plate decrease for higher values of mooring stiffness and flexural rigidity which is because as these two parameters go on increasing the flexible plate becomes stiffer that leads to less deflection.

•
The shear force and bending moment are decreased in nature when there are higher values of mooring stiffness and incident angle. Furthermotr, the reflection coeffi-cients increasing in nature as the mooring stiffness increases whilst the transmission coefficients decrease.

•
The hydroelastic response analysis of different physical parameters suggested that, for higher spring stiffness, the deflection pattern becomes more regular whilst rigidity in the flexible plate creates fewer deformations for higher values of flexural rigidity. This concludes that a suitable choice of springs and flexural rigidity play a vital role in modelling a wave energy conversion device.
As a future step, the present formulation and the integro-differential equation can be generalized for problems arising in civil engineering application to breakwater and WECs.