譚　青，任志湘，魏　來(lái)

（中南大學(xué)　機(jī)電工程學(xué)院，長(zhǎng)沙410083）

球式自動(dòng)平衡裝置數(shù)值仿真系統(tǒng)設(shè)計(jì)

譚青，任志湘，魏來(lái)

（中南大學(xué)機(jī)電工程學(xué)院，長(zhǎng)沙410083）

在球式自動(dòng)平衡裝置的數(shù)值仿真中，傳統(tǒng)的仿真方法是利用Matlab進(jìn)行數(shù)值仿真，得到的結(jié)果不具直觀性和交互性。為了解決上述問(wèn)題，現(xiàn)利用VB平臺(tái)，采用龍格—庫(kù)塔法對(duì)球式自動(dòng)平衡裝置的數(shù)學(xué)模型進(jìn)行數(shù)值求解，設(shè)計(jì)一種能夠使仿真結(jié)果動(dòng)態(tài)可視化的交互式數(shù)值仿真系統(tǒng)。仿真結(jié)果表明，該數(shù)值仿真系統(tǒng)增加仿真結(jié)果的直觀性，實(shí)現(xiàn)滾球運(yùn)動(dòng)狀態(tài)的動(dòng)態(tài)呈現(xiàn)。得到的仿真結(jié)果與自同期理論一致，不平衡轉(zhuǎn)子的振幅得到有效的衰減。

振動(dòng)與波；球式自動(dòng)平衡裝置；數(shù)學(xué)模型；數(shù)值仿真系統(tǒng)；減振

對(duì)于回轉(zhuǎn)機(jī)械而言，轉(zhuǎn)子不平衡將導(dǎo)致回轉(zhuǎn)機(jī)械產(chǎn)生振動(dòng)和噪聲，甚至?xí)l(fā)災(zāi)難性破壞，由此造成的損失不容忽視。目前，解決這一問(wèn)題的一般方法是對(duì)其轉(zhuǎn)子進(jìn)行高精度的動(dòng)靜平衡［1］。由此，學(xué)術(shù)界展開了自動(dòng)平衡課題的研究。球式自動(dòng)平衡裝置是根據(jù)柔性轉(zhuǎn)子和彈性支撐的特性，不需要外部提供任何能量，只利用系統(tǒng)響應(yīng)所形成的能量來(lái)驅(qū)動(dòng)滾球的移動(dòng)和分布，從而自動(dòng)地消除轉(zhuǎn)子的不平衡。球式自動(dòng)平衡裝置對(duì)于運(yùn)動(dòng)在過(guò)臨界轉(zhuǎn)速下的旋轉(zhuǎn)機(jī)械有較好的減振效益［2］。

以往研究中，借助數(shù)值分析方法求解球式自動(dòng)平衡裝置運(yùn)動(dòng)微分方程的結(jié)果不夠直觀，并且無(wú)法直接觀察到其滾球的運(yùn)動(dòng)狀態(tài)及圓盤角速度和圓盤中心軌跡等。

根據(jù)球式自動(dòng)平衡裝置的數(shù)學(xué)模型，借助Visual Basic 6.0平臺(tái)設(shè)計(jì)數(shù)值仿真系統(tǒng)，實(shí)現(xiàn)了球式自動(dòng)平衡裝置仿真結(jié)果的動(dòng)態(tài)可視化?？紤]到在仿真過(guò)程中滾球會(huì)發(fā)生碰撞，將采用經(jīng)過(guò)碰撞調(diào)整的適用于本數(shù)學(xué)模型的4階龍格—庫(kù)塔法求解。

1　球式自動(dòng)平衡裝置的數(shù)學(xué)模型

球式自動(dòng)平衡裝置的工作原理是轉(zhuǎn)子在過(guò)臨界轉(zhuǎn)速下旋轉(zhuǎn)時(shí)，轉(zhuǎn)子內(nèi)滾球的位置將自動(dòng)的受轉(zhuǎn)子機(jī)械振動(dòng)相位的影響，移動(dòng)到偏心質(zhì)量的另一邊，一定程度上抵消偏心質(zhì)量所產(chǎn)生的離心慣性力，從而起到減振的作用［2］。

球式自動(dòng)平衡裝置的主要結(jié)構(gòu)如圖1所示，一個(gè)轉(zhuǎn)盤固定在一根轉(zhuǎn)軸上，轉(zhuǎn)盤滾道內(nèi)裝有若干滾球，其中轉(zhuǎn)盤中心位于O點(diǎn)，G處為不平衡質(zhì)量，偏心距為e，滾球在半徑為R的滾道上運(yùn)行，剛度和阻尼簡(jiǎn)化到水平和豎直方向上。在建立數(shù)學(xué)模型時(shí)，將其結(jié)構(gòu)簡(jiǎn)化為轉(zhuǎn)盤和滾球兩部分，并做如下假設(shè)：

