高云峰
(清華大學(xué)航天航空學(xué)院,北京 100084)
科學(xué)研究離不開系統(tǒng)與模型。系統(tǒng)是根據(jù)研究目的而人為界定,由一些有關(guān)聯(lián)的物體或成份組成的有序整體。簡(jiǎn)單地說,系統(tǒng)就是我們關(guān)注或研究的對(duì)象,包含邊界、成份和相互作用。模型是經(jīng)過處理的簡(jiǎn)化系統(tǒng),能體現(xiàn)原系統(tǒng)的本質(zhì)特征,是描述和理解系統(tǒng)的有效工具。真實(shí)世界很復(fù)雜,但可以用簡(jiǎn)化的模型來(lái)描述。模型超越了我們可觀察的范圍,將不可見的、抽象的現(xiàn)象或特征形象化,有利于理解和研究。
模型可以直觀反映系統(tǒng)的特征,特別是系統(tǒng)的內(nèi)在特征,因此有時(shí)模型比實(shí)物更能展示出系統(tǒng)的特點(diǎn),利用模型也更容易進(jìn)行溝通和交流。模型與系統(tǒng)之間存在著某種相似性,除了可以解釋一部分系統(tǒng)的現(xiàn)象,還可以通過模型預(yù)測(cè)系統(tǒng)更多的現(xiàn)象,甚至是還沒有被發(fā)現(xiàn)的現(xiàn)象,從而在驗(yàn)證模型可靠性的同時(shí),增加對(duì)系統(tǒng)的認(rèn)識(shí)。例如,利用太陽(yáng)系模型預(yù)測(cè)海王星的位置是模型預(yù)測(cè)的典型案例。
建立簡(jiǎn)化的模型來(lái)研究復(fù)雜的系統(tǒng)是一種科學(xué)研究的方法。模型包括物理模型(通常是實(shí)物或示意圖)、數(shù)學(xué)模型(通常是數(shù)學(xué)公式)和概念模型(通常是相關(guān)理論)。
在目前的力學(xué)教學(xué)中,平時(shí)的習(xí)題或作業(yè)可以認(rèn)為是已經(jīng)簡(jiǎn)化好的模型,學(xué)生缺乏如何從真實(shí)系統(tǒng)簡(jiǎn)化為模型的訓(xùn)練(在本系列文章(2)中,從“工程問題”到“物理模型”,實(shí)際也是建模的過程)。
在教學(xué)過程中,筆者會(huì)把少量傳統(tǒng)的習(xí)題進(jìn)行改造,從實(shí)際問題或現(xiàn)象出發(fā),給學(xué)生展示如何簡(jiǎn)化建模,如何分析得到結(jié)果并最終解釋所關(guān)心的問題。這樣的內(nèi)容不用多,給學(xué)生開一個(gè)窗口,讓學(xué)生了解一種重要的科學(xué)研究的方法和思路,而不僅僅是會(huì)做作業(yè)。
吉他、小提琴能演奏出優(yōu)美的旋律,本質(zhì)是琴弦振動(dòng)的頻率和幅度通過空氣傳到耳中(我們聽到的旋律與振動(dòng)頻率相關(guān),聲音的強(qiáng)弱與振動(dòng)幅度相關(guān))。
以吉他為例,吉他由共鳴箱、琴頸、琴頭及琴弦等部件組成,每根琴弦粗細(xì)不同,固定在琴橋與琴頭之間,琴頸上有19根長(zhǎng)品柱,第i個(gè)品柱與琴橋的距離為li(i=0,1,...,18) ,可以注意到這些距離從密到稀排列著(圖1)。
圖1 吉他照片
可以對(duì)琴弦振動(dòng)進(jìn)行建模,通過分析來(lái)解釋與吉他相關(guān)的一些問題,包括上述現(xiàn)象。
首先將某根琴弦簡(jiǎn)化為無(wú)質(zhì)量的彈簧,將琴弦質(zhì)量簡(jiǎn)化為質(zhì)點(diǎn)集中在彈簧中點(diǎn),靜止時(shí)琴弦已繃緊,設(shè)張力為T0。當(dāng)質(zhì)點(diǎn)有位移x時(shí),單側(cè)琴弦長(zhǎng)度為d,張力為T(圖2)。
圖2 琴弦振動(dòng)的簡(jiǎn)化模型
下面尋找系統(tǒng)參數(shù)之間的關(guān)系。單側(cè)琴弦長(zhǎng)度為
通常彈撥吉他時(shí),琴弦的位移很小,因此設(shè)x是一階小量。式(1)級(jí)數(shù)展開后保留一階量,有
即琴弦振動(dòng)時(shí)在忽略高階小量的前提下長(zhǎng)度不變。由于張力與琴弦伸長(zhǎng)量有關(guān),因此在振動(dòng)過程中有
忽略重力(張力遠(yuǎn)大于琴弦重量),列寫質(zhì)點(diǎn)振動(dòng)的動(dòng)力學(xué)方程,有
由于位移x是小量,角度θ也是小量,有
將式(3)和式(5)代入式(4),有
設(shè)琴弦截面積為S,單位長(zhǎng)度質(zhì)量為ρ,則振動(dòng)的固有頻率為
下面利用式(7)來(lái)解釋一些現(xiàn)象。
(1)琴頭上有調(diào)緊琴弦的旋鈕,有經(jīng)驗(yàn)的吉他手演奏前要調(diào)整一下旋鈕,通過改變琴弦的張力調(diào)整頻率(圖3)。式(7)表明頻率ω與成正比,合理。
圖3 琴頭中的琴弦及旋鈕
(2)琴弦有粗細(xì)之分(圖3),粗弦用于演奏低音(低頻),細(xì)弦用于演奏高音(高頻)。式(7)表明頻率與成反比,合理。相反如果6根琴弦一樣粗細(xì),高音對(duì)應(yīng)的琴弦張力就明顯增加,容易導(dǎo)致高音琴弦疲勞損壞。
(3)手按住琴弦某處,相當(dāng)于琴頭處的固定端移到手指處的品柱,改變了琴弦的長(zhǎng)度,從而改變了振動(dòng)頻率。我們知道世界上把一組音(八度)分成十二個(gè)半音音程的律制,各相鄰兩音的振動(dòng)頻率之比完全相等,亦稱“十二等程律”。式(7)表明頻率與l成反比,因此圖1中距離的比值也應(yīng)該滿足“十二等程律”,即前后距離之比為實(shí)際品柱與琴橋的距離見表1,相鄰的距離比值見圖4,合理。
表1 吉他各品柱的距離(根據(jù)實(shí)際照片測(cè)量)
圖4 相鄰品柱距離的比值
對(duì)音樂熟悉的讀者可能會(huì)注意到,國(guó)際標(biāo)準(zhǔn)音的頻率(表2)與表1中的部分品柱距離(黑體)竟然十分接近,我自己也感到有些不可思議(該吉他34寸(1寸 = 3.333 33 cm),全長(zhǎng)866 mm,如果是其他品牌,品柱距離會(huì)變化,但比值不變)。
表2 C大調(diào)國(guó)際標(biāo)準(zhǔn)音頻率
(4)實(shí)際吉他彈起來(lái)音色很豐富,而模型只有一個(gè)頻率。原因是模型把質(zhì)量集中為一個(gè)質(zhì)點(diǎn),這與實(shí)際不同。如果采用連續(xù)質(zhì)量的模型,就更為接近實(shí)際的了。
大家可能注意到下雨時(shí)毛毛細(xì)雨與暴雨有一個(gè)明顯區(qū)別:毛毛細(xì)雨是“飄”下來(lái)的(速度很?。?,而暴雨是“砸”下來(lái)的(速度很大)。如何用相關(guān)的力學(xué)知識(shí),建立模型簡(jiǎn)要地解釋這一現(xiàn)象呢?
建立坐標(biāo)系,確定廣義坐標(biāo)y,分析運(yùn)動(dòng)和受力,畫出受力圖(圖5)。
圖5 雨滴受力分析
根據(jù)一般的經(jīng)驗(yàn),物體在空氣中阻力為
其中CD為阻力系數(shù)(常數(shù)),S為物體橫截面積,ρa(bǔ)為空氣密度。
如果假設(shè)雨滴是球形,半徑為r,水的密度為ρw,則有
根據(jù)圖5,速度逐漸增加導(dǎo)致阻力增加,當(dāng)滿足阻力等于重力后,雨滴就一直以勻速運(yùn)動(dòng)(嚴(yán)格說是趨近于勻速運(yùn)動(dòng))。由式(8)和式(9)獲得雨滴的極限速度為
式(10)表明雨滴的速度與其半徑有關(guān),這就解釋了暴雨速度大,毛毛雨速度小的現(xiàn)象。當(dāng)然上述模型假設(shè)雨滴半徑不變化,這不太合理,如果要更精細(xì),可以修正模型。
考慮雨滴在下降過程中不斷吸收周圍的水汽,設(shè)周圍大氣相對(duì)地面靜止,采用變質(zhì)量模型[1],有
雨滴在單位時(shí)間內(nèi)經(jīng)過的空間是一個(gè)圓柱體,該圓柱內(nèi)的水汽有部分被雨滴吸收(圖6),設(shè)吸收率ζ為常數(shù)(ζ = 1表示全部吸收,ζ = 0表示不吸收;吸收率ζ通常會(huì)很小,否則幾滴雨就把運(yùn)動(dòng)路徑上的水汽全吸收了),則
圖6 考慮雨滴的變質(zhì)量問題
將式(12)取極限后與式(11)聯(lián)立,有
將式(9)改寫,有
設(shè)y1=y,y2=y˙,y3=r,由式(13)和式(14)得到標(biāo)準(zhǔn)的一階微分方程組,為
方程(15)耦合在一起,給出參數(shù)及初始條件 (ρw=1000kg/m3,ρa(bǔ)=1kg/m3,r(0)=0.02mm ,H=2000m ,CD=1 , 估 計(jì) 暴 雨 時(shí)ζ=0.005 ,中 雨 時(shí)ζ=0.001 , 小 雨時(shí)ζ=0.0006 ), 可 以 算出相關(guān)的數(shù)值解。
與前面雨滴質(zhì)量和半徑為常數(shù)的模型相比,變質(zhì)量模型中雨滴不存在極限速度,而是隨著高度減低,速度不斷增加(圖7)。雨滴變質(zhì)量模型的計(jì)算結(jié)果與網(wǎng)上查到的數(shù)據(jù)很接近[2]:一般而言,毛毛雨的雨滴直徑為0.5 mm, 而暴雨雨滴最大直徑約為5.5 mm;毛毛雨的收尾速度為2 m/s,而暴雨最大收尾速度為8~9 m/s(圖8)。且同一高度落下的雨滴沿相同的“速度-直徑”曲線變化,這一點(diǎn)還有待深入研究。
圖7 雨滴高度與速度的關(guān)系
圖8 雨滴直徑與速度的關(guān)系
系統(tǒng)與模型是跨學(xué)科概念,可以應(yīng)用在很多學(xué)科,力學(xué)模型主要研究系統(tǒng)的受力與運(yùn)動(dòng)的關(guān)系,可以遷移到其他領(lǐng)域:任何系統(tǒng)都可以建立類似“力與運(yùn)動(dòng)”的關(guān)系,只不過這里的“力”和“運(yùn)動(dòng)”更為廣泛。例如在生物領(lǐng)域,“力”可以是生存壓力,“運(yùn)動(dòng)”可以是種群數(shù)量的變化。
實(shí)際上,通過食物鏈關(guān)系,不同生物之間存在著某種制約關(guān)系(相當(dāng)于分析力學(xué)中的約束),導(dǎo)致種群數(shù)目保持動(dòng)態(tài)平衡。為了更好地理解這種制約及平衡關(guān)系,設(shè)想在一個(gè)理想的封閉草原環(huán)境中,存在野兔與野狼兩種生物。很明顯,野兔多則野狼食物充分,導(dǎo)致野狼數(shù)量增加;而野狼數(shù)量太多,會(huì)通過捕食使野兔數(shù)量減少,反過來(lái)又導(dǎo)致部分野狼不容易找到食物而餓死,繁殖的數(shù)量也會(huì)下降。因此在這種生存競(jìng)爭(zhēng)中,兩種生物要么最終滅亡,要么經(jīng)過一段時(shí)間后兩者數(shù)量達(dá)到穩(wěn)定的狀態(tài)。
這一現(xiàn)象可以用Lotka-Volterra模型來(lái)描述,也稱作獵食者-獵物(predator-prey)模型[3]。改寫后,設(shè)野兔數(shù)量為R(rabbit),野狼的數(shù)量為W(wolf)。上述競(jìng)爭(zhēng)關(guān)系可以表示為
式中λ1是野兔的自然增長(zhǎng)率(假定草場(chǎng)資源充沛),λ1R表示在單位時(shí)間出生的野兔;λ2表示單位時(shí)間內(nèi)野兔被野狼捕獲的比例,λ2RW表示在單位時(shí)間被吃掉的野兔;λ3是野狼的自然死亡率,λ3W表示在單位時(shí)間死亡的野狼;λ4表示轉(zhuǎn)換系數(shù),表示平均吃λ4只野兔才生出1只小野狼。
給出適當(dāng)參數(shù)和初始條件后,計(jì)算模擬發(fā)現(xiàn)野兔和野狼的數(shù)量的確呈現(xiàn)周期性的變化。其典型的情況是:從時(shí)間變化角度看,兩者的數(shù)量振蕩起伏(圖9(a))。如果用水平軸和豎直軸分別表示它們的數(shù)量,會(huì)出現(xiàn)一個(gè)封閉的曲線(圖9(b));當(dāng)然,如果參數(shù)不合適,會(huì)出現(xiàn)野狼把野兔全部吃完,(在封閉的環(huán)節(jié)下)野狼也會(huì)全部餓死的特殊“平衡”解。
圖9 野兔和野狼數(shù)量變化的關(guān)系
由此看出,在不同物種生存競(jìng)爭(zhēng)這一問題中,不同物種的數(shù)目隨時(shí)間在變化著,但是由于存在食物鏈關(guān)系,每個(gè)物種對(duì)另外的物種生存都會(huì)產(chǎn)生影響,從而導(dǎo)致最終達(dá)到一種“穩(wěn)定”的狀態(tài),或動(dòng)態(tài)的平衡。實(shí)際生物界有著復(fù)雜的食物鏈關(guān)系,模型更為復(fù)雜,但是大致仍有簡(jiǎn)單模型的特征:不同物種相互制約,長(zhǎng)期趨于動(dòng)態(tài)的平衡關(guān)系。
系統(tǒng)和模型是重要跨學(xué)科概念,可以應(yīng)用于不同的學(xué)科領(lǐng)域。一般科研都要對(duì)系統(tǒng)進(jìn)行簡(jiǎn)化,獲得模型后,再進(jìn)行深入研究。很多學(xué)生總感到理論與實(shí)際有巨大的鴻溝,最主要原因是缺乏建模的訓(xùn)練。有建模經(jīng)驗(yàn)的學(xué)生可以將復(fù)雜問題(實(shí)際)簡(jiǎn)化后得到相對(duì)容易的模型(理論),然后就可以用不同的手段進(jìn)行處理。
從以上案例中可以得出幾點(diǎn)結(jié)論或啟示。
(1)對(duì)琴弦的建模及分析并不復(fù)雜,得到的結(jié)果卻可以解釋很多相關(guān)現(xiàn)象,并且意外地發(fā)現(xiàn)某種吉他(34寸)的品柱距離與國(guó)際標(biāo)準(zhǔn)音的頻率數(shù)值十分接近。因此簡(jiǎn)單的建模就可以解決很多問題,同時(shí)注意到模型是簡(jiǎn)化的系統(tǒng),該模型解釋不了實(shí)際琴弦音色豐富的問題。
(2)雨滴問題可以簡(jiǎn)單也可以復(fù)雜,均可以解釋雨滴速度與尺寸的關(guān)系,這表明模型的復(fù)雜程度可以根據(jù)情況而定。注意到分析過程中有大量假設(shè),如形狀是球形、阻力與速度平方成比例、吸收率是常數(shù)等,這些可能與實(shí)際存在差距,也是模型的重要特點(diǎn):模型是簡(jiǎn)化的系統(tǒng),在簡(jiǎn)化時(shí)盡量有依據(jù)。對(duì)于變質(zhì)量雨滴,計(jì)算中的參數(shù)參考了網(wǎng)上的數(shù)據(jù),如云層高度、初始雨滴大小等,最終得到的雨滴速度和直徑與網(wǎng)上資料十分接近。
(3)生物競(jìng)爭(zhēng)的案例擴(kuò)展了學(xué)生的視野,讓學(xué)生了解建模的思想可以用于各個(gè)領(lǐng)域。同時(shí)有利于學(xué)生了解更廣泛的約束概念,以及“力與運(yùn)動(dòng)”的關(guān)系。
如何根據(jù)系統(tǒng)的特點(diǎn)簡(jiǎn)化建模,是大學(xué)生應(yīng)該掌握的重要能力。目前教學(xué)中給學(xué)生的往往是簡(jiǎn)化好的模型,學(xué)生缺少建模的訓(xùn)練。對(duì)學(xué)生來(lái)說,很多知識(shí)在考試后不久就可能忘了,而處理問題的方法和思維模式則可以一輩子不忘記。因此在教學(xué)中少量介紹一些建模的方法,讓學(xué)生自行對(duì)一些有趣的現(xiàn)象進(jìn)行簡(jiǎn)化分析,將使學(xué)生受益終身。