忍者ブログ
20080511~ 13と7と11の倍数の論理積は13と7と11の積の倍数である。 和ァ・・・
[779] [4508] [4507] [4506] [4505] [4504] [4503] [4502] [4501] [4500]

前回、粘性抵抗と慣性抵抗両方を受ける落下運動について運動方程式を解いた。

以下のような結果になったが
v=(4D^2-B^2)/{4Dcoth(-Dt/G)+2B} 


G=m/c2、B=c1/c2、A=mg/c2

このようにB,A,G,Dを定義したことを踏まえて、元の格好に戻すと

v=2mg/{√(c1^2+4mgc2)coth{t√(c1^2+4mgc2)/(2m)}+c1} 

このような姿になる。

ここで、c1やc2をそれぞれゼロとしたとき、

元の微分方程式の解と一致しているかどうかを確かめてみよう。

c1≠0、c2=0の場合は

v=2mg/{√(4mgc2)coth{t√(4mgc2)/(2m)}}=2mg/[c1{1+coth(c1t/(2m))}] 

となり、

運動方程式

mv'=mg-c1v 

の解と一致しているはずだ。

早速解いてみると

-(m/c1)v'=v-(mg/c1) 

dv{v-(mg/c1)}=-(c1/m)dt 

ln{v-(mg/c1)}=-(c1/m)t+G2 

v-(mg/c1)=G2exp(-c1t/m) 

初期条件としてt=0でv=0を与えると任意の定数G2=-mg/c1となるので

v-(mg/c1)=-(mg/c1)exp(-c1t/m) 

v=(mg/c1)-(mg/c1)exp(-c1t/m) 

v=(mg/c1){1-exp(-c1t/m)} 

一方、v=2mg/[c1{1+coth(c1t/(2m))}] こちらはcothの中身を展開し

大きい分母の分母分子にexp(+)-exp(-)をかけると


となって見事一致する。

ちなみに、終端速度もちゃんと一致していることを確かめておこう。

t→∞でのvはv→mg/c1に収束するが

これは実は、元の運動方程式

mv'=mg-c1vのv'=0の条件で解けば、微分方程式ではなくただの代数方程式として
0=mg-c1v
c1v=mg
v=mg/c1

と、簡単に解けてしまう。

拍手[0回]

PR

コメント


コメントフォーム
お名前
タイトル
文字色
メールアドレス
URL
コメント
パスワード
  Vodafone絵文字 i-mode絵文字 Ezweb絵文字


忍者ブログ [PR]
カレンダー
03 2024/04 05
S M T W T F S
1 2 3 4 5 6
7 8 9 10 11 12 13
14 15 16 17 18 19 20
21 22 23 24 25 26 27
28 29 30
ブログランキング
ブログランキング参戦中
にほんブログ村 アニメブログ 深夜アニメへ
にほんブログ村 漫画ブログ SF・ファンタジー漫画へ
にほんブログ村 科学ブログ 自然科学へ
よかったらポチッとお願いします^^
最新CM
[12/30 buy steroids credit card]
[09/26 Rositawok]
[03/24 hydraTep]
[03/18 Thomaniveigo]
[03/17 Robertaverm]
最新TB
プロフィール
HN:
量子きのこ
年齢:
43
性別:
男性
誕生日:
1981/04/04
職業:
WinDOS.N臣T
趣味:
妄想・計算・測定・アニメ
自己紹介:
日記タイトルの頭についてるアルファベットは日記の番号です
26進数を右から読みます
例:H→7番目、XP→15(P)×26+23(X)=413番目。
A=0とする仕様につき一番右の桁はAにできませんのでご了承くださいズコー
バーコード
ブログ内検索
アクセス解析