杜燕茹，李翔宇，韓宏偉，喬 棟

（1.北京理工大學(xué)宇航學(xué)院，北京 100081；2.深空自主導(dǎo)航與控制工信部重點(diǎn)實(shí)驗(yàn)室，北京 100081）

引 言

作為小行星的重要存在形式，雙體小行星在近地小行星、主帶小行星與特洛伊小行星中都已被觀(guān)測(cè)到[1]。在近地小行星中，雙體小行星的數(shù)量約占到其總數(shù)的15±4%[2-3]。由于雙體小行星系統(tǒng)中天體間所具有的獨(dú)特運(yùn)動(dòng)特性可為研究行星系統(tǒng)的演化提供重要線(xiàn)索，因而成為目前小行星探測(cè)的熱點(diǎn)目標(biāo)。雙體小行星模型的不確定性及其附近復(fù)雜多變的動(dòng)力學(xué)環(huán)境對(duì)雙體小行星探測(cè)任務(wù)的設(shè)計(jì)提出了挑戰(zhàn)。

近年來(lái)雙體小行星的穩(wěn)定特性及其附近復(fù)雜的軌道動(dòng)力學(xué)問(wèn)題引起了學(xué)者們的廣泛關(guān)注。Scheeres[4]基于全二體模型在雙體小行星能量與角動(dòng)量守恒條件下，研究了雙體小行星系統(tǒng)旋轉(zhuǎn)運(yùn)動(dòng)與平移運(yùn)動(dòng)中能量與角動(dòng)量的轉(zhuǎn)換關(guān)系，并給出了雙體小行星系統(tǒng)Hill穩(wěn)定性準(zhǔn)則和發(fā)生碰撞的條件。基于球和橢球構(gòu)成的雙體小行星系統(tǒng)，Scheeres提出了評(píng)估平衡狀態(tài)下能量穩(wěn)定性與譜穩(wěn)定性的方法[5]。Bellerose 等基于橢球-球構(gòu)成的全二體問(wèn)題[6-7]，研究了相對(duì)平衡狀態(tài)與穩(wěn)定性問(wèn)題，推導(dǎo)了穩(wěn)定性的解析判據(jù)，并針對(duì)不同構(gòu)型下的雙體小行星系統(tǒng)求解出了穩(wěn)定區(qū)域，研究發(fā)現(xiàn)：與限制性三體問(wèn)題中的穩(wěn)定區(qū)域相比，橢球體-球體系統(tǒng)的穩(wěn)定范圍更小[8]。

本文在前人研究的基礎(chǔ)上，通過(guò)將橢球體與球體模型拓展至雙橢球體模型構(gòu)成的全二體問(wèn)題，研究了雙體小行星系統(tǒng)的平衡態(tài)及穩(wěn)定性，并分析了不同物理參數(shù)對(duì)雙體小行星平衡態(tài)穩(wěn)定性的影響，將為未來(lái)小行星探測(cè)提供理論基礎(chǔ)。

1 雙體小行星動(dòng)力學(xué)與運(yùn)動(dòng)方程

1.1 全二體問(wèn)題

