CN112257268A - 一种液压多路阀可靠性分析方法 - Google Patents

Abstract

本发明提供了一种液压多路阀可靠性分析方法，包括：步骤1：对于液压多路阀的典型失效模式建立树状图；步骤2：进行参数查阅；步骤3：根据最小割集算法，算出最小割集；步骤4：绘制出多路阀失效的逻辑框图；步骤5：计算各部件或模块的可靠度并根据树状图计算液压多路阀系统可靠度；步骤6：对于液压多路阀系统，分析其可靠度与平均寿命；步骤7：根据液压多路阀失效的模式的不同，分析可靠度与平均寿命；步骤8：对比分析不同失效模式下可靠度与平均寿命，得出多路阀系统受不同故障形式影响程度的大小；步骤9：求各基本事件的概率重要度系数I_i；步骤10：设定仿真参数；步骤11：进行模拟仿真并分析仿真结果。

Description

一种液压多路阀可靠性分析方法

技术领域

本发明涉及机械液压可靠性分析领域，特别涉及一种液压多路阀可靠性分析方法。

背景技术

作为工程机械中常见的方向控制阀，液压多路阀具有结构紧凑，压力损失小，寿命长，制造简单等特点。工程机械工作环境恶劣，换向频繁，负载变动剧烈的特点对多路阀的可靠性提出了严苛的要求，因此，对于多路阀进行可靠性分析具有重要意义。

但是，现有对于多路阀可靠性的研究，大多给出多路阀可靠性试验台，而用试验台研究多路阀可靠性成本较高，耗时较长。

发明内容

发明目的：本发明的目的在于克服现有多路阀可靠性分析测试方法成本较高，耗时较长及分析不准确的问题，提供一种更加简洁，结果准确性较高的可靠性分析方法。

本发明要解决的一个技术问题是，在液压多路阀失效模式具有隐蔽性，多样性及相关性的的前提下，提出一种液压多路阀可靠性分析方法。

为实现以上技术目的，本发明提供了一种液压多路阀可靠性分析方法，包括以下步骤：

步骤1：对于液压多路阀系统的典型失效模式建立树状图，树状图包括顶事件、中间事件和基本事件；

步骤2：进行参数查阅，查出各个单元的故障率参数；在进行参数查阅时，对于同一液压多路阀系统中出现i次的故障率参数为λ_i的模块，单次取其故障率参数为λ_i/i(本发明中给出的多路阀失效模式树状图中包含溢流阀弹簧故障和溢流阀阻尼孔堵塞两种故障形式，已知溢流阀故障率参数为λy时，取溢流阀弹簧故障和溢流阀阻尼孔堵塞这两种故障形式的故障率参数分别为

和

)；各个单元具体指的是：溢流阀弹簧，溢流阀阻尼孔，单向阀阀芯，端盖，密封圈，工艺螺堵，滑阀，单向阀螺堵，安全阀泄漏，滑阀阀体间隙，节流槽，阀芯与阀体磨损，安全阀损毁，过载阀损毁，油液污染。故障率参数对应的是各个单元，而某些单元的故障并不一定会导致整体多路阀系统的失效。因此，对于各个单元用故障率参数，整体用失效率参数。

步骤3：根据最小割集算法，算出最小割集；

步骤4：根据树状图计算液压多路阀系统可靠度；

步骤5：对于液压多路阀系统，计算其可靠度与平均寿命；

步骤6：根据液压多路阀系统四种失效模式，分析各模块的可靠度与平均寿命；各模块指的是四种失效模式所分别对应的单元集合。

步骤7：对比分析不同失效模式下各模块的可靠度与平均寿命，得出液压多路阀系统受不同故障形式影响程度的大小；

步骤8：求各基本事件的概率重要度系数I_i；

步骤9：设定仿真参数；

步骤10：进行模拟仿真并分析仿真结果。

步骤1中，树状图的顶事件为T，表示多路阀失效；中间事件分别为：Y1，Y2，Y3，Y4，Y5，Y6，分别表示压力波动及噪声、外泄漏、内泄露及功率损失、阀芯换向异常、换向冲击大、换向缓慢或无法换向；

