This study investigates the hydroelastic interaction of flexural gravity waves with multiple porous elastic plates of varying lengths in finite-depth water, employing an integral equation approach. The floating ice sheet is modelled as a flexible plate of uniform thickness, governed by the Euler–Bernoulli beam equation. The primary objective is to evaluate the effectiveness of porous elastic plates as wave barriers for shoreline protection in ice-covered regions. Within the framework of linearized theory, the problem is formulated as a boundary value problem (BVP) and solved using an eigenfunction expansion method with nonorthogonal eigenfunctions. The mode-coupling relation is utilized to transform the BVP into a system of Fredholm-type integral equations, which is subsequently solved using the multi-term Galerkin approximation technique with Chebyshev polynomials. The numerical analysis evaluates the reflection and transmission coefficients, hydrodynamic forces, and wave energy dissipation, with a particular focus on the influence of the permeability and flexibility of the submerged plates, along with other relevant parameters. Validation is conducted by comparing the results with those of previous studies under specific conditions. This research underscores the practical benefits of incorporating porosity and flexibility into the model, demonstrating improved wave reflection and energy dissipation. Additionally, the findings reveal that the thickness of the ice sheet plays a crucial role in optimizing breakwater performance. The research delivers key insights into mitigating wave-induced forces and offers a reliable framework for designing effective and sustainable coastal protection systems that safeguard shorelines from high waves.