摘要：“子彈射擊木塊”是高中物理涉及牛頓運(yùn)動(dòng)定律知識(shí)、能量 、動(dòng)量知識(shí)等的綜合題，本文結(jié)合圖形對(duì)其中兩個(gè)木塊用彈簧鏈接的問(wèn)題進(jìn)行深入分析，并結(jié)合高中數(shù)學(xué)微分知識(shí)建模與求解，利用MATLAB編程實(shí)現(xiàn)，最后用一個(gè)實(shí)例用圖形闡述了兩個(gè)木塊的運(yùn)動(dòng)情況，結(jié)果可靠、形象生動(dòng).

關(guān)鍵詞：子彈；木塊；物理模型；數(shù)學(xué)模型；建模

作者簡(jiǎn)介：程舒揚(yáng)（2001-），女，山東青島，高中二年級(jí)，理科生.

子彈射擊木塊問(wèn)題是高三綜合復(fù)習(xí)時(shí)的經(jīng)典物理題型，這類(lèi)題型往往涉及動(dòng)量定理、動(dòng)量知識(shí)、能量守恒和轉(zhuǎn)化的知識(shí)的綜合運(yùn)用，研究對(duì)象多、過(guò)程復(fù)雜、綜合性強(qiáng)；尤其是可與彈簧或細(xì)繩結(jié)合，組成涉及到力與運(yùn)動(dòng)、彈簧形變等更加復(fù)雜的問(wèn)題.不少學(xué)生理解起來(lái)非常困難.高三同學(xué)已經(jīng)完成了高中數(shù)學(xué)中有關(guān)簡(jiǎn)單微分和積分知識(shí)的學(xué)習(xí)，因此本文利用圖形和數(shù)學(xué)建模方法來(lái)對(duì)其進(jìn)行闡述.

1問(wèn)題描述

如圖1所示，即為該題型一類(lèi)常見(jiàn)的習(xí)題：中間用彈簧相連接的兩個(gè)木塊靜止于光滑的水平面上，被一質(zhì)量為m、速度為v0的子彈擊中后，子彈留在第一個(gè)木塊內(nèi)，試分析彈簧壓縮到最短或拉伸至最長(zhǎng)時(shí)，兩木塊的運(yùn)行速度關(guān)系.

該類(lèi)題目一般會(huì)假設(shè)在子彈擊中后，在兩個(gè)木塊的運(yùn)動(dòng)過(guò)程中彈簧始終處于彈性限度內(nèi)，且木塊 A、B 不會(huì)發(fā)生直接碰撞.本文針對(duì)這種情況，通過(guò)物理模型分析和數(shù)學(xué)建模求解，詳細(xì)闡述當(dāng)彈簧壓縮到最短或拉伸至最長(zhǎng)時(shí)，兩木塊的運(yùn)行速度相等的關(guān)系.

2物理模型

第一種情況，無(wú)論子彈質(zhì)量、速度、兩木塊大小，運(yùn)動(dòng)始終在彈簧的彈性范圍內(nèi).

這種情況下，子彈射入第一個(gè)木塊后留在木塊內(nèi)，可視為質(zhì)量為（m+M1）的一個(gè)整體，忽略子彈在木塊內(nèi)的摩擦（以下同），以子彈初速度方向?yàn)檎较颍琩x為彈簧實(shí)際長(zhǎng)度與自由長(zhǎng)度之差，被壓縮時(shí)dx<0，被拉伸時(shí)dx>0，分析如下.

（1）在初始時(shí)刻，第一個(gè)木塊（包括子彈，以下同）的速度由機(jī)械能守恒定律獲得，為最大v1max；彈簧不發(fā)生形變dx=0，兩個(gè)木塊均不受彈簧力F1=0，F(xiàn)2=0，第二個(gè)木塊的初始速度v2=0，最小，系統(tǒng)的運(yùn)動(dòng)情況如圖2中的①所示.

（2）第一塊木塊向前滑動(dòng)，壓縮彈簧發(fā)生形變，給第二個(gè)木塊一個(gè)推動(dòng)力，第二個(gè)木塊也發(fā)生滑動(dòng)；由于v1大于v2，相同時(shí)間內(nèi)，第一木塊運(yùn)行的位移大于第二個(gè)木塊的位移，二者間距變小，彈簧被壓縮，發(fā)生壓縮形變dx，彈簧進(jìn)入壓縮階段；在該階段內(nèi)，第一個(gè)木塊受到與運(yùn)動(dòng)方向相反的彈簧力為F1，運(yùn)行速度由最大v1max逐漸變慢，用v1↓表示，第二個(gè)木塊受到的彈簧力與第一個(gè)木塊受到的彈簧力大小相等、方向相反，為F2=-F1，運(yùn)行速度逐漸變快v2↑，系統(tǒng)的運(yùn)動(dòng)情況如圖2②所示.

