Vibració d'una membrana circular

Una membrana elàstica sota tensió pot admetre vibracions transversals. Les propietats d'un pedaç de tambor ideal poden modelitzar-se amb les vibracions d'una membrana circular de gruix uniforme i subjecta a un marc rígid. A causa de les ressonàncies, a una certa freqüència de vibració (la freqüència de ressonància) la membrana pot emmagatzemar energia vibracional, amb la superfície oscil·lant en un patró característic d'ones estacionàries anomenat mode normal. Una membrana té infinits modes normals i el de menor freqüència s'anomena mode fonamental.

Un dels possibles modes de vibració d'un pedaç de tambor ideal (mode $u_{12}$ seguint la notació descrita més endavant). Al final de l'article es mostren altres modes.

Una membrana pot vibrar d'infinites maneres, cadascuna depenent de la forma de la membrana en un cert instant inicial i la velocitat transversal de cada punt de la membrana en aquest instant. Les vibracions de la membrana venen donades per les solucions de l'equació d'ona en dues dimensions amb condicions de vora de Dirichlet que representen la restricció al marc. Es pot demostrar que qualsevol vibració arbitrària de la membrana es pot descompondre en una sèrie, possiblement infinita, de modes normals de la membrana. Aquest resultat és anàleg a la descomposició d'una senyal temporal en sèries de Fourier

Motivació

L'anàlisi de la vibració d'un pedaç de tambor explica el comportament d'instruments de percussió com els tambors o les timbales. Tanmateix, també té aplicacions en biologia en el funcionament del timpà. Des d'un punt de vista education, els modes d'un objecte bidimensional són una forma visual d'explicar el significat dels modes, nodes, antinodes i fins i tot els nombres quàntics. Aquests conceptes són importants per entendre l'estructura de l'àtom.

Descripció del problema

Consideri's un disc obert $\Omega$ de radi $a$ centrat a l'origen, que representarà la membrana en repòs. A qualsevol temps $t,$ l'alçada de la membrana a un punt $(x,y)\in \Omega$ mesurat des de la posició de repòs es pot denotar com $u(x,y,t),$ que pot prendre valors tan positius com negatius. Sigui $\partial \Omega$ la frontera de $\Omega ,$ o sigui, la circumferència de radi $a$ centrada a l'origen, que representa el marc rígid que subjecta la membrana.

L'equació matemàtica que governa la vibració de la membrana és l'equació d'ona amb condició de vora de Dirichlet homogènia,

${\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{ per }}(x,y)\in \Omega \,$

$u=0{\text{ a }}\partial \Omega .\,$

Degut a la geometria circular de la membrana, serà convenient utilitzar coordenades cilíndriques, $(r,\theta ,z).$ En aquest cas, les equacions anteriors es poden escriure com

${\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{ per }}0\leq r<a,0\leq \theta \leq 2\pi \,$

$u=0{\text{ a }}r=a.\,$

Ara, $c$ és una constant positiva, que dona la velocitat a la qual una ona transversal es propaga per la membrana. En termes dels paràmetres físics del sistema, la velocitat d'ona ve donada per

$c={\sqrt {\frac {N_{rr}^{*}}{\rho h}}}$

on $N_{rr}^{*}$, és la força radial resultant a la vora de la membrana ($r=a$), $h$, és el gruix de la membrana, i $\rho$ és la densitat de la membrana. Si la membrana té una tensió uniforme, la força de tensió a un cert radi $r$ ve donada per

$F=rN_{rr}^{r}=rN_{\theta \theta }^{r}$

on $N_{\theta \theta }^{r}=N_{rr}^{r}$ és la resultant de a la membrana en la direcció azimutal.

Cas radialment simètric

Primer estudiarem els modes de vibració d'una membrana circular radialment simètrica. En aquest cas, la funció $u$ no depèn de l'angle $\theta ,$ i l'equació d'ona es redueix a

${\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).$

Buscarem solucions utilitzant separació de variables, $u(r,t)=R(r)T(t).$ Substituint aquesta expressió a l'equació anterior i dividint ambdós costats de la igualtat per $c^{2}R(r)T(t)$ s'obté

