Pendel

:(

Alla accelerationer (a, g + komposanter) är skalade med en faktor 5 för att det ska gå att se, varsågod.

Newton

Okej, den enda kraften som verkar på pendeln är gravitationen.
Vi kan enkelt splitta up den i komposanter (radiell och tangentell [checka gravity]) med lite trig.
Armen är stel och kommer att motverka all kraft i radiell riktning (gr) så vi blir bara kvar med den tangentella kraften (gt).
Vilket ger:

f=ma=mg·sin(θ)

a=g·sin(θ)

För att hitta den svepta sträckan kan vi använda formeln för båglängden i polära koordinater:

s=drdθ+r2dθ

Här är armens längd konstant vilket innebär drdθ=0 och ger:

s=r2dθ=rdθ=rθ

Då kan vi enkelt säga att:

a=s¨=rθ¨

Vilket leder till:

rθ¨=g·sin(θ)

θ¨=gr·sin(θ)

Kaboom! Där har vi vår diff-ekv.

Lagrange

Denna metod bygger på Euler-Lagrange ekvationen:

ddtLθ˙=Lθ

Med

L=TU

Där T är den kinesiska energin och U är den potentiella energin.

I pendeln blir L :

L=12ms˙2mgr(1cos(θ))

Där s˙=rθ˙ precis som i Newton-fallet.

Så då löser vi bara Euler-Lagrange:

Lθ=mgr·sin(θ)


ddtLθ˙=ddtmr2θ˙=mr2θ¨

Nu sätter vi dom lika med varanda och får:

θ¨=gr·sin(θ)

Samma som med Newton.
Det finns många andra sätt att härleda pendeln men att formatera MathML är ett helvete så det får räcka där.
Testa med Hamiltonisk mekanik om du känner dig uttråkad.
Där är H=T+U så man räknar snarare på systemets totala energi.