HOME > natural science Laboratory > コンピュータ・シミュレーション講座 > 仮想物理実験室

VisualC++ と OpenGL を利用した仮想物理実験室
【2-5-2】重力による斜方投射運動２のシミュレーション

文責：遠藤 理平 （2010年8月29日） カテゴリ：仮想物理実験室(325)

【2-5-1】重力による斜方投射運動における初速度の分解にて導入した初期値を用いて、斜方投射運動のシミュレーションを行います。初速度 v_0[m/s]を固定したまま、投射角度θを変更します。 θ= 0°～ 90°まで、10°刻みで変えて投射するシミュレーションを行ないます。 当然、角度によって v_{x0}と v_{z0}とが変わり運動の様子が変化しますが、 それぞれ z軸方向の等加速度直線運動と、x軸方向の等速直線運動であることは変わりません。

VisualC++ ＋ OpenGL プログラミング

【９日目】重力による運動：斜方投射運動の初期値を変更することで、シミュレーションを行ないます。

仮想物理実験室変数の定義

//--------------------------------------------------------
// 仮想物理実験室変数の定義
//--------------------------------------------------------
double t = 0.0;     //時刻
double dt= 0.03;    //時間刻み
int tn = 0;         //ステップ数
double g = 9.80665; //重力加速度 g
double degree;      //角度：度（°）
double radian;      //角度：ラジアン(rad)
double v0=40.0;     //初速度

// ボールの定義
struct BALL {
  double x, y, z;
  double vx, vy, vz;
  double ax, ay, az;
};
const int N = 10;
BALL ball[N];
void SetUp(void){
  for(int i=0; i<N; i++ ){
    degree = 10.0 * double(i);        //度数法で表現
    radian = 2.0 * PI * degree/360.0; //弧度法への変換
    ball[i].x  = -50.0;
    ball[i].y  = 0.0;
    ball[i].z  = 0.0;
    ball[i].vx = v0 * sin(radian);
    ball[i].vy = 0.0;
    ball[i].vz = v0 * cos(radian);
    ball[i].ay = 0.0;
    ball[i].az = -g;
  }
}

ボールを１０個用意して、それぞれの投射角度を１０°ずつ変更します。

計算と物体の描画

void Calculate(){
  t = dt * double(tn);
  for(int i=0; i<N; i++){
    //速度の算出
    ball[i].vx = ball[i].vx + ball[i].ax * dt;
    ball[i].vy = ball[i].vy + ball[i].ay * dt;
    ball[i].vz = ball[i].vz + ball[i].az * dt;
    //位置の算出
    ball[i].x = ball[i].x + ball[i].vx * dt;
    ball[i].y = ball[i].y + ball[i].vy * dt;
    ball[i].z = ball[i].z + ball[i].vz * dt;
  }
 tn++;
}
void DrawStructure(){
  for(int i=0; i<N; i++){
    glPushMatrix();
      glMaterialfv(GL_FRONT, GL_AMBIENT, ms_ruby.ambient);
      glMaterialfv(GL_FRONT, GL_DIFFUSE, ms_ruby.diffuse);
      glMaterialfv(GL_FRONT, GL_SPECULAR, ms_ruby.specular);
      glMaterialfv(GL_FRONT, GL_SHININESS, &ms_ruby.shininess);
      glTranslated(ball[i].x, ball[i].y, ball[i].z); //平行移動値の設定
      glutSolidSphere(4.0, 20, 20);            //引数：(半径, Z軸まわりの分割数, Z軸に沿った分割数)
    glPopMatrix();
  }
}

i 番目のボールの 位置「ball[i].x」「ball[i].y」「ball[i].z」
速度「ball[i].vx」「ball[i].vy」「ball[i].vz」
加速度「ball1.ax」「ball[i].ay」「ball[i].az」
を逐次計算します。水平方向には等速直線運動、鉛直方向（z軸方向）には、重力による等加速度直線運動を行ないます。初期値として、v_{z0}と v_{x0} を設定することでシミュレーションすることができます。

20100127-01.gif

VisualC++ と OpenGL を利用した仮想物理実験室

第０章　仮想物理実験室の構築

第１章　様々な運動

第２章　ニュートンの運動方程式

第３章　剛体の運動（エネルギー保存則と運動量保存則）

付録

  • 【Ａ－１】参考文献
    ・（A-1-1）OpenGL について
    ・（A-1-2）VisualC++ について
    ・（A-1-3）物理シミュレーション
    ・（A-1-4）数値計算

未分類

力学

量子力学

波動論



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

関連記事

仮想物理実験室

【第８話】有限の高さの障壁へ照射アニメーション【Pythonコピペで量子力学完全攻略マニュアル】 2021.09.29 【遠藤 理平｜仮想物理実験室】

【第７話】無限に高い障壁に向けた電子パルスの照射アニメーション【Pythonコピペで量子力学完全攻略マニュアル】 2021.09.24 【遠藤 理平｜仮想物理実験室】


【第６話】無限に高い障壁へ照射アニメーション【Pythonコピペで量子力学完全攻略マニュアル】 2021.09.23 【遠藤 理平｜仮想物理実験室】

【第５話】電子パルスの運動アニメーション【Pythonコピペで量子力学完全攻略マニュアル】 2021.09.21 【遠藤 理平｜仮想物理実験室】


【第４話】電子パルスの作り方【Pythonコピペで量子力学完全攻略マニュアル】 2021.09.20 【遠藤 理平｜仮想物理実験室】

【第３話】自由粒子の運動アニメーション【Pythonコピペで量子力学完全攻略マニュアル】 2021.09.19 【遠藤 理平｜仮想物理実験室】


【第２話】自由粒子のスナップショット【Pythonコピペで量子力学完全攻略マニュアル】 2021.09.18 【遠藤 理平｜仮想物理実験室】

【第１話】プログラムの動作確認【Pythonコピペで量子力学完全攻略マニュアル】 2021.09.17 【遠藤 理平｜仮想物理実験室】


１０日で作る量子ドットコンピュータ～動作原理を数値シミュレーション！～【５日目】１量子ビット万能量子ゲートを設計しよう！ 2019.09.06 【遠藤 理平｜仮想物理実験室】

１０日で作る量子ドットコンピュータ～動作原理を数値シミュレーション！～【４日目】量子井戸の電子状態の基礎実験３：２重量子井戸で量子ビットを改良！ 2019.08.27 【遠藤 理平｜仮想物理実験室】



▲このページのトップ ｜ 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との接続に失敗しました