苗正 牛放 花少震

DOI：10.16661/j.cnki.1672-3791.2105-5042-6439

摘? 要：為了提高等效法求非均質(zhì)或者不規(guī)則幾何體轉(zhuǎn)動慣量的精度，該文以非均質(zhì)發(fā)動機搖臂為研究對象，研究了拉格朗日插值法、牛頓插值法、拋物線插值法在根據(jù)等效原理求非均質(zhì)發(fā)動機搖臂轉(zhuǎn)動慣量的精度與截斷誤差。通過對比分析發(fā)現(xiàn)拋物線插值法的精度最低，拉格朗日插值法與牛頓插值法的截斷誤差較小。綜合精度和誤差分析可以得出，牛頓插值法在等效法求非均質(zhì)發(fā)動機搖臂轉(zhuǎn)動慣量上具有較高精度，該研究可以為等效法求非均質(zhì)或者不規(guī)則幾何體的轉(zhuǎn)動慣量提供參考。

關(guān)鍵詞：轉(zhuǎn)動慣量? ?等效原理? ?插值方法? ?截斷誤差

中圖分類號：O302? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻標(biāo)識碼：A文章編號：1672-3791（2021）08（a）-0164-04

Study on Interpolation Method of Moment of Inertia Based on Equivalent Principle

MIAO Zheng? ?NIU Fang? ?HUA Shaozhen*

（School of Mechanical Engineering， Henan Institute of Technology， Xinxiang， Henan Province， 453003 China）

Abstract： In order to improve the accuracy of calculating the moment of inertia of heterogeneous or irregular geometry based on the equivalent method， heterogeneous engine rocker arm has been taken as the research object， then the Lagrange interpolation method， Newton interpolation method and parabola interpolation method were studied based on the equivalent principle to calculate the accuracy and truncation error of the moment of inertia of. Through comparative analysis， it is found that the accuracy of parabola interpolation method is the lowest， and the truncation error of Lagrange interpolation method and Newton interpolation method is smaller. Comprehensive accuracy and error analysis can be concluded that Newton interpolation method has high accuracy in the equivalent method to calculate the moment of inertia of heterogeneous engine rocker arm. This study can provide a reference for the equivalent method to calculate the moment of inertia of heterogeneous or irregular geometry.

Key Words： Moment of inertia; Equivalent principle; Interpolation method; Truncation error

轉(zhuǎn)動慣量是剛體繞定軸轉(zhuǎn)動時慣性的度量[1]，是航空航天、儀器儀表、汽車、電力等行業(yè)中重要的參量。由于轉(zhuǎn)動慣量只和轉(zhuǎn)軸的位置、剛體的質(zhì)量分布以及港田的形狀有關(guān)，所以對于一些非均質(zhì)的幾何體或者形狀復(fù)雜的幾何體，一般通過實驗手段，基于等效原理，采用插值手段測定。吉丹霞等人[2]利用等效法研究了變速箱中不同輪系等效轉(zhuǎn)動慣量的計算方法;仝哲旭等人[3]通過設(shè)計非接觸平面渦卷彈簧，測量了正弦扭矩標(biāo)準(zhǔn)裝置氣浮軸系的轉(zhuǎn)動慣量;趙亞燊等人[4]利用等效法計算了窄輪距拖拉機的轉(zhuǎn)動慣量。周艷明等人[5]研究了三線擺實驗中大擺角對轉(zhuǎn)動慣量的影響，徐錢欣等人[6]、陳坤波等人[7]、陳慶東等人[8]利用改進的三線擺降低了實驗誤差。目前，已有研究主要側(cè)重于等效法的測量應(yīng)用以及對三線擺實驗裝置的改進，對等效原理中用插值法的研究較少。該文以非均質(zhì)發(fā)動機搖臂為研究對象，采用三線擺實驗，根據(jù)等效原理，研究不同插值方法對等效原理的適用性，通過精度及誤差分析，優(yōu)化對插值方法的選擇。

1? 實驗裝置及實驗原理

1.1 實驗裝置

該研究所采用的實驗裝置為ZME-Ⅰ型實驗臺，主要由柜體、風(fēng)振模型支架、三線擺、不銹鋼頂罩等部分組成，如圖1（a）所示。圖1（b）為實驗中所采用均質(zhì)圓柱體和非均質(zhì)發(fā)動機搖臂。運用此設(shè)備的三線擺可以測得規(guī)則和不規(guī)則物體的轉(zhuǎn)動周期，并可以調(diào)節(jié)三線擺的吊線長度，從而得到不同吊線長度下物體的轉(zhuǎn)動周期。