全二體問(wèn)題研究?jī)蓚€(gè)主天體在考慮形狀因素以及質(zhì)量分布的情況下，受相互引力作用而進(jìn)行的運(yùn)動(dòng)，且運(yùn)動(dòng)被限制在一個(gè)平面內(nèi)，如圖1所示。其中，O點(diǎn)為系統(tǒng)的質(zhì)心，P1和P2為構(gòu)成雙體小行星系統(tǒng)的兩個(gè)形狀任意、質(zhì)量分別為M1與M2的主天體。在該雙體小行星系統(tǒng)中，分別引入慣性坐標(biāo)系OXYZ與旋轉(zhuǎn)坐標(biāo)系Oxyz。慣性坐標(biāo)系OXYZ的坐標(biāo)原點(diǎn)為雙體小行星系統(tǒng)質(zhì)心O，X軸指向慣性空間的固定方向，Z軸指向主天體軌道運(yùn)動(dòng)的角動(dòng)量方向，X、Y、Z軸構(gòu)成右手坐標(biāo)系。旋轉(zhuǎn)坐標(biāo)系Oxyz的坐標(biāo)原點(diǎn)為系統(tǒng)質(zhì)心O，x軸由主天體P1質(zhì)心指向主天體P2質(zhì)心，z軸指向主天體軌道運(yùn)動(dòng)的角動(dòng)量方向，x、y、z構(gòu)成右手坐標(biāo)系。此外，分別以小行星的質(zhì)心為原點(diǎn)建立小行星本體坐標(biāo)系，選取本體系的xi、yi和zi(i=1,2)軸與小行星主軸方向一致，其中zi軸垂直于軌道面[9]。假設(shè)雙體小行星系統(tǒng)質(zhì)心O的運(yùn)動(dòng)相對(duì)于小行星間的相對(duì)運(yùn)動(dòng)可忽略，且小行星僅圍繞慣性主軸進(jìn)行自旋，因此，雙體小行星間的相對(duì)運(yùn)動(dòng)可簡(jiǎn)化為用參數(shù)r，?i(i=1,2)，θ表示的四自由度問(wèn)題，其中r為兩小行星質(zhì)心間的距離，?i(i=1,2)為小行星本體系xi(i=1,2)軸與小行星質(zhì)心連線(xiàn)的夾角，θ為質(zhì)心連線(xiàn)相對(duì)于慣性系的旋轉(zhuǎn)角度。

圖1 平面全二體問(wèn)題示意圖Fig.1 Planar full two-body problem

1.2 雙體小行星系統(tǒng)勢(shì)能與動(dòng)能

對(duì)于空間全二體問(wèn)題，慣性系中雙體小行星系統(tǒng)的勢(shì)能為[5]

方程(1)給出了將慣性矩陣展開(kāi)至二階形式的系統(tǒng)勢(shì)能。其中：Tr 表示矩陣的跡；為兩個(gè)小天體質(zhì)心在慣性系中的相對(duì)位置矢量；1,2)為單位化后的慣性矩陣，Mi(i=1,2)為兩個(gè)小天體的質(zhì)量；Ai(i=1,2)為由慣性坐標(biāo)系轉(zhuǎn)換至小天體本體系的旋轉(zhuǎn)矩陣。對(duì)于橢球體模型，假設(shè)質(zhì)量均勻分布且橢球體的半長(zhǎng)軸分別為αi>βi>γi(i=1,2)，則繞主軸的慣性矩可以表示為

根據(jù)系統(tǒng)動(dòng)能的一般性表達(dá)形式，在慣性系OXYZ中可以推導(dǎo)得到雙體小行星系統(tǒng)的動(dòng)能，在假設(shè)主天體的旋轉(zhuǎn)僅圍繞主慣性矩發(fā)生的前提下，系統(tǒng)動(dòng)能在旋轉(zhuǎn)坐標(biāo)系下的表達(dá)式為

系統(tǒng)的能量可以表示為E=T+V。

1.3 雙體小行星系統(tǒng)運(yùn)動(dòng)方程

對(duì)于雙體小行星系統(tǒng)等經(jīng)典物理系統(tǒng)，定義拉格朗日量L為系統(tǒng)動(dòng)能與系統(tǒng)勢(shì)能之差，即

選取兩主天體質(zhì)心間的相對(duì)距離r、雙體小行星系統(tǒng)角速度以及主天體相對(duì)于旋轉(zhuǎn)坐標(biāo)系的旋轉(zhuǎn)角度?i用于表示平面全二體問(wèn)題中的運(yùn)動(dòng)方程。

對(duì)于q=[r,θ,?1,?2]，根據(jù)拉格朗日方程?L/?qi=可得雙體小行星系統(tǒng)的運(yùn)動(dòng)方程為

2 雙橢球體系統(tǒng)平衡態(tài)與穩(wěn)定性分析

2.1 雙橢球體系統(tǒng)平衡態(tài)的狀態(tài)約束

在求解雙體小行星系統(tǒng)平衡態(tài)對(duì)應(yīng)的約束方程時(shí)，角動(dòng)量是其中的一個(gè)關(guān)鍵參數(shù)。分析式(1)與式(3)中的雙體小行星系統(tǒng)動(dòng)能表達(dá)式與勢(shì)能表達(dá)式可發(fā)現(xiàn)，系統(tǒng)勢(shì)能與動(dòng)能的值與雙體小行星系統(tǒng)旋轉(zhuǎn)角度θ無(wú)關(guān)，即雙體小行星系統(tǒng)角動(dòng)量守恒，角動(dòng)量K的表達(dá)式為

