王洪濤 李艷

(中國礦業(yè)大學(xué)物理學(xué)院 江蘇 徐州 221116) (2018-04-04)

大學(xué)物理課程是本科教育階段的一門重要的公共基礎(chǔ)課程，課程具有知識量大、理論性強、結(jié)論抽象等特點．電磁學(xué)部分由于理論抽象、計算過程需要用到相對復(fù)雜的微積分知識等特征，致使學(xué)生在理解掌握知識上存在一定的困難．實際教學(xué)過程中經(jīng)常性會遇到部分學(xué)生對傳統(tǒng)的定理推導(dǎo)過程抱有疑問，他們更傾向于采用“笨一點”的方法來思考問題．比如電場的高斯定理，該定理是計算規(guī)則對稱分布電場的有力工具，但是總有部分學(xué)生提出能不能利用前面學(xué)過的點電荷的電場強度疊加方法得到均勻帶電球體或球面在周圍任意一點的場強？學(xué)生的疑問其實暴露出了一定的問題，即學(xué)習(xí)知識的慣性思維，每當(dāng)學(xué)到新知識的時候往往還希望嘗試用已熟練掌握的知識去解釋、推導(dǎo)及驗證．基于此，在講授新知識的時候可以根據(jù)教學(xué)進(jìn)度分層次地引導(dǎo)學(xué)生進(jìn)行驗證探索，這對學(xué)生深入理解掌握物理知識、復(fù)習(xí)應(yīng)用數(shù)學(xué)及計算機知識都有很大的益處．

本文以靜電場中的高斯定理為例，采用“笨方法”——點電荷場強的積分方法和“小高招”——數(shù)值編程模擬方法對高斯定理進(jìn)行驗證，使學(xué)生從不同角度不同層次上理解掌握知識．

1 高斯定理計算場強

采用高斯定理[1]計算均勻?qū)ΨQ分布電場的電場強度是大學(xué)物理中計算場強的重要方法之一，均勻帶電的圓球面在球面[2]及其周圍激發(fā)的電場具有均勻?qū)ΨQ分布的特點，是高斯定理應(yīng)用的常見情況．設(shè)球面半徑為R，電荷密度為σ，寫出高斯定理的具體表達(dá)式為[3]

(1)

由此易解出距球心為x處的一點的電場強度為

(2)

2 矢量疊加法計算場強

通過點電荷電場的矢量疊加方法理論上可以直接積分求解任意帶電體周圍激發(fā)的電場．對于具有圓對稱性的帶電球面，可以通過合理劃分微分單元把一般的三維體積分簡化為一維線積分來處理．均勻帶電圓環(huán)在軸線上距圓心x的點的電場強度為

(3)

其中q為圓環(huán)所帶電荷，R為圓環(huán)半徑，x為場點距帶電圓環(huán)圓心的距離．可以把帶電球面沿垂直x軸方向分割為若干個半徑連續(xù)變化的同軸圓環(huán)，如圖1所示．以球心為坐標(biāo)原點，水平向右為x軸方向建立坐標(biāo)系．

圖1 帶電球面微分單元劃分示意圖

用夾角α表示微分元圓環(huán)的位置，則圓環(huán)所帶電荷為

dq=σRdα2πRsinα

其在x點的場強為

通過對α積分可以得到帶電球面在距球心x處的電場強度表達(dá)式

(4)

對式(4)進(jìn)行化簡，令β=x2+R2-2xRcosα則可得到

(5)

其中

β1=(x-R)2

β1=(x+R)2

通過“笨一點”的積分方法得到結(jié)果與高斯定理得到的結(jié)果完全一致，不僅驗證了高斯定理的正確性打消了學(xué)生的疑慮，使學(xué)生對微積分方法在大學(xué)物理學(xué)習(xí)中的靈活運用有了更深入的理解和體會，同時也說明點電荷電場強度的疊加原理對任意帶電體的電場計算都是適用的，區(qū)別只在于積分計算的難易．

3 數(shù)值編程模擬

