[发明专利]先导式安全阀贮存寿命加速试验方法

摘要

本发明一种先导式安全阀贮存寿命加速试验方法，包括以下步骤：步骤1、设计先导式安全阀的加速试验方案；步骤2、将所述样品放入温度试验箱内进行试验；步骤3、在每个加速应力水平下，按照所设定的测试时间对样品进行性能测试；测试项目包括样品中导阀弹簧负荷损失率σ_D(t)和主阀弹簧负荷损失率σ_Z(t)，以及O形圈压缩永久变形率G_M(t)；步骤4、对导阀弹簧负荷损失率σ_D(t)值、主阀弹簧负荷损失率σ_Z(t)值和O形圈压缩永久变形率G_M(t)值进行处理，从而对其贮存寿命进行预测。本发明提供了一套基于加速试验实现先导式安全阀贮存寿命预测的方法。该方法能够反映先导式安全阀库房常规贮存的性能退化进程，以较大的加速系数进行加速试验，缩短了试验时间，为先导式安全阀库房贮存寿命研究提供了一种可行的技术方法。

主权项

1.一种先导式安全阀贮存寿命加速试验方法，其特征是包括以下步骤： 步骤1、设计先导式安全阀的加速试验方案： 1-1以温度作为试验的加速应力，以不同的温度作为试验的加速应力水平； 1-2加速应力水平数等于或大于3，其中，最低加速应力水平应高于库房常规贮存温度，最高加速应力水平的设置必须保证先导式安全阀的退化机理与库房常规贮存时保持一致； 1-3每个加速应力水平下的先导式安全阀样品数等于或大于4； 1-4在每个加速应力水平下对每个样品进行性能测试，测试次数至少为4次，采用先密后疏的原则设定测试时间； 步骤2、将所述样品放入温度试验箱内进行试验； 步骤3、在每个加速应力水平下，按照所设定的测试时间对样品进行性能测试；测试项目包括样品中导阀弹簧负荷损失率σ_D(t)和主阀弹簧负荷损失率σ_Z(t)，以及O形圈压缩永久变形率G_M(t)；对于每次测试，步骤如下： 3-1到达所设定的测试时间时停止试验，并从试验箱中取出样品，等待样品在常温下充分冷却； 3-2从样品中拆卸出导阀弹簧、主阀弹簧和O形圈； 3-3导阀弹簧负荷损失率σ_D(t)和主阀弹簧负荷损失率σ_Z(t)的测试步骤： ①将导阀弹簧和主阀弹簧分别在室内环境下自然恢复1小时； ②将导阀弹簧和主阀弹簧分别压缩至它们在阀体内的安装高度，采用弹力测试仪器（如弹簧试验机）对此时导阀弹簧和主阀弹簧的弹力分别进行测试并记录，每个导阀弹簧和主阀弹簧分别测量5次，分别取它们的平均值作为相应的弹力测试值； ③分别根据式（1）和式（2）计算导阀弹簧负荷损失率σ_D(t)和主阀弹簧负荷损失率σ_Z(t) 其中，t为试验时间，F_D(0)、F_D(t)分别为每个加速应力水平下试验初始时刻和试验至t时 刻时，将导阀弹簧压缩至阀体内安装高度时的弹力； 其中，F_Z(0)、F_Z(t)分别为每个加速应力水平下试验初始时刻和试验至t时刻时，将主阀弹簧压缩至阀体内安装高度时的弹力； 3-4O形圈压缩永久变形率G_M(t)的测试步骤： ①将O形圈在室内环境下自然恢复1小时； ②采用非接触式测量仪器（如数字显微镜）测量O形圈的截面直径，并记录；对于每个O形圈，将其绕中心每旋转一个角度测一次，至少测量6次，取平均值作为O形圈截面直径测试值； ③根据式（3）计算O形圈的压缩永久变形率G_M(t) 其中，H_M(0)、H_M(t)分别为每个加速应力水平下试验初始时刻和试验至t时刻时的O形圈截面直径测试值；H_M1为O形圈在样品中被压缩后的截面直径测试值； 3-5将导阀弹簧、主阀弹簧和O形圈装回样品中，并放入温度试验箱继续试验，直至步骤1中所设计的试验方案全部执行完毕后结束试验； 步骤4、将测试过程中由步骤3的式（1）~式（3）计算得到的所有导阀弹簧负荷损失率σ_D(t)值、主阀弹簧负荷损失率σ_Z(t)值和O形圈压缩永久变形率G_M(t)值分别记为y_D、y_Z和y_M，y_D={y_Dij(t_i，k)|i＝1,...,I;j=1,...,J;k=1,...,K}、y_Z={y_Zij(t_i,k)|i＝1,...,I;j=1,...,J;k=1,...,K}、y_M={y_Mij(t_i,k)|i＝1,...,I;j=1,...,J;k=1,...,K}，其中i表示I个加速应力水平从低到高的编号，j表示每个加速应力水平下J个样品的编号，k表示每个加速应力水平下K次测试的编号，y_Dij(t_i,k)、y_Zij(t_i,k)和y_Mij(t_i,k)分别表示在第i个加速应力水平下对第j个样品进行第k次测试获得的σ_D(t)值、σ_Z(t)值和G_M(t)值； 将y_D、y_Z和y_M作为测试数据进行处理，步骤如下： 4-1对y_D进行分析，得出导阀弹簧在库房常规贮存下的贮存可靠度①采用式（4）描述导阀弹簧负荷损失率与时间的关系（即退化模型）： 其中为模型参数，α_Dij和β_Dij分别为待估计的参数；②令取不同的尝试值，并代入由式（5）表示的误差平方和使取得最小值时所对应的尝试值即为的估计值，记为i＝1,...,I;j=1,...,J      （5）③将代入式（4）中的得出结合厂家给出的导阀弹簧弹性退化失效阈值σ_Df，通过式（6）求取的反函数获得伪失效寿命时间t_Dij：所有样品的导阀弹簧伪失效寿命时间记为t_D，t_D={t_Dij|i＝1,...,I;j=1,...,J}； ④导阀弹簧贮存可靠度R_Di(t)采用式（7）所示的Weibull分布进行描述： 其中，m_D和η_Di分别为导阀弹簧失效分布的形状参数和第i个加速应力水平下的尺度参数； 加速模型采用式（8）所示的Arrhenius模型进行描述： ln η_Di=γ_D0+γ_D1/S_i      （8） 其中γ_D0和γ_D1是模型参数，S_i为第i个加速应力水平； 将式（8）中的η_Di代入式（7），得到式（9）： 将式（9）中的未知参数γ_D0、γ_D1、m_D记为ψ_D，ψ_D=(γ_D0,γ_D1,m_D)； ⑤通过式（10）计算t_D的似然函数L_D(ψ_D|t_D)： 其中R_Di(t_Dij)和h_Di(t_Dij)分别为将t_Dij代入式（7）和式（11）所取得的值， ⑥对L_D(ψ_D|t_D)求对数得到对数似然函数ln L_D(ψ_D|t_D)，令ψ_D取不同的尝试值并代入ln L_D(ψ_D|t_D)，使ln L_D(ψ_D|t_D)取得最大值时所对应的尝试值即为ψ_D的估计值，记为⑦将代入式（13）得到导阀弹簧在库房常规贮存下的贮存可靠度其中S₀为库房常规贮存时的温度； 4-2对y_Z进行分析，得出主阀弹簧在库房常规贮存下的贮存可靠度①采用式（14）描述导阀弹簧负荷损失率与时间的关系（即退化模型）： 其中为模型参数，α_Zij和β_Zij分别为待估计的参数；②令取不同的尝试值，并代入由式（15）表示的误差平方和使取得最小值时所对应的尝试值即为的估计值，记为i＝1,...,I;j=1,...,J      （15）③将代入式（14）中的得出结合厂家给出的主阀弹簧弹性退化失效阈值σ_Zf，通过式（16）求取的反函数获得伪失效寿命时间t_Zij：所有样品的主阀弹簧伪失效寿命时间记为t_Z，t_Z={t_Zij|i＝1,...,I;j=1,...,J}； ④主阀弹簧贮存可靠度R_Zi(t)采用式（17）所示的Weibull分布进行描述： 其中，m_Z和η_Zi分别为主阀弹簧失效分布的形状参数和第i个加速应力水平下的尺度参数； 加速模型采用式（18）所示的Arrhenius模型进行描述： ln η_Zi=γ_Z0+γ_Z1/S_i      （18） 其中γ_Z0和γ_Z1是模型参数，S_i为第i个加速应力水平； 将式（18）中的η_Zi代入式（17），得到式（19）： 式（19）中的未知参数γ_Z0、γ_Z1、m_Z记为ψ_Z，ψ_Z=(γ_Z0,γ_Z1,m_Z)； ⑤通过式（20）计算t_Z的似然函数L_Z(ψ_Z|t_Z)： 其中R_Zi(t_Zij)和h_Zi(t_Zij)分别为将t_Zij代入式（17）和式（21）所取得的值， ⑥对L_Z(ψ_Z|t_Z)求对数得到对数似然函数ln L_Z(ψ_Z|t_Z)，令ψ_Z取不同的尝试值并代入ln L_Z(ψ_Z|t_Z)，使ln L_Z(ψ_Z|t_Z)取得最大值时所对应的尝试值即为ψ_Z的估计值，记为⑦将代入式（23）得到主阀弹簧在库房常规贮存下的贮存可靠度4-3对y_M进行分析，得出O形圈在库房常规贮存下的贮存可靠度①采用式（24）描述O形圈压缩永久变形率与时间的关系（即退化模型）： 其中为模型参数，K_Mij为速度参数，与试验温度有关且样本间存在差异；a_M为取值在0至1区间的常数，a_M∈(0,1)；②计算的参数估计值先求常数a_M的最优值，令a_M在(0,1)区间等间隔取N_a个离散值（N_a在计算量允许情况下应尽量大，如99），a_M的第l个离散值记为a_Ml，a_Ml=l/(N_a+1),l=1,...,N_a；逐一选取a_M的 第l个离散值a_M1，l=1,...,N_a，然后令K_Mij取不同的尝试值并代入式（25）： i＝1,...,I;j=1,...,J      （25）使得式（25）取得最小值时所对应的尝试值即为K_Mij的估计值，记为将和a_Ml代入式（25）得到计算所有样品的和，记为I(a_Ml)，则在I(a_Ml),l=1,...,N_a中取得最小值时所对应的a_Ml即为最优值a_M^*：a_Ml=l/(N_a-1),l=1,...,N_a      （26）然后令K_Mij取不同的尝试值，并与a_M^*一起代入式（25），使得式（25）取得最小值时所对应的尝试值即为K_Mij的最终估计值K_Mij^*，由此得出的参数估计值③将代入式（24）中的得出结合厂家给出的O形圈密封失效阈值G_Mf，通过式（27）求取的反函数获得伪失效寿命时间t_Mij：所有样品的O形圈伪失效寿命时间记为t_M，t_M={t_Mij|i＝1,...,I;j=1,...,J}； ④O形圈贮存可靠度采用式（28）所示的Weibull分布进行描述： 其中，m_M和η_Mi分别为O形圈失效分布的形状参数和第i个加速应力水平下的尺度参数； 加速模型采用式（29）所示的Arrhenius模型进行描述： ln η_Mi=γ_M0+γ_M1/S_i      （29） 其中γ_M0和γ_M1是模型参数，S_i为第i个加速应力水平； 将式（29）中的η_Zi代入式（28），得到式（30）： 式（30）中的未知参数γ_M0、γ_M1、m_M记为ψ_M，ψ_M=(γ_M0,γ_M1,m_M)； ⑤通过式（31）计算t_M的似然函数L_M(ψ_M|t_M)： 其中R_Mi(t_Mij)和h_Mi(t_Mij)分别为将t_Mij代入式（28）和式（32）所取得的值， ⑥对L_M(ψ_M|t_M)求对数得到对数似然函数ln L_M(ψ_M|t_M)，令ψ_M取不同的尝试值并代入ln L_M(ψ_M|t_M)，使ln L_M(ψ_M|t_M)取得最大值时所对应的尝试值即为ψ_M的估计值，记为⑦将代入式（34）得到O形圈在库房常规贮存下的贮存可靠度4-4将4-1~4-3中得出的和库房常规贮存时的温度S₀代入式（35），得出先导式安全阀在库房常规贮存下的贮存可靠度4-5结合厂家给出的贮存可靠度额定值RT，通过式（36）求取的反函数得到先导式安全阀的贮存寿命从而对其贮存寿命进行预测：