求解雙體小行星系統(tǒng)的平衡狀態(tài)通常將全二體問(wèn)題運(yùn)動(dòng)學(xué)方程中關(guān)于時(shí)間求導(dǎo)數(shù)的變量設(shè)為零，可求解平衡狀態(tài)所滿(mǎn)足的約束。因此雙體小行星系統(tǒng)處于平衡狀態(tài)的條件為，系統(tǒng)能量E=T+V需滿(mǎn)足對(duì)各變量求偏導(dǎo)的值均為零。

由于式(1)的勢(shì)能表達(dá)式僅適用于平面全二體問(wèn)題，因此系統(tǒng)的能量表達(dá)式僅適用于平面全二體問(wèn)題。平面全二體問(wèn)題的平衡態(tài)需滿(mǎn)足的約束為

通過(guò)式(6)中的4個(gè)約束可以確定雙體小行星系統(tǒng)平衡狀態(tài)應(yīng)滿(mǎn)足的3個(gè)約束方程為

條件③：Er=0。

本文僅考慮?1=?2=0 的條件下，雙體小行星系統(tǒng)的相對(duì)運(yùn)動(dòng)情況。此時(shí)相對(duì)平衡態(tài)下對(duì)應(yīng)的雙體小行星系統(tǒng)角動(dòng)量平方的值K02與系統(tǒng)能量的大小E0可以分別從式(6)計(jì)算得出。

2.2 雙橢球體系統(tǒng)平衡態(tài)穩(wěn)定性的分析

在進(jìn)行雙體小行星系統(tǒng)平衡態(tài)穩(wěn)定性分析之前，首先需要引入系統(tǒng)零速度狀態(tài)的概念，也即兩天體相對(duì)運(yùn)動(dòng)為零時(shí)狀態(tài)變量構(gòu)成的等式約束。將相對(duì)平衡狀態(tài)r0處的雙體小行星系統(tǒng)能量E0與角動(dòng)量的平方K02帶入零速度狀態(tài)公式，且將平衡狀態(tài)對(duì)應(yīng)的距離r0作為參數(shù)，可以獲得雙體小行星系統(tǒng)運(yùn)動(dòng)可達(dá)范圍的控制方程。由零速度狀態(tài)公式推導(dǎo)得到的控制方程為

其中，A～G為該五次方程的各階系數(shù)，受小行星質(zhì)量分布以及平衡狀態(tài)時(shí)兩小行星的距離r0的影響。計(jì)算驗(yàn)證式(7)的五次方程有二重根r=r0，故可因式分解得到簡(jiǎn)化的控制方程為

利用勞斯-赫爾維茨判據(jù)可得結(jié)論：當(dāng)式(8)中A> 0 時(shí)，該方程的根均分布在復(fù)平面的左半平面；當(dāng)A<0 時(shí)該方程在復(fù)平面的右半平面有一個(gè)根，即存在一個(gè)正實(shí)根。對(duì)該三次不等式進(jìn)行討論研究，可獲得雙體小行星系統(tǒng)平衡態(tài)穩(wěn)定性的3種情況：①運(yùn)動(dòng)可達(dá)范圍無(wú)上界的不穩(wěn)定平衡態(tài)；②運(yùn)動(dòng)可達(dá)范圍有上界的不穩(wěn)定平衡態(tài)；③穩(wěn)定的平衡態(tài)。

由于平衡態(tài)的穩(wěn)定性本質(zhì)上關(guān)注的是能量穩(wěn)定性與Hill 穩(wěn)定性。能量穩(wěn)定性為一類(lèi)基礎(chǔ)的穩(wěn)定性類(lèi)型，決定了系統(tǒng)能否通過(guò)能量耗散的方式達(dá)到能量較低的狀態(tài)。雙體小行星系統(tǒng)具有能量穩(wěn)定性意味著該雙體小行星系統(tǒng)已達(dá)到在給定角動(dòng)量下能量最低的構(gòu)型狀態(tài)。Hill穩(wěn)定性描述雙體小行星系統(tǒng)在相互引力作用下是否具有發(fā)生逃逸能力的穩(wěn)定性，與能量穩(wěn)定性相比，Hill穩(wěn)定性研究對(duì)象不是平衡狀態(tài)附近的運(yùn)動(dòng)，而是關(guān)注雙體小行星系統(tǒng)運(yùn)動(dòng)的絕對(duì)界限。所以以上3 種平衡態(tài)穩(wěn)定性根據(jù)其參數(shù)特性又可歸結(jié)為：①Hill 不穩(wěn)定平衡態(tài)；②Hill 穩(wěn)定、能量不穩(wěn)定平衡態(tài)；③能量穩(wěn)定平衡態(tài)。下面對(duì)這3類(lèi)穩(wěn)定平衡態(tài)分別進(jìn)行分析。

