跳到主要内容

大阪大学 工学研究科 機械工学専攻 2016年8月実施 機械力学・機構学 問題1

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

問題1(解答用紙:赤)

図1に示すような平面連成振り子の微小運動について考える。2つの振り子は、それぞれ長さ ℓ\ell の棒と質量 mm のおもりによって構成され、支点から長さ aa の位置でばね定数 kk のばねによって互いに接続されている。2つの振り子が鉛直下方向を向いた状態ではばねは自然長を持ち、振り子の支点および棒とばねの接続点は回転支持とする。図2に示すように、左と右の振り子の鉛直方向からの傾斜角をそれぞれ θ1\theta_1、θ2\theta_2 とし、いずれも反時計回りを正とする。また、水平方向に速度 VV の風が吹いており、左の振り子のおもりには c(V−ℓθ˙1)2c(V-\ell\dot\theta_1)^2 の水平右方向の力が作用するものとする。このとき、以下の問いに答えなさい。ただし、棒とばねの質量およびおもりの大きさは無視でき、左の振り子のおもり以外は風の力を受けないものとする。θi\theta_i(i=1,2i=1,2)は微小とし、sin⁡θi≒θi\sin\theta_i\fallingdotseq\theta_i、cos⁡θi≒1−θi2/2\cos\theta_i\fallingdotseq1-\theta_i^2/2 の近似が成り立つものとする。θ˙i\dot\theta_i は θi\theta_i の時間微分、

ℓ∣θ˙1∣≪V<m(ka2+mgℓ)(2cℓ)2,\ell|\dot\theta_1|\ll V<\sqrt{\frac{m(ka^2+mg\ell)}{(\sqrt{2}c\ell)^2}},

V>0V>0 および c>0c>0 は定数、重力加速度の大きさを gg とする。

(1) この系に関する運動エネルギー TT およびポテンシャルエネルギー UU を求めなさい。

(2) この系に関する運動方程式を求めなさい。

次に、左右の振り子の風を受けて静止した傾斜角をそれぞれ θ10\theta_{10}、θ20\theta_{20} とする。左右の振り子の θ10\theta_{10}、θ20\theta_{20} からの相対変位角をそれぞれ ϕ1\phi_1、ϕ2\phi_2 とする。ϕ1\phi_1、ϕ2\phi_2 は微小とし、反時計回りを正とする。

(3) この系に関する運動方程式を ϕ1\phi_1、ϕ2\phi_2 を用いて表しなさい。

(4) 右の振り子に、θ20\theta_{20} を中心とした角振動数 ω\omega、振幅一定の正弦波状の強制変位を与え、定常振動状態に達したものとする。ϕ1\phi_1 の振幅が最大となる ω\omega を求めなさい。

原題の図1

原題の図2

题目描述​

问题1(使用红色答题纸)

考虑图1所示平面耦合摆的微小运动。两个摆各由长为 ℓ\ell 的杆和质量为 mm 的重物组成,在距支点 aa 的位置通过弹簧常数为 kk 的弹簧连接。当两个摆均竖直向下时,弹簧处于自然长度;摆的支点以及杆与弹簧的连接点均为转动支承。如图2,左右两摆相对于竖直方向的倾角分别为 θ1\theta_1、θ2\theta_2,均以逆时针为正。沿水平方向有速度为 VV 的风吹来,左摆重物受到向右的水平力 c(V−ℓθ˙1)2c(V-\ell\dot\theta_1)^2。忽略杆和弹簧的质量以及重物的大小,只有左摆重物受到风力。假定 θi\theta_i(i=1,2i=1,2)很小,近似满足 sin⁡θi≒θi\sin\theta_i\fallingdotseq\theta_i、cos⁡θi≒1−θi2/2\cos\theta_i\fallingdotseq1-\theta_i^2/2。θ˙i\dot\theta_i 是 θi\theta_i 的时间导数,且

