|
|
(24 intermediate revisions by 3 users not shown) |
Line 1: |
Line 1: |
− | = Introduction =
| + | This page has moved to [[:Category:Floating Elastic Plate|Floating Elastic Plate]]. |
| | | |
− | The floating elastic plate is one of the best studied problems in hydroelasticity. It can be used to model a range of
| + | Please change the link to the new page |
− | 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_1957a|Stoker 1957]]. The problem can
| |
− | be divided into the two and three dimensional formulations which are closely related. The plate is assumed to be isotropic while the water motion is irrotational and inviscid.
| |
− | | |
− | = Two Dimensional Problem =
| |
− | | |
− | == Equations of Motion ==
| |
− | | |
− | When considering a two dimensional problem, the <math>y</math> variable is dropped and the plate is regarded as a beam. There are various beam theories that can be used to describe the motion of the beam. The simplest theory is the Bernoulli-Euler beam which is commonly used in the two dimensional hydroelastic analysis. Other beam theories include the Timoshenko beam theory and Reddy-Bickford beam theory where shear deformation of higher order is considered.
| |
− | | |
− | | |
− | For a Bernoulli-Euler beam on the surface of the water, the equation of motion is given
| |
− | by the following
| |
− | | |
− | <math>D\frac{\partial^4 \eta}{\partial x^4} + \rho_i h \frac{\partial^2 \eta}{\partial t^2} = p</math>
| |
− | | |
− | where <math>D</math> is the flexural rigidity, <math>\rho_i</math> is the density of the beam,
| |
− | <math>h</math> is the thickness of the beam (assumed constant), <math> p</math> is the pressure
| |
− | and <math>\eta</math> is the beam vertical displacement.
| |
− | | |
− | The edges of the plate satisfy the natural boundary condition (i.e. free-edge boundary conditions).
| |
− | | |
− | <math>\frac{\partial^2 \eta}{\partial x^2} = 0, \,\,\frac{\partial^3 \eta}{\partial x^3} = 0</math>
| |
− | | |
− | at the edges of the plate.
| |
− | | |
− | The pressure is given by the linearised Bernoulli equation at the wetted surface (assuming zero
| |
− | pressure at the surface), i.e.
| |
− | | |
− | <math>p = \rho g \frac{\partial \phi}{\partial z} + \rho \frac{\partial \phi}{\partial t}</math>
| |
− | | |
− | where <math>\rho</math> is the water density and <math>g</math> is gravity, and <math>\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 beam (or possible
| |
− | multiple beams) by <math>P</math> and the free surface by <math>F</math> the equations of motion for the
| |
− | [[Frequency Domain Problem]] with frequency <math>\omega</math> for water of
| |
− | [[Finite Depth]] are the following. At the surface
| |
− | we have the dynamic condition
| |
− | | |
− | <math>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>0=
| |
− | \rho g \frac{\partial \phi}{\partial z} + i\omega \rho \phi, \, x\in F</math>
| |
− | | |
− | and the kinematic condition
| |
− | | |
− | <math>\frac{\partial\phi}{\partial z} = i\omega\eta</math>
| |
− | | |
− | | |
− | | |
− | The equation within the fluid is governed by [[Laplace's Equation]]
| |
− | | |
− | <math>\nabla^2\phi =0 </math>
| |
− | | |
− | and we have the no-flow condition through the bottom boundary
| |
− | | |
− | <math>\frac{\partial \phi}{\partial z} = 0, \, z=-h</math>
| |
− | | |
− | (so we have a fluid of constant depth with the bottom surface at <math>z=-h</math> and the
| |
− | free surface or plate covered surface are at <math>z=0</math>).
| |
− | <math> g </math> is the acceleration due to gravity, <math> \rho_i </math> and <math> \rho </math>
| |
− | are the densities of the plate and the water respectively, <math> h </math> and <math> D </math>
| |
− | the thickness and flexural rigidity 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.
| |
− | | |
− | = Three Dimensional Problem
| |