長(zhǎng)距離管道化輸送用隔膜泵活塞桿的穩(wěn)定性分析與計(jì)算

馮智睿

（中國(guó)有色（沈陽）泵業(yè)有限公司，沈陽 110141）

1 引言

長(zhǎng)距離管道化輸送系統(tǒng)要求整條輸送管道保證連續(xù)不間斷運(yùn)行，通常單條管道設(shè)置多個(gè)泵站，對(duì)每臺(tái)隔膜泵的壓力和流量要求也較高。活塞桿作為連接隔膜泵動(dòng)力端和液力端的重要零件，將曲軸旋轉(zhuǎn)產(chǎn)生的動(dòng)力傳遞到活塞，進(jìn)行水平推動(dòng)完成料漿的一次輸送，其可靠的質(zhì)量是保證隔膜泵連續(xù)運(yùn)行的重要條件［1］。據(jù)不完全統(tǒng)計(jì)，往復(fù)泵的典型機(jī)械事故中，活塞桿的疲勞斷裂居于首位［2］，因此，在活塞桿的設(shè)計(jì)過程中需進(jìn)行活塞桿與活塞連接的強(qiáng)度校核及穩(wěn)定性校核。

當(dāng)結(jié)構(gòu)所受載荷達(dá)到某一值時(shí)，若增加一微小的增量，則結(jié)構(gòu)的平衡位形將發(fā)生很大的改變，這種情況叫做結(jié)構(gòu)失穩(wěn)或屈曲。屈曲分析是一種用于確定結(jié)構(gòu)開始變得不穩(wěn)定時(shí)的臨界載荷和屈曲模態(tài)形狀的穩(wěn)定性校核技術(shù)［3］。

圖1 活塞桿結(jié)構(gòu)圖

2 ANSYS 軟件屈曲分析方法介紹

有限元分析軟件ANSYS 提供了特征值（線性）屈曲分析和非線性屈曲分析兩種方法［3，4］。

特征值屈曲分析屬于線性分析，它對(duì)結(jié)構(gòu)臨界失穩(wěn)力的預(yù)測(cè)往往要高于結(jié)構(gòu)實(shí)際的臨界失穩(wěn)力，但特征值屈曲分析能夠預(yù)測(cè)臨界失穩(wěn)力的大致所在，可以為非線性屈曲分析時(shí)所加力的大小提供依據(jù)。特征值屈曲分析所預(yù)測(cè)的結(jié)果我們只取最小的第一階，所以特征值臨界失穩(wěn)力的大小應(yīng)為F＝實(shí)際施加力×第一階頻率。

非線性屈曲分析是在大變形效應(yīng)開關(guān)打開的情況下的一種非線性靜力學(xué)分析，因?yàn)橹挥屑せ钤撨x項(xiàng)才能使得幾何剛度矩陣保存下來，該分析過程一直進(jìn)行到結(jié)構(gòu)的極限荷載或最大荷載。非線性屈曲分析的方法是：逐步地施加一個(gè)恒定的荷載增量，直到解開始發(fā)散為止。尤其重要的是，要用一個(gè)足夠小的荷載增量來使荷載達(dá)到預(yù)期的臨界屈曲荷載。若荷載增量太大，則屈曲分析所得到的屈曲荷載就可能不準(zhǔn)確，在這種情況下可以打開自動(dòng)時(shí)間步長(zhǎng)功能，ANSYS 程序?qū)⒆詣?dòng)尋找屈曲荷載。

3 桿件臨界載荷解析解計(jì)算

工程中桿件基本分為細(xì)長(zhǎng)桿、中長(zhǎng)桿和粗短桿三類，不同受壓桿件的失效形式為：細(xì)長(zhǎng)桿（λ≥λp）發(fā)生彈性失穩(wěn)（屈曲）；中長(zhǎng)桿（λo≤λ≤λp）發(fā)生彈塑性失穩(wěn)（屈曲）；粗短桿（λ≤λo）發(fā)生強(qiáng)度破壞。

其中對(duì)于碳鋼及部分合金鋼，λo=40，λp＝100。

桿件穩(wěn)定性計(jì)算方法見表1。

表1 壓桿穩(wěn)定性校核公式

屈曲計(jì)算多用于細(xì)長(zhǎng)桿件，但實(shí)際工程中多數(shù)活塞桿的長(zhǎng)徑比較小，因此，為驗(yàn)證解析解計(jì)算與數(shù)值解計(jì)算的相互對(duì)照關(guān)系，暫取柔度＞100的受壓桿件進(jìn)行分析。

細(xì)長(zhǎng)桿承受壓力，兩端鉸支，桿的截面為矩形，高度h和寬度b 均為30mm，桿的長(zhǎng)度l=2000mm，材料為Q235A，彈性模量E=2×105MPa。

桿橫截面的面積：A=bh=900mm2

式中：μ為受壓桿的長(zhǎng)度系數(shù)，與桿的邊界條件有關(guān)，設(shè)基本情況為兩端鉸支且μ=1，則兩端固定時(shí)μ=0.5，一端固定一端鉸支時(shí)μ=0.7。