ℓ∣θ˙1∣≪V<m(ka2+mgℓ)(2cℓ)2.\ell|\dot\theta_1|\ll V<\sqrt{\frac{m(ka^2+mg\ell)}{(\sqrt{2}c\ell)^2}}.

V>0V>0、c>0c>0 为常数,重力加速度大小为 gg。回答下列问题。

(1) 求系统的动能 TT 和势能 UU。

(2) 求系统的运动方程。

接着,设左右两摆在风力作用下静止时的倾角分别为 θ10\theta_{10}、θ20\theta_{20},各摆相对于该静止角度的角位移分别为 ϕ1\phi_1、ϕ2\phi_2。假定 ϕ1\phi_1、ϕ2\phi_2 很小,以逆时针为正。

(3) 用 ϕ1\phi_1、ϕ2\phi_2 表示系统的运动方程。

(4) 对右摆施加以 θ20\theta_{20} 为中心、角频率为 ω\omega、振幅恒定的正弦强制位移,并假定系统已达到稳态振动。求使 ϕ1\phi_1 振幅最大的 ω\omega。

Kai​

(1)​

微小角の二次まで残すと、ばねの伸びは a(θ2−θ1)a(\theta_2-\theta_1) である。よって

T=mℓ22(θ˙12+θ˙22),\boxed{T=\frac{m\ell^2}{2}(\dot\theta_1^2+\dot\theta_2^2)},
U=mgℓ2(θ12+θ22)+ka22(θ1−θ2)2.\boxed{U=\frac{mg\ell}{2}(\theta_1^2+\theta_2^2) +\frac{ka^2}{2}(\theta_1-\theta_2)^2}.

(2)​

右向きの風力は左振り子の正方向にモーメントを与える。速度の一次までとれば

Q1≃cℓ(V−ℓθ˙1)2≃cV2ℓ−2cVℓ2θ˙1,Q2=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.

K=mgℓ+ka2K=mg\ell+ka^2 とおくと、Lagrange の方程式は

mℓ2θ¨1+2cVℓ2θ˙1+Kθ1−ka2θ2=cV2ℓ,\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θ¨2+Kθ2−ka2θ1=0.\boxed{m\ell^2\ddot\theta_2+K\theta_2-ka^2\theta_1=0}.

(3)​

静止角は

Kθ10−ka2θ20=cV2ℓ,Kθ20−ka2θ10=0K\theta_{10}-ka^2\theta_{20}=cV^2\ell, \qquad K\theta_{20}-ka^2\theta_{10}=0

を満たす。参考として

θ10=cV2ℓKK2−k2a4,θ20=cV2ℓka2K2−k2a4.\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}.

θi=θi0+ϕi\theta_i=\theta_{i0}+\phi_i を (2) に代入して静的つり合いを差し引くと

mℓ2ϕ¨1+2cVℓ2ϕ˙1+Kϕ1−ka2ϕ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ϕ¨2+Kϕ2−ka2ϕ1=0.\boxed{m\ell^2\ddot\phi_2+K\phi_2-ka^2\phi_1=0}.

(4)​

右振り子の運動は外部から指定されるので、ϕ2=Bcos⁡ωt\phi_2=B\cos\omega t として左側の式を解く。定常振幅を A1A_1 とすると

A1B=ka2(K−mℓ2ω2)2+(2cVℓ2ω)2.\frac{A_1}{B} =\frac{ka^2}{\sqrt{(K-m\ell^2\omega^2)^2 +(2cV\ell^2\omega)^2}}.

分母の二乗を ω2\omega^2 の関数として最小化すれば

−2mℓ2(K−mℓ2ω2)+4c2V2ℓ4=0.-2m\ell^2(K-m\ell^2\omega^2)+4c^2V^2\ell^4=0.

よって

ω=ka2+mgℓmℓ2−2c2V2m2.\boxed{\omega=\sqrt{\frac{ka^2+mg\ell}{m\ell^2} -\frac{2c^2V^2}{m^2}}}.

原題の VV の上限により根号内は正であり、有限の正の角振動数で振幅が最大となる。