忍者ブログ
20080511～ 13と7と11の倍数の論理積は13と7と11の積の倍数である。 和ァ・・・
[4509] [4508] [4507] [4506] [4505] [4504] [4503] [4502] [4501] [4500] [4499]
粘性抵抗、慣性抵抗、落下運動、両方、運動方程式
質量をm、重力加速度をg、速度v、時間t
速度の1乗に比例する空気抵抗の係数をc1、速度の2乗に比例する空気抵抗の係数をc2とすると

運動方程式は以下のように書ける。
速度は高さy方向の下向きをプラスに定義する。

mv'=mg-c1v-c2v^2 

ここで、v^2の係数を1にしたいので、両辺を-c2で割ると

-(m/c2)v'=v^2+(c1/c2)v-(mg/c2) 

ここで、式の簡単化のためG=m/c2、B=c1/c2、A=mg/c2とすると

-Gv'=v^2+Bv-A 

dvの分母に右辺を持ってくると以下のようになる。

dv/(v^2+Bv-A)=-dt/G 

両辺を積分してvの式にすると、速度vについての運動が解けるのだが、vの2乗が邪魔でうまく積分できないので、部分分数分解を用いる。


1/(v^2+Bv-A)=A1/(v+B/2-D)+A2/(v^2+B/2+D) 

という恒等式を作り、A1、A2、Dがいくつになるのか定めることにする。

まず、(v+B/2-D)(v+B/2+D)=v^2+Bv-A になるようにDを定めたい。

(v+B/2)^2-D^2=v^2+Bv+B^2/4-D^2=v^2+Bv-A 

なので、D^2=B^2/4+A であることがわかる。なぜいきなりDという文字を仮定したのかはあとでわかるが、これは2次方程式の判別式Dの意味である。


次にA1とA2を定める。
A1(v+B/2+D)+A2(v+B/2-D)=1 という恒等式からA1とA2を定めたいので

v(A1+A2)=0
からA1=-A2が導かれるので
B/2(A1+A2)もゼロである。

D(A1-A2)=1にA1=-A2を大入して、

A1=-A2=1/(2D) を得る。


そうすると元の微分方程式は

1/(2D){dv/(v+B/2-D)-dv/(v+B/2-D)}=-dt/G 
なので積分ができ

1/(2D){ln(v+B/2-D)-ln(v+B/2+D)}=-t/G+G2 

となる。G2は任意の積分定数。

対数同士の引き算は対数の中の割り算として引っ込めることができるため

1/(2D)ln{(v+B/2-D)/(v+B/2+D)}=-t/G+G2 

両辺に2Dを掛け算して
ln{(v+B/2-D)/(v+B/2+D)}=-2Dt/G+G2 
(ここでG2はまだ任意の定数なので、2DG2を新たなG2として上書きする)


両辺の指数をとって

(v+B/2-D)/(v+B/2+D)=G2exp(-2Dt/G) 
(ここでもexp(G2)を新たなG2として上書きしている)

左辺の分母分子に2を掛け算してやると
(2v+B-2D)/(2v+B+2D)=G2exp(-2Dt/G) 


初期条件
t=0でv=0を与えて任意の定数G2を定める。

(B-2D)/(B+2D)=G2 

微分方程式の解は
(2v+B-2D)/(2v+B+2D)=(B-2D)/(B+2D)×exp(-2Dt/G) 

となる。


これをvの式に直す。
両辺に(B+2D)と(2v+B+2D)を掛け算する。

(2v+B-2D)(B+2D)=(2v+B+2D)(B-2D)×exp(-2Dt/G) 

展開すると
2v(B+2D)+(B-2D)(B+2D)=2vexp(-2Dt/G)(B-2D)+(B-2D)(B+2D)exp(-2Dt/G)

全部左辺に移行し、vを含むか含まないかで整理すると

2v(B+2D)-2vexp(-2Dt/G)(B-2D)+(B-2D)(B+2D)-(B-2D)(B+2D)exp(-2Dt/G)=0

