考慮穩態液動力的錐閥閥芯可靠性分析

王海芳, 褚天爭, 王繼強, 張 恒

(東北大學秦皇島分?？刂乒こ虒W院,河北秦皇島 066004)

考慮穩態液動力的錐閥閥芯可靠性分析

王海芳, 褚天爭, 王繼強, 張 恒

(東北大學秦皇島分校控制工程學院,河北秦皇島 066004)

以單向錐閥閥芯為例,應用蒙特卡羅數值模擬法,考慮了錐閥閥芯在工作過程中所受到的穩態液動力,將其作為面載荷并與其他工作因素相結合,利用ANSYS建模及其PDS(Probabilistic Design System)模塊進行可靠性分析,獲得了錐閥閥芯的應力分布圖、可靠度以及靈敏度圖.在置信度為95%的情形下,計算得到閥芯可靠度為99.82%.分析結果表明:穩態液動力、工作壓力以及最大截面圓半徑是影響可靠性的最主要因素.為其結構設計優化提供了定性及定量的依據,同時證明了穩態液動力對錐閥閥芯可靠性有較大影響,在設計過程中不可忽略.

可靠性; 液壓閥; 靈敏度; 穩態液動力

0 引言

單向錐閥作為液壓系統中的常用元件,應用廣泛.錐閥閥芯的傳統可靠性設計主要采用常規的安全系數法,形式簡單明了,但也具有一定的經驗性和保守性[1-3].另外,考慮到液壓閥所用材料的強度具有隨機離散性,零件的應力也因尺寸的誤差及表面加工粗糙度的不同而隨機波動,使安全系數法的定量概念變得模糊不清[4-7].而基于數值分析的現代可靠性設計,將非線性的功能函數進行線性化,解決了上述因材料數據及制造誤差帶來的隨機離散化的問題.論文基于現代可靠性設計理論,應用ANSYS/PDS設計分析模塊,將閥芯的材料參數、幾何尺寸、載荷等作為隨機變量參數,并考慮了錐閥閥芯所受的穩態液動力,將其以面載荷的形式進行加載,選擇蒙特卡羅法,應用應力-強度干涉理論,對閥芯的可靠性進行計算,計算過程為閥芯的可靠性定量設計提供了分析基礎.

1 蒙特卡羅可靠性分析法

1.1 蒙特卡羅法基本原理

機械可靠性分析和設計的主要目標是求解研究對象的可靠度R[2]：

(1)

式中：fX(x)為研究對象包含的基本隨機參數向量X=(x1,x2,x3,…,xn)T的聯合概率密度函數,狀態函數Z=g(x)用以表示研究對象的兩種狀態:g(x)≤0為失效狀態,g(x)>0為可靠狀態,g(x)=0是一個n維極限狀態曲面.

蒙特卡羅法是一種用數值模擬來求解與隨機變量有關的實際工程問題的方法,對隨機變量的數值模擬相當于一種“試驗”,所以蒙特卡羅法又稱為統計試驗法.由式(1)顯示的可靠度R是基本隨機變量的聯合概率密度函數fX(x)在可靠域g(x)>0中的多重積分.失效概率Pf可以表示為失效域F的指示函數IF(x)的數學期望形式:

(2)

式中:RMC為蒙特卡羅法數值模擬所獲得的可靠度.蒙特卡羅法數值模擬法進行結構系統可靠度分析的過程如圖1所示.

1.2 蒙特卡羅可靠性靈敏度設計

可靠性靈敏度定義為狀態函數或者是響應量k階矩αkg對基本變量分布參數θx的偏導數?αkg/?θx,可充分反映各設計參數對零部件失效的影響程度,即敏感性,在可靠性優化設計和可靠性穩健設計等方面均有重要的應用[2,8-10].對于相關正態變量情況下的可靠性靈敏度分析,蒙特卡羅直接法以數學期望形式表示的可靠性靈敏度如下:

圖1 蒙特卡羅法進行可靠性分析的流程圖Fig.1 Flow chart of reliability analysis by Monte Carlo method

(3)

(4)

(5)

(6)

2 錐閥的建模與可靠性分析

2.1 問題描述及液動力計算

本文以RVP型單向閥閥芯為例,進行錐閥閥芯的可靠性靈敏度分析計算.圖2為RVP單向閥的整體裝配示意圖.

圖2 RVP單向閥的整體裝配示意圖Fig.2 Schematic diagram of the whole assembly of check valve of RVP

其工作條件為:功率為15 kW,沖次為120次·min-1,流量為145 L·min-1,壓力為20 MPa,通徑為20 mm.錐閥閥芯材料為馬氏體不銹鋼2Cr13,條件屈服極限δ0.2為440 MPa,密度為7.75 g·cm-3.

