CN111914386A - 一种基于退化模型不确定分析的可靠性评估方法及系统 - Google Patents
一种基于退化模型不确定分析的可靠性评估方法及系统 Download PDFInfo
- Publication number
- CN111914386A CN111914386A CN202010332096.7A CN202010332096A CN111914386A CN 111914386 A CN111914386 A CN 111914386A CN 202010332096 A CN202010332096 A CN 202010332096A CN 111914386 A CN111914386 A CN 111914386A
- Authority
- CN
- China
- Prior art keywords
- analysis
- uncertain
- product
- reliability
- model
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
为解决现有技术的不足,本发明提供了一种基于退化模型不确定分析的可靠性评估方法,包括以下步骤:S1.确定产品性能参数的失效阈值的不确定分布。S2.采用时变不确定分布函数建立目标产品的退化模型。S3.建立产品的可靠性模型。S4.建立产品的性能裕量方程,获取性能裕量的逆不确定分布函数,进行可靠性评估。相比现有基于概率论的退化数据分析方法,本发明在分析评估过程中引入了失效阈值的不确定性,对于更关键的产品寿命前期可进行更为准确的评估,从而获得更为准确的评估结果。此外本发明还提供了一种基于退化模型不确定分析的可靠性评估系统,其中分析模块采用基于退化模型不确定分析的可靠性评估方法进行分析。
Description
技术领域
本发明属于可靠性工程技术领域,具体涉及一种基于退化模型不确定分析的可靠性评估方法及系统。
背景技术
传统的产品可靠性评估方法是基于对产品失效时间数据的建模与统计分析,然而对一些高可靠性产品,其失效时间数据很难获取。随着产品失效物理分析技术与性能测试方法的进步,基于性能退化数据分析的可靠性评估方法被越来越多应用。根据测试的产品性能退化数据能够建立起产品性能的随使用时间的退化模型,不需要产品失效,即可根据退化模型计算出产品退化到失效阈值的时间,从而实现高效的可靠性评估。
目前为止,绝大多数退化数据分析方法都是基于概率论的方法,难以表征退化数据分析过程中的认知不确定性,造成分析结果不可信。特别是对于很多新研的产品,产品的失效阈值往往是不确定的,此外,测量到的产品性能退化数据也未必准确,不可避免带有不确定性。
发明内容
本发明针对现有技术存在的不足,提供了一种基于退化模型不确定分析的可靠性评估方法,包括以下步骤:
S1.确定产品性能参数的失效阈值Pth的不确定分布。
S2.采用时变不确定分布函数建立目标产品的退化模型。
S3.建立产品的可靠性模型。
S4.建立产品的性能裕量方程,获取性能裕量的逆不确定分布函数,进行可靠性评估。
进一步的,步骤S1包括:
S1-2基于Pth的预估范围:Pth值小于y1,y2,…,yk的信度为α1,α2,…,αk,其中 y1<y2,…,<yk,α1<α2,…,<αk,建立如下最小二乘估计模型
进一步的,步骤S2包括:
S2-1设Pi,j,ti,j分别为第i个产品的第j次性能退化测量值与测试时刻, i=1,2,…,N1,j=1,2,…,N2,假定各测量时刻ti,j对应的Pi,j都服从不确定分布估计得ti,j对应的参数估计值
S2-4建立如下似然方程
进一步的,步骤S3包括:
S3-1建立产品可靠性模型为
R(t)=M(E(t)>0)
式中,M(·)表示不确定测度函数。
S3-2由不确定理论的对偶公理,得
R(t)=M(E(t)>0)=1-M(E(t)≤0)
S3-3根据不确定分布的定义,得
R(t)=1-M(E(t)≤0)=1-Ψt(0)。
进一步的,步骤S4包括:
S4-1根据P值望大或望小,确定性能裕量方程E(t)。望小是指P值越小对产品性能越有利,望大是指P值越大对产品性能越有利。
S4-2根据如下运算法则推导Ψt(x)对应的逆分布函数
S4-4:将Ψt(0)代入步骤S3建立的可靠性模型中,得到
S4-5:计算产品的可靠寿命为T(α)={t|R(t)=α}={t|Ψt(0)=1-α}。
进一步的,步骤S4-1所述性能裕量方程为:
式中,E(t)为产品的性能裕量,P(t)为性能退化量,Pth为失效阈值。
此外,本发明还提供了一种基于退化模型不确定分析的可靠性评估系统,包括:输入模块、分析模块和输出模块。
所述输入模块向分析模块输入数据。所述分析模块基于输入模块输入的数据,采用任一上述基于退化模型不确定分析的可靠性评估方法进行可靠性评估分析。所述输出模块将分析模块的分析结果输出反馈。
本发明的优点在于:相比现有基于概率论的退化数据分析方法,本发明在分析评估过程中引入了失效阈值的不确定性,对于更关键的产品寿命前期可进行更为准确的评估,从而获得更为准确的评估结果。
附图说明
图1所示为本发明实施例2各测量时刻的ej估计值。
图2所示为本发明实施例2各测量时刻的σj估计值。
图3所示为本发明实施例2可靠度曲线对比示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下对本发明进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例1
一种基于退化模型不确定分析的可靠性评估方法,包括以下步骤:
S1.确定产品性能参数的失效阈值Pth的不确定分布。
S2.采用时变不确定分布函数建立目标产品的退化模型。
S3.建立产品的可靠性模型。
S4.建立产品的性能裕量方程,获取性能裕量的逆不确定分布函数,进行可靠性评估。
步骤S1包括:
S1-2虽然产品性能参数的失效阈值Pth是不确定的,但是产品设计师能够给出失效阈值的大概范围,例如产品设计师分别给出了Pth不小于值y1,y2,…,yk的信度为α1,α2,…,αk,其中y1<y2,…,<yk,α1<α2,…,<αk。根据信息,建立如下最小二乘估计模型
步骤S2包括:
S2-1设Pi,j,ti,j分别为第i个产品的第j次性能退化测量值与测试时刻, i=1,2,…,N1,j=1,2,…,N2,假定各测量时刻ti,j对应的Pi,j都服从不确定正态分布估计得ti,j对应的参数估计值
S2-4建立如下似然方程
步骤S3包括:
S3-1建立产品可靠性模型为
R(t)=M(E(t)>0)
式中,M(·)表示不确定测度函数。
S3-2由不确定理论的对偶公理,得
R(t)=M(E(t)>0)=1-M(E(t)≤0)
S3-3根据不确定分布的定义,得
R(t)=1-M(E(t)≤0)=1-Ψt(0)。
步骤S4包括:
S4-1根据P值望大或望小,确定性能裕量方程E(t)。望小是指P值越小对产品性能越有利,望大是指P值越大对产品性能越有利。为了度量产品的可靠度需要建立产品的裕量方程,裕量为正E(t)>0意味着产品可靠。
S4-2根据如下运算法则推导Ψt(x)对应的逆分布函数
S4-4:将Ψt(0)代入步骤S3建立的可靠性模型中,得到
S4-5:计算产品的可靠寿命为T(α)={t|R(t)=α}={t|Ψt(0)=1-α}。
步骤S4-1所述性能裕量方程为:
式中,E(t)为产品的性能裕量,P(t)为性能退化量,Pth为失效阈值。
实施例2
基于实施例1所述基于退化模型不确定分析的可靠性评估方法,以GaAs 激光器在80℃下的退化实验数据为例,性能退化参数为工作电流相对于原始值的百分比增量,共有10个样品在16个时刻的退化量观测数据。
表1 GaAs激光器的性能参数退化量Pi,j
基于表1数据进行GaAs激光器的可靠性评估:
首先,基于从产品设计工程师出获取的信息:Pth值小于8的信度为0.3,小于10的信度为0.7,小于12的信度为0.95。首先假定Pth服从不确定线性分布,利用最小二乘法估计参数值,得到Pth~L(6.000,12.154),其残差平方和为∑1=3.7E-3。然后假定Pth服从不确定正态分布,利用最小二乘法估计参数值,得到Pth~N(8.987,2.509),其残差平方和为∑2=3.6E-4。由于∑2值更小,据此相信Pth更好得服从不确定正态分布Pth~N(8.987,2.509)。
通过步骤S3建立可靠性模型后,将参数估计值代入步骤S4-4式中,确定出产品可靠性预测模型为
可靠度曲线如图3中的实线,计算出T(0.9)=1.766×103h。
由Ω-1(α)可推导出
进而确定出可靠性预测模型为
式中,设值为不确定变量Pth的期望值,可靠度曲线如图3 中的虚线,计算出T(0.9)*=2.668×103h。从R(t)与R*(t)的可靠度曲线可见,是否考虑失效阈值的不确性导致的可靠度评估结果具有显著不同。如果没有考虑失效阈值的不确定性,在更为关键的产品寿命前期获得的可靠度评估结果明显偏大,得出的可靠寿命值也明显偏大,例如T(0.9)*比T(0.9)大50%以上。因此,为了更为准确地获得可靠度评估结果,应该充分考虑失效阈值的不确定性。
实施例3
一种基于退化模型不确定分析的可靠性评估系统,包括:输入模块、分析模块和输出模块。
所述输入模块向分析模块输入数据。所述分析模块基于输入模块输入的数据,采用实施例1基于退化模型不确定分析的可靠性评估方法进行可靠性评估分析。所述输出模块将分析模块的分析结果输出反馈。
以实施例2中GaAs激光器在80℃下的退化实验数据为例,将表1数据通过输入模块输入至分析模块中,分析模块进行如实施例2所述分析后,所得到的数据通过输出模块反馈给用户。以图形化反馈为例,可产生如图3实线所示图形结果。
应该注意到并理解,在不脱离本发明权利要求所要求的精神和范围的情况下,能够对上述详细描述的本发明做出各种修改和改进。因此,要求保护的技术方案的范围不受所给出的任何特定示范教导的限制。
Claims (8)
1.一种基于退化模型不确定分析的可靠性评估方法,其特征在于,包括以下步骤:
S1.确定产品性能参数的失效阈值Pth的不确定分布;
S2.采用时变不确定分布函数建立目标产品的退化模型;
S3.建立产品的可靠性模型;
S4.建立产品的性能裕量方程,获取性能裕量的逆不确定分布函数,进行可靠性评估。
3.根据权利要求1所述基于退化模型不确定分析的可靠性评估方法,其特征在于,步骤S2包括:
S2-1设Pi,j,ti,j分别为第i个产品的第j次性能退化测量值与测试时刻,i=1,2,…,N1,j=1,2,…,N2,假定各测量时刻ti,j对应的Pi,j都服从不确定分布估计得ti,j对应的参数估计值
S2-4建立如下似然方程
5.根据权利要求1所述基于退化模型不确定分析的可靠性评估方法,其特征在于,步骤S3包括:
S3-1建立产品可靠性模型为
R(t)=M(E(t)>0)
式中,M(·)表示不确定测度函数;
S3-2由不确定理论的对偶公理,得
R(t)=M(E(t)>0)=1-M(E(t)≤0)
S3-3根据不确定分布的定义,得
R(t)=1-M(E(t)≤0)=1-Ψt(0)。
8.一种基于退化模型不确定分析的可靠性评估系统,其特征在于,包括:输入模块、分析模块和输出模块;
所述输入模块向分析模块输入数据;所述分析模块基于输入模块输入的数据,采用权利要求1-7任一基于退化模型不确定分析的可靠性评估方法进行可靠性评估分析;所述输出模块将分析模块的分析结果输出反馈。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010332096.7A CN111914386A (zh) | 2020-04-24 | 2020-04-24 | 一种基于退化模型不确定分析的可靠性评估方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010332096.7A CN111914386A (zh) | 2020-04-24 | 2020-04-24 | 一种基于退化模型不确定分析的可靠性评估方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111914386A true CN111914386A (zh) | 2020-11-10 |
Family
ID=73237409
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010332096.7A Pending CN111914386A (zh) | 2020-04-24 | 2020-04-24 | 一种基于退化模型不确定分析的可靠性评估方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111914386A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112834370A (zh) * | 2021-01-07 | 2021-05-25 | 北京航空航天大学 | 航天机构产品确信可靠性退化方程的建立方法 |
CN114580118A (zh) * | 2022-03-25 | 2022-06-03 | 兰州空间技术物理研究所 | 一种离子推力器寿命及可靠性定量评估方法 |
-
2020
- 2020-04-24 CN CN202010332096.7A patent/CN111914386A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112834370A (zh) * | 2021-01-07 | 2021-05-25 | 北京航空航天大学 | 航天机构产品确信可靠性退化方程的建立方法 |
CN114580118A (zh) * | 2022-03-25 | 2022-06-03 | 兰州空间技术物理研究所 | 一种离子推力器寿命及可靠性定量评估方法 |
CN114580118B (zh) * | 2022-03-25 | 2023-12-15 | 兰州空间技术物理研究所 | 一种离子推力器寿命及可靠性定量评估方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111046564B (zh) | 两阶段退化产品的剩余寿命预测方法 | |
CN110851980B (zh) | 一种设备剩余寿命预测方法及系统 | |
Jia et al. | Quality-related fault detection approach based on dynamic kernel partial least squares | |
CN107153874B (zh) | 水质预测方法及系统 | |
CN108647272B (zh) | 一种基于数据分布的小样本扩充对脱丁烷塔底丁烷浓度进行预测的方法 | |
CN111123188A (zh) | 基于改进最小二乘法的电能表综合检定方法及系统 | |
CN109739214B (zh) | 工业过程间歇故障的检测方法 | |
CN109446189A (zh) | 一种工业参数离群点检测系统及方法 | |
CN109917777B (zh) | 基于混合多采样率概率主成分分析模型的故障检测方法 | |
CN111638707B (zh) | 基于som聚类和mpca的间歇过程故障监测方法 | |
CN108829878B (zh) | 一种工业实验数据异常点检测方法及装置 | |
CN111914386A (zh) | 一种基于退化模型不确定分析的可靠性评估方法及系统 | |
Zhou et al. | A control chart based on likelihood ratio test for detecting patterned mean and variance shifts | |
CN109670549B (zh) | 火电机组的数据筛选方法、装置以及计算机设备 | |
Chen et al. | Two-phase degradation data analysis with change-point detection based on Gaussian process degradation model | |
CN110084301B (zh) | 一种基于隐马尔可夫模型的多工况过程工况辨识方法 | |
CN113791351B (zh) | 基于迁移学习和差值概率分布的锂电池寿命预测方法 | |
CN108537249B (zh) | 一种密度峰值聚类的工业过程数据聚类方法 | |
CN112149054B (zh) | 基于时序扩展的正交邻域保持嵌入模型的构建与应用 | |
CN110750848A (zh) | 一种考虑软件运行的软-硬件退化系统剩余寿命估计的方法 | |
CN112069045A (zh) | 基于隐马尔可夫的云平台软件性能预测方法 | |
WO2024103470A1 (zh) | 自动扶梯生产用性能检测系统及其方法 | |
Rizvi et al. | Real-time zip load parameter tracking using sensitivity-based adaptive window and variable elimination with realistic synchrophasor data | |
CN116048036A (zh) | 基于动态非平稳投影结构的故障检测方法 | |
CN114064760A (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 |