大功率往復(fù)泵曲軸轉(zhuǎn)子動(dòng)力學(xué)特性研究*

□ 張俊峰 □ 張春光 □ 劉光恒 □ 劉佳慧 □ 褚 帥 □ 安 磊

沈陽(yáng)儀表科學(xué)研究院有限公司 沈陽(yáng) 110043

1 研究背景

隨著社會(huì)的發(fā)展與基建規(guī)模的擴(kuò)大,大功率往復(fù)泵的需求不斷提升。大功率往復(fù)泵利用高壓水射流對(duì)鋼筋混凝土破拆,不損壞鋼筋,廣泛應(yīng)用于道路施工、橋梁改建、高速公路維護(hù)、市政工程等領(lǐng)域,成為近幾年工程領(lǐng)域的一大新熱點(diǎn)。曲軸是大功率往復(fù)泵動(dòng)力端的關(guān)鍵零部件,其質(zhì)量往往直接影響大功率往復(fù)泵的性能。隨著大功率往復(fù)泵功率的提高,曲軸性能的要求也越來(lái)越高。曲軸是一種典型的旋轉(zhuǎn)機(jī)械,對(duì)其進(jìn)行轉(zhuǎn)子動(dòng)力學(xué)設(shè)計(jì)是重要環(huán)節(jié),臨界轉(zhuǎn)速、轉(zhuǎn)子不平衡響應(yīng)和轉(zhuǎn)子啟動(dòng)軸心軌跡是研究的重點(diǎn)。近年來(lái),隨著計(jì)算機(jī)計(jì)算性能的顯著提高,利用有限元法來(lái)分析轉(zhuǎn)子動(dòng)力學(xué)特性已成為一種主流工程方法[1]。在有限元模型中考慮轉(zhuǎn)動(dòng)慣量、陀螺力矩、軸向載荷等因素的影響,使計(jì)算相較于其它方法精度更高。筆者采用有限元法對(duì)大功率往復(fù)泵曲軸轉(zhuǎn)子進(jìn)行動(dòng)力學(xué)特性研究。

2 泵動(dòng)力端三維模型

筆者對(duì)轉(zhuǎn)子軸承系統(tǒng)建立動(dòng)力學(xué)模型,一般將轉(zhuǎn)子系統(tǒng)簡(jiǎn)化為彈性軸段、離散剛性圓盤及軸承支撐單元[2],各單元之間由節(jié)點(diǎn)連接,進(jìn)而建立系統(tǒng)的總運(yùn)動(dòng)方程。

根據(jù)大功率往復(fù)泵的性能,對(duì)泵動(dòng)力端整體傳動(dòng)方案和關(guān)鍵零部件進(jìn)行設(shè)計(jì),并利用計(jì)算機(jī)輔助設(shè)計(jì)軟件進(jìn)行三維模型構(gòu)建。泵動(dòng)力端三維模型如圖1所示,其中曲軸采用五拐四支撐布局。曲軸轉(zhuǎn)子三維模型如圖2所示,主要尺寸參數(shù)見表1。

3 剛性圓盤運(yùn)動(dòng)方程

設(shè)圓盤中心節(jié)點(diǎn)的廣義坐標(biāo){u1d}=[x,θy]T,{u2d}=[y,-θx]T,剛性圓盤質(zhì)心與軸心重合,根據(jù)拉格朗日定理,剛性圓盤的運(yùn)動(dòng)方程為:

(1)

(2)

(3)

式中:Ω為回轉(zhuǎn)體的角速度;{Q1d}、{Q2d}為廣義力;m為剛性圓盤質(zhì)量;Jd為過(guò)軸心的徑向轉(zhuǎn)動(dòng)慣量;Jp為極轉(zhuǎn)動(dòng)慣量。

4 彈性軸段運(yùn)動(dòng)方程

設(shè)彈性軸段單元廣義坐標(biāo){u1s}=[xA,θyA,xB,θyB]T,{u2s}=[yA,-θxA,yB,-θxB]T,根據(jù)拉格朗日定理,軸段單元的運(yùn)動(dòng)方程為:

(4)

[Ms]=[MsR]+[MsT]

(5)

式中:[Ms]為包含轉(zhuǎn)動(dòng)慣性矩陣與移動(dòng)慣性矩陣的質(zhì)量矩陣;[MsR]為轉(zhuǎn)動(dòng)慣性矩陣;[MsT]為移動(dòng)慣性矩陣;[Js]為回轉(zhuǎn)矩陣;[Ks]為彈性軸段剛度矩陣。

根據(jù)彈性力學(xué)有限元法,曲軸轉(zhuǎn)子在笛卡兒坐標(biāo)系中的運(yùn)動(dòng)微分方程為:

(6)

