謝善娟

(寧德師范高等專科學(xué)校物理與電氣工程系,福建寧德 352100)

拋體運動是日常生活中一種常見的運動現(xiàn)象,也是力學(xué)教材中經(jīng)常討論的典型問題.在考慮地球自轉(zhuǎn)時,拋體除受重力作用外,還受科里奧利力作用.在一般的理論力學(xué)教材中[1],科里奧利力的效應(yīng)多以自由落體為例進(jìn)行討論,有關(guān)文獻(xiàn)[2-4]也對拋體運動作了分析.但在方法上通常都是使用傳統(tǒng)的解析法,采用了近似處理.本文將借助數(shù)學(xué)軟件MATLAB,利用數(shù)值計算的方法,求得精確、直觀的結(jié)果,分析地球自轉(zhuǎn)對拋體運動產(chǎn)生的影響.

1 拋體運動方程

1.1 運動微分方程

在考慮地球自轉(zhuǎn)時,質(zhì)量為 m的拋體,其動力學(xué)方程為：

其中ω為地球自轉(zhuǎn)角速度,v是拋體的速度.以北半球緯度為 r的某點O為坐標(biāo)原點,在地球上建立直角坐標(biāo)系 O-xyz,x軸沿地面水平向南,y軸沿地面水平向東,z軸垂直地面指向天頂.于是可得(1)式在O-xyz坐標(biāo)系內(nèi)的投影方程為[1]為：

(2)式是二階線性微分方程組,解析求解過程比較繁瑣,可以將其變換為一階線性微分方程組,然后利用MATLAB軟件進(jìn)行數(shù)值求解.

1.2 拋體運動軌跡

設(shè)在坐標(biāo)原點 O處,以初速 v0、拋射角θ向正東拋出某物體,則初始條件為 t=0時,x=y=z=0,=0,=v0cosθ,=v0sinθ.取緯度 r=40°、v0=400 m/s、θ=30°,用MATLAB求解微分方程組 (3),繪制出忽略地球自轉(zhuǎn)和考慮地球自轉(zhuǎn)兩種情況下的拋體運動軌跡圖 (如圖1所示).該圖說明在考慮地球自轉(zhuǎn)時,由于科里奧利力的作用,拋體運動發(fā)生了偏移,以下針對拋體運動特征進(jìn)行具體分析.

2 拋體運動的偏移量分析

地球自轉(zhuǎn)使拋體的落地點發(fā)生了偏移.偏移量的大小與緯度、拋射角和初始速度的大小有關(guān).

圖1 忽略地球自轉(zhuǎn)和考慮地球自轉(zhuǎn)的拋體運動軌跡

2.1 緯度對偏移量的影響

設(shè) v0=500 m/s,分別取拋射角θ=50°和75°運行程序,可得到圖2和圖3.

圖2 南北偏移量與緯度的關(guān)系

圖3 東西偏移量與緯度的關(guān)系

圖2給出了南北偏移量隨緯度的變化曲線.從圖2中可看出,在兩極南北偏移最大,而在赤道南北偏移為零.并且往正東方向發(fā)射的拋體,在北半球其落地點南偏,在南半球則為北偏.

圖3給出了東西偏移量隨緯度的變化曲線.圖3說明,往正東方向發(fā)射的拋體,其落地點可能出現(xiàn)東偏或西偏.這種偏移的結(jié)果相當(dāng)于拋體的射程增大或減小.該效應(yīng)在赤道最顯著,而在兩極為零.

2.2 拋射角對偏移量的影響

圖4 南北偏移量與拋射角的關(guān)系

圖5 東西偏移量與拋射角的關(guān)系

圖4給出了南北偏移量隨拋射角的變化曲線.從圖中可看出,當(dāng)拋射角θ≈0.95 rad時,即θ≈54°時,南偏值最大.該最大的南偏值和對應(yīng)的拋射角也可以利用MATLAB編程精確求出.

圖5給出了東西偏移量隨拋射角的變化曲線.從圖中可看出,當(dāng)拋射角θ≈1.04 rad時,即θ≈60°時,東西偏移量為零;當(dāng)θ＞60°時,存在西偏現(xiàn)象 (射程減小);當(dāng)θ＜60°時,存在東偏現(xiàn)象 (射程增大).由圖可知,在θ≈0.52rad時,即θ≈30°時,東偏現(xiàn)象最顯著.該結(jié)果與文獻(xiàn) [2]理論分析的結(jié)果相同.

2.3 初始速度大小對偏移量的影響

表1對不同初始速度大小情況下,地球自轉(zhuǎn)對拋體運動的影響進(jìn)行了比較.從表1可以看出,隨著發(fā)射速度的增大,拋體的南偏、東偏以及射高均增大,而且當(dāng)初始速度較大時,偏移相當(dāng)顯著.這也說明了,當(dāng)發(fā)射速度較大時,地球自轉(zhuǎn)是一個必須考慮的因素.從表1還可以看出,地球自轉(zhuǎn)使拋體的飛行時間增長了,但增長量很小.

表1 不同初速大小情況下地球自轉(zhuǎn)對拋體運動的影響 (緯度40°,拋射角30°)

3 結(jié)論

(1)地球自轉(zhuǎn)使拋體的運動軌跡發(fā)生偏移,這種現(xiàn)象在初始速度較大時,相當(dāng)顯著.

(2)拋體的南偏效應(yīng)在兩極最顯著,在赤道為零.往正東方向發(fā)射的拋體,在北半球落地點南偏,在南半球落地點北偏.偏移量隨拋射角的不同而變化,拋射角約54°時,偏移量最大.

(3)拋體的東偏效應(yīng)在赤道最顯著,在兩極為零.往正東方向發(fā)射的拋體,東西偏移的結(jié)果將使射程增大或減小,當(dāng)拋射角約為60°時,射程的凈偏移為零;當(dāng)拋射角大于60°時,射程減小 (西偏);拋射角小于60°時,射程增大 (東偏).且拋射角約為30°時,東偏量最大.

[1]周衍柏.理論力學(xué)教程：第2版 [M].北京：高等教育出版社,1986：254-258.

[2]吳永懋.考慮地球自轉(zhuǎn)時拋體運動的偏移 [J].大學(xué)物理,1990,(2)：5-6.

[3]冷水根.關(guān)于拋射體偏移問題的商榷 [J].力學(xué)與實踐,1998,(1)：46-47.

[4]李拓.自由落體和拋體偏離的另一種解法 [J].大學(xué)物理,2005,(8)：61-63.

[5]胡靜,彭芳麟,管靖.理論力學(xué)計算機模擬 [M].北京：清華大學(xué)出版社,2002：136-139.

[6]張志涌.精通MATLAB6.5版 [M].北京：北京航空航天大學(xué)出版社,2003：247-249.