Ondas: Membrana

Ecuación de onda para membranas

En la sección Ondas: Cuerdas y Corrientes- Transversalesse mencionan algunos casos de ondas en una dimensión. Para el caso de dos dimensiones se considera una membrana estirada y flexible siendo esta la contraparte en dos dimensiones de la cuerda estirada y flexible.

Se considera que la membrana se encuentra en un principio en reposo en el plano xy con una tensión superficial ($\tau$) a lo largo de cualquier línea recta que cruce la membrana, el análogo bidimensional de la tensión del a cuerda, y que la membrana posee una densidad por unidad de área ($\sigma$).

Supongamos que tenemos un elemento rectangular $\Delta x$,$\Delta y$ de la membrana, cuando la membrana sufre un desplazamiento surge una fuerza neta en dirección del eje z por cada pareja de fuerzas de tensión $\tau \Delta x$ y $\tau \Delta y$ que actúan en las cuatro orillas del elemento rectangular. Cada par puede considerarse el equivalente a la tensión en una cuerda hipotética que consiste de una tira de la membrana donde su ancho mide $\Delta x$ y se extiendo a lo largo del eje y, o de ancho $\Delta y$ y se extiende a lo largo del eje x. Para la primera tira la fuerza neta en la dirección del eje z es:

$(1)\qquad \tau \Delta x\left({\frac {\partial \Psi }{\partial y}}_{y+\Delta y}-{\frac {\partial \Psi }{\partial y}}_{y}\right)=\tau \Delta x{\frac {\partial ^{2}\Psi }{\partial y^{2}}}\Delta y$

Para la segunda tira, la fuerza neta es:

$(2)\qquad \tau \Delta y\left({\frac {\partial \Psi }{\partial x}}_{x+\Delta x}-{\frac {\partial \Psi }{\partial x}}_{x}\right)=\tau \Delta y{\frac {\partial ^{2}\Psi }{\partial x^{2}}}\Delta x$

La suma de estas fuerzas debe de ser igual a la masa de la membrana $\sigma \Delta x\Delta y$ por su aceleración, eso es,

$(3)\qquad \Delta x\Delta y\tau \left({\frac {\partial ^{2}\Psi }{\partial x^{2}}}+{\frac {\partial ^{2}\Psi }{\partial y^{2}}}\right)=\Delta x\Delta y\sigma {\frac {\partial ^{2}\Psi }{\partial t^{2}}}$

Dividiendo la ecuación entre $\Delta x\Delta y$ se obtiene la ecuación de onda tal que,

$(4)\qquad {\frac {\partial ^{2}\Psi }{\partial x^{2}}}+{\frac {\partial ^{2}\Psi }{\partial y^{2}}}={\frac {1}{C_{m}}}{\frac {\partial ^{2}\Psi }{\partial t^{2}}}$

Donde,

$(5)\qquad C_{m}=\left({\frac {\tau }{\sigma }}\right)^{1/2}$

Es la velocidad de propagación para una onda transversal en la membrana.

--Alejandro Angel Galvan Garcia 21:50 26 mar 2012 (UTC)

---

Modos normales de una membrana circular

Vamos a estudiar ahora los modos normales de vibración de una membrana con límite fijo circular, por lo que resulta conveniente transformar a coordenadas polares la ecuación (4) para simplificar el tratamiento que a continuación se hace.

Usando las relaciones $x=r\cos(\theta )$, $y=r\mathrm {sen} (\theta )$, y por la regla de la cadena obtenemos

$(6)\qquad C_{m}^{2}({\frac {\partial ^{2}\Psi }{\partial r^{2}}}+{\frac {1}{r}}{\frac {\partial \Psi }{\partial r}}+{\frac {1}{r^{2}}}{\frac {\partial ^{2}\Psi }{\partial \theta ^{2}}})={\frac {\partial ^{2}\Psi }{\partial t^{2}}}$

Ahora, por el método de separación de variables trataremos de encontrar una solución de la forma producto a la ecuación anterior que satisfaga la condición de no desplazamiento en la frontera de la membrana. Como veremos en seguida, las variables de la ecuación (6) son separables en el sistema coordenado que estamos usando. Sea

$(7)\qquad \Psi (r,\theta ,t)=R(r)\Theta (\theta )T(t)$

Derivando y sustituyendo en la ecuación (6), obtenemos

