Lagrange Mekaniği: Koordinat Seçiminden Hareket Denklemlerine

Serbestlik derecesini doğru seçmek: $q$ koordinatları

Bir sistemin serbestlik derecesi, konumunu belirlemek için gereken bağımsız değişken sayısıdır. Genelleştirilmiş koordinat $q_i$, bu değişkenlerden herhangi biri olabilir; uzunluk, açı veya başka bir uygun parametre seçilebilir. Koordinatın mutlaka metre cinsinden olması gerekmez. Örneğin tek bir düzlemsel sarkaç için ipin uzunluğu sabitse, konumu belirlemek üzere açı $\theta$ yeterlidir ve sistemin bir serbestlik derecesi vardır.

Koordinat seçiminin temel ölçütü şudur: Seçilen $q_i$ değişkenleri, sistemin izin verilen konfigürasyonlarını bağımsız biçimde tanımlamalıdır. Fazladan koordinat seçilirse aralarında kısıt denklemleri oluşur; eksik koordinat seçilirse sistemin hareketi tam tanımlanamaz. Bu nedenle önce fiziksel kısıtlar belirlenir, sonra bağımsız koordinatlar seçilir.

Kartezyen koordinatlar kullanıldığında bir parçacığın konumu $x$ ve $y$ ile yazılabilir. Ancak parçacık bir çember üzerinde hareket ediyorsa $x^2+y^2=\ell^2$ kısıtı eklenir. Bunun yerine $x=\ell\sin\theta$ ve $y=-\ell\cos\theta$ gibi ifadeler kullanılırsa kısıt doğrudan sağlanır. Bu durumda ip gerilmesi denklemin içinde ayrıca bilinmeyen olarak yer almaz. Bu kolaylık, yalnızca seçilen koordinatların kısıtı gerçekten içermesi ve kısıtın uygun nitelikte olması durumunda geçerlidir; her kısıt kuvveti otomatik olarak ortadan kalkmaz.

Lagrangian nasıl kurulur: $T$ ve $V$ ifadelerini aynı koordinatta yazmak

Standart mekanik problemlerde Lagrangian çoğunlukla

$L(q_i,\dot q_i,t)=T(q_i,\dot q_i,t)-V(q_i,t)$

şeklinde kurulur. Burada $\dot q_i=dq_i/dt$ koordinatın zamana göre türevidir. Kinetik enerji hızlara, potansiyel enerji ise çoğunlukla konuma bağlıdır. Önemli nokta, $T$ ve $V$ ifadelerinin aynı koordinat sistemi ve aynı referans seçimiyle yazılmasıdır.

Potansiyel enerji kullanılabilmesi için kuvvetin ilgili bölgede bir potansiyelden türetilebilmesi gerekir. Yerçekimi gibi korunumlu kuvvetlerde potansiyel enerjiye geçmek doğrudan mümkündür. Sürtünme gibi enerji dağıtan kuvvetler ise standart $L=T-V$ içine yalnızca potansiyel enerji gibi eklenemez; bunlar için genelleştirilmiş kuvvetler veya uygun bir sönüm modeli gerekir.

Potansiyel enerjinin sıfır seviyesi değiştirilebilir. $V$'ye sabit bir sayı eklemek $L$'yi sabit kadar değiştirir, fakat $\partial L/\partial q_i$ ve $\partial L/\partial\dot q_i$ değişmediğinden hareket denklemleri değişmez. Buna karşılık yanlış işaret kullanmak sonucu değiştirir. Örneğin düşey eksen yukarı yönlü ve sarkaç açısı en alt konumdan ölçülüyorsa yerçekimi potansiyelinin açıya bağlı ifadesi dikkatle kurulmalıdır.

Durağan etkiden Euler-Lagrange denklemine geçiş

Etki, Lagrangian'ın zaman aralığında integralidir:

$S=\int_{t_1}^{t_2}L\,dt$.

