Lagrange mekaniği

**Lagrange mekaniği**, klasik mekaniğin yeniden formüle edilmesidir. İtalyan-Fransız matematikçi ve astronom Joseph-Louis Lagrange tarafından 1788'de geliştirilmiştir.

Lagrange mekaniğinde, parçacıkların ve parçacık sistemlerinin hareketleri ve gezingeleri *Lagrange denklemlerinin* iki formundan birinin çözülmesiyle bulunur. *Birinci tip Lagrange denklemleri*, parçacıkların hareketlerindeki kısıtlamaların, Lagrange'in kalkülüste de sık kullanılan Lagrange çarpanı yöntemi kullanılarak ifade edilmesidir. *İkinci tip Lagrange denklemleri*nde ise hareket kısıtlarını, sistemin yapısına uyacak makul genelleştirilmiş koordinat seçimi ile direkt olarak sistem denklemlerine dahil etmek üzerine kuruludur. İki durumda da matematiksel bir fonksiyon olarak tanımlanan, Lagrange'in sistem hakkında verdiği bilgilerden faydalanılır. **Lagrangian**, genelleştirilmiş koordinatların, onların türevlerinin (genelleştirilmiş hızlar) ve zamanın bir fonksiyonudur.

Lagrange mekaniğinde dikkat edilmesi gereken husus; yeni bir fiziğin ortaya atılmıyor oluşudur. Keza, Lagrange mekaniği Newton mekaniğinden daha sığ bir kapsama sahiptir çünkü Newton mekaniği korunumlu olmayan sürtünme gibi kuvvetleri de içerir. Ancak bu tip kuvvetlerin doğru ifade edilmesi için hareketi kısıtlayan kuvvetlerin de sistem denklemlerinde ifade edilmesi gerekir ve kartezyen koordinat sistemi haricindeki koordinatlar bu yönteme pek uygun değildir. Ancak, Lagrange mekaniğindeki yöntemler ile kısıtlayıcı kuvvetlerin ve korunumsuz kuvvetlerin ekstra karmaşık denklemler yaratması önlenebilir. Lagrange mekaniği bu açıdan, korunumlu kuvvetlerin etki ettiği ve iç kuvvetler ile ilgilenmediğimiz sistemlerde çok etkilidir ve kartezyen koordinat harici istediğimiz koordinat sistemini kullanmamızda bize kolaylık sağlar. Ayrıca bazı korunumsuz kuvvetleri de göz ardı etmemize olanak tanıyabilir. Dolayısıyla, kontak noktalarındaki ve yapısal sistemlerle ilgilenilmeyen ve sadece genel dinamiğin ve hareketin önemli olduğu sistemlerde Lagrange mekaniği kullanımı çok etkilidir.

Genelleştirilmiş koordinatlar sistemin serbestlik derecesine göre istenildiği gibi seçilebilir. Ancak Lagrange mekaniğinin kolaylıklarından etkin bir şekilde faydalanmak için, sistemdeki simetrileri ve geometrik kısıtları ve özellikleri dikkate alarak uygun seçimler yapmak gerekir. Böylece matematiksel model ve denklemler önemli oranda basitleşir.

Lagrange mekaniği uygulama açısından kolaylıklarının yanında, teorik olarak doğadaki korunum yasaları yani fiziksel işleyiş hakkında da bize önemli bilgiler verir. Korunan her fiziksel değişken ve onların simetrisi üzerine kurulu olan bu bilgi Noether teoremi ile daha da genelleştirilmiştir. Bu teori aynı zamanda minimum eylem prensibi ile de yakından ilişkilidir. Yine, Lagrange mekaniği, sadece denge (dinamik-statik-kaotik) problemlerine odaklandığından daha sığ kalmaktadır.  Ayrıca, Lagrange mekaniği, sadece holonomik kısıtları olan sistemlere uygulanabilmektedir çünkü formülasyon holonomik olmayan kısıtlar için çalışmamaktadır. Örneğin, integrallenebilir olmayan kısıtlayıcı denklemler (hareket kısıtları hakkındaki denklemler), Coulomb sürtünmesi gibi doğrusal olmayan karmaşık korunumsuz kuvvetler ve aralıklar ile eşitsizliklerle tanımlanmış hareket kısıtlarını tanımlayan denklemlere sahip sistemlerde Lagrange mekaniği etkin olamamaktadır. Holonomik olmayan sınırlamalar özel analiz gerekmektedir. Bu durumlarda Newton mekaniği veya başka yöntemler kullanılabilir zira bu durumlarda Lagrangian'in temelini oluşturan enerji terimlerinin ifadelerinde sıkıntılar yaşanmaktadır. Çünkü bu terimler kuvvet ve yola bağlı terimlerin integralleri ile ifade edilmektedir.

Daha önce de bahsedildiği gibi, Lagrange formülasyonu pratikteki uygulama kolaylıklarının yanı sıra, mekanik fiziği ve doğanın mekanik davranışını anlamamızda bize derin bir bakış açısı sağlar. Kuantum mekaniğinde, Planck sabiti ile ilişkili olan etki ve kuantum fazlarda, minimum eylem prensibi dalga yayılımı ve dalga yayılım denklemleri ile anlaşılabilmekte ve benzer mekanik yaklaşımlar bu gibi apayrı alanlarda da türetilebilmektedir. Lagrangian eğer simetri altında değişmez ise, sonuçta ortaya çıkan hareket denklemi de simetri altında değişmezdir. Bu karakteristik de özel görelilik ve genel görelilik teorisindeki açıklamalar için çok önemlidir. Minimum eylem prensibi, Lagrange formülasyonları gibi konular Noether teoremi ile yakından ilişkilidir ve bu teorilerin sonucunda daha önce sadece gözleme veya basit derivasyonlara dayalı kalmış olan korunum yasaları ve korunan nicelikler; inceleme altındaki fiziksel sistemin sürekli simetrileri ile ilişkilendirilebilmiştir. Dolayısıyla, bu algı; momentum, kütle gibi temel korunum sahibi niceliklerin ötesinde her sistem için farklı ve daha genel korunum ifadeleri yazabileceğimizi ortaya atmış ve bunların nasıl türetileceği, hangi durumlar altında kullanılabileceği hakkında bize bilgi vermiştir. Lagrange mekaniği ve Noether teoremi ayrıca komutatörleri içeren birinci kuantizasyonlar için doğal bir biçimcilik ortaya koymuştur ve bunu belli terimler ile Lagrangian ile türetilmiş hareket denklemleri arasında göstermiştir.

Lagrange mekaniği, Newton mekaniğinin uygun olmadığı, çok vakit aldığı veya etkili sonuç çıkartamadığı yerlerde hem teorik fizikte hem de pratik mekanik uygulamalarda çok etkilidir. Ayrıca minimum eylem prensibi ile yakından ilişkili olduğundan optimizasyon problemlerinde de kullanılabilmektedir. Mekanik biliminde, ikinci tip Lagrange denklemleri birinci tipe göre çok daha sık kullanılmaktadır.

## Giriş

Lagrange mekaniğinin gücü, hareketi kısıtlanmış mekanik sistemlerde, bu kısıtları etkili ifade ederek sistem çözümüne dahil edebilmesinden gelir. Aşağıdaki örnekler, konseptleri ve terminolojiyi anlamada etkili olacaktır.

Rijit bir eğimli yolu ifade eden bir telde hareket eden bir küreyi düşündüğümüzde, ağırlığın etkili olduğu 2 boyutlu uzayda; hareket kısıtları *f*(**r**) = 0 formunda ifade edilebilir. Burada kürenin hareketi **r** = (*x*(*s*), *y*(*s*)) olarak yazılabilir ve *s* eğrisel uzunluk olarak tanımlanır (her eğri, yay parçaları olarak ifade edilebileceğinden, *ds*i diferansiyel yay uzunluğu, *s*i de yay boyuncaki koordinat olarak tanımlayabiliriz). Eğride her noktadaki *x* ve *y* birbirine bağımlı olduğundan, aslında *x* ve *y* gibi iki koordinat tanımlamak yerine tek bir *s* koordinatı hareketimizi tanımlamamıza ve kısıtlarımızı ifade etmemize yetecektir. Çünkü hareket kısıtlarını ifade eden denklemler *x* ve *y* arasında bir bağıntı vermekte ve serbestlik derecesini düşürmektedir. Böylece *x* ve *y* den biri bilindiğinde öbürü bulunabilir.

Hareket kısıtlarının bir sonucu olan kısıtlayıcı kuvvetler ise, hareketin kısıtlandığı yönde uygulanan kuvvetlerdir veya başka bir açıdan, kürenin istenilen yolda kalmasını sağlayan kuvvetlerdir. Ayrıca ağırlık da, kısıtlayıcı olmayan bir kuvvet olarak küreye etki etmektedir.

Eğer telin şekli zamanla değişiyorsa (örneğin esneme ile), o zaman kısıtlayıcı denklemler ve pozisyon zamana bağlı olarak da değişmeye başlar. Çünkü koordinatlar zamanla değişmektedir. Zaman, koordinat tanımlarında örtülü olarak kalırken, kısıtlayıcı denklemlerde açık olarak gelir.

:f(\mathbf{r}, t) = 0 \,,\quad \mathbf{r} = (x(s, t), y(s, t))

