名古屋大学 工学研究科 機械航空系 2017年8月実施 基礎部門 問題5
Author
祭音Myyura (co-authored with GPT 6 Astra)
Description
問題5
図1のような長さ l l l の質量が無視できる棒がある。棒の一端は支持点 O に取り付けられ,棒はそのまわりで自由に回転できる。棒の他端は質量 m m m の物体の重心 G に取り付けられており,その重心 G まわりの物体の慣性モーメントは I I I である。この棒を回転ばね定数(*)が k k k のばねで支持する。以下,これを振り子と呼ぶ。振り子の角度は鉛直上方向から反時計まわりを正とする角度 θ \theta θ で表す。ばねが自然長(変形していない状態)のときの振り子の角度を θ 0 \theta_0 θ 0 とする。鉛直下方向に重力が作用するものとし,重力加速度を g g g とする。以下の問いに答えよ。
(*)回転ばね定数:加えたモーメントを振り子の角変位 θ − θ 0 \theta-\theta_0 θ − θ 0 で割った比例定数
(1) 支持点 O まわりの振り子の慣性モーメント I A I_A I A を求めよ。以下では,慣性モーメントの記号 I A I_A I A をそのまま用いよ。
(2) 振り子の運動方程式を求めよ。
(3) 振り子の平衡状態(つりあいの状態)を θ = θ e q \theta=\theta_{\mathrm{eq}} θ = θ eq とする。平衡状態 θ = θ e q \theta=\theta_{\mathrm{eq}} θ = θ eq を求める条件式を示せ。
(4) ばねの自然長の状態の角度を θ 0 = 0 \theta_0=0 θ 0 = 0 とする。
1) そのときの平衡状態を示す関係式を求めよ。
2) θ e q = 0 \theta_{\mathrm{eq}}=0 θ eq = 0 はその平衡状態の1つである。いま,振り子をそのごく近傍 θ ≈ 0 \theta\approx0 θ ≈ 0 から角速度 θ ˙ = 0 \dot\theta=0 θ ˙ = 0 で離す。平衡状態 θ e q = 0 \theta_{\mathrm{eq}}=0 θ eq = 0 まわりの微小範囲の運動の特性を回転ばね定数 k k k の値に対して分類し,それぞれ説明せよ。
(5) 平衡状態が θ e q = π / 2 \theta_{\mathrm{eq}}=\pi/2 θ eq = π /2 であるとする。
1) そのときのばねの自然長の状態の角度 θ 0 \theta_0 θ 0 を求めよ。
2) その平衡状態のごく近傍 θ ≈ π / 2 \theta\approx\pi/2 θ ≈ π /2 で振り子を離すと,その平衡状態 θ e q = π / 2 \theta_{\mathrm{eq}}=\pi/2 θ eq = π /2 まわりで振り子が微小振幅で単振動する。運動方程式を用いて,その固有角振動数を求めよ。
题目描述
问题5
如图1,有一根长度为 l l l 、质量可忽略的杆。一端安装于支点 O,可绕该点自由转动;另一端连接质量为 m m m 的物体的质心 G。物体绕质心 G 的转动惯量为 I I I 。杆由转动弹簧常数(*)为 k k k 的弹簧支承,下文称此系统为摆。摆角 θ \theta θ 从竖直向上方向起算,以逆时针为正。设弹簧处于自然状态(未变形)时的摆角为 θ 0 \theta_0 θ 0 。重力竖直向下,重力加速度为 g g g 。回答下列问题。
(*)转动弹簧常数:施加的力矩除以摆的角位移 θ − θ 0 \theta-\theta_0 θ − θ 0 所得的比例常数。
(1) 求摆绕支点 O 的转动惯量 I A I_A I A 。以下各问须直接使用转动惯量符号 I A I_A I A 。
(2) 求摆的运动方程。
(3) 设摆的平衡状态为 θ = θ e q \theta=\theta_{\mathrm{eq}} θ = θ eq ,写出确定该平衡状态的条件式。
(4) 令弹簧自然状态对应的角度为 θ 0 = 0 \theta_0=0 θ 0 = 0 。
1) 求此时表示平衡状态的关系式。
2) θ e q = 0 \theta_{\mathrm{eq}}=0 θ eq = 0 是其中一个平衡状态。现将摆从其极近处 θ ≈ 0 \theta\approx0 θ ≈ 0 以角速度 θ ˙ = 0 \dot\theta=0 θ ˙ = 0 释放。按转动弹簧常数 k k k 的取值,对平衡状态 θ e q = 0 \theta_{\mathrm{eq}}=0 θ eq = 0 附近微小范围内的运动特征进行分类,并分别说明。
(5) 设平衡状态为 θ e q = π / 2 \theta_{\mathrm{eq}}=\pi/2 θ eq = π /2 。
1) 求此时弹簧自然状态对应的角度 θ 0 \theta_0 θ 0 。
2) 从该平衡状态极近处 θ ≈ π / 2 \theta\approx\pi/2 θ ≈ π /2 释放摆,则摆绕 θ e q = π / 2 \theta_{\mathrm{eq}}=\pi/2 θ eq = π /2 作微小振幅的简谐运动。利用运动方程求其固有角频率。
Kai
(1)–(3)
平行軸の定理と、重力・ばねのモーメントより
I A = I + m l 2 , I A θ ¨ = m g l sin θ − k ( θ − θ 0 ) . \boxed{I_A=I+ml^2},\qquad
\boxed{I_A\ddot\theta=mgl\sin\theta-k(\theta-\theta_0)}. I A = I + m l 2 , I A θ ¨ = m g l sin θ − k ( θ − θ 0 ) .
従ってつり合い条件は
m g l sin θ e q = k ( θ e q − θ 0 ) . \boxed{mgl\sin\theta_{\rm eq}=k(\theta_{\rm eq}-\theta_0)}. m g l sin θ eq = k ( θ eq − θ 0 ) .
(4)
m g l sin θ e q = k θ e q . \boxed{mgl\sin\theta_{\rm eq}=k\theta_{\rm eq}}. m g l sin θ eq = k θ eq .
特に θ e q = 0 \theta_{\rm eq}=0 θ eq = 0 は常につり合いである。
θ ( 0 ) = θ i \theta(0)=\theta_i θ ( 0 ) = θ i 、θ ˙ ( 0 ) = 0 \dot\theta(0)=0 θ ˙ ( 0 ) = 0 とする。線形化すると
I A θ ¨ + ( k − m g l ) θ = 0. I_A\ddot\theta+(k-mgl)\theta=0. I A θ ¨ + ( k − m g l ) θ = 0.
k > m g l k>mgl k > m g l :安定な微小振動で、θ = θ i cos ( ω t ) \theta=\theta_i\cos(\omega t) θ = θ i cos ( ω t ) 、ω = ( k − m g l ) / I A \omega=\sqrt{(k-mgl)/I_A} ω = ( k − m g l ) / I A 。
k < m g l k<mgl k < m g l :θ = θ i cosh ( λ t ) \theta=\theta_i\cosh(\lambda t) θ = θ i cosh ( λ t ) 、λ = ( m g l − k ) / I A \lambda=\sqrt{(mgl-k)/I_A} λ = ( m g l − k ) / I A 。微小変位が増大するため不安定である。この式は微小角の範囲で有効。
k = m g l k=mgl k = m g l :線形化では θ ¨ = 0 \ddot\theta=0 θ ¨ = 0 となり、静かに放せば θ = θ i \theta=\theta_i θ = θ i 。ただし、これは一次近似の結論である。元の式では I A θ ¨ = m g l ( sin θ − θ ) = − m g l θ 3 / 6 + O ( θ 5 ) I_A\ddot\theta=mgl(\sin\theta-\theta)=-mgl\theta^3/6+O(\theta^5) I A θ ¨ = m g l ( sin θ − θ ) = − m g l θ 3 /6 + O ( θ 5 ) であり、ポテンシャルは U − U ( 0 ) = m g l ( cos θ − 1 + θ 2 / 2 ) = m g l θ 4 / 24 + O ( θ 6 ) U-U(0)=mgl(\cos\theta-1+\theta^2/2)=mgl\theta^4/24+O(\theta^6) U − U ( 0 ) = m g l ( cos θ − 1 + θ 2 /2 ) = m g l θ 4 /24 + O ( θ 6 ) 。従って原点は非線形には安定で、小さい有限振幅では非調和振動となる。振幅をゼロに近づけると周期は発散する。
(5)
θ e q = π / 2 \theta_{\rm eq}=\pi/2 θ eq = π /2 を代入して
θ 0 = π 2 − m g l k . \boxed{\theta_0=\frac\pi2-\frac{mgl}{k}}. θ 0 = 2 π − k m g l .
θ = π / 2 + η \theta=\pi/2+\eta θ = π /2 + η とおく。sin ( π / 2 + η ) = 1 + O ( η 2 ) \sin(\pi/2+\eta)=1+O(\eta^2) sin ( π /2 + η ) = 1 + O ( η 2 ) とつり合い条件から
I A η ¨ + k η = 0 , ω = k I A . I_A\ddot\eta+k\eta=0,\qquad
\boxed{\omega=\sqrt{\frac{k}{I_A}}}. I A η ¨ + k η = 0 , ω = I A k .
Reference