N88-BASICで終端速度 (1回目)

2024/6/16(金)
N88-BASICで終端速度 (1回目)
 
(Terminal velocity)
 
速度に比例する空気抵抗のある投げ上げ
 
■ 前提
▼ 参照
https://ulprojectmail.blogspot.com/2024/08/terminal-1.html
終端速度 (1回目)
より
 
▼ 定義
g:重力加速度(m/s²)
v₀:初速度(m/s)
v:速度(m/s)
k:空気抵抗の比例定数(N･s/m) … 1m/s毎の力(N)
m:質量(kg)
F:力(N)
a:加速度(m/s²)
y:位置(m)
y₀:初期位置(m)
v_∞:終端速度(m/s)
 
▼ 定数
v_∞ = -mg/k, λ = k/m
 
▼ 加速度
a(t) = -(v₀ - v_∞)λexp(-λt)
a(0) = -(v₀ - v_∞)λ = -(k/m)v₀ - g
a(t) = 0  (t→∞)
 
▼ 速度
v(t) = v₀exp(-λt) + v_∞{1 - exp(-λt)}
v(0) = v₀ 
v(t) = v_∞  (t→∞)
 
▼ 位置
y(t) = y₀ + (v₀ - v_∞)(1/λ){1 - exp(-λt)} + v_∞t
y(0) = y₀ 
y(t) = y₀ + (v₀ - v_∞)(1/λ) + v_∞t  (t→∞)
 
■ 解説
上向きを正として
空気抵抗-kvがある場合とない場合の
投げ上げの
加速度a、速度v、位置yの時間t変化を
描画しました
 
 
VL,NL,XL-BASICとdlg～.zip(term001.bas)は
このブログ(以下のリンク)から
ダウンロードできます

https://ulprojectmail.blogspot.com

Readme.txtを読んで遊んで下さい