以該截面尺寸的桿件解析解計(jì)算過程為例，均以歐拉公式校核桿件的穩(wěn)定性，建立Excel 計(jì)算表格，其余截面尺寸的計(jì)算結(jié)果見表2。

4 ANSYS 特征值求解

建立APDL 程序，對(duì)不同截面形狀及尺寸的桿件分別進(jìn)行特征值屈曲分析，得出該桿件的屈曲特征值，由于設(shè)定施加的載荷均為單位載荷，因而ANSYS 求解所得屈曲特征值即為桿件的臨界載荷。不同桿件的相應(yīng)參數(shù)及屈曲特征值見表2。

表2 不同桿件的參數(shù)及ANSYS 屈曲臨界載荷

由表2 可知，對(duì)于截面高度和寬度均為30mm、長(zhǎng)度為2000mm的桿件，其歐拉公式所得的臨界壓力與ANSYS特征值計(jì)算所得臨界載荷相同，其余截面尺寸桿件的結(jié)果相差也不大，但通過柔度項(xiàng)可知，ANSYS 特征值計(jì)算以細(xì)長(zhǎng)桿的歐拉公式為底層運(yùn)算依據(jù)，因此，對(duì)于中小柔度壓桿的穩(wěn)定性計(jì)算以泰特麥爾-雅辛斯基公式及強(qiáng)度校核公式為準(zhǔn)，ANSYS 屈曲分析校核為輔。

5 隔膜泵活塞桿的穩(wěn)定性校核及特征值求解

對(duì)部分實(shí)際生產(chǎn)的隔膜泵活塞桿零件分別進(jìn)行壓桿穩(wěn)定性校核及ANSYS 特征值分析校核。以表3 案例1中活塞桿為例進(jìn)行穩(wěn)定性校核計(jì)算，活塞桿材料為42CrMoA，其機(jī)械性能σb=800MPa，σs=550MPa。

活塞桿估算等截面直徑d=135mm，活塞桿計(jì)算長(zhǎng)度l=880mm

表3 不同型號(hào)隔膜泵活塞桿參數(shù)及安全系數(shù)計(jì)算表

其中，Pmax-最大活塞力，N，此處Pmax=1200000N；［σb］-許用抗拉極限，MPa，此處?。郐襜］＝σs/2=275MPa。

其中，計(jì)算案例1～3為三缸單作用隔膜泵活塞桿，其安全系數(shù)為按表1中小柔度壓桿穩(wěn)定性校核公式計(jì)算的結(jié)果，計(jì)算案例4、5為雙缸雙作用隔膜泵活塞桿，其中案例4 按小柔度計(jì)算，案例5 按中柔度計(jì)算，對(duì)于碳鋼及合金鋼材料，泰特麥爾-雅辛斯基公式中a=1000，b=5.4［5］。

由表2 可知，以歐拉公式計(jì)算細(xì)長(zhǎng)桿件的臨界載荷與ANSYS 特征值屈曲分析所得結(jié)果一致，但是對(duì)于中小柔度桿件（活塞桿多為中小柔度桿件），不能應(yīng)用歐拉公式，折減彈性模量理論是研究小柔度壓桿臨界應(yīng)力的理論之一，即表1 所示壓桿穩(wěn)定性校核中的泰特麥爾-雅辛斯基公式，但在壓桿柔度很小時(shí)，該理論仍存在不足，故取σCr=σS，即表1 所示小柔度壓桿校核公式。而ANSYS分析所得安全系數(shù)為歐拉公式所得結(jié)果，安全系數(shù)結(jié)果多數(shù)偏大。

6 結(jié)論

由于隔膜泵活塞桿多為中小柔度桿件，很少有細(xì)長(zhǎng)桿，因此對(duì)于中柔度活塞桿可以按泰特麥爾-雅辛斯基公式進(jìn)行穩(wěn)定性校核，而對(duì)于小柔度的活塞桿，應(yīng)以靜強(qiáng)度、疲勞強(qiáng)度及小柔度桿件的穩(wěn)定性校核為準(zhǔn)，不必進(jìn)行有限元屈曲分析。

［1］凌學(xué)勤.往復(fù)式活塞隔膜泵與尾礦存儲(chǔ)技術(shù)［J］.金屬礦山，2010（4）：112-115.

［2］陳鵬霏，孫志禮.壓縮機(jī)活塞桿虛擬可靠性試驗(yàn)方法研究［J］.兵工學(xué)報(bào)，2010，31（9）：1235-1240.

［3］高耀東，郭喜平.ANSYS 機(jī)械工程應(yīng)用25 例［M］.北京：電子工業(yè)出版社，2007.

［4］顏高社，劉光明，蔣桂林.基于ANSYS的活塞桿校核與模態(tài)分析［J］.機(jī)械工程師，2010（9）：119-120.

［5］往復(fù)泵設(shè)計(jì)編寫組.往復(fù)泵設(shè)計(jì)［M］.北京：機(jī)械工業(yè)出版社，1987.