跳到主要内容

大阪大学 工学研究科 機械工学専攻 2013年度 機械力学・機構学 問題1

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

図のように、質量 mm、長さ ℓ\ell の一様な棒の上端に、左右からばね定数 k/2k/2 の水平ばねを取り付け、各ばねの他端を固定する。棒の上端は摩擦なく水平方向にのみ移動でき、棒は上端のまわりに鉛直面内で自由に回転できる。重力加速度を gg とする。棒の重心を通り紙面に垂直な軸に関する慣性モーメントは mℓ2/12m\ell^2/12 である。

上端のつり合い位置からの水平変位を xx、棒の鉛直下向きからの傾きを θ\theta とする。図のように右向き変位と右側への傾きを正とする。つり合い時の重心を原点とし、右向きに XX、上向きに YY をとる。

  1. 棒の重心の位置 (X,Y)(X,Y) と速度 (X˙,Y˙)(\dot X,\dot Y) を x,θx,\theta を用いて表せ。
  2. 系の運動エネルギーとポテンシャルエネルギーを x,θx,\theta を用いて表せ。
  3. 微小運動を仮定し、ラグランジアンを求めよ。
  4. 微小運動の運動方程式を求めよ。
  5. k/m=g/ℓ=ω02k/m=g/\ell=\omega_0^2 のとき、2 つの固有角振動数 ω1,ω2\omega_1,\omega_2(ω1<ω2\omega_1<\omega_2)と、低い固有角振動数で振動するときの xx と θ\theta の振幅比を求めよ。

原題の図

题目描述​

如图,质量为 mm、长为 ℓ\ell 的均匀杆,上端分别连接左右两根水平弹簧,每根刚度为 k/2k/2,另一端固定。杆上端可无摩擦地沿水平方向移动,杆可绕上端在竖直面内自由转动。重力加速度为 gg,杆绕质心且垂直图面的轴的转动惯量为 mℓ2/12m\ell^2/12。

上端相对平衡位置的水平位移为 xx,杆相对竖直向下方向的倾角为 θ\theta,向右位移及向右倾斜为正。以平衡时杆质心为原点,XX 向右、YY 向上。

  1. 用 x,θx,\theta 表示质心位置和速度。
  2. 求系统的动能和势能。
  3. 在小振动近似下求拉格朗日量。
  4. 求小振动方程。
  5. 当 k/m=g/ℓ=ω02k/m=g/\ell=\omega_0^2 时,求固有角频率 ω1<ω2\omega_1<\omega_2 及低频模态下 xx 与 θ\theta 的振幅比。

Kai​

(1)​

重心は上端から棒に沿って ℓ/2\ell/2 の位置にある。したがって

X=x+ℓ2sin⁡θ,Y=ℓ2(1−cos⁡θ),\boxed{X=x+\frac\ell2\sin\theta},\qquad \boxed{Y=\frac\ell2(1-\cos\theta)},
X˙=x˙+ℓ2cos⁡θ θ˙,Y˙=ℓ2sin⁡θ θ˙.\boxed{\dot X=\dot x+\frac\ell2\cos\theta\,\dot\theta},\qquad \boxed{\dot Y=\frac\ell2\sin\theta\,\dot\theta}.

(2)​

重心の並進と重心まわりの回転を加えると

T=m2(X˙2+Y˙2)+12mℓ212θ˙2=m2x˙2+mℓ2cos⁡θ x˙θ˙+mℓ26θ˙2.\begin{aligned} T&=\frac m2(\dot X^2+\dot Y^2) +\frac12\frac{m\ell^2}{12}\dot\theta^2\\ &=\boxed{\frac m2\dot x^2 +\frac{m\ell}{2}\cos\theta\,\dot x\dot\theta +\frac{m\ell^2}{6}\dot\theta^2}. \end{aligned}

つり合い位置を位置エネルギーの基準に選ぶ。左右のばねのエネルギーを加えれば

U=k2x2+mgℓ2(1−cos⁡θ).\boxed{U=\frac k2x^2+\frac{mg\ell}{2}(1-\cos\theta)}.

(3)​

変位・速度の二次まで残すと

L=m2x˙2+mℓ2x˙θ˙+mℓ26θ˙2−k2x2−mgℓ4θ2.\boxed{\mathcal L= \frac m2\dot x^2+\frac{m\ell}{2}\dot x\dot\theta +\frac{m\ell^2}{6}\dot\theta^2 -\frac k2x^2-\frac{mg\ell}{4}\theta^2}.

(4)​

Euler–Lagrange 方程式より

mx¨+mℓ2θ¨+kx=0,\boxed{m\ddot x+\frac{m\ell}{2}\ddot\theta+kx=0},
mℓ2x¨+mℓ23θ¨+mgℓ2θ=0.\boxed{\frac{m\ell}{2}\ddot x+\frac{m\ell^2}{3}\ddot\theta +\frac{mg\ell}{2}\theta=0}.

(5)​

x=X0cos⁡ωtx=X_0\cos\omega t、θ=Θ0cos⁡ωt\theta=\Theta_0\cos\omega t とおく。固有値条件は

det⁡(k−mω2−mℓω2/2−mℓω2/2mgℓ/2−mℓ2ω2/3)=0.\det\begin{pmatrix} k-m\omega^2&-m\ell\omega^2/2\\ -m\ell\omega^2/2&mg\ell/2-m\ell^2\omega^2/3 \end{pmatrix}=0.

k/m=g/ℓ=ω02k/m=g/\ell=\omega_0^2 を用いて整理すると

ω4−10ω02ω2+6ω04=0.\omega^4-10\omega_0^2\omega^2+6\omega_0^4=0.

よって

ω1=ω05−19,ω2=ω05+19.\boxed{\omega_1=\omega_0\sqrt{5-\sqrt{19}}},\qquad \boxed{\omega_2=\omega_0\sqrt{5+\sqrt{19}}}.

低い方のモードについて第1行から

X0Θ0=ℓω122(ω02−ω12)=1+196ℓ.\frac{X_0}{\Theta_0} =\frac{\ell\omega_1^2}{2(\omega_0^2-\omega_1^2)} =\boxed{\frac{1+\sqrt{19}}{6}\ell}.

比は正なので、xx と θ\theta は同位相で振動する。