邱智宇 肖桂娜 蘇桂鋒

(上海師范大學物理系 上海 200234)

轉動慣量是剛體繞軸轉動時慣性量度的重要特性，在航空事業(yè)、風力發(fā)電、汽車零件等領域中起重要作用[1]．剛體的轉動慣量不僅取決于構件總質量，還與其形狀、質量分布以及轉軸位置有關．無論剛體形狀規(guī)則與否，利用扭擺法測量剛體轉動慣量是精度較高的方法，是大學物理的重要必修實驗，更是被廣泛地應用于測量制造業(yè)零件的轉動慣量[1,2]．由此，教師應重視剛體轉動慣量測量的實驗教學，充分展現現代信息技術時代的特征，合理利用多媒體技術實現教育現代化，調動學生對物理模型的感性認知，喚醒物理觀念的理性思維；信息、編程技術與物理實驗教學相結合以實現教育信息化，已經成為高校物理教育改革目前的熱門研究課題和教育發(fā)展的優(yōu)良趨勢[5]．

本文以質量均勻物體對其幾何中心轉軸轉動慣量的測量為例，利用MATLAB編程與Curve Fitting Tool工具，不僅能夠實現轉動慣量理論值可視化、實驗數據擬合精確化，更能將實驗過程以動畫形式呈現．

1 實驗工具與實驗原理

1.1 實驗原理

將待測物體置于螺旋彈簧載物臺上，物體與螺旋彈簧同步運動．根據胡克定律，當物體在水平面上轉動，轉過某一角度θ后，螺旋彈簧產生的恢復力矩M為

M=-Kθ

(1)

式中，K為螺旋彈簧的扭轉系數．再根據轉動定理，物體的力矩M為

M=Iα

(2)

式中，I表示物體繞轉軸的轉動慣量，α為角加速度．物體在螺旋彈簧恢復力矩的作用下與彈簧共同繞垂直于水平面的轉軸做往復運動，因此式(1)、(2)聯立得

(3)

忽略軸承的摩擦力矩，由式(3)和角加速度的定義可得

(4)

令

可得關于θ的二階常微方程

(5)

由此可知置于扭擺上物體的運動具有簡諧振動的特性，此方程的解為

θ=Acos(ωt+φ)

(6)

式中A為簡諧振動的振幅，φ為相位角，皆由初始條件決定；ω為角速度，扭擺法過程的簡諧振動周期為

(7)

所以待測剛體的轉動慣量為

(8)

圖1 扭擺法測量物體轉動慣量實驗原理圖

(9)

式中由于T0和T1能夠通過實驗測量

是已知量，則可得到彈簧的扭轉系數K的值為

(10)

由此，只要測得待測物體隨著螺旋彈簧轉動的總周期，即可得到待測物體和支架的轉動慣量之和[3,4]．綜上所述，待測物體的轉動慣量為

(11)

1.2 MATLAB介紹

MATLAB被應用于算法、數據分析與計算、數據可視化的高級計算語言和交互式環(huán)境，主要包括MATLAB和Simulink兩大部分．MATLAB的主界面如圖2所示.

圖2 MATLAB軟件主界面

圖2中編輯器為主體，功能為程序編入與調試，能以顏色區(qū)別代碼顯示的功能．運行編輯器，工作區(qū)以矩陣形式存儲數據，并得出數據計算結果和繪制相關函數圖像或動態(tài)圖像．命令行窗口(Command Window)為主要交互式窗口，位于編輯器下方，輸入指令后顯示除圖形以外的執(zhí)行結果．命令窗口中“?”符號為命令提示符，表示MATLAB正在處于準備狀態(tài)．在命令提示符后輸入命令并按回車鍵，MATLAB就會解釋執(zhí)行所輸入的命令，并在命令后面給出計算結果．

MATLAB語言具有解釋性，程序簡潔易懂．不僅能夠直接在命令行窗口輸入命令，也可以在編輯器內編寫應用程序，然后在MATLAB環(huán)境下對程序命令進行處理，處理的對象為矩陣或調用函數，最后返回運算結果．因此，其程序編寫較易上手，可以通過直觀的數學表達式實現對數據的處理．在扭擺法測量剛體轉動慣量的實驗中需要運用其主體與工具包Curve Fitting，前者用于實驗過程動畫繪制與剛體轉動慣量理論值函數圖像繪制；后者于工具欄中APP一欄，能直接讀取實驗實測中時間與位移的數據，通過傅里葉函數分析擬合，得函數表達式，可直接讀出其周期值．在扭擺法測量剛體轉動慣量的實驗教學中， MATLAB以其語言易懂化、圖像可視化、數值分析能力優(yōu)異性等特點，不論是實驗過程、實驗理論值，還是實驗測量值的函數擬合，均起到極佳的教學輔助作用．

2 實驗方法

根據實驗原理需要測定各類規(guī)則物體轉動慣量的理論值和扭轉常數K，根據規(guī)則物體轉動慣量的積分計算過程，通過編輯MATLAB程序，可獲得如圖3所示的不同厚度的均勻圓筒和均勻球殼、不同長度均勻細桿繞中心軸或轉軸，以及其他形狀規(guī)則剛體的轉動慣量理論值[5]．點擊“工具”欄中的數據提示，便可選取理論曲線上任意一點獲得單位厚度或距離細桿軸線一定距離的物體的轉動慣量數值,如圖3中均勻圓筒對柱體軸線的轉動慣量圖像中，當圓筒中空心部分半徑R0與圓筒總體半徑R之比為某值時，圖3(a)中可直接讀得數據與圓柱體的質量和半徑平方的乘積即為已知質量和幾何參數圓柱體的轉動慣量理論值．理論值的獲得不僅有助于分析實際測量物理量數據的精確程度，而且在扭擺法實驗數據處理過程中具有輔助作用．