${\frac {T''(t)}{c^{2}T(t)}}={\frac {1}{R(r)}}\left(R''(r)+{\frac {1}{r}}R'(r)\right).$

El costat esquerre de l'equació no depèn de $r,$ i el costat dret no depèn de $t,$ de manera que ambdós costats han de ser iguals a una constant $K.$ D'aquesta manera, obtenim equations separates per $T(t)$ i $R(r)$:

$T''(t)=Kc^{2}T(t)\,$

$rR''(r)+R'(r)-KrR(r)=0.\,$

L'equació per $T(t)$ té solucions que creixen o decauen exponencialment per $K>0,$ són lineals o constants per $K=0$ i són periòdiques per $K<0$. Físicament, és d'esperar que una solució al problema de la vibració d'una membrana circular serà oscil·latòria en el temps, per tant l'únic cas possible és el tercer, $K<0,$ on $K=-\lambda ^{2}.$ Llavors, $T(t)$ és combinació lineal de sinus i cosinus

$T(t)=A\cos c\lambda t+B\sin c\lambda t.\,$

Estudiant ara l'equació per $R(r),$ i utilitzant que $K=-\lambda ^{2},$ totes les solucions d'aquesta equació diferencial de segon ordre són combinació lineal de funcions de Bessel d'ordre 0, ja que aquest és un cas particular de l'equació diferencial de Bessel:

$R(r)=c_{1}J_{0}(\lambda r)+c_{2}Y_{0}(\lambda r).\,$

La funció de Bessel $Y_{0}$ és no acotada quan $r\to 0,$ el que donaria un resultat físicament impossible per la vibració d'una membrana, de manera que la constant $c_{2}$ ha de ser nul·la. Podem suposar que $c_{1}=1,$ ja que en cas contrari aquesta constant pot ser absorbida per les constants $A$ i $B$ de $T(t).$ D'aquesta manera s'obté

$R(r)=J_{0}(\lambda r).$

La condició que l'alçada $u$ ha de ser zero a la vora de la membrana dona la condició

$R(a)=J_{0}(\lambda a)=0.$

La funció de Bessel $J_{0}$ té un nombre infinit d'arrels positives,

$0<\alpha _{01}<\alpha _{02}<\cdots$

així doncs $\lambda a=\alpha _{0n},$ per $n=1,2,\dots ,$ i d'aquesta manera

$R(r)=J_{0}\left({\frac {\alpha _{0n}}{a}}r\right).$

Per tant, les solucions radialment simètriques $u$ de la vibració d'una membrana circular es poden representar mitjançant separació de variables com

$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{ per }}n=1,2,\dots ,\,$

on $\lambda _{0n}=\alpha _{0n}/a.$

Cas general

El cas general, on $u$ també pot dependre de l'angle $\theta ,$ es resol de manera similar. Suposem que la solució admet separació de variables,

$u(r,\theta ,t)=R(r)\Theta (\theta )T(t).\,$

Substituint en l'equació d'ona i separant les variables s'obté

${\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$

on $K$ és una constant. Igual que abans, de l'equació de $T(t)$ es dedueix que $K=-\lambda ^{2}$ amb $\lambda >0$ i

$T(t)=A\cos c\lambda t+B\sin c\lambda t.\,$

De l'equació

${\frac {R''(r)}{R(r)}}+{\frac {R'(r)}{rR(r)}}+{\frac {\Theta ''(\theta )}{r^{2}\Theta (\theta )}}=-\lambda ^{2}$

s'obté, multiplicant ambdós costant per $r^{2}$ i separant les variables, que

$\lambda ^{2}r^{2}+{\frac {r^{2}R''(r)}{R(r)}}+{\frac {rR'(r)}{R(r)}}=L$

i

$-{\frac {\Theta ''(\theta )}{\Theta (\theta )}}=L,$

per alguna constant $L.$ Com que $\Theta (\theta )$ és periòdica, amb període $2\pi ,$ on $\theta$ és una variable angular, se segueix que

$\Theta (\theta )=C\cos m\theta +D\sin m\theta ,\,$

on $m=0,1,\dots$ i $C$ i $D$ són constants. Això també implica que $L=m^{2}.$

Tornant a l'equació per $R(r),$ la seva solució és combinació lineal de funcions de Bessel $J_{m}$ i $Y_{m}.$ Utilitzant un raonament similar al de l'apartat anterior obtenim

$R(r)=J_{m}(\lambda _{mn}r),\,$ $m=0,1,\dots ,$ $n=1,2,\dots ,$

on $\lambda _{mn}=\alpha _{mn}/a,$ amb $\alpha _{mn}$ la $n$-èsima arrel positiva de $J_{m}.$

Per tant, hem demostrat que totes les solucions en variables separades del problema de vibració d'una membrana circular tenen 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 )$

per $m=0,1,\dots ,n=1,2,\dots$

Animacions dels modes de vibració

A continuació es mostren alguns dels primers modes juntament amb els seus nombres quàntics. Les funcions d'ona anàlogues per l'àtom d'hidrogen també s'indiquen així com les freqüències angulars associades $\omega =\lambda _{mn}c={\frac {c}{a}}\alpha _{mn}$.

• 
  Mode $u_{01}$ (1s) amb $\alpha _{01}=2.40483$
• 
  Mode $u_{02}$ (2s) amb $\alpha _{02}=5.52008$
• 
  Mode $u_{03}$ (3s) amb $\alpha _{03}=8.65373$

• 
  Mode $u_{11}$ (2p) amb $\alpha _{11}=3.83171$
• 
  Mode $u_{12}$ (3p) amb $\alpha _{12}=7.01559$
• 
  Mode $u_{13}$ (4p) amb $\alpha _{13}=10.1735$

• 
  Mode $u_{21}$ (3d) amb $\alpha _{21}=5.13562$
• 
  Mode $u_{22}$ (4d) amb $\alpha _{22}=8.41724$
• 
  Mode $u_{23}$ (5d) amb $\alpha _{23}=11.6198$

Enllaços externs

• Applet Java que il·lustra els modes d'una membrana circular