基于Simulink的圓形限制性三體運動模型可視化

石云鳳，張洪銀，胡宇，楊洪明，梁攀

（1.西南科技大學數學與應用數學系，綿陽621010；2.西南科技大學信息與計算科學系，綿陽 621010）

0 引言

3個質點，在只有萬有引力相互作用下的運動問題，就是三體問題。由于三體運動方程的非線性，使得其運動軌跡非常依賴于初值，初值改變那怕微乎其微，運動也會從周期性轉化為混沌[1]。因此本文旨在通過建立三體問題的數學模型，再利用Simulink仿真軟件，實現對地球、太陽和月球的圓形限制性三體運動的仿真，從而克服由于初值改變使得三體運動軌跡有較大變動的特點，便于觀察三體運動的軌跡，得到三體運動的規律。

1 三體運動模型的建立

1.1 坐標系建立

可視為質點的三個天體相互作圓周運動,取質心旋轉坐標系,且旋轉速度等于兩天體的圓周運動角速度。則質心在此坐標系中處于靜止,三體在直角坐標系中的相對位置示意圖如圖1所示。

1.2 數學模型建立

當所討論的三個天體中，有一個天體的質量與其他兩個天體的質量相比小到可以忽略時，這樣的三體問題稱為限制性三體問題。一般地把小質量的天體稱為無限小質量體，或簡稱小天體，把兩個大質量的天體稱為有限質量體。由于小天體的質量被視為無限小，則它不影響兩個有限質量體的運動。

圖1 三體相對位置示意圖

根據牛頓萬有引力定理和牛頓第二定律，可以視為質點的兩個天體之間產生的萬有引力關系為：

式中M為天體一的質量，m為天體二的質量，G為萬有引力常數，r為兩個天體的距離。故在空間直角坐標系中，天體二在x軸方向上的加速度分量為：

式中，表示天體二與天體一的所在向量與x軸的夾角。同理可得天體二在y軸方向上的加速度分量與z軸方向上的加速度分量為：

各加速度的分量與天體二質心投影在x、y、z軸上的距離關系為：

本文在建立仿真系統時討論太陽、地球和月球三體問題。由于太陽、月球和地球距離較遠且質量分布較均勻，故可將它們視為質點。在仿真系統中，不考慮其他天體對該三體運動問題的影響。以太陽中心作為坐標原點，月球和地球的運動平面為xOz平面，根據右手法則確定y軸，建立空間直角坐標系。設太陽、月球和地球坐標分別為：(0 ,0,0)，(x2,y2,z2),(x3,y3,z3)，質量分別為M1,M2,M3。由于萬有引力常量實際值太小，為了提高仿真效果，在建立仿真系統時可取G=1。

（1）子系統Gravity accelerations的數學模型：

月球的加速度在x，y和z軸方向的分量分別為：

地球的加速度在x，y和z軸方向的分量分別為：

由于地球的質量和月球的質量相對太陽質量來說太小，故太陽的加速度在x，y和z軸方向的分量均可取0，即：

（2）子系統Distance scaling中數學關系：

考慮到月球繞地球公轉，同時為使仿真效果逼真，故對月球的空間坐標加以修正。根據實際情況取地球與月球質量的比例系數為k=49，故可得月球修正后的空間坐標：

2 仿真系統的建立

根據三體運動中建立的微分方程模型，本文利用MATLAB中的Simulink仿真軟件建立仿真系統，如下：

圖2 仿真總圖

圖3 Gravity accelerations子系統中的組件

圖4 Distance scaling子系統中的組件

3 仿真結果與結果分析

通過查閱相關資料，結合實際，在仿真系統中設定初值如下表1。

表1 仿真系統中的初值

最終通過組件的運行計算，可在VR組件部分得到三體運動的仿真軌跡：以下是仿真結果的動態圖形部分時刻的截圖，如圖5。

由圖可清晰看到地球、太陽、月球三個天體在不同時刻的運動狀態，實現了該三體問題的運動軌跡可視化。通過改變相應初值，還可方便看到該三體的運動狀態的改變，非常地直觀形象。

圖5 三體仿真在不同時刻的運動圖形

4 結語

一般三體問題的運動方程為十八階方程,必須得到18個積分才能得到完全解。然而,目前還只能得到三體問題的10個初積分,還遠不能解決三體問題[2]。通過利用Simulink仿真軟件對太陽、地球、月球三體運動問題進行仿真，得到較為直觀的三體運動模型，這種方法對于三體問題的研究有很好的效果。利用這樣的方法，還可以對其他三體問題進行研究,從而得到更為直觀逼真的三體運動過程。

參考文獻：

[1]楊潔，姜付錦，邱為鋼.三體問題的一種初等解[J].大學物理，2017，36（11）：70-72.

[2]徐衛青，陳朋，連磊,陳亮，謝欣欣.三體問題初探[J].中國高新技術企業，2007（05）：243+245.