鉆井泵活塞缸套密封系統密封性能流-固耦合分析*

侯勇俊 都 磊 張芳芳 蒲宗珉 席建秋

(1.西南石油大學機電工程學院 四川成都 610500；2.宏華集團 四川成都 610036；3.寶雞石油機械有限公司成都裝備制造公司 四川成都 610052)

我國對石油天然氣的開采隨鉆井工藝技術的進步而逐漸向深地層發展。隨著鉆井深度增加，地層情況越來越復雜，對鉆井泵的工作壓力要求也越來越高。在鉆井泵液力端中，活塞缸套密封性能將直接影響到整個鉆井生產的效率[1]。由于高壓作業時，活塞的使用壽命和密封性能均急速下降，所以研究活塞的力學特性，不僅能分析出活塞失效的原因，還能為活塞后期的優化設計提供相應的指導依據。近年來，國內很多學者對活塞的使用性能進行了深入研究。王曉偉等[2]分析了活塞缸套的失效形式，并提出了合理使用缸套、活塞的方法。王新怡[3]通過有限元軟件對活塞皮碗及連桿進行了力學分析，并針對其薄弱環節提出了改進措施。文獻[4-6]通過觀察生活中動植物形態，設計了多種仿生形活塞，并通過有限元軟件ANSYS建立了仿生活塞三維軸對稱結構，研究了仿生活塞受力特性。苑惠娟等[7]通過實驗分析，對活塞的失效原因進行了歸納總結，并提出了提高活塞壽命的建議。

現有的相關研究大多都忽略了鉆井泵工作過程中活塞缸套與鉆井液的流-固耦合作用，得出的結論與實際的受力特性有一定的差別。本文作者以某2000HP型鉆井泵活塞為研究對象，活塞缸套示意圖如圖1所示，基于Workbench下FSI流-固耦合分析模塊，研究活塞在吸排液過程中不同工作壓力下活塞缸套與鉆井液流-固耦合作用的力學特征，同時探究過盈量和皮碗唇角等參數對活塞缸套密封性能的影響。

圖1 2000HP型活塞缸套示意

1 數學模型

1.1 丁腈橡膠本構方程

活塞皮碗材料為丁腈橡膠，是一種高度非線性、不可壓縮的超彈性體，一般的彈塑性理論不能準確地描述其力學特性[8]。當前針對橡膠材料的模型多樣，文中選用Mooey-Rivlin模型來描述橡膠的力學特性[9-11]。超彈性本構模型通過應變能密度函數來定義，關系如下：

(1)

式中：W為應變能密度；Cij為材料常數；I1、I2、I3為變形張量；D為材料不可壓縮參數，J為彈性體積比。

通過主伸長比λ1、λ2、λ3可得到I1、I2、I3，具體如下：

(2)

在橡膠材料單軸拉伸實驗過程中，主伸長比λ1、λ2、λ3的關系如下：

(3)

式中：λμ為沿載荷方向上的主拉伸比；λ2、λ3為垂直于載荷方向上的主伸長比。

文中皮碗橡膠材料的本構模型采用兩參數Mooey-Rivlin模型，則有：

W=C10(I1-3)+C01(I2-3)

(4)

式中的C10、C01通過橡膠材料的拉伸實驗數據及ANSYS仿真軟件擬合得到，最終得到C10=2.688 MPa，C01=4.021 MPa。

1.2 流-固耦合有限元理論

流-固耦合除遵循最基本的守恒原則，滿足流體控制方程(質量、動量守恒方程)和結構域控制方程，還需在流固耦合交界面處，應滿足流體和固體應力(τ)、位移(d)、熱流量(q)、溫度(T)等變量相等或守恒。但文中未考慮溫度對活塞缸套的影響，因此在未考慮能量傳遞情況下的流-固耦合在交界處的控制方程應滿足以下條件[12-13]：

(5)

式中：d為耦合邊界位移；τ為應力邊界(下標f、s分別代表流體和固體)。

2 模型及設置

2.1 幾何模型及參數

根據2000HP型鉆井泵活塞缸套尺寸，在三維建模軟件中建立活塞缸套的幾何模型，其中皮碗唇角θ=20°，過盈量λ=1.5 mm，缸套內徑為114 mm。由于在鉆井液的壓力作用下，活塞缸套密封系統中卡簧和壓板只是起到傳遞載荷的作用，對密封性能的影響較小，為避免模型過于復雜導致結果不收斂，將結構中非關鍵部位進行簡化，其最終模型結構如圖2所示。