1）轉(zhuǎn)盤為均質(zhì)圓形薄盤；

2）轉(zhuǎn)速視為按給定的函數(shù)變化；

3）球與滾道之間無(wú)相對(duì)滑動(dòng)；

4）轉(zhuǎn)盤處在水平方向，不考慮軸向運(yùn)動(dòng)及重力影響。

圖1　球式自動(dòng)平衡裝置結(jié)構(gòu)圖

所采用的符號(hào)表示如下（單位均采用國(guó)際標(biāo)準(zhǔn)單位）：

XOY：系統(tǒng)整體坐標(biāo)系。以X和Y表示轉(zhuǎn)盤中心O的坐標(biāo)。

Kx、KY：X方向的剛度系數(shù)和Y方向的剛度系數(shù)；CX、CY：X方向的粘性阻力系數(shù)和Y方向的粘性阻力系數(shù)。

θ：不平衡質(zhì)量相對(duì)于X軸的轉(zhuǎn)角；

φi：第i個(gè)滾球的轉(zhuǎn)角。

β0：滾球的滾動(dòng)摩擦系數(shù)；β1:滾球的粘性阻尼系數(shù)。

M：系統(tǒng)總質(zhì)量，m1、m2、m3、m分別為偏心質(zhì)量、轉(zhuǎn)盤質(zhì)量、電機(jī)質(zhì)量和滾球質(zhì)量。

I、I1、I2：滾球的轉(zhuǎn)動(dòng)慣量、偏心質(zhì)量的轉(zhuǎn)動(dòng)慣量、圓盤的轉(zhuǎn)動(dòng)慣量。

R、r：滾球轉(zhuǎn)動(dòng)半徑和球徑。

利用拉格朗日運(yùn)動(dòng)微分方程建立球式自動(dòng)平衡裝置的數(shù)學(xué)模型，其推導(dǎo)過(guò)程參照文獻(xiàn)［3］第17至19頁(yè)，得到系統(tǒng)微分方程如下

2　數(shù)值仿真系統(tǒng)的實(shí)現(xiàn)

采用四階龍格—庫(kù)塔法對(duì)上節(jié)給出的數(shù)學(xué)模型進(jìn)行數(shù)值分析求解。數(shù)值分析需要借助編程平臺(tái)，對(duì)于一般的數(shù)學(xué)計(jì)算，首選Matlab，但Matlab在可視化上不如通用編程語(yǔ)言靈活，因此本文的計(jì)算將采用VB，利用VB強(qiáng)大的圖形化功能增加仿真結(jié)果的直觀性，實(shí)現(xiàn)滾球運(yùn)動(dòng)狀態(tài)的動(dòng)態(tài)呈現(xiàn)。仿真程序流程如圖2所示。

圖2　仿真程序流程圖

2.1系統(tǒng)微分方程化簡(jiǎn)

將通用運(yùn)動(dòng)方程（1）化為便于仿真的通用格式1階微分方程組，以3球?yàn)槔?，?/p>

將定義的廣義坐標(biāo)代入式（1）得到式（2），觀察可知，式（2）為已知初值的通用格式1階微分方程組，可以采用龍格—庫(kù)塔法求解。

2.2仿真參數(shù)

仿真參數(shù)依據(jù)實(shí)驗(yàn)室原有的實(shí)驗(yàn)設(shè)備確定，見(jiàn)表1。系統(tǒng)中仿真參數(shù)為變量，可通過(guò)改變仿真參數(shù)進(jìn)行不同情況的仿真。

2.3仿真算法

仿真算法為4階龍格—庫(kù)塔法［4］。對(duì)于通用格式的1階微分方程組，如果初值可以確定，借助微機(jī)可以采用此算法求解。對(duì)于初值問(wèn)題

4階龍格—庫(kù)塔法求解下一步的方法為

上式中h為步長(zhǎng)，ki為

2.4碰撞調(diào)整

球式自動(dòng)平衡裝置運(yùn)動(dòng)過(guò)程中的滾球可能會(huì)相撞，因此在龍格-庫(kù)塔法一個(gè)步的求解后，需檢查球的位置是否干涉，如果干涉須做碰撞處理。對(duì)于在滾道中運(yùn)動(dòng)的滾球，由于存在摩擦和粘性阻尼，可視為完全非彈性碰撞，即指滾球碰撞結(jié)束后，這一瞬間內(nèi)兩物體的運(yùn)動(dòng)狀況完全相同，擁有相同的速度［5］。

