馬曉偉，劉 健，孫延迪，劉忠硯，肖柳勝

(中國石油大學(xué)(華東)機(jī)電工程學(xué)院，山東青島266580)

0 引言

近幾年來隨著油田增產(chǎn)措施的推進(jìn)，壓裂采油工藝需要的壓力等級越來越高［1-2］，壓裂車也朝著大功率、大型化快速發(fā)展，而隨之產(chǎn)生的振動問題也日益突出。一些新研制的大型壓裂車，作業(yè)工況下某些檔位出現(xiàn)了較大振動，已經(jīng)嚴(yán)重影響到整車的使用性能、安全及壽命。從已有的分析可以明確，三缸泵的振動激勵是引起壓裂車作業(yè)工況振動的主要原因，而如何減小三缸泵的振動激勵向車架的傳遞也成為目前制約壓裂車壓力等級提高的技術(shù)瓶頸。因此，明確三缸泵對車架的實際激勵大小及形式，進(jìn)而采取合理的隔振措施，對于壓裂車的大功率、大型化發(fā)展具有重要的意義。

三缸泵結(jié)構(gòu)復(fù)雜、激勵源多，目前國內(nèi)對其振動進(jìn)行全面分析研究的較少，僅有少量文獻(xiàn)針對其局部振動的研究，例如文獻(xiàn)［3-4］仿真分析了三缸泵的曲軸和滑塊部分的受力情況，模型簡單，忽略因素較多。

本研究綜合考慮三缸泵的慣性力和液力端壓力波動等多激勵源工況，其能夠真實地反映壓裂車作業(yè)時三缸泵對車架的振動激勵。

1 振動激勵源分析

三缸泵可以看成是由3 個單缸泵按相位角相差120°排列而成［5］。因此分析三缸泵的一個缸的曲柄滑塊機(jī)構(gòu)慣性力激勵及液力端壓力波動激勵，將其按相位角相差120°進(jìn)行疊加即可得到三缸泵的相應(yīng)的振動激勵。

1.1 曲柄滑塊機(jī)構(gòu)的慣性力

三缸泵的單個缸的原理圖如圖1所示。

其主體為典型的曲柄滑塊機(jī)構(gòu)，慣性力主要包括曲軸的旋轉(zhuǎn)慣性力、連桿的旋轉(zhuǎn)慣性力和旋轉(zhuǎn)慣性力矩、柱塞的往復(fù)慣性力。

圖1 曲柄滑塊機(jī)構(gòu)圖

設(shè)曲柄AB 長度為R，曲柄角速度為ω，連桿長度為L，連桿瞬時角速度和角加速度分別為ω2和ε2，連桿和柱塞的質(zhì)量分別為m2和m3，由文獻(xiàn)［6-7］可求得連桿AB 的角速度ω2和角加速度ε2及柱塞的加速度a，此處不再做推導(dǎo)，因此可以求得:

(1)曲軸的慣性力F=m×R×ω2;

(2)連桿的慣性力F2=m2×r2×ω22，慣性力矩M=J×ε2;

(3)柱塞的慣性力F3=m3×a。

1.2 液力端的液壓波動力

三缸泵柱塞上激勵力:

FP=S×P

式中:S—柱塞端面積;P—柱塞缸內(nèi)壓力。

由于三缸泵液力端存在泵入、泵出液體時壓力的波動，加載在柱塞端的液力激勵為隨著液力端壓力波動的變量，泵入液體時的壓力較小，本研究在此忽略［8］，僅考慮泵出液體時的壓力。

參考圖1，假設(shè)其泵出液體時的壓力恒為液力端壓力，當(dāng)曲軸轉(zhuǎn)過角度0≤φ≤π 時，泵入液體，π≤φ≤2π時，泵出液體，則加載在柱塞上的激勵力變化如圖2所示。

圖2 柱塞端激勵力

2 基于ADAMS 三缸泵的虛擬樣機(jī)

由于三缸泵結(jié)構(gòu)復(fù)雜，而ADAMS 軟件三維建模能力一般［9］，本研究采用Pro/E 來建立三缸泵的各部分模型并裝配，然后導(dǎo)入ADAMS 軟件中，添加相應(yīng)的材料屬性、連接約束、驅(qū)動和激勵。

2.1 Pro/E 中三缸泵模型的建立

