CN110889187A - 基于等退化量时间Pearson系数的加速贮存与自然贮存退化数据一致性检验法 - Google Patents

基于等退化量时间Pearson系数的加速贮存与自然贮存退化数据一致性检验法 Download PDF

Info

Publication number
CN110889187A
CN110889187A CN201811049215.7A CN201811049215A CN110889187A CN 110889187 A CN110889187 A CN 110889187A CN 201811049215 A CN201811049215 A CN 201811049215A CN 110889187 A CN110889187 A CN 110889187A
Authority
CN
China
Prior art keywords
data
degradation
storage
accelerated
test
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
Application number
CN201811049215.7A
Other languages
English (en)
Other versions
CN110889187B (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.)
Hunan Gingko Reliability Technology Research Institute Co Ltd
Original Assignee
Hunan Gingko Reliability Technology Research Institute Co Ltd
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 Hunan Gingko Reliability Technology Research Institute Co Ltd filed Critical Hunan Gingko Reliability Technology Research Institute Co Ltd
Priority to CN201811049215.7A priority Critical patent/CN110889187B/zh
Publication of CN110889187A publication Critical patent/CN110889187A/zh
Application granted granted Critical
Publication of CN110889187B publication Critical patent/CN110889187B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Abstract

一种基于等退化量时间Pearson相关系数的加速贮存与自然贮存退化数据一致性检验方法。方法是对来自加速贮存试验的数据进行回归拟合，得到与自然贮存数据等退化量水平下所需经历的试验时间序列。通过计算各加速应力水平与自然贮存下所对应的试验时间间隔序列的Pearson相关系数，判断两组序列的相关性，从而推断加速贮存与自然贮存下的失效机理是否保持一致。该一致性检验方法对加速试验的实施方案及其试验数据的有效性进行验证，以保证产品寿命预测与验证的有效性。

Description

