CN101847834A - 电站小样本及无失效数据部件失效率的可靠性评估方法 - Google Patents

电站小样本及无失效数据部件失效率的可靠性评估方法 Download PDF

Info

Publication number
CN101847834A
CN101847834A CN201010192847A CN201010192847A CN101847834A CN 101847834 A CN101847834 A CN 101847834A CN 201010192847 A CN201010192847 A CN 201010192847A CN 201010192847 A CN201010192847 A CN 201010192847A CN 101847834 A CN101847834 A CN 101847834A
Authority
CN
China
Prior art keywords
lambda
failure rate
centerdot
power station
alpha
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
CN201010192847A
Other languages
English (en)
Other versions
CN101847834B (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.)
China General Nuclear Power Corp
Suzhou Nuclear Power Research Institute Co Ltd
Original Assignee
China General Nuclear Power Corp
Suzhou Nuclear Power 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 China General Nuclear Power Corp, Suzhou Nuclear Power Research Institute Co Ltd filed Critical China General Nuclear Power Corp
Priority to CN201010192847.6A priority Critical patent/CN101847834B/zh
Publication of CN101847834A publication Critical patent/CN101847834A/zh
Application granted granted Critical
Publication of CN101847834B publication Critical patent/CN101847834B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • General Factory Administration (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种电站小样本及无失效数据部件失效率的可靠性评估方法,该方法包括如下步骤:a)判断可用失效事件数据数量是否大于等于n,当判断结果为否时进入步骤b),否则根据双参数威布尔分布确定电站部件的失效率;b)判断可用失效事件数据数量是否大于0,当判断结果为否时进入步骤c),否则利用bayes方法计算部件失效率;c)判断是否存在通用失效率,当判断结果为否时进入步骤d),否则利用电站运行数据对通用失效率进行修正;d)利用减函数方法计算部件失效率。本发明可以快速有效的进行部件失效评估,有助于提高管理水平和促进产品质量提高,为电站的生产运行建立安全可靠的评估机制。

Description

