- 締切済み
放物線の運動
放物線の運動で空気抵抗を考える場合のプログラムってどう作ったらいいんですか?fortran初心者です。
- fortran
- お礼率0% (0/3)
- その他(プログラミング・開発)
- 回答数3
- ありがとう数1
- みんなの回答 (3)
- 専門家の回答
みんなの回答
- ka_8511
- ベストアンサー率0% (0/0)
1年ほど前に、放物線運動を解析したことがあるので、コメントさせていただきます。 プログラムのことを考える前に、物体の動きを良く理解する必要があります。物体の動きを理解していないと、作ったプログラムの動作が正しいかどうか確認することが出来なくなってしまいます。今回に限らず、技術計算の場合、作ったプログラムが正しく動作しているかどうかは、そのプログラムの基本となる数式通りに計算が行われているかどうかで判断することになります。 さて、物体が空気中で運動するときは、速度の二乗の抵抗を受けます。物体の運動はx軸とy軸の成分に分けて計算しますが、どちらの運動も速度の二乗の抵抗を受けます。空気の抵抗を無視した場合は、x軸方向の力は常に0でしたが、今回は空気の抵抗を考慮していますので、x方向の運動に対しても速度の二乗の抵抗を受ける運動として計算します。 y方向は重力と速度の二乗の力を受けます。重力は常に大地に向かう向きで働きますが、空気抵抗は物体の速度が0に近づくように働きます。 以上が物体の運動を考える上で必要な力です。これらの力を使って物体の運動方程式を作ります。 m(dVx/dt)=-kVx^2 ・・・・・・・・・・(1) m(dVy/dt)=-kVy^2 -mg (上昇時)・・・(2) m(dVy/dt)=kVy^2 -mg (降下時)・・・(3) 但し、kは粘性係数 y方向は上昇と降下でそれぞれ違っの式になってしまうことに注意してください。 この式をオイラー法またはルンゲ・クッタ法で数値的に計算することになります。上昇時は式1と式2を使って計算します。物体が頂点に達したら、即ち、Vy=0になったら式(2)を式(3)代えて計算を続けます。 物体が降下するとき、Vyは一定の速度に達します。これは重力と空気の粘性抵抗の大きさが等しくなるためです。そのときの速度は以下のとおりです。 Vy(終端速度)=sqrt(mg/k) さて、プログラムで上記の計算を行わせて得られた結果が正しいことを証明するには、上記の運動の解析解と比較するのが一番です。 上記の運動の解析解は長い計算を必要としますので割愛します。しかし、計算の一部は、例えば、書籍「力学 新しい視点にたって」戸田ほか共訳 培風館に有ります。 なお、空気抵抗が速度に比例すると考えた場合は、上記の運動方程式が少し違ってきます。(簡単になります)この計算の一部は、例えば、書籍「物理入門コース1 力学」戸田著 岩波書店 に有ります。 放物線運動は解析解が得られるので、それと比較することで、数値計算が正しいかどうか確認できます。私は、表計算でこれらの運動方程式を計算したことがありますが、あまりお勧めしません。但し、このような技術計算が初めての場合は、運動方程式を理解するためには、必要です。しかし、今後も、運動方程式の数値解析を行うのであれば、簡単なシミュレーションソフト(市販)を用意した方が良いと思います。
- Interest
- ベストアンサー率31% (207/659)
Fortran云々の前に、数式で空気抵抗を考慮した物体の運動を記述できますか?(言い換えれば、紙と鉛筆で運動方程式を立てて解けますか?)
- Hayashi_Trek
- ベストアンサー率44% (366/818)
いわゆる微分方程式を解く必要があります。 コンピュータで解くには微小時間ごとの状態を計算することになります。 回転していないボールとして ある瞬間のボールの X座標,Y座標, X方向の速度,Y方向の速度, X方向の空気抵抗,Y方向の空気抵抗,Y方向の引力, X方向の移動量,Y方向の移動量 等を例えば1ミリ秒単位で計算してボールの軌道を求めます。 エクセルでも可能です。
関連するQ&A
- 放物運動の速度の符号について
放物運動なんですが、空気抵抗はなしとして初速V0で打ち上げた時、再び地面に到達する時の速度もV0なのは分かりますが、符号は正でいいのですよね? y軸を上に正と取ってあるので負なのかと思いましたが、それだと運動方向とは逆になりおかしくなりますよね? 考え方が間違っているのか教えてください。
- ベストアンサー
- 物理学
- 鉛直投げ上げ運動&放物運動
真上に物体を投げたとして 最高点の座標Xを導き出したい場合、 鉛直投げ上げ運動だと、 公式は、 x=xo+vot+(1/2)at^2を使うのに対し、 放物運動で最高点の座標を導き出す場合、 公式は、 (v^2)-(vo^2)=2a(x-xo) を使うみたいです。 鉛直と放物で投げ上げ方に違いはありますが、 投げた後どちらも上方向の運動のみを考慮して 物体の速度が0になった状態から計算するのに関わらず、 使う公式が違うのは何故でしょうか?
- ベストアンサー
- 物理学