大阪大学 工学研究科 機械工学専攻 2016年8月実施 機械力学・機構学 問題1
Author
祭音Myyura (refined with GPT 6 Astra)
Description
問題1(解答用紙:赤)
図1に示すような平面連成振り子の微小運動について考える。2つの振り子は、それぞれ長さ ℓ \ell ℓ の棒と質量 m m m のおもりによって構成され、支点から長さ a a a の位置でばね定数 k k k のばねによって互いに接続されている。2つの振り子が鉛直下方向を向いた状態ではばねは自然長を持ち、振り子の支点および棒とばねの接続点は回転支持とする。図2に示すように、左と右の振り子の鉛直方向からの傾斜角をそれぞれ θ 1 \theta_1 θ 1 、θ 2 \theta_2 θ 2 とし、いずれも反時計回りを正とする。また、水平方向に速度 V V V の風が吹いており、左の振り子のおもりには c ( V − ℓ θ ˙ 1 ) 2 c(V-\ell\dot\theta_1)^2 c ( V − ℓ θ ˙ 1 ) 2 の水平右方向の力が作用するものとする。このとき、以下の問いに答えなさい。ただし、棒とばねの質量およびおもりの大きさは無視でき、左の振り子のおもり以外は風の力を受けないものとする。θ i \theta_i θ i (i = 1 , 2 i=1,2 i = 1 , 2 )は微小とし、sin θ i ≒ θ i \sin\theta_i\fallingdotseq\theta_i sin θ i ≒ θ i 、cos θ i ≒ 1 − θ i 2 / 2 \cos\theta_i\fallingdotseq1-\theta_i^2/2 cos θ i ≒ 1 − θ i 2 /2 の近似が成り立つものとする。θ ˙ i \dot\theta_i θ ˙ i は θ i \theta_i θ i の時間微分、
ℓ ∣ θ ˙ 1 ∣ ≪ V < m ( k a 2 + m g ℓ ) ( 2 c ℓ ) 2 , \ell|\dot\theta_1|\ll V<\sqrt{\frac{m(ka^2+mg\ell)}{(\sqrt{2}c\ell)^2}}, ℓ ∣ θ ˙ 1 ∣ ≪ V < ( 2 c ℓ ) 2 m ( k a 2 + m g ℓ ) ,
V > 0 V>0 V > 0 および c > 0 c>0 c > 0 は定数、重力加速度の大きさを g g g とする。
(1) この系に関する運動エネルギー T T T およびポテンシャルエネルギー U U U を求めなさい。
(2) この系に関する運動方程式を求めなさい。
次に、左右の振り子の風を受けて静止した傾斜角をそれぞれ θ 10 \theta_{10} θ 10 、θ 20 \theta_{20} θ 20 とする。左右の振り子の θ 10 \theta_{10} θ 10 、θ 20 \theta_{20} θ 20 からの相対変位角をそれぞれ ϕ 1 \phi_1 ϕ 1 、ϕ 2 \phi_2 ϕ 2 とする。ϕ 1 \phi_1 ϕ 1 、ϕ 2 \phi_2 ϕ 2 は微小とし、反時計回りを正とする。
(3) この系に関する運動方程式を ϕ 1 \phi_1 ϕ 1 、ϕ 2 \phi_2 ϕ 2 を用いて表しなさい。
(4) 右の振り子に、θ 20 \theta_{20} θ 20 を中心とした角振動数 ω \omega ω 、振幅一定の正弦波状の強制変位を与え、定常振動状態に達したものとする。ϕ 1 \phi_1 ϕ 1 の振幅が最大となる ω \omega ω を求めなさい。
题目描述
问题1(使用红色答题纸)
考虑图1所示平面耦合摆的微小运动。两个摆各由长为 ℓ \ell ℓ 的杆和质量为 m m m 的重物组成,在距支点 a a a 的位置通过弹簧常数为 k k k 的弹簧连接。当两个摆均竖直向下时,弹簧处于自然长度;摆的支点以及杆与弹簧的连接点均为转动支承。如图2,左右两摆相对于竖直方向的倾角分别为 θ 1 \theta_1 θ 1 、θ 2 \theta_2 θ 2 ,均以逆时针为正。沿水平方向有速度为 V V V 的风吹来,左摆重物受到向右的水平力 c ( V − ℓ θ ˙ 1 ) 2 c(V-\ell\dot\theta_1)^2 c ( V − ℓ θ ˙ 1 ) 2 。忽略杆和弹簧的质量以及重物的大小,只有左摆重物受到风力。假定 θ i \theta_i θ i (i = 1 , 2 i=1,2 i = 1 , 2 )很小,近似满足 sin θ i ≒ θ i \sin\theta_i\fallingdotseq\theta_i sin θ i ≒ θ i 、cos θ i ≒ 1 − θ i 2 / 2 \cos\theta_i\fallingdotseq1-\theta_i^2/2 cos θ i ≒ 1 − θ i 2 /2 。θ ˙ i \dot\theta_i θ ˙ i 是 θ i \theta_i θ i 的时间导数,且
ℓ ∣ θ ˙ 1 ∣ ≪ V < m ( k a 2 + m g ℓ ) ( 2 c ℓ ) 2 . \ell|\dot\theta_1|\ll V<\sqrt{\frac{m(ka^2+mg\ell)}{(\sqrt{2}c\ell)^2}}. ℓ ∣ θ ˙ 1 ∣ ≪ V < ( 2 c ℓ ) 2 m ( k a 2 + m g ℓ ) .
V > 0 V>0 V > 0 、c > 0 c>0 c > 0 为常数,重力加速度大小为 g g g 。回答下列问题。
(1) 求系统的动能 T T T 和势能 U U U 。
(2) 求系统的运动方程。
接着,设左右两摆在风力作用下静止时的倾角分别为 θ 10 \theta_{10} θ 10 、θ 20 \theta_{20} θ 20 ,各摆相对于该静止角度的角位移分别为 ϕ 1 \phi_1 ϕ 1 、ϕ 2 \phi_2 ϕ 2 。假定 ϕ 1 \phi_1 ϕ 1 、ϕ 2 \phi_2 ϕ 2 很小,以逆时针为正。
(3) 用 ϕ 1 \phi_1 ϕ 1 、ϕ 2 \phi_2 ϕ 2 表示系统的运动方程。
(4) 对右摆施加以 θ 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 の上限により根号内は正であり、有限の正の角振動数で振幅が最大となる。