基本事件为：X1，X2，X3，X4，X5，X6，X7，X8，X9，X10，X11，X12，X13，X14，X15，分别表示溢流阀弹簧故障、溢流阀阻尼孔故障、单向阀关闭不严、端盖泄漏、密封圈失效、工艺螺堵失效、滑阀磨损失效、单向阀螺堵漏油、安全阀泄漏、滑阀阀体间隙过大、节流槽与系统不匹配、阀芯与阀体磨损、油液污染、安全阀损毁、过载阀损毁。

步骤3包括：根据最小割集算法推出液压多路阀系统的失效模式为串联模式，根据如下公式计算最小割集：

T＝Y1+Y2+Y3+Y4

＝X1+X2+...+X9+Y5+Y6

＝X1+X2+...+X13+(X14×X15)

计算得到最小割集为：{X1}，{X2}，{X3}，{X4}，{X5}，{X6}，{X7}，{X8}，{X9}，{X10}，{X11}，{X12}，{X13}，{X14，X15}。

步骤4包括：计算各单元的可靠度，公式为：

其中，R_i(t)表示第i个单元的可靠度函数，λ_i(t)表示第i个单元的失效率函数；

根据上述公式计算出各单元的可靠度分别记为R1，R2，R3.....R14，R15；

计算液压多路阀系统可靠度R：

R＝R1×R2×......×R13×[1-(1-R14)(1-R15)]。

步骤5包括：并联的液压多路阀系统的可靠度R(t)计算公式为：

串联的液压多路阀系统的可靠度R(t)计算公式为：

计算液压多路阀系统可靠度R(t)，设运行时间为t：

液压多路阀系统失效率参数λ为：

进而得到液压多路阀系统平均寿命θ为：

步骤6包括：

根据液压多路阀系统四种失效模式：压力波动及噪声失效模式、外泄漏失效模式、内泄漏及功率损失失效模式、阀芯换向异常失效模式，分为四个系统：系统1：压力波动及噪声失效模式的相关单元；系统2：外泄漏失效模式的相关单元；系统3：内泄漏及功率损失失效模式的相关单元；系统4：阀芯换向异常失效模式的相关单元；

计算系统1的可靠度Ra：Ra＝R1×R2×R3；

对于压力波动及噪声失效模式，计算其相关单元的失效率参数λa与平均寿命θa：

计算系统2的可靠度Rb：Rb＝R4×R5×R6；

对于外泄漏失效模式，计算其相关单元的失效率参数λb与平均寿命θb；

计算系统3的可靠度Rc：Rc＝R7×R8×R9×R10；

对于内泄漏及功率损失失效模式，计算其相关单元的失效率参数λc与平均寿命θc；

计算系统4的可靠度Rd：Rd＝R11×R12×R13×[1-(1-R14)(1-R15)]；

对于阀芯换向异常失效模式，计算其相关单元的失效率参数λd与平均寿命θd。

步骤8包括：顶上事件发生函数g(t)为：

g(t)＝1-R1×R2×……×R13×[1-(1-R14)(1-R15)]，

即：g(t)＝1-(1-F1)×(1-F2)×……×(1-F13)×(1-F14×F15)；

F_i表示第i个单元的不可靠度；

各基本事件的概率重要度系数为：

其中i取值为1～15；

步骤9包括：规定蒙特卡罗仿真时总工作时间为T，再将其分割为m个时间段，即[0，T1]，(0，T2]，......，(Tm-1，Tm]，令仿真运行总次数为N，仿真次数累加变量为j，基本事件个数为n，基本事件累加数为i；Tm数值上就等于工作总时间T，T1,T2,T3......Tm是时间跨度确定后，每一时段的右区间端点。此处可以按需要求统计任意时间区间跨度内液压多路阀系统发生故障次数。

步骤10包括：如果取基本事件所对应的部件或模块的故障规律均符合指数分布，同时取仿真系统运行时间为t，第i个单元的失效概率函数F_i(t)为：F_i(t)＝1-e^-λit；

选取[0，1]区间上均匀分布的随机数ξ_ij，再根据F_i(t)的反函数产生第i个基本部件失效时间T_ij：T_ij＝F_i ^-1(ξ_ij)；