三線擺吊線線長選用60 cm，等重圓柱體直徑d=18 cm，高h=20.7 mm，圓柱體密度ρ=7.8 g/cm3。分別取兩圓柱體中心距s為30、40、50、60 mm。

1.2 實驗原理

該文通過使用三線擺測量圓盤的轉(zhuǎn)動慣量方法，使用了等效法來測量非均質(zhì)發(fā)動機搖臂的轉(zhuǎn)動慣量。在實驗中同時放置兩個同樣三線擺，使它們下降到相同的位置，一個三線擺放置非均質(zhì)發(fā)動機搖臂，另一個三線擺放置兩個相同的均質(zhì)圓柱，兩個均質(zhì)圓柱的重量與非均質(zhì)發(fā)動機搖臂的重量一致，通過不斷調(diào)節(jié)兩個等效均質(zhì)圓柱體之間的中心距s，向放置有非均質(zhì)發(fā)動機搖臂的擺動周期逼近。當(dāng)兩個三線擺系統(tǒng)的擺動周期相同時，則可以大致認為均質(zhì)圓柱與發(fā)動機搖臂的轉(zhuǎn)動慣量是一致的，而放置有均質(zhì)圓柱體的三線擺系統(tǒng)，可以計算得出非均質(zhì)發(fā)動機搖臂的等效轉(zhuǎn)動慣量。

1.3 實驗方法

在已知等重圓柱體的物理參數(shù)后，兩圓柱對中心軸O的轉(zhuǎn)動慣量的理論計算公式為：

式中，s為兩均質(zhì)圓柱體的中心距，具體情況見圖2所示。

通過多次實驗測得轉(zhuǎn)動周期，并通過式（1）計算轉(zhuǎn)動慣量的理論值，如表1所示。

為了使測量的非均質(zhì)發(fā)動機搖臂的轉(zhuǎn)動周期精度較高，多次實驗后刨除明顯有差異的，選取5組實驗，測試5次搖臂轉(zhuǎn)動10次的時間，得到搖臂扭動周期，并取平均值得到平均周期。實驗測量數(shù)據(jù)如表2所示。

2? 實驗結(jié)果與分析

該文將運用牛頓插值法、拉格朗日插值法、拋物線插值法估算非均質(zhì)發(fā)動機搖臂的轉(zhuǎn)動慣量。

2.1 牛頓插值法

2.1.1 使用牛頓插值法對非均質(zhì)發(fā)動機搖臂轉(zhuǎn)動慣量進行插值

通過表1的得到的實驗數(shù)據(jù)，將數(shù)據(jù)代入得到牛頓插值多項式：

其中：

將代入式（2），得到，此時的就是由牛頓插值法估計的非均質(zhì)發(fā)動機搖臂的等效轉(zhuǎn)動慣量。

即非均質(zhì)發(fā)動機搖臂的轉(zhuǎn)動慣量。

2.1.2 牛頓插值法的誤差分析

使用牛頓插值法分析轉(zhuǎn)動慣量，不在插值節(jié)點上的函數(shù)值會有截斷誤差，牛頓插值法的插值余項為：

其中，，得到：

即牛頓插值法計算轉(zhuǎn)動慣量的截斷誤差為：

2.2 拉格朗日插值法

2.2.1 使用拉格朗日插值法對非均質(zhì)發(fā)動機搖臂轉(zhuǎn)動慣量進行插值

通過表1得到的實驗數(shù)據(jù)，得到拉格朗日插值多項式：

式（3）中，，，，。

將代入式（3），得到，此時就是由拉格朗日插值法擬合的非均質(zhì)發(fā)動機搖臂的等效轉(zhuǎn)動慣量。

即非均質(zhì)發(fā)動機搖臂的轉(zhuǎn)動慣量。

2.2.2 拉格朗日插值法截斷誤差分析

因為使用拉格朗日多項式擬合離散數(shù)據(jù)，不在插值節(jié)點上的值會有誤差，用截斷誤差表示，截斷誤差也稱為拉格朗日插值余項，而拉格朗日插值函數(shù)的截斷誤差為：

