Hareketin Katı-Gövde Denklemleri - 2

Daha önce rotasyon matrisi bazlı hareket denklemlerini vermiştik. Bu yazıda cisim duruşu ve dönüş mekaniğini kuaterniyon kavramı üzerinden hesaplayacağız, daha önce görülen bazı kavramların da tekrar üzerinden geçeceğiz.

Lineer Hız, Momentum

Cisim yeri ve hızı şöyle alakalı,

$$ \frac{\mathrm{d} x(t)}{\mathrm{d} t} = v(t) $$

Katı gövdenin momentumu

$$ p(t) = m v(t) $$

Newton'un ikinci kanununa gore,

$$ \frac{\mathrm{d} p(t)}{\mathrm{d} t} = F(t) $$

Yani momentumun zamansal türevi kuvvettir. Hareket hesaplarında $F$ entegre edilerek $p$ hesaplanır, $p$'yi $m$ ile bölerek $v$ elde ederiz, ve onu entegre ederek cisim yeri $x$ elde edilir [1, sf. 466].

Açısal Momentum

Lineer momentum kavramı anlaşılması kolay, onu bir tür atalet olarak görüyoruz, hareket eden bir objenin düz bir çizgi üzerindeki hareketini devam ettirmeye olan meyili. Tabii ki bu meyil, hareketin devamlılığı o anki hız ve cismin kütlesiyle orantılı, $mv$ buradan geliyor. Açısal momentum da benzer bir kavram, tek fark bir eksen etrafında dönmekte olan bir cismin dönmeye devam etme meyili. O zaman üç boyutta, orijin etrafındaki açısal momentum $\vec{L}$, $m$ kütleli cismin orijine olan uzaklığı $\vec{r}$ ile $m\vec{v}$ çapraz çarpımıdır [1, sf. 42].

$$ L = r \times p = r \times mv $$

Aynen kuvvetin momentumun zamansal türevi olması gibi, tork benzer şekilde açısal momentumun türevidir, bunu görmek için üstteki $L$'nin türevini alalım,

$$ \frac{\mathrm{d} L}{\mathrm{d} t} = \frac{\mathrm{d} (r \times p)}{\mathrm{d} t} $$

Eşitliğin sağındaki türev zincirleme kuralı ile şöyle açılır,

$$ = r \times \frac{\mathrm{d} p}{\mathrm{d} t} + \frac{\mathrm{d} r}{\mathrm{d} t} \times p $$

$\mathrm{d} p / \mathrm{d} t = F$, $\mathrm{d} r / \mathrm{d} t = v$ olduğunu biliyoruz, üstte yerine geçirelim,

$$ = r \times F + v \times p $$

Fakat $v$ ile $p$ aynı yöne işaret eden vektörler, onların çapraz çarpımı sıfır, geri kalanlar,

$$ \frac{\mathrm{d} L}{\mathrm{d} t} = r \times F = \tau $$

Bir diğer eşitlik

$$ L(t) = J(t) w(t) $$

Bu formül açısal momentum ile açısal hız $w$ (diğer kaynaklarda $\omega$) ile ilişkilendirir, lineer momentum'daki kütle yerine burada $J$ (diğer kaynaklarda $I$ diye geçer) var, ki $J$ bir atalet tensoru, objenin şekline, ağırlık dağılımına göre değişir, ve bu değerin kendisi de obje döndükçe değişime uğrar.

Hesaplarda gereken uygulanan tork $\tau(t)$'un sebep olduğu cisim dönüşünü bulmak, bunun için $\tau$ entegre edilerek $L$ elde edilir, $L$ değeri $J$ ile bölünerek (daha doğrusu $J$ bir matris olduğu için onun tersi alıp çarpılarak) $w$ elde edilir. Kuaterniyon durumunda

$$ \frac{\mathrm{d} q(t)}{\mathrm{d} t} = \frac{1}{2} \omega(t) q(t) $$

entegre edilerek objenin dönüş sonrası yeni işaret ettiği yer bulunur (burada $\omega$ açısal hızı temsil eden kuaterniyon). Üstteki denklemi türetmek gerekirse [4, sf. 264] kaynağına başvurulabilir. Herhangi iki kuaterniyonun bir geçiş kuaterniyonu ile birbirine bağlanabildiğini biliyoruz. O zaman ufak bir zaman dilimi $\Delta t$ içinde $q(t)$ ve $q(t + \Delta t)$ kuaterniyonlarının nasıl ilişkili olduğuna bakarız, yani aradığımız bir $\Delta r(t)$ öyle ki

$$ q(t + \Delta t) = q(t) \Delta r(t) $$