1）Hill不穩(wěn)定平衡態(tài)

A>0 時(shí)，不等式的根均分布于復(fù)平面的左半平面，則式(8)中控制方程始終成立，即r的取值范圍沒(méi)有約束，r的取值范圍無(wú)上界，對(duì)應(yīng)的雙體小行星系統(tǒng)在受擾動(dòng)的情況下可發(fā)生逃逸，且r0處對(duì)應(yīng)的自由能始終為正。

A=0時(shí)，若2A+B≠0，通過(guò)推導(dǎo)分析可知簡(jiǎn)化后控制不等式始終成立。若2A+B=0，則三次不等式可以簡(jiǎn)化為2Gr/r0+G≥0，顯然控制不等式同樣恒成立。

圖2 運(yùn)動(dòng)可達(dá)范圍無(wú)上界的不穩(wěn)定平衡態(tài)Fig.2 Unstable equilibrium，the motion is allowable on all parts of the curve beneath the equilibrium energy

由以上分析可知，當(dāng)A>0與A=0時(shí)相應(yīng)的控制方程始終成立，即對(duì)應(yīng)的雙體小行星系統(tǒng)在受擾動(dòng)的情況下可發(fā)生逃逸，運(yùn)動(dòng)范圍r無(wú)上界，該情況為運(yùn)動(dòng)可達(dá)范圍無(wú)上界的不穩(wěn)定平衡態(tài)。圖2所示為系統(tǒng)能量與相對(duì)距離的關(guān)系，平衡狀態(tài)位于r/r0=1，運(yùn)動(dòng)可達(dá)區(qū)域?yàn)槲挥谄胶鉅顟B(tài)能量E0之下且零能量曲面之上對(duì)應(yīng)的r的取值范圍。

2）Hill穩(wěn)定、能量不穩(wěn)定平衡態(tài)

在A(yíng)<0 的情況下，控制方程有一實(shí)數(shù)根均分布于復(fù)平面的右半平面，其余兩根分布在復(fù)平面的左半平面。假設(shè)控制方程的正實(shí)根為X，若X> 1 即r >r0，此時(shí)r有上界M，且上界M=Xr0應(yīng)滿(mǎn)足M>r0。如果對(duì)平衡態(tài)r/r0=1處的雙體小行星系統(tǒng)施加擾動(dòng)，則運(yùn)動(dòng)可達(dá)區(qū)域可能達(dá)到上界M，即雙體小行星系統(tǒng)內(nèi)兩個(gè)主天體的相對(duì)距離的值可能在r0和M之間變動(dòng)，運(yùn)動(dòng)可達(dá)區(qū)域?yàn)槲挥谄胶鉅顟B(tài)能量之下的曲線(xiàn)部分對(duì)應(yīng)的距離。在該情況下，r/r0=1處的雙體小行星系統(tǒng)的系統(tǒng)總能量為負(fù)值，但為不穩(wěn)定狀態(tài)。雙體小行星系統(tǒng)能量與相對(duì)距離的關(guān)系如圖3所示。

圖3 雙橢球體系統(tǒng)下運(yùn)動(dòng)受限的平衡態(tài)Fig.3 Equilibrium with limited movement in the double ellipsoid system

3）能量穩(wěn)定平衡態(tài)

在A(yíng)< 0 且X≤ 1 即r

圖4 雙橢球體系統(tǒng)能量最小的平衡態(tài)Fig.4 Minimal energy equilibrium in the double ellipsoid system

3 雙體小行星系統(tǒng)物理參數(shù)對(duì)平衡態(tài)穩(wěn)定性影響分析