（3）彈簧繼續(xù)被壓縮，彈簧形變繼續(xù)增大，兩個(gè)木塊受到相反的彈簧力繼續(xù)增大，v1更快得變小，v2更快地增加，當(dāng)v1=v2時(shí)，兩木塊之間沒(méi)有相對(duì)位移，彈簧壓縮形變最大即彈簧最短，系統(tǒng)的運(yùn)動(dòng)情況如圖2③所示.

（4）然后v1小于v2，相同時(shí)間內(nèi)，第一木塊運(yùn)行的位移小于第二個(gè)木塊的位移，二者間距拉大，彈簧由最短狀態(tài)逐漸被拉伸，但依然處于壓縮狀態(tài)，兩木塊的受力方向不變，v1繼續(xù)變小，v2繼續(xù)增加，系統(tǒng)的運(yùn)動(dòng)情況如圖2④所示.

（5）當(dāng)v1=0時(shí)，v2達(dá)到最大，彈簧恢復(fù)到最初的自由狀態(tài)，兩木塊受到的彈簧力為0，系統(tǒng)的運(yùn)動(dòng)情況如圖2中的⑤所示.

（6）第一個(gè)木塊不動(dòng)，第二個(gè)木塊速度最大，二者發(fā)生相對(duì)運(yùn)動(dòng)，間距變大，彈簧進(jìn)入拉伸階段：此時(shí)，彈簧的拉力使第一個(gè)木塊的速度v1增加，使第二個(gè)木塊的速度v2減少，系統(tǒng)的運(yùn)動(dòng)情況如圖2中的⑥所示.

（7）當(dāng)?shù)诙蝪1=v2時(shí)，彈簧拉伸至最長(zhǎng)，系統(tǒng)的運(yùn)動(dòng)情況如圖2中的⑦所示.

（8）此后v1>v2，彈簧形變由最長(zhǎng)慢慢減小，但依然處于拉伸階段，v1受到拉力作用v1繼續(xù)增大，v2繼續(xù)減小，系統(tǒng)的運(yùn)動(dòng)情況如圖2中的⑧所示.

（9）當(dāng)v2減小到0時(shí)，v1達(dá)到最大，彈簧再次恢復(fù)到自由狀態(tài)，系統(tǒng)所有的能量轉(zhuǎn)變?yōu)榈谝粋€(gè)木塊的動(dòng)能，v1達(dá)到初始的v1max，系統(tǒng)的運(yùn)動(dòng)情況如圖2中的⑨所示.

由圖2上述分析可見(jiàn)，狀態(tài)⑨與初始態(tài)①完全相同，兩木塊的一個(gè)運(yùn)動(dòng)周期結(jié)束，下一個(gè)運(yùn)動(dòng)周期開(kāi)始.在每個(gè)周期內(nèi)，彈簧都經(jīng)過(guò)如下的周期變化：自由狀態(tài)→壓縮→壓縮至最短→壓縮（恢復(fù)）→自由狀態(tài)→拉伸→拉伸至最長(zhǎng)→拉伸（恢復(fù)）→自由狀態(tài).彈簧形變方向（壓縮還是拉伸）不同，兩木塊的受力方向不同，加速度方向也不同，可根據(jù)彈簧形變將周期分為壓縮階段和拉伸階段.

3數(shù)學(xué)模型建立及求解

3.1數(shù)學(xué)模型

根據(jù)題意，豎直方向上的重力與支持力是一對(duì)平衡力，由于在光滑水平面上運(yùn)動(dòng)，因此可以忽略，只考慮水平方向的速度（動(dòng)能）和彈簧作用力（勢(shì)能）相互變化的關(guān)系.

根據(jù)能量守恒定律，在忽略了子彈在木塊內(nèi)的運(yùn)行情況及能量損失后，系統(tǒng)的能量守方程如下：

E=E1+EP+E2（1）

其中，E為系統(tǒng)的總能量，E1為第一個(gè)木塊（含子彈）的動(dòng)能，E2為第二個(gè)木塊的動(dòng)能，EP為彈簧的勢(shì)能.

E=12mv20（2）

E1=12（M1+m）v21（3）

E2=12M2v22（4）

Ep=12K（Δx）2（5）

根據(jù)速度公式可得兩木塊的速度

v1=v10+a1t（6）

v2=v20+a2t（7）