olsun, ve

$$ \Delta r(t) = \cos (\Delta \alpha) + v(t) \sin(\theta \alpha) $$

Üstteki formülü kuaterniyon tanımından [3] biliyoruz, dönüş açısı $2 \Delta \alpha$, ve etrafında dönüşün yapıldığı eksen vektör $v(t)$.

Dönüş açısı $\Delta \alpha$ çok ufak, bu sebeple küçük açı yaklaşıklaması [5] üzerinden $\cos(\Delta \alpha) \approx 1$, $\sin(\Delta \alpha) \approx \Delta \alpha$, o zaman

$$ \Delta r(t) = 1 + v(t) (\Delta \alpha) $$

Bu demektir ki

$$ q(t + \Delta t) = q(t) [1 + v(t) (\Delta \alpha)] $$

$$ = q(t) + q(t) v(t) (\Delta \alpha) $$

$$ \Rightarrow q(t + \Delta t) - q(t) = q(t) v(t) (\Delta \alpha) $$

Şimdi her iki tarafı $\Delta t$ ile bölüp limit uygularsak,

$$ \frac{\mathrm{d} q}{\mathrm{d} t} = \lim_{\Delta t \to 0} \frac{q(t + \Delta t) - q(t)}{\Delta t} $$

$$ = \lim_{\Delta t \to 0} \frac{q(t) v(t) (\Delta \alpha)}{\Delta t} $$

Dikkat edersek skalar $\Delta \alpha$ değişimi $\Delta t$ ile bölünüyor, bu bize skalar açısal değişim $\omega(t)$'yı verir. Bu açısal hızın vektör $v(t)$ ile çarpılması ise $v(t)$ etrafındaki skalar açısal değişim $\vec{\omega}(t)$'yı verir, yani

$$ = q(t) v(t) \omega(t) = q(t) \vec{\omega}(t) $$

Son formül dönüş açısı $2 \Delta \alpha$ içindi, eğer $\Delta \alpha$ istiyorsak, formülü $1/2$ ile çarparız. Böylece

$$ \frac{\mathrm{d} q(t)}{\mathrm{d} t} = \frac{1}{2} \vec{\omega}(t) q(t) $$

formülüne erişmiş olduk. $q(t)$ için gereken entegrasyon üstteki formül üzerinden yapılır.

Not: Hesapsal hataların birikmemesi için $q(t)$'yi sürekli normalize etmek iyi bir fikirdir, bu normalizasyon

$$ q(t) = \frac{\hat{q(t)}}{|\hat{q(t)}|} $$

ile gerçekleştirilebilir [1, sf. 468].

Hepsini bir araya koyarsak entegrasyonun uygulanacağı nihai yapı şu şekilde gösterilebilir,

$$ \frac{\mathrm{d}}{\mathrm{d} t} \left[\begin{array}{c} x \\ q \\ p \\ L \end{array}\right] = \left[\begin{array}{c} \dot{x} \\ \dot{q} \\ \dot{p} \\ \dot{L} \end{array}\right] = \left[\begin{array}{ccc} m^{-1} p \\ \omega q / 2 \\ F \\ \tau \end{array}\right] $$

Başlangıç Kuvvetleri

Şimdi "sıfırıncı anda" yani ilk başlangıçta uygulanan kuvvetleri, lineer, açısal, hesaplamak lazım. Obje üzerinde uygulanan noktaya bildiğimizi farzedelim, sonu o noktada başlangıcı nesne ağırlık merkezinde olan bir vektör ile kuvvet vektörü arasında çapraz çarpım yapıyoruz, bu bize torku veriyor.

Benzer şekilde sonu nesne merkezinde başı o noktada olan bir vektör daha var, lineer kuvvet bu doğrultuda uygulanacak, o vektör üzerine iki üstte görülen kırmızı vektörü yansıtıyoruz, bu da lineer kuvvet oluyor. Bir üstteki resim üzerinde gösterirsek,

Daha önce söylediğimiz gibi her iki kuvvet de ilk anda lineer ve açısal momentumu ekileyen faktörler, sonraki adımlarda etkileri yok.

Başlangıç torku ve kuvvetin nasıl başlangıç hızına çevirilecebileceğini daha önceki derste gördük.

[devam edecek]

Kaynaklar

[1] Eberly, Game Physics 2nd Ed

[2] Bayramlı, Diferansiyel Denklemler Ders 2

[3] Bayramlı, Döndürme (Rotation) - 2

[4] Kuipers, Quaternions and Rotation Sequences

[5] Bayramlı, Normal Diferansiyel Denklemler, Trigonometri