Başka bir ilginç 2 boyutlu örnek ise kaotik/düzensiz bir dinamiğe sahip olan çift sarkaçtır (yine yer çekimi altında). Bir sarkaçın uzunluğu *L*1 ve diğerinin uzunluğu ise *L*2 olsun. Her bir kütle hareketlerinin kısıtlanmasını ifade eden denklemlere sahiptir.
Her sarkacın kısıtlayıcı denklemi vardır,

:f(\mathbf{r}_1) = x_1^2 + y_1^2 - L_1^2 = 0 \,,\quad f(\mathbf{r}_2) = x_2^2 + y_2^2 - L_2^2 = 0 \,, 

Kütlelerin pozisyonları ise,

:\mathbf{r}_1 = (x_1(\theta_1), y_1(\theta_1)) \,,\quad \mathbf{r}_2 = (x_2(\theta_2), y_2(\theta_2)) \,, 

olarak ifade edilebilir.

*θ* 1 ve *θ* 2 belli bir referansa göre sarkaçların açısını ifade eder. Her sarkaç, kısıtlardan ötürü birer koordinatla ifade edilebilir. Çünkü bu kısıtlar mekansal koordinatları (2 tane) birbirine bağlar.

3 boyutlu bir örnek olarak, küresel eksenlerde hareket eden sarkaç örnek verilebilir. Sarkacın uzunluğu *l* ve sarkaç her yönde serbestçe sallanacak şekilde ağırlığının etkisine bırakılmış durumdadır. Bu durumda, hareket kısıtlanması şu şekilde ifade edilebilir;

:f(\mathbf{r}) = x^2 + y^2 + z^2 - l^2 = 0 \,, 

Kütlenin pozisyonu ise daha önce s ve θ değişkenlerinin kullanıldığı gibi genel koordinatlarda şu şekilde ifade edilebilir;

:\mathbf{r} = (x(\theta,\phi),y(\theta,\phi),z(\theta,\phi)) \,, 

ve burada (*θ*, *φ*) çifti küresel kutup açılarını ifade eder çünkü kütle bir kürenin yüzeyinde hareket edermişçesine davranır aynı zamanda küresel koordinatlarda tanımlanabilir ve aynı bu durumda olduğu gibi, yarıçapın sabit olduğu bir kürenin denklemi de kısıtımızı ifade eder. Böylece 3 koordinat yerine 2 koordinat bilgisi yeterli olmaktadır.

Özellikle sarkaç gibi sistemler, salınım hareketi yaparlar çünkü yükleme, pozisyon değişkeni olan açının sinüsel bir fonksiyonudur. Böylece basit harmonik osilatörler Lagrange mekaniği ile yine rahatça analiz edilebileceği anlaşılabilir. Ayrıca, bu tip enerji metotları denk sistemler bulmak, çoklu serbestlik dereceye sahip sistemlerin analizi gibi problemlerde ciddi işe yaramaktadırlar çünkü yazılımsal olarak adapte edilebilmeleri kolaydır. Dolayısıyla, bağlanmış harmonik osilatörler (çift sarkaç, üçlü sarkaç gibi) özellikle Lagrange dinamiği metotları ile incelenir.

3 boyutlu uzayda yer alan *N* partikül için, pozisyon vektörleri 3 değişkenli vektör fonksiyonları olarak şöyle ifade edilebilir;

:\mathbf{r}_1 = (x_1,y_1,z_1) \,, \quad \mathbf{r}_2 = (x_2,y_2,z_2) \,, \ldots \,, \mathbf{r}_N = (x_N,y_N,z_N)\,, 

Dolayısıyla, her bir parçanın net pozisyonunu bilmek için 3*N* kadar koordinata ihtiyacımız vardır. Zaten bir vektör uzayında da bir noktanın genel ifadesi **r** = (*x*, *y*, *z*). fonksiyonu ile yapılır. Bu partiküllerden herhangi birisi veya birkaçı holonom kısıtlara tabi tutulursa, kısıt ifadesi *f*(**r**, *t*) = 0 şeklinde yazılabilir. Bu, herhangi bir anda partiküllerin pozisyonları birbirlerine bağlanmaya başlar ve bağımsız olmaktan çıkarlar. Eğer sistemde *C* kadar kısıtlayıcı var ise, bunların her birinin fonksiyonel ifadesi için,

:f_1(\mathbf{r}, t) = 0\,,\quad f_2(\mathbf{r}, t) = 0 \,, \ldots \,, f_C(\mathbf{r}, t) = 0 \,, 

Fonksiyon eşitlikleri kullanılır ve her bir denklem bir koordinat ihtiyacını eler. Dolayısıyla bağımsız koordinatların sayısı *n* = 3*N* − *C* kadardır. Unutulmamalıdır ki bağımlı koordinatlar kullanılarak yapılan hesaplar modeli karmaşıklaştıracaktır.

Dolayısıyla, partikülleri veya sistemi tanımlamamız için *n* tane koordinata ihtiyacımız bulunmaktadır. Her pozisyon vektörü istenilen *n* elemanlı genelleştirilmiş koordinat setine dönüştürülebilir. Ve uygun bir şekilde **q** = (*q*1, *q*2, ... *q<sub>n</sub>*) olarak ifade edilebilir. Dolayısıyla pozisyonlar şu şekilde ifade edilebilir;

:\mathbf{r}_k = \mathbf{r}_k(\mathbf{q}, t) = (x_k(\mathbf{q}, t), y_k(\mathbf{q}, t), z_k(\mathbf{q}, t), t)\,.

Böylece, pozisyonlar, genelleştirilmiş koordinatların fonksiyonu olarak karşımıza çıkar çünkü pozisyon değişkenlerini herhangi bir sete dönüştürebiliriz.

Genelleştirilmiş koordinatların birinci türevine genelleştirilmiş hızlar denir ve her partikül için hız vektörlerinin dönüşümleri;

:\mathbf{v}_k = \dot{\mathbf{r}}_k = \frac{d\mathbf{r}_k}{dt} = \left(\frac{dx_k}{dt},\frac{dy_k}{dt},\frac{dz_k}{dt}\right)\,, \quad \dot{q}_j = \frac{dq_j}{dt} \,, \quad \mathbf{v}_k = \sum_{j=1}^n \frac{\partial \mathbf{r}_k}{\partial q_j}\dot{q}_j +\frac{\partial \mathbf{r}_k}{\partial t}\,, 

(Noktalar Newton notasyonu olarak zamana göre türevi ifade eder.)

Burada unutulmamalıdır ki hız, pozisyonun toplam diferansiyelidir. Veya basitçe *r*k olarak tanımladığımız fonksiyonun zamana göre türevidir ve *q* da zamana bağlı olduğu için zincir kuralı kullanılır.

Önceki örneklerde, biri eğer her kütlenin hareketini partikül olarak inceleyecek basitleştirmeler yapabiliyorsa, Newton mekaniği zamana bağlı olarak değişen kısıtlayıcı kuvvetlerin (partikülleri kısıtlanmış hareket patikasında tutmaya yarayan kuvvetler) zamana göre değişimi ile uğraşmayı gerektirir (telin küreye uyguladığı kuvvet, çekme kuvvetinin sarkaca aktardığı kuvvet gibi).Aynı problem için Lagrange mekaniği, bizim partikülün aldığı yola bakarak uygun bir genelleştirilmiş koordinat seçmemize dayanır. Bu seçimlerdeki koordinatlar *bağımsız* olmalıdır ve hareketi tamamen karakterize edebilmelidir. Hareketin tamamen serbest olduğu koordinatlar seçilirse, kısıtlayıcı denklem ve kuvvetlerin işin içine girmesine gerek kalmaz. Örneğin *x*-*y* düzleminde tel üzerindeki bir kürenin hareketi her an tanjantlara dik olarak kısıtlanmıştır ama tek bir *s* koordinatı seçilirse, zaten o koordinat yolun ve hareketin kendisini tanımladığı için ve hareketin olmadığı yerlerde sıfır olduğu için ekstra kısıt ifadelerine gerek kalmaz. Böylece, kısıtların partiküller üzerindeki direkt etkisi hesaplanmadan daha az denklem çözülerek sistem modellenebilir.

## Lagrangian'ın tanımı
Lagrange mekaniğindeki temel nicelik **Lagrangian fonksiyonu**dur. Bu fonksiyon, sistemin dinamik davranışı hakkında önemli bilgiler barındırır. Ancak, her fiziksel sistem için tek bir Lagrangian ifadesi türetmek mümkün değildir. Esasen, Euler-Lagrange denklemlerini fiziksel yasalarla uyumlu bir şekilde türetebilecek her fonksiyon Lagrangian olarak ele alınabilir. Ancak yine de, bazı gruplanabilecek sistemler için genel Lagrangian ifadeleri oluşturmak mümkündür.

Örneğin, görelilik etkilerin ihmal edilebileceği sistemler için Lagrangian;L = T - V \,,

olarak ifade edilebilir. (mekanikte en sık karşılaşılan ifade budur, ancak bunun tüm sistemler için genel bir ifade olmadığına dikkat edilmesi gerekir. Farklı fiziksel yasaların da işin içine girmesiyle değişmesi gerektiğini unutmamak gerekir.) Burada *T*, sistemin toplam kinetik enerjisini ve *V* ise potansiyel enerjisini ifade eder. Lagrangian'ın birimi, seçilen enerji birimleriyle aynıdır.

