CN107702905B - 一种基于威布尔分布的橡胶圈可靠寿命预测方法 - Google Patents

一种基于威布尔分布的橡胶圈可靠寿命预测方法 Download PDF

Info

Publication number
CN107702905B
CN107702905B CN201710897911.2A CN201710897911A CN107702905B CN 107702905 B CN107702905 B CN 107702905B CN 201710897911 A CN201710897911 A CN 201710897911A CN 107702905 B CN107702905 B CN 107702905B
Authority
CN
China
Prior art keywords
rubber ring
test
sample
ijk
ring sample
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.)
Expired - Fee Related
Application number
CN201710897911.2A
Other languages
English (en)
Other versions
CN107702905A (zh
Inventor
黄洪钟
任彬
李彦锋
刘俊
黄鹏
彭卫文
许焕卫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201710897911.2A priority Critical patent/CN107702905B/zh
Publication of CN107702905A publication Critical patent/CN107702905A/zh
Application granted granted Critical
Publication of CN107702905B publication Critical patent/CN107702905B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/005Sealing 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,根据公式
Figure GDA0002280888500000021
得到多个橡胶圈样本在试验过程中的变形率εik,进而得到各个温度应力水平压缩永久变形率的平均增长速率λi;其中Hx表示限制器的高度;
S1-4、根据公式
Figure GDA0002280888500000022
得到试验时间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:根据公式
Figure GDA0002280888500000031
得到第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,采用公式
Figure GDA0002280888500000032
和橡胶老化经验公式
Figure GDA0002280888500000033
通过逐次逼近法使参数I(α)取最小值得到老化常数α,并利用最小二乘回归分析得到第i个温度应力水平下第j个橡胶圈样本的试验常数Bij和老化反应速率Kij,进而得到各温度应力水平下各橡胶圈样本的老化公式:
ln(1-εijk)=lnBij-Kijtk α
S3-2、将公式ln(1-εijk)=lnBij-Kijtijk α中的参数εijk替换为给定压缩永久变形率失效阈值ε0,得到各温度应力水平下各样本的伪失效寿命tij
S3-3、将各温度应力水平下各样本的伪失效寿命tij数据按从小到大排序,根据威布尔分布的最优线性无偏估计方法:
Figure GDA0002280888500000041
得到各温度应力水平下形状参数:
Figure GDA0002280888500000042
得到尺度参数,即特征寿命:
Figure GDA0002280888500000043
S3-4、根据阿伦尼斯模型:
Figure GDA0002280888500000044
将阿伦尼斯模型线性化,得到线性加速模型:
Figure GDA0002280888500000045
将给定的工作温度应力T0带入线性加速模型得到该工作温度下的预测特征寿命η0
Figure GDA0002280888500000046
S3-5、根据公式:
得到各温度应力水平下形状参数的平均值m0
根据公式:
Figure GDA0002280888500000051
和给定可靠度R预测橡胶圈样本的寿命t(R);
根据公式:
Figure GDA0002280888500000052
得到橡胶圈样本的可靠性R(t);
其中D(v,v,j)和C(v,v,j)为方差系数,η表示特征寿命,A为常数,Ea是激活能,r是玻尔兹曼常数,T为热力学温度;
Figure GDA0002280888500000053
Figure GDA0002280888500000054
均为中间参数;a和b的值根据特征寿命
Figure GDA0002280888500000055
和温度应力水平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,根据公式
Figure GDA0002280888500000071
得到多个橡胶圈样本在试验过程中的变形率εik,进而得到各个温度应力水平压缩永久变形率的平均增长速率λi;其中Hx表示限制器的高度;
S1-4、根据公式
Figure GDA0002280888500000072
得到试验时间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:根据公式
Figure GDA0002280888500000081
得到第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,采用公式
Figure GDA0002280888500000082
和橡胶老化经验公式
Figure GDA0002280888500000083
通过逐次逼近法使参数I(α)取最小值得到老化常数α,并利用最小二乘回归分析得到第i个温度应力水平下第j个橡胶圈样本的试验常数Bij和老化反应速率Kij,进而得到各温度应力水平下各橡胶圈样本的老化公式:
ln(1-εijk)=lnBij-Kijtk α
S3-2、将公式ln(1-εijk)=lnBij-Kijtijk α中的参数εijk替换为给定压缩永久变形率失效阈值ε0,得到各温度应力水平下各样本的伪失效寿命tij
S3-3、将各温度应力水平下各样本的伪失效寿命tij数据按从小到大排序,根据威布尔分布的最优线性无偏估计方法:
Figure GDA0002280888500000084
得到各温度应力水平下形状参数:
得到尺度参数,即特征寿命:
Figure GDA0002280888500000091
S3-4、根据阿伦尼斯模型:
Figure GDA0002280888500000092
将阿伦尼斯模型线性化,得到线性加速模型:
Figure GDA0002280888500000093
将给定的工作温度应力T0带入线性加速模型得到该工作温度下的预测特征寿命η0
Figure GDA0002280888500000094
S3-5、根据公式:
Figure GDA0002280888500000095
得到各温度应力水平下形状参数的平均值m0
根据公式:
Figure GDA0002280888500000096
和给定可靠度R预测橡胶圈样本的寿命t(R);
根据公式:
得到橡胶圈样本的可靠性R(t);
其中D(v,v,j)和C(v,v,j)为方差系数,η表示特征寿命,A为常数,Ea是激活能,r是玻尔兹曼常数,T为热力学温度;
Figure GDA0002280888500000098
Figure GDA0002280888500000099
均为中间参数;a和b的值根据特征寿命
Figure GDA0002280888500000101
和温度应力水平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试验过程中样品压缩永久变形率(%)
Figure GDA0002280888500000102
Figure GDA0002280888500000111
80℃下压缩永久变形变化趋势图见图3;100℃下压缩永久变形变化趋势图见图4;115℃下压缩永久变形变化趋势图见图5。
步骤3:试验数据处理。
由试验数据,应用逐次逼近法计算得到α=0.5,通过老化方程ln(1-εijk)=lnBij-Kijtk α结合表1数据进行回归分析,可以得到各应力水平下个样本的老化方程系数,见表2。
表2依据试验样品数据的老化方程系数回归分析
Figure GDA0002280888500000112
依据老化方程,代入O型氟橡胶圈压缩永久变形率失效阈值ε0=0.5,计算出各应力水平下各样本的伪失效寿命,见表3。
表3各试验样本的伪失效寿命
Figure GDA0002280888500000122
对伪失效寿命进行威布尔参数估计,计算中间过程参数见表4~表6。
表4 80℃下试验样本伪失效寿命最佳线性无偏估计相关系数表
Figure GDA0002280888500000132
Figure GDA0002280888500000141
形状参数
Figure GDA0002280888500000142
尺度参数
Figure GDA0002280888500000143
表5 100℃下试验样本伪失效寿命最佳线性无偏估计相关系数表
序号 t<sub>ij</sub>(h) lnt<sub>ij</sub> C(v,v,j) C(v,v,j)lnt<sub>ij</sub> D(v,v,j) D(v,v,j)lnt<sub>ij</sub>
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
形状参数
Figure GDA0002280888500000144
尺度参数
Figure GDA0002280888500000145
表6 115℃下试验样本伪失效寿命最佳线性无偏估计相关系数表
Figure GDA0002280888500000146
Figure GDA0002280888500000151
形状参数
Figure GDA0002280888500000152
尺度参数
将各个应力水平下的威布尔参数估计值汇于表7中。
表7各应力水平下试验样本伪失效寿命的威布尔参数估计值
温度(℃) 形状参数m 尺度参数η
80 9.61 2553
100 8.18 1259
115 9.09 745
利用线性加速模型的加速方程,结合表7进行回归分析可得到加速方程的参数的估计值,见表8。
表8加速方程的参数的估计值
a b R-square
-5.772 4811 0.9995
得到加速方程
Figure GDA0002280888500000154
将O型圈工作温度60℃,即333.15K,代入可得到η的估计值:η0=5819。
形状参数m的估计值为:
Figure GDA0002280888500000155
步骤4:O型氟橡胶密封圈寿命预测与可靠性分析。
O型氟橡胶圈可靠度为0.9的可靠寿命为
Figure GDA0002280888500000156
O型氟橡胶圈可靠度为0.8的可靠寿命为
Figure GDA0002280888500000157
可靠度函数为:
Figure GDA0002280888500000158
O型氟橡胶圈的可靠度随时间变化曲线见图6。
综上所述,本发明将橡胶圈的压缩永久变形率作为性能退化指标,用温度作为加速应力进行恒定应力加速退化试验,利用试验退化数据,结合橡胶材料经验老化方程推算出各应力水平下各样本的伪失效寿命,利用威布尔分布的最佳线性无偏估计得到分布参数估计值,结合加速方程外推正常工作温度下的分布参数,利用分布函数进行可靠寿命估计,利用退化数据外推产品失效寿命,节省了大量时间且具有良好的可行性;相较传统橡胶寿命预测方法中对测量老化性能参数的平均处理,考虑的产品个体之间的差异性,引入了伪失效寿命分布的概念,充分利用了每一个的样本的试验信息;相较传统的橡胶寿命预测方法得到单一寿命预测值,利用寿命分布函数得到了可靠寿命,使评估结果具有更高的可信性。

Claims (1)

1.一种基于威布尔分布的橡胶圈可靠寿命预测方法,其特征在于:包括以下步骤:
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,根据公式
Figure FDA0002280888490000011
得到多个橡胶圈样本在试验过程中的变形率εik,进而得到各个温度应力水平压缩永久变形率的平均增长速率λi;其中Hx表示限制器的高度;
S1-4、根据公式
Figure FDA0002280888490000021
得到试验时间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:根据公式
Figure FDA0002280888490000022
得到第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,采用公式
Figure FDA0002280888490000031
和橡胶老化经验公式通过逐次逼近法使参数I(α)取最小值得到老化常数α,并利用最小二乘回归分析得到第i个温度应力水平下第j个橡胶圈样本的试验常数Bij和老化反应速率Kij,进而得到各温度应力水平下各橡胶圈样本的老化公式:
ln(1-εijk)=lnBij-Kijtk α
S3-2、将公式ln(1-εijk)=lnBij-Kijtijk α中的参数εijk替换为给定压缩永久变形率失效阈值ε0,得到各温度应力水平下各样本的伪失效寿命tij
S3-3、将各温度应力水平下各样本的伪失效寿命tij数据按从小到大排序,根据威布尔分布的最优线性无偏估计方法:
Figure FDA0002280888490000033
得到各温度应力水平下形状参数:
Figure FDA0002280888490000034
得到尺度参数,即特征寿命:
Figure FDA0002280888490000035
S3-4、根据阿伦尼斯模型:
Figure FDA0002280888490000036
将阿伦尼斯模型线性化,得到线性加速模型:
Figure FDA0002280888490000037
将给定的工作温度应力T0带入线性加速模型得到该工作温度下的预测特征寿命η0
Figure FDA0002280888490000041
S3-5、根据公式:
Figure FDA0002280888490000042
得到各温度应力水平下形状参数的平均值m0
根据公式:
Figure FDA0002280888490000043
和给定可靠度R预测橡胶圈样本的寿命t(R);
根据公式:
得到橡胶圈样本的可靠性R(t);
其中D(v,v,j)和C(v,v,j)为方差系数,η表示特征寿命,A为常数,Ea是激活能,r是玻尔兹曼常数,T为热力学温度;
Figure FDA0002280888490000045
Figure FDA0002280888490000046
均为中间参数;a和b的值根据特征寿命
Figure FDA0002280888490000047
和温度应力水平Ti对线性加速模型进行最小二乘回归分析得到。
CN201710897911.2A 2017-09-28 2017-09-28 一种基于威布尔分布的橡胶圈可靠寿命预测方法 Expired - Fee Related CN107702905B (zh)

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 CN107702905A (zh) 2018-02-16
CN107702905B true 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)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110207909B (zh) * 2018-02-28 2023-02-03 全球能源互联网研究院有限公司 一种管路系统密封结构的老化寿命预测方法及压缩工装
CN109145502B (zh) * 2018-09-17 2023-05-12 中国人民解放军海军工程大学 一种威布尔型单元寿命分布参数估计方法
CN111380772A (zh) * 2019-04-30 2020-07-07 中国人民解放军63961部队 长期贮存产品内部的橡胶类密封件贮存寿命的测定方法
CN111751093A (zh) * 2020-06-16 2020-10-09 江苏高聚识别技术有限公司 用于身份识别研发的系统寿命检验系统
CN112560237B (zh) * 2020-12-03 2023-05-09 中国市政工程华北设计研究总院有限公司 一种燃气胶管的寿命快速测试方法及其测试系统
CN112686397B (zh) * 2020-12-18 2024-06-11 广东机电职业技术学院 空调器故障率计算方法、装置、计算机设备及存储介质
CN114065548B (zh) * 2021-11-30 2024-05-07 深蓝汽车科技有限公司 一种预测电池包箱体盖密封寿命的方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102494992B (zh) * 2011-12-13 2013-07-24 北京航空航天大学 基于步进应力的丁腈橡胶o型密封圈的加速退化试验方法
CN104236889B (zh) * 2014-09-11 2017-02-22 航天科工防御技术研究试验中心 橡胶密封圈的贮存寿命优化设计方法
CN105352815B (zh) * 2015-10-22 2018-04-10 广州机械科学研究院有限公司 一种静密封氟橡胶o形圈剩余使用寿命的预测方法
CN105844079B (zh) * 2016-03-15 2018-08-28 北京航空航天大学 压缩载荷下氟橡胶密封产品加速老化模型与寿命评估方法
CN106124191A (zh) * 2016-08-31 2016-11-16 中国兵器工业第五九研究所 一种橡胶密封圈的剩余寿命评估方法

