科学のおもちゃ箱 @wiki

ペットボトルの口を水平にはじいて，再度立たせるゲーム，「ヒットペット」。
ペットボトルの起き上がり
および
Interactive Physics
で紹介した，ペットボトルの運動のモデルシミュレーションに挑戦していたが，概ね満足できる結果を得たので紹介する。

本来，おもりとして水をわずかに入れたペットボトルが対象物体であるが，簡単のため次のようにモデル化する。図のように質量およびの質点が長さの軽い棒でつながれていて，鉛直に立っている状態から上の質点に水平な初速度を与えるものとする。

座標として，重心の座標および鉛直からの角度をとる。下の質点は床から離れないものとすれば，運動の自由度は2であるから，これらの3つの座標は独立ではない。さて，それぞれの座標変数に対する運動方程式および拘束条件は，次のようになるだろう。

　　　　　　　　　　---(1)
　　　　　　　---(2)
　　---(3)
　　　　　　　　　　　　---(4)

ここに，は物体と床の動摩擦係数，は未知の垂直抗力である。なお物体と床の間は常にすべりの状態にあるものとし，静止摩擦は考慮しない。また，慣性モーメントは，

である。
(4)を時間で微分すると，


となるから，(2)によりを消去した上で(3)に代入して整理すると，に関する二階微分方程式

を得る。これをまず，初期条件

のもとで積分することになる。見るからに解析的な解は望むべくもない。
床との衝突は，はねかえり定数3/4の摩擦のない衝突とし，衝突後の位置と速度を初期条件にして同様に数値積分すれば，起き上がりまでの運動の過程を得るだろう。
手持ちの数学ソフトMathcadによって数値積分した結果を下図に示す。

赤丸が重心座標の積分結果である。ソフトに描かせたグラフの上に物体の位置と向きを
描き込んでみた。衝突までの軌跡はほとんどそのまま倒れる場合に同じで，衝突後はさっと足を引き寄せる感じの動きがおもしろい。次の図は，エネルギーを追跡してみたグラフである。はじめ持っていた運動エネルギーが衝突と摩擦によって散逸する力学的エネルギーに等しければ，直立して止まることになる。最後鉛直に立ったとき運動エネルギーがゼロにならなかったのは，ご愛嬌？^^;