Kütleleri *m*1, *m*2, ..., *m<sub>N</sub>* olan *N* parçacıktan oluşan bir sistem için toplam kinetik enerji her bir parçacığın kinetik enerjisi toplanarak bulunabilir. (Enerji metodlarındaki toplanabilme özelliği bu açıdan çok güçlüdür.) Her parçacığın enerjisi sistemin hareketine göre şöyle tanımlanabilir;

:T = \frac {1}{2} \sum_{k=1}^N m_k v^2_k\,,\quad v^2_k  = \mathbf{v}_k\cdot \mathbf{v}_k = \dot{\mathbf{r}}_k \cdot \dot{\mathbf{r}}_k\,,

Burada, işaretlerin doğru verilmesi için noktasal çarpım kullanılmıştır. Yani kinetik enerji ne olursa olsun pozitif bir ifadedir. Bu ifadeye göre, kinetik enerji sadece hızların yani **v***k* terimlerinin bir fonksiyonudur. Kinetik enerji **r***k*'ya yani parçacıkların pozisyonlarına ve zamana bağlı değildir. Dolayısıyla *T*(**v**1, **v**2, ...) olarak ifade edilebilir. Ancak, yukarıdaki genelleştirilmiş koordinat ve hız ifadeleri göz önüne alındığında, özellikle zamana bağlı bir değişim söz konusu ise; genelleştirilmiş koordinatlar çerçevesinde ifade edilen kinetik enerji genelleştirilmiş hızlara, koordinatlara ve pozisyon vektörleri zamana bağlı çarpanlar içeriyorsa zamana bağlıdır. Yani *T* = *T*(**q**, *d***q**/*dt*, *t*) olarak daha doğru bir şekilde ifade edilebilir.

Potansiyel enerji *V* ise sistemdeki parçacıkların birbirleri ile olan ilişkilerini yansıtır. Başka bir deyişle, bir parçacığın diğer tüm parçacıklara ve dış etkilere göre ne kadar enerjisi olduğunu ifade eder. korunumlu kuvvetlerin etkisinden ötürü oluşan potansiyel enerji, sadece parçacıkların pozisyonlarına bağlıdır yani *V* = *V*(**r**1, **r**2, ...)’dir. Korunumlu olmayan kuvvetler için de uygun potansiyel enerji fonksiyonları türetilebilir. Ancak bu durumda işin içine hızlar da girmek zorundadır. Yani korunumsuz kuvvetler için *V* = *V*(**r**1, **r**2, ..., **v**1, **v**2, ...) olarak ifade edilebilir. Ayrıca dış bir fiziksel alan veya kuvvet zamana bağlı olarak değişiyorsa potansiyel enerji de bundan etkilenecektir. Dolayısıyla en genel haliyle; *V* = *V*(**r**1, **r**2, ..., **v**1, **v**2, ..., *t*) olacaktır.

Eğer potansiyel ve kinetik enerjilerden biri veya ikisi de zamana açık bir bağlılık gösteriyorsa (zamanla değişen kısıtlar, etkiler, çarpanlar veya sistem parametreleri) Lagrangian da zamana açık bir şekilde bağlı olacaktır. Yani *L*[**q**,*d***q**/*dt*,*t*) olacaktır. Eğer potansiyel enerji ve kinetik enerji zamana göre değişmiyorsa bu durumda, Lagrangian, *L*(**q**, *d***q**/*dt*) formunda olacaktır ve zamana açık bir bağlılık göstermeyecektir. Ancak zamana göre türev alındığı için zamana örtülü bir şekilde bağlı olacaktır. İki durumda da Lagrangian’ın zamana örtülü olarak bağlı olacağı kesindir (Genelleştirilmiş koordinatlar kullanılırsa).

Lagrangian’a ek olarak ayrıca, disipatif yani hareket süresince enerjiyi o an sistem içinde dönüştürülemez şekilde harcayan sistem elemanlarını göz önüne alarak çeşitli fonksiyonlar türetmek de mümkündür.

Göreli-mekanikde, *L* daha derin bir şekilde analiz edilip doğru ifade edilmelidir. Bunun için relativistik Lagrange mekaniğine bakılabilir. (Enerji ifadeleri artık daha farklı olacak ve dönüşümler, zamana bağımlılık farklı bir yapı gösterecektir.)

## Hareket denklemleri
Newton mekaniğinde, hareket Newton Yasaları ile ifade edilir. Newton'un ikinci yasası **F** = *m***a** her parçacığa uygulanır ve 3 boyutta hareket sahibi olan*N* parçacıklı bir sistem için 3*N* tane adi diferansiyel denklem pozisyonların bulunması için çözülür(ivmenin iki kere integrali alınarak).

Lagrange mekaniğinde ise, Lagrangian fonksiyonu hareket denklemlerini türetir. Hareket denklemleri, kısmi türev  içerse bile, denklemler yine de genelleştirilmiş pozisyon koordinatlarında çalışan adi diferansiyel denklemlerdir. Ayrıca zaman göre toplam türev ifadesi,*d*/*dt* örtülü türev alınmasını gerektirebilir.

İki çeşit hareket denklemi vardır. Biri kısıtlayıcı kuvvetleri ve hareketlerin kısıtlarını ifade eden denklemlerini işin içine katarken diğeri bunları modellemenin dışında tutarak sadece kısıtlayıcı olmayan ve hareket sağlayan kuvvetleri genelleştirilmiş koordinatlarda inceler. Bu ifadeleri Lagrangian’a koyup gerekli denklemleri çözerek, sınır ve ilk değerleri katarak (hız ve pozisyon) parçacıkların pozisyonlarını zamana göre çözebilir ve sistemin zamana göre hareketini, değişimini bulabiliriz.

### Birinci tip Lagrange denklemleri
3 boyutlu uzayda *N* parçacıktan oluşan ve *C* tane holonomik kısıta maruz kalan bir sistemin kısıtları yukarıda belirtildiği üzere, *f*2,..., *f<sub>C</sub>* gibi fonksiyonlarla ifade edilebilir ve sistemin dinamiği, Lagrangian *L*(**r**, *d***r**/*dt*, *t*)'a bağlı olarak **birinci tip Lagrange denklemleri** ile şöyle ifade edilebilir;

\frac{\partial L}{\partial \mathbf{r}_k} - \frac{d}{dt}\frac{\partial L}{\partial \dot{\mathbf{r}}_k} + \sum_{i=1}^C \lambda_i\frac{\partial f_i}{\partial \mathbf{r}_k }=0 

burada *k* = 1, 2, ..., *N* her bir parçacığı ifade eder ve bu denklemler sisteme dahil olan her eleman veya parçacık için çözülmek üzere bir denklem sistemi oluşturur. Ayrıca, her bir kısıtlayıcı denklem *f<sub>i</sub>* için birer Lagrange çarpanı *λ<sub>i</sub>* bulunur. Ve bu indisler her bir vektörün her bir ilgili değişkene göre türevini ifade etmede kolaylık sağlar.

:\frac{\partial}{\partial \mathbf{r}_k} \equiv \left(\frac{\partial}{\partial x_k},\frac{\partial}{\partial y_k},\frac{\partial}{\partial z_k}\right)\,,\quad \frac{\partial}{\partial \dot{\mathbf{r}}_k} \equiv \left(\frac{\partial}{\partial \dot{x}_k},\frac{\partial}{\partial \dot{y}_k},\frac{\partial}{\partial \dot{z}_k}\right)

Burada dikkat edilmesi gereken, Lagrangian'ın her bir vektöre göre türevi alınmasından ziyade yine mekansal koordinatlara göre türev alındığı ve buradaki gösterimin kolaylık olduğudur.

:\frac{\delta}{\delta \mathbf{r}_k} \equiv \frac{\partial}{\partial \mathbf{r}_k} - \frac{d}{dt}\frac{\partial}{\partial \dot{\mathbf{r}}_k} \,

Bu prosedür Newton yasalarıyla çözülebilecek 3*N* denklem sayısını kısıtlayıcı denklemlerin de işin içine girmesiyle 3*N* + *C*ye çıkartır. Çünkü *N* parçacık için 3 boyuttta – 3 pozisyonda 3*N* denklem vardır ve *C* kadar hareket kısıtlarını ifade eden denklem vardır. Ancak burada etkili olan nokta Lagrange çarpanlarıdır ve bu çarpanlar kısıtlayıcı kuvvetler hakkında bize çözüm sırasında bilgi verir. Bu metotta genelleştirilmiş koordinatlar veya kısıt denklemlerinin temel hareket denklemlerine katılması kullanılarak denklemleri elimine etme ihtiyacı yoktur. Kısıt denklemleri ekstra denklemler olarak Lagrange çarpanları ile ifade edilir ve sistem çözülür.

## İkinci tip Lagrange denklemleri

**Euler-Lagrange denklemleri** veya Lagrange'ın **ikinci tip hareket denklemleri** olarak adlandırılır. Değişim kalkülüsünün matematiksel bir sonucu olarak ortaya konmuştur. Bu sonuç fiziksel – mekanik sistemlerde de hem doğalarını anlamada hem de hareketlerini hesaplamada kullanılır.

