CN107702905A - 一种基于威布尔分布的橡胶圈可靠寿命预测方法 - Google Patents
一种基于威布尔分布的橡胶圈可靠寿命预测方法 Download PDFInfo
- Publication number
- CN107702905A CN107702905A CN201710897911.2A CN201710897911A CN107702905A CN 107702905 A CN107702905 A CN 107702905A CN 201710897911 A CN201710897911 A CN 201710897911A CN 107702905 A CN107702905 A CN 107702905A
- Authority
- CN
- China
- Prior art keywords
- mrow
- rubber ring
- msub
- sample
- ijk
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/005—Sealing rings
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种基于威布尔分布的橡胶圈可靠寿命预测方法,其包括以下步骤:S1、通过摸底试验确定橡胶圈样本的压缩永久变形率εik的平均增长速率λi,进而确定橡胶圈样本退化试验方案的试验时间t;S2、根据橡胶圈样本退化试验,得到橡胶圈样本的轴向初始厚度H0、各检测时间点橡胶圈样本的轴向厚度Hijk和各个试验过程中的压缩永久变形率εijk;S3、根据各个试验过程中的压缩永久变形率εijk,利用橡胶老化公式和基于威布尔分布的加速模型预测橡胶圈寿命并分析橡胶圈的可靠性。本发明相较传统的橡胶寿命预测方法得到单一寿命预测值,利用寿命分布函数得到了可靠寿命,使评估结果具有更高的可信性。
Description
技术领域
本发明涉及橡胶圈可靠性研究领域,具体涉及一种基于威布尔分布的橡胶圈可靠寿命预测方法。
背景技术
橡胶圈由于其优异的密封性能和较低的成本在密封行业得到了广泛应用。但是橡胶材料容易老化,从而影响其密封性能。著名的“挑战者”号航天飞机爆炸事故就是由于O型密封圈失效引起的。因此橡胶圈的可靠性评估和寿命预测对工程应用具有十分重要的意义。传统的橡胶圈寿命预测方法只是将各温度水平下测得样本老化性能参数进行了平均处理,没有考虑产品的分布特征,导致预测的精度不高,同时由于没有引入统计分布,这种方法也没有将橡胶密封圈的寿命与可靠性结合起来,预测的结果仅仅是一个寿命值,没有体现产品的可靠性特征。
发明内容
针对现有技术中的上述不足,本发明提供的一种基于威布尔分布的橡胶圈可靠寿命预测方法不仅预测了寿命,还分析了橡胶圈的可靠性。
为了达到上述发明目的,本发明采用的技术方案为:
提供一种基于威布尔分布的橡胶圈可靠寿命预测方法,其包括以下步骤:
S1、通过摸底试验确定橡胶圈样本的压缩永久变形率εik的平均增长速率λi,进而确定橡胶圈样本退化试验方案的试验时间t;
S2、根据橡胶圈样本退化试验的试验时间t进行橡胶圈样本退化试验,得到橡胶圈样本的轴向初始厚度H0、各检测时间点橡胶圈样本的轴向厚度Hijk和各个试验过程中的压缩永久变形率εijk;
S3、根据各个试验过程中的压缩永久变形率εijk,利用橡胶老化公式和基于威布尔分布的加速模型预测橡胶圈的可靠寿命。
进一步地,通过摸底实验确定橡胶圈样本的压缩永久变形率εik的平均增长速率λi,进而确定橡胶圈样本退化试验的试验时间t的方法包括:
S1-1、在室温下将橡胶圈样本装入夹具进行变形压缩至少24小时,取出橡胶圈样本并在静置至少一小时后,测量橡胶圈样本至少四个不同地方的轴向厚度,取其平均值作为摸底试验橡胶圈样本的轴向初始厚度h0;
S1-2、将该橡胶圈样本再次装入夹具中静置至少半小时后放入老化试验箱中进行温度应力水平试验,在检测时间点取出橡胶圈样本并在室温下静置至少一小时,测量橡胶圈样本至少四个不同地方的轴向厚度,取其平均值作为摸底试验橡胶圈样本的检测时间点橡胶圈样本的轴向厚度h1;
S1-3、改变摸底实验中的温度应力水平并重复步骤S1-2,根据公式得到多个橡胶圈样本在试验过程中的变形率εik,进而得到各个温度应力水平压缩永久变形率的平均增长速率λi;其中Hx表示限制器的高度;
S1-4、根据公式
得到试验时间t的范围,其中ε0为压缩永久变形率失效阈值。
进一步地,橡胶圈样本退化试验的方法为:
S2-1、设置试验参数:设置温度应力水平总数u、各应力水平下橡胶圈试验样本量v、试验时间t和测量次数w;
S2-2、测量v个橡胶圈样本轴向初始厚度H0:将橡胶圈样本以实际工作压缩量分别安装在夹具内,室温下至少压缩橡胶圈样本24小时后取出,并在室温下静置橡胶圈样本至少1小时,测量每个橡胶圈样品至少4个不同地方的轴向厚度,取其平均值作为每组橡胶圈样品轴向初始厚度H0;
S2-3、测量各检测时间点橡胶圈样品的轴向厚度Hijk:将每组橡胶圈样品以实际工作压缩量装入夹具并装入温度老化试验箱中;在检测时间点tijk取出橡胶圈样品并放置在室温下恢复弹性变形至少1小时,测量同一橡胶圈样品至少4个不同地方的轴向厚度,取其平均值作为退化试验后的轴向厚度Hijk;
S2-4、计算各个试验过程中的压缩永久变形率εijk:根据公式得到第j个橡胶圈试验样本在第i个温度应力水平、第k次试验时的压缩永久变形率εijk;
其中Hx表示限制器的高度;j=1,2,…,v;i=1,2,…,u;k=1,2,…,w。
进一步地,根据各个试验过程中的压缩永久变形率εijk、利用橡胶老化公式和基于威布尔分布的加速模型预测橡胶圈的可靠寿命方法为:
S3-1、根据第j个橡胶圈试验样本在第i个温度应力水平、第k次检测时的压缩永久变形率εijk和第k次试验检测时间tk,采用公式和橡胶老化经验公式通过逐次逼近法使参数I(α)取最小值得到老化常数α,并利用最小二乘回归分析得到第i个温度应力水平下第j个橡胶圈样本的试验常数Bij和老化反应速率Kij,进而得到各温度应力水平下各橡胶圈样本的老化公式:
ln(1-εijk)=ln Bij-Kijtk α
S3-2、将公式ln(1-εijk)=ln Bij-Kijtijk α中的参数εijk替换为给定压缩永久变形率失效阈值ε0,得到各温度应力水平下各样本的伪失效寿命tij;
S3-3、将各温度应力水平下各样本的伪失效寿命tij数据按从小到大排序,根据威布尔分布的最优线性无偏估计方法:
得到各温度应力水平下形状参数:
得到尺度参数,即特征寿命:
S3-4、根据阿伦尼斯模型:
将阿伦尼斯模型线性化,得到线性加速模型:
将给定的工作温度应力T0带入线性加速模型得到该工作温度下的预测特征寿命η0:
S3-5、根据公式:
得到各温度应力水平下形状参数的平均值m0;
根据公式:
和给定可靠度R预测橡胶圈样本的寿命t(R);
根据公式:
得到橡胶圈样本的可靠性R(t);
其中D(v,v,j)和C(v,v,j)为方差系数,η表示特征寿命,A为常数,Ea是激活能,r是玻尔兹曼常数,T为热力学温度;和均为中间参数;a和b的值根据特征寿命和温度应力水平Ti对线性加速模型进行最小二乘回归分析得到。
本发明的有益效果为:本发明将橡胶圈的压缩永久变形率作为性能退化指标,用温度作为加速应力进行恒定应力加速退化试验,利用试验退化数据,结合橡胶材料经验老化方程推算出各应力水平下各样本的伪失效寿命,利用威布尔分布的最佳线性无偏估计得到分布参数估计值,结合加速方程外推正常工作温度下的分布参数,利用分布函数进行可靠寿命估计,利用退化数据外推产品失效寿命,节省了大量时间且具有良好的可行性;相较传统橡胶寿命预测方法中对测量老化性能参数的平均处理,考虑的产品个体之间的差异性,引入了伪失效寿命分布的概念,充分利用了每一个的样本的试验信息;相较传统的橡胶寿命预测方法得到单一寿命预测值,利用寿命分布函数得到了可靠寿命,使评估结果具有更高的可信性。
附图说明
图1为本发明的流程图;
图2为本发明实施例摸底实验压缩永久变形变化趋势图;
图3为本发明实施例80℃下压缩永久变形变化趋势图;
图4为本发明实施例100℃下压缩永久变化趋势图;
图5为本发明实施例115℃下压缩永久变化趋势图;
图6为本发明实施例在工作温度下可靠度函数曲线图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
如图1所示,该基于威布尔分布的橡胶圈可靠寿命预测方法包括以下步骤:
S1、通过摸底试验确定橡胶圈样本的压缩永久变形率εik的平均增长速率λi,进而确定橡胶圈样本退化试验方案的试验时间t;
S2、根据橡胶圈样本退化试验的试验时间t进行橡胶圈样本退化试验,得到橡胶圈样本的轴向初始厚度H0、各检测时间点橡胶圈样本的轴向厚度Hijk和各个试验过程中的压缩永久变形率εijk;
S3、根据各个试验过程中的压缩永久变形率εijk,利用橡胶老化公式和基于威布尔分布的加速模型预测橡胶圈的可靠寿命。
通过摸底实验确定橡胶圈样本的压缩永久变形率εik的平均增长速率λi,进而确定橡胶圈样本退化试验的试验时间t的方法包括:
S1-1、在室温下将橡胶圈样本装入夹具进行变形压缩至少24小时,取出橡胶圈样本并在静置至少1小时后,测量橡胶圈样本至少四个不同地方的轴向厚度,取其平均值作为摸底试验橡胶圈样本的轴向初始厚度h0;
S1-2、将该橡胶圈样本再次装入夹具中静置至少半小时后放入老化试验箱中进行温度应力水平试验,在检测时间点取出橡胶圈样本并在室温下静置至少1小时,测量橡胶圈样本至少四个不同地方的轴向厚度,取其平均值作为摸底试验橡胶圈样本的检测时间点橡胶圈样本的轴向厚度h1;
S1-3、改变摸底实验中的温度应力水平并重复步骤S1-2,根据公式得到多个橡胶圈样本在试验过程中的变形率εik,进而得到各个温度应力水平压缩永久变形率的平均增长速率λi;其中Hx表示限制器的高度;
S1-4、根据公式
得到试验时间t的范围,其中ε0为压缩永久变形率失效阈值。
橡胶圈样本退化试验的方法为:
S2-1、设置试验参数:设置温度应力水平总数u、各应力水平下橡胶圈试验样本量v、试验时间t和测量次数w;
S2-2、测量v个橡胶圈样本轴向初始厚度H0:将橡胶圈样本以实际工作压缩量分别安装在夹具内,室温下至少压缩橡胶圈样本24小时后取出,并在室温下静置橡胶圈样本至少1小时,测量每个橡胶圈样品至少4个不同地方的轴向厚度,取其平均值作为每组橡胶圈样品轴向初始厚度H0;
S2-3、测量各检测时间点橡胶圈样品的轴向厚度Hijk:将每组橡胶圈样品以实际工作压缩量装入夹具并装入温度老化试验箱中;在检测时间点tijk取出橡胶圈样品并放置在室温下恢复弹性变形至少1小时,测量同一橡胶圈样品至少4个不同地方的轴向厚度,取其平均值作为退化试验后的轴向厚度Hijk;
S2-4、计算各个试验过程中的压缩永久变形率εijk:根据公式得到第j个橡胶圈试验样本在第i个温度应力水平、第k次试验时的压缩永久变形率εijk;
其中Hx表示限制器的高度;j=1,2,…,v;i=1,2,…,u;k=1,2,…,w。
根据各个试验过程中的压缩永久变形率εijk、利用橡胶老化公式和基于威布尔分布的加速模型预测橡胶圈的可靠寿命方法为:
S3-1、根据第j个橡胶圈试验样本在第i个温度应力水平、第k次检测时的压缩永久变形率εijk和第k次试验检测时间tk,采用公式和橡胶老化经验公式通过逐次逼近法使参数I(α)取最小值得到老化常数α,并利用最小二乘回归分析得到第i个温度应力水平下第j个橡胶圈样本的试验常数Bij和老化反应速率Kij,进而得到各温度应力水平下各橡胶圈样本的老化公式:
ln(1-εijk)=ln Bij-Kijtk α
S3-2、将公式ln(1-εijk)=ln Bij-Kijtijk α中的参数εijk替换为给定压缩永久变形率失效阈值ε0,得到各温度应力水平下各样本的伪失效寿命tij;
S3-3、将各温度应力水平下各样本的伪失效寿命tij数据按从小到大排序,根据威布尔分布的最优线性无偏估计方法:
得到各温度应力水平下形状参数:
得到尺度参数,即特征寿命:
S3-4、根据阿伦尼斯模型:
将阿伦尼斯模型线性化,得到线性加速模型:
将给定的工作温度应力T0带入线性加速模型得到该工作温度下的预测特征寿命η0:
S3-5、根据公式:
得到各温度应力水平下形状参数的平均值m0;
根据公式:
和给定可靠度R预测橡胶圈样本的寿命t(R);
根据公式:
得到橡胶圈样本的可靠性R(t);
其中D(v,v,j)和C(v,v,j)为方差系数,η表示特征寿命,A为常数,Ea是激活能,r是玻尔兹曼常数,T为热力学温度;和均为中间参数;a和b的值根据特征寿命和温度应力水平Ti对线性加速模型进行最小二乘回归分析得到。
在本发明的一个实施例中,某机械密封O型氟橡胶圈,其工作温度为60℃,压缩永久变形失效阈值为0.5,按上述具体实施方式对其进行试验并进行寿命预测和可靠性分析:
步骤1:开展O型氟橡胶圈摸底试验。
选取4个样本,温度水平为80℃、100℃、115℃进行加速试验,每24小时测量一次厚度。其性能退化轨迹见图2。
选取4个加速应力水平分别为80℃、100℃、115℃,各温度应力水平下的样本量为8个,由步骤1的性能退化轨迹,设定试验时间为16天,每48小时测量一次压缩永久变形率。
步骤2:开展O型氟橡胶圈加速退化试验。
试验过程中测得的数据见表1。
表1试验过程中样品压缩永久变形率(%)
80℃下压缩永久变形变化趋势图见图3;100℃下压缩永久变形变化趋势图见图4;115℃下压缩永久变形变化趋势图见图5。
步骤3:试验数据处理。
由试验数据,应用逐次逼近法计算得到α=0.5,通过老化方程ln(1-εijk)=ln Bij-Kijtk α结合表1数据进行回归分析,可以得到各应力水平下个样本的老化方程系数,见表2。
表2依据试验样品数据的老化方程系数回归分析
依据老化方程,代入O型氟橡胶圈压缩永久变形率失效阈值ε0=0.5,计算出各应力水平下个样本的伪失效寿命,见表3。
表3各试验样本的伪失效寿命
利用式(4)~(6)对伪失效寿命进行威布尔参数估计,计算中间过程参数见表4~表6。
表4 80℃下试验样本伪失效寿命最佳线性无偏估计相关系数表
形状参数尺度参数
表5 100℃下试验样本伪失效寿命最佳线性无偏估计相关系数表
序号 | tij(h) | lntij | C(v,v,j) | C(v,v,j)lntij | D(v,v,j) | D(v,v,j)lntij |
1 | 1000 | 6.9078 | -0.1019 | -0.7039 | 0.0365 | 0.2521 |
2 | 1021 | 6.9285 | -0.1081 | -0.7490 | 0.0561 | 0.3887 |
3 | 1049 | 6.9556 | -0.1027 | -0.7143 | 0.0759 | 0.5279 |
4 | 1153 | 7.0501 | -0.0872 | -0.6148 | 0.0971 | 0.6846 |
5 | 1240 | 7.1229 | -0.0589 | -0.4195 | 0.1212 | 0.8633 |
6 | 1309 | 7.1770 | -0.0111 | -0.0797 | 0.1502 | 1.0780 |
7 | 1315 | 7.1816 | 0.0758 | 0.5444 | 0.1894 | 1.3602 |
8 | 1412 | 7.2528 | 0.3942 | 2.8590 | 0.2735 | 1.9836 |
Σ | — | — | — | 0.1222 | — | 7.1384 |
形状参数尺度参数
表6 115℃下试验样本伪失效寿命最佳线性无偏估计相关系数表
形状参数尺度参数
将各个应力水平下的威布尔参数估计值汇于表7中。
表7各应力水平下试验样本伪失效寿命的威布尔参数估计值
温度(℃) | 形状参数m | 尺度参数η |
80 | 9.61 | 2553 |
100 | 8.18 | 1259 |
115 | 9.09 | 745 |
利用式(8)的加速方程,结合表7进行回归分析可得到加速方程的参数的估计值,见表8。
表8加速方程的参数的估计值
a | b | R-square |
-5.772 | 4811 | 0.9995 |
得到加速方程将O型圈工作温度60℃,即333.15K,代入可得到η的估计值:η0=5819。
形状参数m的估计值为:
步骤4:O型氟橡胶密封圈寿命预测与可靠性分析。
O型氟橡胶圈可靠度为0.9的可靠寿命为
O型氟橡胶圈可靠度为0.8的可靠寿命为
可靠度函数为:
O型氟橡胶圈的可靠度随时间变化曲线见图6。
综上所述,本发明将橡胶圈的压缩永久变形率作为性能退化指标,用温度作为加速应力进行恒定应力加速退化试验,利用试验退化数据,结合橡胶材料经验老化方程推算出各应力水平下各样本的伪失效寿命,利用威布尔分布的最佳线性无偏估计得到分布参数估计值,结合加速方程外推正常工作温度下的分布参数,利用分布函数进行可靠寿命估计,利用退化数据外推产品失效寿命,节省了大量时间且具有良好的可行性;相较传统橡胶寿命预测方法中对测量老化性能参数的平均处理,考虑的产品个体之间的差异性,引入了伪失效寿命分布的概念,充分利用了每一个的样本的试验信息;相较传统的橡胶寿命预测方法得到单一寿命预测值,利用寿命分布函数得到了可靠寿命,使评估结果具有更高的可信性。
Claims (4)
1.一种基于威布尔分布的橡胶圈可靠寿命预测方法,其特征在于:包括以下步骤:
S1、通过摸底试验确定橡胶圈样本的压缩永久变形率εik的平均增长速率λi,进而确定橡胶圈样本退化试验的试验时间t;
S2、根据橡胶圈样本退化试验的试验时间t进行橡胶圈样本退化试验,得到橡胶圈样本的轴向初始厚度H0、各检测时间点橡胶圈样本的轴向厚度Hijk和各个试验过程中的压缩永久变形率εijk;
S3、根据各个试验过程中的压缩永久变形率εijk,利用橡胶老化公式和基于威布尔分布的加速模型预测橡胶圈的可靠寿命。
2.根据权利要求1所述的基于威布尔分布的橡胶圈可靠寿命预测方法,其特征在于,通过摸底实验确定橡胶圈样本的压缩永久变形率εik的平均增长速率λi,进而确定橡胶圈样本退化试验的试验时间t的方法包括:
S1-1、在室温下将橡胶圈样本装入夹具进行变形压缩至少24小时,取出橡胶圈样本并在静置至少一小时后,测量橡胶圈样本至少四个不同地方的轴向厚度,取其平均值作为摸底试验橡胶圈样本的轴向初始厚度h0;
S1-2、将该橡胶圈样本再次装入夹具中静置至少半小时后放入老化试验箱中进行温度应力水平试验,在检测时间点取出橡胶圈样本并在室温下静置至少一小时,测量橡胶圈样本至少四个不同地方的轴向厚度,取其平均值作为摸底试验橡胶圈样本的检测时间点橡胶圈样本的轴向厚度h1;
S1-3、改变摸底实验中的温度应力水平并重复步骤S1-2,根据公式得到多个橡胶圈样本在试验过程中的变形率εik,进而得到各个温度应力水平压缩永久变形率的平均增长速率λi;其中Hx表示限制器的高度;
S1-4、根据公式
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo>{</mo>
<mfrac>
<mrow>
<mn>0.1</mn>
<msub>
<mi>&epsiv;</mi>
<mn>0</mn>
</msub>
</mrow>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>}</mo>
<mo>&le;</mo>
<mi>t</mi>
<mo>&le;</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo>{</mo>
<mfrac>
<msub>
<mi>&epsiv;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>}</mo>
</mrow>
得到试验时间t的范围,其中ε0为压缩永久变形率失效阈值。
3.根据权利要求2所述的基于威布尔分布的橡胶圈可靠寿命预测方法,其特征在于,橡胶圈样本退化试验的方法为:
S2-1、设置试验参数:设置温度应力水平总数u、各应力水平下橡胶圈试验样本量v、试验时间t和测量次数w;
S2-2、测量v个橡胶圈样本轴向初始厚度H0:将橡胶圈样本以实际工作压缩量分别安装在夹具内,室温下至少压缩橡胶圈样本24小时后取出,并在室温下静置橡胶圈样本至少1小时,测量每个橡胶圈样品至少4个不同地方的轴向厚度,取其平均值作为每组橡胶圈样品轴向初始厚度H0;
S2-3、测量各检测时间点橡胶圈样品的轴向厚度Hijk:将每组橡胶圈样品以实际工作压缩量装入夹具并装入温度老化试验箱中;在检测时间点tijk取出橡胶圈样品并放置在室温下恢复弹性变形至少1小时,测量同一橡胶圈样品至少4个不同地方的轴向厚度,取其平均值作为退化试验后的轴向厚度Hijk;
S2-4、计算各个试验过程中的压缩永久变形率εijk:根据公式得到第j个橡胶圈试验样本在第i个温度应力水平、第k次试验时的压缩永久变形率εijk;
其中Hx表示限制器的高度;j=1,2,…,v;i=1,2,…,u;k=1,2,…,w。
4.根据权利要求3所述的基于威布尔分布的橡胶圈可靠寿命预测方法,其特征在于,所述根据各个试验过程中的压缩永久变形率εijk、利用橡胶老化公式和基于威布尔分布的加速模型预测橡胶圈的可靠寿命方法为:
S3-1、根据第j个橡胶圈试验样本在第i个温度应力水平、第k次检测时的压缩永久变形率εijk和第k次试验检测时间tk,采用公式和橡胶老化经验公式通过逐次逼近法使参数I(α)取最小值得到老化常数α,并利用最小二乘回归分析得到第i个温度应力水平下第j个橡胶圈样本的试验常数Bij和老化反应速率Kij,进而得到各温度应力水平下各橡胶圈样本的老化公式:
ln(1-εijk)=ln Bij-Kijtk α
S3-2、将公式ln(1-εijk)=ln Bij-Kijtijk α中的参数εijk替换为给定压缩永久变形率失效阈值ε0,得到各温度应力水平下各样本的伪失效寿命tij;
S3-3、将各温度应力水平下各样本的伪失效寿命tij数据按从小到大排序,根据威布尔分布的最优线性无偏估计方法:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>&mu;</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>v</mi>
</munderover>
<mo>&lsqb;</mo>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>,</mo>
<mi>v</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mi>ln</mi>
<mi> </mi>
<msub>
<mi>t</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>v</mi>
</munderover>
<mo>&lsqb;</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>,</mo>
<mi>v</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mi>ln</mi>
<mi> </mi>
<msub>
<mi>t</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>u</mi>
</mrow>
得到各温度应力水平下形状参数:
<mrow>
<msub>
<mover>
<mi>m</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
</mfrac>
</mrow>
得到尺度参数,即特征寿命:
<mrow>
<msub>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>&mu;</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
S3-4、根据阿伦尼斯模型:
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<msup>
<mi>Ae</mi>
<mrow>
<mo>-</mo>
<mfrac>
<msub>
<mi>E</mi>
<mi>a</mi>
</msub>
<mrow>
<mi>r</mi>
<mi>T</mi>
</mrow>
</mfrac>
</mrow>
</msup>
</mrow>
将阿伦尼斯模型线性化,得到线性加速模型:
<mrow>
<mi>l</mi>
<mi>n</mi>
<mi>&eta;</mi>
<mo>=</mo>
<mi>a</mi>
<mo>+</mo>
<mi>b</mi>
<mo>&CenterDot;</mo>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
</mrow>
将给定的工作温度应力T0带入线性加速模型得到该工作温度下的预测特征寿命η0:
<mrow>
<msub>
<mi>ln&eta;</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mi>a</mi>
<mo>+</mo>
<mi>b</mi>
<mo>&CenterDot;</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
</mfrac>
</mrow>
S3-5、根据公式:
<mrow>
<msub>
<mi>m</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mi>i</mi>
<mi>u</mi>
</munderover>
<msub>
<mover>
<mi>m</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
</mrow>
<mi>u</mi>
</mfrac>
</mrow>
得到各温度应力水平下形状参数的平均值m0;
根据公式:
<mrow>
<mi>t</mi>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&eta;</mi>
<mn>0</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>ln</mi>
<mi> </mi>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mn>0</mn>
</msub>
</mfrac>
</msup>
</mrow>
和给定可靠度R预测橡胶圈样本的寿命t(R);
根据公式:
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>t</mi>
<msub>
<mi>&eta;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>m</mi>
<mn>0</mn>
</msub>
</msup>
<mo>)</mo>
</mrow>
</mrow>
得到橡胶圈样本的可靠性R(t);
其中D(v,v,j)和C(v,v,j)为方差系数,η表示特征寿命,A为常数,Ea是激活能,r是玻尔兹曼常数,T为热力学温度;和均为中间参数;a和b的值根据特征寿命和温度应力水平Ti对线性加速模型进行最小二乘回归分析得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710897911.2A CN107702905B (zh) | 2017-09-28 | 2017-09-28 | 一种基于威布尔分布的橡胶圈可靠寿命预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710897911.2A CN107702905B (zh) | 2017-09-28 | 2017-09-28 | 一种基于威布尔分布的橡胶圈可靠寿命预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107702905A true CN107702905A (zh) | 2018-02-16 |
CN107702905B CN107702905B (zh) | 2020-01-10 |
Family
ID=61175899
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710897911.2A Expired - Fee Related CN107702905B (zh) | 2017-09-28 | 2017-09-28 | 一种基于威布尔分布的橡胶圈可靠寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107702905B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109145502A (zh) * | 2018-09-17 | 2019-01-04 | 中国人民解放军海军工程大学 | 一种威布尔型单元寿命分布参数估计方法 |
CN110207909A (zh) * | 2018-02-28 | 2019-09-06 | 全球能源互联网研究院有限公司 | 一种管路系统密封结构的老化寿命预测方法及压缩工装 |
CN111380772A (zh) * | 2019-04-30 | 2020-07-07 | 中国人民解放军63961部队 | 长期贮存产品内部的橡胶类密封件贮存寿命的测定方法 |
CN111751093A (zh) * | 2020-06-16 | 2020-10-09 | 江苏高聚识别技术有限公司 | 用于身份识别研发的系统寿命检验系统 |
CN112560237A (zh) * | 2020-12-03 | 2021-03-26 | 中国市政工程华北设计研究总院有限公司 | 一种燃气胶管的寿命快速测试方法及其测试系统 |
CN112686397A (zh) * | 2020-12-18 | 2021-04-20 | 广东机电职业技术学院 | 空调器故障率计算方法、装置、计算机设备及存储介质 |
CN114065548A (zh) * | 2021-11-30 | 2022-02-18 | 重庆长安新能源汽车科技有限公司 | 一种预测电池包箱体盖密封寿命的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102494992A (zh) * | 2011-12-13 | 2012-06-13 | 北京航空航天大学 | 基于步进应力的丁腈橡胶o型密封圈的加速退化试验方法 |
CN104236889A (zh) * | 2014-09-11 | 2014-12-24 | 航天科工防御技术研究试验中心 | 橡胶密封圈的贮存寿命优化设计方法 |
CN105352815A (zh) * | 2015-10-22 | 2016-02-24 | 广州机械科学研究院有限公司 | 一种静密封氟橡胶o形圈剩余使用寿命的预测方法 |
CN105844079A (zh) * | 2016-03-15 | 2016-08-10 | 北京航空航天大学 | 压缩载荷下氟橡胶密封产品加速老化模型与寿命评估方法 |
CN106124191A (zh) * | 2016-08-31 | 2016-11-16 | 中国兵器工业第五九研究所 | 一种橡胶密封圈的剩余寿命评估方法 |
-
2017
- 2017-09-28 CN CN201710897911.2A patent/CN107702905B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102494992A (zh) * | 2011-12-13 | 2012-06-13 | 北京航空航天大学 | 基于步进应力的丁腈橡胶o型密封圈的加速退化试验方法 |
CN104236889A (zh) * | 2014-09-11 | 2014-12-24 | 航天科工防御技术研究试验中心 | 橡胶密封圈的贮存寿命优化设计方法 |
CN105352815A (zh) * | 2015-10-22 | 2016-02-24 | 广州机械科学研究院有限公司 | 一种静密封氟橡胶o形圈剩余使用寿命的预测方法 |
CN105844079A (zh) * | 2016-03-15 | 2016-08-10 | 北京航空航天大学 | 压缩载荷下氟橡胶密封产品加速老化模型与寿命评估方法 |
CN106124191A (zh) * | 2016-08-31 | 2016-11-16 | 中国兵器工业第五九研究所 | 一种橡胶密封圈的剩余寿命评估方法 |
Non-Patent Citations (2)
Title |
---|
刘炜等: "基于加速老化试验的某型战车橡胶密封圈寿命预测", 《质量与可靠性》 * |
肖坤等: "基于恒定应力加速退化试验的某引信用O型橡胶密封圈可靠性评估", 《机械工程学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110207909A (zh) * | 2018-02-28 | 2019-09-06 | 全球能源互联网研究院有限公司 | 一种管路系统密封结构的老化寿命预测方法及压缩工装 |
CN110207909B (zh) * | 2018-02-28 | 2023-02-03 | 全球能源互联网研究院有限公司 | 一种管路系统密封结构的老化寿命预测方法及压缩工装 |
CN109145502A (zh) * | 2018-09-17 | 2019-01-04 | 中国人民解放军海军工程大学 | 一种威布尔型单元寿命分布参数估计方法 |
CN111380772A (zh) * | 2019-04-30 | 2020-07-07 | 中国人民解放军63961部队 | 长期贮存产品内部的橡胶类密封件贮存寿命的测定方法 |
CN111751093A (zh) * | 2020-06-16 | 2020-10-09 | 江苏高聚识别技术有限公司 | 用于身份识别研发的系统寿命检验系统 |
CN112560237A (zh) * | 2020-12-03 | 2021-03-26 | 中国市政工程华北设计研究总院有限公司 | 一种燃气胶管的寿命快速测试方法及其测试系统 |
CN112686397A (zh) * | 2020-12-18 | 2021-04-20 | 广东机电职业技术学院 | 空调器故障率计算方法、装置、计算机设备及存储介质 |
CN114065548A (zh) * | 2021-11-30 | 2022-02-18 | 重庆长安新能源汽车科技有限公司 | 一种预测电池包箱体盖密封寿命的方法 |
CN114065548B (zh) * | 2021-11-30 | 2024-05-07 | 深蓝汽车科技有限公司 | 一种预测电池包箱体盖密封寿命的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107702905B (zh) | 2020-01-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107702905A (zh) | 一种基于威布尔分布的橡胶圈可靠寿命预测方法 | |
CN104020401B (zh) | 基于云模型理论的变压器绝缘热老化状态的评估方法 | |
CN102494992B (zh) | 基于步进应力的丁腈橡胶o型密封圈的加速退化试验方法 | |
CN111458225B (zh) | 一种质子交换膜燃料电池密封材料的寿命预测方法 | |
CN105242155A (zh) | 一种基于熵权法和灰色关联分析的变压器故障诊断方法 | |
CN107436963A (zh) | 一种基于Copula函数多元退化的O型橡胶密封圈寿命预测方法 | |
CN111024528B (zh) | 一种沥青混合料粘弹性参数测量方法 | |
CN103279657B (zh) | 一种基于工程经验的产品加速退化试验方案设计方法 | |
CN103246821A (zh) | 一种基于仿真的多应力小样本加速寿命试验方案设计优化方法 | |
CN105352815B (zh) | 一种静密封氟橡胶o形圈剩余使用寿命的预测方法 | |
CN105844079A (zh) | 压缩载荷下氟橡胶密封产品加速老化模型与寿命评估方法 | |
CN104778354A (zh) | 定期抽检型产品贮存寿命评估方法 | |
CN103678938B (zh) | 一种面向空间形状和误差范围的退化模型一致性检验方法 | |
CN106021685A (zh) | 一种考虑测量误差的退化可靠性分析方法 | |
CN105158147B (zh) | 一种用于测试密封圈材料老化的装置及方法 | |
CN103970129B (zh) | 控制阀粘滞检测方法 | |
Johlitz et al. | Chemical ageing of elastomers: experiments and modelling | |
CN107561452A (zh) | 一种蓄电池组健康度的评估方法 | |
CN109359375B (zh) | 一种数据产品贮存寿命预测方法 | |
CN116773239A (zh) | 一种智能燃气表控制器可靠性寿命预测方法 | |
CN105956977A (zh) | 采用全观测指标方案的灰色关联分析健康分级评价方法 | |
CN113792446B (zh) | 一种基于晶体塑性理论的蠕变疲劳剩余寿命评定方法 | |
CN110991001B (zh) | 一种基于单调回归理论的卷簧寿命评估方法 | |
Baussaron et al. | Degradation test plan for Wiener degradation processes | |
CN113361025B (zh) | 一种基于机器学习的蠕变疲劳概率损伤评定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200110 Termination date: 20200928 |