Durağan etki ilkesi, başlangıç ve bitiş durumları sabit tutulduğunda fiziksel yolun, küçük yol değişimleri karşısında $S$'yi durağan yapan yol olduğunu ifade eder. Durağan değer minimum, maksimum veya daha genel bir eyer davranışı olabilir; bu nedenle yalnızca 'en küçük etki' demek her durumda doğru değildir.

Bu ilkeden her bağımsız koordinat için şu denklem elde edilir:

$\frac{d}{dt}\left(\frac{\partial L}{\partial\dot q_i}\right)-\frac{\partial L}{\partial q_i}=0$.

Uygulama sırası nettir: Önce $L$ yazılır, sonra $\partial L/\partial\dot q_i$ alınır, bunun zamana göre türevi hesaplanır ve son olarak $\partial L/\partial q_i$ çıkarılır. Sonuç, $q_i(t)$ için diferansiyel hareket denklemidir.

$p_i=\partial L/\partial\dot q_i$ ifadesi genelleştirilmiş momentum olarak adlandırılır. Böylece denklem $\dot p_i=\partial L/\partial q_i$ biçiminde de yazılabilir. Ancak $p_i$ her durumda basitçe 'kütle çarpı hız' değildir; bunun biçimi seçilen koordinata ve Lagrangian'a bağlıdır. Kartezyen koordinatta bu tanım çoğu basit sistemde alışılmış momentuma karşılık gelirken açısal koordinatta açısal momentum benzeri bir büyüklük verir.

Sabit açılı sarkaçta hareket denklemi: tek koordinatlı uygulama

Uzunluğu $\ell$, kütlesi $m$ olan ve sürtünmesiz kabul edilen basit sarkaçta tek bağımsız koordinat açı $\theta$'dır. Kütlenin hızı, açısal hızla $v=\ell\dot\theta$ ilişkisine sahip olduğundan kinetik enerji

$T=\frac{1}{2}m\ell^2\dot\theta^2$

olur. En alt konuma göre yerçekimi potansiyeli

$V=mg\ell(1-\cos\theta)$

şeklinde seçilebilir. Dolayısıyla

$L=\frac{1}{2}m\ell^2\dot\theta^2-mg\ell(1-\cos\theta)$.

Şimdi $\theta$ için Euler-Lagrange denklemi uygulanır. Türevler

$\frac{\partial L}{\partial\dot\theta}=m\ell^2\dot\theta$

ve

$\frac{\partial L}{\partial\theta}=-mg\ell\sin\theta$

olduğundan

$m\ell^2\ddot\theta+mg\ell\sin\theta=0$

elde edilir. $m\ell$ ortak çarpanı sıfır olmayan bir sistem için sadeleştirilirse

$\ddot\theta+\frac{g}{\ell}\sin\theta=0$

sonucuna ulaşılır.

Bu denklem küçük açı yaklaşımı yapılmadıkça doğrusal değildir. Yalnızca $|\theta|$ küçük olduğunda $\sin\theta\approx\theta$ alınarak $\ddot\theta+(g/\ell)\theta=0$ biçimindeki basit harmonik salınım denklemi elde edilir. Dolayısıyla küçük salınım sonucu genel sarkaç hareketine değil, bu özel yaklaşım koşuluna aittir.

Simetri, korunum yasası ve yöntemin sınırları

Lagrangian'ın belirli bir değişkene açıkça bağlı olmaması önemli bir karar kuralı verir. Eğer $L$ bir koordinata değil, o koordinatın hızına bağlıysa, yani $\partial L/\partial q_i=0$ ise Euler-Lagrange denklemi $d(\partial L/\partial\dot q_i)/dt=0$ sonucunu verir. Bu durumda karşılık gelen genelleştirilmiş momentum korunur.

Benzer biçimde Lagrangian zamanın açık fonksiyonu değilse, uygun koşullarda enerji benzeri büyüklük korunur. Standart $L=T-V$ biçiminde ve zaman bağımlı olmayan korunumlu sistemlerde bu büyüklük genellikle $E=T+V$ olarak yorumlanır. Ancak zamanla değişen dış etkiler veya enerji dağıtan kuvvetler varsa bu korunum sonucu doğrudan kullanılamaz.