基于上述關(guān)于雙體小行星系統(tǒng)平衡態(tài)穩(wěn)定性的結(jié)論，研究雙體小行星系統(tǒng)平衡態(tài)受小行星形狀參數(shù)、質(zhì)量分?jǐn)?shù)以及雙體小行星系統(tǒng)平衡態(tài)時(shí)小行星間相對(duì)距離等物理參數(shù)的影響。

假設(shè)小行星Pi為密度均勻一致的雙橢球體模型，Pi的橢球體三軸滿(mǎn)足αi>βi>γi。根據(jù)2.2小節(jié)中對(duì)雙體小行星系統(tǒng)平衡狀態(tài)穩(wěn)定性的分類(lèi)，分別選取形狀參數(shù)一致的雙體小行星系統(tǒng)，并通過(guò)調(diào)整雙體小行星主天體的構(gòu)型與質(zhì)心間距離，分別獲得3種不同穩(wěn)定性情況的系統(tǒng)平衡狀態(tài)。各系統(tǒng)主天體形狀參數(shù)如表1所示。

表1 3種不同穩(wěn)定性情況參數(shù)選取Table 1 Parameters for three different kinds of equilibrium m

隨后對(duì)主天體P2質(zhì)心沿徑向施加5 N·s 的沖量，并以施加沖量后的雙體小行星系統(tǒng)狀態(tài)參數(shù)為狀態(tài)初值，依照第2節(jié)中提供的雙體小行星系統(tǒng)運(yùn)動(dòng)方程沿時(shí)間正向遞推，從而研究雙體小行星受到小擾動(dòng)后的運(yùn)動(dòng)變化。遞推時(shí)長(zhǎng)選取為25 h。

圖5為分別針對(duì)具有能量穩(wěn)定性、Hill 穩(wěn)定且能量不穩(wěn)定性、Hill不穩(wěn)定性的平衡狀態(tài)雙體小行星系統(tǒng)，施加小擾動(dòng)力后主天體質(zhì)心間的運(yùn)動(dòng)變化規(guī)律。由仿真結(jié)果可知，對(duì)于處于平衡狀態(tài)的雙體小行星系統(tǒng)受到擾動(dòng)后主天體質(zhì)心間的距離變化可依據(jù)不同的穩(wěn)定特性分為3 類(lèi)。針對(duì)能量穩(wěn)定系統(tǒng)施加小擾動(dòng)時(shí)，雙體小行星系統(tǒng)主天體質(zhì)心間距離具有保持原平衡狀態(tài)不改變的特性。在平衡狀態(tài)的雙體小行星系統(tǒng)具有Hill穩(wěn)定的情況中，施加小擾動(dòng)使得兩個(gè)主天體質(zhì)心間的相對(duì)距離在一定范圍內(nèi)改變，系統(tǒng)呈現(xiàn)有界的相對(duì)運(yùn)動(dòng)。對(duì)于平衡狀態(tài)下Hill不穩(wěn)定的雙體小行星系統(tǒng)，施加極其微小的擾動(dòng)即可使系統(tǒng)主天體質(zhì)心間的相對(duì)距離達(dá)到無(wú)窮遠(yuǎn)，即發(fā)生逃逸。

圖5 能量穩(wěn)定、Hill穩(wěn)定但能量不穩(wěn)定、Hill不穩(wěn)定雙體小行星系統(tǒng)質(zhì)心距離變化Fig.5 Centroid distance of energy stable，Hill stable but energy unstable，Hill unstable binary asteroid systems

由2.2節(jié)中關(guān)于能量穩(wěn)定性與Hill穩(wěn)定性的結(jié)論，繪制質(zhì)量分?jǐn)?shù)的值在[0,1]之間變化的不同橢球體形狀下曲線(xiàn)。圖6～8 中的距離長(zhǎng)度已相對(duì)于兩橢球體的半長(zhǎng)軸之和α1+α2單位化。將雙球體模型的曲線(xiàn)與不同形狀參數(shù)下雙橢球體的曲線(xiàn)結(jié)果對(duì)比，顯示出不同形狀對(duì)應(yīng)結(jié)果的差異。

