Doble pèndol

Un exemple de doble pèndol.

En general, un doble pèndol és un sistema compost per dos pèndols, amb el segon penjant de l'extrem del primer. En el cas més simple, es tracta de dos pèndols simples, amb l'inferior penjant de la massa pendular del superior.

Normalment se sobreentén que ens referim a un doble pèndol pla, amb dos pèndols plans coplanaris. Aquest sistema físic posseeix dos graus de llibertat i exhibeix un ric comportament dinàmic. El seu moviment està governat per dues equacions diferencials ordinàries acoblades. Per sobre de certa energia, el seu moviment és caòtic.

Anàlisi del moviment del pèndol doble pla

Cinemàtica

A la cinemàtica només estem interessats a trobar les expressions de la posició, la velocitat, l'acceleració i en termes de les variables que especifiquen l'estat del doble pèndol, sense interessar-nos per les forces actuants. Ens servirem de les següents coordenades:

• x, i = posició horitzontal i vertical de la massa d'un pèndol
• θ = angle d'un pèndol respecte a la vertical (0 = vertical cap avall, antihorari és positiu)
• l = longitud de la vareta (constant)

Associarem al pèndol superior el subíndex 1, i al de baix el subíndex 2. Posarem l'origen de coordenades en el punt de pivot del pèndol superior. El sentit de les ordenades creixents es pren cap amunt.

A partir de consideracions trigonomètriques escrivim les expressions de les posicions x ₁ , i ₁ , x ₂ , i ₂ en termes dels angulos θ ₁ , θ ₂ :

$x_1=l_1\sin {\theta }_1$

$y_1=-l_1\cos {\theta }_1$

$x_2=x_1+l_2\sin {\theta }_2$

$y_2=y_1-l_2\cos {\theta }_2$

Derivant respecte al temps obtenim:

${\dot{x}}_1={\dot{\theta }}_1{l}_1\cos {\theta }_1$

${\dot{y}}_1={\dot{\theta }}_1{l}_1\sin {\theta }_1$

${\dot{x}}_2={\dot{x}}_1+{\dot{\theta }}_2{l}_2\cos {\theta }_2$

${\dot{y}}_2={\dot{y}}_1+{\dot{\theta }}_2{l}_2\sin {\theta }_2$

I derivant una segona vegada:

${\ddot{x}}_1=-{\dot{\theta }}_1^2{l}_1\sin {\theta }_1+{\ddot{\theta }}_1{l}_1\cos {\theta }_1$

${\ddot{y}}_1={\dot{\theta }}_1^2{l}_1\cos {\theta }_1+{\ddot{\theta }}_1{l}_1\sin {\theta }_1$

${\ddot{x}}_2={\ddot{x}}_1-{\dot{\theta }}_2^2{l}_2\sin {\theta }_2+{\ddot{\theta }}_2{l}_2\cos {\theta }_2$

${\ddot{y}}_2={\ddot{y}}_1+{\dot{\theta }}_2^2{l}_2\cos {\theta }_2+{\ddot{\theta }}_2{l}_2\sin {\theta }_2$

Forces

Definim les variables:

• T = tensió en la vareta
• M = massa del pèndol
• G = acceleració de la gravetat

Farem servir la llei de Newton $F=ma$, escrivint per separat les equacions de les components verticals i horitzontals de les forces.

Sobre la massa $m_1$ actuen la tensió a la part superior de la vareta $T_1$, la tensió en la part inferior de la vareta $T_2$, i la gravetat -m ₁ g :

$M_1{\ddot{x}}_1=-T_1\sin {\theta }_1+T_2\sin {\theta }_2$

$M_1{\ddot{i}}_1=T_1\cos {\theta }_1-T_2\cos {\theta }_2-m_{1}g$

Sobre la massa $m_2$, actuen la tensió $T_2$ i la gravetat -m ₂ g :

$M_2{\ddot{x}}_2=-T_2\sin {\theta }_2$

$M_2{\ddot{i}}_2=T_2\cos {\theta }_2-m_{2}g$

Equacions de moviment

