CN106202725B - 基于多源数据的成败型产品可靠性试验设计和决策方法 - Google Patents
基于多源数据的成败型产品可靠性试验设计和决策方法 Download PDFInfo
- Publication number
- CN106202725B CN106202725B CN201610544423.9A CN201610544423A CN106202725B CN 106202725 B CN106202725 B CN 106202725B CN 201610544423 A CN201610544423 A CN 201610544423A CN 106202725 B CN106202725 B CN 106202725B
- Authority
- CN
- China
- Prior art keywords
- test
- success
- source data
- failure type
- distribution
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/32—Circuit design at the digital level
- G06F30/33—Design verification, e.g. functional simulation or model checking
- G06F30/3323—Design verification, e.g. functional simulation or model checking using formal methods, e.g. equivalence checking or property checking
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于多源数据的成败型产品可靠性试验设计和决策方法,包括以下几个步骤:步骤一、获取多源数据参数的修正因子;步骤二、建立试验评价模型;步骤三、试验的费用权衡;步骤四、试验方案及规划决策;本发明能够综合利用成败型产品的多源数据信息,并以此为依据制定合理的试验方案;由于能够扩充样本量,使得设计得到的试验方案能够以相对低廉的代价,取得预期的效果。
Description
技术领域
本发明是一种利用多源成败型数据,通过引入参数折合系数,对成败型产品的可靠性试验进行设计的方法,为产品的试验实施提供决策依据。以此充分利用、合理规划产品开发过程中的各个来源的数据,在降低试验成本的同时确保评估的精度,属于质量与可靠性领域。
背景技术
在生产生活中,成败型产品被广泛使用,如火箭、安全气囊、灭火器等等,对它们进行可靠性评估,主要是利用“成功”、“失败”等成败型数据;除成败型产品外,在考核复杂大***(如飞机)的任务可靠度等指标时,亦将任务“成功”与“失败”的次数数据作为评估的依据。因此就需要设计相应的可靠性试验,在要求的可信范围内获得并验证产品的可靠性指标,验证满足一定置信度的可靠性指标通常需要投入一定量的样本方可实现,而在产品的实际研制过程中,由于经费、时间等因素的限制,往往难以满足样本量的要求。另一方面,产品有时存在着其它来源的样本量,如历史数据、相似产品数据、其它类型试验数据等。若能在可靠性试验设计中考虑上述多源数据的作用,在此基础上所设计的试验及其试验结果可以协同多源数据综合考量指标,那么就可以在降低试验成本的同时确保评估的精度,为产品的开发提供有力支持。
发明内容
本发明的目的是为了解决上述问题,提出一种能够有效的利用多源数据,设计并规划相关试验,以此获得准确有效的评估结果。
本发明的具体步骤为:
步骤一、获取多源数据参数的修正因子;
步骤二、建立试验评价模型;
步骤三、试验的费用权衡;
步骤四、试验方案及规划决策。
本发明的优点在于:
(1)本发明能够综合利用成败型产品的多源数据信息,并以此为依据制定合理的试验方案;
(2)由于能够扩充样本量,使得设计得到的试验方案能够以相对低廉的代价,取得预期的效果。
附图说明
图1是本发明的流程图;
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
方法的流程图如图1所示,包括以下几个步骤:
步骤一、获取多源数据参数的修正因子
已知在开展可靠性试验前,搜集到m个来源的产品成败型数据,每个来源的成败型数据可表示为(ri,ni),i=1…m,其中ni为样本数,ri为成功数,它服从二项分布为:
ri~B(ni,pi) (1)
若投入可靠性试验的样本数为n0,成功数为r0,那么其二项分布可表示为:
r0~B(n0,p0) (2)
因此引入修正因子ki来描述多源数据间的关系,可表示为:
pi=kip0(3)
修正因子ki可以通过贝叶斯方法进行求解,也可以通过专家评分法、类比分析法等手段获取。
步骤二、建立试验评价模型
根据二项分布的性质可知,
r0=n0·p0 (4)
亦即
ri=ni·pi=ni·ki·p0 (5)
也就是说,
ri/ni=ki·r0/n0,ri=r0·ki·ni/n0
根据二项分布的性质,在置信度γ已知的情况下,可靠度的置信下限RL为
其中
二项分布的共轭先验分布为Beta分布,a,b为Beta分布中的超参数。根据共轭先验分布的性质:
a=r,b=n-r
那么在p0条件下,超参数可以写作:
a=(ri·n0)/(ki·ni),b=ni-a,
因此对于多源数据而言,可以得到:
那么公式(7)可以写作:
那么公式(6)和(8)即是试验评价模型,可根据试验评价模型来计算满足一定评估精度时所需投入的样本量n0。
步骤三、试验的费用权衡
通过公式(6)和(8)可以预计满足一定评估精度时所需投入的样本量n0。但在设计试验时,除了要考虑样本量还应考虑经费C这一约束条件,
C=(Cs+Ct)·n0+CLoss·Prisk
其中Cs是每个样本的费用,Ct是每个样本所需的试验费用,Prisk是试验方案存在的风险,CLoss是试验方案失败所造成的额外损失。
Prisk是一个同样本量n0,评估结果RL,置信度γ以及修正因子ki相关的函数,可根据上述参数的情况权衡确定。
步骤四、试验方案及规划决策
试验方案以评估结果RL误差最小化为目标,以试验费用为C为约束,再根据需要考虑其它因素,采用优化算法进行优化设计。
在设计多次试验的方案时,可以根据可靠性要求的增长与产品实际的可靠性增长情况,设计每次试验需投入的样本量nj,每次RL的计算公式为:
实施例:
现通过对某***的应用,对本专利的方法进行说明。
步骤一、获取多源数据参数的修正因子
已知某产品在进行可靠性试验前开展了3次试验,试验数据分别为(81,84),(115,125),(101,113),通过分析得到其对应的修正因子分别为k1=0.894,k2=0.956,k3=0.994。
步骤二、建立试验评价模型
通过公式(6)和(8)可以计算得到要满足γ=0.9,RL=0.865的条件,需要最少投入的样本量为23,其中需成功数为21。
注:若采用常规方法,在高风险的情况下也要最少投入51个样本。
步骤三、试验的费用权衡
已知每个产品的费用Cs为2万元,每次试验费用Ct为1万元,CLoss为100万,计算得到不同样本量对应的风险和费用分别为:
表1试验方案权衡
试验方案 | A | B | C | D | E | F | G | H | I |
样本量n<sub>0</sub> | 23 | 24 | 25 | 26 | 27 | 28 | 29 | 30 | 31 |
风险值P<sub>risk</sub> | 25% | 20% | 16% | 13% | 9% | 8% | 6% | 5% | 5% |
费用C(万) | 94 | 92 | 91 | 91 | 90 | 92 | 93 | 95 | 98 |
步骤四、试验方案及规划决策
在满足评估精度的方案中(样本量大于23个),若将试验费用作为权衡约束,那么可确定选择试验方案E。
Claims (1)
1.基于多源数据的成败型产品可靠性试验设计和决策方法,其特征在于,包括以下几个步骤:
步骤一、获取多源数据参数的修正因子
已知在开展可靠性试验前,搜集到m个来源的产品成败型数据,第i个来源的成败型数据可表示为(ri,ni),i=1…m,其中ni为样本数,ri为成功数,它服从二项分布为:
ri~B(ni,pi), (1)
若投入可靠性试验的样本数为n0,成功数为r0,那么其二项分布可表示为:
r0~B(n0,p0), (2)
因此引入修正因子ki来描述第i个来源数据与可靠性试验间的成功率pi和p0的关系,可表示为:
pi=kip0, (3)
修正因子ki可以通过贝叶斯方法求解,也可以通过专家评分法、类比分析法获取;
步骤二、建立试验评价模型
根据二项分布的性质可知,
r0=n0·p0 (4)
亦即
ri=ni·pi=ni·ki·p0 (5)
也就是说,
ri/ni=ki·r0/n0,ri=r0·ki·ni/n0
在置信度γ已知的情况下,可靠度的置信下限RL为
根据二项分布的共轭先验分布性质,先验分布为Beta分布,a,b为Beta分布中的超参数;那么:
且:
a=r,b=n-r,
那么在p0条件下,超参数可以写作:
a=(ri·n0)/(ki·ni),b=ni-a,
因此对于m个来源数据而言,可以得到:
那么公式(7)可以写作:
那么公式(6)和(8)即是试验评价模型,可根据试验评价模型来计算满足一定评估精度时所需投入的样本量n0;
步骤三、试验的费用权衡
通过公式(6)和(8)可以预计满足一定评估精度时所需投入的样本量n0;但在设计试验时,除了要考虑样本量还应考虑经费C这一约束条件,
C=(Cs+Ct)·n0+CLoss·Prisk
其中Cs是每个样本的费用,Ct是每个样本所需的试验费用,Prisk是试验方案存在的风险,CLoss是试验方案失败所造成的额外损失;
Prisk是一个同样本量n0,评估结果RL,置信度γ以及修正因子ki相关的函数,可根据上述参数的情况权衡确定;
步骤四、试验方案及规划决策
试验方案以评估结果RL误差最小化为目标,以试验费用为C为约束,再根据需要考虑其它因素,采用优化算法进行优化设计;
在设计多次试验的方案时,可以根据可靠性要求的增长与产品实际的可靠性增长情况,设计第j次试验需投入的样本量nj,第j次的计算公式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610544423.9A CN106202725B (zh) | 2016-07-12 | 2016-07-12 | 基于多源数据的成败型产品可靠性试验设计和决策方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610544423.9A CN106202725B (zh) | 2016-07-12 | 2016-07-12 | 基于多源数据的成败型产品可靠性试验设计和决策方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106202725A CN106202725A (zh) | 2016-12-07 |
CN106202725B true CN106202725B (zh) | 2021-03-12 |
Family
ID=57476432
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610544423.9A Active CN106202725B (zh) | 2016-07-12 | 2016-07-12 | 基于多源数据的成败型产品可靠性试验设计和决策方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106202725B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109086290A (zh) * | 2018-06-08 | 2018-12-25 | 广东万丈金数信息技术股份有限公司 | 基于多源数据决策树的注册信息真实性判断方法及*** |
CN109508455B (zh) * | 2018-10-18 | 2021-11-19 | 山西大学 | 一种GloVe超参数调优方法 |
CN110134090B (zh) * | 2019-05-23 | 2020-09-29 | 电子科技大学 | 融合多源信息的工业机器人控制***可靠性评估方法 |
CN112883497B (zh) * | 2021-03-22 | 2023-06-16 | 中国人民解放军国防科技大学 | 基于多源信息融合的航天阀门可靠性评估方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8024128B2 (en) * | 2004-09-07 | 2011-09-20 | Gene Security Network, Inc. | System and method for improving clinical decisions by aggregating, validating and analysing genetic and phenotypic data |
CN102360332A (zh) * | 2011-09-28 | 2012-02-22 | 北京航空航天大学 | 一种软件可靠性加速测试与评估方法及其计算机辅助工具 |
CN103646138A (zh) * | 2013-12-03 | 2014-03-19 | 北京航空航天大学 | 基于贝叶斯理论的定时截尾加速验收抽样试验优化设计方法 |
CN104077445A (zh) * | 2014-07-01 | 2014-10-01 | 北京航空航天大学 | 基于模糊理论的加速寿命试验统计分析方法 |
CN104133994A (zh) * | 2014-07-24 | 2014-11-05 | 北京航空航天大学 | 融合多源成败型数据的可靠度评估方法 |
WO2015014295A1 (en) * | 2013-07-31 | 2015-02-05 | Tencent Technology (Shenzhen) Company Limited | Multi-source download scheduling |
CN104933323A (zh) * | 2015-07-10 | 2015-09-23 | 北京航空航天大学 | 融合产品成败型数据和故障时间数据的可靠性评估方法 |
-
2016
- 2016-07-12 CN CN201610544423.9A patent/CN106202725B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8024128B2 (en) * | 2004-09-07 | 2011-09-20 | Gene Security Network, Inc. | System and method for improving clinical decisions by aggregating, validating and analysing genetic and phenotypic data |
CN102360332A (zh) * | 2011-09-28 | 2012-02-22 | 北京航空航天大学 | 一种软件可靠性加速测试与评估方法及其计算机辅助工具 |
WO2015014295A1 (en) * | 2013-07-31 | 2015-02-05 | Tencent Technology (Shenzhen) Company Limited | Multi-source download scheduling |
CN103646138A (zh) * | 2013-12-03 | 2014-03-19 | 北京航空航天大学 | 基于贝叶斯理论的定时截尾加速验收抽样试验优化设计方法 |
CN104077445A (zh) * | 2014-07-01 | 2014-10-01 | 北京航空航天大学 | 基于模糊理论的加速寿命试验统计分析方法 |
CN104133994A (zh) * | 2014-07-24 | 2014-11-05 | 北京航空航天大学 | 融合多源成败型数据的可靠度评估方法 |
CN104933323A (zh) * | 2015-07-10 | 2015-09-23 | 北京航空航天大学 | 融合产品成败型数据和故障时间数据的可靠性评估方法 |
Non-Patent Citations (1)
Title |
---|
A Bayesian reliability evaluation method with integrated accelerated degradation testing and field information;lizhi wang等;《Reliability Engineering and System Safety》;20121129;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN106202725A (zh) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106202725B (zh) | 基于多源数据的成败型产品可靠性试验设计和决策方法 | |
CN102629300A (zh) | 一种基于灰色预测模型的步进应力加速退化数据评估方法 | |
CN104133994A (zh) | 融合多源成败型数据的可靠度评估方法 | |
CN105608263A (zh) | 一种面向涡轮叶盘结构寿命概率分析的自适应处理方法 | |
Luo et al. | Accelerated reliability demonstration under competing failure modes | |
Gail et al. | Calibration and seasonal adjustment for matched case–control studies of vitamin D and cancer | |
Lehtinen | Derivational robustness and indirect confirmation | |
CN103955892A (zh) | 一种目标跟踪方法及扩展截断无迹卡尔曼滤波方法、装置 | |
CN113627098A (zh) | 一种cfd模型确认方法及产品设计方法 | |
CN104820771A (zh) | 一种航天工程制造成熟度等级确定方法 | |
Zheng et al. | Robust estimation in joint mean–covariance regression model for longitudinal data | |
CN103106332B (zh) | 一种测量不确定度的分析方法 | |
CN110850358A (zh) | 基于逐步回归算法的电能表综合检定方法及*** | |
David et al. | Resistant measures in assessing the adequacy of regression models | |
Lin et al. | Estimation of P (Y< X) for location-scale distributions under joint progressively type-II right censoring | |
Cao et al. | Dynamic prediction and compensation of aerocraft assembly variation based on state space model | |
CN110096820B (zh) | 单元寿命服从威布尔分布时k/n(G)***的剩余寿命预测方法 | |
Jin et al. | Transfer learning with quantile regression | |
CN109598059B (zh) | 一种基于代理模型的热防护***优化设计方法及设计*** | |
CN105389406B (zh) | 基于故障的单位加权累计数的整车设计可靠性评估方法 | |
CN113722961A (zh) | 一种基于广义协同高斯过程模型的结构不确定性量化解析方法 | |
CN102306353A (zh) | 一种仿真***可信性评估方法及*** | |
CN102306249B (zh) | 一种设备效能仿真方法及*** | |
Kim et al. | Testing for an excessive number of zeros in time series of bounded counts | |
Plotkin et al. | Detailed weather and terrain analysis for aircraft noise modeling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |