# Vibraciones de una membrana circular

Las vibraciones de una membrana circular idealizada, esencialmente una membrana elástica de espesor uniforme fijada a un marco circular rígido, existen soluciones de la ecuación de onda con condiciones de contorno nulas.
Existe un número infinito de formas en las cuales la membrana puede vibrar, dependiendo de la forma de la deformación de la membrana en un cierto instante de tiempo inicial y de la derivada de la forma de la membrana en el instante inicial. Utilizando el método de separación de variables, es posible encontrar un conjunto de modos de vibración simples, y se puede demostrar que cualquier vibración compleja arbitraria de una membrana puede ser descompuesta en una serie de vibraciones simples (análogas a las serie de Fourier).

## Motivación

El análisis del problema de la membrana vibrante permite explicar el funcionamiento de instrumentos de percusión tales como los tambores y timbales. Sin embargo, existe también una aplicación biológica en explicar el funcionamiento del tímpano. Desde un punto de vista educativo los modos de un objeto bidimensional son una forma conveniente de visualmente demostrar el significado de modos, nodos, antinodo y aun los números cuánticos. Estos conceptos son importantes para la comprensión de la estructura del átomo.

## El problema

Considérese un disco de $\Omega$ de radio $a$ centrado en el origen, el cual representa la forma quieta de la membrana del tambor. En todo instante $t,$ la elevación de la forma de la membrana en un punto $(x,y)$ en $\Omega$ medido a partir de la forma quieta de la membrana se escribirá como $u(x,y,t),$ el cual puede tomar tanto valores positivos como negativos. Llamemos $\partial \Omega$ la frontera de $\Omega ,$ es decir, el círculo de radio $a$ centrado en el origen, que representa el marco rígido al cual se encuentra fijada la membrana.
La ecuación matemática que gobierna la vibración de la membrana es la ecuación de onda con condiciones de contorno nulas,
    ${\frac {\partial ^{2}u}{\partial t^{2}}}=c^{2}\left({\frac {\partial ^{2}u}{\partial x^{2}}}+{\frac {\partial ^{2}u}{\partial y^{2}}}\right){\text{ para }}(x,y)\in \Omega \,$
    $u=0{\text{ en }}\partial \Omega .\,$
A causa de la geometría circular de $\Omega$, será conveniente utilizar coordenadas cilíndricas, $(r,\theta ,t).$ Entonces, la ecuación previa se expresa como
    ${\frac {\partial ^{2}u}{\partial t^{2}}}=c^{2}\left({\frac {\partial ^{2}u}{\partial r^{2}}}+{\frac {1}{r}}{\frac {\partial u}{\partial r}}+{\frac {1}{r^{2}}}{\frac {\partial ^{2}u}{\partial \theta ^{2}}}\right){\text{ para }}0\leq r<a,0\leq \theta \leq 2\pi \,$
    $u=0{\text{ para }}r=a.\,$
Donde, $c$ es una constante positiva, que es la velocidad a la cual las ondas de vibración transversal se propagan en la membrana. En función de los parámetros físicos, la velocidad de onda c, queda expresada como
$c={\sqrt {\frac {N_{rr}^{*}}{\rho h}}}$ :
| Símbolo | Nombre |
| $N_{rr}^{*}$ | Resultante radial de la membrana en la frontera de la membrana ($r=a$) |
| $h$ | Espesor de la membrana |
| $\rho$ | Densidad de la memebrana |
Si la membrana posee una tensión uniforme, la fuerza de la tensión uniforme en un radio dado, $r$ queda expresada como
    $F=rN_{rr}^{r}=rN_{\theta \theta }^{r}$
donde $N_{\theta \theta }^{r}=N_{rr}^{r}$ es la resultante de la membrana en la dirección azimutal.

## El caso radial simétrico

Primero se analizan los modos de vibración posibles de una membrana que son simétricos de manera radial. O sea que la función $u$ no depende del ángulo $\theta ,$ y la ecuación de onda se simplifica quedando
    ${\frac {\partial ^{2}u}{\partial t^{2}}}=c^{2}\left({\frac {\partial ^{2}u}{\partial r^{2}}}+{\frac {1}{r}}{\frac {\partial u}{\partial r}}\right).$
Se buscan soluciones utilizando separación de variables, $u(r,t)=R(r)T(t).$ Substituyendo esto en la ecuación anterior y dividiendo ambos lados por $c^{2}R(r)T(t)$ se obtiene
    ${\frac {T''(t)}{c^{2}T(t)}}={\frac {1}{R(r)}}\left(R''(r)+{\frac {1}{r}}R'(r)\right).$
El lado izquierdo de la igualdad no depende de $r,$ y el lado derecho no depende de $t,$ por lo tanto se infiere que ambos lados deben ser iguales a una constante $K.$ Se obtienen de este modo ecuaciones separadas para $T(t)$ y $R(r)$:
    $T''(t)=Kc^{2}T(t)\,$
    $rR''(r)+R'(r)-KrR(r)=0.\,$
La ecuación de $T(t)$ tiene soluciones que crecen o decaen exponencialmente para $K>0,$ son lineales o constante para $K=0,$ y son periódicas si $K<0.$ Desde un punto de vista físico se espera que la solución al problema de la membrana que vibra sea oscilatorio en el tiempo, por lo que queda el tercer caso, $K<0,$ cuando $K=-\lambda ^{2}.$ Entonces, $T(t)$ es una combinación lineal de funciones seno y coseno,
    $T(t)=A\cos c\lambda t+B\sin c\lambda t.\,$
Analizando la ecuación para $R(r),$ notando que $K=-\lambda ^{2},$ todas las soluciones de esta ecuación diferencial de segundo orden son una combinación lineal de funciones de Bessel de orden 0,
    $R(r)=c_{1}J_{0}(\lambda r)+c_{2}Y_{0}(\lambda r).\,$
La función de Bessel $Y_{0}$ no se encuentra acotada para $r\to 0,$ lo que daría lugar a una solución del problema de la membrana que carece de sentido físico, por lo tanto la constante $c_{2}$ debe ser nula. También se supone que $c_{1}=1,$ ya que de todas formas esta constante será absorbida posteriormente en las constantes $A$ y $B$ provenientes de $T(t).$ De donde resulta que
    $R(r)=J_{0}(\lambda r).$
El requerimiento que la elevación $u$ sea cero en la frontera de la membrana determina la condición
    $R(a)=J_{0}(\lambda a)=0.$
La función de Bessel $J_{0}$ posee una cantidad infinita de raíces positivas,
    $0<\alpha _{01}<\alpha _{02}<\cdots$
Se obtiene que $\lambda a=\alpha _{0n},$ para $n=1,2,\dots ,$ por lo tanto
    $R(r)=J_{0}\left({\frac {\alpha _{0n}}{a}}r\right).$
Entonces, las soluciones simétricas radiales $u$ del problema de la membrana vibrante pueden ser expresadas en variables separadas como
    $u_{0n}(r,t)=\left(A\cos c\lambda _{0n}t+B\sin c\lambda _{0n}t\right)J_{0}\left(\lambda _{0n}r\right){\text{ for }}n=1,2,\dots ,\,$
donde $\lambda _{0n}=\alpha _{0n}/a.$

## El caso general

El caso general, cuando $u$ también puede depender del ángulo $\theta ,$ es tratada de manera similar. Se supone existe una solución mediante separación de variables,
    $u(r,\theta ,t)=R(r)\Theta (\theta )T(t).\,$
Substituyendo esto en la ecuación de onda y separando las variables, se obtiene
    ${\frac {T''(t)}{c^{2}T(t)}}={\frac {R''(r)}{R(r)}}+{\frac {R'(r)}{rR(r)}}+{\frac {\Theta ''(\theta )}{r^{2}\Theta (\theta )}}=K$
donde $K$ es una constante. Donde al igual que antes, de la ecuación para $T(t)$ se concluye que $K=-\lambda ^{2}$ con $\lambda >0$ y
    $T(t)=A\cos c\lambda t+B\sin c\lambda t.\,$
De la ecuación
    ${\frac {R''(r)}{R(r)}}+{\frac {R'(r)}{rR(r)}}+{\frac {\Theta ''(\theta )}{r^{2}\Theta (\theta )}}=-\lambda ^{2}$
resulta, multiplicando ambos lados por $r^{2}$ y separando variables, que
    $\lambda ^{2}r^{2}+{\frac {r^{2}R''(r)}{R(r)}}+{\frac {rR'(r)}{R(r)}}=L$
y
    $-{\frac {\Theta ''(\theta )}{\Theta (\theta )}}=L,$
para alguna constante $L.$ Dado que $\Theta (\theta )$ es periódica, con período $2\pi ,$ $\theta$ es una variable angular, se deduce que
    $\Theta (\theta )=C\cos m\theta +D\sin m\theta ,\,$
donde $m=0,1,\dots$ y $C$ y $D$ son ciertas constantes. Esto implica que $-L=m^{2}.$
Analizando ahora la ecuación para $R(r),$ su solución es una combinación lineal de funciones de Bessel $J_{m}$ y $Y_{m}.$ Con un argumento similar al utilizado con anterioridad, se obtiene que
    $R(r)=J_{m}(\lambda _{mn}r),\,$ $m=0,1,\dots ,$ $n=1,2,\dots ,$
donde $\lambda _{mn}=\alpha _{mn}/a,$ con $\alpha _{mn}$ la $n$-ésima raíz positiva de $J_{m}.$
Se ha demostrado que todas las soluciones en variables separadas del problema de la membrana circular vibrante son de la forma
    $u_{mn}(r,\theta ,t)=\left(A\cos c\lambda _{mn}t+B\sin c\lambda _{mn}t\right)J_{m}\left(\lambda _{mn}r\right)(C\cos m\theta +D\sin m\theta )$
para $m=0,1,\dots ,n=1,2,\dots$

## Animaciones de varios modos de vibración

A continuación se muestran algunos modos de vibración junto con sus números cuánticos. Las funciones de onda análogas del átomo de hidrógeno también se indican como la frecuencia angular asociada $\omega =\lambda _{mn}$.
* Modo $u_{01}$ (1s) con $\lambda _{01}=2.40483$
* Modo $u_{02}$ (2s) con $\lambda _{02}=5.52008$
* Modo $u_{03}$ (3s) con $\lambda _{03}=8.65373$
* Modo $u_{11}$ (2p) con $\lambda _{11}=3.83171$
* Modo $u_{12}$ (3p) con $\lambda _{12}=7.01559$
* Modo $u_{13}$ (4p) con $\lambda _{13}=10.1735$
* Modo $u_{21}$ (3d) con $\lambda _{21}=5.13562$
* Modo $u_{22}$ (4d) con $\lambda _{22}=8.41724$
* Modo $u_{23}$ (5d) con $\lambda _{23}=11.6198$