王艷海 魏江南 孟義昌 張玉 馬乾偉 擺永強

摘 要：為了解決交通運輸和工業(yè)行業(yè)中移動重物的起重機動臂和支撐點不穩(wěn)定的擺動幅度問題，采用一個懸掛點在水平面內(nèi)做圓周運動的擺球系統(tǒng)模型來模擬這一問題的動力學(xué)過程。應(yīng)用拉格朗日方程推導(dǎo)了擺球運動的微分方程，采用四階龍格-庫塔算法進行了數(shù)值求解，得到了給定懸掛點圓周運動半徑、頻率和擺長的條件下，擺球做穩(wěn)定圓周運動的擺角，計算了擺球的穩(wěn)定運動軌跡，并將計算結(jié)果與實驗進行了對比，理論計算和實驗結(jié)果精確吻合。研究內(nèi)容在科研和工程技術(shù)領(lǐng)域有著廣泛的應(yīng)用，研究結(jié)果對起重機的機械吊臂設(shè)計有重要的參考價值。同時，文中所用到的理論建模、計算機編程、數(shù)值計算、實驗設(shè)計等也可以作為一個非常好的教學(xué)案例，培養(yǎng)學(xué)生應(yīng)用所學(xué)理論解決實際問題的能力。

關(guān)鍵詞：理論力學(xué)；圓周運動；拉格朗日方程；四階龍格-庫塔算法

中圖分類號：O311.3 文獻標(biāo)志碼：A

文章編號：1008-1542（2018）03-0238-05

在一個力學(xué)擺的系統(tǒng)中，如果擺球的懸掛點在水平面內(nèi)做圓周運動，在某些給定的初始條件下，擺球?qū)霭霃礁〉姆€(wěn)定的圓周運動，這一研究題目來自第二十九屆國際青年物理學(xué)家錦標(biāo)賽IYPT第二題：滯后的擺。對該問題的研究，無論在科研、工程技術(shù)領(lǐng)域和教學(xué)中都有著廣泛的應(yīng)用，有助于對起重機吊臂的設(shè)計[1-3]，也可以作為理論力學(xué)教學(xué)的補充，增強學(xué)生對理論力學(xué)相關(guān)理論知識的理解，提高學(xué)生將物理理論應(yīng)用到實際問題中去的能力[4-6]。

1 理論模型

一根長為R的硬質(zhì)桿，桿的一端固定在o點，另一端s繞o點作水平圓周運動，并在硬質(zhì)桿的s端連接一根長度為l的細繩，細繩的下端懸掛一個金屬擺球，繩子的伸縮和質(zhì)量不計。s端開始運動后，小球開始擺動，如圖1所示（擺球的三維坐標(biāo)系，xoy平面為水平面）。給定不同的初始條件，小球會呈現(xiàn)出不同的運動軌跡，研究擺球做半徑小于R的穩(wěn)定的圓周運動的條件。該系統(tǒng)所受的約束是完整約束，可以應(yīng)用拉格朗日方程進行求解[1-2，7-14]。

3 實驗驗證

實驗器材包括：交流變壓器（0～300 V，可調(diào)），風(fēng)扇（55 W），單擺小球，轉(zhuǎn)速測量儀，米尺，視頻拍攝儀，細繩。

實驗過程如下：在風(fēng)扇葉片上距離轉(zhuǎn)軸15 cm處鉆一個小孔，把半徑為1 cm的擺球用細繩系在風(fēng)扇小孔處，把風(fēng)扇轉(zhuǎn)動面放置為水平并固定好，調(diào)節(jié)交流變壓器電壓控制風(fēng)扇轉(zhuǎn)速，用轉(zhuǎn)速測速儀測量風(fēng)扇的實際轉(zhuǎn)速，用視頻拍攝儀拍攝小球的運動軌跡。根據(jù)數(shù)值模擬部分設(shè)定的初始條件進行實驗，采用運動軌跡分析軟件Tracker恢復(fù)擺球的軌跡，將實驗結(jié)果與計算結(jié)果進行對比，如圖5所示。實線為實驗數(shù)據(jù)，擺球軌跡非常接近于圓（離心率很小的橢圓）。偏離圓軌跡的原因有二：其一是由于在實驗中初始條件非常難以控制，很難達到其理想值；其二是拍攝的角度稍微偏離了豎直方向（這是為了避開風(fēng)扇葉片對擺球的遮擋），導(dǎo)致用Tracker軟件恢復(fù)擺球軌跡時帶來誤差。圖5中虛線為初始條件θ=11.6°，其他值與數(shù)值模擬計算部分相同時的擺球的理論計算軌跡，從圖中可以看出數(shù)值模擬計算的軌跡與實驗軌跡符合的很好。擺球做圓周運動時擺線與豎直方向平均夾角的實驗值為11.02°，與理論值的平均夾角（11.6°）的相對誤差約為6.2%，誤差很小，理論計算與實驗結(jié)果吻合的很好。懸掛點和擺球的運動步調(diào)如圖6所示，圖中實線和虛線分別表示懸掛點和擺球的運動步調(diào)。兩點的運動步調(diào)剛好相反，說明擺球的運動步調(diào)滯后于懸掛點，滯后的相位為π，實驗中也確實如此。

4 結(jié) 語

