N88-BASICで剛体力学 (1回目)

2024/1/28(日)
N88-BASICで剛体力学 (1回目)
 
剛体力学
(Rigid body mechanics)
 
斜面の転がりと滑り(移動)
 
■ 定義
https://ulprojectmail.blogspot.com/2024/01/rigid-2.html
剛体力学 (2回目)
より
 
▼ 座標
θ:回転角(rad)
ω:回転角速度(rad/s)
α:回転角加速度(rad/s2)
φ:傾斜角(rad)
m :剛体の質量(kg)
r :球(円柱)の半径(m)
h :円柱の高さ(m)
g :重力加速度(m/s2)
a :剛体回転軸位置の加速度(m/s2)
v :剛体回転軸位置の速度(m/s)
x :剛体回転軸位置(m)
I :慣性モーメント(kg･m2)
F :静止摩擦力(N)
 
I = (2/5)mr2  … 球
I = (1/2)mr2  … 円柱
 
▼ 運動方程式
ma = mgsinφ – F
a = (mgsinφ)/(m + I/r2)
 
I = 0      if F = 0
a = gsinφ if F = 0
 
a = (5/7)gsinφ  … 球
a = (2/3)gsinφ  … 円柱
 
a = const.(一定)
v = v0 + at
x = x0 + v0t + (1/2)at2 
 
■ 動作
斜面の傾斜角と初めの位置の高さを入力し
高さ0m地点までの移動距離、経過時間
加速度
高さ0m地点での速度を
 
円柱、球、滑る(摩擦なし)
について
それぞれ表示します
 
VL,NL,XL-BASICとdlg～.zip(rigid001.bas)は
以下のリンクからダウンロードできます
https://ulprojectmail.blogspot.com
Readme.txtを読んで遊んで下さい
















 

このブログの人気の投稿

NEWS

2024/12/27(金 ) N88-BASICでハノイの塔 (2回目)(再帰なし) を公開しました    2024/12/25(水) VL,NL,XL-BASIC(ver～28w3) の?LBasic01.txt修正版公開しました    2024/12/24(火) VL,NL,XL-BASIC(ver～28w3)  を公開しました    2024/12/21(土 ) ハノイの塔 (3回目)  を公開しました
続きを読む

N88-BASICでゲーム (1回目)

イメージ
2021/8/16(月) N88-BASICでゲーム  (1回目) by ULproject   球同士の衝突と跳返り   大文字ベクトル、小文字スカラー 文字の後の、 'は衝突後、数字は球番号   球同士の反発係数、衝突面に水平 ,垂直u,e 球の中心位置 P,半径r,質量m,速度V 衝突面の球 1方向の法線N = (P 1 -P 2 )/|P 1 -P 2 |   球同士の衝突は、 |P 1 -P 2 | ≦ r 1 +r 2   なので √{(P 1 -P 2 )･(P 1 -P 2 )} ≦ r 1 +r 2   (P 1 -P 2 )･(P 1 -P 2 ) ≦ (r 1 +r 2 )(r 1 +r 2 ) で判断できます。   ここで、ベクトルの計算ですが、 P(i) (i=0,1,2) をPのx,y,z成分と すると、 P = P 1  - P 2 は、 P(i) = P1(i) - P2(i) (i=0,1,2) P･Pは ΣP(i)P(i) (i=0,1,2) となります   BASICでは、 P(0) = P1(0) - P2(0) P(1) = P1(1) - P2(1) P(2) = P1(2) - P2(2) D = P(0)*P(0)+P(1)*P(1)+P(2)*P(2) R = R1 + R2 IF D <= R*R THEN 重なり又は接触 となります。   反発させるときは、球が接触している 間は再び反発しないようにしないと 球が離れなくなります。   次に反発についてですが、衝突面に 水平な成分は、球の回転に関係します が、ここでは球の回転は考えずに 垂直成分と同じように反発係数で 減速させる事にします。   衝突面に垂直な成分 (v 1 は球 1の速度、v 1 'は衝突後) 運動量 p = mv 運動量保存則 p 1 +p 2 =p 1 '+ p 2 'より p 1 +p 2  = m 1 v 1 '+ m 2 v 2 ' 反発係数 e = -(v 1 '- v 2 ')/(v 1  - v 2 )より ev 1  - ev 2  = -v 1 '+ v ...
続きを読む