HOME > natural science Laboratory > コンピュータ・シミュレーション講座 > 計算物理学

【物理シミュレーションに挑戦!】古典力学
斜面を転がる剛体球運動の計算アルゴリズム8:斜面との衝突する剛体球(回転なし)

文責:遠藤 理平 (2016年9月21日) カテゴリ:仮想物理実験室(325)計算物理学(165)

本項は古典力学の様々な系の物理現象を解析的に扱うのではなく、数値計算による物理シミュレーションを実行するために必要な計算アルゴリズムを示すことを目的とします。 様々な初期条件に対する物理シミュレーションを実現するために、最も汎用的な直交座標系を用います。
斜面を転がる剛体球運動の計算アルゴリズム1:摩擦力無しの場合
斜面を転がる剛体球運動の計算アルゴリズム2:摩擦力無しの場合の例(曲面上の運動)
斜面を転がる剛体球運動の計算アルゴリズム3:静止摩擦力のみの場合
斜面を転がる剛体球運動の計算アルゴリズム4:静止摩擦力のみの場合の例(曲面上の運動)
斜面を転がる剛体球運動の計算アルゴリズム5:滑りながら転がる剛体球
斜面を転がる剛体球運動の計算アルゴリズム6:滑りながら転がる剛体球(曲面上の運動)
斜面を転がる剛体球運動の計算アルゴリズム7:空気抵抗力と転がり摩擦抵抗力
斜面を転がる剛体球運動の計算アルゴリズム8:斜面との衝突する剛体球(回転なし)
斜面を転がる剛体球運動の計算アルゴリズム9:斜面との衝突する剛体球(回転あり)
斜面を転がる剛体球運動の計算アルゴリズム10:剛体球同士の衝突(回転なし)
斜面を転がる剛体球運動の計算アルゴリズム11:複数剛体球の同時衝突(回転なし)
斜面を転がる剛体球運動の計算アルゴリズム12:接触中の剛体球に衝突した場合(回転なし)
斜面を転がる剛体球運動の計算アルゴリズム13:剛体球同士の衝突(回転あり)
斜面を転がる剛体球運動の計算アルゴリズム14:固定された剛体球に衝突した場合(回転あり)
斜面を転がる剛体球運動の計算アルゴリズム15:力学的エネルギー保存則を満たした斜面との衝突(回転あり)
斜面を転がる剛体球運動の計算アルゴリズム16:力学的エネルギー保存則を満たした剛体球同士の衝突(回転あり)

前項までは指定した曲面上を運動する剛体球の運動をシミュレーションしました。 本項では斜面と衝突する剛体球をシミュレーションします。まずは剛体球の回転を考慮しない場合の衝突時の計算アルゴリズムを示します。

衝突力の計算アルゴリズム

剛体同士の衝突は一瞬の出来事です。計算時間間隔を \Delta t とし、衝突前と衝突後の運動量ベクトルを \mathbf{ p' }\mathbf{ p } と表した場合、両者の関係は運動量保存則から

\mathbf{ p' } =  \mathbf{ p } + \mathbf{ F }_{\rm coll}\,\Delta t

と表されます。 \mathbf{F}_{\rm coll} が衝突で剛体に加えられた衝突力です。この衝突力は運動エネルギー保存則から得ることができます。 本項では、衝突対象の物体が完全に固定されている系(衝突による運動なし)における衝突力は剛体の速度ベクトルを \mathbf {v}、 質量をM、衝突面の法線ベクトルを \mathbf {n} と表した場合、

\mathbf{F}_{\rm coll}  = \frac{ -M(1+e)(\mathbf {v} \cdot \mathbf {n} )}{ \Delta t}\, \mathbf{n}

となります(詳細な導出はこちらの書籍を参照ください)。eは反発係数と呼ばれ、通常 0\leq e \leq1 となります。e=1で弾性衝突(エネルギー減衰なし)、e=0で衝突後のなし非弾性衝突(エネルギー減衰のある場合)を表現することができます。なお、反発係数eは、衝突前後の速度ベクトルと衝突面の法線ベクトルを用いて

|\mathbf{v'}\cdot \mathbf{n}| = e|\mathbf{v}\cdot \mathbf{n}|

で定義されます。

重心運動に対する運動方程式

衝突力を加えた運動方程式は次のとおりです。衝突の瞬間のみに適用されます。

M\,\frac{d^2\mathbf{R}}{dt^2} =M\mathbf{g}+\mathbf{f} _{\rm coll}+ \mathbf{f}_{\gamma}

ただし、\mathbf{f}_{\gamma}空気抵抗力です。

斜面と衝突する剛体球のシミュレーション

次の物理シミュレーションは、若干傾いた斜面に衝突する剛体球です。衝突時に剛体球の回転は考慮されません。 斜面との接触時には従来の斜面を滑りながら転がります。

曲面と衝突する剛体球のシミュレーション

次の物理シミュレーションは、曲面に衝突する剛体球です。衝突時に剛体球の回転は考慮されません。 斜面との接触時には従来の斜面を滑りながら転がります。



▲このページのトップNPO法人 natural science トップ

関連記事

仮想物理実験室







計算物理学

▲このページのトップNPO法人 natural science トップ




Warning: mysqli_connect(): (28000/1045): Access denied for user 'xsvx1015071_ri'@'sv102.xserver.jp' (using password: YES) in /home/xsvx1015071/include/natural-science/include_counter-d.php on line 8
MySQL DBとの接続に失敗しました