跳到主要内容

京都大学 理学研究科 物理学・宇宙物理学専攻 2017年度 問題I-1

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

質量 M1,M2M_1,M_2 の二天体が、距離 aa を保ち、共通重心を原点とする慣性系 S∗S^* の x∗y∗x^*y^* 平面内で角速度

ω=G(M1+M2)a3\omega=\sqrt{\frac{G(M_1+M_2)}{a^3}}

で円運動する。両天体を常に xx 軸上に置く回転系 SS を導入し、

(x∗y∗)=(cos⁡ωt−sin⁡ωtsin⁡ωtcos⁡ωt)(xy)\begin{pmatrix}x^*\\y^*\end{pmatrix} =\begin{pmatrix}\cos\omega t&-\sin\omega t\\\sin\omega t&\cos\omega t\end{pmatrix} \begin{pmatrix}x\\y\end{pmatrix}

とする。μi=Mi/(M1+M2)\mu_i=M_i/(M_1+M_2) とすると、二天体の座標は (−μ2a,0),(μ1a,0)(-\mu_2a,0),(\mu_1a,0) である。これらに比べ質量 mm が無視できる第三天体の運動を同じ平面内で考える。座標系と各天体の配置は図のとおりである。

慣性系と回転系における三天体の配置(原題図)

  1. SS 系における第三天体のラグランジアンを書け。ただし r1=(x+μ2a)2+y2r_1=\sqrt{(x+\mu_2a)^2+y^2}、r2=(x−μ1a)2+y2r_2=\sqrt{(x-\mu_1a)^2+y^2} を用いてよい。

  2. x,yx,y 成分の運動方程式を書け。

  3. SS 系で速度・加速度がともに零になる平衡点のうち、xx 軸外の y>0y>0 の点 L4L_4 と y<0y<0 の点 L5L_5 の座標を求めよ。

  4. 平衡点からの微小変位 ξ,η\xi,\eta の線形方程式を

    ξ¨−2ωη˙+aξ+bη=0,η¨+2ωξ˙+bξ+cη=0\ddot\xi-2\omega\dot\eta+a\xi+b\eta=0,\qquad \ddot\eta+2\omega\dot\xi+b\xi+c\eta=0

    とする。ここで a,b,ca,b,c は定数であり、上の天体間距離 aa とは別の係数である。ξ=Aeλt,η=Beλt\xi=Ae^{\lambda t},\eta=Be^{\lambda t} が非自明解を持つ条件を示せ。

  5. L4L_4 では

    a=−34ω2,b=−334(1−2μ1)ω2,c=−94ω2a=-\frac34\omega^2,\quad b=-\frac{3\sqrt3}{4}(1-2\mu_1)\omega^2,\quad c=-\frac94\omega^2

    である。微小運動が不安定とならない条件を示せ。

  6. (4) の方程式から速度に比例する二項を除き、a,b,ca,b,c は (5) のままとすると、微小運動が常に不安定となることを示せ。

  7. L4L_4 で重力ポテンシャルは上に凸であり、重力だけなら常に不安定になる。それでも (5) で不安定とならない場合がある理由を簡潔に述べよ。

题目描述​

质量 M1,M2M_1,M_2 的两天体间距为 aa,绕共同质心以 ω=G(M1+M2)/a3\omega=\sqrt{G(M_1+M_2)/a^3} 作平面圆周运动。旋转系 SS 相对惯性系 S∗S^* 逆时针转过 ωt\omega t,坐标变换为上式。令 μi=Mi/(M1+M2)\mu_i=M_i/(M_1+M_2),两天体在 SS 中位于 (−μ2a,0),(μ1a,0)(-\mu_2a,0),(\mu_1a,0)。第三天体质量 mm 可忽略,也在同一平面内运动。坐标系与天体位置见上图。

  1. 用到两天体的距离 r1,r2r_1,r_2 写出旋转系拉格朗日量。
  2. 写出 x,yx,y 方向运动方程。
  3. 求轴外两个静止平衡点 L4L_4(y>0y>0)、L5L_5(y<0y<0)。
  4. 微扰满足题中两条线性方程,求指数形式非零解的特征方程;其中 a,b,ca,b,c 是线性系数。
  5. 代入给定的 L4L_4 系数,求不发生不稳定运动的条件。
  6. 保持这些系数,去掉速度项,证明平衡点总不稳定。
  7. 解释为何重力单独作用时不稳定,完整旋转系方程却可能稳定。