2vB{1-exp(-2Dt/G)}+4vD{1+exp(-2Dt/G)}+(B-2D)(B+2D){1-exp(-2Dt/G)}=0


2vB{1-exp(-2Dt/G)}+4vD{1+exp(-2Dt/G)}=-(B^2-4D^2){1-exp(-2Dt/G)}  

2v[B{1-exp(-2Dt/G)}+2D{1+exp(-2Dt/G)}]=(4D^2-B^2){1-exp(-2Dt/G)} 

ここで1/(2v)の式に着目してみることにする。
1/(2v)=[[B{1-exp(-2Dt/G)}+2D{1+exp(-2Dt/G)}]]/[(4D^2-B^2){1-exp(-2Dt/G)}]



第1項は約分できるので
1/(2v)=[B/(4D^2-B^2)+2D{1+exp(-2Dt/G)}]/[(4D^2-B^2){1-exp(-2Dt/G)}] 

問題は第2項だが

{1+exp(-2x)}/{1-exp(-2x)}の分母分子にexp(x)を掛け算すると

{exp(x)+exp(-x)}/{exp(x)-exp(-x)}となって、これは双曲線関数ハイパボリックタンジェントtanhの逆数ハイパボリックコタンジェントcothであることがわかる。
{exp(x)+exp(-x)}/{exp(x)-exp(-x)}=coth(x) 

1/(2v)=B/(4D^2-B^2)+2D/(4D^2-B^2)×coth(-Dt/G)



ここでようやく2vの式に戻すと
v=(4D^2-B^2)/{4Dcoth(-Dt/G)+2B}

拍手[0回]

PR
【2023/05/27 11:39 】
CATEGORY[物理] COMMENT[0]

コメント


コメントフォーム
お名前
タイトル
文字色
メールアドレス
URL
コメント
パスワード
  Vodafone絵文字 i-mode絵文字 Ezweb絵文字
<<粘性抵抗、慣性抵抗、落下運動、両方、運動方程式2 HOME トンネル障壁にパルス波をぶち込む>>


忍者ブログ [PR]
カレンダー
01 2025/02 03
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
リンク

