Oscillation of an Elastic Sheet

A straightforward generalization of the analysis of Section 4.3 reveals that the transverse oscillation of a uniform elastic sheet, stretched over a rigid frame, is governed by the two-dimensional wave equation (Pain 1999)

$\displaystyle \frac{\partial^{\,2}\psi}{\partial t^{\,2}} = v^{\,2}\left(\frac{...
...\,2}\psi}{\partial x^{\,2}}+\frac{\partial^{\,2}\psi}{\partial y^{\,2}}\right).$ (7.14)

Here, $\psi(x,y,t)$ is the sheet's transverse (i.e., in the $z$-direction) displacement, $v=\sqrt{T/\rho}$ the characteristic speed of elastic waves on the sheet, $T$ the tension, and $\rho$ the mass per unit area. In equilibrium, the sheet is assumed to lie in the $x$-$y$ plane. The boundary condition is that $\psi=0$ at the rigid frame.

Figure 7.4: Density plot illustrating the spatial variation of the $m=2$, $n=3$ normal mode of a rectangular elastic sheet with $a=2\,b$. Dark/light regions indicate positive/negative displacements.
\includegraphics[width=1\textwidth]{Chapter07/fig7_04.eps}

Suppose that the frame is rectangular, extending from $x=0$ to $x=a$, and from $y=0$ to $y=b$. Let us search for a normal mode of the form

$\displaystyle \psi(x,y,t)=\hat{\psi}(x,y)\,\cos(\omega\,t-\phi).$ (7.15)

Substitution into Equation (7.14) yields

$\displaystyle \frac{\partial^{\,2}\hat{\psi}}{\partial x^{\,2}} + \frac{\partia...
...2}\hat{\psi}}{\partial y^{\,2}} + \frac{\omega^{\,2}}{v^{\,2}}\,\hat{\psi} = 0,$ (7.16)

subject to the boundary conditions $\hat{\psi}(0,y)=\hat{\psi}(a,y)=\hat{\psi}(x,0)=\hat{\psi}(x,b)=0$. Let us search for a separable solution of the form

$\displaystyle \hat{\psi}(x,y)= X(x)\,Y(y).$ (7.17)

Such a solution satisfies the boundary conditions provided $X(0)=X(a)=Y(0)=Y(b)=0$. It follows that

$\displaystyle \frac{1}{X}\,\frac{d^{\,2}X}{dx^{\,2}} + \frac{1}{Y}\,\frac{d^{\,2}Y}{dy^{\,2}}+\frac{\omega^{\,2}}{v^{\,2}} = 0.$ (7.18)

The only way that the preceding equation can be satisfied at all $x$ and $y$ is if

$\displaystyle \frac{1}{X}\,\frac{d^{\,2}X}{dx^{\,2}}$ $\displaystyle = -\alpha^{\,2},$ (7.19)
$\displaystyle \frac{1}{Y}\,\frac{d^{\,2}Y}{dy^{\,2}}$ $\displaystyle = -\beta^{\,2},$ (7.20)
$\displaystyle \omega^{\,2}$ $\displaystyle = (\alpha^{\,2}+\beta^{\,2})\,v^{\,2},$ (7.21)

where $\alpha^{\,2}$ and $\beta^{\,2}$ are positive constants. (The constants have to be positive, rather than negative, to give oscillatory solutions that are capable of satisfying the boundary conditions.) Appropriate solutions of Equations (7.19) and (7.20) are

$\displaystyle X(x)$ $\displaystyle = X_0\,\sin(\alpha\,x),$ (7.22)
$\displaystyle Y(y)$ $\displaystyle = Y_0\,\sin(\beta\,y),$ (7.23)

where $X_0$ and $Y_0$ are arbitrary constants. These solutions automatically satisfy the boundary conditions $X(0)=Y(0)=0$. The boundary conditions $X(a)=Y(b)=0$ are satisfied provided

$\displaystyle \alpha$ $\displaystyle = \frac{m\,\pi}{a},$ (7.24)
$\displaystyle \beta$ $\displaystyle = \frac{n\,\pi}{b},$ (7.25)

where $m$ and $n$ are positive integers. Thus, the normal modes of a rectangular elastic sheet, which are indexed by the mode numbers $m$ and $n$, take the form

$\displaystyle \psi_{m,n}(x,y,t) = A_{m,n}\,\sin\left(m\,\pi\,\frac{x}{a}\right)\,\sin\left(n\,\pi\,\frac{y}{b}\right)\,\cos(\omega_{m,n}\,t-\phi_{m,n}),$ (7.26)

where

$\displaystyle \omega_{m,n} = \pi\left(\frac{m^{\,2}}{a^{\,2}}+ \frac{n^{\,2}}{b^{\,2}}\right)^{1/2}v.$ (7.27)

Here, $A_{m,n}$ and $\phi_{m,n}$ are arbitrary constants. Because Equation (7.14) is linear, its solutions are superposable. Hence, the most general solution is a superposition of all of the normal modes; that is,

$\displaystyle \psi(x,y,t)=\sum_{m=1,\infty}\sum_{n=1,\infty}A_{m,n}\,\sin\left(...
...right)\,\sin\left(n\,\pi\,\frac{y}{b}\right)\,\cos(\omega_{m,n}\,t-\phi_{m,n}).$ (7.28)

The amplitudes $A_{m,n}$, and the phase angles $\phi_{n,m}$, are determined by the initial conditions. (See Exercise 4.) Figure 7.4 illustrates the spatial variation of the $m=2$, $n=3$ normal mode of a rectangular elastic sheet with $a=2\,b$.

Figure 7.5: The Bessel function $J_0(z)$.
\includegraphics[width=0.8\textwidth]{Chapter07/fig7_05.eps}

Suppose that an elastic sheet is stretched over a circular frame of radius $a$. Defining the radial cylindrical coordinate $\rho= (x^{\,2}+y^{\,2})^{1/2}$, where $\rho=a$ corresponds to the location of the frame, the axisymmetric oscillations of the sheet are governed by the cylindrical wave equation (see Section 7.4)

$\displaystyle \frac{\partial^{\,2}\psi}{\partial t^{\,2}}=v^{\,2}\left(\frac{\p...
...}{\partial\rho^{\,2}}+\frac{1}{\rho}\,\frac{\partial\psi}{\partial\rho}\right),$ (7.29)

subject to the boundary condition $\psi=0$ at $\rho=a$. Let us search for a normal mode of the form

$\displaystyle \psi(\rho,t) = \hat{\psi}(\rho)\,\cos(\omega\,t-\phi).$ (7.30)

Substitution into Equation (7.29) yields

$\displaystyle \frac{d^{\,2}\hat{\psi}}{d\rho^{\,2}} + \frac{1}{\rho}\,\frac{d\hat{\psi}}{d\rho}+ \frac{\omega^{\,2}}{v^{\,2}}\,\hat{\psi} = 0,$ (7.31)

subject to the boundary condition $\hat{\psi}(a)=0$. Let us define the scaled radial coordinate $z=(\omega/v)\,\rho$. When written in terms of this new coordinate, the previous equation transforms to

$\displaystyle \frac{d^{\,2}\hat{\psi}}{dz^{\,2}} + \frac{1}{z}\,\frac{d\hat{\psi}}{dz}+\hat{\psi}=0.$ (7.32)

This well-known equation is called Bessel's equation (of degree zero), and has the standard solution (that is well-behaved at $z=0$) (Abramowitz and Stegun 1965)

$\displaystyle J_0(z)= \frac{1}{\pi}\int_0^\pi\cos(z\,\sin\nu)\,d\nu.$ (7.33)

Here, $J_0(z)$ is termed a Bessel function (of degree zero), and is plotted in Figure 7.5. It can be seen that the function is broadly similar to a cosine function, except that its zeros are not quite evenly spaced, and its amplitude gradually decreases as $z$ increases. The first few values of $z$ at which $J_0(z)=0$ are listed in Table 7.1. Let the $j$th zero be located at $z=z_j$. In order to satisfy the boundary condition, we require that $(\omega/v)\,a=z_j$. Hence, the axisymmetric normal modes of an elastic sheet, stretched over a circular frame of radius $a$, are indexed by the mode number $j$ (which is a positive integer), and take the form

$\displaystyle \psi_j(\rho,t) = A_j \,J_0\left(\frac{z_j\,\rho}{a}\right)\cos(\omega_j\,t-\phi_j),$ (7.34)

where

$\displaystyle \omega_j = z_j\,\frac{v}{a},$ (7.35)

and $A_j$, $\phi_j$ are arbitrary constants. Figure 7.6 illustrates the spatial variation of the $j=5$ normal mode.


Table: 7.1 First few zeros of the Bessel function $J_0(z)$. Source: Abramowitz and Stegun 1965.
$j$ $z_j$ $j$ $z_j$
       
1 2.40482 6 18.07106
2 5.52007 7 21.21163
3 8.65372 8 24.35247
4 11.79153 9 27.49347
5 14.93091 10 30.63460


Figure 7.6: Density plot illustrating the spatial variation of the $j=5$ normal mode of a circular elastic sheet of radius $a$. Dark/light regions indicate large/small displacement amplitudes.
\includegraphics[width=0.8\textwidth]{Chapter07/fig7_06.eps}