圖6～8 展示了考慮兩天體形狀為橢球時(shí)，因形狀參數(shù)的改變而產(chǎn)生的能量穩(wěn)定極限、Hill穩(wěn)定極限以及角動(dòng)量相同時(shí)對(duì)應(yīng)的共軛解的變化。由圖6～8可知，在改變形狀參數(shù)的過(guò)程中穩(wěn)定性界限的特性未發(fā)生質(zhì)的改變，比如對(duì)于Hill穩(wěn)定性的極限，存在某一部分使得對(duì)于任意天體質(zhì)心間距離，雙體小行星系統(tǒng)均具有Hill穩(wěn)定性，對(duì)于能量穩(wěn)定性則不存在這種情況；Hill穩(wěn)定極限對(duì)應(yīng)的兩天體間距離始終小于能量穩(wěn)定極限對(duì)應(yīng)的距離。微小的橢球率可一定程度改變系統(tǒng)的穩(wěn)定性極限，且兩主天體三軸比例不同的情況會(huì)使得對(duì)應(yīng)的能量穩(wěn)定性極限、Hill穩(wěn)定性極限以及共軛解的曲線(xiàn)不再具有對(duì)稱(chēng)性。下面分析小行星形狀參數(shù)變化對(duì)雙體小行星系統(tǒng)平衡態(tài)穩(wěn)定性的影響。假設(shè)主天體形狀參數(shù)滿(mǎn)足βi=γi，選取小行星歸一化后的的長(zhǎng)半軸為α1=6，α2=4，兩個(gè)小行星質(zhì)心間的初始平衡態(tài)距離為r*=10。

圖6 不同形狀參數(shù)天體的能量穩(wěn)定質(zhì)心間距離極限隨質(zhì)量分?jǐn)?shù)的變化Fig.6 Centroid distance limit of energy stable system for a range of different shape bodies

圖7 不同形狀參數(shù)天體Hill穩(wěn)定質(zhì)心間距離極限隨質(zhì)量分?jǐn)?shù)的變化Fig.7 Centroid distance limit of Hill stable system for a range of different shape bodies

圖8 與接觸系統(tǒng)具有相同角動(dòng)量的不同形狀天體的共軛解隨質(zhì)量分?jǐn)?shù)的變化Fig.8 Conjugate stable solutions for a range of different shape bodies.The conjugate solution has the same angular momentum as the solution for a touching system

圖9給出了控制方程的三次項(xiàng)系數(shù)A與控制方程實(shí)數(shù)解Xr的變化曲線(xiàn)。由圖9可知，每一條曲線(xiàn)均具有相似的變化趨勢(shì)，即當(dāng)主天體P1短半軸與長(zhǎng)半軸的比值β1/α1極小時(shí)，控制方程三次項(xiàng)系數(shù)A>0，此時(shí)雙體小行星系統(tǒng)處于可能發(fā)生逃逸的不穩(wěn)定狀態(tài)；隨著β1/α1略微增大，存在小部分區(qū)間可滿(mǎn)足A<0且Xr>1，該情況下雙體小行星系統(tǒng)運(yùn)動(dòng)可達(dá)范圍具有上界M=Xr r*，為能量不穩(wěn)定的相對(duì)平衡狀態(tài)；主星形狀參數(shù)β1/α1繼續(xù)增大，雙體小行星系統(tǒng)的結(jié)果為A<0且Xr<1，即平衡狀態(tài)能量穩(wěn)定。而針對(duì)不同取值的β2/α2可發(fā)現(xiàn)，主天體P2短半軸與長(zhǎng)半軸的比值的增大將使得雙體小行星系統(tǒng)達(dá)到穩(wěn)定狀態(tài)所需的小行星P1的形狀參數(shù)比例β1/α1更小。

