Difference between revisions of "Two Semi-Infinite Elastic Plates of Identical Properties"

From WikiWaves
Jump to navigationJump to search
Line 1: Line 1:
 +
= Introduction =
 +
 
We present here the solution of [[Evans and Porter 2005]] for the simple
 
We present here the solution of [[Evans and Porter 2005]] for the simple
 
case of a single crack with waves incident from normal (they also considered multiple cracks
 
case of a single crack with waves incident from normal (they also considered multiple cracks
Line 5: Line 7:
 
<math>\phi</math> in terms of a linear combination of the incident wave and certain source functions located at the crack.
 
<math>\phi</math> in terms of a linear combination of the incident wave and certain source functions located at the crack.
 
Along with satisfying the field and boundary conditions, these source functions satisfy the jump conditions in the displacements and gradients across the crack.
 
Along with satisfying the field and boundary conditions, these source functions satisfy the jump conditions in the displacements and gradients across the crack.
 +
 +
= Governing Equations =
 +
 +
We consider the entire free surface to be occupied by a [[Floating Elastic Plate]] with a single discontinuity at <math>z=0</math>.
 +
The equations are the following
 +
<center><math>
 +
\nabla^2 \phi = 0, -H<z<0,
 +
</math></center>
 +
<center><math>
 +
\frac{\partial \phi}{\partial z} =0, z=-H,
 +
</math></center>
 +
<center><math>
 +
{\left( \beta \frac{\partial^4}{\partial x^4} -
 +
\gamma\alpha + 1\right)\frac{\partial \phi}{\partial z} - \alpha \phi}
 +
= 0, z=0, x\neq 0.
 +
</math></center>
 +
 +
 
They first define <math>\chi(x,z)</math> to be the Two-Dimensional solution to the [[Free-Surface Green Function for a Floating Elastic Plate]]
 
They first define <math>\chi(x,z)</math> to be the Two-Dimensional solution to the [[Free-Surface Green Function for a Floating Elastic Plate]]
 
given by
 
given by
Line 22: Line 42:
 
-\alpha \cos(kH). </math>
 
-\alpha \cos(kH). </math>
 
</center>
 
</center>
 +
This satisfies the following equation
 +
The Green function <math>G(x,z)</math> representing
 +
outgoing waves as <math>|x|\rightarrow \infty</math> satisfies
 +
<center><math>
 +
\nabla^2 G = 0, -H<z<0,
 +
</math></center>
 +
<center><math>
 +
\frac{\partial G}{\partial z} =0, z=-H,
 +
</math></center>
 +
<center><math>
 +
{\left( \beta \frac{\partial^4}{\partial x^4} -
 +
\gamma\alpha + 1\right)\frac{\partial G}{\partial z} - \alpha G}
 +
= \delta(x), z=0,
 +
</math></center>
 
Consequently, the source functions for a single crack at <math>x=0</math> can be defined as
 
Consequently, the source functions for a single crack at <math>x=0</math> can be defined as
 
<center>
 
<center>

Revision as of 22:34, 5 April 2007

Introduction

We present here the solution of Evans and Porter 2005 for the simple case of a single crack with waves incident from normal (they also considered multiple cracks and waves incident from different angles). The solution of Evans and Porter 2005 expresses the potential [math]\displaystyle{ \phi }[/math] in terms of a linear combination of the incident wave and certain source functions located at the crack. Along with satisfying the field and boundary conditions, these source functions satisfy the jump conditions in the displacements and gradients across the crack.

Governing Equations

We consider the entire free surface to be occupied by a Floating Elastic Plate with a single discontinuity at [math]\displaystyle{ z=0 }[/math]. The equations are the following

[math]\displaystyle{ \nabla^2 \phi = 0, -H\lt z\lt 0, }[/math]
[math]\displaystyle{ \frac{\partial \phi}{\partial z} =0, z=-H, }[/math]
[math]\displaystyle{ {\left( \beta \frac{\partial^4}{\partial x^4} - \gamma\alpha + 1\right)\frac{\partial \phi}{\partial z} - \alpha \phi} = 0, z=0, x\neq 0. }[/math]


They first define [math]\displaystyle{ \chi(x,z) }[/math] to be the Two-Dimensional solution to the Free-Surface Green Function for a Floating Elastic Plate given by

[math]\displaystyle{ \chi(x,z) = -i\sum_{n=-2}^\infty\frac{\sin{(k_n h)}\cos{(k_n(z-h))}}{2\alpha C_n}e^{-k_n|x|},\,\,\,(1) }[/math]

where

[math]\displaystyle{ C_n=\frac{1}{2}\left(h + \frac{(5\beta k_n ^4 + 1 - \alpha\gamma)\sin^2{(k_n h)}}{\alpha}\right), }[/math]

and [math]\displaystyle{ k_n }[/math] are the solutions of the Dispersion Relation for a Floating Elastic Plate.

[math]\displaystyle{ \beta k^5 \sin(kH) + k \left(1 - \alpha \gamma \right) \sin(kH) = -\alpha \cos(kH). }[/math]

This satisfies the following equation The Green function [math]\displaystyle{ G(x,z) }[/math] representing outgoing waves as [math]\displaystyle{ |x|\rightarrow \infty }[/math] satisfies

[math]\displaystyle{ \nabla^2 G = 0, -H\lt z\lt 0, }[/math]
[math]\displaystyle{ \frac{\partial G}{\partial z} =0, z=-H, }[/math]
[math]\displaystyle{ {\left( \beta \frac{\partial^4}{\partial x^4} - \gamma\alpha + 1\right)\frac{\partial G}{\partial z} - \alpha G} = \delta(x), z=0, }[/math]

Consequently, the source functions for a single crack at [math]\displaystyle{ x=0 }[/math] can be defined as

[math]\displaystyle{ \psi_s(x,z)= \beta\chi_{xx}(x,z),\,\,\, \psi_a(x,z)= \beta\chi_{xxx}(x,z),\,\,\,(2) }[/math]

It can easily be shown that [math]\displaystyle{ \psi_s }[/math] is symmetric about [math]\displaystyle{ x = 0 }[/math] and [math]\displaystyle{ \psi_a }[/math] is antisymmetric about [math]\displaystyle{ x = 0 }[/math].

Substituting (1) into (2) gives

[math]\displaystyle{ \psi_s(x,z)= { -\frac{\beta}{\alpha} \sum_{n=-2}^\infty \frac{g_n\cos{(k_n(z+h))}}{2k_{xn}C_n}e^{k_n|x|} }, \psi_a(x,z)= { {\rm sgn}(x) i\frac{\beta}{\alpha}\sum_{n=-2}^\infty \frac{g_n'\cos{(k_n(z+h))}}{2k_{xn}C_n}e^{k_n|x|}}, }[/math]

where

[math]\displaystyle{ g_n = ik_n^3 \sin{k_n h},\,\,\,\, g'_n= -k_n^4 \sin{k_n h}. }[/math]

We then express the solution to the problem as a linear combination of the incident wave and pairs of source functions at each crack,

[math]\displaystyle{ \phi(x,z) = e^{-k_0 x}\frac{\cos(k_0(z+h))}{\cos(k_0h)} + (P\psi_s(x,z) + Q\psi_a(x,z))\,\,\,(3) }[/math]

where [math]\displaystyle{ P }[/math] and [math]\displaystyle{ Q }[/math] are coefficients to be solved which represent the jump in the gradient and elevation respectively of the plates across the crack [math]\displaystyle{ x = a_j }[/math]. The coefficients [math]\displaystyle{ P }[/math] and [math]\displaystyle{ Q }[/math] are found by applying the edge conditions and to the [math]\displaystyle{ z }[/math] derivative of [math]\displaystyle{ \phi }[/math] at [math]\displaystyle{ z=0 }[/math],

[math]\displaystyle{ \frac{\partial^2}{\partial x^2}\left. \frac{\partial \phi}{\partial z}\right|_{x=0,z=0}=0,\,\,\, {\rm and}\,\,\,\, \frac{\partial^3}{\partial x^3}\left. \frac{\partial \phi}{\partial z}\right|_{x=0,z=0}=0. }[/math]

The reflection and transmission coefficients, [math]\displaystyle{ R }[/math] and [math]\displaystyle{ T }[/math] can be found from (3) by taking the limits as [math]\displaystyle{ x\rightarrow\pm\infty }[/math] to obtain

[math]\displaystyle{ R = {- \frac{\beta}{\alpha} (g'_0Q + ig_0P)} }[/math]

and

[math]\displaystyle{ T= 1 + {\frac{\beta}{\alpha}(g'_0Q - ig_0P)} }[/math]