电站小样本及无失效数据部件失效率的可靠性评估方法
技术领域
本发明涉及一种电站小样本及无失效数据部件失效率的可靠性评估方法。
背景技术
金属构件(包括机械零部件、金属结构件及其他金属制品)在使用过程中失去原有功能(或达不到原有功能)的现象称为金属构件失效,包括金属构件变形、断裂(开裂)、腐蚀与磨损,失效分析通常是指对失效产品为寻找失效原因和预防措施所进行的一切技术活动,就是研究失效现象的特征和规律,从而找出失效的模式和原因。对于各类电站来说,工作部件一旦失效会产生巨大的经济损失乃至社会损失,因此,电站部件失效分析是一门综合性的质量系统工程,是一门解决材料、工程结构、系统组元等质量问题的工程学。它的任务是既要揭示产品功能失效的模式和原因,弄清失效的机理和规律,又要找出纠正和预防失效的措施,能够对现有金属构件的可靠性进行评估。失效分析具有巨大的社会经济效益,有助于提高管理水平和促进产品质量提高,有助于分清责任和保护生产者利益,是修订产品技术规范及标准的依据,为新产品开发提供依据,促进材料、力学、机械等学科的发展。
发明内容
为克服上述缺点,本发明的目的在于提供一种快速有效的电站小样本及无失效数据部件失效率的可靠性评估方法。
为了达到以上目的,本发明提供了一种电站小样本及无失效数据部件失效率的可靠性评估方法,该方法包括如下步骤:
a)判断可用失效事件数据数量是否大于等于n,当判断结果为否时进入步骤b),否则根据双参数威布尔分布确定电站部件的失效率;
b)判断可用失效事件数据数量是否大于0,当判断结果为否时进入步骤c),否则利用bayes方法计算部件失效率;
c)判断是否存在通用失效率,当判断结果为否时进入步骤d),否则利用电站运行数据对通用失效率进行修正;
d)利用减函数方法计算部件失效率。
优选地,在步骤a)中当失效事件数据数量大于等于n时,利用双参数威布尔分布模型计算部件失效率函数、故障分布函数、可靠度函数的点估计:
λ ( t ) = βt β - 1 α β , t ≥ 0 F ( t ) = 1 - exp [ - t β α ] , t ≥ 0 R ( t ) = exp [ - t β α ] , t ≥ 0 .
优选地,在步骤b)中当失效事件数据数量大于0时,
试验(r,T)似然函数:
λ取共轭型先验:f0(λ)=Γ(λ/r0,T0)
后验:f0(λ/r,T)=Γ(λ/f0+r,T0+T)
失效率单侧置信上限:
Figure BSA00000145400100025
优选地,在步骤c)中当存在通用失效率时,将通用失效率作为较保守的上限λ0,部件的寿命服从指数分布,有n个产品进行定时截尾试验,截尾时间段[0,ti],没有部件失效,则无失效时的修正似然函数为:
Figure BSA00000145400100026
C为常数
按Jeffreys准则,λ的先验分布取πλ∝λ-1,0<λ<λ0,在平方损失下的bayes估计为:
Figure BSA00000145400100027
优选地,在步骤d)中,采用无失效数据处理方法获得部件可靠度对于寿命服从指数分布的产品进行n次定时截尾试验,结果所有试验品无一失效,记试验总时间为T=nt0
取π(λ)∝λα-1,α∈(0,1)为先验分布,则:
h ( λ / T ) = L ( t 1 , t 2 , · · · t r / λ ) · π ( λ ) ∫ 0 ∞ L ( t 1 , t 2 , · · · t r / λ ) · π ( λ ) dλ
= λ α - 1 · exp ( - n t 0 λ ) ∫ 0 ∞ λ a - 1 · exp ( - n t 0 λ ) dλ = ( n t 0 ) a Γ ( α ) λ α - 1 exp ( - n t 0 λ )
λ的条件期望值:
λ ^ = ∫ 0 ∞ λh ( λ / nt 0 ) dλ = Γ ( a + 1 ) T 0 = a nt 0 a = 1 2 χ α 2 ( 2 ) = - ln α
因此,取一定的置信度求得部件失效率:
λ ^ = - ln α T 0 .
更进一步地,n大于等于5。
由于采用了以上技术方案,使得本发明可以快速有效的进行部件失效率可靠性评估,有助于提高管理水平和促进产品质量提高,有助于分清责任和保护生产者利益,是修订产品技术规范及标准的依据,为电站的生产运行建立安全可靠的评估机制。
附图说明
附图1为本发明电站小样本及无失效数据部件失效率的可靠性评估方法的流程图。
具体实施方式
下面结合附图对本发明的较佳实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
参见附图1所示,一种电站小样本及无失效数据部件失效率的可靠性评估方法,该方法包括如下步骤:
a)判断可用失效事件数据数量是否大于等于n,当判断结果为否时进入步骤b),否则根据双参数威布尔分布确定电站部件的失效率;
b)判断可用失效事件数据数量是否大于0,当判断结果为否时进入步骤c),否则利用bayes方法计算部件失效率;
c)判断是否存在通用失效率,当判断结果为否时进入步骤d),否则利用电站运行数据对通用失效率进行修正;
d)利用减函数方法计算部件失效率。
在步骤a)中当失效事件数据数量大于等于n时,也就是可用失效的数据较多情形下利用双参数威布尔分布模型计算部件失效率函数、故障分布函数、可靠度函数的点估计:
λ ( t ) = βt β - 1 α β , t ≥ 0 F ( t ) = 1 - exp [ - t β α ] , t ≥ 0 R ( t ) = exp [ - t β α ] , t ≥ 0 .
若取ln[-lnR(t)]为纵轴,lnt为横轴,则以上方程可表达为一条直线,斜率为β,纵截距为(-βlnα),这样用做图法求出α与β,MTTF与可靠度R(t)点估计为:
MTTF = ∫ 0 ∞ tf ( t ) dt = ∫ 0 ∞ R ( t ) dt = αΓ ( 1 + β β )
R ( t 0 ) = exp ( - ( t 0 α ) β )
在步骤b)中当失效事件数据数量大于0时,也就是说当无通用失效率且经验的反馈中可用失效事件较少时,外部经验反馈同样得到较少数据,在此基础上根据电站内部数对外部数据据进行Bayes修正:
试验(r,T)似然函数:f0(λ)=Γ(λ/r0,T0)
λ取共轭型先验:
Figure BSA00000145400100047
则后验:f0(λ/r,T)=Γ(λ/r0+r,T0+T)
采用无信息先验可以推出失效率单侧置信上限:
f 0 ( λ ) = Γ ( λ / 0,0 ) ∝ λ - 1 ⇒ λ u = χ 2 ( 2 r , C ) 2 T .
在步骤c)中当存在通用失效率时,对于部件有通用失效率但是经验反馈中无可用失效事情情形,对于指数分布在无失效的情形下,失效率不会很大,将通用失效率作为较保守的上限λ0,部件的寿命服从指数分布,有n个产品进行定时截尾试验,截尾时间段[0,ti],没有部件失效,则无失效时的修正似然函数为:
Figure BSA00000145400100049
C为常数
按Jeffreys准则,λ的先验分布取πλ∝λ-1,0<λ<λ0,在平方损失下的bayes估计为:
Figure BSA00000145400100051
在步骤d)中,也就是无任何失效数据部件的可靠性评估主要采用减函数方法采用无失效数据处理方法获得部件可靠度对于寿命服从指数分布的产品进行n次定时截尾试验,结果所有试验品无一失效,记试验总时间为T=nt0:取π(λ)∝λα-1,α∈(0,1)为先验分布,则:
h ( λ / T ) = L ( t 1 , t 2 , · · · t r / λ ) · π ( λ ) ∫ 0 ∞ L ( t 1 , t 2 , · · · t r / λ ) · π ( λ ) dλ
= λ α - 1 · exp ( - n t 0 λ ) ∫ 0 ∞ λ a - 1 · exp ( - n t 0 λ ) dλ = ( n t 0 ) a Γ ( α ) λ α - 1 exp ( - n t 0 λ )
λ的条件期望值:
λ ^ = ∫ 0 ∞ λh ( λ / nt 0 ) dλ = Γ ( a + 1 ) T 0 = a nt 0 a = 1 2 χ α 2 ( 2 ) = - ln α
因此,取一定的置信度求得部件失效率:
λ ^ = - ln α T 0 .
在这里需要指出的是,在实际操作过程中n的选择一般大于等于5,而在不同的情况下,n值的大小也可以根据需要选择不同的数值。
以上结合实施方式对本发明做了详细说明,只为说明本发明的技术构思及特点,其目的在于让熟悉此项技术的人了解本发明的内容并加以实施,并不能以此限定本发明的保护范围,凡根据本发明精神实质所做的等效变化或修饰,都应涵盖在本发明的保护范围之内。