物理學(xué)知識理論性較強，結(jié)論定理大都比較抽象．可以通過計算機編程把抽象的公式定理用直觀的圖線演示出來，增強學(xué)生對知識的掌握[4]．積分的數(shù)值編程方法有很多，學(xué)生在數(shù)學(xué)及計算機課程上都有所涉及，通過計算機編程模擬不僅可以把抽象的知識形象化，而且能夠借此讓其了解數(shù)學(xué)及計算機編程手段的應(yīng)用，可謂一舉多得．編程語言有很多，常見的有Matlab，C++，F(xiàn)ortran，Python等等，雖然在代碼書寫上存在一些差異，但是算法是共通的．Python 是一種不受局限、跨平臺的開源編程語言，它功能強大且簡單易學(xué)，代碼簡潔易懂，目前在很多行業(yè)已得到廣泛應(yīng)用．?dāng)?shù)值編程模擬的主要目的是驗證定理結(jié)論或?qū)Τ橄笾R進(jìn)行形象化展示，因此對運算速度及精度要求不高，可以采用經(jīng)典的梯形法計算積分[5]， Python程序積分代碼如下：

#------------------

import math

sigma = 1e-8

epsilon = 8.854187817e-12

R = 0.2; x = 0.4; a = -1; b = 1; h = 2.0/16000; s = 0; m = a; f0 = 0

for i in range(16000):

m += h

f1 = (sigma*R**2/2/epsilon)*(L-R*m)/(L**2+R**2-2*L*R*m)**1.5

s += h*(f0+f1)/2

f0 = f1

print(x,s)

#-------------------

為了考察積分區(qū)間[-1,1]劃分子區(qū)間的數(shù)目與計算精度的關(guān)系，程序連續(xù)計算了子區(qū)間數(shù)目分別為1 000, 2 000，…，16 000共16組數(shù)據(jù)，并把計算結(jié)果與公式(2)進(jìn)行了對比，求出其相對誤差，如圖2所示，其中橫軸為等分區(qū)間數(shù)目N，縱軸為模擬值與理論值的相對誤差．由圖可以看出，隨著等分區(qū)間數(shù)目的增加，計算結(jié)果很快收斂于理論值．當(dāng)?shù)确謪^(qū)間數(shù)目為N=7 000時，相對誤差已降至5×10-5以下．梯形積分法編程簡單，計算精度較高，能夠滿足教學(xué)驗證及學(xué)生課下實踐的要求，在增加對物理知識的理解及對抽象知識的直觀展示方面具有積極意義，適合在大學(xué)物理教學(xué)過程中推廣應(yīng)用．

圖2 相對誤差與積分區(qū)間等分?jǐn)?shù)目的關(guān)系

取N=16 000，R=0.2 m固定不變，場點距球心坐標(biāo)由近及遠(yuǎn)連續(xù)變化，x=1.2R，1.4R，…，4.0R，通過積分程序分別計算得到對應(yīng)x距離下的電場強度，把x及對應(yīng)電場強度數(shù)據(jù)存入數(shù)據(jù)文件“E.data”中，然后利用Gnuplot進(jìn)行繪圖．Gnuplot是免費開源的科學(xué)繪圖軟件，具有很強的數(shù)據(jù)繪圖及參數(shù)擬合能力[6]，同時支持交互和腳本兩種工作模式，在輸出圖片的高質(zhì)量和高可控性上有較強的優(yōu)勢．在Gnuplot中寫入下述命令：

sigma=1e-8; R=0.2; epsilon=

8.854 187 817e-12

f(x)=sigma*R**2/(epsilon*x**2)

plot f(x) with line, ‘E.data’ using 1:2 with point

圖3 理論結(jié)果與數(shù)值模擬結(jié)果對比圖

4 小結(jié)

本文針對電場高斯定理教學(xué)實踐中的問題，通過點電荷場強積分方法計算了均勻帶電球面周圍的電場分布，并且采用開源、高效的程序語言Python及科學(xué)繪圖軟件Gnuplot對電場分布進(jìn)行了數(shù)值模擬．采用多管齊下的方法使學(xué)生從不同角度對高斯定理展開了全面深入的學(xué)習(xí)和理解，不僅促進(jìn)了學(xué)生對物理知識的掌握，而且讓他們體會到了微積分及數(shù)值編程知識在物理學(xué)習(xí)過程中的具體應(yīng)用．采用Python及Gnuplot開源軟件實現(xiàn)抽象物理知識的直觀形象化在教師教學(xué)實踐及學(xué)生課下探索過程中具有很高的可執(zhí)行性，值得在大學(xué)物理教學(xué)中進(jìn)行推廣和應(yīng)用．