「本日午後からにわか雨」
ブログ村
プロヒ
ツイ太
管理画面
新しい記事を書く
検索エンジン登録 SEM UG
人気ブログランキング
bloog
カテゴリー
未選択 ( 2 )
雑念 ( 1101 )
アニメ・漫画 ( 1406 )
節約 ( 19 )
哲学？ ( 7 )
科学？ ( 505 )
日常？ ( 74 )
ニュス ( 38 )
夢 ( 34 )
オカルト ( 77 )
ストーリ ( 43 )
今日の、定理。 ( 44 )
照明節電 ( 4 )
マクロでやらないマクローリン展開シリーズ ( 5 )
量子力学～井戸型ポテンシャル～ ( 31 )
量子力学～水素原子～ ( 31 )
DX全波整流シリーズ ( 5 )
次元のバーゲンセール ( 22 )
ピタゴラ・ジュース・メーカー ( 8 )
対角化シリーズ ( 53 )
シュタインズゲート ( 90 )
キロクログ ( 6 )
けいおん ( 33 )
シャンゼリオン ( 12 )
数学 ( 359 )
物理 ( 131 )
相対論 ( 32 )
電気 ( 35 )
地球防衛企業ダイガード ( 0 )
複素行列をエクセルで ( 1 )
パウリ行列指数関数 ( 49 )
量子力学 ( 71 )
ブログランキング
ブログランキング参戦中
にほんブログ村 アニメブログ 深夜アニメへ
にほんブログ村 漫画ブログ SF・ファンタジー漫画へ
にほんブログ村 科学ブログ 自然科学へ
よかったらポチッとお願いします＾＾
最新CM
無題
[12/30 buy steroids credit card]
мега darknet market
[09/26 Rositawok]
Новая ссылка на гидру v3 hydra onion
[03/24 hydraTep]
Wisdоm is а mixturе оf experiеnсе, cоuragе аnd intеlligenсe
[03/18 Thomaniveigo]
ван вин
[03/17 Robertaverm]
最新記事
はじめまして
(01/01)
4次の特性方程式の実装
(01/03)
SU(4)の特性方程式の実数解の解析
(09/23)
特殊ユニタリSU(4)と4次方程式の4つの実数解
(09/23)
調和振動子の昇降演算子
(02/11)
粘性抵抗、慣性抵抗、落下運動、両方、運動方程式4
(05/30)
粘性抵抗、慣性抵抗、落下運動、両方、運動方程式3
(05/28)
粘性抵抗、慣性抵抗、落下運動、両方、運動方程式2
(05/28)
粘性抵抗、慣性抵抗、落下運動、両方、運動方程式
(05/27)
トンネル障壁にパルス波をぶち込む
(08/04)
Excelでマイクラ風フルル
(10/24)
2次元井戸型ポテンシャル。解析解につきこのような形のポテンシャルに
(06/08)
4次元空間＋1次元時間状態図の熱力学「カルノーサイクル」
(05/22)
3Dカルノーサイクル状態図(Excelおもちゃ)
(01/13)
身勝手の極意～井戸型ポテンシャル編～
(11/04)
最新TB
プロフィール
ＨＮ：
量子きのこ
年齢：
43
ＨＰ：
希望的観測日記　ブログ枝部
性別：
男性
誕生日：
1981/04/04
職業：
WinDOS.N臣T
趣味：
妄想・計算・測定・アニメ
自己紹介：
日記タイトルの頭についてるアルファベットは日記の番号です
26進数を右から読みます
例：H→7番目、XP→15(P)×26＋23(X)=413番目。
A=0とする仕様につき一番右の桁はAにできませんのでご了承くださいｽﾞｺｰ
バーコード
RSS
RSS 0.91
RSS 1.0
RSS 2.0
ブログ内検索
アーカイブ
2038 年 01 月 ( 1 )
2025 年 01 月 ( 1 )
2024 年 09 月 ( 2 )
2024 年 02 月 ( 1 )
2023 年 05 月 ( 4 )
2022 年 08 月 ( 1 )
2021 年 10 月 ( 1 )
2021 年 06 月 ( 1 )
2021 年 05 月 ( 1 )
2021 年 01 月 ( 1 )
2020 年 11 月 ( 2 )
2020 年 10 月 ( 7 )
2020 年 06 月 ( 8 )
2020 年 05 月 ( 20 )
2020 年 03 月 ( 1 )
最古記事
Bこんなデスノートは嫌だ。てんこもり編(BはBファクトリー)
(05/11)
C：今、お笑いが寒い！(Cは炭素言語)
(05/11)
D：いまさらながら大長編ドラえもんの考察(Dはデューテリウム・ドライブ)
(05/13)
E関西黒メディアの娯楽以外の実用性(Eはエネルギー空っぽ)
(05/13)
F：ダイGシーン、主に見ない永続CM、花より団子(Fは昔のレーサーアニメ)
(05/13)
G：1ヶ月アイスですごすときに残り日数を片手で数える方法ほか(Gは重力定数と地上での加速度)
(05/13)
H：ハンドメイドラッキー大王仮面(Hはエッチなハンドメイド・ハミルトーニャン)
(05/13)
I：土産にもらった地元のう・どん(Iは行列の虚数単位)
(05/13)
J：乾燥(Jは色んなものの流れ密度。虚数単位)
(05/14)
K：仮面ノリダー072(Kは仮面ノリダーの運動エネルギー)
(05/14)
L：中国自身(Lは昔の僕と死んだ男)
(05/14)
M：上戸で愛撫、青き白木(Mは相互誘導定数、宇宙の理論、頭の形、女王様)
(05/14)
N：車載プレーヤーついにバグってきた最終局面(Nは北極、でもそこはS。中性(子))
(05/16)
O：ケータイメモの予想外な文字節約復活(Oは酸素リング)
(05/16)
P：デカ・ルートの最終結論について(Pは運動量)
(05/16)
アクセス解析
コガネモチ