圖2 幾何模型

皮碗材料為丁腈橡膠，墊圈材料為MC尼龍，其余結構材料采用結構鋼，具體材料參數如表1。

表1 材料參數

2.2 網格劃分及無關性分析

根據建立的活塞缸套模型提取的流體域模型采用四面體網格劃分，并合理地分布網格的疏密，如圖3(a)所示；固體域模型網格劃分采用以四面體網格為主的方式，并對皮碗、墊圈、缸套內表面等重要的部位進行網格細化，如圖3(b)所示。

圖3 幾何模型網格

為保證仿真結果的準確性和穩定性，需驗證仿真模型的網格無關性。采用相同的物理參數和邊界條件對整個耦合區域進行網格無關性驗證，選取5組不同網格，計算得到皮碗在最大鉆井液壓力51.7 MPa作用下皮碗最大變形量，如表2所示。根據表2，最終選擇網格節點為29 518，其中皮碗網格為2 574，此時最大變形量變化較小，滿足計算精度和時間成本要求。

表2 網格無關性驗證

2.3 邊界條件

由于分析的是活塞在排吸液過程中，鉆井液壓力變化過程中不同壓力下的活塞缸套的力學特征，而缸套內鉆井液壓力變化過程中，在不同壓力值的各個瞬時狀態下，活塞在缸套中都是處于平衡狀態,故采用單向流-固耦合方法，流體計算采用穩態分析，結構計算采用靜力學分析。

穩態流體分析：選定鉆井液作為流體域材料，活塞排液過程中無入口邊界條件，出口選用壓力出口，由參考文獻[14]通過理論計算和實驗測試得到鉆井泵缸套內壓力變化規律，如圖4所示。根據圖4分別選取工作壓力為0、10、20、30、40和51.7 MPa進行計算；選用k-ω湍流模型[15-16]；采用分離式求解器，壓力速度耦合方式選用simple算法，采用二階迎風格式進行求解；通過設置計算殘差等參數來判斷計算是否收斂。初始化后，迭代次數設置為2 000次，計算時間為1 s；計算完成后，將鉆井液產生載荷導入固體域作為活塞受到的外載。

圖4 缸套內壓力變化曲線

穩態靜力學分析：為得到正確的活塞缸套受力情況，文中根據鉆井泵活塞缸套實際運行狀態來設置接觸和約束。由于皮碗和墊圈在鉆井液壓力作用下會沿軸向發生移動變形，所以將皮碗和缸套接觸面、墊圈和缸套接觸面都設置為摩擦接觸(Frictional)，摩擦因數為0.124[17]；由于討論的是活塞某一時刻靜態受力特性，所以鋼芯與缸套之間設置為固定約束(Fix Support)；活塞皮碗和鋼芯接觸面、墊圈和鋼芯接觸面之間沒有間隙，同時又允許相互滑動，所以將其接觸類型設置為不分離接觸(No Separate)；活塞皮碗與墊圈是整體硫化黏結而成，皮碗與墊圈之間的接觸面之間不允許相互滑動，所以將其接觸類型設置為綁定接觸(Bonded)。定義缸套固定不動，設置皮碗端面和缸套內表面為流-固耦合面。

3 結果分析與討論

3.1 一個周期內皮碗受力變化趨勢

根據缸套內壓力變化規律，分別選取缸套內壓力節點0、10、20、30、40和51.7 MPa為耦合計算壓力，再根據凸輪轉角所對應的缸套內壓力繪制一個周期內各壓力下皮碗與缸套接觸面之間的最大等效應力、接觸壓力和摩擦應力曲線。

圖5所示是皮碗與缸套接觸面之間的最大等效應力、接觸壓力和摩擦應力隨凸輪轉角的變化曲線?？梢钥闯?，其變化規律同缸套內壓力變化曲線相似，在凸輪轉角0°～22°范圍內，也就是排液階段，皮碗與缸套間的最大等效應力、接觸壓力以及摩擦應力都隨缸套內壓力增大而增大；在凸輪轉角為22°時，壓力達到最大51.7 MPa，最大等效應力、接觸壓力和摩擦應力也達到最大；在凸輪轉角為180°～202°時，此時進入吸液階段，皮碗受力也隨缸套內壓力減小而減??；在凸輪轉角為202°時，缸套內壓力減小為0后，皮碗與缸套間的受力也減小到最小。