由公式得到第j次仿真运行中，所有基本部件或模块的失效时间为T_1j,T_2j,..T_nj，将他们按照由小到大的顺序排列，即：P_1j≤P_2j≤...P_kj≤...Pnj；P_1j指的是，第j次仿真中，所有单元失效时间T_1j,T_2j,...T_nj的最小值，P_2j指的是其中的次小值......依此类推，Pnj指的是其中的最大值；

每次仿真中，当系统仿真时间达到P_kj时，相对应的基本部件或模块发生故障，如果此时顶事件发生即液压多路阀系统失效，则此次仿真停止，如果顶事件未发生，仿真继续直至顶事件发生则此次仿真停止；

仿真结束后，会生成m×n的矩阵A，用来记录此次仿真中哪个基本部件或模块的失效引起了系统的失效，此时系统运行时间记为P_sj，则T_j＝P_sj；T_j表示第j次仿真中，系统失效时的运行时间，数值上等于P_sj；

仿真N次结束后得到N×n的矩阵B，矩阵B记录每次仿真中基本部件或模块的状态情况；

根据仿真所得数据，分别统计落[0,T_μ]内顶事件的累积发生次数C_μ，μ＝1,2,...,m，则此时液压多路阀系统的失效概率F(T_μ)为：F(T_μ)＝C_μ/N，由此得到液压多路阀系统失效曲线；

步骤10中，根据仿真得到的数据得到平均无故障运行时间MTEP：

根据仿真结果计算系统的结构重要度W(X_i)：

W(X_i)＝z/N

其中，z表示Xi失效次数。

有益效果：本发明所提供的解析法与基于蒙特卡罗方法的液压多路阀可靠性分析方法可以针对不同复杂程度的液压多路阀系统灵活的采取方法。本发明提供的方法不仅包含对整体系统可靠性的分析，也包含了对于不同失效模式所对应模块，和底层单元重要度的分析，从而判断出液压多路阀因何种失效形式失效的可能性更大，评估众多导致液压多路阀系统失效的基本事件中，哪一个单元的故障对液压多路阀系统的影响最大。对以后多路阀状态监测及检修维护具有指导意义。

附图说明

下面结合附图和具体实施方式对本发明做更进一步的具体说明，本发明的上述和/或其他方面的优点将会变得更加清楚。

图1是对于多路阀失效模式进行分析之后，得到的多路阀失效模式树状图。

图2是根据最小割集的结果，绘制出的多路阀失效逻辑框图。

图3是基于蒙特卡罗方法的液压多路阀可靠性分析方法流程图。

具体实施方式

所述液压多路阀可靠性分析方法应当包括以下步骤：

步骤S1，对于液压多路阀的典型失效模式建立树状图。如图1所示。顶事件为T：多路阀失效。中间事件分别为：Y1，Y2，Y3，Y4，Y5，Y6。基本事件为：X1，X2，X3，X4，X5，X6，X7，X8，X9，X10，X11，X12，X13，X14，X15。本发明所述的液压多路阀的失效模式树状图涵盖了典型的多路阀失效模式，从直接导致多路阀失效的四种失效模式向下逐层追因，得到了十五个导致多路阀失效的基本事件。所述失效模式的基本事件描述到某种或某些部件的某种故障形式，方便后面对基本事件所对应的失效率进行查找。考虑到液压多路阀种类繁多以及液压失效的隐蔽性，多样性。所述液压多路阀的失效模式树状图可以根据需要进行添加或删减，再进行后面的可靠性计算。这种图形推导演绎方法，能直观形象地表示目标分析系统失效原因与故障之间的逻辑关系，分析全面，分析效果好。

对照液压多路阀失效模式树状图基本事件并参考机械设计手册或厂家提供数据进行失效率参数查找。考虑到多路阀故障具有相关性，故在同一次可靠性分析中，对于不同基本事件中出现同一种部件或模块失效但失效形式不同的，在查阅到的失效率参数后，取λ_i/i(失效率参数/部件或模块出现次数)代入计算。基本事件X1，X2，X3，X4，X5，X6，X7，X8，X9，X10，X11，X12，X13，X14，X15所对应的失效率参数分别为…λ1，λ2，λ3，λ4，λ5，λ6，λ7，λ8，λ9，λ10，

λ11，λ12，λ1＝3，λ14，λ15。