兩球碰撞的判斷條件為兩球夾角φ滿足

完全非彈性碰撞后的結(jié)果為兩球分別往各自兩側(cè)移動(dòng)α

兩球速度均變?yōu)?/p>

2.5數(shù)值仿真系統(tǒng)界面

基于VB設(shè)計(jì)的數(shù)值仿真系統(tǒng)主要輸出有X和Y方向的振幅、轉(zhuǎn)盤中心的運(yùn)動(dòng)軌跡、各滾球相對(duì)不平衡質(zhì)量的轉(zhuǎn)角、滾球在轉(zhuǎn)盤內(nèi)運(yùn)動(dòng)的動(dòng)畫等。

該仿真系統(tǒng)具備計(jì)算0個(gè)滾球（視為不平衡轉(zhuǎn)子）和1～3個(gè)滾球（視為不同的球式自動(dòng)平衡裝置）系統(tǒng)各項(xiàng)參數(shù)的能力。

軟件界面分三列，第一列四幅圖顯示振幅，分別為無(wú)滾球和有滾球情況下X和Y方向的振幅輸出，在無(wú)滾球計(jì)算完成后，保留振幅輸出圖像，便于同有滾球的時(shí)候形成振幅對(duì)比。中間列第一幅圖為角速度輸入，中間列后三幅圖為各滾球相對(duì)不平衡質(zhì)量的轉(zhuǎn)角。最右邊列第一幅圖動(dòng)態(tài)顯示滾球在轉(zhuǎn)盤內(nèi)的運(yùn)動(dòng)情況，本模塊可以根據(jù)計(jì)算得到的轉(zhuǎn)角，以動(dòng)畫的方式顯示滾球運(yùn)動(dòng)情況，不平衡質(zhì)量在X軸正向，滾球以轉(zhuǎn)盤為參考系運(yùn)行，最右邊列下面一幅顯示轉(zhuǎn)盤中心運(yùn)動(dòng)軌跡，也就是振幅在XY平面內(nèi)的點(diǎn)云圖像。右下角的四個(gè)按鈕可以根據(jù)滾球數(shù)點(diǎn)擊，用來(lái)啟動(dòng)仿真。

表1　仿真參數(shù)（未注單位為國(guó)際單位）

3　仿真系統(tǒng)的仿真結(jié)果

根據(jù)自同期理論［6］預(yù)測(cè)滾球在不同轉(zhuǎn)速下的位置：當(dāng)轉(zhuǎn)速較低時(shí)，滾球向偏心質(zhì)量一側(cè)移動(dòng)，此時(shí)會(huì)造成轉(zhuǎn)子振幅增大；當(dāng)轉(zhuǎn)速大于1階臨界轉(zhuǎn)速時(shí)，如果只是略大于，可能會(huì)產(chǎn)生自激振蕩，如果大得多，則滾球開始遠(yuǎn)離偏心質(zhì)量，達(dá)到平衡的目的［7］。下面分析四種不同滾球數(shù)目下的輸出，與自同期理論形成對(duì)比，驗(yàn)證仿真系統(tǒng)結(jié)果的可信性。

殘留振幅是衡量過(guò)臨界轉(zhuǎn)速時(shí)移動(dòng)物體對(duì)不平衡轉(zhuǎn)子減振能力的標(biāo)志，最大振幅是決定轉(zhuǎn)子啟動(dòng)時(shí)能否安全通過(guò)共振區(qū)的一項(xiàng)指標(biāo)，因此最大振幅、殘留振幅作為評(píng)價(jià)自動(dòng)平衡回轉(zhuǎn)機(jī)械工作機(jī)理影響因素的指標(biāo)［3］。由于X、Y方向振幅規(guī)律相近，本文僅分析X方向振幅。

為了使仿真盡量接近實(shí)際情況，角速度輸入采用從0以5 s時(shí)間加速到25 Hz后勻速轉(zhuǎn)動(dòng)。仿真時(shí)間持續(xù)30 s，每0.008 s取一個(gè)步，總共3 750個(gè)計(jì)算步。這里指出，仿真步的大小要選擇合適，步太大影響結(jié)果的準(zhǔn)確性，步太小會(huì)導(dǎo)致仿真時(shí)間過(guò)長(zhǎng)。

