****正方形枠の回転2 鉛直面内に立てて,重力を考慮する。ひとまず,姿勢保持に必要なトルクの導出を考えてみたい。 ---- 重力を考慮した系のラグランジアンは, $$L = 2ma^2\left\{\frac{5}{3}({\dot\theta}^2+{\dot\phi}^2) + 2\dot\theta\dot\phi\cos(\phi - \theta)\right\} + \tau_\theta\theta + \tau_\phi\phi - 4mga(\sin\theta + \sin\phi)$$ となる。ラグランジュ方程式をつくり,2階微分について解いて整理すると $$\ddot\theta = \frac{3\sin(\phi - \theta)\{5{\dot\phi}^2 + 3\cos(\phi - \theta)\cdot{\dot\theta}^2\}+\displaystyle\frac{3}{4ma^2}\{5\tau_\theta - 3\cos(\phi - \theta)\tau_\phi\} - \displaystyle\frac{3g}{a}\{5\cos\theta - 3\cos\phi\cos(\phi - \theta)\}}{25 - 9\cos^2(\phi - \theta)}$$ $$\ddot\phi = \frac{-3\sin(\phi - \theta)\{5{\dot\theta}^2 + 3\cos(\phi - \theta)\cdot{\dot\phi}^2\}-\displaystyle\frac{3}{4ma^2}\{3\cos(\phi - \theta)\tau_\theta - 5\tau_\phi\}+\displaystyle\frac{3g}{a}\{3\cos\theta\cos(\phi - \theta) - 5\cos\phi\}}{25 - 9\cos^2(\phi - \theta)}$$ を得る。 まず手始めに,角 $$\theta , \phi$$ への姿勢保持において,左下の固定回転軸まわりに必要なトルクを求めてみる。$$\ddot\theta=\ddot\phi=\dot\theta=\dot\phi=0$$ とおいて, $$5\tau_\phi - 3\cos(\phi - \theta)\cdot\tau_\phi = 4mga\{5\cos\theta - 3\cos\phi\cos(\phi - \theta)\}$$ $$3\cos(\phi - \theta)\cdot\tau_\theta - 5\tau_\phi = 4mga\{3\cos\theta\cos(\phi - \theta) - 5\cos\phi\}$$ となる。これを $$\tau_\theta , \tau_\phi$$ について解けば, $$\tau_\theta = 4mga\cos\theta\qquad \tau_\phi = 4mga\cos\phi$$ を得る。 Algodoo の設定は,$$m=1.0{\rm kg} , \theta=\pi/6 , \phi=\pi/3$$ である。Algodooのモーターの強制停止時トルクでは静止させられなかったので,静止させるための支点を与えて力を見ることにした。 #ref(http://www14.atwiki.jp/yokkun?cmd=upload&act=open&pageid=319&file=Square3.bmp) こうして見ると,姿勢保持に要するトルクを求めるのにラグランジアンから導くのは全く大げさで,初歩的な剛体のつりあいを考えれば十分である。ラグランジアンの強みは,運動にこそある。しかし,ラグランジュの方法が,束縛における「[[仮想仕事の原理]]」をその手法として内に含んでいるということをあらためて学んだ例となった。 ---- Algodoo シーン >http://www14.atwiki.jp/yokkun?cmd=upload&act=open&pageid=319&file=Square3.phz ----