以某公司的3000 型壓裂車車載三缸泵為研究原型，為節(jié)省建模時間和簡化分析，本研究對于三缸泵的模型做了適當(dāng)?shù)暮喕?

(1)建模時忽略對分析影響不大的連接螺栓，螺栓連接處全部用剛性連接代替;

(2)建模時對曲軸輸入端進(jìn)行簡化，忽略曲軸輸入端減速器，僅建立位于泵中的曲軸;

(3)為便于在ADAMS 中添加各種約束副，建模時將十字頭、小連桿及柱塞當(dāng)成一個零件處理。

本研究建立的三缸泵三維模型如圖3所示。

圖3 三缸泵三維模型

2.2 ADAMS 中三缸泵虛擬樣機(jī)的構(gòu)建

筆者把Pro/E 導(dǎo)出的三缸泵中性體文件導(dǎo)入ADAMS，設(shè)置各部分的材料屬性，三缸泵主體曲柄滑塊部分材質(zhì)為42CrMo，其彈性模量為2.12 ×1011N/m2，泊松比為0.28，密度為7.85 ×103kg/m3，添加完材料屬性之后，添加各部件之間的連接約束。

本研究更真實地反應(yīng)三缸泵工況時對車架的振動激勵，考慮了各運動部件間的摩擦力，由于篇幅有限，在此不詳細(xì)列出各部分約束及摩擦系數(shù)，完成相應(yīng)設(shè)置的虛擬樣機(jī)主體部分如圖4所示。

圖4 主體曲柄滑塊部分約束

3 虛擬樣機(jī)仿真實驗及振動激勵分析

3.1 驅(qū)動及振動激勵的設(shè)置

本研究假定發(fā)動機(jī)在大泵達(dá)到穩(wěn)定工況時的輸出功率是恒定值，且到三缸泵輸入端的傳遞率不隨其他因素變化而變化。

由:

P=T×ω

式中:P—大泵輸入功率，T—轉(zhuǎn)矩，ω—轉(zhuǎn)速。

即可求出大泵輸入扭矩和轉(zhuǎn)速的關(guān)系如表1所示。

表1 輸入轉(zhuǎn)速、扭矩關(guān)系

激勵部分主要包括慣性力和液力端壓力［10-11］，慣性力通過在虛擬樣機(jī)中設(shè)置約束和材料屬性，進(jìn)行仿真時虛擬樣機(jī)會自動計算。而對于液力端壓力，為了縮短仿真時間，本研究假設(shè)壓裂車的大泵在每個檔位工作2 s，由于在吸液沖程壓力可忽略不計，液力端壓力波動主要考慮排液沖程。由柱塞作用力為柱塞面積(11 432.6 mm2)和最大排出壓力(100 MPa)可知，加載柱塞端的液體壓力值為F=P×S=1 143 260 N。其第一個活塞端的加載函數(shù)為(其余兩個活塞加載函數(shù)與第一個活塞的加載函數(shù)相差120°和240°相位角):

3 個柱塞的壓力加載曲線如圖5所示。

圖5 3 個柱塞壓力加載曲線

3.2 振動激勵的分析

本研究設(shè)置仿真分析時間為16 s，steps 為160。由于三缸泵引起的壓裂車的振動主要表現(xiàn)為橫向振動，此處主要關(guān)心三缸泵工作時對車架的橫向激勵，仿真分析車架對三缸泵的橫向支反力，并對支反力進(jìn)行FFT 變換。

支反力及支反力FFT 變換如圖6所示。

圖6 泵座橫向支反力及支反力FFT 變換

由圖6(b)可知，壓裂車在剛啟動不加壓作業(yè)時，工作平穩(wěn)，無明顯橫向振動產(chǎn)生;加壓后，出現(xiàn)明顯橫向振動，且在四檔、五檔時振動幅度較大。由仿真數(shù)據(jù)可以看出，車架對三缸泵的橫向支反力幅值只受液力端的壓力波動影響，且峰值出現(xiàn)在5.3 Hz 處，7 s～9 s(6 s～8 s 為四檔，8 s～10 s 為五檔)時，這與壓裂車實際作業(yè)工況時在四檔和五檔處出現(xiàn)低頻橫向振動吻合。因此，由仿真分析可以得出，引起壓裂車作業(yè)工況整車橫向振動的激勵源主要為三缸泵加壓時液力端的壓力波動引起的，且其激勵主要成分為低頻激勵。

