跳到主要内容

東京大学 工学系研究科 機械工学専攻 2021年度 機械力学・制御 問題7

Author​

祭音Myyura (Based on zzy's answer refined with GPT 6 Astra)

Description​

水平面内で微小回転する一様な棒状剛体 A を、固定軸受 B と2本のばね S1、S2 で支える。棒の全長は LL、質量は mm、重心回りの慣性モーメントは II。ばね取付位置は初め重心から左右に距離 aa とする。

  1. 軸受が重心にあり、ばね定数が k1,k2k_1,k_2、回転粘性減衰係数が cc のとき、回転角 θ\theta の運動方程式を求めよ。
  2. 平衡位置から変位させて静かに放すと、振動せずに元へ戻った。このような2種類の応答を図で説明せよ。

次に軸受だけを重心から右へ bb 移し、軸受の回転抵抗を0、ばね定数を k1=k2=kk_1=k_2=k とする。

  1. 運動方程式を求めよ。
  2. b=0b=0 の場合と比べた固有振動数の変化を説明せよ。
  3. 固有振動数を b=0b=0 のときと等しくするため、S2 の取付位置だけを変える方法と、S2 のばね定数だけを変える方法をそれぞれ求めよ。

重心支持と偏心支持の棒・ばね系

题目描述​

均匀刚杆在水平面内绕固定轴承微小转动,两弹簧起初安装在重心两侧距离 aa 处。先求重心支承时带粘性阻尼的运动方程,并比较过阻尼与临界阻尼自由响应。然后将轴承向右移 bb、忽略阻尼,求方程和固有频率变化;分别仅调右侧弹簧位置或刚度,使频率恢复到轴承位于重心时的值。

Kai​

(1)​

ばね変位はそれぞれ ±aθ\pm a\theta なので、軸回りのモーメントの釣合いより

Iθ¨+cθ˙+(k1+k2)a2θ=0.\boxed{I\ddot\theta+c\dot\theta+(k_1+k_2)a^2\theta=0}.

(2)​

ωn=(k1+k2)a2I,ζ=c2I(k1+k2)a2.\omega_n=\sqrt{\frac{(k_1+k_2)a^2}{I}},\qquad \zeta=\frac{c}{2\sqrt{I(k_1+k_2)a^2}}.

初期条件を θ(0)=θ0, θ˙(0)=0\theta(0)=\theta_0,\ \dot\theta(0)=0 とする。

臨界減衰 ζ=1\zeta=1 では

θ=θ0(1+ωnt)e−ωnt.\boxed{\theta=\theta_0(1+\omega_nt)e^{-\omega_nt}}.

過減衰 ζ>1\zeta>1 では、λ1,2=−ωn(ζ∓ζ2−1)\lambda_{1,2}=-\omega_n(\zeta\mp\sqrt{\zeta^2-1}) として

θ=θ0−λ2eλ1t+λ1eλ2tλ1−λ2.\boxed{\theta=\theta_0\frac{-\lambda_2e^{\lambda_1t}+\lambda_1e^{\lambda_2t}}{\lambda_1-\lambda_2}}.

両者とも振動しないが、同じ慣性・剛性で静かに放したとき、臨界減衰の方が速く平衡位置へ戻る。

臨界減衰と過減衰の自由応答

(3)​

軸回りの慣性モーメントは I+mb2I+mb^2、ばねの腕の長さは a+b,a−ba+b,a-b。よって

(I+mb2)θ¨+2k(a2+b2)θ=0.\boxed{(I+mb^2)\ddot\theta+2k(a^2+b^2)\theta=0}.

(4)​

固有角振動数と固有振動数は

ωn(b)=2k(a2+b2)I+mb2,fn(b)=ωn(b)2π.\boxed{\omega_n(b)=\sqrt{\frac{2k(a^2+b^2)}{I+mb^2}},\qquad f_n(b)=\frac{\omega_n(b)}{2\pi}}.

b=0b=0 との差は

ωn2(b)−ωn2(0)=2kb2(I−ma2)I(I+mb2).\omega_n^2(b)-\omega_n^2(0) =\frac{2kb^2(I-ma^2)}{I(I+mb^2)}.

したがって b>0b>0 では、a<I/ma<\sqrt{I/m} なら増加、a>I/ma>\sqrt{I/m} なら減少、a=I/ma=\sqrt{I/m} なら不変である。

(5)​

目標値は ωn2(0)=2ka2/I\omega_n^2(0)=2ka^2/I とする。

S2 の軸受からの距離を dd に変えると

k[(a+b)2+d2]I+mb2=2ka2I.\frac{k[(a+b)^2+d^2]}{I+mb^2}=\frac{2ka^2}{I}.

よって

d2=D:=2a2(1+mb2I)−(a+b)2.\boxed{d^2=D:=2a^2\left(1+\frac{mb^2}{I}\right)-(a+b)^2}.

D≥0D\ge0 かつ棒上に取付可能な場合、重心基準の位置を x2=b±Dx_2=b\pm\sqrt D のいずれかにする。

位置を変えずに S2 のばね定数を k2′k_2' とすると

k(a+b)2+k2′(a−b)2I+mb2=2ka2I\frac{k(a+b)^2+k_2'(a-b)^2}{I+mb^2}=\frac{2ka^2}{I}

より、a≠ba\ne b に対して

k2′=kD(a−b)2.\boxed{k_2'=k\frac{D}{(a-b)^2}}.

通常の正のばね定数で実現するには D>0D>0 が必要である。a=ba=b では S2 の腕が0なので、そのばね定数を変えても固有振動数は変わらない。

Reference​