傳統液壓閥可靠性設計未考慮由于液體動量的改變而使閥芯受到的液動力.液動力分為穩態液動力和瞬態液動力兩種,本文研究了閥芯移動完畢閥口開度固定之后,液流流經閥口時因動量改變而附加作用在閥芯上的穩態液動力[11-14].當錐閥穩定工作時,作用在錐閥閥芯上的力平衡方程為

PA=Ft±Ff+Fs=k(x0+x)±Ff+Fs

(7)

式中:p為錐閥閥口壓力;A為閥芯底部有效作用面積;Ft為錐閥調壓彈簧力;Ff為錐閥閥芯所受的摩擦力;Fs為錐閥閥芯所受的液動力;k為彈簧剛度;x0為彈簧預壓縮量;x為閥的開口量.

一般錐閥在工作時錐閥閥芯所受的摩擦力Ff和閥的開口量x比較小,所以在以下的分析中將其忽略.當錐閥閥芯的錐角為2α(半錐角α一般為12°～20°),液體在壓力p的作用下以流量q流經錐閥,液體流入速度為v1,出流速度為v2.則可沿液流方向列出動量方程[15-16]:

(8)

又因為v1?v2,忽略v1,θ2=α,θ1=0°,帶入后得液動力Fs:

(9)

根據以上分析帶入相應的數據便可以求出液壓錐閥閥芯工作時的穩態液動力,并將其加入到錐閥閥芯可靠性的分析過程中,以此提高閥芯可靠性及靈敏度分析的可信度.

2.2 錐閥閥芯的建模及應力分析

根據錐閥閥芯的工作受力及幾何尺寸情況進行結構簡化,建立錐閥閥芯的三維模型,選用Solid45計算單元進行劃分,Solid45單元用于三維實體結構,單元通過8個節點來定義,每個節點有3個沿著xyz方向平移的自由度[17,18].由于錐閥閥芯的不規則性,本文采用自由網格劃分,有限元劃分網格模型如圖3所示.

將穩態液動力作為面載荷進行加載,施加各類邊界工作條件和其它工作載荷進行靜力分析,得到閥芯的VonMises等效應力圖,如圖4與圖5所示,在圖中可以看到最大節點等效應力的位置為閥芯的根部,并利用命令“get,maxstr,sort,,max”提取最大節點等效應力,其數值為45.65MPa.

2.3 錐閥閥芯的可靠性分析

ANSYS提供的PDS系統又稱概率設計或可靠性設計系統,是用來評估輸入參數的不確定性對于系統輸出的影響行為及其特性的一個模塊[18],在ANSYS中,蒙特卡羅模擬技術可以選擇直接抽樣法或拉丁抽樣法進行處理.文章利用該模塊對單向閥的可靠性進行了分析計算,得到了影響其可靠性定性因素的定量關系,PDS的可靠性分析過程如圖6所示.

圖3 錐閥閥芯有限元模型Fig.3 The finite element model of cone valve spool

圖4 Von Mises等效應力分布圖Fig.4 Von Mises equivalent stress distribution

圖5 Von Mises等效應力分布圖Fig.5 Von Mises equivalent stress distribution

錐閥閥芯在工作中承受的最大應力不允許超過選用材料的條件屈服極限,如果應力超過屈服極限則認為失效.定義失效準則為[19]:δmax≥δ0.2,其中δmax為閥芯在使用過程中出現的最大應力值,其值與圖4和圖5提取的最大節點等效應力相等,δ0.2為材料的條件屈服極限.根據應力-強度干涉模型定義極限狀態函數g(x)=δ0.2-δmax,當g(x)<0時,錐閥閥芯功能失效,其中X為上式中所有不確定量組成的向量[20-21],錐閥閥芯的可靠度就是g(x)>0的概率,最后生成概率分析文件.

圖6 PDS可靠性分析過程Fig.6 Process of PDS reliability analysis

讀入概率分析文件后,根據錐閥的實際工作情況,以錐閥閥芯直徑(D1,D2,D3分別為圖2中錐閥閥芯從左至右的不同橫截圓面的直徑)、壓力載荷(P)、液動力(Fs)、調壓彈簧力(Ft)為隨機輸入變量,定義最大等效應力(δmax)和極限狀態函數g(x)為隨機輸出變量.分析中,假設上述隨機輸入變量均服從正態分布,隨機輸入變量統計值如表1所示.選擇蒙特卡羅法中拉丁抽樣方法來進行概率設計,模擬樣本數為700,即產生700組基本隨機變量的隨機樣本Xj(j=1,2,…,700),對進入可靠域g(x)>0的樣本組Nr進行統計,用安全發生頻率Nr/700來近似代替可靠度R,這樣就可以得可靠度R的近似估計值.

表1 隨機輸入變量統計值