$(8)\qquad R\Theta T''=C_{m}^{2}(R''\Theta T+{\frac {1}{r}}R'\Theta T+{\frac {1}{r^{2}}}R\Theta ''T)$

donde las primas denotan la diferenciación ordinaria. Dividiendo por $R\Theta T$ en ambos lados de la ecuación anterior nos queda

$(9)\qquad {\frac {C_{m}^{2}}{R}}(R''+{\frac {1}{r}}R')+{\frac {C_{m}^{2}}{r^{2}\Theta }}\Theta ''={\frac {1}{T}}T''$

Como el lado izquierdo no depende de la variable $t$ así como el derecho no depende de las variables $r$ y $\theta$, podemos concluir que para que sea cierta la igualdad ambas partes de la ecuación deben de ser iguales a una constante, que llamaremos de separación, por ejemplo $-\omega ^{2}$. Introduciendo esta constante en la ecuación anterior obtenemos las dos siguientes

$(10)\qquad T''+\omega ^{2}T=0$

$(11)\qquad {\frac {C_{m}^{2}}{R}}(R''+{\frac {1}{r}}R')+\omega ^{2}=-{\frac {C_{m}^{2}}{r^{2}\Theta }}\Theta ''$

La segunda ecuación podemos a su vez separarla en dos ecuaciones diferenciales ordinarias introduciendo una segunda variable de separación $n^{2}$. Después de multiplicar por $r^{2}/C_{m}^{2}$ a (11), obtenemos

$(12)\qquad \Theta ''+n^{2}\Theta =0$

$(13)\qquad R''+{\frac {1}{r}}R'+[{\frac {\omega ^{2}}{C_{m}^{2}}}-{\frac {n^{2}}{r^{2}}}]R=0$

Con el cambio de variable $u=\omega r/C_{m}^{2}$ en (13), obtenemos la ecuación de Bessel

$(14)\qquad R''+{\frac {1}{u}}R'+(1-{\frac {n^{2}}{u^{2}}})R=0$

En este caso las primas representan la diferenciación con respecto a la variable $u$. Siendo de segundo orden, tiene dos soluciones linealmente independientes para cada valor del parámetro m. Las dos soluciones de la ecuación (14) normlmente se presentan como $J_{n}(u)$ y $N_{n}(u)$. La solución $J_{n}(u)$ es conocida como función de Bessel de primer tipo de orden n. La otra solución $N_{n}(u)$, es conocida como función de Neumann ( o función de Bessel de segundo tipo). Como la función de Neumann no puede representar el desplazamiento de una membrana circular , solo se utiliza la función de Bessel de primer tipo, y esta función puede ser expresada en una serie infinita de la manera:

$(15)\qquad J_{n}(u)={\frac {u^{n}}{2^{n}n!}}\left[1-{\frac {u^{2}}{1!2^{2}(n+1)}}+{\frac {u^{4}}{2!2^{4}(n+1)(n+2)}}-...\right]$

El coeficiente numerico $1/2^{n}n!$ es puramente convencional. Todas las funciones de Besse a exepcion de $J_{0}$ se desvanecen en el origen y $J_{0}(0)=1$. La gráfica de la función de Bessel obedece un comportamiento similar al de una onda sinusoidal amortiguada.

Las funciones de Bessel obedecen las siguientes relaciones de recursión:

$(16)\qquad J_{n+1}(u)={\frac {2n}{u}}J_{n}(u)-J_{n-1}(u)$

$(17)\qquad {\frac {dJ_{n}(u)}{du}}=-{\frac {n}{u}}J_{n}(u)+J_{n-1}(u)$

Estas relaciones pueden establecerse directamente de la serie infinita dada en (15). Muestran que solo es necesario tener los valores numericos de $J_{0}$ y $J_{1}$ para conocer los valores de las funciones de Bessel de orden mayor, asi como las primeras derivadas.

Por otra parte, las soluciones generales de (10) y (11) son, respectivamente,

$(18)\qquad T(t)=b_{1}\cos(\omega t)+b_{2}\mathrm {sen} (\omega t)$

$(19)\qquad \Theta (\theta )=c_{1}\cos(n\theta )+c_{2}\mathrm {sen} (n\theta )$

Consideremos ahora algunos ejemplos de los modos normales de la membrana que estamos tratando. Cuando $m=0$, podemos escribir a (7) como

$(20)\qquad \Psi (r,t)=AJ_{0}({\frac {\omega }{C_{m}}}r)\cos(\omega t)$

Nótese que se uso también la condición inicial ${\frac {\partial \Psi }{\partial t}}=0$, es decir la membrana parte del reposo, de donde la constante $b_{2}=0$ en (18). Sea a el radio de las membranas, para satisfacer ahora la condición de frontera

$(21)\qquad \Psi (a,\theta ,t)=0$

es decir, que los límites estén fijos, tiene que cumplirse

$(22)\qquad J_{0}({\frac {\omega }{C_{m}}}r)=0\Rightarrow \omega _{0l}={\frac {C_{m}u_{0l}}{a}},l=1,2,...$

donde las $u_{0l}$ son las raíces de $J_{0}(u)=0$. Las $\omega _{0l}$ son pues las frecuencias en las que se excitan los modos normales de las membranas circulares.

Figura 1.

Figura 2.

Figura 3.

Figura 4.

--Alejandro Angel Galvan Garcia 19:29 29 mar 2012 (UTC)

--Roberto Verdel Aranda 05:03 28 mar 2012 (UTC)