下面分析不同小行星質(zhì)量分?jǐn)?shù)對(duì)雙體小行星系統(tǒng)平衡態(tài)穩(wěn)定性的影響。令主天體的三軸比例滿(mǎn)足1:0.75:0.5，選取小行星P2的長(zhǎng)半軸為α2=1，兩個(gè)小行星質(zhì)心間的距離為r*=10，給出控制方程的三次項(xiàng)系數(shù)A與控制方程實(shí)數(shù)解Xr的變化曲線(xiàn)如圖10所示。圖中變化曲線(xiàn)與圖9(b)中變化趨勢(shì)相近，出現(xiàn)該現(xiàn)象的原因?yàn)樵诩僭O(shè)密度相同的條件下，β1/α1比值的增加將導(dǎo)致質(zhì)量分?jǐn)?shù)的增大。并且在質(zhì)量分?jǐn)?shù)極大時(shí)，即小行星P1的質(zhì)量M1遠(yuǎn)大于M2時(shí)，雙體小行星平衡狀態(tài)對(duì)應(yīng)的A與Xr增大，雙體小行星系統(tǒng)開(kāi)始由能量穩(wěn)定的相對(duì)平衡狀態(tài)向Hill 穩(wěn)定且能量不穩(wěn)定狀態(tài)轉(zhuǎn)變。

圖9 β1/α1與β2/α2對(duì)控制方程的三次項(xiàng)系數(shù)A與控制方程實(shí)數(shù)解Xr的影響Fig.9 Effect of β1/α1 and β2/α2 on oefficient A and real solution Xr in the control equation

圖10 質(zhì)量分?jǐn)?shù)對(duì)雙體小行星系統(tǒng)平衡態(tài)穩(wěn)定性的影響Fig.10 Influence of the mass ratio on the stability of equilibrium

最后分析小行星間相對(duì)距離變化對(duì)雙體小行星系統(tǒng)平衡態(tài)穩(wěn)定性的影響。假設(shè)兩個(gè)小行星具有相同的形狀參數(shù)，選取小行星的長(zhǎng)半軸α1=α2=1，為保證兩個(gè)小行星初始狀態(tài)不為碰撞狀態(tài)，質(zhì)心間的距離取值需r*>α1+α2，即在[2,10]的區(qū)間內(nèi)改變雙體小行星系統(tǒng)的質(zhì)心距離。

控制方程的三次項(xiàng)系數(shù)A與控制方程實(shí)數(shù)解Xr的隨r*變化的曲線(xiàn)如圖11所示。由圖11可知，隨著雙星間距離的增加，A為負(fù)值且持續(xù)減小，Xr始終為小于1的正數(shù)，即雙體小行星系統(tǒng)的相對(duì)平衡狀態(tài)將在r*增大的過(guò)程中始終保持穩(wěn)定狀態(tài)。

圖11 小行星間距離r0對(duì)雙體小行星系統(tǒng)平衡態(tài)穩(wěn)定性的影響Fig.11 Influence of the centroid distance on the stability of equilibrium

4 結(jié) 論

本文研究了雙體小行星系統(tǒng)在考慮形狀參數(shù)與質(zhì)量分布的情況下，在相互引力作用運(yùn)動(dòng)的平衡態(tài)與穩(wěn)定性問(wèn)題。文中基于雙橢球構(gòu)成的全二體模型從分析雙體小行星系統(tǒng)的運(yùn)動(dòng)可達(dá)范圍的角度給出了確定雙體小行星系統(tǒng)平衡狀態(tài)的穩(wěn)定性的方法。研究發(fā)現(xiàn)小行星形狀參數(shù)β1/α1的增大可使得雙體小行星系統(tǒng)由具有逃逸潛能的不穩(wěn)定平衡態(tài)，經(jīng)歷運(yùn)動(dòng)范圍存在上界的平衡狀態(tài)，最終轉(zhuǎn)變?yōu)榉€(wěn)定的平衡狀態(tài)；β2/α2增大將使得雙體小行星系統(tǒng)最終達(dá)到穩(wěn)定狀態(tài)對(duì)應(yīng)的β1/α1臨界值更小。雙體小行星系統(tǒng)質(zhì)量分?jǐn)?shù)M1/(M1+M2)的增大對(duì)穩(wěn)定性的影響與形狀參數(shù)產(chǎn)生的影響類(lèi)似，但在M1>>M2時(shí)系統(tǒng)由穩(wěn)定平衡狀態(tài)趨近于有上界的不穩(wěn)定狀態(tài)。系統(tǒng)平衡態(tài)時(shí)小行星質(zhì)心間的距離r0的改變不會(huì)對(duì)系統(tǒng)的穩(wěn)定性產(chǎn)生質(zhì)的影響。這些研究結(jié)論對(duì)于未來(lái)雙體小行星附近探測(cè)器的運(yùn)動(dòng)研究具有重要的理論指導(dǎo)意義。