跳到主要内容

九州大学 工学府 機械系専攻 2018年度 機械力学 問3

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

長さ ll、張力 TT、線密度 ρ\rho の弦の両端が壁に固定されている。左端を原点、長手方向を xx とし、横変位を u(x,t)u(x,t) とする。弦は空気中で単位長さ当たり ν(∂u/∂t)\nu(\partial u/\partial t) の抵抗力を受け、張力は一定である。c=T/ρc=\sqrt{T/\rho}、μ=ν/ρ\mu=\nu/\rho を用いてよい。(50 点)

両端固定の弦(公式原図)

  1. 運動方程式を導出せよ。
  2. 自由振動解を u(x,t)=U(x)W(t)u(x,t)=U(x)W(t) と仮定し、U,WU,W の微分方程式を求めよ。
  3. UU の一般解を求めよ。
  4. 振動するための ν\nu または μ\mu の条件と、そのときの WW の一般解を求めよ。
  5. 境界条件から rr 次の固有モード UrU_r を求めよ。
  6. 自由振動の一般解と rr 次の減衰固有角振動数 ωdr\omega_{dr} を求めよ。

题目描述​

两端固定弦长 ll、张力 TT、线密度 ρ\rho,单位长度受到与横向速度相反的黏性阻力 νut\nu u_t。推导阻尼波动方程,以乘积分离变量,求空间与时间方程、振荡条件、固有模态、一般自由响应及各模态阻尼频率。

Kai​

(1)​

微小要素の力のつり合いから

ρutt dx=Tuxx dx−νut dx,utt+μut=c2uxx.\rho u_{tt}\,dx=T u_{xx}\,dx-\nu u_t\,dx, \qquad\boxed{u_{tt}+\mu u_t=c^2u_{xx}.}

(2)​

分離定数を −ω2-\omega^2 とすると

W¨+μW˙W=c2U′′U=−ω2.\frac{\ddot W+\mu\dot W}{W}=c^2\frac{U''}{U}=-\omega^2.

よって

U′′+ω2c2U=0,W¨+μW˙+ω2W=0.\boxed{U''+\frac{\omega^2}{c^2}U=0,\qquad \ddot W+\mu\dot W+\omega^2W=0.}

(3)​

U(x)=Acos⁡(ωx/c)+Bsin⁡(ωx/c).\boxed{U(x)=A\cos(\omega x/c)+B\sin(\omega x/c).}

(4)​

時間方程式の特性根は λ=−μ/2±μ2/4−ω2\lambda=-\mu/2\pm\sqrt{\mu^2/4-\omega^2}。したがって振動条件は

μ<2ω(ν<2ρω),\boxed{\mu<2\omega\quad(\nu<2\rho\omega),}

その場合

W(t)=e−μt/2[Ccos⁡(ω2−μ2/4 t)+Dsin⁡(ω2−μ2/4 t)].\boxed{W(t)=e^{-\mu t/2} \left[C\cos\left(\sqrt{\omega^2-\mu^2/4}\,t\right) +D\sin\left(\sqrt{\omega^2-\mu^2/4}\,t\right)\right].}

(5)​

U(0)=U(l)=0U(0)=U(l)=0 と非自明性から A=0A=0、ωrl/c=rπ\omega_rl/c=r\pi(r=1,2,…r=1,2,\ldots)。よって

Ur(x)=sin⁡rπxl,ωr=rπcl.\boxed{U_r(x)=\sin\frac{r\pi x}{l},\qquad \omega_r=\frac{r\pi c}{l}.}

r=0r=0 は恒等的にゼロとなるため固有モードではない。

(6)​

各モードが不足減衰である場合、

u(x,t)=∑r=1∞e−μt/2(Arcos⁡ωdrt+Brsin⁡ωdrt)sin⁡rπxl,\boxed{u(x,t)=\sum_{r=1}^{\infty} e^{-\mu t/2}(A_r\cos\omega_{dr}t+B_r\sin\omega_{dr}t) \sin\frac{r\pi x}{l},}
ωdr=(rπcl)2−μ24.\boxed{\omega_{dr}=\sqrt{\left(\frac{r\pi c}{l}\right)^2-\frac{\mu^2}{4}}.}

全モードを不足減衰とする十分必要条件は μ<2πc/l\mu<2\pi c/l。より大きい減衰では、μ=2ωr\mu=2\omega_r のモードの時間因子を (Ar+Brt)e−μt/2(A_r+B_rt)e^{-\mu t/2}、μ>2ωr\mu>2\omega_r のモードでは Areλr+t+Breλr−tA_re^{\lambda_{r+}t}+B_re^{\lambda_{r-}t} に置き換える。定数は初期変位・初速度で定まる。