3.3 數(shù)據(jù)可靠性的試驗驗證

為驗證仿真分析數(shù)據(jù)的可靠性，本研究通過將仿真得到的載荷數(shù)據(jù)加載在壓裂車虛擬樣機(jī)上，測得不同部位的加速度響應(yīng)值，與現(xiàn)場實驗實際工況加速度傳感器測得的數(shù)據(jù)進(jìn)行對比。

壓裂車虛擬樣機(jī)如圖7所示。

圖7 壓裂車虛擬樣機(jī)

壓裂車現(xiàn)場實驗及加速度傳感器4464、4459、4458 在車架上的安裝位置如圖8所示。

圖8 現(xiàn)場實驗

由于篇幅所限，本研究著重對比分析壓裂車四檔時，4464、4459 及4458 位置的試驗測得的加速度峰值和仿真測得的加速度峰值對比。

數(shù)據(jù)對比如表2所示。

從表2 可以看出，4464、4459 及4458 位置仿真測得的加速度峰值與試驗數(shù)據(jù)均較為接近，其數(shù)據(jù)誤差在20%以內(nèi)，滿足工程誤差要求，可知本研究仿真的

表2 數(shù)據(jù)對比

三缸泵對車架的激勵數(shù)據(jù)能夠真實反應(yīng)壓裂車實際作業(yè)工況三缸泵對車架的振動激勵情況，具有較好的工程應(yīng)用價值，確保了壓裂車三缸泵隔振方案設(shè)計所需數(shù)據(jù)的正確性與可靠性。

4 結(jié)束語

(1)本研究綜合考慮三缸泵的慣性力和液力端壓力波動等多激勵源，建立了貼近實際工況的三缸泵的虛擬樣機(jī)模型，仿真得到的結(jié)果與壓裂車的實際工作特性吻合;

(2)通過加載仿真激勵數(shù)據(jù)，本研究得到與實際試驗相同位置點的加速度響應(yīng)值，其誤差均在20%以內(nèi)，證明了虛擬樣機(jī)仿真結(jié)果是可靠的;

(3)本研究明確了三缸泵振動引起壓裂車橫向振動的激勵的幅值及形式，為解決壓裂車三缸泵的振動問題及隔振方案的設(shè)計提供了依據(jù)。

［1］劉 健，劉忠硯，龐 罕，等.大型壓裂泵車車架實驗應(yīng)變模態(tài)分析方法研究［J］.機(jī)電工程，2014，31(2):150-153.

［2］劉 健，肖柳勝，吳漢川，等.油田壓裂車車架的多工況疲勞壽命分析［J］.石油機(jī)械，2014，42(1):75-78.

［3］桑曙東，鐘功祥，曹炳軍.基于ADAMS 的3NB-1300 泥漿泵運動學(xué)分析［J］.裝備制造技術(shù)，2014(5):27-29.

［4］宋保華，王立梅，賀建隆.三缸單作用泥漿泵曲軸受力分析［J］.化學(xué)工程與裝備，2011(6):131-134.

［5］綦耀光，石臨嵩，周鳳石.鉆井泵設(shè)計中的兩個理論問題［J］.石油機(jī)械，1990，18(9):13-14.

［6］黃錫愷，鄭文緯.機(jī)械原理［M］.北京:人民教育出版社，1981.

［7］盧其宜，謝菊平.曲柄滑塊機(jī)構(gòu)中連桿加速度瞬心位置的確定［J］.中國水運:學(xué)術(shù)版，2006，8(7):50-51.

［8］朱永有，葉永彪，周 容，等.泥漿泵運動學(xué)、動力學(xué)分析［J］.研發(fā)與制造，2005(5):65-67.

［9］丁乾坤，劉正士.Pro-E 與ADAMS 聯(lián)合對發(fā)動機(jī)振動的研究［J］.現(xiàn)代車用動力，2008(4):34-38.

［10］安 平，趙文宏，謝良江，等.基于LabVIEW 的磨球機(jī)振動信壓分析系統(tǒng)［J］.輕工機(jī)械，2014，32(1):42-44.

［11］趙俊龍.往復(fù)式壓縮機(jī)振動信號特征分析及故障診斷方法研究［D］.大連:大連理工大學(xué)機(jī)械工程學(xué)院，2010.