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