從ANSYS中捕獲到如圖7所示的95%置信度下的閥芯最大等效應力抽樣曲線圖.圖中的中間曲線表示最大等效應力均值的變化曲線；上面曲線表示置信區間上限的變化曲線;下面曲線表示置信區間下限的變化曲線.置信區間的寬度隨著抽樣次數的增加而減小,且最大等效應力均值收斂于21.19 MPa.應力均值變化趨近于水平,說明循環次數足夠多.

2.4 可靠性及靈敏度分析結果

2.4.1 累計分布函數及可靠度

通過累積分布函數圖像可以獲得零件的可靠度或失效概率,其函數上任意一點值等于數據出現在該點的概率值.圖8為置信度為95%的極限函數g(x)的累積分布.并且經過計算分析后,可以得到置信度為95%時極限方程小于0的概率.對于本算例的單向錐閥閥芯在置信度為95%的情形下,極限方程小于零的概率平均為0.18%,則可得出在該工作環境下該錐閥閥芯的可靠度為99.82%.

2.4.2 靈敏度分析

對錐閥閥芯概率靈敏度分析可以得到影響其失效的最主要因素.圖9是錐閥閥芯所受最大等效應力在重要水平為2.5%時的靈敏度圖,從圖中可以看出,對錐閥閥芯最大等效應力影響較大的主要參數是錐閥閥口壓力(P)、錐閥閥芯所受的液動力(Fs)、錐閥閥芯的最大截面圓半徑(D2).因此如要提高錐閥閥芯的可靠度,在設計過程中應嚴格控制這些參數.靜態液動力(Fs)對錐閥閥芯可靠性的靈敏度是負值,閥芯直徑(D2)和流體壓力對其靈敏度是正值,表明錐閥閥芯所受的應力隨錐閥閥芯直徑和流體壓力的增加而增加,隨靜態液動力的增加而減小,同時也表明穩態液動力(Fs)對錐閥閥芯的可靠性有一定影響,在分析設計過程中不能忽略.

圖9 靈敏度分析結果示意圖Fig.9 Theresults of the sensitivity analysis

3 結語

在算例中,單向錐閥閥芯在置信度為95%的情形下,其可靠度計算為99.82%,由可靠性靈敏度的結果分析可知穩態液動力、工作壓力以及最大截面圓半徑是影響錐閥閥芯可靠度的主要因素,驗證了ANSYS/PDS概率設計模塊分析錐閥閥芯可靠性的可行性.分析結果反映了以往在對錐閥閥芯可靠性進行分析時將液動力忽略的做法存在缺陷,為錐閥閥芯可靠性設計提供了新的修改依據.

[1] 黃天成,張炎,王宏麗,等. 基于ANSYS/PDS的抽油機支架可靠性分析[J].制造業自動化. 2011,33(7):119-121.

HUANG Tiancheng,ZHANG Yan,WANG Hongli,et al. Reliable analysis for pumping unit framework based on ANSYS/PDS[J].Manufacturing Automation,2011,33(7): 119-121.

[2] 呂震宇,宋述芳,李洪雙,等. 結構機構可靠性及可靠性靈敏度分析[M].北京: 科學出版社,2009.

LV Zhengyu,SONG Shufang,LI Hongshuang,et al.Strucural reliability and reliability sensitvity analysis[M].Beijing:Scinece Press ,2009.

[3] 杜永英,孫志禮,呂春梅,等.基于ANSYS某型電連接器插拔過程的可靠性分析[J].中國工程機械學報,2015,13(6): 543-549.

DU Yongying,SUN Zhili,LV Chunmei,et al. Reliability analysis on terminal contact parts for electrical connector based on ANSYS[J].Chinese Journal of Construction Machinery,2015,13(6):543-549.

[4] 楊周,張義民. 非正態隨機參數的車輛零件的可靠性靈敏度設計[J].中國工程機械學報,2007,5(3):282-287.

YANG Zhou,ZHANG Yimin. Sensitivity design for reliability of vehicle components based on non-normal random parameters[J].Chinese Journal of Construction Machinery,2007,5(3):282-287.

[5] 李民,舒歌群,梁興雨. 4100QB柴油機機體的形狀優化及靈敏度分析[J].農業工程學報,2010,26(2):160-164.

LI Min,SHU Gequn,LIANG Xinyu. Shape optimization and sensitivity analysis for cylinder block of 4100QB diesel engine[J].Transactions of the Chinese Society of Agricultural Engineering,2010,26(2): 160-164.

[6] ZHANG Yimin,ZHU Lisha,WANG Xingang. Advanced method to estimate reliability-based sensitivity of mechanical components with strongly nonlinear performance function[J].Applied Mathematics and Mechanics,2010,31(10): 1325-1336.

[7] SHIRAISHI F,TOMITA T,IWATA M,et al. A reliable taylor series-based computational method for the calculation of dynamic sensitivities in large-scale metabolic reaction systems:algorithm and software evaluation[J].Mathematical Biosciences,2009,222(2): 73-85.