基于等退化量时间Pearson系数的加速贮存与自然贮存退化 数据一致性检验法
一、技术领域
本发明涉及产品在不同试验应力水平下的退化数据一致性检验方法，具体是基于等退化量时间Pearson系数的加速贮存与自然贮存退化数据一致性检验法，属于可靠性建模技术与寿命预测分析领域。
二、背景技术
对于长期贮存产品，其有效贮存寿命是重要设计和使用指标之一。但长贮产品在正常贮存应力下状态变化极其缓慢，为了尽快获取产品的贮存失效规律并预测贮存寿命，通常采用提高应力的方式加速产品的贮存失效。对于某些长贮产品，即使在加速应力下，也很难观测到失效寿命数据，只能通过监测产品某些关键性能参数的退化失效规律，预测产品在正常应力下的贮存寿命。为了保证这种统计推断的可信性，必须论证产品在正常应力和加速应力下贮存时具有相同的失效机理，即提高应力水平仅仅增大失效速率而不改变失效机理。这是进行加速贮存退化试验设计的重要前提，也是解决长贮产品贮存寿命预测的关键问题和基础问题。目前，国内外文献中关于加速试验失效机理一致性检验方法的研究，主要分为三类：第一类是突发失效产品失效机理的一致性检验，如正态分布方差的一致性；第二类是失效机理明确的退化失效产品，一般采用参数方法，检验两种应力下的退化轨道是否属于同一族而参数各异的随机过程；第三类是失效机理不明确的退化失效产品，一般凭借专家经验进行定性判断，此时结论具有一定的主观性，难以推广应用。一方面由于现代长贮产品的工艺复杂度提高，明确产品退化失效机理的难度相应增大，另一方面由于产品检测水平的提高，收集到的产品检测数据的种类和数量都相对充足。因此，基于数据驱动的非参数检验方法可有效提高试验数据一致性检验的水平。Pearson 系数是在统计学方法下得到的，可以定量地衡量变量之间线性相关关系的相关系数，因此是用来检验模型线性化程度的系数。它需要满足以下条件：（1）两变量均为测量得到的连续变量；（2）变量必须是成对的数据；（3）两变量间为线性关系。因此，对于已知退化模型为非线性模型的情形，对其进行线性转换后再进行线性相关检验。本发明是基于等退化量时间Pearson系数的加速贮存与自然贮存退化数据一致性检验可以对加速试验的实施方案及其试验退化数据的有效性进行初步地验证。
三、发明内容
（一）目的
本发明的目的是进行加速贮存退化数据和自然贮存退化数据一致性检验，它能够检验加速贮存退化数据的有效性：一方面确保产品在试验过程中退化机理一致，验证加速试验的有效性；另一方面，提高产品寿命预测与验证的可信度和精度。
（二）技术方案
本发明一种基于等退化量时间Pearson系数的加速贮存与自然贮存退化数据一致性检验法，该发明的基本思想是：若自然贮存与加速贮存下的失效机理一致，则各性能参数的退化趋势基本保持不变，呈现极强的正相关性，加速可视为于一个时间压缩的过程。
本发明一种基于等退化量时间Pearson系数的加速贮存与自然贮存退化数据一致性检验法，其步骤如下：
步骤1，分别采集自然贮存退化数据和加速贮存退化数据，所述自然贮存退化数据的采集方式通常是专门的自然贮存试验或产品现场检测，加速贮存退化数据的采集于加速贮存试验中同一应力（通常为温度、湿度、盐雾等）下的不同水平的退化数据。
步骤2，对退化数据进行回归分析，得到表征退化量与贮存时间关系的回归方程，其中自然贮存退化数据可得到自然贮存环境下的回归方程，利用不同加速应力水平下的退化数据可以得到不同应力水平下的回归方程。
步骤3，给定若干退化量水平，通过解回归方程，可以得到各给定退化量水平下所需经历的试验时间。
步骤4，计算自然贮存与各应力水平下退化时间数据的Pearson系数。
步骤5，利用Pearson系数对自然贮存和所对应应力水平下的退化失效过程进行判定，判定规则包括不具有相关性、具有一定相关性、暂时无法下结论，需采用其他相关性检验方法进一步判别。其中Pearson系数越接近1表征相关性越强。
其中，在步骤1所述的加速贮存试验应为恒定应力加速退化试验。恒定应力加速退化试验为工程上最常见最便利开展的加速试验类型，若实际进行的加速试验为步进应力或序进应力，则需采用一定的统计方法，将数据等效转化为恒定应力加速试验下的数据。
其中，在步骤1所述的加速贮存试验中，在参试样品量的设置上，工程上比较普遍的做法是各加速贮存应力水平下分别投入一个或几个样品进行性能试验，并获得各样品的性能监测数据。
其中，在步骤1所述的分别采集自然贮存退化数据和加速贮存退化数据，要求：
（1）自然贮存环境下至少获得一个同型产品的性能监测数据，在自然贮存环境下，若同时获得了多个同型产品的性能监测数据，则需先根据各产品监测时刻点，采用插值方法对性能数据进行插值，将各产品的测试时刻对齐，然后得到各测试时刻的样本均值，再将样本均值随贮存时间变化的序列看作是单样品性能变化数据，这样就将自然贮存环境下的多样本数据转化成了单样本退化序列数据。
（2）加速贮存环境下至少获得一个同型产品的性能监测数据，在加速贮存环境下，若同时获得了多个同型产品的性能监测数据，则需先根据各产品监测时刻点，采用插值方法对性能数据进行插值，将各产品的测试时刻对齐，然后得到各测试时刻的样本均值，再将样本均值随贮存时间变化的序列看作是单样品性能变化数据，这样就将加速贮存应力下的多样本数据转化成了单样本退化序列数据。
四、附图说明
图1为本发明的数据一致性检验方法的流程图。
图2为不同线性相关程度的数据点特征图。
五、具体实施方式
本发明一种基于等退化量时间Pearson系数的加速贮存与自然贮存退化数据一致性检验法基于如下假设：
（1）产品所进行的加速贮存试验为恒定应力加速退化试验， 若需要验证一致性的产品不满足此条件，即实际进行的加速试验为步进应力或序进应力，则需采用一定的统计方法，将数据等效转化为恒定应力加速试验下的数据；
（2）各加速贮存应力水平下分别投入一个或几个样品进行性能试验，并获得各样品的性能监测数据；
（3）自然贮存环境下至少获得一个同型产品的性能监测数据。在自然贮存环境下，若同时获得了多个同型产品的性能监测数据，则需先根据各产品监测时刻点，采用插值方法对性能数据进行插值，将各产品的测试时刻对齐；然后得到各测试时刻的样本均值；再将样本均值随贮存时间变化的序列看作是单样品性能变化数据，这样就将自然贮存环境下的多样本数据转化成了单样本退化序列数据。对于加速贮存试验下的多样本数据也需作类似处理。
如图1所示，本发明一种基于等退化量时间Pearson系数的加速贮存与自然贮存退化数据一致性检验法，其步骤如下：
步骤1，分别采集自然贮存退化数据和加速贮存退化数据，所述自然贮存退化数据的采集方式通常是专门的自然贮存试验或产品的现场检测，加速贮存退化数据的采集于加速贮存试验中同一应力（通常为温度、湿度、盐雾等）的不同水平下的退化数据。
步骤2，根据自然贮存下的试验数据，进行回归分析，得到表征自然贮存环境下退化量与贮存时间关系的回归方程
Figure DEST_PATH_IMAGE002
；根据应力水平
Figure DEST_PATH_IMAGE004
下的试验数据，进行回归分析，得到表征应力水平
Figure 981793DEST_PATH_IMAGE004
下退化量与贮存时间关系的回归方程
Figure DEST_PATH_IMAGE006
，
Figure DEST_PATH_IMAGE008
，
Figure DEST_PATH_IMAGE010
为进行恒定应力加速退化试验的应力水平数。
步骤3，给定若干退化量水平
Figure DEST_PATH_IMAGE012
，令
Figure DEST_PATH_IMAGE014
，
Figure DEST_PATH_IMAGE016
，
Figure 148201DEST_PATH_IMAGE008
，
Figure DEST_PATH_IMAGE018
。解回归方程，得到达到给定退化量水平
Figure DEST_PATH_IMAGE020
所需经历的试验时间
Figure DEST_PATH_IMAGE022
，
Figure DEST_PATH_IMAGE024
，
Figure 52572DEST_PATH_IMAGE018
，其中
Figure DEST_PATH_IMAGE026
表示自然贮存。
步骤4，计算自然贮存
Figure DEST_PATH_IMAGE028
与应力水平
Figure 634732DEST_PATH_IMAGE004
下退化时间数据的相关系数
Figure DEST_PATH_IMAGE030
，
Figure DEST_PATH_IMAGE032
其中，
Figure DEST_PATH_IMAGE034
，
Figure DEST_PATH_IMAGE036
，
Figure 505561DEST_PATH_IMAGE008
。
步骤5，制定判别规则。若
Figure DEST_PATH_IMAGE038
，则可判定自然贮存
Figure 514974DEST_PATH_IMAGE028
与应力水平
Figure 719691DEST_PATH_IMAGE004
下的退化失效过程具有一定的相关性；若
Figure DEST_PATH_IMAGE040
，则可初步判定自然贮存
Figure 105542DEST_PATH_IMAGE028
与应力水平
Figure 774420DEST_PATH_IMAGE004
下的退化失效过程不具有相关性，可能存在失效机理不一致的问题，需对加速退化试验条件进行进一步的物理审查，是否由于应力的增加而引入了新的失效机理；若
Figure DEST_PATH_IMAGE042
，则暂时无法下结论，需采用其他相关性检验方法进一步判别。
兹举实施案例如下：
本案例以发射药为例，陈述本发明一种基于等退化量时间Pearson系数的加速贮存与自然贮存退化数据一致性检验法的应用。
本案例的加速贮存试验设计如下。
（1）采用温度单一恒定应力，应力分别为50℃、60℃、70℃、80℃、90℃。
（2）试验样品全部采用新出厂的发射药。
（3）加热设备采用恒温水浴炉。
（4）使用液相色谱仪定时取样测定安定剂含量。
通过以上加速贮存试验，记录不同安定剂剩余含量下的试验时间数据如表1所示。
案例实施流程为上述五个步骤。针对本案例，通过步骤一至步骤四得到加速应力为363K、353K、343K、333K、323K与自然贮存（293K）下的Pearson相关系数如表2所示；,从表1可知，相关性都非常接近1，由步骤五制定的判定规则，可以判定上述各个应力下的失效机理与自然贮存时的失效机理相同。
表1不同温度应力下的试验时间
Figure DEST_PATH_IMAGE044
Figure DEST_PATH_IMAGE046
表2 发射药的Pearson系数表
温度/K 363 353 343 333 323
相关系数 0.999 0. 998 0.998 0.996 0.994