圖5 一個周期內皮碗最大受力曲線

3.2 最大鉆井液壓力下活塞缸套受力情況

通過分析發現工作壓力越大皮碗受力也越大，因此研究活塞缸套的密封性能時，需重點分析活塞缸套在51.7 MPa最大鉆井液作用下的受力情況。通過計算得到在51.7 MPa鉆井液作用下，缸套與活塞接觸內表面、墊圈表面、活塞皮碗與缸套接觸面的等效應力、接觸壓力和摩擦應力等。

缸套與活塞接觸內表面等效應力分布如圖6所示?？梢钥闯?，缸套內表面等效應力分布分為3個區域，活塞左端既沒有與活塞接觸，也不受鉆井液壓力作用，等效應力較小，為40～45 MPa；缸套與活塞接觸處即為等效應力最大的區域，其變化規律與活塞皮碗尺寸的變化規律相同，由皮碗根部接觸處至唇部逐漸增大，為60～109 MPa，在皮碗的唇部接觸處等效應力最大；在活塞右端，僅有鉆井液壓力作用，等效應力分布較均勻，為55～60 MPa?？傮w來看，缸套內表面最大等效應力小于結構鋼的屈服應力，表明缸套在51.7 MPa鉆井液作用下，不會產生明顯的變形，能夠正常工作。

圖6 缸套內表面等效應力分布

圖7所示為墊圈表面等效應力分布?？梢钥闯?，由于皮碗屬于非線性材料，所以墊圈的等效應力為非均勻分布。墊圈的等效應力主要集中在與皮碗接觸的內圈表面，這是由于通過皮碗傳遞的鉆井液壓力首先作用在墊圈與皮碗接觸的表面，同時皮碗受壓會使其徑向體積增加，導致與皮碗接觸的墊圈內圈受到的等效應力大于外圈，為100～138 MPa。但由于墊圈材料為MC尼龍，綜合性能良好，墊圈也不易失效。

圖7 墊圈表面等效應力分布

活塞皮碗，相較于缸套和墊圈是活塞密封中最易失效的部位，也是文中研究的重點對象。皮碗上等效應力大小反映了皮碗表面各主應力的差值。當等效應力增大時，會加劇皮碗橡膠材料的松弛程度，降低其剛度，因而皮碗容易產生裂紋。由于橡膠材料的松弛程度和裂紋對活塞的密封性能有很大的影響，所以等效應力可間接反映活塞的密封性能[19]。圖8(a)所示為活塞皮碗與缸套接觸面的等效應力云圖，圖8(b)所示為皮碗沿軸向取一定間隔點繪制的等效應力曲線?？梢钥闯?，由于活塞皮碗是非線性材料，其等效應力是非均勻分布的。在皮碗的唇部和根部與墊圈接觸的位置等效應力值較大，且唇部等效應力大于根部，為9～12 MPa。這是由于活塞在唇口處為過盈配合，在壓力作用下橡膠產生形變使活塞達到密封效果；由于活塞皮碗直徑從根部逐漸向唇部增大，其背部等效應力也由根部到唇部逐漸增大，為1～9 MPa；同時皮碗根部會出現擠入墊圈與缸套之間的趨勢，所以皮碗根部與墊圈接觸位置應力較大，為6～8 MPa。

圖8 皮碗等效應力分布

活塞缸套密封為接觸型密封，靠外部壓力作用緊貼在密封面上，密封介質壓力小于彈性體對表面的接觸壓力時就能形成良好的密封，所以接觸壓力的大小一定程度上能反映其密封性能。圖9所示為活塞皮碗與缸套接觸面間接觸壓力云圖和曲線圖?？梢钥闯?，2000HP型往復泵活塞皮碗接觸壓力從皮碗根部至唇部逐漸增大，這正與皮碗尺寸的變化規律相吻合?；钊ね氡砻娼佑|壓力也是非均勻分布的，皮碗唇部過盈量最大，所以此處的接觸壓力也最大，為60～72 MPa，且大于鉆井液壓力值。該仿真結果表明活塞在鉆井液工作壓力作用下密封性良好；皮碗背部接觸壓力值變化較為平緩，為25～45 MPa；而由于皮碗根部會擠入墊圈與缸套之間的間隙中，但根部徑向變形后不足以與缸套接觸，所以此處的接觸壓力為0。