Claims (6)

1.一种电站小样本及无失效数据部件失效率的可靠性评估方法,其特征在于:该方法包括如下步骤:
a)判断可用失效事件数据数量是否大于等于n,当判断结果为否时进入步骤b),否则根据双参数威布尔分布确定电站部件的失效率;
b)判断可用失效事件数据数量是否大于0,当判断结果为否时进入步骤c),否则利用bayes方法计算部件失效率;
c)判断是否存在通用失效率,当判断结果为否时进入步骤d),否则利用电站运行数据对通用失效率进行修正;
d)利用减函数方法计算部件失效率。
2.根据权利要求1所述的电站小样本及无失效数据部件失效率的可靠性评估方法,其特征在于:在步骤a)中当失效事件数据数量大于等于n时,则利用双参数威布尔分布模型计算部件失效率函数、故障分布函数、可靠度函数的点估计:
λ ( t ) = βt β - 1 α β , t ≥ 0 F ( t ) = 1 - exp [ - t β α ] , t ≥ 0 R ( t ) = exp [ - t β α ] , t ≥ 0 .
3.根据权利要求1所述的电站小样本及无失效数据部件失效率的可靠性评估方法,其特征在于:在步骤b)中当失效事件数据数量大于0时,
试验(r,T)似然函数:
λ取共轭型先验:f0(λ)=Γ(λ/r0,T0)
后验:f0(λ/r,T)=Γ(λ/r0+r,T0+T)
失效率单侧置信上限:
Figure FSA00000145400000015
4.根据权利要求1所述的电站小样本及无失效数据部件失效率的可靠性评估方法,其特征在于:在步骤c)中当存在通用失效率时,将通用失效率作为较保守的上限λ0,部件的寿命服从指数分布,有n个产品进行定时截尾试验,截尾时间段[0,ti],没有部件失效,则无失效时的修正似然函数为:
Figure FSA00000145400000016
C为常数
按Jeffreys准则,λ的先验分布取πλ∝λ-1,0<λ<λ0,在平方损失下的bayes估计为:
Figure FSA00000145400000017
5.根据权利要求1所述的电站小样本及无失效数据部件失效率的可靠性评估方法,其特征在于:在步骤d)中,采用无失效数据处理方法获得部件可靠度对于寿命服从指数分布的产品进行n次定时截尾试验,结果所有试验品无一失效,记试验总时间为T=ntx:
取π(λ)∝λα-1,α∈(0,1)为先验分布,则:
h ( λ / T ) = L ( t 1 , t 2 , · · · t r / λ ) · π ( λ ) ∫ 0 ∞ L ( t 1 , t 2 , · · · t r / λ ) · π ( λ ) dλ
= λ α - 1 · exp ( - n t 0 λ ) ∫ 0 ∞ λ a - 1 · exp ( - n t 0 λ ) dλ = ( n t 0 ) a Γ ( α ) λ α - 1 exp ( - n t 0 λ )
λ的条件期望值:
λ ^ = ∫ 0 ∞ λh ( λ / nt 0 ) dλ = Γ ( a + 1 ) T 0 = a nt 0 a = 1 2 χ α 2 ( 2 ) = - ln α
因此,取一定的置信度求得部件失效率:
λ ^ = - ln α T 0 .
6.根据权利要求1至5所述的电站小样本及无失效数据部件失效率的可靠性评估方法,其特征在于:n大于等于5。
CN201010192847.6A 2010-06-07 2010-06-07 电站小样本及无失效数据部件失效率的可靠性评估方法 Active CN101847834B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201010192847.6A CN101847834B (zh) 2010-06-07 2010-06-07 电站小样本及无失效数据部件失效率的可靠性评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010192847.6A CN101847834B (zh) 2010-06-07 2010-06-07 电站小样本及无失效数据部件失效率的可靠性评估方法

