Settes føringene inn i L, forsvinner føringskreftene fra regningen. Det er en fordel helt til du trenger dem. Snorkraften i en pendel og normalkraften på en kloss som glir på en kuleflate, er føringskrefter, og for klossen er det normalkraften som avgjør når den forlater flaten. Da må føringen holdes utenfor L og tas med som en egen ligning.
Behold alle koordinatene q1,…,qn og skriv føringen som f(q,t)=0. Legg den til Lagrangefunksjonen ganget med en ny, ukjent funksjon λ(t):
L′=L+λf(q,t)
λ kalles en Lagrange-multiplikator og behandles som en koordinat til. L′ inneholder ikke λ˙, så Lagranges ligning for λ er ∂L′/∂λ=f=0. Føringen kommer tilbake som en av bevegelsesligningene. Ligningene for qk får et ledd fra λf på høyresida:
dtd∂q˙k∂L−∂qk∂L=α∑λα∂qk∂fαLagranges ligninger med multiplikatorer, for m føringer fα(q,t)=0. Sammen med føringene er det n+m ligninger for de n+m ukjente funksjonene qk(t) og λα(t). Høyresida er føringskraften i retning qk.
Venstresida er bevegelsesligningen uten føringer, så høyresida er den generaliserte kraften som holder føringen oppfylt. Den peker langs gradienten til f, altså vinkelrett på flaten f=0 som bevegelsen er bundet til. Det er grunnen til at føringskraften ikke gjør arbeid langs de tillatte forskyvningene, og til at den kunne sløyfes da føringene ble satt inn i L.
Ta pendelen, ei kule i ei snor med lengde l, med avstanden r fra opphenget og vinkelen θ fra loddlinja som koordinater. Føringen er f=r−l=0, og
L′=21m(r˙2+r2θ˙2)+mgrcosθ+λ(r−l)
Ligningen for λ gir r=l. Ligningen for r har ∂f/∂r=1 på høyresida:
mr¨−mrθ˙2−mgcosθ=λ
Sett inn r=l og r¨=0, så er λ=−(mlθ˙2+mgcosθ). f har dimensjon lengde, så λ har dimensjon kraft: den er kraften fra snora i retning r, og fortegnet sier at den peker inn mot opphenget. Snorkraften er S=−λ=mlθ˙2+mgcosθ, sentripetalkraften pluss komponenten av tyngden langs snora, det samme som Newton gir. Ligningen for θ har ingen λ, for f er uavhengig av θ, og med r=l er den ml2θ¨=−mglsinθ som før.
Energien gir S som funksjon av vinkelen alene. Har kula farten v0 i bunnpunktet, er 21ml2θ˙2−mglcosθ=21mv02−mgl, og
S=lmv02−2mg+3mgcosθ
S er størst i bunnpunktet og minker oppover. Ei snor kan bare trekke, så der S blir null, blir den slakk, og kula går i en kastebane. Det skjer ved cosθ=(2−v02/gl)/3, over vannrett. Med v02≤2gl snur kula før den når vannrett, og S rekker aldri null. Med v02≥5gl er S≥0 også i toppunktet, og kula går rundt. Ei stiv stang kan dytte like godt som trekke, og med stang skifter λ fortegn uten at noe skjer.
Snorkraften i pendelen
Kula henger i ei snor og får farten glidebryteren setter i bunnpunktet. Den oransje pila er snorkraften S, den grå er tyngden mg, tegnet i samme målestokk. Der S når null, blir snora slakk, og kula går i en kastebane til snora strammes igjen. De oransje prikkene er stedene der det skjer. Over 5gl i bunnpunktet holder snora seg stram hele veien rundt.
Fortegnet på λ er en konvensjon. Skriver du L−λf i stedet, bytter λ fortegn og ingenting annet. Føringskraften er den samme.
Hvorfor de frie koordinatene aldri ser multiplikatoren
Bytt koordinater fra q1,…,qn til et sett der de m føringsfunksjonene fα selv er koordinater, sammen med n−m frie koordinater q~i. Det går så lenge føringene er uavhengige av hverandre. Lagranges ligninger har samme form i de nye koordinatene, og høyresida i ligningen for q~i er ∑αλα∂fα/∂q~i. Men fα og q~i er uavhengige koordinater, så ∂fα/∂q~i=0. Ligningene for de frie koordinatene er Lagranges ligninger for L med føringene satt inn, uten spor av λ. Multiplikatorene lever bare i ligningene for koordinatene føringene låser, og det er de ligningene som gir føringskreftene.
Regneeksempel Kloss som glir av en kuleflate
En liten kloss med masse m glir friksjonsfritt fra ro på toppen av en fast kule med radius R. Ved hvilken vinkel forlater den overflaten?
Vis løsningSkjul løsning
Bruk avstanden r fra sentrum og vinkelen θ fra toppen som koordinater, med føringen f=r−R=0 og V=mgrcosθ:
L′=21m(r˙2+r2θ˙2)−mgrcosθ+λ(r−R)
Ligningen for r er
mr¨−mrθ˙2+mgcosθ=λ
Med r=R og r¨=0 blir
λ=mgcosθ−mRθ˙2
Det er normalkraften, positiv når flaten dytter utover. Energien fjerner θ˙: 21mR2θ˙2+mgRcosθ=mgR gir mRθ˙2=2mg(1−cosθ), og
λ=mg(3cosθ−2)
Normalkraften er mg på toppen og synker med vinkelen. Flaten kan ikke trekke, så klossen forlater den der λ=0:
cosθ=32⇒θ=48,2°
Det er i høyden 32R over sentrum, med farten v2=R2θ˙2=32gR. Derfra går klossen i en kastebane.
Svar
cosθ=2/3, altså θ=48,2° fra toppen, uavhengig av m og R.
Consider a uniform ladder of length l and total mass m, as shown in
Fig. 1.
Figure 1: Ladder leaning against a wall.
The ladder starts to slide without friction. The initial conditions are
θ(0)=θ0 and θ˙(0)=0.
a) Show that the Lagrangian of the systems [sic] is
L=21m(x˙cm2+y˙cm2)+21Icmθ˙2−mgycm,(1)
where Icm=121ml2.
b) The system is described in terms of the three coordinates
xcm, ycm, and θ. As long as the ladder is in
contact with the wall, the system is constrained, with a single degree of
freedom. Write down the two constraint equations and classify the
constraints.
c) Use the Lagrange multiplier method and write down the Lagrangian for the
constrained system. Derive the equations of motion.
d) Find the constraint forces as functions of m, g, and θ. For
which value of θ does the ladder lose contact with the wall?
Vis løsningSkjul løsning
a) Den kinetiske energien til et stivt legeme er massesenterets
translasjonsenergi pluss rotasjonsenergien om massesenteret,
21mvcm2+21Icmθ˙2, og for en
jevn stang er Icm=ml2/12. Den potensielle energien er
mgycm.
b) Med veggen langs y-aksen og gulvet langs x-aksen ligger massesenteret
midt på stigen: f1=xcm−2lcosθ=0 og
f2=ycm−2lsinθ=0. Begge er holonome, og
tiden inngår ikke.
c) L′=L+λ1f1+λ2f2. Lagranges ligninger for de tre
koordinatene:
λ1 er normalkraften fra veggen og λ2 normalkraften fra
gulvet.
d) Deriver føringene to ganger og sett inn:
λ1=−21ml(cosθθ˙2+sinθθ¨) og
λ2=mg+21ml(cosθθ¨−sinθθ˙2).
I θ-ligningen faller θ˙2 bort, og igjen står
31ml2θ¨=−21mglcosθ. Gang med θ˙ og
integrer fra start: θ˙2=l3g(sinθ0−sinθ).
Sett θ¨ og θ˙2 inn i uttrykkene over:
Stigen mister kontakten med veggen når λ1=0, altså når
sinθ=32sinθ0.
Svar
Kontakten med veggen forsvinner ved sinθ=32sinθ0.
Generalisert impuls og sykliske koordinater
Lagranges ligning for qi kan leses som en bevegelsesligning for én størrelse:
pi=∂q˙i∂L,p˙i=∂qi∂LDen generaliserte impulsen til koordinaten qi, også kalt kanonisk eller konjugert impuls. Lagranges ligning sier at den endres med den generaliserte kraften ∂L/∂qi.
For en kartesisk koordinat er px=mx˙, den vanlige impulsen. For pendelvinkelen er pθ=ml2θ˙, dreieimpulsen om opphenget. For den ladde partikkelen er p=mr˙+qA, og den er ikke den mekaniske impulsen. Den generaliserte impulsen er det L sier den er.
Mangler en koordinat i L, er ∂L/∂qi=0, og da sier ligningen at pi er konstant. En slik koordinat kalles syklisk. Det er den raskeste veien til en bevaringslov. ϕ mangler i Lagrangefunksjonen til den sfæriske pendelen, så pϕ=ml2ϕ˙sin2θ, dreieimpulsen om loddlinja, er bevart. Om en symmetri viser seg som en manglende koordinat, avhenger av koordinatvalget. Et sentralpotensial gjør ϕ syklisk i polarkoordinater og ingen koordinat syklisk i kartesiske. Velg koordinatene etter symmetrien.
Energifunksjonen
Tida kan også mangle i L. Størrelsen som hører til det, er
HdtdH=i∑q˙i∂q˙i∂L−L=−∂t∂LEnergifunksjonen. Den er bevart når L ikke avhenger eksplisitt av tida.
Deriver H langs en bane og bruk pi=∂L/∂q˙i og p˙i=∂L/∂qi:
Leddene med q¨i opphever hverandre fordi pi=∂L/∂q˙i, og leddene med q˙i fordi p˙i=∂L/∂qi.
H er energien når posisjonene ra(q) ikke avhenger av tida og V ikke avhenger av hastighetene. Da er T=21∑ijMij(q)q˙iq˙j kvadratisk i hastighetene, så ∑iq˙i∂T/∂q˙i=2T, og H=2T−(T−V)=T+V. Avhenger føringene av tida, får T ledd av lavere grad i q˙, og H og E skiller lag. For perla på en ring som roterer med fast ω er H=E−ma2ω2sin2θ. H er bevart, E er det ikke, for motoren som holder ω fast, gjør arbeid.
Ei perle glir friksjonsfritt på en rett stang som roterer i horisontalplanet med fast vinkelfart ω, med avstanden r fra aksen som koordinat. Hva er bevart?
Vis svarSkjul svar
Posisjonen er (rcosωt,rsinωt), så L=T=21m(r˙2+r2ω2). Tida forsvant da fartene ble kvadrert, så ∂L/∂t=0, og H=r˙mr˙−L=21mr˙2−21mr2ω2 er bevart. Energien E=T er ikke bevart. Ligningen for r er r¨=ω2r, perla slynges utover, og motoren som holder ω fast, gjør arbeidet. Leddet −21mr2ω2 i H er sentrifugalpotensialet.
Regneeksempel Den sfæriske pendelen mellom to breddesirkler
Ei kule henger i en stang med lengde l=1,0 m og slippes fra θ=60° fra loddlinja med farten v=1,4 m/s vannrett, vinkelrett på planet gjennom stanga og loddlinja. Hvor nær bunnpunktet kommer den?
Vis løsningSkjul løsning
Med θ fra loddlinja og ϕ rundt den er L=21ml2(θ˙2+ϕ˙2sin2θ)+mglcosθ. ϕ er syklisk, og ∂L/∂t=0, så to størrelser er bevart:
Kula beveger seg der E≥Veff, altså mellom to vendepunkter. I vendepunktene er θ˙=0, og startpunktet er det ene, for farten er vannrett der: pϕ=mlvsin60° og E=21mv2−mglcos60°. Det andre er den andre løsningen av Veff(θ)=E. Med u=cosθ, u1=cos60°=0,5 og k=v2/2gl=1,96/19,62=0,0999 blir ligningen, etter multiplikasjon med sin2θ/mgl,
k(1−u12)−u(1−u2)=(k−u1)(1−u2)
Flytt alt over på én side og del ut den kjente roten u=u1:
(u−u1)(u2+ku+ku1−1)=0
Andregradsligningen u2+0,0999u−0,950=0 har den fysiske roten
u=2−0,0999+0,0100+3,800=0,926
altså θ2=22,2°. Kula svinger mellom 22,2° og 60° og når aldri bunnpunktet, for Veff→∞ når θ→0 så lenge pϕ=0. En sirkelbane ved 60° krever v2=glsin2θ/cosθ, altså 3,8 m/s. Med 1,4 m/s faller kula innover.
Svar
Til θ=22,2° fra loddlinja, det andre vendepunktet av Veff(θ)=E.
Fiktive krefter fra Lagranges ligninger
De fiktive kreftene kom fra å derivere vektorer i et roterende system. De kommer også rett ut av Lagrangefunksjonen, uten noe kraftregnskap. La K′ rotere med konstant ω om en akse gjennom origo, og bruk komponentene av r i K′ som koordinater. Hastigheten i inertialsystemet er v=v′+ω×r, så
L=21m(v′+ω×r)2−V(r)Lagrangefunksjonen i det roterende systemet K′. v′ er hastigheten målt i K′, og ω er konstant.
Den generaliserte impulsen er p=∂L/∂v′=m(v′+ω×r)=mv, impulsen målt i K, ikke mv′. Den deriverte med hensyn på posisjonen er
∂r∂L=mv′×ω−mω×(ω×r)−∇V
og dp/dt=m(a′+ω×v′). Lagranges ligninger gir
ma′=−∇V−2mω×v′−mω×(ω×r)
Corioliskraften og sentrifugalkraften, som før. Et ledd −mω˙×r kommer til når ω ikke er konstant.
Lagrangefunksjonen inneholder ikke tida, så energifunksjonen er bevart:
H=v′⋅p−L=21mv′2+V−21m∣ω×r∣2
Det siste leddet er sentrifugalpotensialet. H er energien sett fra det roterende systemet, med sentrifugalkraften regnet som en konservativ kraft, og den er ikke den samme størrelsen som 21mv2+V i K.
Noethers teorem
En syklisk koordinat er det enkleste tilfellet av noe mer generelt. At ϕ mangler i L, betyr at forskyvningen ϕ→ϕ+ϵ lar L stå uendret, og at systemet ser likt ut etter en rotasjon om loddlinja. Det er en kontinuerlig symmetri: en familie av transformasjoner med en parameter ϵ som kan gjøres så liten du vil. Noethers teorem sier at hver slik symmetri gir en bevart størrelse, også når symmetrien ikke er en manglende koordinat i de koordinatene du bruker.
δL=ϵdtdF⇒Q=i∑piFi−F=konstantNoethers teorem for forskyvningen qi→qi+ϵFi. δL er endringen i L til første orden i ϵ.
Bevis
Tida endres ikke, så δq˙i=dtdδqi=ϵF˙i. Endringen i L til første orden er
der Lagranges ligninger ∂L/∂qi=p˙i er brukt i det andre trinnet. Er transformasjonen en symmetri, er den samme endringen også ϵdF/dt. Trekk fra:
dtd(i∑piFi−F)=0
Kravet δL=ϵdF/dt må gjelde for enhver bane q(t), ikke bare for den som løser bevegelsesligningene. Langs den virkelige banen er δL alltid en total tidsderivert, som beviset viser, så der sier kravet ingenting. Symmetrien er en egenskap ved L, bevaringen en egenskap ved løsningene.
En syklisk koordinat qj er tilfellet Fi=δij og F=0, med Q=pj. Bevaringslovene fra Newton kommer på samme måte. Se på et lukket system av partikler med L=∑a21mar˙a2−V, der V bare avhenger av avstandene ∣ra−rb∣. Forskyv alle partiklene like mye, ra→ra+ϵn^. Avstandene og hastighetene er uendret, så δL=0 og F=0, med Fa=n^ for hver partikkel. Den bevarte størrelsen er Q=∑apa⋅n^, den totale impulsen langs n^, og siden n^ er vilkårlig, hele P. Roter i stedet alle partiklene en liten vinkel ϵ om aksen n^, ra→ra+ϵn^×ra. Lengder og avstander er uendret til første orden i ϵ, og hastighetene roterer likt, så δL=0 igjen, og
Q=a∑pa⋅(n^×ra)=n^⋅a∑ra×pa
er den totale dreieimpulsen langs n^. At rommet er likt overalt, gir impulsbevaring, og at det er likt i alle retninger, gir dreieimpulsbevaring. Tidsforskyvning t→t+ϵ passer ikke rett inn i formen qi→qi+ϵFi, men resultatet er det samme: er L uavhengig av tida, er H bevart.
I disse eksemplene er δL=0. Teoremet er laget for å tåle mer, for en transformasjon kan være en symmetri selv om L endres, så lenge endringen er en total tidsderivert.
En fri partikkel har L=21mx˙2. Vis at x→x+ϵt er en symmetri, og finn den bevarte størrelsen.
Vis svarSkjul svar
Her er Fx=t, så δx˙=ϵ og δL=mx˙ϵ=ϵdtd(mx). Det er en total tidsderivert med F=mx, så transformasjonen er en symmetri selv om L endres. Noethers teorem gir Q=pxFx−F=mx˙t−mx. Konstant Q betyr x=x˙t−Q/m: partikkelen går med konstant fart, og −Q/m er der den var ved t=0. Transformasjonen er overgangen til et inertialsystem som beveger seg med farten −ϵ, og at den er en symmetri, er Galileis relativitetsprinsipp.
Speiling r→−r er også en symmetri av mange systemer, men den har ingen liten parameter ϵ, og gir ingen bevart størrelse i klassisk mekanikk.
b) Consider a system described by a Lagrangian that is a function of a set
of generalized coordinates qi, their time derivatives q˙i, and
time t, L=L(qi,q˙i,t). A set of coordinate transformations are
given by
qi→qi′=qi+δqi=qi+ϵFi(qj,q˙j,t).(6)
where ϵ is a small dimensionless parameter. Show that the
corresponding change in the Lagrangian can be written as
δL=ϵdtdF(7)
for a suitable function F=F(qj,q˙j,t).
c) Consider a system consisting of two particles of mass m1 and m2 in
three dimensions. The potential energy V is a function of the distance
∣r12∣=∣r1−r2∣, i.e. V=V(∣r12∣). Prove
that the total angular momentum of the system is conserved. Hint: Use
Eq. (7). Are there other conserved quantities? (proof not required).
Vis løsningSkjul løsning
a) Til hver kontinuerlig symmetri i Lagrangefunksjonen hører en bevart
størrelse: endrer L seg bare med en total tidsderivert,
δL=ϵdF/dt, under transformasjonen
qi→qi+ϵFi, så er Q=∑ipiFi−F konstant langs
bevegelsen.
b) Til første orden i ϵ er
δL=i∑∂qi∂LϵFi+i∑∂q˙i∂LϵF˙i.
Langs en løsning gjelder Lagranges ligninger,
∂L/∂qi=p˙i med pi=∂L/∂q˙i,
så summen er ϵ∑i(p˙iFi+piF˙i)=ϵdtd∑ipiFi.
Altså F=∑ipiFi.
Med denne F er Q=∑ipiFi−F identisk null langs løsningen, så
regningen i b) gjelder for enhver transformasjon og gir ingen bevart størrelse
alene. En symmetri krever at δL=ϵdF/dt holder for alle
baner, uten å bruke Lagranges ligninger, og det er dette c) bruker.
c) Roter begge partiklene en liten vinkel ϵ om en akse n^:
δra=ϵn^×ra. Avstanden ∣r12∣ og
fartene endrer seg ikke, så δL=0. Etter b) er da
dtd∑apa⋅(n^×ra)=0, og
pa⋅(n^×ra)=n^⋅(ra×pa), så
n^⋅L er bevart for enhver retning n^. Hele
dreieimpulsen L=r1×p1+r2×p2 er
bevart. L endrer seg heller ikke under en felles forskyvning av begge
partiklene eller under en forskyvning i tid, så total impuls og energi er
også bevart.
Svar
c) L er bevart fordi V(∣r12∣) er rotasjonsinvariant; total impuls og energi er også bevart.