式中:[M]為總體質(zhì)量矩陣;[K]為總體剛度矩陣;[C]為由阻尼矩陣和陀螺力矩陣組成的矩陣;{F(t)}為外部激振力。

5 曲軸轉(zhuǎn)子有限元模型

將曲軸轉(zhuǎn)子支撐系統(tǒng)作為研究對(duì)象。曲軸轉(zhuǎn)子為五拐四支撐結(jié)構(gòu)布局,采用SOLID185單元類型,單元數(shù)為503 243,節(jié)點(diǎn)數(shù)為195 929,網(wǎng)格質(zhì)量達(dá)到計(jì)算精度要求。軸承采用COMBI214單元類型模擬,這一單元類型可以模擬交叉剛度和阻尼,還能模擬拉壓剛度,但不能模擬彎曲或扭轉(zhuǎn)剛度。曲軸材料彈性模量為210 GPa,密度為7.8×103kg/m3,泊松比為0.3。模型Z方向?yàn)榍S轉(zhuǎn)子旋轉(zhuǎn)方向,假設(shè)軸承各向同性,且不受阻尼影響。

在曲軸轉(zhuǎn)子四個(gè)支撐位的截面中心建立主節(jié)點(diǎn),在截面外圓周建立從節(jié)點(diǎn),創(chuàng)建剛性區(qū)域。在每個(gè)支撐位置通過(guò)創(chuàng)建COMBI214單元來(lái)模擬軸承,單元某一端與剛性區(qū)域主節(jié)點(diǎn)相連,另一端施加全約束。曲軸轉(zhuǎn)子有限元模型如圖3所示。

▲圖3 曲軸轉(zhuǎn)子有限元模型

6 模態(tài)分析

大功率往復(fù)泵曲軸轉(zhuǎn)子是典型的旋轉(zhuǎn)零件,相較于一般零部件,模態(tài)分析時(shí)需考慮陀螺效應(yīng)對(duì)轉(zhuǎn)子的影響[3-5],應(yīng)將陀螺效應(yīng)設(shè)置為打開。筆者使用QR阻尼法對(duì)轉(zhuǎn)子系統(tǒng)進(jìn)行模態(tài)分析,這一方法可以很好地對(duì)低阻尼系統(tǒng)進(jìn)行求解,支持不對(duì)稱阻尼和不對(duì)稱剛度。計(jì)算得到曲軸轉(zhuǎn)子前六階固有頻率見表2,前六階模態(tài)振型如圖4所示。

▲圖4 曲軸轉(zhuǎn)子前六階模態(tài)振型

由于考慮到陀螺效應(yīng),曲軸轉(zhuǎn)子轉(zhuǎn)速對(duì)特征頻率有一定影響。計(jì)算不同轉(zhuǎn)速所對(duì)應(yīng)的固有頻率,并繪制固有頻率隨轉(zhuǎn)速變化曲線,這一曲線圖稱為坎貝爾圖[6-8]。在坎貝爾圖上,從縱軸出發(fā)曲線斜率為正的直線代表正進(jìn)動(dòng),斜率為負(fù)的直線代表反進(jìn)動(dòng),從原點(diǎn)出發(fā)的轉(zhuǎn)速直線與正進(jìn)動(dòng)直線相交點(diǎn)所對(duì)應(yīng)的橫坐標(biāo)代表臨界轉(zhuǎn)速。一般計(jì)算臨界轉(zhuǎn)速時(shí)只考慮正進(jìn)動(dòng)轉(zhuǎn)速。

表2 曲軸轉(zhuǎn)子前六階固有頻率

載荷步分為40步,轉(zhuǎn)速為0～600 r/min,利用模態(tài)分析繪制曲軸轉(zhuǎn)子坎貝爾圖,如圖5所示。

在圖5中可以看到,從原點(diǎn)出發(fā)的轉(zhuǎn)速直線并沒有與正進(jìn)動(dòng)直線相交,也就是說(shuō)曲軸轉(zhuǎn)子工作轉(zhuǎn)速?zèng)]有達(dá)到臨界轉(zhuǎn)速,旋轉(zhuǎn)結(jié)構(gòu)不會(huì)產(chǎn)生共振,運(yùn)行安全,滿足設(shè)計(jì)要求。

▲圖5 曲軸轉(zhuǎn)子坎貝爾圖

7 不平衡響應(yīng)分析

由力學(xué)知識(shí)可知,大功率往復(fù)泵曲軸轉(zhuǎn)子由于質(zhì)量偏心引起的離心力F為:

F=meω2

(7)

式中:m為偏心質(zhì)量;e為偏心距;ω為曲軸轉(zhuǎn)子轉(zhuǎn)速。

在理想情況下,無(wú)需考慮軸向力,離心力被分解為X方向與Y方向的減諧力,為:

(8)

式中:t為時(shí)間。