Claims (4)

1.一种基于等退化量时间Pearson系数的加速贮存与自然贮存退化数据一致性检验法，该发明的基本思想是：若自然贮存与加速贮存下的失效机理一致，则各性能参数的退化趋势基本保持不变，呈现极强的正相关性，加速可视为于一个时间压缩的过程；其特征在于，该方法的具体步骤如下：
步骤1，分别采集自然贮存退化数据和加速贮存退化数据，所述自然贮存退化数据的采集方式通常是专门的自然贮存试验或产品现场检测，加速贮存退化数据的采集于加速贮存试验中同一应力（通常为温度、湿度、盐雾等）下的不同水平的退化数据；
步骤2，对退化数据进行回归分析，得到表征退化量与贮存时间关系的回归方程，其中自然贮存退化数据可得到自然贮存环境下的回归方程，利用不同加速应力水平下的退化数据可以得到不同应力水平下的回归方程；
步骤3，给定若干退化量水平，通过解回归方程，可以得到各给定退化量水平下所需经历的试验时间；
步骤4，计算自然贮存与各应力水平下退化时间数据的Pearson系数；
步骤5，利用Pearson系数对自然贮存和所对应应力水平下的退化失效过程进行判定，判定规则包括不具有相关性、具有一定相关性、暂时无法下结论，需采用其他相关性检验方法进一步判别；
其中Pearson系数越接近1表征相关性越强。
2.根据权利要求1所述的一种基于等退化量时间Pearson系数的加速贮存与自然贮存退化数据一致性检验法，其特征在于：在步骤1所述的加速贮存试验应为恒定应力加速退化试验，
恒定应力加速退化试验为工程上最常见最便利开展的加速试验类型，若实际进行的加速试验为步进应力或序进应力，则需采用一定的统计方法，将数据等效转化为恒定应力加速试验下的数据。
3.根据权利要求1所述的一种基于等退化量时间Pearson系数的加速贮存与自然贮存退化数据一致性检验法，其特征在于：在步骤1所述的加速贮存试验中，在参试样品量的设置上，工程上比较普遍的做法是各加速贮存应力水平下分别投入一个或几个样品进行性能试验，并获得各样品的性能监测数据。
4.根据权利要求1所述的一种基于等退化量时间Pearson系数的加速贮存与自然贮存退化数据一致性检验法，其特征在于：在步骤1所述的分别采集自然贮存退化数据和加速贮存退化数据，要求（1）自然贮存环境下至少获得一个同型产品的性能监测数据，在自然贮存环境下，若同时获得了多个同型产品的性能监测数据，则需先根据各产品监测时刻点，采用插值方法对性能数据进行插值，将各产品的测试时刻对齐，然后得到各测试时刻的样本均值，再将样本均值随贮存时间变化的序列看作是单样品性能变化数据，这样就将自然贮存环境下的多样本数据转化成了单样本退化序列数据；（2）加速贮存环境下至少获得一个同型产品的性能监测数据，在加速贮存环境下，若同时获得了多个同型产品的性能监测数据，则需先根据各产品监测时刻点，采用插值方法对性能数据进行插值，将各产品的测试时刻对齐，然后得到各测试时刻的样本均值，再将样本均值随贮存时间变化的序列看作是单样品性能变化数据，这样就将加速贮存应力下的多样本数据转化成了单样本退化序列数据。
CN201811049215.7A 2018-09-10 2018-09-10 基于Pearson系数的贮存退化数据一致性检验法 Active CN110889187B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811049215.7A CN110889187B (zh) 2018-09-10 2018-09-10 基于Pearson系数的贮存退化数据一致性检验法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811049215.7A CN110889187B (zh) 2018-09-10 2018-09-10 基于Pearson系数的贮存退化数据一致性检验法

