大阪大学 工学研究科 機械工学専攻 2018年度 力学 問題2
Author
祭音Myyura (refined with GPT 6 Astra)
Description
液滴が空気中を鉛直に落下しながら、静止している水滴を取り込み、質量を増していく。液滴の質量の時間微分は落下の速さに比例し、比例定数を α > 0 \alpha>0 α > 0 とする。時刻 t t t の質量を m m m 、落下の速さを v v v 、運動量の大きさを p p p とする。初期条件は t = 0 t=0 t = 0 で m = m 0 > 0 , v = 0 m=m_0>0,v=0 m = m 0 > 0 , v = 0 。重力加速度の大きさを g g g とし、空気抵抗は無視する。
運動方程式を p p p の時間微分で表せ。
d p / d m dp/dm d p / d m を p , m p,m p , m を用いて表せ。
落下の速さ v v v を m m m の関数として表せ。
十分に時間が経ち m ≫ m 0 m\gg m_0 m ≫ m 0 となったとする。(3) で m 0 / m → 0 m_0/m\to0 m 0 / m → 0 とすることにより、このときの落下加速度の大きさを求めよ。
题目描述
液滴在空气中竖直下落,吸收原本静止的水滴而增重,质量增长率与下落速率成正比,比例常数为 α > 0 \alpha>0 α > 0 。记质量为 m m m 、下落速率为 v v v 、动量大小为 p p p ,初始时 m = m 0 > 0 , v = 0 m=m_0>0,v=0 m = m 0 > 0 , v = 0 。重力加速度为 g g g ,忽略空气阻力。
用动量的时间导数写出运动方程。
用 p , m p,m p , m 表示 d p / d m dp/dm d p / d m 。
求速率 v v v 关于质量 m m m 的函数。
在长时间后 m ≫ m 0 m\gg m_0 m ≫ m 0 时,利用上一问中 m 0 / m → 0 m_0/m\to0 m 0 / m → 0 的极限求下落加速度。
Kai
(1)
下向きを正とする。取り込まれる水滴の運動量は 0 0 0 なので、液滴の運動量変化は重力の力積に等しい。
d p d t = m g , p = m v , d m d t = α v = α p m . \boxed{\frac{dp}{dt}=mg},\qquad
p=mv,\qquad \frac{dm}{dt}=\alpha v=\frac{\alpha p}{m}. d t d p = m g , p = m v , d t d m = αv = m α p .
(2)
m > m 0 m>m_0 m > m 0 では p > 0 p>0 p > 0 であり、連鎖律より
d p d m = d p / d t d m / d t = g m 2 α p . \boxed{\frac{dp}{dm}=\frac{dp/dt}{dm/dt}
=\frac{gm^2}{\alpha p}}. d m d p = d m / d t d p / d t = α p g m 2 .
(3)
p d p = ( g / α ) m 2 d m p\,dp=(g/\alpha)m^2\,dm p d p = ( g / α ) m 2 d m を初期状態 ( m , p ) = ( m 0 , 0 ) (m,p)=(m_0,0) ( m , p ) = ( m 0 , 0 ) から積分すると
p 2 2 = g 3 α ( m 3 − m 0 3 ) . \frac{p^2}{2}=\frac{g}{3\alpha}(m^3-m_0^3). 2 p 2 = 3 α g ( m 3 − m 0 3 ) .
p = m v p=mv p = m v 、v ≥ 0 v\ge0 v ≥ 0 より
v ( m ) = 2 g 3 α ( m − m 0 3 m 2 ) = 2 g m 3 α [ 1 − ( m 0 m ) 3 ] . \boxed{v(m)=\sqrt{\frac{2g}{3\alpha}
\left(m-\frac{m_0^3}{m^2}\right)}
=\sqrt{\frac{2gm}{3\alpha}
\left[1-\left(\frac{m_0}{m}\right)^3\right]}}. v ( m ) = 3 α 2 g ( m − m 2 m 0 3 ) = 3 α 2 g m [ 1 − ( m m 0 ) 3 ] .
(4)
v 2 v^2 v 2 を m m m で微分し、m ˙ = α v \dot m=\alpha v m ˙ = αv を用いると
v ˙ = α 2 d ( v 2 ) d m = g 3 [ 1 + 2 ( m 0 m ) 3 ] . \dot v=\frac{\alpha}{2}\frac{d(v^2)}{dm}
=\frac g3\left[1+2\left(\frac{m_0}{m}\right)^3\right]. v ˙ = 2 α d m d ( v 2 ) = 3 g [ 1 + 2 ( m m 0 ) 3 ] .
従って
lim m 0 / m → 0 v ˙ = g 3 . \boxed{\lim_{m_0/m\to0}\dot v=\frac g3}. m 0 / m → 0 lim v ˙ = 3 g .