\frac{d}{dt} \left ( \frac {\partial  L}{\partial \dot{q}_j} \right ) =  \frac {\partial L}{\partial q_j} 

Lagrangian'ın *L*(**q**, *d***q**/*dt*, *t*) formulasyonu bu denkleme konduğunda her bir genelleştirilmiş koordinata karşılık bir denklem bulunur. Ve genelleştirilmiş koordinatlar eksenel koordinatlar – kısıtlar sayısı kadar seçildiğinden toplamda 3*N*-*C* kadar 2. Dereceden adi diferansiyel denklemi çözmemiz gerekir. Bu denklemler hareket kısıtlarını içermezler ve sadece kısıtlayıcı olmayan – hareketi gerçekleştiren kuvvetlerin modele katılması gerekir.

Bu metot, değişim kalkülüsündeki minimum değişim ve hareket prensibinin fiziksel bağlamda enerji olarak ifade edilebilecek bir değişim fonksiyonuyla etkili bir şekilde kullanımı ile ortaya çıkar.

## Newton mekaniğinden Lagrange mekaniğine

### Newton yasaları

Basitçe, Newton yasaları bir parçacık için en genel bir şekilde;

: \mathbf{F} = m \mathbf{a} \,, 

olarak ifade edilir. Burada *a*, ivme; pozisyonun 2. Türevidir. Bu sayede *N* parçacığa sahip bir sistemde her parçacığa bu prensip uygulanarak hareket denklemleri çıkartılır. *F* ise, parçacığa etkiyen toplam kuvvettir.

Üç mekansal koordinatta her bir parçacık için üç hareket denklemi ortaya çıkar. Ve bu denklemlerin her biri adi diferansiyel denklemdir. 3 boyutlu uzayda üç değişken içerecek sonuçlar, birer vektör olarak ifade edilebilir. Elbette, her bir çözüm ilk hız, ilk pozisyon gibi başlangıç değerlerine tabidir.

Newton yasalarının kartezyen koordinatlarda kullanımı çok rahattır. Ancak çoğu harekette (özellikle de eğrisel hareketlerde) kartezyen koordinat kullanımı zordur. Bu durumda kartezyen koordinatların uygun koordinatlara çevrilmesi ve denklemlerin bu şekilde ifade edilmesi gerekir. Ancak bu durumda da pozisyonların türevlenmesi gerektiği için hareket denklemleri karmaşıklaşır. Örneğin üç eğrisel koordinatta  **ξ** = (*ξ*1, *ξ*2, *ξ*3), Newton yasası *"Lagrange formunda"* şöyle ifade edilebilir;

: F^a = m \left( \frac{d^2 \xi^a}{dt^2} + \Gamma^a {}_{bc} \frac{d\xi^b}{dt}\frac{d\xi^c}{dt} \right) = \frac{d}{dt} \frac{\partial T}{\partial \dot{\xi}^a} - \frac{\partial T}{\partial \xi^a} \,, \quad \dot{\xi}^a \equiv \frac{d \xi^a }{dt} \,, 

Burada *F<sup>a</sup>*, kartezyen uzaydan eğrisel uzaya geçişte, uzayların taban vektörlerinin değişmesinden ötürü değişen kuvvetlerinkarşıdeğişkeni *a* bileşenidir. Γ*<sup>a</sup>* ise ikinci tip Christoffel sembolleridir. Burada bir parçacığın kinetik enerjisi;

:T = \frac{1}{2} m g_{bc} \frac{d \xi^b}{dt} \frac{d \xi^c}{dt} \,, 

olarak ifade edilebilir. *g<sub>bc</sub>*, eğrisel koordinat sisteminin *[[metrik tensör]]*ün eşdeğişken bileşenidir. Her  *a*, *b*, *c*, indisi 1,2,3 gibi değerler alır. Burada dikkat edilmesi gereken husus, eğrisel hareketin önemi vurgulandığından; bu koordinat sisteminin genelleştirilmiş koordinatlarla karıştırılmaması gerektiğidir.

Burada yapılan türetmeler Newton yasalarının fazlaca karmaşıklaştırılması olarak düşünülebilir. Ancak bu yöntemin de çeşitli avantajları bulunmaktadır. Christoffel sembolleri cinsinden ifade edilen ivme bileşenleri, yukarıda tanımlanan kinetik enerjinin türevleri alınarak ifade edilebilir ve böylece  *T* gibi bir her referans yapısında aynı kalacak skalerin kolaylığı kullanılabilir. Zaten Lagrangian gibi enerji metodları da vektörlerin referanslara bağımlılığını elimine ederek işlemleri kolaylaştırma ihtiyacına dayanır (Bunu yaparken ise sistemin detaylarından feragat etmek durumunda kalırlar).

Eğer sistem üzerine etkiyen bir kuvvet yoksa, yani **F** = **0** ise, sistem ivmelenemez. Ya durur ya da hareketini değiştirmeden devam eder. Hareket durumu daha genel bir ifadeyle; doğru üzerinde sabit hızlı harekettir. Ancak bu, bildiğimiz anlamda doğru olmak zorunda değildir. Örneğin eğrisel koordinatlara geçtiğimizde (jeodezik koordinatlar) doğrular artık keseller olarak karşımıza çıkar. Dolayısıyla, eğrisel yüzeyler üzerinde hareket etmek üzere kısıtlanmış veya o şekilde dengeye gelmiş parçacıkların hareket denklemlerinin çözümü keselleri verir. *Keseller* uzayda iki genel nokta arasındaki en kısa uzunluğa karşılık gelen eğrilere verilen genel addır. Düzlemde ise bunlar basitçe doğrulardır. Dolayısıyla 3 boyutlu uzayda da keseller düz doğrulara karşılık gelirler. Dolayısıyla, serbest bir parçacık için, Newton'un 2. Yasası geometrik jeodezik denklemleri ile uyum halindedir ve keselleri takip eden serbest parçacıkların hareketlerini ifade etmede kullanılır. Burada en önemli nokta, kesellerin 3 boyutlu uzayda parçacıkların iki nokta arasında izleyeceği en kısa yolu ifade etmesidir. Eğer parçacıklar belli bir kuvvet altında kalırlarsa ivmelenirler ve bu nedenle jeodezik hareketten sapmaya başlarlar.