利用有限元法進(jìn)行不平衡分析[9],打開陀螺效應(yīng),設(shè)置不平衡響應(yīng)的頻率為0～1 500 Hz,子步數(shù)為40。經(jīng)分析得到曲軸轉(zhuǎn)子不平衡響應(yīng)結(jié)果,如圖6所示,橫坐標(biāo)為掃頻頻率,縱坐標(biāo)為曲軸轉(zhuǎn)子軸心的振動(dòng)幅度。曲軸轉(zhuǎn)子在525 Hz、637.5 Hz、900 Hz、1 200 Hz、1 425 Hz頻率處振動(dòng)幅度較大,即在這些頻率處產(chǎn)生較強(qiáng)共振,與模態(tài)分析的結(jié)果基本相同,由此驗(yàn)證了模態(tài)分析的準(zhǔn)確性。

▲圖6 曲軸轉(zhuǎn)子不平衡響應(yīng)結(jié)果

8 啟動(dòng)過(guò)程響應(yīng)分析

大功率往復(fù)泵曲軸轉(zhuǎn)子由電機(jī)驅(qū)動(dòng),啟動(dòng)時(shí)間較快,1 s時(shí)間從0加速到39.76 rad/s,隨后進(jìn)入穩(wěn)定轉(zhuǎn)速狀態(tài)。曲軸轉(zhuǎn)子啟動(dòng)時(shí)角速度變化如圖7所示。

▲圖7 曲軸轉(zhuǎn)子啟動(dòng)時(shí)角速度變化

整個(gè)加速過(guò)程中角加速度為39.76 rad/s2。曲軸轉(zhuǎn)子動(dòng)不平衡等級(jí)選擇為G6.3,由此得到曲軸轉(zhuǎn)子加速過(guò)程中的激振力為:

(9)

對(duì)曲軸轉(zhuǎn)子進(jìn)行瞬態(tài)動(dòng)力學(xué)響應(yīng)分析,得到曲軸轉(zhuǎn)子軸心的運(yùn)動(dòng)軌跡,如圖8所示。從圖8中可以看到,曲軸轉(zhuǎn)子軸心的最大響應(yīng)振幅為4.72×10-8m,振幅很小,沒有產(chǎn)生共振,并且軸心的運(yùn)動(dòng)軌跡經(jīng)過(guò)啟動(dòng)時(shí)螺旋線后達(dá)到穩(wěn)定運(yùn)行狀態(tài),近似圓形運(yùn)動(dòng)軌跡,符合預(yù)期要求。

▲圖8 曲軸轉(zhuǎn)子軸心運(yùn)動(dòng)軌跡

9 結(jié)束語(yǔ)

為了驗(yàn)證大功率往復(fù)泵曲軸轉(zhuǎn)子設(shè)計(jì)的合理性,筆者對(duì)曲軸轉(zhuǎn)子進(jìn)行了模態(tài)分析、坎貝爾圖分析、不平衡響應(yīng)分析與啟動(dòng)過(guò)程響應(yīng)分析。

對(duì)曲軸轉(zhuǎn)子進(jìn)行模態(tài)分析,得到曲軸轉(zhuǎn)子的前六階固有頻率及振型,充分了解了曲軸轉(zhuǎn)子的振動(dòng)情況,為設(shè)計(jì)驗(yàn)證提供依據(jù)。

對(duì)曲軸轉(zhuǎn)子進(jìn)行坎貝爾圖分析,發(fā)現(xiàn)從原點(diǎn)出發(fā)的轉(zhuǎn)速直線并沒有與正進(jìn)動(dòng)直線相交,也就是說(shuō)曲軸轉(zhuǎn)子工作轉(zhuǎn)速?zèng)]有達(dá)到臨界轉(zhuǎn)速,旋轉(zhuǎn)結(jié)構(gòu)不會(huì)產(chǎn)生共振,運(yùn)行安全,滿足設(shè)計(jì)要求。

根據(jù)計(jì)算的不平衡量對(duì)曲軸轉(zhuǎn)子進(jìn)行不平衡響應(yīng)分析,得到不平衡響應(yīng)結(jié)果及共振頻率點(diǎn),確認(rèn)與模態(tài)分析的結(jié)果基本相同,驗(yàn)證了模態(tài)分析的準(zhǔn)確性。

對(duì)曲軸轉(zhuǎn)子進(jìn)行啟動(dòng)過(guò)程響應(yīng)分析,得到啟動(dòng)瞬間與穩(wěn)定運(yùn)行時(shí)的運(yùn)動(dòng)軌跡,確認(rèn)曲軸轉(zhuǎn)子軸心的振動(dòng)幅度很小,沒有出現(xiàn)共振情況,由此驗(yàn)證了設(shè)計(jì)的合理性。