3.1不裝滾球

不裝滾球下的球式自動(dòng)平衡裝置可視為普通不平衡轉(zhuǎn)子，其水平方向的振幅如圖3所示，從圖中的振幅走向可以看出，剛開始由于轉(zhuǎn)速不大，轉(zhuǎn)子振幅較小，之后轉(zhuǎn)子轉(zhuǎn)速達(dá)到一階臨界轉(zhuǎn)速后通過(guò)共振區(qū)，達(dá)到最大振幅，最后逐漸衰減為穩(wěn)定值。最大振幅值和殘留振幅值可從圖上獲得（水平方向最大振幅由Xmax表示，殘留振幅由Xbal表示）。

圖3　未裝滾球的計(jì)算結(jié)果

3.2 1個(gè)滾球

球式自動(dòng)平衡裝置內(nèi)裝有1個(gè)滾球的情況如圖4所示，滾球的原始位置相對(duì)不平衡質(zhì)量的轉(zhuǎn)角為0，其X方向的振幅如圖4（a）所示。從圖中的振幅走向可以看出，剛開始由于轉(zhuǎn)速不大，轉(zhuǎn)子振幅較小，之后轉(zhuǎn)子轉(zhuǎn)速達(dá)到1階臨界轉(zhuǎn)速后通過(guò)共振區(qū)，達(dá)到最大振幅，該振幅超過(guò)了不裝滾球時(shí)的最大振幅，說(shuō)明此時(shí)滾球加大了不平衡量，最后逐漸衰減為一個(gè)穩(wěn)定值，該穩(wěn)定值低于未裝滾球時(shí)的殘留振幅，說(shuō)明滾球已經(jīng)起到平衡不平衡質(zhì)量的作用，符合自同期現(xiàn)象。圖4（b）顯示轉(zhuǎn)盤的中心運(yùn)動(dòng)軌跡。最大振幅值和殘留振幅值具體數(shù)值如下

Xmax=1.699×10-3m，Xbal=2.270×10-6m滾球最后相對(duì)不平衡質(zhì)量的轉(zhuǎn)角為φ1'=-53.396 rad，化為角度，得滾球與不平衡質(zhì)量最后的夾角為180.63°，接近理想角度180°。圖4（c）中的區(qū)域動(dòng)態(tài)顯示滾球在轉(zhuǎn)過(guò)8圈后，緩慢移動(dòng)到最后停留位置。

3.3 2個(gè)滾球

球式自動(dòng)平衡裝置內(nèi)裝有2個(gè)滾球，滾球初始轉(zhuǎn)角為0和π，基本情況同裝有1個(gè)滾球時(shí)相似，區(qū)別是兩個(gè)滾球可能會(huì)發(fā)生碰撞，上述初值下未發(fā)生碰撞，最大振幅同樣大于不裝滾球時(shí)的振幅。振幅穩(wěn)定值低于未裝滾球時(shí)的振幅穩(wěn)定值，減振效益十分明顯，說(shuō)明滾球已經(jīng)起到平衡不平衡質(zhì)量的作用。最大振幅值和殘留振幅值及滾球停留位置具體數(shù)據(jù)如下

化為角度，得滾球1和滾球2與不平衡質(zhì)量最后的夾角分別為118.34°和238.78°，分別接近理想角度120°和240°。球式自動(dòng)平衡裝置裝有2個(gè)滾球的計(jì)算結(jié)果如圖5所示。

圖4　裝有1個(gè)滾球的計(jì)算結(jié)果

3.4 3個(gè)滾球

球式自動(dòng)平衡裝置內(nèi)裝有3個(gè)滾球的情況如圖

6所示，滾球初始轉(zhuǎn)角為0、2 π/3和4 π/3，其X方向的振幅如圖6（a）所示。從圖中的振幅走向可以看出，剛開始由于轉(zhuǎn)速不大，轉(zhuǎn)子振幅較小，之后轉(zhuǎn)子轉(zhuǎn)速達(dá)到1階臨界轉(zhuǎn)速后通過(guò)共振區(qū)，達(dá)到最大振幅，由于初始位置比較理想，該振幅略低于不裝滾球時(shí)的最大振幅，實(shí)際上此時(shí)滾球加大了不平衡量。最后逐漸衰減為穩(wěn)定值，該穩(wěn)定值低于未裝滾球時(shí)的殘留振幅，說(shuō)明滾球已起到平衡不平衡質(zhì)量的作用。這里指出在上述初始轉(zhuǎn)角下，滾球之間有相互碰撞。最大振幅值和殘留振幅值及滾球轉(zhuǎn)角具體數(shù)值如下

