吳高峰,周思柱,王峻喬,吳漢川

(1.長(zhǎng)江大學(xué)機(jī)械工程學(xué)院,湖北荊州434023; 2.中國(guó)石化集團(tuán)江漢石油管理局第四石油機(jī)械廠,湖北荊州434023) *

基于Ansys/Ls-Dyna的壓裂泵泵閥強(qiáng)度分析

吳高峰1,周思柱1,王峻喬2,吳漢川2

(1.長(zhǎng)江大學(xué)機(jī)械工程學(xué)院,湖北荊州434023; 2.中國(guó)石化集團(tuán)江漢石油管理局第四石油機(jī)械廠,湖北荊州434023)*

泵閥是壓裂泵中最重要的易損部件之一,其強(qiáng)度關(guān)系到壓裂泵的工作特性。運(yùn)用Ansys有限元分析軟件中Ls-Dyna模塊,模擬100 MPa高壓環(huán)境下閥盤(pán)以一定速度沖擊閥座的整個(gè)過(guò)程,對(duì)泵閥進(jìn)行靜力以及顯示動(dòng)力下應(yīng)力、應(yīng)變分析,找出泵閥的失效原因,為改進(jìn)泵閥提供了理論依據(jù)。

壓裂泵;泵閥;Ansys/Ls-Dyna;靜力分析;動(dòng)力分析

高壓是壓裂泵不同于鉆井泵的重要方面。某地區(qū)4 000 m井深要求壓裂泵工作壓力100～140 MPa,導(dǎo)致泵閥難以完成一次作業(yè)(2.5 h)就損壞[1-3]。本文采用Ansys/Ls-Dyna軟件對(duì)泵閥關(guān)閉過(guò)程進(jìn)行靜力和顯示動(dòng)力分析,分析了泵閥在高壓條件下的承載能力,找出失效原因,為提高泵閥使用壽命提供了理論基礎(chǔ)。

1 靜力分析

1.1 結(jié)構(gòu)及受力

壓裂泵泵閥結(jié)構(gòu)如圖1。當(dāng)泵閥關(guān)閉后,閥盤(pán)受到重力 G、彈簧拉力 Ft和液體的壓力差作用力F1-F2,由力學(xué)平衡原理可知

式中,FFP為閥盤(pán)受力,N。

1.2 建立分析模型

泵閥在關(guān)閉的狀態(tài)之下,閥盤(pán)工作面與閥座工作面接觸,是典型的接觸問(wèn)題。因此,泵閥工作面的應(yīng)力、應(yīng)變是分析重點(diǎn)。為使分析結(jié)論更接近實(shí)際的惡劣工況,在此沒(méi)有考慮橡膠密封圈。閥盤(pán)在不考慮導(dǎo)向爪時(shí),其結(jié)構(gòu)和受載均是對(duì)稱的,所以取泵閥的 1/4結(jié)構(gòu)進(jìn)行分析。建立的分析模型如圖2～3。材料參數(shù)如表1。

圖1 壓裂泵泵閥結(jié)構(gòu)

表1 某壓裂泵泵閥材料參數(shù)

圖2 閥盤(pán)和閥座有限元模型

圖3 閥盤(pán)和閥座接觸定義

1.3 施加載荷和邊界條件及計(jì)算結(jié)果

鑒于所建泵閥分析模型結(jié)構(gòu)的特殊性,對(duì)閥盤(pán)截面位置進(jìn)行徑向約束,閥座與箱體接觸,y方向位移為零。根據(jù)前面靜力分析理論可知,閥盤(pán)上部受液體壓力作用,壓力為100 MPa。閥盤(pán)、閥座選擇Plane42單元[4-6],施加相應(yīng)的邊界條件和載荷后進(jìn)行分析,閥盤(pán)、閥座變形如圖4,應(yīng)力云圖如圖5。

圖4 閥盤(pán)和閥座變形云圖

圖5 閥盤(pán)和閥座應(yīng)力云圖

2 顯示動(dòng)力分析

2.1 理論基礎(chǔ)

Ls-Dyna采用中心差分法來(lái)進(jìn)行時(shí)間積分,在已知1,…,tn時(shí)間步的情況下,求解tn+1時(shí)間步的解,運(yùn)動(dòng)方程為

式中,M為集中質(zhì)量矩陣;¨U為節(jié)點(diǎn)加速度列陣;?U為節(jié)點(diǎn)速度列陣;P(t)為外力向量列陣;H(tn)為砂漏阻力;Fint(t)為內(nèi)力矢量,為單元內(nèi)力與接觸內(nèi)力之和。Fint表達(dá)式為

單元內(nèi)力由當(dāng)前構(gòu)型的應(yīng)力場(chǎng)求得。

這樣可以求得在 tn+1時(shí)刻的位移,更新 tn時(shí)刻的系統(tǒng)幾何構(gòu)型,得到tn+1時(shí)刻系統(tǒng)新的幾何構(gòu)型。由于采用集中質(zhì)量矩陣M,運(yùn)動(dòng)方程的求解是非耦合的,不需要組集成總體剛度矩陣,并采用中心單點(diǎn)積分,因此極大節(jié)省存儲(chǔ)空間和求解時(shí)間[7-8]。

2.2 建立模型

導(dǎo)向爪在關(guān)閉過(guò)程中對(duì)閥盤(pán)、閥座強(qiáng)度的影響較小,同時(shí)橡膠圈的緩沖作用增強(qiáng)了泵閥的承受能力。采用 solid164單元對(duì)泵閥進(jìn)行實(shí)體建模,如圖6。