本文研究了懸掛點在水平面內(nèi)做勻速圓周運動時擺球的運動規(guī)律，應(yīng)用保守力系的拉格朗日方程推導(dǎo)了擺球運動的微分方程，并采用python語言進行數(shù)值模擬計算，理論計算結(jié)果與實驗結(jié)果精確吻合。研究表明當(dāng)滿足一定的初始條件時，小球也將在水平面內(nèi)做穩(wěn)定的半徑更小的勻速圓周運動，而且相位滯后懸掛點為π，該運動對初值非常敏感，如果偏離特定的初始條件，擺球的運動軌跡將變得混亂，其緊鄰的軌道是多個自由度簡諧振動的疊加。研究結(jié)果可以作為理論力學(xué)教學(xué)的一個極佳的案例，將物理學(xué)原理應(yīng)用到具體的實際問題中去，涉及到理論推導(dǎo)、計算機編程模擬計算和實驗驗證，將單純枯燥的理論與計算和實驗相結(jié)合，有助于激發(fā)學(xué)生的學(xué)習(xí)興趣，提高學(xué)生的學(xué)術(shù)研究能力。同時，本文的理論結(jié)果和數(shù)值模擬程序也可以應(yīng)用到起重機的吊臂設(shè)計等相關(guān)機械設(shè)計的實際問題當(dāng)中[18-20]。

參考文獻/References：

[1] GHIGLIAZZA R M， HOLMES P. On the dynamics of cranes， or spherical pendula with moving supports[J]. International Journal of Non-Linear Mechanics， 2002， 37（7）：1211-1221.

[2] NELSON R A， OLSSON M G. The pendulum：Rich physics from a simple system[J]. American Journal of Physics， 1986， 54（2）：112-121.

[3] TRUEBA J L， BALTANS J P， SANJUN M A F. A generalized perturbed pendulum[J]. Chaos Solitons & Fractals， 2003， 15（5）：911-924.

[4] 周衍柏. 理論力學(xué)教程[M]. 北京：高等教育出版社， 1986.

[5] 張偶利， 胡其圖， 張小靈，等. 傅科擺運動軌跡的計算機動態(tài)模擬及其教學(xué)應(yīng)用[J]. 物理與工程， 2006， 16（2）：36-40.

ZHANG Ouli， HU Qitu， ZHANG Xiaoling， et al. The dynamic computer simulation of the motion track of foucault pendulum and its application in teaching[J]. Physics and Engineering， 2006， 16（2）：36-40.

[6] 劉樹青， 吳洪濤. 一種新型起重機器人的運動學(xué)分析與應(yīng)用[J]. 河北科技大學(xué)學(xué)報， 2004， 25（2）：58-61.

LIU Shuqing， WU Hongtao. Kinematic analysis and application of a new style crane robot[J]. Journal of Hebei University of Science and Technology， 2004， 25（2）：58-61.

[7] PAULA A S D， SAVI M A. Chaos and transient chaos in an experimental nonlinear pendulum[J]. Journal of Sound & Vibration， 2006， 294（3）：585-595.

[8] DOOREN R V. Chaos in a pendulum with forced horizontal support motion： A tutorial[J]. Chaos Solitons & Fractals， 1996， 7（1）：77-90.

[9] HORTON B， SIEBER J， THOMPSON J M T， et al. Dynamics of the nearly parametric pendulum[J]. International Journal of Non-Linear Mechanics， 2011， 46（2）：436-442.

[10]BELYAKOV A O. On rotational solutions for elliptically excited pendulum[J]. Physics Letters A， 2011， 375（25）：2524-2530.

[11]KHOLOSTOVA O V. Some problems of the motion of a pendulum with horizontal vibrations of the point of suspension[J]. Journal of Applied Mathematics & Mechanics， 1995， 59（4）：553-561.

[12]BUTIKOV E I. Subharmonic resonances of the parametrically driven pendulum[J]. Journal of Physics A Mathematical General， 2002， 35（30）：6209-6231.

[13]LAHIRI A， SINHA R M. The Hamiltonian Hopf bifurcation： An elementary perturbative approach[J]. International Journal of Non-Linear Mechanics， 2001， 36（5）：787-802.

[14]JR W T G， SCHCK M. Simulations of nonlinear pivot-driven pendula[J]. American Journal of Physics， 1997， 65（5）：376-381.

[15]LANDAU R H， P M J， BORDEIANU C C. Computational physics： Problem solving with python[J]. John Wiley & Sons， 2015， 51（36）： 12039-12048.

[16]WILLIAM H. C++數(shù)值算法[M]. 2版. 北京：電子工業(yè)出版社， 2005.

[17]LUTZ M， 李軍， 劉紅偉. Python 學(xué)習(xí)手冊[M]. 北京：機械工業(yè)出版社， 2011.

[18]王雷， 劉志虎， 李震，等. 六自由度裝配機械手設(shè)計及運動學(xué)求解[J]. 河北科技大學(xué)學(xué)報， 2014， 35（5）：417-427.

WANG Lei， LIU Zhihu， LI Zhen， et al. Design and kinematics solution of 6-DOF assembly manipulator[J]. Journal of Hebei University of Science and Technology， 2014， 35（5）： 417-427.

[19]吳學(xué)禮， 劉浩南， 許晴. 機器人手臂控制系統(tǒng)的設(shè)計與研究[J]. 河北科技大學(xué)學(xué)報， 2014， 35（4）：361-365.

WU Xueli， LIU Haonan， XU Qing. Design and research of robot arm control system[J]. Journal of Hebei University of Science and Technology， 2014， 35（4）：361-365.

[20]梁建術(shù)， 蔡建軍. 一類含間隙滯回動力系統(tǒng)周期運動和分岔[J]. 機械強度， 2007， 29（3）：390-393.

LIANG Jianshu， CAI Jianjun. Periodic motion and bifurcation of a class of hysteretic dynamics with gaps[J]. Journal of Mechanical Strength， 2007， 29（3）： 390-393.