步骤S2，根据最小割集算法推出系统的失效模式为串联模式，下面计算其最小割集：

T＝Y1+Y2+Y3+Y4

＝X1+X2+...+X9+Y5+Y6

＝X1+X2+...+X13+(X14×X15)

计算得其最小割集为：{X1}，{X2}，{X3}，{X4}，{X5}，{X6}，{X7}，{X8}，{X9}，{X10}，{X11}，{X12}，{X13}，{X14，X15}。

表1是中间事件Y1，Y2，Y3，Y4，Y5，Y6所分别对应的多路阀故障或失效模式。

表2是基本事件X1，X2，......X15，所分别对应的多路阀故障形式。

表1

表2

步骤S3，根据最小割集的结果，绘制出多路阀失效的逻辑框图。如图2所示。

步骤S4，计算各部件或模块的失效度。其公式为：

根据公式计算出各部件或模块的可靠度分别记为R1，R2，R3.....R14，R15；

步骤S5，计算液压多路阀系统可靠度：

R＝R1×R2×......×R13×[1-(1-R14)(1-R15)]。

所述液压多路阀系统可靠度计算公式仅适合上面所述的多路阀失效模式。

对于并联的液压多路阀系统的可靠度计算方法，其公式为：

对于串联的液压多路阀系统的可靠度计算方法，其公式为：

步骤S6，计算所述液压多路阀系统的失效率参数与平均寿命。设运行时间为t,

液压多路阀系统失效率参数为：

进而得到系统平均寿命为：

根据所述液压多路阀失效的四种主要模式，将其失效模式分为四个系统。系统1：压力波动及噪声失效模式的相关单元；系统2：外泄漏失效模式的相关单元；系统3：内泄漏及功率损失失效模式的相关单元；系统4：阀芯换向异常失效模式的相关单元。

步骤S7，计算系统1的可靠度Ra：Ra＝R1×R2×R3，

步骤S8，对于压力波动及噪声失效模式，计算其相关部件或模块的失效率参数λa与平均寿命θa：

步骤S9，计算系统2的可靠度Rb：Rb＝R4×R5×R6；

步骤S10，对于外泄漏失效模式，计算其相关部件或模块的失效率参数λb与平均寿命θb。步骤及公式见步骤S8。

步骤S11，计算系统3的可靠度Rc：Rc＝R7×R8×R9×R10；

步骤S12，对于内泄漏及功率损失失效模式，计算其相关部件或模块的失效率参数λc与平均寿命θc。步骤及公式见步骤S8。

步骤S13，计算系统4的可靠度Rd：Rd＝R11×R12×R13×[1-(1-R14)(1-R15)]；

步骤S14，对于阀芯换向异常失效模式，计算其相关部件或模块的失效率参数λd与平均寿命θd。步骤及公式见步骤S8。

步骤S15，对比分析四种失效模式相关部件或模块的可靠度大小及平均寿命大小，可以判断出液压多路阀因何种失效形式失效的可能性更大，对以后多路阀状态监测及检修具有指导意义。

步骤S16，求各基本事件的概率重要度系数I_i，i＝1，2......15。

顶上事件发生函数为：g(t)＝1-R1×R2×……×R13×[1-(1-R14)(1-R15)]；

即：g(t)＝1-(1-F1)×(1-F2)×……×(1-F13)×(1-F14×F15)；

F_i表示第i个单元的不可靠度；

各基本事件的概率重要度系数：

所述基本事件概率重要度系数可以评估众多导致液压多路阀系统失效的基本事件中，哪一种部件或模块的失效对液压多路阀系统的影响最大。对以后多路阀状态监测及检修具有指导意义。

上述是对于简单液压多路阀故障树模型的解析法定量计算，先对液压多路阀进行总体的可靠性分析，能够对液压多路阀的可靠性有一个定性的分析。然后，根据四种不同的失效模式：压力波动及噪声，外泄漏，内泄漏及功率损失和阀芯换向异常把多路阀失效分为四个系统。对比分析四种失效模式相关部件或模块的可靠度大小及平均寿命大小，判断出液压多路阀因何种失效形式失效的可能性更大，对以后多路阀状态监测及检修具有指导意义。最后，分析导致多路阀失效的基本事件概率重要度系数，可以评估众多导致液压多路阀系统失效的基本事件中，哪一种部件或模块的失效对液压多路阀系统的影响最大。