Publications (2)

Publication Number Publication Date
CN110889187A true CN110889187A (zh) 2020-03-17
CN110889187B CN110889187B (zh) 2020-12-11

Family

ID=69744877

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811049215.7A Active CN110889187B (zh) 2018-09-10 2018-09-10 基于Pearson系数的贮存退化数据一致性检验法

Country Status (1)

Country Link
CN (1) CN110889187B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114091627A (zh) * 2022-01-20 2022-02-25 华谱科仪(北京)科技有限公司 色谱仪试验数据验证方法、装置、存储介质及电子设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101894221A (zh) * 2010-08-02 2010-11-24 北京航空航天大学 基于退化量分布非平稳时序分析的加速退化试验产品寿命预测方法
CN102033182A (zh) * 2010-12-10 2011-04-27 北京航空航天大学 一种固体钽电解电容器寿命预测方法
WO2012165300A1 (ja) * 2011-05-31 2012-12-06 株式会社日立製作所 蓄電システム
CN104778354A (zh) * 2015-04-03 2015-07-15 北京航空航天大学 定期抽检型产品贮存寿命评估方法
CN108280608A (zh) * 2017-12-27 2018-07-13 中国人民解放军63908部队 产品寿命分析方法及终端设备
CN108303317A (zh) * 2018-01-23 2018-07-20 中国兵器工业第五九研究所 一种橡胶密封圈失效检测方法
CN108446523A (zh) * 2018-05-11 2018-08-24 北京航天自动控制研究所 一种电子整机贮存寿命评估与预测方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101894221A (zh) * 2010-08-02 2010-11-24 北京航空航天大学 基于退化量分布非平稳时序分析的加速退化试验产品寿命预测方法
CN102033182A (zh) * 2010-12-10 2011-04-27 北京航空航天大学 一种固体钽电解电容器寿命预测方法
WO2012165300A1 (ja) * 2011-05-31 2012-12-06 株式会社日立製作所 蓄電システム
CN104778354A (zh) * 2015-04-03 2015-07-15 北京航空航天大学 定期抽检型产品贮存寿命评估方法
CN108280608A (zh) * 2017-12-27 2018-07-13 中国人民解放军63908部队 产品寿命分析方法及终端设备
CN108303317A (zh) * 2018-01-23 2018-07-20 中国兵器工业第五九研究所 一种橡胶密封圈失效检测方法
CN108446523A (zh) * 2018-05-11 2018-08-24 北京航天自动控制研究所 一种电子整机贮存寿命评估与预测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
LIU TIAN-YU，SUN QUAN，PAN ZHENG-QIANG，FENG JING: "Methods for Ｒeliability Assessment under Irregular Time-Varying-Stress", 《JOURNAL OF DONGHUA UNIVERSITY ( ENG．ED.)》 *
冯静: "基于秩相关系数的加速贮存退化失效机理一致性检验", 《航空动力学报》 *
李文华，周露等: "基于统计分析的继电器贮存寿命神经网络预测", 《航天控制》 *
王长安，牛跃听: "某型陀螺仪自然贮存寿命评估与加速寿命试验时间预估", 《中国惯性技术学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114091627A (zh) * 2022-01-20 2022-02-25 华谱科仪(北京)科技有限公司 色谱仪试验数据验证方法、装置、存储介质及电子设备
CN114091627B (zh) * 2022-01-20 2022-03-25 华谱科仪(北京)科技有限公司 色谱仪试验数据验证方法、装置、存储介质及电子设备