Also Published As

Publication number Publication date
CN107702905A (zh) 2018-02-16

Similar Documents

Publication Publication Date Title
CN107702905B (zh) 一种基于威布尔分布的橡胶圈可靠寿命预测方法
CN110851980B (zh) 一种设备剩余寿命预测方法及系统
CN103630343B (zh) 隔振器贮存寿命加速试验方法
CN107436963B (zh) 一种基于Copula函数多元退化的O型橡胶密封圈寿命预测方法
CN107885928B (zh) 考虑测量误差的步进应力加速性能退化可靠性分析方法
CN111929514B (zh) 一种电气元器件可靠性评价及寿命预测方法
CN111751199A (zh) 基于eifs分布的疲劳寿命预测方法
CN112331281B (zh) 基于环境大数据和机器学习的高分子材料服役寿命预测方法
CN113567795A (zh) 基于威布尔分布的步退应力加速可靠性试验方法
Qiu et al. Probability analysis for failure assessment of electric energy metering equipment under multiple extreme stresses
CN109030548B (zh) 基于变活化能的聚合物材料热老化寿命评估方法
CN114970157B (zh) 电子产品在电压应力作用下的小样本试验寿命预测方法
CN109342310B (zh) 一种基于老化损伤的橡胶密封制品寿命预测方法
CN108846239A (zh) 一种基于温湿度的弹性环氧树脂加速贮存试验及评估方法
CN109766518B (zh) 考虑样本个体差异的不确定加速退化建模和分析方法
CN112304478A (zh) 一种基于蠕变轮廓法的残余应力测试方法
CN111380772A (zh) 长期贮存产品内部的橡胶类密封件贮存寿命的测定方法
CN111307483B (zh) 一种机械产品长周期疲劳试验数据处理与趋势预判方法
CN113761755A (zh) 考虑认知不确定性的温湿度双应力下的加速寿命分析方法
CN110889077B (zh) 一种加速贮存与自然贮存退化数据一致性检验方法
Hao et al. A novel quality requirement design method for the quality characteristic of rubber products based on the reliability constraint
CN113361025B (zh) 一种基于机器学习的蠕变疲劳概率损伤评定方法
CN113916763A (zh) 一种甲基乙烯基硅橡胶湿热老化寿命预测的方法
CN110972562B (zh) 一种橡胶材料环境适应性实验室快速表征方法
CN113960298A (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