則通過拉格朗日多項式計算的截斷誤差為：

其中，，得到：

即拉格朗日插值法計算轉(zhuǎn)動慣量的截斷誤差。

2.3 拋物線插值法

2.3.1 使用拋物線插值法對非均質(zhì)發(fā)動機搖臂轉(zhuǎn)動慣量進行插值

利用表1的得到的實驗數(shù)據(jù)，使用逐次線性插值法分析轉(zhuǎn)動慣量。因為測得的非均質(zhì)發(fā)動機搖臂的平均周期為1.052s，故取與1.052最近三點作為插值節(jié)點，即：，，三點。

將3個離散數(shù)據(jù)點代入得到擬合二次多項式：

（5）

式（5）中，

將代入式（5），得到，此時的L（1.052）就是通過拋物線插值法擬合的非均質(zhì)發(fā)動機搖臂的轉(zhuǎn)動慣量。

即非均質(zhì)發(fā)動機搖臂的轉(zhuǎn)動慣量。

2.3.2 物線插值法截斷誤差分析

拋物線插值法的插值余項為：

則拋物線插值法計算的插值余項為：

其中，，于是有：

即拋物線插值法計算轉(zhuǎn)動慣量的截斷誤差≤。

通過牛頓插值法、拉格朗日插值法、拋物線插值法估算非均質(zhì)發(fā)動機搖臂的轉(zhuǎn)動慣量，并計算了使用3種方法的誤差，得出的分析，見圖3。

3? 結(jié)論

基于等效法原理測量了與非均質(zhì)發(fā)動機搖臂等重的均質(zhì)圓柱的轉(zhuǎn)動周期，為了找到一種基于等效法求轉(zhuǎn)動慣量方面精度較高、計算方便的插值法，分別使用拉格朗日插值法、牛頓插值法、拋物線插值法對非均質(zhì)發(fā)動機搖臂的轉(zhuǎn)動慣量進行插值，通過文中分析，主要得出以下結(jié)論。

（1）牛頓插值法與拉格朗日插值法的插值余項相同，這兩種插值方法的截斷誤差為≤。

（2）拋物線插值法得到的轉(zhuǎn)動慣量值，由于拋物線為二次多項式，得到的截斷誤差≤，誤差相對于前兩個插值方法較大。

（3）而對于截斷誤差相同的拉格朗日插值法與牛頓插值法，牛頓插值法相對于拉格朗日插值法計算更加簡便、快捷。拉格朗日插值法在新增數(shù)據(jù)點時只能再次計算拉格朗日插值多項式，而牛頓插值法則不需要。

（4）3種插值法轉(zhuǎn)動慣量進行插值計算時，使用牛頓插值法更加方便。

參考文獻

[1] 哈爾濱工業(yè)大學(xué)理論力學(xué)教研室.理論力學(xué)[M].8版.哈爾濱：高等教育出版社，2016：284.

[2] 吉丹霞，陳瀟，嚴(yán)思敏.變速箱中不同輪系等效轉(zhuǎn)動慣量的計算方法[J].汽車實用技術(shù)，2021，46（6）：93-95.

[3] 仝哲旭，彭軍.轉(zhuǎn)動慣量在線測量用渦卷彈簧的設(shè)計[J].計測技術(shù)，2021，41（1）：53-56.

[4] 趙亞燊，陳俊杰，劉江輝，等.窄輪距拖拉機轉(zhuǎn)動慣量與質(zhì)心高度試驗臺研究[J].拖拉機與農(nóng)用運輸車，2021，48（1）：20-22.

[5] 周艷明，翦知漸，謝中.三線擺實驗中用大擺角測量轉(zhuǎn)動慣量[J].物理實驗，2020，40（10）：22-27，46.

[6] 徐錢欣，丁益民，范兵，等.利用智能手機測物體轉(zhuǎn)動慣量的居家DIY實驗[J].物理與工程，2020，30（6）：95-99.

[7] 陳坤波，朱瑜，鄭偉，等.測量三線擺上下懸盤間垂直距離的改進方案[J].大學(xué)物理實驗，2019，32（3）：25-29.

[8] 陳慶東，王俊平.三線擺測量剛體轉(zhuǎn)動慣量實驗的改進[J].物理通報，2016（10）：79-81.