邱為鋼

(湖州師范學(xué)院理學(xué)院, 浙江 湖州　313000)

?

半球之間的引力

邱為鋼

(湖州師范學(xué)院理學(xué)院, 浙江 湖州313000)

摘要兩個(gè)不同半球沿對(duì)稱軸方向(第1種情形)或者接觸面分開(第2種情形)，球心間的距離是h.在情形1中，計(jì)算得到了展開到h一次項(xiàng)的引力勢能表達(dá)式.由這個(gè)引力勢能的解析表達(dá)式，得到了兩個(gè)半球間的引力表達(dá)式.在情形2中，計(jì)算得到了展開到h二次項(xiàng)的引力勢能表達(dá)式.當(dāng)大的半球固定時(shí)，小的半球作簡諧振動(dòng)，并給出了周期公式.以上的解析式得到了數(shù)值計(jì)算值的驗(yàn)證.

關(guān)鍵詞萬有引力；半球；周期

基金資助: 國家自然科學(xué)基金(11275067，11475062).

第一顆原子彈的研發(fā)過程中, 有這樣一個(gè)故事: 某次試驗(yàn), 不小心兩個(gè)鈾半球合在一起, 馬上就要重合超過臨界質(zhì)量,將發(fā)生核反應(yīng).有位物理學(xué)家眼捷手快, 把兩個(gè)半球分開,阻止了一次嚴(yán)重事故. 撇開故事，現(xiàn)在只考慮半球之間的引力(平方反比關(guān)系), 掰開兩個(gè)半球的力量有多大? 這就是本文主要討論的問題.

設(shè)兩個(gè)半球的半徑分別為R1,R2，不失一般性，設(shè)R1>R2.起始時(shí)球心重合，重合面為x-y平面.設(shè)大球固定，小球拉開.為計(jì)算方便，只考慮兩種情況，一是沿z軸方向拉開一段距離；二是沿著重合面(假設(shè)接觸面完全光滑)拉開一段距離.第二種情況下，如果拉開距離是小量的話，小球就在平衡位置作小幅度的簡諧振動(dòng).那么，這個(gè)振動(dòng)周期是多少？為便于計(jì)算，假設(shè)球心距離h遠(yuǎn)小于兩球的半徑，計(jì)算過程中只展開到h的一階或兩階小量.

圖1　半球坐標(biāo)示意圖

設(shè)大球在下面，其中一點(diǎn)的球坐標(biāo)是(r1,θ1,φ1)，參數(shù)滿足π/2<θ1<π；小球在上面，其中一點(diǎn)的球坐標(biāo)是(r2,θ2,φ2)，參數(shù)滿足0<θ2<π/2.設(shè)兩個(gè)球的質(zhì)量密度分別是ρ1,ρ2，由萬有引力勢能公式, 兩球之間的引力勢能是

(1)

式(1)中的積分因子有積分變換式

(2)

其中，波矢量k的球坐標(biāo)是(k,θ3,φ3).式(2)中積分因子的分子有以下展開式[1,2]

(3)

(4)

只考慮到h的一階，那么有

(5)

把式(3)～式(5)代入式(2)， 計(jì)算得到(具體計(jì)算過程請(qǐng)參閱文獻(xiàn)[3])引力勢能中包含h一次方項(xiàng)為

(6)

由積分公式[1,2]

(7)

計(jì)算得到

(8)

由此得到沿對(duì)稱軸方向掰開兩個(gè)半球的力是

(9)

下面討論沿接觸面的滑動(dòng)，設(shè)沿x軸方向移動(dòng)h距離，考慮到h的一階和二階:

(10)

把式(3)，式(4)和式(10)代入式(2), 由對(duì)稱性，平移距離h一次方貢獻(xiàn)為零，即沿切合面方向掰開半球所需力為零.引力勢能中包含h二次方項(xiàng)分為V2,0和V2,2(具體計(jì)算過程請(qǐng)參閱文獻(xiàn)[3])，其中V2,0的表達(dá)式是

(11)

其中，

(12)

V2,2的表達(dá)式是

(13)

其中，

(14)

(15)

(16)

實(shí)際計(jì)算中，只取前4項(xiàng)求和，即l=0,1,2,3這4項(xiàng), 這時(shí)兩半球引力勢能的二次項(xiàng)是

(17)

式(17)表明小的半球所感受到的引力勢能是二次項(xiàng)勢能，也意味著小球沿切面作簡諧運(yùn)動(dòng).由二次方勢能中簡諧運(yùn)動(dòng)的周期公式，計(jì)算得到

(18)

表1　兩球引力勢能的二次項(xiàng)比例系數(shù)

綜合利用數(shù)學(xué)物理的積分求和和特殊函數(shù)公式，給出了兩個(gè)半球垂直拉開和平行拉開時(shí)引力勢能的近似解析表達(dá)式，結(jié)果與數(shù)值求值一致.也給出了小球沿切面作簡諧運(yùn)動(dòng)周期的表達(dá)式.

參考文獻(xiàn)

[1]王竹溪，郭敦仁.特殊函數(shù)概論[M].北京: 北京大學(xué)出版社，2000.

[2]Gradshteyn I S， Ryzhik I M. Table of integrals, series, and products[M].7版. 北京: 世界圖書出版公司，2008.

[3]邱為鋼.數(shù)學(xué)物理專題選講[Z]. 湖州： 湖州師范學(xué)院，2014.

審稿意見摘錄：

文中的兩個(gè)例子，吸引力和簡諧頻率可以有解析解，對(duì)于多數(shù)大學(xué)物理的學(xué)生和教師來說，是一件新鮮事.可以作為復(fù)雜的習(xí)題、大作業(yè)之類.文中討論的內(nèi)容，不是當(dāng)下教學(xué)的急所，但可以拓寬大學(xué)物理師生的視野，本身也有點(diǎn)新意.

■

GRAVITATION BETWEEN TWO HEMISPHERES

Qiu Weigang

(School of Science, Huzhou Teachers College, Huzhou, Zhejiang 313000)

AbstractTwo different hemispheres are separated from each other along symmetric axis (case one) or a contact plane (case two). The distance between their centers is denoted as h. In the case one, the gravitational potential energy containing first order of h is derived. The attraction force between two hemispheres is obtained from the analytical expression of the gravitational potential energy. In the case two, the gravitational potential energy containing second order of h is derived. When the bigger hemisphere is fixed, the smaller one’s motion is a simple harmonic vibration. The period time is given. All of those analytical expressions are confirmed by numerical calculation results.

Key wordsgravitation; hemisphere; period

作者簡介:2014-12-23 邱為鋼, 男,副教授，主要從事數(shù)學(xué)物理的教學(xué)和研究.wgqiu@hutc.zj.cn