圖6 顯示動(dòng)力分析模型

2.3 載荷和邊界條件及計(jì)算結(jié)果

泵閥在關(guān)閉時(shí),受到自身重力、彈簧拉力及較小的流體壓力差的作用,閥盤(pán)以一定的速度(取150 m/s)向閥座運(yùn)動(dòng),與閥座發(fā)生碰撞。閥座下端面軸向位移為零,閥座1/4截面徑向位移為零,利用Ls-Dyna求解器計(jì)算出的閥盤(pán)碰撞閥座的應(yīng)力變化,在碰撞后2×10-3ms時(shí)刻的應(yīng)力云圖如圖7,最大應(yīng)力達(dá)900 MPa,閥座與閥盤(pán)下錐面應(yīng)力為750 MPa。

圖7 閥盤(pán)碰撞閥座2×10-3ms時(shí)刻的應(yīng)力云圖

3 綜合分析

綜合靜力、動(dòng)力分析可知,泵閥在關(guān)閉過(guò)程中,閥盤(pán)對(duì)閥座的沖擊對(duì)泵閥強(qiáng)度的影響遠(yuǎn)大于靜力作用。在靜壓作用下閥盤(pán)關(guān)閉后泵閥最大應(yīng)力為393 MPa,在泵閥材料的屈服極限內(nèi)。在泵閥沖擊過(guò)程中,錐面下端出現(xiàn)的最大局部應(yīng)力為 750 MPa,小于材料的屈服極限,是安全的。錐面上端出現(xiàn)最大局部應(yīng)力達(dá)到900 MPa,即錐面處出現(xiàn)的應(yīng)力集中更易使得橡膠圈磨損,這一點(diǎn)與實(shí)際上橡膠圈出現(xiàn)嚴(yán)重的變形相一致,同時(shí)也是泵閥橡膠圈失效的原因。

4 結(jié)論

1) 利用Ansys/Ls-Dyna軟件對(duì)高壓壓裂泵泵閥進(jìn)行靜態(tài)、動(dòng)態(tài)的仿真計(jì)算,求解出泵閥在靜力作用和沖擊過(guò)程中出現(xiàn)的最大局部應(yīng)力。

2) 壓裂泵泵閥錐面上端局部應(yīng)力過(guò)大是導(dǎo)致泵閥橡膠密封圈過(guò)度磨損以至泵閥失效的原因。

3) 增強(qiáng)橡膠密封圈的性能,選擇強(qiáng)度更高的材料,設(shè)計(jì)新型的泵閥結(jié)構(gòu)是提高壓裂泵泵閥使用壽命的有效方法。

[1] 羅安賢,齊樹(shù)敏,楊啟明.優(yōu)質(zhì)壓裂泵泵閥分析[J].西南石油學(xué)院學(xué)報(bào),1992(3):60-66.

[2] 周錫容,羅安賢,王子瑜,等.壓裂泵泵閥失效分析[J].石油機(jī)械,1990(11):22-24.

[3] 楊國(guó)安,張 冬,黃 聰.鉆井泵閥的沖擊特性分析[J].振動(dòng)與沖擊,2008(12):18-22.

[4] 楊國(guó)安,喬 艦,宋 征,等.鉆井泵閥壽命分析[J].振動(dòng)與沖擊,2010(1):58-61.

[5] 王中輝,戴靜君,孟 波,等.往復(fù)泵泵閥強(qiáng)度計(jì)算[J].石油礦場(chǎng)機(jī)械,2010,39(2):34-37.

[7] 李洪波,劉振龍,周天明,等.F-1600型泥漿泵閥座的接觸分析[J].石油礦場(chǎng)機(jī)械,2010,39(5):26-29.

[8] 往復(fù)泵設(shè)計(jì)編寫(xiě)組.往復(fù)泵設(shè)計(jì)[M].北京:機(jī)械工業(yè)出版社,1987.

[9] 楊啟明,周錫容,譙 英.新型壓裂泵泵閥材料的耐磨損性能研究[J].西南石油學(xué)院學(xué)報(bào),2001(6):73-76.

Strength Analysis of Fracturing Pump Valve Based on Ansys/Ls-Dyna

WU Gao-feng1,ZHOU Si-zhu1,WANGJun-qiao2,WU Han-chuan2
(1.College of Mechanical Engineering,Yangtze University,J ingzhou434023,China; 2.S j Petroleum Machinery Co.,Sinopec,J ingzhou434023,China)

Fracturing pump valve was the most important wearing parts;the strength of pump valve was directly related to working characteristics of fracturing pump valve.With the Ls-Dyna module of Ansys finite element analysis software,the whole impact process in which the valve disc was impacting the seat in 100Mpa high pressure environment was simulated,the static and dynamic stress-strain of pump valve was analyzed and the failure cause of fracturing pump valve was found,then theory foundation was provided for prolonging the service life of fracturing pump valve.

fracturing pump;valve;Ansys/Ls-Dyna;static analysis;dynamic analysis

1001-3482(2011)05-0050-03

TE934.202

A

2010-08-03

國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)項(xiàng)目“2500HP大型數(shù)控成套壓裂裝備研制”(2009AA063601)資助

吳高峰(1986-),男,江西萍鄉(xiāng)人,碩士研究生,研究方向是機(jī)械結(jié)構(gòu)與強(qiáng)度。