Noether teoremi, sürekli simetriler ile korunum yasaları arasındaki ilişkiyi sistematikleştirir: zaman öteleme simetrisi enerji, uzay öteleme simetrisi momentum, dönme simetrisi ise açısal momentum korunumu ile ilişkilidir. Bu ifadelerin kullanılabilmesi için ilgili simetrinin gerçekten Lagrangian tarafından korunması gerekir.

Lagrange yöntemi, Newton yönteminin fiziksel olarak farklı bir yasası değildir; uygun koşullarda aynı hareket denklemlerini verir. İdeal kısıtlar ve uygun koordinatlar için kısıt kuvvetleri denklemden çıkabilir. Fakat sürtünme gibi ideal olmayan kuvvetler, bağımlı koordinatlar veya açık kısıt kuvvetlerini bulma amacı varsa ek yöntemler gerekebilir. Lagrange çarpanları, kısıt denklemlerini koruyarak kısıt kuvvetlerini hesaplamada kullanılabilir.

Adım adım çözümlü örnekler

Örnek 1 — Düşey yay-kütle sistemi:

Verilenler: Kütle $m$, yay sabiti $k$, düşey yer değiştirme koordinatı $x$. Koordinat, yayın doğal uzunluğundan itibaren ve yukarı yön pozitif seçilsin. Yayın uzama enerjisi $V_y=\frac{1}{2}kx^2$; yerçekimi potansiyeli ise sabit bir referans farkı dışında $V_g=mgx$ biçiminde yazılabilir.

1. Kinetik enerjiyi yaz: $T=\frac{1}{2}m\dot x^2$.
2. Potansiyel enerjileri topla: $V=\frac{1}{2}kx^2+mgx$.
3. Lagrangian'ı kur: $L=\frac{1}{2}m\dot x^2-\frac{1}{2}kx^2-mgx$.
4. Türevleri al: $\partial L/\partial\dot x=m\dot x$ ve $\partial L/\partial x=-kx-mg$.
5. Euler-Lagrange denklemini uygula: $m\ddot x-(-kx-mg)=0$.

Sonuç: $m\ddot x+kx+mg=0$. Denge konumu $x_0=-mg/k$'dır. $y=x-x_0$ ile denge etrafındaki koordinata geçilirse denklem $m\ddot y+ky=0$ olur. Bu, yerçekiminin denge noktasını değiştirdiğini; yay-kütle salınımının biçimini ise denge etrafında yazıldığında belirlediğini gösterir.

Örnek 2 — Düzlemde merkezi kuvvet için koordinat seçimi:

Verilenler: Kütlesi $m$ olan parçacık, düzlemde konuma bağlı korunumlu bir merkezi potansiyel $V(r)$ içinde hareket ediyor. Kartezyen koordinatlar yerine kutupsal koordinatlar $r$ ve $\varphi$ seçilsin.

1. Hızı yaz: $v^2=\dot r^2+r^2\dot\varphi^2$.
2. Kinetik enerjiyi kur: $T=\frac{1}{2}m(\dot r^2+r^2\dot\varphi^2)$.
3. Lagrangian'ı yaz: $L=\frac{1}{2}m(\dot r^2+r^2\dot\varphi^2)-V(r)$.
4. Açısal koordinatın durumunu incele: $L$ içinde $\varphi$ açıkça görünmüyor; yalnızca $\dot\varphi$ var. Bu nedenle $\partial L/\partial\varphi=0$.
5. Korunan büyüklüğü belirle: $p_\varphi=\partial L/\partial\dot\varphi=mr^2\dot\varphi$ sabittir.

Sonuç: Merkezi potansiyelde açısal momentumun düzlemsel karşılığı $mr^2\dot\varphi$ korunur. Buradaki çıkarım, potansiyelin yalnızca $r$'ye bağlı olması koşuluna dayanır; $V$ açıya da bağlı olsaydı $\varphi$ çevrimsel koordinat olmayacak ve bu korunum doğrudan elde edilemeyecekti.