Uzayda noktalar arası en kısa yolları bulma fikri mekanikte doğal bir merak konusudur. Özellikle enerjinin tanımından sonra mekanik sistemlerin minimum enerjiyle amaçlarını gerçekleştirmeleri önem kazanmıştır. Bu nedenle, değişimler hesabı (kalkülüs) metotlarının mekaniğe uygulanması ile yukarıda bahsedilen sonuçlar elde edilmiştir. *Brachistochrone problemi* (sabit ivme altında en kısa yol) gibi önemli konular değişimler hesabı sayesinde açıklığa kavuşturulmuştur. [Jean Bernoulli)), [Gottfried Leibniz|Leibniz, Daniel Bernoulli, L’Hôpital gibi bilim insanları bu konu üzerine çalışmış ve farklı zamanlarda benzer çözümler bularak katkı sağlamışlardır. Newton'un kendisi bile değişimler hesabı prensiplerini o zamanlarda düşünmeye başlamıştı. Hatta yukarıda bahsedilen eğri uzunluğu gibi kavramları tanımlayıp, onların türevlerine bakarak düzlemde en kısa yolu ifade eden doğruları eğrilere nasıl adapte edebileceğini düşünmekteydi. Bu fikirler üzerine varyasyonel prensipler Fermat, Maupertuis, Euler ve Hamilton gibi bilim adamları tarafından net bir şekilde ortaya kondu. Ayrıca, uzay-zaman gibi 4 boyutlu uzaylara da uygun şekilde genelleştirilerek Einstein'ın genel görelilik teorisine katkı sağladı. Bahsedildiği gibi uzay eğrisel olarak, zamanda büküldüğünde oluşan uzay zaman boyunca; serbest parçacıklar, düzlemde doğruya karşılık gelecek şekilde jeodezik hareket yaparlar.

Ancak, Newton mekaniğinde, parçacıkların hareketlerini tam olarak modelleyebilmek için üzerlerine etkiyen kuvvetleri **F** bilmek gerekir. Genel bir harekette bu kuvvet, hareketi çeşitli eksenlerde kısıtlayan veya kısıtlamayan kuvvetlerin birleşimidir. Kısıtlayıcı kuvvetleri **C**, kısıtlayıcı olmayan ve hareketi direkt sağlayan kuvvetlere ise **N** dersek;

: \mathbf{F} = \mathbf{C} + \mathbf{N} \,. 

olarak ifade edebiliriz.

Kısıtlayıcı kuvvetler karmaşık olabilir ve yola göre değişkenlik gösterebilir. Kısıtlanan mekanik sistemlerde genelde hareket yolu olabildiğinden çoğunlukla zamana (veya harekete) bağlı olurlar. Bu durumlarda eğrisel koordinatlar bağımsız değillerdir ve çeşitli kısıt denklemleri ile bağlıdırlar.

### D'Alembert prensibi

Analitik mekaniğin ilk temel sonucu olan D’alambert prensibi 1708'de Jacques Bernoulli tarafından ortaya statik dengeyi anlamak için atıldı ve D’Alambert tarafından 1743'te dinamik problemlerini çözmek için geliştirildi.

Bu prensip, N parçacık için şu denklemi tanımlar;

:\sum_{k=1}^N ( \mathbf {N}_k + \mathbf {C}_k - m_k \mathbf{a}_k )\cdot \delta \mathbf{r}_k = 0\,.

Burada, δ**r***k*,sanal yer değiştirmeleri tanımlar. *[[Sanal yer değiştirme]]* ise, sistem konfigurasyonundaki kısıtlar (hareket ve kuvvet bazında) ile uyumlu olan, diferansiyel (çok küçük) değişimlere denir. Değişimler çok küçük olduğundan, *r*, kısıtlanmış hareket eksenlerinde tanımlansa bile kısıtlayıcı kuvvetlerin etkisini bozmaz. Bu sayede her tip hareket incelenebilir ve denklemlere katılabilir. Burada kısıtlayıcı veya değil her kuvvetin sanal yer değiştirmeler boyunca yaptığı işlere sanal iş adı verilir.

Kısıtlayıcı kuvvetler her parçacığın hareketine dik yönde etki eder. (En genel haliyle uzaylarda hareketi tanımlamak için *n* boyutlu uzay ise *n* kadar birbirine dik vektör ifadesi gerektiği için oluşacak tüm kısıtlayıcı kuvvetlerin bileşkesi harekete dik yönde olacaktır. Böylece o yöndeki vektör boyunca gerçekleşecek hareket iptal olacaktır.) Kısıtlayıcı kuvvetler tarafından gerçekleştirilen toplam sanal iş bu nedenle hareket boyunca sıfırdır (Noktasal çarpımdan ötürü).\sum_{k=1}^N \mathbf {C}_k \cdot \delta \mathbf{r}_k = 0\,,

dolayısıyla;

:\sum_{k=1}^N (\mathbf {N}_k - m_k \mathbf{a}_k ) \cdot \delta \mathbf{r}_k = 0\,.

eşitliği elde edilir.

Bu nedenle, D'Alambert prensibi bize kısıtlayıcı olmayan kuvvetlere odaklanmamızı sağlar ve kısıtlayıcı kuvvetleri hareket denklemlerinden çıkartır.  Hareket denklemlerinin bu formu aynı zamanda birbirinden bağımsız olarak sistemle uyumlu genelleştirilmiş koordinat seçimine de direkt uyum sağlar. Ancak, yine de bu denklemleri her parçacık ve genelleştirilmiş koordinat için ifade edip bağıl bir şekilde çözmek gereklidir.

### D’Alambert prensibiden çıkartılan hareket denklemleri

Eğer *k* parçacığının üzerinde kısıtlar varsa, kartezyen uzaydaki  **r***k* = (*x<sub>k</sub>*, *y<sub>k</sub>*, *z<sub>k</sub>*) pozisyonundaki koordinatlar kısıt denklemleri ile birbirine bağlanır. Burada sanal yer değiştirmeler; *δ***r***k* = (*δx<sub>k</sub>*, *δy<sub>k</sub>*, *δz<sub>k</sub>*) olarak tanımlanabilir. Bu koordinatları hareket denklemleri içerisinde kısıt denklemleri ile uğraşmamak için genelleştirilmiş koordinatlara sanal yer değiştirme olarak dönüştürebiliriz. Sanal yer değiştirmeler diferansiyel boyutta olduğu için toplam diferansiyel formunu kullanabiliriz. Yani türevleri kullanarak yine işimizi analitik olarak kolaylaştırabiliriz.\delta \mathbf{r}_k = \sum_{j=1}^n \frac {\partial \mathbf{r}_k} {\partial q_j} \delta q_j \,.

Bu sayede direkt *δ***r***k*yı kullanarak yaşayacağımız sorunları, genelleştirilmiş koordinatları kullanarak ve sanal küçük hareketlerin aslında diferansiyel (türev ile ilintili) olduğu gerçeğini kullanarak hareket denklemlerini çıkarmayı kolaylaştırabiliriz.

Unutulmamalıdır ki, burada zamana göre kısmi türevin yine küçük bir zaman artışıyla çarpımı yoktur. Burada artağan bir zamandan ziyade anlık bir zamandaki ufak-diferansiyel yer değiştirmeler söz konusudur dolayısıyla türevlemeler ve artışlar yer değiştirmeler cinsindendir.

D'alambert prensibinin birinci terimi, kısıtlayıcı olmayan  **N***k* kuvvetlerinin diferansiyel *δ***r***k* boyunca yaptığı sanal işi kapsar. Ve genel olarak bu terim genelleştirilmiş kuvvetler analojisiyle ifade edilebilir (genelleştirilmiş koordinatlar boyunca etkiyen);

:Q_j = \sum_{k=1}^N \mathbf {N}_k \cdot \frac {\partial \mathbf{r}_k} {\partial q_j} \,,

dolayısıyla,

:\sum_{k=1}^N \mathbf{N}_k \cdot \delta \mathbf{r}_k = \sum_{k=1}^N  \mathbf {N}_k \cdot \sum_{j=1}^n \frac {\partial \mathbf{r}_k} {\partial q_j} \delta q_j = \sum_{j=1}^n Q_j \delta q_j \,.

şeklinde ifade edilir.

Burada yapılan genelleştirilmiş koordinatlara yapılan dönüşümün yarısını içerir. Geriye, ivmeleri içeren terimin genelleştirilmiş koordinatlara dönüştürülmesi kalır. Bu hemen görülebilir değildir. Yukarıda bahsedilen Newton Yasalarının Lagrange formu, kinetik enerjinin genelleştirilmiş koordinatlara ve hızlara göre kısmi türevini kullanarak istenilen sonucu verebilir;\sum_{k=1}^N m_k \mathbf{a}_k \cdot \frac {\partial \mathbf{r}_k}{\partial q_j} = \frac{d}{dt}\frac{\partial T}{\partial \dot{q}_j} - \frac{\partial T}{\partial q_j} \,.

Bu sayede D’alambert prensibini istediğimiz şekilde genelleştirilmiş koordinatlarda şu şekilde yazabiliriz;

:\sum_{j=1}^n \left[ Q_j - \left(\frac{d}{dt}\frac{\partial T}{\partial \dot{q}_j} - \frac{\partial T}{\partial q_j} \right) \right] \delta q_j = 0 \,, 

Dolayısıyla artık sanal yer değiştirmeler genelleştirilmiş koordinatlar cinsinden yapılabilir. *δq<sub>j</sub>*ler birbirlerinden bağımsız ve sıfırdan farklıdır. Çünkü herhangi bir zaman anında belirli bir yer değiştirmeden bahsediliyor. Bu nedenle *δq<sub>j</sub>*ların katsayıları sıfıra eşitlenebilir. Bu katsayılar da özellikle statik durumlar için genelleştirilmiş kuvvetler olarak adlandırılabilir. Bu sayede, esasen **Lagrange denklemleri**ni veya **genelleştirilmiş hareket denklemleri**ni elde etmiş oluyoruz.

:Q_j = \frac{d}{dt}\frac{\partial T}{\partial \dot{q}_j} - \frac{\partial T}{\partial q_j} 

Lagrange Denklemleri; Euler-Lagrange Denklemleriyle karıştırılmamalı, burada enerji terimleri yerine genelleştirilmiş kuvvet terimleri yer alır.

Bu denklemler, kısıtlayıcı olmayan kuvvetler altında bize Newton'un yasaları ile aynı sonucu verir. Yukarıdaki denklemlerdeki genelleştirilmiş kuvvetler sadece *kısıtlayıcı olmayan kuvvetler*den türetilmiş ve kısıtlayıcı kuvvetler dışarıda D'Alambert'in sanal iş prensibi sayesinde bırakılmıştır ve ihtiyaç yoktur. D'Alambert prensibi kuvvetlerin korunumluluğu üzerinde geneldir. Bu kuvvetler D'Alambert prensibini sağladığı sürece korunumsuz olabilir.

### Hamilton prensibi ile Euler-Lagrange denklemleri

Korunumlu olmayan kuvvetler, pozisyon ve zamana ek olarak, korunum kavramının tanımı gereği hıza da bağlıdır. Dolayısıyla, pozisyona ve hızlara bağlı bir potansiyel enerji fonksiyonu tanımlamak mümkündür. Eğer genelleştirilmiş kuvvetler bir *V* potansiyelinden şu şekilde türetilebiliyorsa;Q_j = \frac{d}{dt}\frac{\partial V}{\partial \dot{q}_j} - \frac{\partial V}{\partial q_j} \,,

ve Lagrangian'ı *L* = *T* − *V* olarak tanımlayıp bunu Lagrange denklemlerine eşitlersek ikinci tip **Lagrange denklemleri**ni veya **Euler-Lagrange denklemleri**ni  elde ederiz;

: \frac{\partial L}{\partial q_j} - \frac{d}{dt}\frac{\partial L}{\partial \dot{q}_j}  = 0 \,. 

Ancak, Euler-Lagrange denklemlerinde korunumsuz kuvvetler, sadece onları ifade edebilecek bir potansiyel fonksiyonu bulunursa ifade edilebilir. Dolayısıyla Lagrange denklemleri potansiyel içermediği, kuvvet içerdiği için Euler-Lagrange denklemlerinden daha geneldir. Çünkü Euler-Lagrange denklemlerinde kuvvetleri korunumlu veya korunumsuz; potansiyel fonksiyonları olarak ifade etmeye çalışıyoruz. Ve böyle bir potansiyel bulunamazsa EL denklemleri o kuvvetleri kapsayamıyor.

Euler-Lagrange denklemleri değişim hesabı metotlarıyla bulunabilir. *Lagrangian fonksiyonunun değişimi* şöyle ifade edilebilir;

:\delta L = \sum_{j=1}^n \left(\frac{\partial L}{\partial q_j} \delta q_j + \frac{\partial L}{\partial \dot{q}_j} \delta \dot{q}_j \right) \,,\quad \delta \dot{q}_j \equiv \delta\frac{dq_j}{dt} \equiv \frac{d(\delta q_j)}{dt} \,, 

Bu tanım, *L*nin toplam diferansiyeline çok benzerdir. Ancak sanal yer değiştirmeler ve onların zamana göre türevleri sadece diferansiyel terimlerle yer değiştirir (del operatörü ile, belli bir büyüklük olarak ifade edilirler). Ayrıca, sanal yer değiştirme'nin tanımına uyacak şekilde, zaman artışı bu tanımda yer almamaktadır. Onun yerine genelleştirilmiş koordinatların zamanda bir ana göre ifade edilir ve türevlerden ziyade diferansiyel terimlerin çarpımları kullanılır. Zaman göre uygulanan bir parçalı integral yöntemi genelleştirilmiş koordinatların türevlerini yani *δq<sub>j</sub>*ı ∂*L*/∂(*dq<sub>j</sub>*/*dt*) terimine dönüştürür. Diğer bir deyişle *d*(*δq<sub>j</sub>*)/*dt*nın *δq<sub>j</sub>*ya dönüşümü gerçekleştirilir. Buradan, bağımsız sanal yer değiştirmelerin Lagrangian fonksiyonunun türevi ile türetilebileceği ortaya çıkar;

:\int_{t_1}^{t_2} \delta L \, dt = \sum_{j=1}^n\left[\frac{\partial L}{\partial \dot{q}_j}\delta q_j\right]_{t_1}^{t_2} + \int_{t_1}^{t_2} \sum_{j=1}^n \left(\frac{\partial L}{\partial q_j} - \frac{d}{dt}\frac{\partial L}{\partial \dot{q}_j} \right)\delta q_j \, dt  \,. 

Eğer *δq<sub>j</sub>*(*t*1) = *δq<sub>j</sub>*(*t*2) = 0  koşulu her *j* terimi için doğru oluyorsa, integrali alınmayan terimler sıfır olur. Buna ek olarak eğer **δL**in zamana göre integrali sıfır ise, δqj  ‘ler bağımsız olduğu için, integralı sıfıra eşitlemenin tek yolu integrand'ın sıfır olması yani bağımsız genelleştirilmiş koordinat *δq<sub>j</sub>'lerin katsayılarının sıfır olmasıdır. Burada, yukarıdaki figürde bahsedilen durağan etkiler altında, çok küçük değişimler sonucu alınan kırmızı yolun etkileri incelendiği için her durumda konfigurasyon değişimi sıfıra yakın olur (sanal yer değiştirme, diferansiyel yer değiştirme; *δqj' çok küçük). Bu sonuçlardan, sistemin hareket denklemi ortaya konulabilir ve bu, Hamilton Prensibi ile şu şekilde özetlenebilir;

:\int_{t_1}^{t_2}\delta L \, dt = 0 \,. 

Bu duruma göre eğer sistem durağan etkiler altındaysa, Lagrangian'ın zamana göre toplamı sıfır olmaktadır. Eğer sistem üzerinde dinamik etkiler varsa, sistemin yukarıdaki hareket eğrisi kırmızıdan etkiye göre sapmakta ve mavi eğrilere gelmektedir. Bu durumda, Lagrangian'ın zaman göre integralı bize esasen dinamik etkileri verecektir. Dolayısıyla dinamik etkileri tanımlayan bir eylem ifadesini şu şekilde açıklayabiliriz;S = \int_{t_1}^{t_2} L\,dt\,,

*S*, yani eylem; bir fonksiyoneldir (fonksiyonun fonksiyonu). Lagrangian fonksiyonunu *t*1den *t*2ye kadar alır ve skaler bir nicelik döndürür. Açısal momentum, [zaman].[enerji], [momentum].[uzunluk] gibi terimlerin birimiyle aynı birime sahiptir. Bu tanıma göre Hamilton'un minimum eylem prensibi *sanal iş*, sanal-difreansiyel yer değiştirme gibi nicelikler arasındaki ilişki şöyle özetlenebilir;

:\delta S = 0\,. \,\!

Dolayısıyla, Hamilton Prensibi ve Lagrangian yaklaşımı; Newton mekaniğindeki anlayışı tamamen değiştirmiş oluyor. (Metotlar doğru uygulandığında, sonuçların birebir aynı olduğu unutulmamalıdır.) Parçacıklar üzerindeki kuvvetlerin çıktısı olarak ivmeleri düşünmek yerine, modelimizi en az eyleme sahip olan yola ve eylemlerin o yolu nasıl değiştireceğine göre tanımlayabiliriz. Elbette burada eylem ve enerji fonksiyonlarını (Lagrangian) uygun tanımlamak  ve genelleştirilmiş koordinatları veya değişkenleri doğru seçmek gerekmektedir. Parçacıkların konfigurasyon uzayında başlangıç ve bitiş anlık zamanları ile o zamanlara karşılık gelen noktaları tanımlayarak bu yolları ifade edebiliriz. Hamilton'un türettiği ilkeler bazen en az eylem ilkesi ile birebir tutulsa da, en az eylem ilkesinde hareket boyunca uygulanan eylemin durgun olması gerekmektedir. Hamilton prensipleri ile durgun olmayan dinamik eylemleri ve korunumsuz kuvvetleri de sistem modeline dahil etmemiz mümkündür ancak holonomik olmayan kısıtlarda çalışmamız yine mümkün değildir. Holonomik olmayan kısıtlar için ise D'Alambert'in sanal iş kavramından faydalanmamız gerekmektedir.

### Lagrange çarpanları ve kısıtların ifade edilmesi

Lagrangian *L*nin kartezyen **r***k*   koordinatlarında *N* parçacık için değişimini şu şekilde ifade edebiliriz;

:\int_{t_1}^{t_2} \sum_{k=1}^N \left(\frac{\partial L}{\partial \mathbf{r}_k} - \frac{d}{dt}\frac{\partial L}{\partial \dot{\mathbf{r}}_k} \right)\cdot\delta \mathbf{r}_k \, dt = 0 \,. 

Hamilton prensipleri, Lagrangian'daki genelleştirilmiş koordinatlar –sistemin davranışını tamamen tanımlıyor olması koşulu ile- bağımlı olsalar bile geçerlidir. Ancak kısıtların yine de holonomik olması gerekmektedir. Bu açıdan Hamilton prensipleri durum uzayı ifadeleri ile benzerlik gösterir.

Eğer seçilen genelleştirilmiş koordinatlar **r***k*ler bağımsız değil ise sanal işte ve yer değiştirmede yaptığımız gibi δ**r***k*nin katsayılarını sıfıra eşitlememiz mümkün değildir. Çünkü artık bağımsız yer değiştirmelerin sıfır olmayacağından emin olsak bile bağlıların sıfır olma veya onları oluşturacak kuvvetlerin olmaması ya da olması gibi bilinmezlikler ortaya çıkmaktadır. Bu sorunu, Lagrange çarpanları metodunu kısıtları ifade edecek şekilde kullanarak çözebiliriz. Her kısıt denklemi *f<sub>i</sub>*(**r***k*, *t*) = 0yı bir Lagrange çarpanı *λ<sub>i</sub>* for *i* = 1, 2, ...,*C* ile çarpabiliriz. Ve bu ifadeleri orijinal Lagrangian ifadeleri ile birleştirdiğimizde artık yeni bir Lagrangian tanımlamış oluruz;

:L' = L(\mathbf{r}_1,\mathbf{r}_2,\ldots,\dot{\mathbf{r}}_1,\dot{\mathbf{r}}_2,\ldots,t) + \sum_{i=1}^C \lambda_i(t) f_i(\mathbf{r}_k,t) \,. 

Bunu, bağımlı olarak tanımlanmış koordinatları bağımsız koordinatlara çevirmek yerine onları tanımlayan fonksiyonu değiştirmişiz gibi düşünebiliriz. Böylece, özellikle komplike sistemlerdeki bağımsız koordinat seçimi probleminden kurtularak; hesaplamaya daha uygun bir metodoloji ortaya konmuş olur.

Lagrange çarpanları zamanın herhangi birer fonksiyonudur, ancak genelleştirilmiş **r***k* koordinatlarının fonksiyonu değillerdir. Dolayısıyla karmaşık fonksiyonlar oluşturmazlar(Çarpanlar, pozisyon koordinatları ile eşit düzeyde ifade bulmaktadırlar). Bu yeni Lagrangian, sistematik olarak koordinatlarımızın bağımsızlığı ve kısıtların ayrı ifadesi (ki bu iki durum birbirine yakından bağlıdır) sorunlarına çözüm bulabilmektedir. Dolayısıyla, yeni Lagrangian'ı Euler-Lagrange denklemlerine koyarak;

:\frac{\partial L'}{\partial \lambda_i} - \frac{d}{dt}\frac{\partial L'}{\partial \dot{\lambda}_i}  = 0 \quad \Rightarrow \quad f_i(\mathbf{r}_k,t) = 0 

İfadesini kısıt denklemleri olarak elde edebiliriz. Ayrıca **r***k* cinsinden hareket denklemlerini **Lagrange'in  birinci tip denklemleri**ylede şu şekilde ifade edebiliriz;

:\frac{\partial L'}{\partial \mathbf{r}_k} - \frac{d}{dt}\frac{\partial L'}{\partial \dot{\mathbf{r}}_k} = 0 \quad \Rightarrow \quad \frac{\partial L}{\partial \mathbf{r}_k} - \frac{d}{dt}\frac{\partial L}{\partial \dot{\mathbf{r}}_k} + \sum_{i=1}^C \lambda_i \frac{\partial f_i }{\partial \mathbf{r}_k} = 0 \,, 

Bu denklemlerde yeni Lagrangian'ın yukarıda bahsettiğimiz değişimini zaman göre integre edersek;

:\int_{t_1}^{t_2} \sum_{k=1}^N \left(\frac{\partial L}{\partial \mathbf{r}_k } - \frac{d}{dt}\frac{\partial L}{\partial \dot{\mathbf{r}}_k} + \sum_{i=1}^C \lambda_i \frac{\partial f_i}{\partial \mathbf{r}_k } \right)\cdot \delta \mathbf{r}_k \, dt = 0 \,. 

ifadesini buluruz ve temel aldığımız prensiplerden ötürü, uç noktaların ve ona bağlı anların sabit olduğunu her *k* için unutmamamız gerekir. (*δ***r***k*(*t*1) = *δ***r***k*(*t*2) = 0) Çarpanlar, daha önceki metotlarda olduğu gibi *δ***r***k*ların katsayılarını sıfırlayacak şekilde seçilirse, karşılaştığımız çarpan problem sorunu çözülür ve kısıtlar ile bağlılık problemini içsel bir şekilde ifade etmiş oluruz. Böylece bağımsız olmayan genelleştirilmiş koordinatlar için uygun çarpanlarla (kısıtları ifade eden çarpanlar) hareket denklemini bulabiliriz.

Eğer bir korunumlu kuvvetin gradyanı sisteme bir potansiyel enerji *V* bazında etkiyorsa ve bu V sadece **r**kların birer fonksiyonu ise Lagrangian *L* = *T* − *V*nin denklemlere konulması ile;

:\frac{\partial T}{\partial \mathbf{r}_k} - \frac{d}{dt}\frac{\partial T}{\partial \dot{\mathbf{r}}_k} - \frac{\partial V}{\partial \mathbf{r}_k} + \sum_{i=1}^C \lambda_i \frac{\partial f_i}{\partial \mathbf{r}_k} = 0 \,, 

elde edilir. Newton’un 2. Yasasının bu Lagrangian formunda, kinetik enerjinin türevi(negatif olarak denklemde yer bulmakta); *k* parçacığına etkiyen toplam kuvvete; potansiyel enerjilerin türevi de kısıtlayıcı olmayan kuvvetlere karşılık gelmektedir. Yani;

:-\mathbf{F}_k = \frac{\partial T}{\partial \mathbf{r}_k} - \frac{d}{dt}\frac{\partial T}{\partial \dot{\mathbf{r}}_k} \,, \quad \mathbf{N}_k = - \frac{\partial V}{\partial \mathbf{r}_k} \,, 

Ve bunu yukarıdaki temel denkleme koyduğumuzda, kısıtlayıcı kuvvetleri şöyle buluruz;

:\mathbf{C}_k = \sum_{i=1}^C \lambda_i \frac{\partial f_i}{\partial \mathbf{r}_k}  \,, 

ve en temel denklemimize dönersek; birinci terim (pozisyona göre türevlenmiş) momentumun türevini ifade ederken diğer terimler ise kuvvetleri ifade eder ve sistemi çözmemizi sağlar. Burada hem kısıtlayıcı (ve belki sistemimizden boş yere enerji çekecek) kuvvetleri hem de diğer kuvvetleri ayrı olarak bulmamız da önemlidir. Bu sayede sistem optimizasyonu da yapılabilir.

## Lagrange-Euler denklemlerinin özellikleri

Bazı durumlarda; Lagrangian, hareket denklemlerini çözmeye gerek kalmadan bize sistem hakkında bilgi verebilir. Bunu Lagrange’in ikinci tip denklemleri sağlar.

### Çözümlerin tekilliği

Bir sistem için Lagrangian fonksiyonu tekil değildir. Zira, sistemlerin özelliklerine göre de Lagrangian fonksiyonu farklı şekillerde tanımlanmalıdır. Lagrangian *L* sıfır olmayan bir çarpan *a* ile çarpılabilir ve rastgele bir *b* sabiti ile toplanabilir. Dolayısıyla yeni Lagrangian *aL* + *b* olacaktır ve bu fonksiyon, *L* ile aynı hareket tanımını gerçekleştirmektedir. Aynı zamanda Lagrangian, genelleştirilmiş koordinatlardan oluşan fonksiyonların kısmi olmayan, normal türevlerinden de etkilenmez.
 
:L' = L + \frac{df(\mathbf{q},t)}{dt}\,.

Burada *L*&prime; ve *L* tamamen aynı hareketi tanımlarlar.

### Noktasal değişimlerden bağımsızlık

Lagrangian’ın tanımlandığı genelleştirilmiş koordinat setini yani *q*yu başka bir genelleştirilmiş koordinat seti *s*e çevirmek için kullanılacak nokta değişimi **q** = **q**(**s**, *t*) uygulanırsa yeni Lagrangian *L*&prime; (yeni koordinat setinin bir fonksiyonu);

:L(\mathbf{q}(\mathbf{s},t), \dot{\mathbf{q}}(\mathbf{s},\dot{\mathbf{s}},t), t ) = L'(\mathbf{s}, \dot{\mathbf{s}},t) \,,

şeklinde ifade edilir. Ayrıca, zincir kuralı kısmi diferansiyellere uygulanırsa, Lagrange denklemlerinin bu tip bir noktasal dönüşüm altında değişmez olduğu görülür;\frac{d}{dt}\frac{\partial L'}{\partial \dot{s}_i} = \frac{\partial L'}{\partial s_i} \,.

Bu sayede birçok denklem daha basit çalışabileceğimiz koordinat sistemlerine dönüştürülerek; hesaplama süresi azaltılabilir. Bu prosedür, Hamilton mekaniğindeki kanonik dönüşümler ile aynı mantığı içerir. Hamilton
fonksiyonları da benzer şekilde, kanonik dönüşümlerden etkilenmezler.

### Devrî koordinatlar ve korunum yasaları

Lagrangian’ın başka bir önemli özelliği ise, tanımladığı sistem içinde; dinamik davranış boyunca korunan özellikleri ifade etmesidir. Örneğin genelleştirilmiş momentum, *kanonik olarak  genelleştirilmiş* koordinat *q<sub>i</sub>*a eşleniktir ve şu şekilde tanımlanır;
:p_i =\frac{\partial L}{\partial\dot q_i}.

Eğer, Lagrangian *L*  *q<sub>i</sub>*a bağlı değilse Euler-Lagrange denklemleri bize genelleştirilmiş momentumun *korunan* bir değer olduğunu söyler. Çünkü zamana bağlı türevi sıfırdır. Hareket de zamanın değişimiyle gerçekleştiğine göre momentumun hareket boyunca sabit kalması gerektiği sonucuna ulaşırız;

:\dot{p}_i = \frac{d}{dt}\frac{\partial L}{\partial\dot q_i} = \frac{\partial L}{\partial q_i}=0\,.

Bu, Noether teoreminin özel bir halidir. Bu tip koordinatlar *devrî* veya *ihmal* edilebilir koordinatlar olarak tanımlanır. Örneğin sistem   şeklinde bir Lagrangian’a sahip olabilir. *r* ve *z* burada doğru boyuncaki uzunluklar, *s* eğrisel uzunluk, * θ*  ve *φ* ise açılardır. Lagrangian’da *z*,*s* ve *φ*nin bulunmadığına dikkat etmek gerekir. Ancak bu konumların hızları bulunmaktadır. Bu nedenle momentum değerleri;

:p_z =\frac{\partial L}{\partial\dot z}\,,\quad p_s =\frac{\partial L}{\partial \dot s}\,,\quad p_\phi =\frac{\partial L}{\partial \dot \phi}\,,

olarak bulunur ve bunların hepsi korunan değerlerdir. Bu genelleştirilmiş momentumların birimleri tanımlandıkları koordinatlara bağlıdır. Örneğin *p<sub>z</sub>* yönündeki doğrusal hareket momentumu, *p<sub>s</sub>*de *s* boyuncaki doğrusal momentum, *p<sub>φ</sub>* ise ölçüldüğü düzlemdeki açısal momentumu ifade eder. Ancak karmaşık hareketlere sahip olan sistemlerde bu korunumlar hareketi kompleksleştirir.

### Enerji korunumu

Lagrangian’ın zamana göre toplam diferansiyeli alındığında çok genel bir sonuç elde ederiz;

:\frac{\partial L}{\partial t} = \frac{d }{d t}\left(\sum_{i=1}^n \dot{q}_i\frac{\partial L}{\partial \dot{q}_i} - L\right)\,.
 
Eğer bütün Lagrangian zamandan harici olarak bağımsız ise, Lagrangian’ın zamana göre kısmî türevi sıfırdır. () Dolayısıyla parantez içindeki terimi şöyle tanımlayabiliriz;

:E = \sum_{i=1}^n \dot{q}_i\frac{\partial L}{\partial \dot{q}_i} - L 

Bu terim, bütün hareket boyunca sabit kalmalıdır. Buradan da kinetik enerjinin genelleştirilmiş koordinatlara bağlı ikinci dereceden homojen fonksiyon olduğu sonucuna varabiliriz. Ayrıca, ek olarak potansiyel *V* de sadece genelleştirilmiş koordinatların birer fonksiyonu olup hızlardan bağımsız ise, direkt hesaplamayla veya Euler’in homojen fonksiyonlar için geliştirdiği teoremler ile;

:\sum_{i=1}^n \dot{q}_i\frac{\partial L}{\partial \dot{q}_i} = \sum_{i=1}^n \dot{q}_i\frac{\partial T}{\partial \dot{q}_i}  = 2T \,.

sonucuna varılır. Bu durumlar altında, *E* sabiti şöyle tanımlanabilir;

:E = 2T - L = T + V 

Bu ise sistemin *[[Enerjinin korunumu|toplam korunan enerjisine]]* karşılık gerekebilir. Sistem değiştikçe ve dinamik davranış gösterdikçe kinetik ve potansiyel enerjiler değişebilir ancak sistemin hareketi üzerindeki kısıt; toplam enerjisinin sıfır olmasıdır. Bu özellik hesaplamaları oldukça kolaylaştırmaktadır çünkü enerji yani E sanki sabit bir integrasyondaki herhangi bir sabitmiş gibi davranır. Ayrıca bu durumda, enerji ilişkilerinden hızların da integralini almak mümkündür. Ve bu, koordinatlar için yapılacak çözümlere kolaylık sağlar. İki durumda da (kinetik enerji veya hız) zamana bağlı olabilir; bu durumda simetri bozulur ve enerji korunmaz. Burada dikkat edilmesi gereken başka bir konu ise, Hamiltonian’ın Lagrangian ile Legendre dönüşümüne bağlı olan ilişkisinden bu sonuç yine türetilebilir;

: H = \sum_{i=1}^n \dot{q}_i\frac{\partial L}{\partial \dot{q}_i} - L \,. 

Aynı durumlar geçerli ise, Hamilton fonksiyonu sistemin korunan toplam enerjisine eşittir.

### Mekanik benzerlik

Eğer potansiyel enerji, genelleştirilmiş koordinatların zamandan bağımsız homojen fonksiyonu ise ve her pozisyon vektörü sıfır olmayan bir sabit ile büyütülürse, yani;

:V(\alpha\mathbf{r}_1,\alpha\mathbf{r}_2,\ldots, \alpha\mathbf{r}_N)=\alpha^N V(\mathbf{r}_1,\mathbf{r}_2,\ldots, \mathbf{r}_N)

Ve zaman; bir beta faktörü ile büyütülürse, *t*&prime; = *βt*, bu durumda hızlar **v***k*da *α*/*β* faktörü ile büyümüş olur(kinetik enerji de bunun karesi kadar). Tüm Lagrangian aslında aynı faktör ile büyütülür;

:\frac{\alpha^2}{\beta^2}=\alpha^N \quad\Rightarrow\quad \beta = \alpha^{1-N/2}\,.

Uzunluklar ve zaman büyütüldüğü için, sistem parçacıklarının islediği yollar geometrik olarak benzer ama büyüklük olarak farklıdır. *t* zamanda alınmış *l* uzunluğu eğer orijinal yol ise, bu dönüşümlerle artık yeni bir *l&prime;* olacak, zaman ise *t&prime;* olacaktır ve bunlar şu şekilde tanımlanır;

:\frac{t'}{t}= \left(\frac{l'}{l}\right)^{1-N/2}\,.

### İlişkili parçacıklar

Verilen bir sistem için iki alt sistem *A* ve *B* eğer ilişkili değilse sistemin toplam Lagrangian’ı;A* and *LB* for the subsystems:

:<math>L = L_A + L_B \,. </math>

olarak ifade edilir. Ancak parçacıklar ilişkiliyse bu mümkün değildir. Bazı durumlarda sistemin toplam Lagrangian'ını alt parçaların Lagrangian'ları ile ifade etmek çok işe yarar. Bazı daha genel sistemlerde ise ilişkinin olmadığı Lagrangianlar ile ilişki hakkında bilgi taşıyan başka bir Lagrangian *LAB* ile toplayarak sistemi modellemek mümkündür.

:<math>L = L_A + L_B + L_{AB}\,.</math>

Özellikle mekanik ilişkilerin bulunduğu sistemlerde bu doğrusal özellik çok kullanışlıdır. Bu uygulama, fiziksel olarak da etkileşim göstermeyen Lagrangian'ların kinetik enerjiyle ilişkili; etkileşim içinde olanların ise potansiyel enerjiye karşılık geleceği gibi gerçeklenebilir. Önemsenmeyen etkileşimler için Lab sıfır alınabilir. Burada, potansiyel enerjinin daha önceki denklemlerde, sistemlerin koordinatlar boyunca nasıl kısıtlandıkları hakkında bilgi verirken kinetik enerjinin hareket hakkında bilgi verdiği konusu da tekrarlanmış oluyor. Elbette, iki parçacıklı bir sistem için tanımlanan bu yöntem, çok parçacıklı sistemlere genelleştirilebilir. Bu durumda her parçacık arası ilişkilerin ifade edilmesi gerekir.

==Lagrange mekaniğinin başka konulardaki uygulamaları==

Hamiltonian, Lagrangian'a [[Legendre dönüşümü]] uygulanarak bulunur. Dolayısıyla, orijinal değişkenler ile kanonik olarak eşlenik yeni değişkenler elde edilmiş olur. Dolayısıyla Hamilton mekaniğinde değişken sayısı ve diferansiyel denklem sayısı iki katına çıkar ancak dereceler bire indirgenir. Bu, durum uzayı ifadeleri ile de esasen oldukça paraleldir. Hamilton fonksiyonu aynı Lagrangian fonksiyonunun Lagrange denklemleriyle olan ilişkisi gibi [[Hamilton mekaniği]]nin temelidir. Klasik görüşe göre zaman bağımsız bir değişkendir ve hız ile pozisyon zamana bağlıdır. Bu değişkenler ise genelde durum-faz uzayında incelenir ve çevrimler, dengeler ve stabilite bu şekilde incelenir.

[[Analitik Mekanik#Routhian Mekaniği|Routhian mekaniği]]; Lagrange ve Hamilton mekaniklerinin bir hibrit formülasyonudur. Ancak çevrimsel koordinatlar haricinde hesaplama açısından pek etkili değildir.

===Mühendislikteki kullanımları===

Circa 1963'te Lagrangian fonksiyonlarının mühendislikte kullanımını önermiş ancak birkaç yıl sonra, [[dinamik sistem]]lerdeki üstünlüklerine rağmen mühendislik müfredatlarından çıkarılarak teorik dinamik alanına kaydırılmıştır. Zira, Lagrangian sadece makina mühendisliğinde değil, Elektrik Elektronik Mühendisliği; finans, ekonomi, biyoloji ve birçok diferansiyel değişimin olduğu alanda kendine yer bulmuştur. Zira bu alanlarda birçok yol-integral ifadeleri ve parabolik kısmî diferansiyel denklemler bulunmaktadır.

Lagrangian fonksiyonları birçok mühendislik problemine hızlı, yazılıma uygun ve rahat çözümler sunmaktadır. Ayrıca, normal mekanikle formülasyonu bile çok zor olan sistemleri hem felsefi olarak hem de dinamik olarak daha rahat formüle etmemizi sağlar. [[Robotik]], [[akışkanlar mekaniği]]nde akış modelleri, ([[Akışın Euler ve Lagrange spesifikasyonları]]), [[sinyal işleme]], [[nanoteknoloji]], [[süper bilgisayar]]lar, [[kimya mühendisliği]], [[inşaat mühendisliği]] gibi alanlarda kendine yer bulur. Lagrangian gibi enerji fonksiyonları özellikle *akış* kavramının etkin olduğu diferansiyel denklemlerde kendine yer bulur. Özellikle doğrusal olmayan sistemlerde akış konsepti önemlidir ve bunları çözmek için geliştirilen [[Lyapunov fonksiyonları]] da enerji temellidir.

==Ayrıca bakınız==
* [[Lagrange noktası|Lagrange noktaları]]

==Notlar==


==Kaynakça==
*
*
*
*
*
*
*
*
*
*
**The Principle of Least Action*, R. Feynman
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*

== Konuyla ilgili yayınlar ==
* Gupta, Kiran Chandra, *Classical mechanics of particles and rigid bodies'' (Wiley, 1988).
* Cassel, Kevin W.: Variational Methods with Applications in Science and Engineering, Cambridge University Press, 2013.