随着液压多路阀的发展，多路阀种类繁多，结构复杂，所包含的各种阀结构及其他部件越来越多，用解析法定量计算是比较困难的。例如当部件故障规律服从任意分布时，就难以采用常规解析法精确求解。

针对上述问题，本发明提供了基于蒙特卡罗方法的液压多路阀可靠性分析方法。

所述基于蒙特卡罗方法的液压多路阀可靠性分析方法应当包含以下步骤：构造模型，定义随机变量，通过模型获得子样，仿真计算。

所述构造模型，即针对所要求解问题的特点，建立一个实用而简便的概率统计模型，使要求得的解恰好等于所建概率统计模型的数学期望、分布的分位点或事件发生的概率。

所述定义随机变量，即跟据问题的要求定义一个随机变量，使它的数字特征值恰好等于所求问题的解。

所述通过模拟获得子样，即根据概率统计模型和概率分布类型找出对应的抽样方法，在计算机上进行数字模拟，随机产生一些样本观察值，它近似具有简单随机子样的性质。

所述仿真计算，即对上述子样进行统计分析，得到有关的数字特征值、分布的分位点或某事件发生的频率，以此作为所求问题的解。

所述基于蒙特卡罗方法的液压多路阀可靠性分析方法流程图如图3所示。

所述基于蒙特卡罗方法的液压多路阀可靠性分析方法应具体步骤为：

步骤1：规定蒙特卡罗仿真总工作时间为T，再将其分割为m个时间段，即[0，T1]，(0，T2]，......，(Tm-1，Tm]，令仿真运行总次数为N，N的取值要待多次试验后，仿真结果趋于稳定后再确定。仿真次数累加变量为j。基本事件个数为n，基本事件累加数为i。

步骤2：如果取基本事件所对应的部件或模块的故障规律均符合指数分布，同时取仿真系统运行时间为t，第i个单元的失效概率函数F_i(t)为：F_i(t)＝1-e^-λit。

当部件故障规律服从其他分布式，应采用其他公式。此处单元故障规律符合正态分布，有些其他单元的故障率不符合正态分布。具体实现需要结合具体单元，查阅后可知其故障规律分布方式。

选取[0，1]区间上均匀分布的随机数ξ_ij，再根据F_i(t)的反函数产生第i个基本部件失效时间。公式为：T_ij＝F_i ^-1(ξ_ij)；

由公式得到第j次仿真运行中，所有基本部件或模块的失效时间为T_1j,T_2j,..T_nj。

将他们按照由小到大的顺序排列，即：P_1j≤P_2j≤...P_kj≤...Pnj；

每次仿真中，当系统仿真时间达到P_kj时，相对应的那个基本部件或模块发生故障；如果此时顶事件发生即液压多路阀系统失效，则此次仿真停止，如果顶事件未发生，仿真继续直至顶事件发生则此次仿真停止。仿真结束后，会生成m×n的矩阵A。用来记录此次仿真中哪个基本部件或模块的失效引起了系统的失效。此时系统运行时间记为P_sj。则T_j＝P_sj。

仿真N次结束后得到N×n的矩阵B，矩阵B记录每次仿真中基本部件或模块的状态情况。

根据仿真实验所得数据，分别统计落[0,T_μ](μ＝1,2,...,m)内顶事件的累积发生次数C_μ，则此时液压多路阀系统的失效概率为：F(T_μ)＝C_μ/N，由此得到液压多路阀系统失效曲线。用以评估液压多路阀系统的可靠度。

根据仿真得到的数据得到平均无故障运行时间MTEP：

根据仿真结果计算系统的结构重要度：

W(X_i)＝z/N

其中z表示Xi失效次数。

上述为基于蒙特卡罗方法的液压多路阀可靠性分析方法，用到MATLAB作为仿真软件。通过上述方法可以近似得到液压多路阀系统的失效概率函数，可以得到系统的平均无故障运行时间和系统的结构重要度。通过以上三个参数，可以精确的衡量液压多路阀系统的可靠性。也可以根据需要，修改程序以获得需要的参数。面对复杂多路阀系统结构复杂时，此方法更方便明了。