A partir de les equacions anteriors, després de realitzar nombroses operacions algebraiques amb la finalitat de trobar les expressions de $\ddot{{\theta }_1}$, $\ddot{{\theta }_2}$ en termes de ${\theta }_1$, $\dot{{\theta }_1}$, ${\theta }_2$, $\dot{{\theta }_2}$, arribaríem a les equacions de moviment per al pèndol doble:

${\ddot{\theta }}_1=\frac{-g(2m_1+m_2)\sin {\theta }_1-m_{2}g\sin ({\theta }_1-2{\theta }_2)-2\sin ({\theta }_1-{\theta }_2)m_2({\dot{\theta }}_2^2{l}_2+{\dot{\theta }}_1^2{l}_1\cos ({\theta }_1-{\theta }_2))}{l_1(2m_1+m_2-m_2\cos (2{\theta }_1-2{\theta }_2))}$

${\ddot{\theta }}_2=\frac{2\sin ({\theta }_1-{\theta }_2)({\dot{\theta }}_1^2{l}_1(m_1+m_2)+g(m_1+m_2)\cos {\theta }_1+{\dot{\theta }}_2^2{l}_2{m}_2\cos ({\theta }_1-{\theta }_2))}{l_2(2m_1+m_2-m_2\cos (2{\theta }_1-2{\theta }_2))}$

Energia

L'energia cinètica ve expressada per:

$T=\frac{1}{2}{m}_1({\dot{x}}_1^2+{\dot{i}}_1^2)+\frac{1}{2}{m}_2({\dot{x}}_2^2+{\dot{i}}_2^2)=\frac{1}{2}{m}_1{l}_1^2{\dot{\theta }}_1^2+\frac{1}{2}{m}_2[L_1^2{\dot{\theta }}_1^2+L_2^2{\dot{\theta }}_2^2+2l_1{l}_2{\dot{\theta }}_1{\dot{\theta }}_2\cos ({\theta }_1-{\theta }_2)]$

L'energia potencial:

$V=m_{1}g{y}_1+m_{2}g{y}_2=-(m_1+m_2)g{l}_1\cos {\theta }_1-m_{2}g{l}_2\cos {\theta }_2$.

Per tant, el moviment es regirà per la lagrangiana

$ℒ=T-V=\frac{1}{2}(m_1+m_2)L_1^2{\dot{\theta }}_1^2+\frac{1}{2}{m}_2{l}_2^2{\dot{\theta }}_2^2+M_2{l}_1{l}_2{\dot{\theta }}_1{\dot{\theta }}_2\cos ({\theta }_1-{\theta }_2)+(m_1+m_2)g{l}_1\cos {\theta }_1+m_{2}g{l}_2\cos {\theta }_2$

Equacions de moviment de Lagrange

Usant les equacions de Lagrange en aquest cas particular són: Plantilla:Equació Calculant explícitament les derivades de l'expressió anterior s'arriba a: Plantilla:Equació Aquestes són les equacions de Lagrange per a un pèndol doble on hem escollit com coordenades generalitzades les polars i en el qual hi ha dues lligadures ($L_1$ i $L_2$ constants)

Vegeu també

• Pèndol
• Pèndol balístic
• Pèndol cicloïdal
• Pèndol cònic
• Pèndol de Foucault
• Pèndol de Newton
• Pèndol de Pohl
• Pèndol de torsió
• Pèndol esfèric
• Pèndol físic
• Pèndol simple
• Pèndol simple equivalent

Referències

Bibliografia

• Plantilla:Ref-llibre

Enllaços externs

• Eric W. Weisstein, Double pendulum, ScienceWorld, 2005 (en anglès, conté detalls sobre les equacions involucrades) i Rob Morris, Double Pendulum, Wolfram Demonstrations Project, 2007 (animacions).
• Animació que mostra moviment del pèndol doble i repartiment d'energia entre un i altre pèndol Plantilla:Webarchive
• Double pendulum physics simulation, MyPhysicsLab (en anglès, codi Javascript lliure d'ús disponible a Github)