Publications (2)

Publication Number Publication Date
CN101847834A true CN101847834A (zh) 2010-09-29
CN101847834B CN101847834B (zh) 2016-02-03

Family

ID=42772338

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010192847.6A Active CN101847834B (zh) 2010-06-07 2010-06-07 电站小样本及无失效数据部件失效率的可靠性评估方法

Country Status (1)

Country Link
CN (1) CN101847834B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102207998A (zh) * 2011-06-09 2011-10-05 河南科技大学 基于无失效数据的评估产品寿命及其可靠性的方法
CN103761421A (zh) * 2013-12-31 2014-04-30 电子科技大学 一种大型矿用挖掘机提升机构可靠性评估方法
CN104899423A (zh) * 2015-05-06 2015-09-09 同济大学 一种动车组子系统关键部件运用可靠性评估方法
CN105023109A (zh) * 2015-07-30 2015-11-04 苏州热工研究院有限公司 一种核电站系统或设备可靠性评估方法和系统
CN105203152A (zh) * 2014-06-27 2015-12-30 国家电网公司 一种光伏发电设备故障风险指数预测装置及预测方法
CN106202776A (zh) * 2016-07-20 2016-12-07 上海宇航系统工程研究所 一种运载火箭结构机构产品可靠性评估方法
CN106251044A (zh) * 2016-07-21 2016-12-21 中国科学院数学与系统科学研究院 多批次成败型试验下产品贮存期评估的Buehler 方法
CN103995970B (zh) * 2014-05-26 2017-04-19 北京航空航天大学 一种离子推力器极小子样可靠性评估方法
CN107358342A (zh) * 2017-06-27 2017-11-17 苏州热工研究院有限公司 一种核电设备需求失效概率的评估方法及系统
CN109376407A (zh) * 2018-09-30 2019-02-22 中国人民解放军92942部队 一种按次使用武器装备的可靠度评估方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
周源泉等: "《可靠性评定》", 31 March 1990, 科学出版社 *
张宝珍等: "威布尔方法在机械产品可靠性设计、分析和评价中的应用(上)", 《航空标准化与质量》 *
王超等: "《机械可靠性工程》", 30 June 1992, 冶金工业出版社 *
韩明: "《基于无失效数据的可靠性参数估计》", 30 June 2005, 中国统计出版社 *
韩明: "《无失效数据的可靠性分析》", 30 April 1999, 中国统计出版社 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102207998B (zh) * 2011-06-09 2014-06-25 河南科技大学 基于无失效数据的评估产品寿命及其可靠性的方法
CN102207998A (zh) * 2011-06-09 2011-10-05 河南科技大学 基于无失效数据的评估产品寿命及其可靠性的方法
CN103761421A (zh) * 2013-12-31 2014-04-30 电子科技大学 一种大型矿用挖掘机提升机构可靠性评估方法
CN103761421B (zh) * 2013-12-31 2017-06-06 电子科技大学 一种大型矿用挖掘机提升机构可靠性评估方法
CN103995970B (zh) * 2014-05-26 2017-04-19 北京航空航天大学 一种离子推力器极小子样可靠性评估方法
CN105203152B (zh) * 2014-06-27 2018-06-19 国家电网公司 一种光伏发电设备故障风险指数预测装置及预测方法
CN105203152A (zh) * 2014-06-27 2015-12-30 国家电网公司 一种光伏发电设备故障风险指数预测装置及预测方法
CN104899423A (zh) * 2015-05-06 2015-09-09 同济大学 一种动车组子系统关键部件运用可靠性评估方法
CN104899423B (zh) * 2015-05-06 2017-12-15 同济大学 一种动车组子系统关键部件运用可靠性评估方法
CN105023109A (zh) * 2015-07-30 2015-11-04 苏州热工研究院有限公司 一种核电站系统或设备可靠性评估方法和系统
CN106202776A (zh) * 2016-07-20 2016-12-07 上海宇航系统工程研究所 一种运载火箭结构机构产品可靠性评估方法
CN106251044A (zh) * 2016-07-21 2016-12-21 中国科学院数学与系统科学研究院 多批次成败型试验下产品贮存期评估的Buehler 方法
CN107358342A (zh) * 2017-06-27 2017-11-17 苏州热工研究院有限公司 一种核电设备需求失效概率的评估方法及系统
CN109376407A (zh) * 2018-09-30 2019-02-22 中国人民解放军92942部队 一种按次使用武器装备的可靠度评估方法
CN109376407B (zh) * 2018-09-30 2023-06-20 中国人民解放军92942部队 一种舰载机弹射器的可靠度评估方法