Kai​

(1)​

慣性系での速度の二乗は x˙2+y˙2+2ω(xy˙−yx˙)+ω2(x2+y2)\dot x^2+\dot y^2+2\omega(x\dot y-y\dot x)+\omega^2(x^2+y^2) なので、

L=m2[x˙2+y˙2+2ω(xy˙−yx˙)+ω2(x2+y2)]+GM1mr1+GM2mr2.\boxed{L=\frac m2\left[\dot x^2+\dot y^2+2\omega(x\dot y-y\dot x)+\omega^2(x^2+y^2)\right] +\frac{GM_1m}{r_1}+\frac{GM_2m}{r_2}.}

(2)​

オイラー・ラグランジュ方程式より

x¨−2ωy˙−ω2x=−GM1(x+μ2a)r13−GM2(x−μ1a)r23,\boxed{\ddot x-2\omega\dot y-\omega^2x =-\frac{GM_1(x+\mu_2a)}{r_1^3}-\frac{GM_2(x-\mu_1a)}{r_2^3},}
y¨+2ωx˙−ω2y=−GM1yr13−GM2yr23.\boxed{\ddot y+2\omega\dot x-\omega^2y =-\frac{GM_1y}{r_1^3}-\frac{GM_2y}{r_2^3}.}

(3)​

y≠0y\ne0 の平衡点では ω2=GM1/r13+GM2/r23\omega^2=GM_1/r_1^3+GM_2/r_2^3。これを xx の式に代入すると r1=r2r_1=r_2、さらに r1=r2=ar_1=r_2=a を得る。したがって二天体と正三角形をなし、

L4,5=(μ1−μ22a, ±32a).\boxed{L_{4,5}=\left(\frac{\mu_1-\mu_2}{2}a,\ \pm\frac{\sqrt3}{2}a\right).}

(4)​

(λ2+ab−2ωλb+2ωλλ2+c)(AB)=0.\begin{pmatrix}\lambda^2+a&b-2\omega\lambda\\b+2\omega\lambda&\lambda^2+c\end{pmatrix} \begin{pmatrix}A\\B\end{pmatrix}=0.

よって非自明解の条件は

λ4+(a+c+4ω2)λ2+ac−b2=0.\boxed{\lambda^4+(a+c+4\omega^2)\lambda^2+ac-b^2=0.}

(5)​

与えられた係数から

λ4+ω2λ2+274μ1μ2ω4=0,λ2=−ω22(1±1−27μ1μ2).\lambda^4+\omega^2\lambda^2+\frac{27}{4}\mu_1\mu_2\omega^4=0, \qquad \lambda^2=-\frac{\omega^2}{2}\left(1\pm\sqrt{1-27\mu_1\mu_2}\right).

正の質量について、異なる純虚数の特性根だけを持ち、すべての線形微小運動が有界になる条件は

μ1μ2<127⟺(M1+M2)2>27M1M2.\boxed{\mu_1\mu_2<\frac1{27} \quad\Longleftrightarrow\quad (M_1+M_2)^2>27M_1M_2.}

指数的に増大する根が存在しない条件だけなら μ1μ2≤1/27\mu_1\mu_2\le1/27 である。ただし等号では λ=±iω/2\lambda=\pm i\omega/\sqrt2 が重根となり、対応する固有空間は各一次元なので、一般解には te±iωt/2te^{\pm i\omega t/\sqrt2} が現れる。有界な線形安定性には等号を含めない。

(6)​

速度項を除くと特性方程式は

λ4−3ω2λ2+274μ1μ2ω4=0,λ2=3ω22(1±1−3μ1μ2).\lambda^4-3\omega^2\lambda^2+\frac{27}{4}\mu_1\mu_2\omega^4=0, \qquad \lambda^2=\frac{3\omega^2}{2}\left(1\pm\sqrt{1-3\mu_1\mu_2}\right).

0<μ1μ2≤1/40<\mu_1\mu_2\le1/4 より両方の λ2\lambda^2 は正である。したがって正の実特性根を持ち、常に不安定である。

(7)​

回転系では遠心力に加えて、速度に依存するコリオリ力が働く。特にコリオリ力による二方向の運動の結合がジャイロ的な安定化をもたらすため、ポテンシャルの曲率だけでは安定性を判定できない。

Reference​