大阪大学 工学研究科 機械工学専攻 2017年度 機械力学・機構学 問題1
Author
祭音Myyura (refined with GPT 6 Astra)
Description
図1の平面連成振り子を考える。2 つの振り子は、それぞれ長さ ℓ \ell ℓ の棒と質量 m m m のおもりからなり、支点から距離 a a a の位置でばね定数 k k k のばねにより接続されている。両振り子が鉛直下向きのとき、ばねは自然長である。振り子の支点と棒・ばねの接続点は回転支持とする。図2のように左右の振り子の鉛直方向からの傾斜角を θ 1 , θ 2 \theta_1,\theta_2 θ 1 , θ 2 とし、いずれも反時計回りを正とする。
水平方向に速さ V V V の風が吹き、左のおもりにのみ、水平方向右向きの力 c ( V − ℓ θ ˙ 1 ) 2 c(V-\ell\dot\theta_1)^2 c ( V − ℓ θ ˙ 1 ) 2 が作用する。棒とばねの質量、おもりの大きさは無視する。各角は微小で、sin θ i ≃ θ i \sin\theta_i\simeq\theta_i sin θ i ≃ θ i 、cos θ i ≃ 1 − θ i 2 / 2 \cos\theta_i\simeq1-\theta_i^2/2 cos θ i ≃ 1 − θ i 2 /2 と近似できる。V , c > 0 V,c>0 V , c > 0 は定数、重力加速度を g g g とし、
∣ ℓ θ ˙ 1 ∣ ≪ V < m ( k a 2 + m g ℓ ) ( 2 c ℓ ) 2 |\ell\dot\theta_1|\ll V
<\sqrt{\frac{m(ka^2+mg\ell)}{(\sqrt2c\ell)^2}} ∣ ℓ θ ˙ 1 ∣ ≪ V < ( 2 c ℓ ) 2 m ( k a 2 + m g ℓ )
を仮定する。
系の運動エネルギー T T T とポテンシャルエネルギー U U U を求めよ。
系の運動方程式を求めよ。
風を受けて静止したときの角度を θ 10 , θ 20 \theta_{10},\theta_{20} θ 10 , θ 20 とし、それらからの相対変位角を ϕ 1 , ϕ 2 \phi_1,\phi_2 ϕ 1 , ϕ 2 とする。ϕ 1 , ϕ 2 \phi_1,\phi_2 ϕ 1 , ϕ 2 も微小で、反時計回りを正とする。
ϕ 1 , ϕ 2 \phi_1,\phi_2 ϕ 1 , ϕ 2 を用いて運動方程式を表せ。
右の振り子に、θ 20 \theta_{20} θ 20 を中心とする振幅一定・角振動数 ω \omega ω の正弦波状の角運動を与え、定常振動に達したとする。ϕ 1 \phi_1 ϕ 1 の振幅が最大となる ω \omega ω を求めよ。
题目描述
图1的两个平面摆各由长为 ℓ \ell ℓ 的无质量杆及质量为 m m m 的物体组成,在距支点 a a a 处由刚度为 k k k 的弹簧连接。两摆均竖直向下时弹簧为自然长度,摆支点及杆与弹簧的连接均可转动。左右摆相对竖直方向的角度为 θ 1 , θ 2 \theta_1,\theta_2 θ 1 , θ 2 ,均以逆时针为正。
水平风速为 V V V ,仅左侧物体受到向右的力 c ( V − ℓ θ ˙ 1 ) 2 c(V-\ell\dot\theta_1)^2 c ( V − ℓ θ ˙ 1 ) 2 。忽略杆和弹簧的质量以及物体尺寸。角度微小,采用 sin θ i ≃ θ i \sin\theta_i\simeq\theta_i sin θ i ≃ θ i 、cos θ i ≃ 1 − θ i 2 / 2 \cos\theta_i\simeq1-\theta_i^2/2 cos θ i ≃ 1 − θ i 2 /2 。V , c > 0 V,c>0 V , c > 0 为常量,重力加速度为 g g g ,且
∣ ℓ θ ˙ 1 ∣ ≪ V < m ( k a 2 + m g ℓ ) ( 2 c ℓ ) 2 . |\ell\dot\theta_1|\ll V
<\sqrt{\frac{m(ka^2+mg\ell)}{(\sqrt2c\ell)^2}}. ∣ ℓ θ ˙ 1 ∣ ≪ V < ( 2 c ℓ ) 2 m ( k a 2 + m g ℓ ) .
求动能 T T T 和势能 U U U 。
求运动方程。
受风作用静止时的角度为 θ 10 , θ 20 \theta_{10},\theta_{20} θ 10 , θ 20 ,相对它们的小角位移为 ϕ 1 , ϕ 2 \phi_1,\phi_2 ϕ 1 , ϕ 2 ,仍以逆时针为正。
用 ϕ 1 , ϕ 2 \phi_1,\phi_2 ϕ 1 , ϕ 2 表示运动方程。
强制右摆绕 θ 20 \theta_{20} θ 20 作振幅固定、角频率为 ω \omega ω 的正弦角运动。系统达到稳态后,求使左摆 ϕ 1 \phi_1 ϕ 1 的振幅最大的 ω \omega ω 。
Kai
(1)
微小角の二次まで残すと、ばねの伸びは a ( θ 2 − θ 1 ) a(\theta_2-\theta_1) a ( θ 2 − θ 1 ) である。よって
T = m ℓ 2 2 ( θ ˙ 1 2 + θ ˙ 2 2 ) , \boxed{T=\frac{m\ell^2}{2}(\dot\theta_1^2+\dot\theta_2^2)}, T = 2 m ℓ 2 ( θ ˙ 1 2 + θ ˙ 2 2 ) ,
U = m g ℓ 2 ( θ 1 2 + θ 2 2 ) + k a 2 2 ( θ 1 − θ 2 ) 2 . \boxed{U=\frac{mg\ell}{2}(\theta_1^2+\theta_2^2)
+\frac{ka^2}{2}(\theta_1-\theta_2)^2}. U = 2 m g ℓ ( θ 1 2 + θ 2 2 ) + 2 k a 2 ( θ 1 − θ 2 ) 2 .
(2)
右向きの風力は左振り子の正方向にモーメントを与える。速度の一次までとれば
Q 1 ≃ c ℓ ( V − ℓ θ ˙ 1 ) 2 ≃ c V 2 ℓ − 2 c V ℓ 2 θ ˙ 1 , Q 2 = 0. Q_1\simeq c\ell(V-\ell\dot\theta_1)^2
\simeq cV^2\ell-2cV\ell^2\dot\theta_1,
\qquad Q_2=0. Q 1 ≃ c ℓ ( V − ℓ θ ˙ 1 ) 2 ≃ c V 2 ℓ − 2 c V ℓ 2 θ ˙ 1 , Q 2 = 0.
K = m g ℓ + k a 2 K=mg\ell+ka^2 K = m g ℓ + k a 2 とおくと、Lagrange の方程式は
m ℓ 2 θ ¨ 1 + 2 c V ℓ 2 θ ˙ 1 + K θ 1 − k a 2 θ 2 = c V 2 ℓ , \boxed{m\ell^2\ddot\theta_1+2cV\ell^2\dot\theta_1
+K\theta_1-ka^2\theta_2=cV^2\ell}, m ℓ 2 θ ¨ 1 + 2 c V ℓ 2 θ ˙ 1 + K θ 1 − k a 2 θ 2 = c V 2 ℓ ,
m ℓ 2 θ ¨ 2 + K θ 2 − k a 2 θ 1 = 0 . \boxed{m\ell^2\ddot\theta_2+K\theta_2-ka^2\theta_1=0}. m ℓ 2 θ ¨ 2 + K θ 2 − k a 2 θ 1 = 0 .
(3)
静止角は
K θ 10 − k a 2 θ 20 = c V 2 ℓ , K θ 20 − k a 2 θ 10 = 0 K\theta_{10}-ka^2\theta_{20}=cV^2\ell,
\qquad K\theta_{20}-ka^2\theta_{10}=0 K θ 10 − k a 2 θ 20 = c V 2 ℓ , K θ 20 − k a 2 θ 10 = 0
を満たす。参考として
θ 10 = c V 2 ℓ K K 2 − k 2 a 4 , θ 20 = c V 2 ℓ k a 2 K 2 − k 2 a 4 . \theta_{10}=\frac{cV^2\ell K}{K^2-k^2a^4},\qquad
\theta_{20}=\frac{cV^2\ell ka^2}{K^2-k^2a^4}. θ 10 = K 2 − k 2 a 4 c V 2 ℓ K , θ 20 = K 2 − k 2 a 4 c V 2 ℓ k a 2 .
θ i = θ i 0 + ϕ i \theta_i=\theta_{i0}+\phi_i θ i = θ i 0 + ϕ i を (2) に代入して静的つり合いを差し引くと
m ℓ 2 ϕ ¨ 1 + 2 c V ℓ 2 ϕ ˙ 1 + K ϕ 1 − k a 2 ϕ 2 = 0 , \boxed{m\ell^2\ddot\phi_1+2cV\ell^2\dot\phi_1
+K\phi_1-ka^2\phi_2=0}, m ℓ 2 ϕ ¨ 1 + 2 c V ℓ 2 ϕ ˙ 1 + K ϕ 1 − k a 2 ϕ 2 = 0 ,
m ℓ 2 ϕ ¨ 2 + K ϕ 2 − k a 2 ϕ 1 = 0 . \boxed{m\ell^2\ddot\phi_2+K\phi_2-ka^2\phi_1=0}. m ℓ 2 ϕ ¨ 2 + K ϕ 2 − k a 2 ϕ 1 = 0 .
(4)
右振り子の運動は外部から指定されるので、ϕ 2 = B cos ω t \phi_2=B\cos\omega t ϕ 2 = B cos ω t として左側の式を解く。定常振幅を A 1 A_1 A 1 とすると
A 1 B = k a 2 ( K − m ℓ 2 ω 2 ) 2 + ( 2 c V ℓ 2 ω ) 2 . \frac{A_1}{B}
=\frac{ka^2}{\sqrt{(K-m\ell^2\omega^2)^2
+(2cV\ell^2\omega)^2}}. B A 1 = ( K − m ℓ 2 ω 2 ) 2 + ( 2 c V ℓ 2 ω ) 2 k a 2 .
分母の二乗を ω 2 \omega^2 ω 2 の関数として最小化すれば
− 2 m ℓ 2 ( K − m ℓ 2 ω 2 ) + 4 c 2 V 2 ℓ 4 = 0. -2m\ell^2(K-m\ell^2\omega^2)+4c^2V^2\ell^4=0. − 2 m ℓ 2 ( K − m ℓ 2 ω 2 ) + 4 c 2 V 2 ℓ 4 = 0.
よって
ω = k a 2 + m g ℓ m ℓ 2 − 2 c 2 V 2 m 2 . \boxed{\omega=\sqrt{\frac{ka^2+mg\ell}{m\ell^2}
-\frac{2c^2V^2}{m^2}}}. ω = m ℓ 2 k a 2 + m g ℓ − m 2 2 c 2 V 2 .
原題の V V V の上限により根号内は正であり、有限の正の角振動数で振幅が最大となる。