Also Published As

Publication number Publication date
CN101847834B (zh) 2016-02-03

Similar Documents

Publication Publication Date Title
CN101847834A (zh) 电站小样本及无失效数据部件失效率的可靠性评估方法
Gao et al. Use of double mass curves in hydrologic benefit evaluations
CN102103658B (zh) 一种实用的电子产品寿命评估模型参数高精度萃取方法
Valor et al. Stochastic approach to pitting-corrosion-extreme modelling in low-carbon steel
Wilshire et al. A new approach to creep data assessment
Chaves et al. Extreme value analysis for assessing structural reliability of welded offshore steel structures
CN110374047A (zh) 基于变形的拱坝运行期实时安全监控阈值确定方法
Zhou Reliability evaluation of corroding pipelines considering multiple failure modes and time-dependent internal pressure
Harris XIMIS, a penultimate extreme value method suitable for all types of wind climate
CN102156820A (zh) 一种机械装配工艺的危害性分析方法
CN109145382A (zh) 一种载机驱动桥小子样可靠性的评估方法
CN103308334B (zh) 一种非线性累积的构件疲劳评估方法
CN105302946A (zh) 一种腐蚀管道可靠度确定方法及装置
Heng et al. GPR-assisted evaluation of probabilistic fatigue crack growth in rib-to-deck joints in orthotropic steel decks considering mixed failure models
Li et al. Reliability analysis of Lan Chengyu corroded pipeline with associated defects
Mahmoodian et al. Stochastic failure analysis of defected oil and gas pipelines
Wu et al. Time-dependent reliability model of components with strength degradation based-on gamma process
Keisman et al. Total maximum daily load criteria assessment using monitoring and modeling data
Prakash Probabilistic model for remaining fatigue life estimation of bridge components
CN110096820B (zh) 单元寿命服从威布尔分布时k/n(G)系统的剩余寿命预测方法
Zhang et al. Temporal change of spatial heterogeneity and its effect on regional trend of annual precipitation heterogeneity indices
Li et al. Extreme value theory-based structural health prognosis method using reduced sensor data
Espósito et al. Application of risk analysis methods on tailings dams
Shadabfar et al. Time-Variant System Reliability Analysis of Concrete Sewer Pipes under Corrosion Considering Multiple Failure Modes
CN113742814A (zh) 大坝安全预警方法、装置、计算机设备和存储介质

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent of invention or patent application
CB02 Change of applicant information

Address after: 215004 West Ring Road, Jiangsu, Suzhou, No. 1788

Applicant after: Suzhou Nuclear Power Research Institute Co., Ltd.

Applicant after: China General Nuclear Power Corporation

Address before: 215004 West Ring Road, Jiangsu, Suzhou, No. 1788

Applicant before: Suzhou Nuclear Power Research Institute Co., Ltd.

Applicant before: China Guangdong Nuclear Power Group Co., Ltd.

C14 Grant of patent or utility model
GR01 Patent grant