運動方程式を解いてみよう!(粘性抵抗力)
こんにちは!
まろです.
今日は粘性抵抗のある系を考えましょう.
粘性抵抗とは,文字通りネバネバと運動の邪魔をする力のことです.
レイノルズ数(流体力学の用語で,慣性項と粘性項の比を表す無次元数)が小さいときにはおおよそ速度に比例することが分かっています.
今は大きさのない質点を考えているのでレイノルズ数は定義しにくいですが,ここではあまり速くなく,粘っこい流体中の運動と考えればいいでしょう.
つまり,今回考える運動方程式は
です.
は「減衰係数」または「ダンパ係数」と呼ばれる定数です.
「粘性係数」は流体の物性値の名前なので,あまり使わない方がいいと思います.
さて,この式をよく見ると・・・
両辺が簡単に積分できますね!
.
さらにこれは,簡単な微分方程式の形で,
と変形できます.
自分で検算してみてください.
両辺を積分すると,
.
整理して,
となります.
またこれを微分して,
です.
この一般解に初期条件を入れてみましょう.
簡単化のためにここからは方向のみを考えて,
とします.
初期条件を,としてみると,
,
となります.
まとめると,
です.
物理でおなじみの指数関数が出てきましたね.
まぁ線形常微分方程式の解の基底が指数関数だから当然と言えば当然ですが.
指数が負なので,指数の項は時間が経つと減衰していきます.
つまり質点の速度はに収束し,位置はに収束します.
速度が小さくなると抵抗力も小さくなるため,いつまで経っても静止することはないのです.
ばねの時は永遠に同じ振動が続いていましたが,この場合はどんどん質点のエネルギーが小さくなっていきます(「エネルギー」はまた今度定義します).
これが減衰のある運動の特徴です.
次回はもう少し複雑な系を考えてみましょう!
ちなみに,初学者には「青山秀明,力学,学術図書出版社」が読みやすくて良いと思います.
演習問題と解説もちゃんと書いてあります.
著者の大学初級向けの講義の教科書ですので,難しすぎることなくまとめられています.
まろはこの本の2010年度版で入門しました.
ではまた!