在空氣阻力下的拋物運動

1、在空氣阻力下的拋物運動空氣阻力一般物體在空氣中運動所受的阻力大小大略與物體對空氣相對速度平方成正比，方向則與相對速度方向相反，因此阻力可表示為：f = -D|v|v其中 v 為物體對空氣的相對速度，如果空氣為靜止，v 就是物體的速度；式中的負號表示阻力與速度方向相反。當 v 為正時，|v|v 為正；當 v 為負時，|v|v 為負（若寫成 v2 則始終為正）。上式中的 D可表示為：其中 ra 為空氣密度，在常溫常壓下約為 1.2 kg/m3；A 為物體面對運動方向上的截面積（cross-sectional area），當物體以較大的截面積面朝向運動方向時，所受的阻力較大，反之則較?。籆 則是一個

2、與物體形狀、表面粗糙度、空氣黏性、等因素相關(guān)的無因次阻力係數(shù)（resistance coefficient）。下表所列為一些狀況下的阻力係數(shù)值。光滑圓球棒球圓柱圓盤人機翼汽車0.50.31.21.10.90.010.4當物體在靜止空氣中做拋物運動時，在水平（x）與鉛垂（y，朝上為正）方向的牛頓第二定律方程式分別為：或(1)(2)其中 ax 與 ay 分別為物體在兩個方向上的加速度；m 為物體質(zhì)量；vx 與 vy 分別為物體在兩個方向上的運動速度，g 為重力加速度。由於物體面對兩個方向的截面積可能不同，故分別以 Dx 與 Dy 區(qū)別?？紤]一個在靜止空氣中做拋物運動的圓球，圓球密度為 r，半徑為

3、r，物體無論面對任何方向的截面積都是於是而物體的質(zhì)量可由體積 V = 4pr3/3 得代入式（1）與（2）可得：(3)(4)顯然這種運動不是等加速度運動，因為加速度會隨速度的大小而改變。所以一般等加速度運動方程式都不適用在這種狀況下。但我們還是可以利用等加速度運動方程式解決這個問題。近似解法如果在一段相當短由 t 到 t+Dt 的時段內(nèi)（意思就是 Dt 很?。?，ax 與 ay 變化不太大，我們可以將原來的加速度方程式稍加修改：(5)(6)(7)(8)上面幾個式子的基本理念如下：1.用 t 到 t+Dt 時段開頭的位置 x(t) 與 y(t) 做為起始位置。2.用 t 到 t+Dt 時段開頭的速

4、度 vx(t) 與 vy(t) 做為起始速度。3.用 t 到 t+Dt 時段開頭的速度 vx(t) 與 vy(t) 計算阻力與加速度 ax(t) 與 ay(t)。4.用 t 到 t+Dt 時段開頭的加速度 ax(t) 與 ay(t) 做為整個時段的等加速度值（只是個近似值，如果 Dt 很小的話，加速度值在這個時段內(nèi)應(yīng)該不會變化太多）。此外我們還需要知道圓球在 t = 0 時的位置 x(0) = x0，y(0) = y0 以及速度 vx(0) = vx0 = v0cosq0，vy(0) = vy0 = v0sinq0，其中 q0 為射角。範例以一個範例做說明：空氣密度 ra：1.2 kg/m3阻

5、力係數(shù) C：0.3圓球密度 r：2000 kg/m3圓球半徑 r：0.025 m時間間距 Dt：0.01 s（希望夠?。┏跛?v0：50 m/s重力加速度 g：9.8 m/s2射角 q0：40 度起始位置 (x0, y0)：(0, 0) m這差不多是一個棒球被擊出的狀況。先將後面會用到的參數(shù)算好：圓球截面積 A： m2D： kg/m圓球體積 V： m3圓球質(zhì)量 m： kgD/m： m-1水平初速 v0x： = 38.30 m/s垂直初速 v0y： = 32.14 m/s這些資料齊備後，就可以逐步計算：t = 0 時在起始時，所有條件均已給定，所以：x(0) = x0 = 0 my(0) = y