本发明提供了液压多路阀可靠性分析方法。解析法与基于蒙特卡罗方法的液压多路阀可靠性分析方法可以针对不同复杂程度的液压多路阀灵活的采取方法。本发明提供的方法不仅包含对整体系统可靠性的分析，也包含了对于不同失效模式所对应的部件或模块，和基本事件底层重要度的分析，从而判断出液压多路阀因何种失效形式失效的可能性更大，评估众多导致液压多路阀系统失效的基本事件中，哪一种部件或模块的失效对液压多路阀系统的影响最大。对以后多路阀状态监测及检修维护具有指导意义。

本发明提供的液压多路阀可靠性分析方法需要基于以下假设：假定所有底事件之间都是独立的；假设顶事件和所有底事件与中间事件都只有两种状态，即发生和不发生；假定故障分布为指数分布。

本发明提供了一种液压多路阀可靠性分析方法，具体实现该技术方案的方法和途径很多，以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。

Claims (10)

1.一种液压多路阀可靠性分析方法，其特征在于，包括以下步骤：

步骤1：对于液压多路阀系统的典型失效模式建立树状图，树状图包括顶事件、中间事件和基本事件；

步骤2：进行参数查阅，查出各个单元的故障率参数；在进行参数查阅时，对于同一液压多路阀系统中出现i次的故障率参数为λ_i的模块，单次取其故障率参数为λ_i/i；

步骤3：根据最小割集算法，算出最小割集；

步骤4：根据树状图计算液压多路阀系统可靠度；

步骤5：对于液压多路阀系统，计算其可靠度与平均寿命；

步骤6：根据液压多路阀系统四种失效模式，分析各模块的可靠度与平均寿命；各模块指的是四种失效模式所分别对应的单元集合；

步骤7：对比分析不同失效模式下各模块的可靠度与平均寿命，得出液压多路阀系统受不同故障形式影响程度的大小；

步骤8：求各基本事件的概率重要度系数I_i；

步骤9：设定仿真参数；

步骤10：进行模拟仿真并分析仿真结果。

2.根据权利要求1所述的方法，其特征在于，步骤1中，树状图的顶事件为T，表示多路阀失效；中间事件分别为：Y1，Y2，Y3，Y4，Y5，Y6，分别表示压力波动及噪声、外泄漏、内泄露及功率损失、阀芯换向异常、换向冲击大、换向缓慢或无法换向；

基本事件为：X1，X2，X3，X4，X5，X6，X7，X8，X9，X10，X11，X12，X13，X14，X15，分别表示溢流阀弹簧故障、溢流阀阻尼孔故障、单向阀关闭不严、端盖泄漏、密封圈失效、工艺螺堵失效、滑阀磨损失效、单向阀螺堵漏油、安全阀泄漏、滑阀阀体间隙过大、节流槽与系统不匹配、阀芯与阀体磨损、油液污染、安全阀损毁、过载阀损毁。

3.根据权利要求2所述的方法，其特征在于，步骤3包括：根据最小割集算法推出液压多路阀系统的失效模式为串联模式，根据如下公式计算最小割集：

T＝Y1+Y2+Y3+Y4

＝X1+X2+...+X9+Y5+Y6

＝X1+X2+...+X13+(X14×X15)

计算得到最小割集为：{X1}，{X2}，{X3}，{X4}，{X5}，{X6}，{X7}，{X8}，{X9}，{X10}，{X11}，{X12}，{X13}，{X14，X15}。

4.根据权利要求3所述的方法，其特征在于，步骤4包括：计算各单元的可靠度，公式为：

其中，R_i(t)表示第i个单元的可靠度函数，λ_i(t)表示第i个单元的失效率函数；

根据上述公式计算出各单元的可靠度分别记为R1，R2，R3.....R14，R15；

计算液压多路阀系统可靠度R：

R＝R1×R2×......×R13×[1-(1-R14)(1-R15)]。

5.根据权利要求4所述的方法，其特征在于，步骤5包括：并联的液压多路阀系统的可靠度R(t)计算公式为：

串联的液压多路阀系统的可靠度R(t)计算公式为：

计算液压多路阀系统可靠度R(t)，设运行时间为t：

