三体同時衝突はどう解いたらいいのでしょうか?
撃力ベースで質量・大きさが同じ正方形の弾性衝突の2Dシミュレーションをしていて、
物体 X の速度を vX、物体 X に対する撃力を fX と書くとして
物体 A と B で衝突が起こったら
fA += vB-vA, fB += vA-vB
という撃力が発生するように実装しました。
ここで、止まっている中央の物体 C に、左右から同じ速度 v で他の2体 L, R が衝突する
状況をシミュレーションすると、
fL = +v
fC = +v-v = 0
fR = -v
となり、撃力適用後の速度は、
vL = 0
vC = 0
vR = 0
となって固まってしまいます。
vL = +v
vC = 0
vR = -v
というのが解だと思うのですが、これを得るにはどう取り扱ったらいいのでしょうか?
そもそも3体が同時に衝突することは確率的にあり得ず、2体同士の3回の衝突に分けて計算する、
というのが一般的な打開案ですが、今回は位置・速度の演算を有理数形式で管理し、
厳密なシミュレーションをするとどうなるか、ということを念頭を置いているので、
衝突時刻の誤差を期待できません。
cf.
検索すると、三体衝突問題は数学的に解けないことが証明されている、と書いてあるページがちらほらみつかりました。
関係があるかもしれません。
(注:3天体の万有引力相互作用がカオス系になる、
というのはまた別で、ここではあくまで三体の同時衝突の問題です。
根は同じかも知れませんが)
よろしくお願いします。
A 回答 (2件)
- 最新から表示
- 回答順に表示
No.2
- 回答日時:
お使いの撃力の式は2体衝突には有効でも、三体衝突には適さないものですね。
3体衝突では、L / R には C だけでなく、R / L からの撃力も間接的に(圧縮波を通じて)作用しますが、これが取り入れられていません。そのために L と R が止まってしまうわけです。2体衝突と3体衝突の両方に使えるような撃力の式が見つかればよいですが、見つからなければ、2体の場合と3体の場合を分けて扱う必要があるでしょう。その際、跳ね返り係数を使う方法を検討されたらどうでしょう。
例えば、跳ね返り係数が 1 である場合、質量が同じ2物体の衝突では速度が交換されます。
3物体の場合には、衝突後の速度を「’」を付けて表すと、
vL' - vC' = vC - VL
vR' - vC' = vC - vR
vL' + vC' + vR' = vL + vC + vR (運動量保存)
より
vL' = (-vL + 2vC + 2vR) / 3
vC' = (2vL - vC + 2vR) / 3
vR' = (2vL + 2vC - vR) / 3
が得られます。例えば
vL = v
vC = 0
vR = -v
に対しては
vL' = -v
vC' = 0
vR' = v
で、望みの結果となります。
ありがとうございます。
vL' - vC' = vC - VL (作用反作用の法則)
vR' - vC' = vC - vR (作用反作用の法則)
vL' + vC' + vR' = vL + vC + vR (運動量保存)
これで目が覚めました。
二体衝突も元はこの式から導出しているのだから、
変数が増えるだけで式も増えてるから同じようにすればいいのですね。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 「反発係数=e」「衝突前の球の速度=V」「平面な床に衝突後の球の速度をv」とした時、eを表す関係式は 7 2022/12/18 18:06
- 物理学 物理の問題(車関係)で質問があります。 東方向に15.5m/sで走行中の車(1250kg)が角度不明 6 2022/12/09 13:17
- 物理学 高1力学の運動量の問題です。問題を一通り解いたのですが、行き詰まってしまったのでご回答頂ければ嬉しい 3 2022/06/29 11:20
- 物理学 衝突後の相対速度=-e×(衝突前の相対速度)という式がありますが、物体AとBがあり衝突前のAの速度を 2 2023/05/25 16:51
- 物理学 台車の上の物体が衝突時に受ける力積について 4 2023/01/19 13:58
- 物理学 質量 M,半径αの円板が1つの直径を固定軸として回転できるようになっている。質量mの物体が速さvで円 2 2022/10/21 20:16
- 物理学 長さaの軽い棒の各端に質量mの物体A,Bを取り付け、なめらかな床の上におき、これを棒の中点Oを中心と 2 2022/10/09 19:16
- 物理学 力学の衝突の事例です。物体Aが物体Bに衝突する時、物体Bが運動する方向は衝突時に受ける撃力の方向に制 1 2022/12/24 15:04
- 物理学 高校物理 二次元の衝突 画像の問題の解答では、静止系での球2の速度v2を -運動エネルギー保存 -運 3 2022/11/12 00:34
- 物理学 力学の問題です。質量m1、速度v1の物体Aと質量m2、速度v2の物体Bがx軸上を等速直線運動していて 2 2022/12/24 13:26
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
等加速度運動の変位の公式でな...
-
【至急】衝突時間ってどうやっ...
-
力学・v-tグラフの面積がなぜ移...
-
物理の教科書で最初の公式 s = ...
-
慣性の法則 走行中に人間が飛び...
-
あるくのがはやい
-
速度は「速くなる」?「大きく...
-
速さは速度の絶対値?
-
物理の問題で速さを求めよとい...
-
100mで10秒と11秒の距離の差は?
-
メートル毎秒m/sについてです
-
重力加速度
-
制動距離が速度の2乗に比例す...
-
12kmを普通の自転車で行ったら
-
Allegretto moderato とは??
-
ドリフト移動度とホール移動度...
-
台車を押す力 (物理?)
-
「速く移動すると時間の進み方...
-
「陽」という言葉の意味が分か...
-
電圧の単位V(ボルト)の表記につ...
おすすめ情報