Floating Elastic Plate

Introduction

The floating elastic plate is one of the best studied problems in hydroelasticity. It can be used to model a range of physical structures such as a floating break water, an ice floe or a VLFS). The equations of motion were formulated more than 100 years ago and a discussion of the problem appears in Stoker 1957. The problem can be divided into the two and three dimensional formulations which are closely related.

Equations of Motion

For a thin plate floating on the surface of the water, the motion of the plate is derived from the Classical Thin Plate Theory which was developed by Kirchhoff. The equation of motion for a thin plate (also known as Kirchhoff's equation) is given by the following

[math]\displaystyle{ D\frac{\partial^4 \eta}{\partial x^4} + \rho_i h \frac{\partial^2 \eta}{\partial t^2} = p }[/math]

where [math]\displaystyle{ D }[/math] is the flexural rigidity, [math]\displaystyle{ \rho_i }[/math] is the density of the plate, [math]\displaystyle{ h }[/math] is the thickness of the plate (assumed constant), [math]\displaystyle{ p }[/math] is the pressure and [math]\displaystyle{ \eta }[/math] is the plate displacement.

The pressure is given by the linearised Bernoulli equation at the wetted surface (assuming zero pressure at the surface), i.e.

[math]\displaystyle{ p = \rho g \frac{\partial \phi}{\partial z} + \rho \frac{\partial \phi}{\partial t} }[/math]

where [math]\displaystyle{ \rho }[/math] is the water density and [math]\displaystyle{ g }[/math] is gravity, and [math]\displaystyle{ \phi }[/math] is the velocity potential. The velocity potential is governed by Laplace's equation through out the fluid domain subject to the free surface condition and the condition of no flow through the bottom surface. If we denote the region of the fluid surface covered in the plate (or possible multiple plates) by [math]\displaystyle{ P }[/math] and the free surface by [math]\displaystyle{ F }[/math] the equations of motion for the Frequency Domain Problem with frequency [math]\displaystyle{ \omega }[/math] for water of Finite Depth are the following. At the surface we have the dynamic condition

[math]\displaystyle{ D\frac{\partial^4 \eta}{\partial x^4} +\left(\rho g- \omega^2 \rho_i h \right)\eta = i\omega \rho \phi, \, z=0, \, x\in P }[/math]

[math]\displaystyle{ 0= \rho g \frac{\partial \phi}{\partial z} + i\omega \rho \phi, \, x\in F }[/math]

and the kinematic condition

[math]\displaystyle{ \frac{\partial\phi}{\partial z} = i\omega\eta }[/math]

the equation within the fluid is Laplace's Equation

[math]\displaystyle{ \nabla^2\phi =0 }[/math]

and we have the no-flow condition through the bottom boundary

[math]\displaystyle{ \frac{\partial \phi}{\partial z} = 0, \, z=-h }[/math]

(so we have a fluid of constant depth with the bottom surface at [math]\displaystyle{ z=-h }[/math] and the free surface or plate covered surface are at [math]\displaystyle{ z=0 }[/math]). [math]\displaystyle{ g }[/math] is the acceleration due to gravity, [math]\displaystyle{ \rho_i }[/math] and [math]\displaystyle{ \rho }[/math] are the densities of the plate and the water respectively, [math]\displaystyle{ h }[/math] and [math]\displaystyle{ D }[/math] the thickness and flexural rigidity of the plate.

Finally we need to include some boundary conditions at the edge of the plate. The most common boundary conditions in pratical applications are that the edges are free, this means that we have the additional conditions that

[math]\displaystyle{ \frac{\partial^2 \eta}{\partial x^2} = 0, \,\,\frac{\partial^3 \eta}{\partial x^3} = 0 }[/math]

at the edges of the plate.

Solution Method

There are many different methods to solve the corresponding equations ranging from highly analytic such as the Wiener-Hopf to very numerical based on Eigenfunction Matching Method which are applicable and have advantages in different situations.