[8] KUMAWAT H L,SINGH O P,MUNSHI P . Uncertainty evaluation of reliability of safety grade decay heat removal system of indian prototype fast breeder reactor[J].Annals of Nuclear Energy,2013,62: 242-250.

[9] WU Y T,MOHANTY S. Variable screening and ranking using sampling-based sensitivity measures[J].Reliability Engineering and System Safety,2006,91(6): 634-647.

[10] SAMADANI M,BEHBAHANI S,NATARAJ C. A reliability-based manufacturing process planning method for the components of a complex mechatronic system[J].Applied Mathematical Modelling,2014,37(24): 9829-9845.

[11] 陳從桂,林國湘,郭平輝,等.鏟運機動臂抗斷裂強度的可靠性設計[J].中國工程機械學報. 2006,4(2):170-173.

CHEN Conggui,LIN Guoxiang,GUO Pinghui,et al.Reliability design on the strength of loader boom resisting gracture[J].Chinese Journal of Construction Machinery,2006,4(2):170-173.

[12] KIM C,WANG S,CHOI K K.Efficient response surface modeling by using moving least-squares method and sensitivity[J].AIAA Journal,2005,43(11): 2404-2411

[13] WEI D,RAHMAN S.Structural reliability analysis by univariate decomposition and numerical integration[J].Probabilistic Engineering Mechanics,2007(22):27-38.

[14] 田勇,高長銀.液壓與氣壓傳動技術與應用[M].北京:電子工業出版社,2011.

TIAN Yong,GAO Changyin. Hydraulic and pneumatic transmission technology and Application[M].Beijing:Publishing House of Electronics Industry,2011.

[15] 許福玲,陳堯明.液壓與氣壓傳動[M].北京: 機械工業出版社,2007.

XU Fuling,CHEN Yaoming. Hydraulic and pneumatic transmission[M].Beijing:China Machine Press,2007.

[16] LI X,GUAN T,LIU C B,et al. Modeling and simulation of urethra valve of bladder power pump[J].Journal of Mechanics,2014,30(3): 255-264.

[17] MELCHERS R E,AHAMMED M.A fast approximate method for parameter sensitivity estimation in Monte Carlo structural reliability[J].Computers and Structures,2004,82(1): 55-61.

[18] QIN H W,HU H M,YE W,et al. Reliability assessment of the water distributing valve of a hydrostatic sediment corer[J].Journal of Pressure Vessel Technology,Transactions of the ASME,2014,136(6): 061601-061607.

[19] 鞠彥忠.材料力學 [M].第2版.武漢:華中科技大學出版社,2014.

JU Yanzhong.Materials mechanics[M].2nd ed.Wuhan:Huangzhong Univesity of Science and Techonolgy Press,2014.

[20] 施衛東,顏品蘭,蔣小平,等.基于ANSYS/PDS的泵軸可靠性分析[J].排灌機械,2008,26(05):1-4.

SHI Weidong,YAN Pinlan,JIANG Xiaoping,et al. Probabilistic analyse of pump shaft based on ANSYS/PDS[J].Drainage and Irrigation Machinery,2008,26(05): 1-4.

[21] 王金龍,王清明,王偉章. ANSYS12.0有限元分析與范例解析[M].北京: 機械工業出版社,2010.

WANG Jinglong,WANG Pingming,WANG Weizhang.ANSYS12.0 Finite element analysis and example analysis[M].Beijing:China Machine Press,2010.

Reliability analysis on check-cone valve spool based on steady flow force

WANG Hai-fang, CHU Tian-zheng, WANG Ji-qiang, ZHANG Heng

(College of Control Engineering, Northeastern University at Qinhuangdao, Qinhuangdao 066004, China)

With regard to the check-cone valve spools,the surface loading, together with other working conditions, is considered in terms of steady flowforce in the working process. By using Monte Carlo numerical simulation method, the reliability analysis is conducted on spool stress distribution, reliability and sensitivity via ANSYSTMmodeling and PDSTMmodule. In details, the reliability of spoolreaches 99.82%, whilst the confidence level is 95%. Therein, it is found from analysis results that such majorimpact factors as steady flowforce,working pressure and radius of maximum section circle set a reference to thedirectional and quantitative analysis on structural design optimization. In addition, the steady flowforce of check-conevalve spooldelivers significant impacts upon pool reliability.

reliability; hydraulic valve; sensitivity; steady flow force

國家自然科學基金資助項目(51475086);河北省自然科學基金資助項目(E2012407010);秦皇島科技支撐項目(201501B011)

王海芳(1976-),男,副教授,博士.E-mail:hfwang0335@126.com

TH 122

A

1672-5581(2016)05-0426-06