圖5　裝有2個(gè)滾球的計(jì)算結(jié)果

圖6　裝有3個(gè)滾球的計(jì)算結(jié)果

化為角度，得滾球1、滾球2和滾球3與不平衡質(zhì)量最后夾角分別為270.23°、89.40°和179.52°，接近理想角度270°、90°和180°。圖6（c）顯示了滾球的最后停留位置。

4　結(jié)語(yǔ)

基于VB設(shè)計(jì)的球式自動(dòng)平衡裝置數(shù)值仿真系統(tǒng)，通過(guò)參照（未裝滾球）和三組帶滾球的仿真，得到的仿真結(jié)果與自同期理論一致，不平衡轉(zhuǎn)子的振幅都得到了有效的衰減，減振效益達(dá)到40倍左右。

仿真過(guò)程中發(fā)現(xiàn)，減振速度與仿真步的大小有關(guān)，仿真步越小，滾球達(dá)到穩(wěn)定位置的時(shí)間越長(zhǎng)；殘留振幅與仿真時(shí)間有關(guān)，仿真時(shí)間越長(zhǎng)，殘留振幅越小。

該數(shù)值仿真系統(tǒng)實(shí)現(xiàn)了球式自動(dòng)平衡裝置仿真結(jié)果的動(dòng)態(tài)可視化，增加仿真結(jié)果的直觀性，實(shí)現(xiàn)了滾球運(yùn)動(dòng)狀態(tài)的動(dòng)態(tài)呈現(xiàn)。通過(guò)修改仿真系統(tǒng)的參數(shù)，可以很便捷的分析仿真結(jié)果的影響因素，包括初始位置、滾道摩擦系數(shù)、滾道阻尼系數(shù)、轉(zhuǎn)子旋轉(zhuǎn)加速度等。

［1］劉文倩，譚青，謝燕琴，等.自動(dòng)平衡裝置減振效益分析［J］.噪聲與振動(dòng)控制，2010，30（4）：153-157.

［2］譚青，周鐵，黃秀祥.球式自動(dòng)平衡裝置的數(shù)值仿真與實(shí)驗(yàn)研究［J］.噪聲與振動(dòng)控制，2008，28（3）：142-145.

［3］周鐵.內(nèi)裝滾球轉(zhuǎn)子自同期現(xiàn)象的理論與實(shí)驗(yàn)研究［D］.長(zhǎng)沙：中南大學(xué)，2008.

［4］Brooks Cole.馮煙利譯.數(shù)值分析［M］.北京：高等教育出版社，2005.

［5］哈爾濱工業(yè)大學(xué)理論力學(xué)教研室.理論力學(xué)（第七版）［M］.北京：高等教育出版社，2010.

［6］井上順吉，荒木嘉昭.振動(dòng)機(jī)械的自己同期化研究（日文）［C］.日本機(jī)械學(xué)會(huì)論文集，1967，33（2）.246.

［7］楊義蛟，譚青，周鐵.自動(dòng)平衡回轉(zhuǎn)機(jī)械的工作機(jī)理研究［J］.機(jī)械設(shè)計(jì)與研究，2009，（4）：33-36.

Design of Numerical Simulation System of Ball-type Automatic Balancers

TANQing，REN Zhi-xiang，WEILai
（College of Mechanical&Electrical Engineering，Central South University，Changsha 410083，China）

The Matlab code is used very often in numerical simulation of ball-type automatic balancers.But the results are neither intuitive nor interactive.In this paper，a new dynamic simulation system with VB platform was designed.The mathematical model of the ball-type automatic balancers was solved numerically using Runge-Kutta method.The simulation results show that the system enables interactive visualization and intuition of the simulation results and realizes the dynamic presentation of the state of dynamic motion of the balls.The simulation results obtained by the numerical simulation system of the ball-type automatic balancers are consistent with those of the self-synchronization theory.The amplitude of the unbalanced rotor is reduced effectively.

vibration and wave；ball-type automatic balancer；mathematical model；numerical simulation system；vibration reduction

TP311；TH113.1

ADOI編碼：10.3969/j.issn.1006-1335.2015.05.009

1006-1355（2015）05-0048-05

2015－02－12

任志湘（1990－），女，湖南省岳陽(yáng)市人，碩士生，主要研究方向：自動(dòng)平衡裝置虛擬實(shí)驗(yàn)平臺(tái)。

Email：jds-share@163.com

譚青（1955－），男，博士生導(dǎo)師。