Difference between revisions of "Two-Dimensional Floating Elastic Plate"

From WikiWaves
Jump to navigationJump to search
 
(12 intermediate revisions by the same user not shown)
Line 1: Line 1:
= Equations of Motion =
+
== 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 [http://en.wikipedia.org/wiki/Euler_Bernoulli_beam_equation 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.
+
{{equations for a beam}}
  
 +
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.
  
For a Bernoulli-Euler beam on the surface of the water, the equation of motion is given
+
{{linear elastic plate on water time domain}}
by the following
 
<center><math>
 
D\frac{\partial^4 \eta}{\partial x^4} + \rho_i h \frac{\partial^2 \eta}{\partial t^2} = p
 
</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 satisfy the natural boundary condition (i.e. free-edge boundary conditions).
+
=== Frequency Domain Equations ===
<center><math>
 
\frac{\partial^2 \eta}{\partial x^2} = 0, \,\,\frac{\partial^3 \eta}{\partial x^3} = 0
 
</math></center>
 
at the edges of the plate.
 
  
The pressure is given by the linearised Bernoulli equation at the wetted surface (assuming zero
+
{{frequency domain equations for a floating plate}}
pressure at the surface), i.e.
 
<center>
 
<math>p = \rho g \frac{\partial \phi}{\partial z} + \rho \frac{\partial \phi}{\partial t}</math>
 
</center>
 
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 beam-covered (or possible multiple beams covered) region of the fluid 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
 
<center>
 
<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>
 
</center>
 
<center><math>0=
 
\rho g \frac{\partial \phi}{\partial z} + i\omega \rho \phi, \, x\in F</math>
 
</center>
 
and the kinematic condition
 
<center>
 
<math>\frac{\partial\phi}{\partial z} = i\omega\eta</math>
 
</center>
 
  
 +
== Computational Methods ==
  
The equation within the fluid is governed by [[Laplace's Equation]]  
+
This wiki contains a number of methods to solve for a floating elastic plate in the frequency domain.
<center>
+
* [[Eigenfunction Matching for a Semi-Infinite Floating Elastic Plate]]
<math>\nabla^2\phi =0 </math>
+
* [[Eigenfunction Matching for a Finite Floating Elastic Plate using Symmetry]]
</center>
+
* [[Eigenfunction Matching Method for Floating Elastic Plates]]
and we have the no-flow condition through the bottom boundary
+
* [[Eigenfunction Matching for a Circular Floating Elastic Plate]]
<center>
+
* [[Green Function Methods for Floating Elastic Plates]]
<math>\frac{\partial \phi}{\partial z} = 0, \, z=-h</math>
 
</center>
 
(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.
 
  
=Energy Balance=
+
== Solution Methods ==
An energy balence relation is derived in
 
[[Evans and Davies 1968]] which is simply a condition that the incident energy is equal to the sum of the radiated energy including both the energy
 
in the water and the energy in the plate.
 
If the properties of the first and last semi-infinite plates were identical, then
 
this would be the familiar requirement that
 
<center><math>  |T_{\Lambda}(0)|^2 + |R_{1}(0)|^2 = |I|^2. </math></center>
 
However, when the first and last plates have different properties, then the
 
energy balance condition becomes the following
 
<center><math>
 
    D|T_{\Lambda}(0)|^2 + |R_{1}(0)|^2 = |I|^2,
 
</math></center>
 
where <math>D</math> is found by applying Green's theorem to <math>\phi</math> and its conjugate,
 
[[Evans and Davies 1968]]
 
and is given by
 
  <center><math>  \begin{matrix}    D & = &  {  \frac{\kappa_{\Lambda}(0)k_1(0)\cosh^2{(k_1(0)h)}}{\kappa_{1}(0)k_\Lambda(0)\cosh^2{(k_\Lambda(0)h)}} \times }\\        &  &  {  \frac{\left(\frac{\beta_\Lambda}{\alpha}4k_\Lambda(0)^3(\kappa_{\Lambda}(0)^2 +k_y^2)\sinh^2{(k_\Lambda(0) h)} +                    \frac{1}{2}{\sinh{(2k_\Lambda(0)h)}}+k_\Lambda(0)h\right)}            {\left(\frac{\beta_1}{\alpha}4k_1(0)^3(\kappa_{1}(0)^2 +k_y^2)\sinh^2{(k_1(0)h)} +                    \frac{1}{2}{\sinh{(2k_1(0)h)}}+k_1(0)h\right)}.  }  \end{matrix}   </math></center> 
 
The energy balance condition is useful to help check that the
 
solution is not incorrect (it does not of course guarantee the solution
 
is correct). The energy balance condition is surprisingly well
 
satisfied by our solutions, for example with <math>M=20</math> we can easily get ten decimal places.
 
 
 
 
 
 
 
= Solution Methods =
 
  
 
There are many different methods to solve the corresponding equations ranging from highly analytic such
 
There are many different methods to solve the corresponding equations ranging from highly analytic such
as the [[:Category:Wiener-Hopf|Wiener-Hopf]] to very numerical based on [[Eigenfunction Matching Method]] which are
+
as the [[:Category:Wiener-Hopf|Wiener-Hopf]] to very numerical based on the [[:Category:Eigenfunction Matching Method|Eigenfunction Matching Method]] which are
 
applicable and have advantages in different situations. We describe here some of the solutions
 
applicable and have advantages in different situations. We describe here some of the solutions
 
which have been developed, grouped by problem
 
which have been developed, grouped by problem
Line 114: Line 54:
 
described in [[Eigenfunction Matching Method for Floating Elastic Plates]]
 
described in [[Eigenfunction Matching Method for Floating Elastic Plates]]
 
The eigenfunction matching method using  
 
The eigenfunction matching method using  
incomplete mode-matching was used by [[Hendrickson and Webb 1963]], [[Wadhams 1973]],[[Squire 1978]], [[Squire 1984]], [[Wadhams 1986]]. [[Fox and Squire 1990]] computed the solution for the full set of eigenfunctions using a conjugate gradient technique. [[Fox and Squire 1994]] used the same method to complete further studies on this problem, investigating the strain in the ice, and also the effect of shore fast ice on an incoming directional wave spectrum of specified structure.
+
incomplete mode-matching was used by [[Hendrickson and Webb 1963]], [[Wadhams 1973]], [[Squire 1978]], [[Squire 1984]], [[Wadhams 1986]]. [[Fox and Squire 1990]] computed the solution for the full set of eigenfunctions using a conjugate gradient technique. [[Fox and Squire 1994]] used the same method to complete further studies on this problem, investigating the strain in the ice, and also the effect of shore fast ice on an incoming directional wave spectrum of specified structure.
  
 
Other authors have extended mode-matching approach of [[Fox and Squire 1990]]. The eigenfunction matching solution was extended  
 
Other authors have extended mode-matching approach of [[Fox and Squire 1990]]. The eigenfunction matching solution was extended  
to plates of different properties by [[Barrett and Squire 1996]]. [[Sahoo et al. 2001]] defined an inner product enabling the solution to be found by using <math>N</math> eigenfunctions and inverting an <math>N\times N</math> matrix, Earlier [[Chakrabarti 2000]] had used an infinite depth mode-matching scheme to set up a singular integral equation, the equivalent problem to a residue calculus problem when the eigenvalue spectrum is continuous. This singular integral equation was solved by transforming it into a Hilbert problem ([[Roos 1969]]), a generalization of the [[Wiener-Hopf]] problem (although in practice both are solved in exactly the same way).
+
to plates of different properties by [[Barrett and Squire 1996]]. [[Sahoo et al. 2001]] defined an inner product enabling the solution to be found by using <math>N</math> eigenfunctions and inverting an <math>N\times N</math> matrix, Earlier [[Chakrabarti 2000]] had used an infinite depth mode-matching scheme to set up a singular integral equation, the equivalent problem to a residue calculus problem when the eigenvalue spectrum is continuous. This singular integral equation was solved by transforming it into a Hilbert problem ([[Roos 1969]]), a generalization of the [[:Category:Wiener-Hopf|Wiener-Hopf]] problem (although in practice both are solved in exactly the same way).
  
 
===[[:Category:Residue Calculus|Residue Calculus]]===
 
===[[:Category:Residue Calculus|Residue Calculus]]===
Line 137: Line 77:
  
 
[[Category:Floating Elastic Plate]]
 
[[Category:Floating Elastic Plate]]
 +
[[Category:Complete Pages]]

Latest revision as of 07:15, 4 April 2009

Equations of Motion

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]

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]

Computational Methods

This wiki contains a number of methods to solve for a floating elastic plate in the frequency domain.

Solution Methods

There are many different methods to solve the corresponding equations ranging from highly analytic such as the Wiener-Hopf to very numerical based on the Eigenfunction Matching Method which are applicable and have advantages in different situations. We describe here some of the solutions which have been developed, grouped by problem

Two Semi-Infinite Plates of Identical Properties

The simplest problem to consider is one where there are only two semi-infinite plates of identical properties separated by a crack. A related problem in acoustics was considered by Kouzov 1963 who used an integral representation of the problem to solve it explicitly using the Riemann-Hilbert technique. Recently the crack problem has been considered by Squire and Dixon 2000 and Williams and Squire 2002 using a Green function method applicable to infinitely deep water and they obtained simple expressions for the reflection and transmission coefficients. Squire and Dixon 2001 extended the single crack problem to a multiple crack problem in which the semi-infinite regions are separated by a region consisting of a finite number of plates of finite size with all plates having identical properties. Evans and Porter 2005 further considered the multiple crack problem for finitely deep water and provided an explicit solution which is described in Two Semi-Infinite Elastic Plates of Identical Properties

Two Semi-Infinite Plates of Different Properties

The next most simple problem is two semi-infinite plates of different properties. Often one of the plates is taken to be open water which makes the problem simpler. In general, the solution method developed for open water can be extended to two plates of different properties, the exception to this is the Residue Calculus solution which applies only when one of the semi-infinite regions is water.

Wiener-Hopf

The solution to the problem of two semi-infinite plates with different properties can be solved by the Wiener-Hopf method. The first work on this problem was by Evans and Davies 1968 but they did not actually develop the method sufficiently to be able to calculate the solution. At the time the infinite products that formed a large part of the solution proved to be too difficult to compute and so they were not able to present any results, except for certain limiting cases (such as the shallow water limit). This was because at the time the infinite products that formed a large part of the solution proved to be too difficult to compute and so they were not able to present general results. Gol'dshtein and Marchenko 1989 presented a later Wiener-Hopf solution for infinite depth, but again only certain limiting situations were discussed. However, in the late 1990's and early [math]\displaystyle{ 21^{st} }[/math] century improvements in computing power enabled other authors to compute results using the Wiener-Hopf solution of Evans and Davies 1968. Balmforth and Craster 1999 turned the required infinite products into integrals which were evaluated by quadrature, while Chung and Fox 2002 showed that the products themselves could be evaluated directly with relatively little effort. Ironically Tkacheva 2001 finally showed that if the inertia term in the thin plate equation was neglected for normally incident waves (as it can for most wavelengths), then [math]\displaystyle{ |R| }[/math] could be calculated by simply using the correct value for [math]\displaystyle{ k_0 }[/math], the wavenumber in the ice, in the formula for Mass Loading Model of Ice given by Keller and Weitz 1953. The Wiener-Hopf solution is described in Wiener-Hopf Elastic Plate Solution

Eigenfunction Matching Method

Essentially the solution is expanded on either side of the crack using separation of variables. The theory is described in Eigenfunction Matching Method for Floating Elastic Plates The eigenfunction matching method using incomplete mode-matching was used by Hendrickson and Webb 1963, Wadhams 1973, Squire 1978, Squire 1984, Wadhams 1986. Fox and Squire 1990 computed the solution for the full set of eigenfunctions using a conjugate gradient technique. Fox and Squire 1994 used the same method to complete further studies on this problem, investigating the strain in the ice, and also the effect of shore fast ice on an incoming directional wave spectrum of specified structure.

Other authors have extended mode-matching approach of Fox and Squire 1990. The eigenfunction matching solution was extended to plates of different properties by Barrett and Squire 1996. Sahoo et al. 2001 defined an inner product enabling the solution to be found by using [math]\displaystyle{ N }[/math] eigenfunctions and inverting an [math]\displaystyle{ N\times N }[/math] matrix, Earlier Chakrabarti 2000 had used an infinite depth mode-matching scheme to set up a singular integral equation, the equivalent problem to a residue calculus problem when the eigenvalue spectrum is continuous. This singular integral equation was solved by transforming it into a Hilbert problem (Roos 1969), a generalization of the Wiener-Hopf problem (although in practice both are solved in exactly the same way).

Residue Calculus

The solution using Residue Calculus was developed Linton and Chung 2003 who effectively showed that the equations Sahoo et al. 2001 had set up could be solved analytically using Residue Calculus. In the process they also demonstrated its equivalence to the Wiener-Hopf solution, and confirmed the formula of Tkacheva 2001.

Single Floating Plate

The problem of a single floating plate in two-dimensions was treated by Newman 1994, Meylan and Squire 1994 and Hermans 2003 using the Free-Surface Green Function, described in Green Function Methods for Floating Elastic Plates

Multiple Floating Plates

The most general problem consists of multiple floating plate. The methods which generalises to this are the Eigenfunction Matching Method for Floating Elastic Plates (Kohout et. al. 2006) and the Green Function Methods for Floating Elastic Plates (Hermans 2004).