6、0 = 0 mvx(0) = v0x = 38.30 m/svy(0) = v0y = 32.14 m/s有了速度後，加速度便可利用式（3）與（4）計算： m/s2 m/s2t = 0.01 s 時利用式（5）與（6）求得： m/s m/s利用式（7）與（8）求得： m m接著利用式（3）與（4）求得 t = 0.01 s 時的加速度，以提供次一個時間（t = 0.02 s）計算使用： m/s2 m/s2t = 0.02 s 以後依前述方法，每隔 0.01 s 計算一次，最後你就會得到位置、速度與加速度對時間的關(guān)係。使用 Excel你或許會說，天呀！照這樣算下去要算到什麼時候才能畫出一個拋物軌

7、跡來。記得二十年前，美國軍方就開發(fā)出一種雷達，可以偵測到敵方射出的砲彈（當然是在空氣中飛行），並預(yù)估砲彈的落點，一方面對己方發(fā)出預(yù)警（當然還得躲得快一點），另一方面可以推估敵方砲位，以便進行反砲擊。如果用剛才的人工計算方法（雖然配合你手邊的計算機），當你還沒算完第一步驟時，敵方砲彈早就落地了，而且敵方的火砲早就跑了。所以這種計算應(yīng)該交給電腦做，你只需要設(shè)計計算的程序來。你可以在所附的 Excel 檔案中看看整個計算方法。在這個檔案中的工作表上方是必須輸入的參數(shù)，輸入的儲存格以底色標示。再下來幾列則是往後計算過程中要用到的參數(shù)，這些參數(shù)是由輸入的數(shù)據(jù)自動推算。再下來近一千列的數(shù)字就是每隔一個 D

8、t 的位置、速度與加速度值，每列的計算結(jié)果都是依據(jù)前面方法由前一列的結(jié)果逐步計算而得。工作表中還有根據(jù)每個時間的 x 與 y 值以散佈圖所繪出的軌跡圖。工作表中的計算還包括不考慮空氣阻力的兩種算法做為比較：一種是純粹套拋物運動公式；另一種方法也是使用前述逐步計算的方法，但不加入空氣阻力相關(guān)的項。這兩種結(jié)果所繪出的拋物軌跡應(yīng)該重疊。這種逐步計算的近似方式在數(shù)學(xué)上稱為數(shù)值方法（numerical method），因為它的結(jié)果必然是數(shù)字，不是代數(shù)公式。如果物體不是圓球呢？前面的例子是針對最簡單的圓球形物體，無論在任何方向截面積都相同。若是像砲彈或飛機之類的物體，還必須知道每個時間彈頭的方向，因為這會

9、影響不同方向的空氣阻力。這個問題涉及更高深的航空動力學(xué)（aerodynamics）知識，根據(jù)物體四周的壓力分布計算物體的指向。微粒在空氣中的運動小微粒（直徑在數(shù)微米或 10-6 m 左右）在空氣中所受的阻力又與一般物體不同。它們所受的阻力與速度成正比：其中 m 為空氣黏滯係數(shù)（viscosity），在常溫常壓下約為 10-4 N-s/m2。小微粒的空氣阻力主要來自於空氣在微粒表面的摩擦；但一般物體空氣阻力主要來自於對其所造成的亂流（因此產(chǎn)生前後壓力差），所以阻力公式不一樣?；蛟S你會想用前面的近似計算法來計算微粒的軌跡，但你會發(fā)現(xiàn)行不通，因為計算幾個步驟後數(shù)字就會大到相當不合理的值（發(fā)散）。這是因為微粒的質(zhì)量太?。ㄅc r 的三次方成正比），雖然微粒所受的阻力也很?。ㄅc r 成正比），但減小的趨勢遠低於質(zhì)量，因此算出來的加速絕對值都超大（可以高達 1010 m/s2），除非你將