跳到主要内容

東京農工大学 工学府 機械システム工学専攻 2021年度 専門科目 2

Author​

祭音Myyura (co-authored with GPT 6 Astra)

Description​

2021年4月入学(2020年10月入学を含む)の問題。

水平面に対し α\alpha 傾く斜面に、半径 aa、質量 mm の円柱を置く。中心軸 O は質量を無視できるばね定数 kk のばねで斜面上方の壁に支えられ、円柱は中心軸まわりに滑らかに回転できる。重心 O まわりの慣性モーメントは I=ma2/2I=ma^2/2、斜面との接線上の点を P とする。重力加速度は gg。

  1. P を通る接線まわりの慣性モーメント IpI_p を求めよ。以降、解答に I,IpI,I_p が含まれる場合はそれらの式を代入せず、記号のまま用いること。
  2. つり合い位置を原点に、斜面下向き(左下)を xx の正方向、反時計回りを回転角 θ\theta の正方向に取る。滑らない場合の xx と θ\theta の関係を示せ。
  3. 初期変位を与えて静かに放し、滑らず転がりながら振動するとき、重心位置 xx の運動方程式を示せ。
  4. (3) の固有角振動数 ωn\omega_n を示せ。
  5. 中心軸に、ばねと同じ方向の力 FF を加えるアクチュエータを取り付ける。FF の正方向は xx の正方向とし、アクチュエータは力を加えるほか運動を拘束しない。初期値をゼロとし、FF から xx への伝達関数をラプラス演算子 ss で求めよ。
  6. x=x1x=x_1 を目標として F=−k1(x−x1)−k2x˙F=-k_1(x-x_1)-k_2\dot x で制御する。制御系の減衰比と固有角振動数を求めよ。
  7. ばねを外す。摩擦がなく回転せず滑り落ちる場合の加速度は、滑らず転がり落ちる場合の何倍か。

题目描述​

质量 mm、半径 aa 的圆柱置于倾角 α\alpha 的斜面上,中心轴由沿斜面的弹簧 kk 连至上方固定壁,圆柱可绕中心轴转动,I=ma2/2I=ma^2/2,接触点 P,重力加速度 gg。① 求接触线轴的转动惯量 IpI_p;后续保留符号 I,IpI,I_p。② 平衡处为零,沿斜面向下为 xx 正向、逆时针为 θ\theta 正向,求无滑动关系。③ 求纯滚动振动方程。④ 求固有角频率。⑤ 在中心轴沿斜面加执行器力 FF,求零初值下 X(s)/F(s)X(s)/F(s)。⑥ 采用 F=−k1(x−x1)−k2x˙F=-k_1(x-x_1)-k_2\dot x 时求阻尼比和固有角频率。⑦ 去掉弹簧,比较无摩擦滑动与无滑动滚动的下落加速度之比。

Kai​

(1)・(2)​

平行軸の定理と転がり条件より

Ip=I+ma2,x=aθ.\boxed{I_p=I+ma^2},\qquad\boxed{x=a\theta}.

(3)・(4)​

静的なばね伸びを δ0\delta_0 とすると kδ0=mgsin⁡αk\delta_0=mg\sin\alpha。接点まわりのモーメントから

Ipθ¨=a{mgsin⁡α−k(δ0+x)}=−kax.I_p\ddot\theta=a\{mg\sin\alpha-k(\delta_0+x)\}=-kax.

よって

Ipa2x¨+kx=0,ωn=akIp.\boxed{\frac{I_p}{a^2}\ddot x+kx=0},\qquad \boxed{\omega_n=a\sqrt{\frac{k}{I_p}}}.

(5)​

運動方程式は (Ip/a2)x¨+kx=F(I_p/a^2)\ddot x+kx=F。初期値ゼロでラプラス変換して

X(s)F(s)=a2Ips2+ka2.\boxed{\frac{X(s)}{F(s)}=\frac{a^2}{I_ps^2+ka^2}}.

(6)​

制御則を代入すると

Ipa2x¨+k2x˙+(k+k1)x=k1x1.\frac{I_p}{a^2}\ddot x+k_2\dot x+(k+k_1)x=k_1x_1.

k+k1>0k+k_1>0 として標準的な二次系の形と比較すれば

ωn=ak+k1Ip,ζ=ak22Ip(k+k1).\boxed{\omega_n=a\sqrt{\frac{k+k_1}{I_p}},\qquad \zeta=\frac{ak_2}{2\sqrt{I_p(k+k_1)}}}.

(7)​

滑る場合は x¨slide=gsin⁡α\ddot x_{\mathrm{slide}}=g\sin\alpha、転がる場合は Ipθ¨=mgasin⁡αI_p\ddot\theta=mga\sin\alpha より x¨roll=(ma2/Ip)gsin⁡α\ddot x_{\mathrm{roll}}=(ma^2/I_p)g\sin\alpha。したがって

x¨slidex¨roll=Ipma2=1+Ima2.\boxed{\frac{\ddot x_{\mathrm{slide}}}{\ddot x_{\mathrm{roll}}} =\frac{I_p}{ma^2}=1+\frac{I}{ma^2}}.

Reference​