京都大学 理学研究科 物理学・宇宙物理学専攻 2017年8月実施 問題I-1
Author
祭音Myyura (refined with GPT 6 Astra)
Description
図のように、質量 M M M の棒が、その重心 G G G から距離 ℓ \ell ℓ の位置で、x x x 軸に平行な回転軸に垂直に取りつけられている。棒の回転軸まわりの慣性モーメントは I I I で、棒は回転軸のまわりに y z yz yz 平面内を滑らかに回転するものとする。z z z 軸は鉛直方向を向き、z z z の値が減少する方向に重力加速度 g g g で重力がはたらいているものとする。図に示したように、鉛直下向きから測った棒の傾きの角度を θ \theta θ とする。
(1) 回転軸の位置が不動であるとして、棒にはたらく回転軸まわりのトルクを θ \theta θ の関数として求めよ。
(2) 前問と同様、回転軸の位置が不動であるとして、θ \theta θ の微小振動を考える。微小振動の角振動数 ω 0 \omega_0 ω 0 を求めよ。
以下では、回転軸の位置が y y y 方向に加速度 a ( t ) a(t) a ( t ) で運動する場合を考える。
(3) このとき、θ \theta θ の運動を表す運動方程式を書き下せ。
[ヒント]θ \theta θ は回転軸の運動とともに並進する非慣性座標系における角度座標であることに注意せよ。
以下の問には、M , ℓ , g , I M,\ell,g,I M , ℓ , g , I を用いずに 、小問 (2) で定義した ω 0 \omega_0 ω 0 を用いて 答えよ。
(4) 回転軸の運動の加速度が
a ( t ) = ϵ g sin ( Ω t ) (A) a(t)=\epsilon g\sin(\Omega t)\tag{A} a ( t ) = ϵ g sin ( Ω t ) ( A )
となるように、y y y 方向に振動させる。θ \theta θ の振幅が小さいとして、θ \theta θ に関して線形化し、Ω ≠ ω 0 \Omega\ne\omega_0 Ω = ω 0 の場合に、角振動数 Ω \Omega Ω で周期振動する θ \theta θ の特解を求めよ。
(5) 小問 (4) において、ϵ ≪ 1 \epsilon\ll1 ϵ ≪ 1 であっても、Ω \Omega Ω が ω 0 \omega_0 ω 0 に十分近い場合、線形近似が破綻する。その場合についても適用できる近似解を得るために、振動の非線形性を取り入れることを考える。運動方程式において θ \theta θ の3次の非線形項までを取り入れる近似をする。ただし、ϵ \epsilon ϵ が乗じられた項に関しては θ \theta θ の高次の補正をすべて無視するものとする。このとき、θ 1 \theta_1 θ 1 、および、θ 3 \theta_3 θ 3 を定数として、近似解が
θ = θ 1 sin ( Ω t ) + θ 3 sin ( 3 Ω t ) (B) \theta=\theta_1\sin(\Omega t)+\theta_3\sin(3\Omega t)\tag{B} θ = θ 1 sin ( Ω t ) + θ 3 sin ( 3Ω t ) ( B )
という高次の高調波を無視する形で与えられるとして、θ 1 \theta_1 θ 1 を決定する代数方程式を導け。ただし、∣ θ 3 ∣ |\theta_3| ∣ θ 3 ∣ は、たかだか、O ( ∣ θ 1 ∣ 3 ) O(|\theta_1|^3) O ( ∣ θ 1 ∣ 3 ) であると仮定してよいものとする。
题目描述
如上图,质量为 M M M 的棒在距重心 G G G 为 ℓ \ell ℓ 的位置,垂直安装于平行 x x x 轴的转轴上。棒绕该转轴的转动惯量为 I I I ,可以在 y z yz yz 平面内绕转轴光滑转动。z z z 轴沿竖直方向,重力以加速度 g g g 沿 z z z 减小的方向作用。如图,将从竖直向下方向测量的棒的倾角记为 θ \theta θ 。
(1) 设转轴位置固定,求作用于棒、绕转轴的力矩,表示为 θ \theta θ 的函数。
(2) 同前问,设转轴位置固定,考虑 θ \theta θ 的微小振动。求微小振动的角频率 ω 0 \omega_0 ω 0 。
以下考虑转轴位置沿 y y y 方向以加速度 a ( t ) a(t) a ( t ) 运动的情况。
(3) 写出描述 θ \theta θ 运动的运动方程。
提示:注意 θ \theta θ 是在随转轴一起平移的非惯性坐标系中定义的角坐标。
以下各问不要使用 M , ℓ , g , I M,\ell,g,I M , ℓ , g , I ,而应使用 第 (2) 问定义的 ω 0 \omega_0 ω 0 作答。
(4) 使转轴沿 y y y 方向振动,加速度为
a ( t ) = ϵ g sin ( Ω t ) . (A) a(t)=\epsilon g\sin(\Omega t).\tag{A} a ( t ) = ϵ g sin ( Ω t ) . ( A )
假设 θ \theta θ 的振幅很小,对 θ \theta θ 线性化。在 Ω ≠ ω 0 \Omega\ne\omega_0 Ω = ω 0 时,求以角频率 Ω \Omega Ω 周期振动的 θ \theta θ 的一个特解。
(5) 第 (4) 问中,即使 ϵ ≪ 1 \epsilon\ll1 ϵ ≪ 1 ,当 Ω \Omega Ω 充分接近 ω 0 \omega_0 ω 0 时,线性近似也会失效。为得到在这种情况下仍适用的近似解,考虑振动的非线性。将运动方程中关于 θ \theta θ 的非线性项保留到三次,但对乘有 ϵ \epsilon ϵ 的项忽略 θ \theta θ 的所有高阶修正。设 θ 1 , θ 3 \theta_1,\theta_3 θ 1 , θ 3 为常数,近似解以忽略更高次谐波的形式给出:
θ = θ 1 sin ( Ω t ) + θ 3 sin ( 3 Ω t ) . (B) \theta=\theta_1\sin(\Omega t)+\theta_3\sin(3\Omega t).\tag{B} θ = θ 1 sin ( Ω t ) + θ 3 sin ( 3Ω t ) . ( B )
导出确定 θ 1 \theta_1 θ 1 的代数方程。可以假设 ∣ θ 3 ∣ |\theta_3| ∣ θ 3 ∣ 至多为 O ( ∣ θ 1 ∣ 3 ) O(|\theta_1|^3) O ( ∣ θ 1 ∣ 3 ) 。
Kai
(1)
重力による回転軸方向のトルクは
τ x = − M g ℓ sin θ . \boxed{\tau_x=-Mg\ell\sin\theta.} τ x = − M g ℓ sin θ .
(2)
I θ ¨ = − M g ℓ sin θ ≃ − M g ℓ θ I\ddot\theta=-Mg\ell\sin\theta\simeq-Mg\ell\theta I θ ¨ = − M g ℓ sin θ ≃ − M g ℓ θ より
ω 0 = M g ℓ I . \boxed{\omega_0=\sqrt{\frac{Mg\ell}{I}}.} ω 0 = I M g ℓ .
(3)
軸とともに動く座標系では、重心に慣性力 − M a ( t ) e y -Ma(t)\mathbf e_y − M a ( t ) e y が働き、そのトルクは − M ℓ a ( t ) cos θ -M\ell a(t)\cos\theta − M ℓ a ( t ) cos θ 。したがって
I θ ¨ = − M g ℓ sin θ − M ℓ a ( t ) cos θ . \boxed{I\ddot\theta=-Mg\ell\sin\theta-M\ell a(t)\cos\theta.} I θ ¨ = − M g ℓ sin θ − M ℓ a ( t ) cos θ .
(4)
線形化すると θ ¨ + ω 0 2 θ = − ϵ ω 0 2 sin Ω t \ddot\theta+\omega_0^2\theta=-\epsilon\omega_0^2\sin\Omega t θ ¨ + ω 0 2 θ = − ϵ ω 0 2 sin Ω t 。よって
θ p ( t ) = ϵ ω 0 2 Ω 2 − ω 0 2 sin Ω t . \boxed{\theta_{\mathrm p}(t)=\frac{\epsilon\omega_0^2}{\Omega^2-\omega_0^2}\sin\Omega t.} θ p ( t ) = Ω 2 − ω 0 2 ϵ ω 0 2 sin Ω t .
(5)
指定の次数では
θ ¨ + ω 0 2 θ − ω 0 2 6 θ 3 = − ϵ ω 0 2 sin Ω t . \ddot\theta+\omega_0^2\theta-\frac{\omega_0^2}{6}\theta^3
=-\epsilon\omega_0^2\sin\Omega t. θ ¨ + ω 0 2 θ − 6 ω 0 2 θ 3 = − ϵ ω 0 2 sin Ω t .
θ 3 = θ 1 3 sin 3 Ω t + O ( θ 1 5 ) \theta^3=\theta_1^3\sin^3\Omega t+O(\theta_1^5) θ 3 = θ 1 3 sin 3 Ω t + O ( θ 1 5 ) とし、sin 3 u = ( 3 sin u − sin 3 u ) / 4 \sin^3u=(3\sin u-\sin3u)/4 sin 3 u = ( 3 sin u − sin 3 u ) /4 を用いる。基頻成分の係数から
( ω 0 2 − Ω 2 ) θ 1 − ω 0 2 8 θ 1 3 = − ϵ ω 0 2 . (\omega_0^2-\Omega^2)\theta_1-\frac{\omega_0^2}{8}\theta_1^3=-\epsilon\omega_0^2. ( ω 0 2 − Ω 2 ) θ 1 − 8 ω 0 2 θ 1 3 = − ϵ ω 0 2 .
したがって
ω 0 2 8 θ 1 3 + ( Ω 2 − ω 0 2 ) θ 1 − ϵ ω 0 2 = 0. \boxed{\frac{\omega_0^2}{8}\theta_1^3+(\Omega^2-\omega_0^2)\theta_1-\epsilon\omega_0^2=0.} 8 ω 0 2 θ 1 3 + ( Ω 2 − ω 0 2 ) θ 1 − ϵ ω 0 2 = 0.
近共鳴でも振幅を決める三次項が残り、Ω = ω 0 \Omega=\omega_0 Ω = ω 0 では θ 1 3 = 8 ϵ \theta_1^3=8\epsilon θ 1 3 = 8 ϵ となる。
Reference