Also Published As

Publication number Publication date
CN110889187B (zh) 2020-12-11

Similar Documents

Publication Publication Date Title
CN109657937B (zh) 一种基于退化数据的产品可靠性评估与寿命预测方法
CN107451392B (zh) 一种含有多个相关退化过程的剩余寿命预测方法
CN108460230B (zh) 基于数据融合的腐蚀油气管道剩余寿命及可靠性预测方法
CN110174261B (zh) 多退化量监测的齿轮实时剩余寿命预测方法
CN111611294B (zh) 星敏感器数据异常检测方法
CN112883550A (zh) 一种考虑多重不确定性的退化设备剩余寿命预测方法
CN115962797A (zh) 一种基于温度应力下的传感器可靠性测试方法及系统
CN110889187B (zh) 基于Pearson系数的贮存退化数据一致性检验法
CN114091320B (zh) 一种天然气管道腐蚀失效时间预测方法及装置
CN110889077B (zh) 一种加速贮存与自然贮存退化数据一致性检验方法
CN110889189B (zh) 一种产品在不同试验应力水平下的退化数据一致性检验法
CN111428360A (zh) 一种不同维修深度分布下海底管道剩余寿命预测方法
CN110895624B (zh) 基于最大熵谱估计的加速贮存与自然贮存退化数据一致性检验法
CN108875169B (zh) 水面舰船装备用数字多用表的退化建模与寿命预测方法
CN110889186B (zh) 灰关联分析的加速贮存与自然贮存退化数据一致性检验法
CN109375143A (zh) 一种确定智能电能表剩余寿命的方法
CN111967140A (zh) 一种考虑混合不确定性的性能退化实验建模与分析方法
CN111783308A (zh) 一种精确预测隧道围岩位移的方法
CN116046415A (zh) 整车功能件可靠性台架验证方法、系统及整车验证台架
CN111967142B (zh) 一种考虑认知不确定性的性能退化实验建模与分析方法
CN110889083B (zh) 基于窗谱估计的退化数据一致性检验法
CN114492074A (zh) 一种概率损伤容限评估分析方法
CN114266006A (zh) 一种加速退化试验测量不确定度的评定方法
Louahem M'Sabah et al. Simulation of Bearing Degradation by the Use of the Gamma Stochastic Process.
CN117929171B (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
PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Consistency test method of storage degradation data based on Pearson coefficient

Effective date of registration: 20211008

Granted publication date: 20201211

Pledgee: Bank of Changsha Co.,Ltd. Yinde sub branch

Pledgor: Hunan gingko Reliability Technology Research Institute Co.,Ltd.

Registration number: Y2021980010281