Difference between revisions of "Category:Floating Elastic Plate"

From WikiWaves
Jump to navigationJump to search
 
(34 intermediate revisions by 2 users not shown)
Line 1: Line 1:
= Introduction =
+
{{complete pages}}
 +
 
 +
== Introduction ==
  
 
The floating elastic plate is one of the best studied problems in hydroelasticity. It can be used to model a range of
 
The floating elastic plate is one of the best studied problems in hydroelasticity. It can be used to model a range of
Line 13: Line 15:
 
We present here the theory of an elastic plate in a vacuum, concentrating of the two-dimensional problem.
 
We present here the theory of an elastic plate in a vacuum, concentrating of the two-dimensional problem.
  
For a [http://en.wikipedia.org/wiki/Euler_Bernoulli_beam_equation Bernoulli-Euler Beam] on the surface of the water, the equation of motion is given
+
{{equations for a uniform free beam}}
by the following
+
 
<center><math>
+
If we assume that the pressure is of the form <math>p(x,t) = e^{i\omega t} \bar{p}(x)</math> then it follows that
D\frac{\partial^4 \eta}{\partial x^4} + \rho_i h \frac{\partial^2 \eta}{\partial t^2} = p
+
<math>\zeta(x,t) = e^{i\omega t} \bar{\zeta}(x) </math> from linearity. In this case the equations reduce to
</math></center>
 
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 can satisfy a range of boundary conditions. The natural boundary condition (i.e. free-edge boundary conditions).
 
 
<center><math>
 
<center><math>
\frac{\partial^2 \eta}{\partial x^2} = 0, \,\,\frac{\partial^3 \eta}{\partial x^3} = 0
+
\partial_x^2\left(\beta \partial_x^2 \bar{\zeta}\right)  -\omega^2 \rho_i h \bar{\zeta} = \bar{p}
 
</math></center>
 
</math></center>
at the edges of the plate.
 
  
If we assume that the pressure is of the form <math>p(x,t) = e^{i\omega t} \bar{p}(x)</math> then it follows that
+
== Linear Elastic Plate on Water ==
<math>\eta(x,t) = e^{i\omega t} \bar{\eta}(x) </math> from linearity. In this case the equations reduce to
 
  
<center><math>
+
We derive here the equations for a plate on a fluid, ignoring boundary conditions at the plate edge and assuming
D\frac{\partial^4 \bar{\eta}}{\partial x^4} -\omega^2 \rho_i h \bar{\eta} = \bar{p}
+
the plate occupies the entire fluid region.
</math></center>
 
  
== Linear Elastic Thin Plate on Water ==
+
{{linear elastic plate on water time domain}}
  
{{linear equations for finite depth in time domain}}
+
== Frequency Domain Equations ==
  
<center>
+
{{frequency domain equations for a floating plate}}
<math>D\frac{\partial^4}{\partial x^4} \frac{\partial \phi}{\partial z}
 
+ \left(\rho g
 
- \omega^2 \rho_i d \right) \frac{\partial \phi}{\partial z} =
 
- \rho \omega^2 \phi, \, z=0</math>
 
</center>
 
plus the equations within  the fluid
 
<center>
 
<math>\nabla^2\phi =0 </math>
 
</center>
 
<center>
 
<math>\frac{\partial \phi}{\partial z} = 0, \, z=-h</math></center>
 
where <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> \,d </math> and <math> \,D </math> are
 
the thickness and flexural rigidity of the plate.
 
  
 +
== Nonlinear Elastic Thin Plate on Water ==
  
<center><math> \frac{\partial\zeta_1}{\partial t} = \frac{\partial\Phi_1}{\partial z} , \ z=0;  </math></center>
+
The nonlinear plate equations can be found by include
which follows from the Kinematic equation and
+
We begin with the [[Conservation Laws and Boundary Conditions | nonlinear equations]] for a fluid.
<center><math> \zeta_1 = -\frac{1}{g} \frac{\partial\Phi_1}{\partial t}, \ z=0; </math></center>
+
The kinematic condition is
which follows from the Dynamic equation. These are the linear free surface conditions.  
+
<center><math>  
 +
\frac{\partial\zeta}{\partial t}+\frac{\partial\Phi}{\partial x} \frac{\partial\zeta}{\partial x}
 +
=\frac{\partial\Phi}{\partial z},  
 +
\ z=\zeta(x,t) </math></center>
 +
the dynamic condition is
 +
<center><math>  
 +
\partial_x^2 \left( \frac{ D \partial_x^2 \zeta}{(1+(\partial_x\zeta^2))^{3/2}} \right)
 +
+ \rho_i h \frac{\partial^2 \zeta}{\partial t^2}
 +
= - \rho\frac{\partial\Phi}{\partial t}-\rho\frac{1}{2}\nabla\Phi \cdot \nabla \Phi - \rho g \zeta, \ z=\zeta(x,t)
 +
</math></center>
 +
There are different versions of the nonlinear boundary condition for a plate and this one is based on
 +
[[Parau and Dias 2002]].
 +
We also have Laplace's equation
 +
<center><math>
 +
\Delta \Phi = 0,\,\,-h<z<\zeta,
 +
</math></center>
 +
and the usual non-flow condition at the bottom surface
 +
<center><math>
 +
\partial_z \Phi = 0,\,\,z=-h.
 +
</math></center>
  
 
[[Category:Linear Hydroelasticity]]
 
[[Category:Linear Hydroelasticity]]

Latest revision as of 09:39, 20 October 2009


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. The plate is assumed to be isotropic while the water motion is irrotational and inviscid.

The solution methods are divided up into those for Two-Dimensional Floating Elastic Plate and those for a Three-Dimensional Floating Elastic Plate.

Linear Thin Elastic Plate Theory

We present here the theory of an elastic plate in a vacuum, concentrating of the two-dimensional problem.

There are various beam theories that can be used to describe the motion of the beam. The simplest theory is the Bernoulli-Euler Beam theory (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, the equation of motion is given by the following

[math]\displaystyle{ \partial_x^2\left(\beta(x)\partial_x^2 \zeta\right) + \gamma(x) \partial_t^2 \zeta = p }[/math]

where [math]\displaystyle{ \beta(x) }[/math] is the non dimensionalised flexural rigidity, and [math]\displaystyle{ \gamma }[/math] is non-dimensionalised linear mass density function. Note that this equations simplifies if the plate has constant properties (and that [math]\displaystyle{ h }[/math] is the thickness of the plate, [math]\displaystyle{ p }[/math] is the pressure and [math]\displaystyle{ \zeta }[/math] is the plate vertical displacement) .

The edges of the plate can satisfy a range of boundary conditions. The natural boundary condition (i.e. free-edge boundary conditions).

[math]\displaystyle{ \partial_x^2 \zeta = 0, \,\,\partial_x^3 \zeta = 0 }[/math]

at the edges of the plate.

The problem is subject to the initial conditions

[math]\displaystyle{ \zeta(x,0)=f(x) \,\! }[/math]
[math]\displaystyle{ \partial_t \zeta(x,0)=g(x) }[/math]

If we assume that the pressure is of the form [math]\displaystyle{ p(x,t) = e^{i\omega t} \bar{p}(x) }[/math] then it follows that [math]\displaystyle{ \zeta(x,t) = e^{i\omega t} \bar{\zeta}(x) }[/math] from linearity. In this case the equations reduce to

[math]\displaystyle{ \partial_x^2\left(\beta \partial_x^2 \bar{\zeta}\right) -\omega^2 \rho_i h \bar{\zeta} = \bar{p} }[/math]

Linear Elastic Plate on Water

We derive here the equations for a plate on a fluid, ignoring boundary conditions at the plate edge and assuming the plate occupies the entire fluid region.

We begin with the linear equations for a fluid. The kinematic condition is the same

[math]\displaystyle{ \frac{\partial\zeta}{\partial t} = \frac{\partial\Phi}{\partial z} , \ z=0; }[/math]

but the dynamic condition needs to be modified to include the effect of the the plate

[math]\displaystyle{ \rho g\zeta + \rho \frac{\partial\Phi}{\partial t} = D \frac{\partial^4 \eta}{\partial x^4} + \rho_i h \frac{\partial^2 \eta}{\partial t^2} , \ z=0; }[/math]

We also have Laplace's equation

[math]\displaystyle{ \Delta \Phi = 0,\,\,-h\lt z\lt 0 }[/math]

and the usual non-flow condition at the bottom surface

[math]\displaystyle{ \partial_z \Phi = 0,\,\,z=-h, }[/math]

where [math]\displaystyle{ \zeta }[/math] is the surface displacement, [math]\displaystyle{ \Phi }[/math] is the velocity potential, and [math]\displaystyle{ \rho }[/math] is the fluid density.

Frequency Domain Equations

If we make the assumption of Frequency Domain Problem that everything is proportional to [math]\displaystyle{ \exp (-\mathrm{i}\omega t)\, }[/math] the equations become

[math]\displaystyle{ \begin{align} -\mathrm{i}\omega\zeta &= \partial_z\phi , &z=0 \\ \rho g\zeta - \mathrm{i}\omega\rho \phi &= D \partial_x^4 \eta -\omega^2 \rho_i h \zeta, &z=0 \\ \Delta \phi &= 0, &-h\lt z\lt 0 \\ \partial_z \phi &= 0, &z=-h, \end{align} }[/math]

where [math]\displaystyle{ \zeta }[/math] is the surface displacement and [math]\displaystyle{ \phi }[/math] is the velocity potential in the frequency domain.

These equations can be simplified by defining [math]\displaystyle{ \alpha = \omega^2/g }[/math], [math]\displaystyle{ \beta = D/\rho g }[/math] and [math]\displaystyle{ \gamma = \rho_i h/\rho }[/math] to obtain

[math]\displaystyle{ \begin{align} \Delta \phi &= 0, &-h \lt z \leq 0 \\ \partial_z \phi &= 0, &z = - h \\ \beta \partial_x^4 \zeta + \left( 1 - \gamma\alpha \right) \zeta &= -\mathrm{i} \sqrt{\alpha}\phi, &z = 0 \\ -\mathrm{i}\omega\zeta &= \partial_z\phi , &z=0 . \end{align} }[/math]

Nonlinear Elastic Thin Plate on Water

The nonlinear plate equations can be found by include We begin with the nonlinear equations for a fluid. The kinematic condition is

[math]\displaystyle{ \frac{\partial\zeta}{\partial t}+\frac{\partial\Phi}{\partial x} \frac{\partial\zeta}{\partial x} =\frac{\partial\Phi}{\partial z}, \ z=\zeta(x,t) }[/math]

the dynamic condition is

[math]\displaystyle{ \partial_x^2 \left( \frac{ D \partial_x^2 \zeta}{(1+(\partial_x\zeta^2))^{3/2}} \right) + \rho_i h \frac{\partial^2 \zeta}{\partial t^2} = - \rho\frac{\partial\Phi}{\partial t}-\rho\frac{1}{2}\nabla\Phi \cdot \nabla \Phi - \rho g \zeta, \ z=\zeta(x,t) }[/math]

There are different versions of the nonlinear boundary condition for a plate and this one is based on Parau and Dias 2002. We also have Laplace's equation

[math]\displaystyle{ \Delta \Phi = 0,\,\,-h\lt z\lt \zeta, }[/math]

and the usual non-flow condition at the bottom surface

[math]\displaystyle{ \partial_z \Phi = 0,\,\,z=-h. }[/math]