圖3 規(guī)則物體轉動慣量的獲得

獲得轉動慣量的理論數據和相關實驗準備工作后，利用扭擺法開始進行實物測量，第一步測量金屬托盤空載時的周期．其具體步驟為：先調整扭擺底座的底角螺絲，使螺旋彈簧上側水準儀中的氣泡居中．然后，裝配金屬載物盤，并調整光電探頭的位置，測定其擺動周期T0，測3次求平均值．扭擺法測量轉動慣量的實驗過程步驟存在先后順序，且口頭描述較為復雜，實驗教學的過程中可通過MATLAB進行動畫展示，分析每一步驟的細節(jié)與注意點．如圖4所示，圖中灰色圓框為金屬底座的俯視圖像，突出部分為擋光軸，黑色線條組成部分為光電門．

圖4 扭擺法測量剛體轉動慣量實驗過程俯視動畫關鍵幀圖像

開始時，如圖4(a)所示，需要將擋光軸與光電門所在軸線垂直，保證初始擺角相同且約為90°，此時螺旋彈簧對于轉動慣量的實驗測量誤差影響最小[6]．釋放螺旋彈簧后，載物盤隨螺旋彈簧旋轉(如圖4所示為先順時針旋轉)，通過光電門之后，周期測定儀會記錄讀數，擋光軸繼續(xù)向順時針方向旋轉，到達彈性限度后再沿著逆時針方向旋轉，以此進行往復運動，最終周期測定儀記錄擋光軸往復10次的時間，該時間除以周期測量次數后，便可得載物盤轉動一次的運動周期，減少扭擺周期測量的誤差．利用動畫演示，可以將初始設置情況與周期這類微小數據測量方法描述得更加細致、準確．通過Mxcel的編程功能，可以再現螺旋彈簧運動情況，由此能夠讓學生更好地理解實驗過程各參數的測量步驟及注意事項．將具象的動畫過程與單一的語言相結合，提升了物理實驗課堂教學的多樣性與趣味性．

此外，實驗中需要捕捉大量的周期數據，按照傳統(tǒng)方法需要花費大量精力進行數值整理與計算．利用Tracker軟件輔助自動追蹤擋光軸，捕獲擋光軸在10個周期內的運動軌跡，可將捕捉的數據軌跡以Excel的文件形式直接導入MATLAB軟件，在命令行窗口轉化為矩陣的形式．通過工具包Curve Fitting中的傅里葉函數分析擬合，不僅能夠直觀地觀察其運動軌跡，還可以得到函數表達式，從而可直接讀出載物臺或其他物體運動的周期值．圖5中實線框處下拉，選擇“Fourier”軟件將自動進行數據擬合處理，虛線框中的數據即為數據處理完畢后，載物盤隨螺旋彈簧運動過程中的角頻率，用2π除以該數值即為扭擺法測量剛體轉動慣量實驗中需要的周期．

圖5 Curve Fitting傅里葉函數擬合數據界面

3 數據處理

本次實驗在螺旋彈簧扭轉系數的確定中選用質量為0.877 kg，直徑為0.089 5 m的圓柱體，所得螺旋彈簧的扭轉系數K如表1所示.

表1 螺旋彈簧扭轉系數K的確定

本文提及實驗中彈簧的扭轉系數K=0.033 1 N/m，已知支架的轉動慣量為I桿支架=3.21×10-3kg·m2，I球底座=1.87×10-3kg·m2，本次實驗測量的對象有：大塑料圓柱體、小塑料圓柱體、塑料球、金屬圓筒和金屬細桿，表2為在MATLAB軟件輔助下利用扭擺法測量上述物體轉動慣量的實驗數據與誤差.

表2 各種物體轉動慣量的測定

利用Curve Fitting工具包中的傅里葉函數分析擬合直接讀取周期數據與傳統(tǒng)實驗教學相比：轉動周期的獲得操作便捷且代入公式得到的轉動慣量值相對誤差較小，實驗數據較可靠，具有參考價值．由實驗數據可見，被測樣品的實驗誤差均低于或在1.5%，金屬物體測得的轉動慣量較理論值的相對誤差較大，可能由于其內部結構不均勻，可成為后續(xù)轉動慣量實驗改進的突破口．

4 結論

本文將MATLAB編程與Curve Fitting工具包應用于扭擺法測量剛體轉動慣量的大學普通物理實驗教學及其數據處理中，不僅能夠科學嚴謹地得到規(guī)則剛體轉動慣量的理論值，也能夠通過傅里葉函數擬合迅速便捷地得出物體轉動慣量的實驗值，在掌握物體轉動慣量實驗原理的基礎上，摒棄了傳統(tǒng)實驗數據冗長的計算步驟，同時實驗數據的相對誤差較??；而對于課堂教學方面，通過MATLAB編程實現圖像的運動與旋轉，能再現扭擺法測量的過程，便于學生的想象與理解，利用多媒體設備實現課堂教學的多元化、可視化和信息化，融教育于信息時代的社會大背景，實現教育形式的進步，貫徹教育信息化的落實．