- 締切済み
微分積分を用いた力学について
物体が空気中で落下運動する時、物体は空気中から抵抗力を受ける。抵抗力の大きさは、物体の速度に比例するとする。すなわち、物体の速度がvのときの抵抗力は-γvと表される。ただしγは正の定数である。 質量mの物体を時刻0において空気中で静かに話し、落下運動させた。鉛直下向きを正の向きにとる。時刻tにおける物体の速度をv、重力加速度をgとすると、物体の運動方程式は m((dv)/(dt))=mg-γv と表されその解はv=Aexp(-λt)+B で与えられる。ただし、A、Bおよびλは、条件によって決まる定数である。次の各問いに答えよ。 (1)十分時間が経過すると、物体は一定の速度に達する。この速度を終端速度と呼ぶ。その速度はいくらか。 (2)初期条件(t=0でv=0)および(1)の結果よりA,Bを定めよ。 (3)解を運動方程式に代入し、任意の時刻tにおいて等式が成立することを考慮すると、λは決まる。λはいくらになるか。
- 物理学
- 回答数3
- ありがとう数2
- みんなの回答 (3)
- 専門家の回答
みんなの回答
- Tann3
- ベストアンサー率51% (708/1381)
No.1です。「お礼」に「よろしくお願いします」とのことですので、今回だけ特別に。 ていねいに書くので、ちょっと長くなります。 No.2さんが、(1)に関しては的確な回答をされていますが、(3)が違うように思います。(指数関数のべき乗係数が「正」では、発散してしまうので) (1)解として、 v = A * exp(-λt) + B (A) が与えられています。 解(A)で、初期条件として t = 0 のとき v = 0 を代入すれば 0 = A + B ∴ A = -B (B) よって、(A)は v = B * [ 1 - exp(-λt) ] (A1) ここで、「t →∞」にすれば、「exp(-λt)→0」ですよね。 ということで、終端速度「v∞」は v∞ = B この「v∞ = B」は、そこでは「速度が一定」になっているということですから、 dv/dt = 0 という条件から、微分方程式で m * g - γ * B = 0 ということです。これから,、γ≠0 として、 B = m * g / γ よって、(A)は v = (m * g / γ) * [ 1 - exp(-λt) ] (A2) これで(2)まで解けてしまいますね。 (注:(1)の模範解答に関しては、No.2さんの方が題意に沿っていて良さそうですね) (2)念のため、解(A1)を運動方程式に代入すると dv/dt = λ * B * exp(-λt) ですから、 m * λ * B * exp(-λt) = m * g - γ * B * [ 1 - exp(-λt) ] 整理すると (m * λ * B - γ * B) * exp(-λt) = m * g - γ * B (C) 「t →∞」にすれば、左辺はゼロになるので、γ≠0 として、 B = m * g / γ (D) (1)の最後の計算とも一致しますね。 (B)より A = - m * g / γ (D)を(A1)に代入して v = (m * g / γ) * [ 1 - exp(-λt) ] (A2) という、(1)と同じ結果が得られます。 (注:これは、t=0 のときに「v = ( λ * m * g / γ) * t 」の直線の接線で始まり、t が増加するとともに徐々にカーブが寝てきて、最終的に「 v = (m * g / γ) 」の一定値に漸近する、という曲線になります) (3)(C)式が任意の時刻 t で成立するためには、左辺で m * λ - γ = 0 であればよいということです。 (注:この仮定がちょっと分かりづらいですが、下の区切り線以下に書いた「微分方程式を解析的に解く」ことと結果を合わせるために、無理やり条件と付けている感じがします) これより、 λ = γ / m 従って、 v = (m * g / γ) * [ 1 - exp(-[γ / m] * t) ] (A3) (解答終わり) **************************** 問題のとおり進めれば上でよいのですが、(3)で「何故こうするの?」というのがちょっと分かりづらいですね。 これは、もともとも微分方程式 m * (dv/dt) = m * g - γ * v (F) をきちんと解析的に解かずに、何となく「解を想定して」解を求めているからです。 質問者さんが、「微分・積分」を勉強しているなら、次のようにきちんと解いてみることをお勧めします。 (F)のままだとちょっと面倒なので、変形します。 dv/dt = g - (γ * v) / m = - (γ / m) * v + g (F2) さらに、積分しやすくするために、γ≠0 として v - m * g / γ = X (G) と置換します。そうすると dX/dt = dv/dt で、(F2)式は dX/dt = - (γ / m) * X (F3) になります。さらに、γ / m = λ とおけば、 (H) dX/dt = - λ * X (F4) この微分方程式は簡単に解けて、 ∫(1/X)dX = ∫(-λ)dt より、C1を積分定数として -ln(X) = -λ * t + C1 従って、 X = C2 * exp(-λ * t) (ただし C2 = exp(C1) (G)式で、Xを元の v に戻して、 v = C2 * exp(-λ * t) + m * g / γ 初期条件 t = 0 のとき v = 0 より C2 = - m * g / γ よって v = (m * g / γ) * [ 1 - exp(-λ * t) ] (H)のλを戻して v = (m * g / γ) * [ 1 - exp(-[γ / m] * t) ] で(A3)と同じ式が得られます。 上の方で、「γ≠0 として」としていますが、では「γ=0」だとどうなるのでしょうか? この場合には、最初の微分方程式は m * dv/dt = m * g つまり dv/dt = g で、簡単に解けて v = g * t + C3 初期条件 t=0 で v=0 から C3=0 で v = g * t という単純な「空気抵抗のない自由落下」になります。 (A3)の曲線も、t=0 のときの「接線の傾き」は「g」になっていますね。 (dv/dtを計算して、t=0 を代入してみてください) ちょっと長くなりました。
(1)は関数の極限値の話ではありません。 物体は一定の速度に達する、ということは加速度が0であるということです。つまり m((dv)/(dt))=mg-γv という運動方程式で、tが十分大きければ dv/dt=0 となる、ということを言っています。これはtが十分大きいところで 0=mg-γv となるということです。つまり v=mg/y が終端速度になります。 (2)を考えるときにはexp(-x)のx→∞の極限値が0であることを知っていなけれはなりません。exp(x)がx→∞で無限大に発散するのですから、exp(-x)は当然0に収束しますね。 このことを知らなかったのなら覚えておかなければなりません。あるいはエクセルでもなんでもいいですから、exp(x)のグラフを書いてこの関数が決まった値に収束しないことを感じておくといいです。 このことを知っていれば終端速度がBと表すことができることがわかります。Aexp(-λt)の項が0に収束するからです。 一方(1)で終端速度をmg/yと求めたのですから B=mg/y です。 また、t=0でv=0なのですから 0=Aexp(-λ・0)+B です。つまり A=-B となります。 A=-mg/y です。 (3)では解を運動方程式に代入してごらんと書いてあるのですから、そうしてみましょう。 m(-λ)Aexp(-λt)=mg-γv Aは上で求めていますから、これも使って書き直すと m(-λ)mg/yexp(-λt)=mg-yv t=0でv=0なのですから m(-λ)mg/y=mg です。"任意の時刻tにおいて等式が成立することを考慮"して、t=0の所で計算すれば十分だと考えました。この式を整理すれば、 λ=-y/m ということになります。 いかが?
- Tann3
- ベストアンサー率51% (708/1381)
それで、どこが分からないのですか? (1)v の「解」がああ得られているのだから、それが t →∞ にしたらどうなるか。 (2)方程式に「解」を代入して整理し、t=0 の初期条件と(1)の条件でA,Bが決まる。 (3)「時間 t を含む項」で、任意の t で成立する条件を考える。具体的には、「時間 t を含む項」の係数がゼロになればよい。 以上のヒントで解けるでしょう? ただの「与えられた「解」の式で微分方程式を解いて、係数の値を確定する」問題ですから。
関連するQ&A
- 微分積分を用いた力学
質量mの物体が一直線上を運動している。この物体には、変位xに比例する復元力F=-kxが働いている。kは正の定数である。時刻t=0のとき物体の変位はx=0、速度はvoであったとする。この物体の運動について次の各問いについて教えてください。 (a)この物体の運動方程式を求めよ (b)(a)の運動方程式の解として時刻tにおける物体の変位xはx=Asinωtと表されるとする。ただしωとAは条件によって決まる定数である。この式を運動方程式に代入し初期条件を用いるとω、Aの値が定まる。ωとAはどんな値になるか。 よろしくお願いします。
- ベストアンサー
- 物理学
- 力学の質問です!
以下の回答をお願いします。 質量mの物体Aが、つりあいのいち(x=0)を中心にx軸上を書く振動数ωで単振動している。 (1)物体Aの一座標x(t)について、運動方程式を示せ。 (2)時刻tにおける物体Aに位置x(t)を一般解で記せ。 (3)この物体に粘性抵抗Fr=-2mrv と時間tに依存する外力Fe(t)=mβt が加わった場合の運動方程式を示せ。 (4)設問(3)でもとめた運動方程式の特解を1つ求めよ。(ヒント:試行関数としてx(t)=at+bを用いて、定数a,bを求める) (5)設問(3)でもとめた運動方程式従う物体Aが南西抵抗Frによって失う単位時間当たりの運動エネルギーを求めよ。
- 締切済み
- 物理学
- 空気抵抗が速度の自乗に比例する場合の運動方程式
自然落下する物体があり、その物体のt秒後の速度は速度が遅いうち(空気抵抗が速度に比例している時)は m・dv/dt=mg-kv mは物体の質量、gは重力加速度、kは空気抵抗定数。 という運動方程式を解いて v=[mg-exp{-k(t+C)/m}]・1/k という一般解を得ます。(Cは積分定数) ここまではできるのですが、この先の空気抵抗が速度の二乗に比例するばあいの運動方程式 m・dv/dt=mg-kv^2 を解くことができないのです。解こうとしてもどうしても、途中で -log|mg-kv^2/C|/2kv=t/m という形になってしまい解くことができません。logの中にvが含まれていて、さらにlogの外にvのある形をv=の形に直せません。強引にv=の形に直せますが、両辺にvが現れてしまい解いたことになりません。 どうすればいいのでしょうか? 余談ではありますが、A=xlogx のx=への変形のしかたや、B=exp(x)+x のx=への変形のしかたを教えてくれると幸いです。
- 締切済み
- 数学・算数
- 力学(抵抗力について)
鉛直上向きにy軸をとり、重力加速度の大きさをgとする。時刻t=0に位置y=y(0)から質量mの物体を鉛直下向きに初速度v(0) (v(0)<0)で放り投げた。物体には、速度vに比例する抵抗力 -bv が働く。 1.運動方程式をvの微分方程式として書き出せ。 2.終端速度の大きさv(t)を求めよ。 3.運動方程式を解いて、t >0におけるv(t).y(t)を求めよ 4.y(0)=0、v(0)=-3v(t)/2 の時、十分時間がたったときのv(t),y(t)の漸近線を求めよ。 問題数多いですが、簡単でいいので途中式もお願いします。
- ベストアンサー
- 物理学
- 慣性抵抗
時刻t=0で十分高い位置にある質量mの物体を速度が0(ベクトル)となるようにそっと放した。この物体には地表付近の重力と物体の速度の大きさの2乗に比例する慣性抵抗が作用する。以下の問に答えよ。 なお、重力加速度の大きさはg、単位質量あたりの慣性抵抗の比例定数はβとし、鉛直方向下向きの単位ベクトルをiベクトル、落下している物体の時刻t≠0での速度をv(t)i(ベクトル)とする。 1)この物体に作用する慣性抵抗Fbを文中の記号であらわせ。 2)この物体に作用する力の総和Fを文中の記号であらわせ。 3)この物体の加速度a(t)を文中の記号であらわせ。 4)この物体の時刻tでの運動方程式(ma(t)=F)を文中の記号であらわせ。 5)4)で求めた運動方程式を微分方程式と考えて解き、この物体の時刻tでの速度を表すv(t)を時間tの関数としてあらわせ。 解く途中経過も示すこと。 6)時刻tが十分経過(t→∞)した場合、この物体の速度を表す関数v(t)はどうなるか、示しなさい。 自分の考えだと 1)はFb=-βv^2i(ベクトル) 2)はF=mgi-mβv^2i(ベクトル) 3)はma(t)=mgi-mβv^2i a(t)=gi-βv^2i までは出来たと思うんですけど、4)~6)がよくわからないので 解き方を教えてください。
- 締切済み
- 物理学
- 雨粒の自由落下の振舞についてです。
雨粒の自由落下の振舞についてです。 ご協力お願いいたします。 変数、定数を適切に定義し、雨粒の自由落下の振舞を微分方程式で近似せよ。また、 2km上空から半径 1mmの球形の雨粒が落下するとき、時刻 tと雨粒の高さ h(t)の関係を求めよ (時刻 tの関数 h(t)を求めよ )。ただし、摩擦 (空気抵抗、抗力 )によって発生する力は、物体の形状に無関係で、物体の速度に比例すると仮定してよい。比例定数は適当な値を選択せよ。 ヒント :微分方程式を解くときは高さ h(t)ではなく、 速度 v(t) := d h(t),速度の微分 dt d v(t)= d2 h(t)を用いて時刻 tと速度 v(t)の関係を先に求めよ。その結果を積分し、積分 dt dt2 定数に適切な値を代入して、関数 h(t)を求めよ。
- ベストアンサー
- 物理学
- (力学)空気抵抗による単振動の減衰運動についての質問
質量mの物体を一端に付けたバネ(バネ定数k自然長k0バネの質量は無視) の他端を天井からつるし、鉛直方向に物体を振動させる。 重力加速度をgとする。その他必要な物理量があれば定義して使ってよい t=0に自然長の位置から物体に初速度を与えずに運動を開始させた。 天井を原点とし、下方を正とするx軸座標を使って以下の問いに答えよ 1.空気抵抗を考える場合(比例定数b)、物体はいわゆる減衰運動を行う。 このとき、この物体の減衰運動の運動方程式をこたえよ。 2.この運動の力学的エネルギーEを求め、ついでその時間変化率をもとめ正負を 判定せよ。この場合の力学的エネルギーとは物体の運動エネルギー、重力の 位置エネルギー、バネの変形エネルギーの和である。 物体に加わる外力が一定の場合は高校の時よく解いた問題なのでわかると思いますが 今回の場合は空気抵抗ということでどのように解いたらよいか良くわかりません。 自由落下時の空気抵抗は解ります。教科書には空気抵抗の記述すらありません。 お願いします
- ベストアンサー
- 物理学
- 抵抗がある場合の自由落下問題を教えてください
空気中を質量mの質点が自由落下する場合に、鉛直下方に働く重力(mg)と、速度に比例した空気抵抗(-cv:cは比例定数)が働くものとする。 1.鉛直下方にx軸の向きをとり、質点の運動方程式を書け。 2.t=0での速さを0とし、運動方程式を解いて、時刻tにおける速度v(t)を求めよ。 3.質点の終端速度 v∞=lim t→∞ v(t)を求めよ。 4.空気抵抗が無くなる極限c→0において、速度はgtになることを示せ。 という問題が分からないです。 合ってるか分からないですが図を書いてみました。図の右側にあるのは3.のlimの表記です。解答例と答えを教えてください。 特に4.が分からなくて困ってます。私が出した1-3の答えは 1. ma = mg - cv 2. v = (mg/c)(1 - exp(-ct/m)) 3. v∞ = mg/c です。
- 締切済み
- 物理学
- 単振動(バネ)
鉛直上向きにy軸をとり、重力加速度の大きさをgとする。ばね定数kのばねの上端を固定し、下端に質量mの物体をつける。ばねが自然長であるときの物体の位置をy=0とする。ばねの質量、空気抵抗は無視できる。物体は鉛直方向のみ運動する。 1、物体の運動方程式を求めよ 2、つりあいの位置y_eを求めよ 3、つりあいの位置からの変位をy_2(t) = y(t) - y_eとする。y_2に関する運動方程式を求めよ 4、運動方程式を解いて、位置y(t)と、v(t) = y '(t)の一般解を求めよ 5、時刻t= 0にy = y_0 の位置で静かに物体を放した(v (0) = 0 )とする。その後の運動y(t),v(t)を求めよ y_eはyの右下に小さいeがあるという意味 よろしくお願いします
- ベストアンサー
- 物理学
お礼
回答ありがとうございます。 現在予習している形で、(1)すら答えを導出できない状態です。 おっしゃられているように(1)の答えがわかれば(2)(3)の答えも導けることは分かるのですが.... よろしくお願いします