圖9 皮碗接觸壓力分布

圖10所示為活塞皮碗與缸套接觸面間摩擦應力云圖和曲線圖?？梢钥闯?，摩擦應力和接觸壓力的變化規律相似。摩擦應力最大值也出現在皮碗唇部邊緣，為6～7 MPa，這是由于在皮碗唇部過盈量較大，皮碗唇部體積會被擠壓，促使唇部邊緣橡膠在有限距離內徑向體積增大，而周圍的直徑會相應減小，導致皮碗唇部向其中心彎曲，邊緣的摩擦應力變大，也在一定程度上反映了活塞皮碗存在向心效應。皮碗唇部由于存在向心效應，鉆井液中磨粒會進入缸套與皮碗的接觸面，可能發生磨損失效。

圖10 皮碗摩擦應力分布

3.3 驗證

為驗證數值模擬方法的可行性，將圖9中的仿真模擬結果與文獻[18]實驗結果進行對比，發現接觸壓力模擬結果與實驗結果變化規律基本吻合。由于模型尺寸和壓力大小有一定差別，所以實驗結果與仿真結果有一定的誤差，但總體的變化趨勢是一致的，表明數值模擬方法的可行。

4 皮碗尺寸對接觸壓力的影響

接觸壓力作為判斷密封性能的重要指標，接觸壓力越大，活塞能密封的介質壓力越大。圖11和圖12分別示出了活塞皮碗在不同唇部過盈量λ和不同唇角θ下皮碗與缸套之間接觸壓力沿軸向的變化情況。

在51.7 MPa鉆井液壓力作用下，活塞唇部過盈量λ從0.5 mm增大到2.0 mm時，活塞皮碗與缸套接觸面間接觸壓力的變化情況如圖11所示?？梢钥闯?，當過盈量λ≤1.0 mm時，皮碗唇部接觸壓力從根部至唇部逐漸增大；過盈量λ>1.0 mm時，由于唇部過盈量大，變形量也更大，橡膠的松弛程度也會加劇，唇部橡膠將向中心彎曲，產生向心效應，導致唇部邊緣后部橡膠在有限距離內徑向體積增大，唇部邊緣的接觸壓力也大于唇部接觸壓力。但總的來說，過盈量λ≤1.0 mm時，沿軸向的接觸壓力變化不大，對密封性能影響也不大；當λ>1.0 mm時，由于皮碗唇部存在向心效應，發生磨損失效的可能性增大，皮碗缸套的密封性能也會受到影響。

圖11 不同過盈量時的接觸壓力曲線

圖12顯示了在51.7 MPa鉆井液壓力作用下皮碗唇角θ為15°～30°時，活塞皮碗與缸套接觸面間接觸壓力變化情況?？梢钥闯?，唇角的變化會影響向心效應，當唇角θ≤15°時皮碗唇部不會產生向心效應，密封性能最好；當θ>15°時，皮碗唇部橡膠將向中心彎曲，產生向心效應，進而影響皮碗的密封性能。

圖12 不同唇角時的接觸壓力曲線

5 結論

(1)一個周期內皮碗與缸套接觸面之間的最大等效應力、接觸壓力和摩擦應力曲線變化規律同缸套內壓力變化曲線相似，都隨鉆井液壓力增大而增大，隨壓力的減小而減小。

(2)在最大鉆井液壓力作用下，缸套內表面和墊圈表面與皮碗接觸部位等效應力都較大。后期改進設計時應重點考慮強化受力較大的部位。

(3)在最大鉆井液壓力作用下，皮碗受到應力都主要集中在活塞皮碗唇部和根部與墊圈接觸的位置，也是皮碗易發生失效的2個部位，所以加強皮碗這2個部位將有助于活塞壽命的提高。

(4)當皮碗與缸套之間過盈量λ≤1.0 mm時，接觸壓力沿皮碗軸向分布規律相似，對密封性能影響不大；當過盈量λ>1.0 mm時，會導致皮碗唇部的接觸壓力小于皮碗唇部邊緣處，產生向心效應，發生磨損失效的可能性增大，皮碗缸套的密封性能也會受到影響。

(5)當皮碗的唇角θ≤15°時，皮碗唇部不會產生向心效應，密封性能最好；當唇角θ>15°時，皮碗唇部橡膠將向中心彎曲，產生向心效應，進而影響皮碗的密封性能。