液压多路阀系统失效率参数λ为：

进而得到液压多路阀系统平均寿命θ为：

6.根据权利要求5所述的方法，其特征在于，步骤6包括：

根据液压多路阀系统四种失效模式：压力波动及噪声失效模式、外泄漏失效模式、内泄漏及功率损失失效模式、阀芯换向异常失效模式，分为四个系统：系统1：压力波动及噪声失效模式的相关单元；系统2：外泄漏失效模式的相关单元；系统3：内泄漏及功率损失失效模式的相关单元；系统4：阀芯换向异常失效模式的相关单元；

计算系统1的可靠度Ra：Ra＝R1×R2×R3；

对于压力波动及噪声失效模式，计算其相关单元的失效率参数λa与平均寿命θa：

计算系统2的可靠度Rb：Rb＝R4×R5×R6；

对于外泄漏失效模式，计算其相关单元的失效率参数λb与平均寿命θb；

计算系统3的可靠度Rc：Rc＝R7×R8×R9×R10；

对于内泄漏及功率损失失效模式，计算其相关单元的失效率参数λc与平均寿命θc；

计算系统4的可靠度Rd：Rd＝R11×R12×R13×[1-(1-R14)(1-R15)]；

对于阀芯换向异常失效模式，计算其相关单元的失效率参数λd与平均寿命θd。

7.根据权利要求6所述的方法，其特征在于，步骤8包括：顶上事件发生函数g(t)为：g(t)＝1-R1×R2×……×R13×[1-(1-R14)(1-R15)]，

即：g(t)＝1-(1-F1)×(1-F2)×……×(1-F13)×(1-F14×F15)；

F_i表示第i个单元的不可靠度；

各基本事件的概率重要度系数为：

其中i取值为1～15。

8.根据权利要求7所述的方法，其特征在于，步骤9包括：规定蒙特卡罗仿真时总工作时间为T，再将其分割为m个时间段，即[0，T1]，(0，T2]，......，(Tm-1，Tm]，令仿真运行总次数为N，仿真次数累加变量为j，基本事件个数为n，基本事件累加数为i。

9.根据权利要求8所述的方法，其特征在于，步骤10包括：如果取基本事件所对应的单元故障规律均符合指数分布，同时取仿真系统运行时间为t，第i个单元的失效概率函数F_i(t)为：F_i(t)＝1-e^-λit；

选取[0，1]区间上均匀分布的随机数ξ_ij，再根据F_i(t)的反函数产生第i个基本部件失效时间T_ij：T_ij＝F_i ^-1(ξ_ij)；

由公式得到第j次仿真运行中，所有单元的失效时间为T_1j,T_2j,..T_nj，将他们按照由小到大的顺序排列，即：P_1j≤P_2j≤...P_kj≤...Pnj，P_1j指的是，第j次仿真中，所有单元失效时间T_1j,T_2j,...T_nj的最小值，Pnj指的是第j次仿真中，所有单元失效时间T_1j,T_2j,...T_nj的最大值；

每次仿真中，当系统仿真时间达到P_kj时，相对应的基本部件或模块发生故障，如果此时顶事件发生即液压多路阀系统失效，则此次仿真停止，如果顶事件未发生，仿真继续直至顶事件发生则此次仿真停止；

仿真结束后，会生成m×n的矩阵A，用来记录此次仿真中哪个基本部件或模块的失效引起了系统的失效，此时系统运行时间记为P_sj，则T_j＝P_sj；T_j表示第j次仿真中，系统失效时的运行时间，数值上等于P_sj；

仿真N次结束后得到N×n的矩阵B，矩阵B记录每次仿真中基本部件或模块的状态情况；

根据仿真所得数据，分别统计落[0,T_μ]内顶事件的累积发生次数C_μ，μ＝1,2,…,m，则此时液压多路阀系统的失效概率函数F(T_μ)为：F(T_μ)＝C_μ/N，由此得到液压多路阀系统失效曲线。

10.根据权利要求9所述的方法，其特征在于，步骤10中，根据仿真得到的数据得到平均无故障运行时间MTEP：

根据仿真结果计算系统的结构重要度W(X_i)：

W(X_i)＝z/N

其中，z表示Xi失效次数。

Images