根據(jù)兩木塊相對(duì)運(yùn)動(dòng)關(guān)系可得彈簧的形變量

Δx=S2-S1=（v2-v1）t（8）

根據(jù)胡克定律可得彈簧的作用力

F1=KΔx（9）

F2=-F1（10）

根據(jù)牛頓第二定律可得彈簧作用力和加速的關(guān)系：

F1=（M1+m）a1（11）

F2=M2a2（12）

聯(lián)立式（1）～（12），即可求得在任意時(shí)間點(diǎn)上兩木塊的運(yùn)動(dòng)速度、加速度、受到的彈簧力，彈簧的形變量以及能量轉(zhuǎn)化情況.

3.2模型求解

以上模型中的變量隨時(shí)間一直變化，矢量還存在方向的改變.用常規(guī)的高中知識(shí)是沒(méi)有辦法求解的.但可以借助高中數(shù)學(xué)中學(xué)過(guò)的微積分知識(shí).

假設(shè)在極短的時(shí)間Δt（比如0.001秒）內(nèi)，各變量的大小和方向不變，根據(jù)機(jī)械能守恒獲得第一個(gè)木塊的初速度有

12mv20=12（M1+m）v210（13）

v10=v0m（M1+m）（14）

即可逐步求得彈簧形變、彈簧作用力、兩個(gè)木塊的加速度，從而求得兩個(gè)木塊的運(yùn)行速度，以及能量轉(zhuǎn)化情況.這些計(jì)算可以借助MATLAB編程實(shí)現(xiàn)，計(jì)算步驟如圖3所示.

4實(shí)例討論與分析

假設(shè)有一個(gè)質(zhì)量為20g的子彈以300m/s的速度射入質(zhì)量為15kg的木塊，該木塊用彈簧與第二個(gè)質(zhì)量為5kg的木塊相連，彈簧的勁度系數(shù)為50N/mm， 用上述模型求解，所得的結(jié)果如圖4、圖5所示.

其中圖4描述了彈簧形變量、兩木塊的速度、加速度、彈簧力的變化趨勢(shì)，圖5描述了彈簧形變、兩木塊動(dòng)能以及總動(dòng)能、彈簧勢(shì)能的變化趨勢(shì).由此可見(jiàn)，木塊的運(yùn)動(dòng)情況與前文的物理分析一致.

對(duì)比圖4描述的各條曲線可見(jiàn)，當(dāng)彈簧形變最大（壓縮至最短或拉伸至最長(zhǎng)）時(shí)，兩木塊的速度相等，加速度最大，彈簧作用力最大；對(duì)比圖5描述的各條曲線可見(jiàn)，當(dāng)彈簧形變最大（壓縮至最短或拉伸至最長(zhǎng)）時(shí)，彈簧勢(shì)能最大，總動(dòng)能最小.

當(dāng)子彈的質(zhì)量m、速度v0、兩木塊的質(zhì)量M1和M2，以及彈簧的勁度系數(shù)K變化時(shí)，假定滿足彈性極限范圍時(shí)，圖4和圖5的振幅和頻率發(fā)生變化，變化趨勢(shì)即運(yùn)動(dòng)規(guī)律不變.

5結(jié)論

本文通過(guò)結(jié)合高中數(shù)學(xué)的微分知識(shí)，利用數(shù)學(xué)建模的思想來(lái)探討高中物理“子彈射擊木塊”問(wèn)題，通過(guò)繪制木塊運(yùn)動(dòng)的物理分析圖以及運(yùn)動(dòng)參數(shù)的變化趨勢(shì)圖，形象生動(dòng)地分析了兩個(gè)被彈簧鏈接的木塊在被子彈射中后，它們的運(yùn)行特點(diǎn)，為分析此類(lèi)習(xí)題提供了新思路.

參考文獻(xiàn)：

[1]明德.子彈打木塊的綜合分析[J].物理教師，2015，36（1）：60-62.

[2]翟燕.對(duì)一道“子彈穿木塊”問(wèn)題的討論[J].物理通報(bào)，2015，（3）：41-42.

[3]許冬保.高三復(fù)習(xí)中習(xí)題教學(xué)思維導(dǎo)圖的應(yīng)用——以“子彈擊木塊”問(wèn)題為例[J].物理教學(xué)，2014，36（5）：73-76.

[4]高尚.“子彈打木塊”模型的非動(dòng)量解法[J].物理教師，2013，34（4）：89-90.

[5]陳新.圖象法在動(dòng)量守恒中的應(yīng)用——再談子彈打木塊問(wèn)題[J].中學(xué)物理，2010，28（11）：63-64.