可靠性数据分析的计算方法

可靠性数据分析的计算方法

PROCEEDINGS,Annual RELIABILITY and MAINTAINABILITY Symposium（1996）

可靠性数据分析的计算方法

Gordon Johnston, SAS Institute Inc., Cary 关键词：寿命数据分析加速试验修复数据分析软件工具

摘要&结论

许多从事组件和系统可靠度研究的专业人员并没有意识到，通过廉价的台式电脑的普及使用，很多用于可靠度分析的功能强大的统计工具已经用于实践中。软件的计算功能还可以将复杂的计算统计和图形技术应用于可靠度分析问题。这大大的便利了工业统计学家和可靠性工程师，他们可以将这些灵活精确的方法应用于在可靠度分析时所遇到的许多不同类型的数据。

在本文中，我们在SAS@系统中将一些最有用的统计数据和图形技术应用到例子的当中，这些例子主要包涵了寿命数据，加速试验数据，以及可修复系统中的数据。随着越来越多的人意识到创新性软件在可靠性数据分析中解决问题的需要，毫无疑问，计算密集型技术在可靠性数据分析中的应用的趋势将会继续扩大。

1.介绍

本文探讨了人们在可靠性数据分析普遍遇到的三个方面：

寿命数据分析

试验加速数据分析

可修复系统数据的分析

在上述各领域，图形和分析的统计方法已被开发用于探索性数据分析，可靠性预测，并用于比较不同的设计系统，供应商等的可靠性性能。

为了体现将现代统计方法用于结合使用高分辨率图形的使用价值，在下面的章节中图形和统计方法将被应用于含有上述三个方面的可靠性数据的例子中。2.寿命数据分析

概率统计图的寿命数据分析中使用的最常见的图形工具之一。Weibull 图是最常见的使用可靠性的概率图的类型，但是当Weibull概率分布并不符合实际数据的时候，类似于对数正态分布和指数分布这一类的概率图在寿命数据分析中也能够起到帮助。

在许多情况下，可用的数据不仅包含故障时间，但也包含在分析时没有发生故障的单位的运行时间。在某些情况下，只能够知道两次故障发生之间的时间间隔。例如，在测试大量的电子元件时，如果记录每一个发生故障的元件的故障时间，那么这可能不经济。相反，在固定的时间间隔内

检查这些电子元件，并在每个时间间隔内记录下故障的数量。概率图，可以构建右删失，间隔删失这两种不完整数据的类型，在这两种类型中所有的元件都有相同的检测时间间隔。

通常，这些概率图都是通过计算绘制位置来构造，并且将它们绘制在专业的可靠度论文上。这是一项繁琐而耗时的过程。Hazard图含有类似的信息和更容易绘制，但手工绘制这些图仍然是耗时的。利用计算机程序绘制概率图可以更快速，更精准。而且它也更容易与不同的概率分布和模型，从而提高结果的准确性。

Nelson在下文给出了柴油发动机风扇装置出现故障的例子。在观测的70个风扇单元中，在进行数据分析时有12个发生了故障。其他的能够正常工作，分析的目的是为了预测在8000个小时的保修期内故障的发生率。图1显示了风扇数据的Weibull概率分布图。

上图所标示的点代表一个风扇的寿命数据的概率分布函数的估计值。由于图像是呈Weibull分布，所以当所标示的点大致呈直线分布的时候，那么Weibull概率分布就能与观测数据大致吻合。在这种情况下，Weibull 分布比较适用。顶端的点代表没有发生故障的风扇的检查时间和运行时间

在图中的位置。

图1发动机风扇数据的Weibull概率分布图

确定了Weibull概率分布与数据完全拟合之后，下一步就是进行分布参数的估计。许多统计学家和工程师更倾向于用最大似然法来估计参数。

最大似然法涉及数似然函数的最大化。

这是一个计算密集型程序，并且在没有计算机程序时难以执行，除非在某些特殊情况下。这个软件还可进行分布参数、可靠度函数以及分布百分比的最大似然估计。

最大似然方法有以下几个优点：

1.它提供了具有良好的统计特性的客观估计值（参考文献1）

2.它提供常见性的错误，这有利于评估参数估计和相关的概率分布百分比和可靠性函数估计的精确性。

3.它适用于任何数据类型的组合。

4.它适用于任何概率分布

对于发动机风扇的数据，最大似然Weibull的规模和形状参数估计如图1的插图所示。实直线代表拟合的Weibull分布，实曲线代表拟合的Weibull分布的95%置信区间。由于分析的目的是预测未来在8000小时的保修期内会失效的部件的概率，所以在8000小时的时间范围内绘制了一条参考线。这使得预测在8000小时内发生失效的概率更容易在Y轴上读取约25%。相同的信息可以用表格的形式表达，但是通过拟合线将用于报告或演示的所需的信息和结果表示在图上是显示数据的一种简便方法。

3.加速寿命试验模型

产品在设计的温度，荷载，电压，或其他应力的水平上操作时，可能有很长的寿命。加速试验主要是对产品进行过度加载，让产品快速产生故障，然后收集可靠度信息。统计模型已应用于许多类型的产品加速寿命试验分析，并用来显示产品可靠度和应力之间的关系。最大似然法通常用于拟合结合了加速试验的统计模型因为这类方法可用于多种模型和数据类型。除了分析方法，利用图形的方法也可用于显示数据，评估模型拟合程度，显示拟合的模型。这些都是计算密集型的程序，并且需要专业软件来完成它们。

在下面改编自Nelson的例子中，在电力变压器绝缘液体中做加速寿命试验来显示电压和时间之间的关系。在图2中并列的显示了两张图，用于显示加速寿命试验的数据和分析结果。

图2.绝缘液体数据的Weibull分布图和关系图

数据的Weibull概率图绘制在左侧，不同的电压等级对应不同的标准。每条对电压的拟合线代表拟合的统计模型。对于所有电压等级，模型都有相等的Weibull形状参数，在这种情况下，它是拟合的。这表现为在Weibull图上的平行拟合线。相同形状的模型的拟合度可以直观的从图2的绘图中进行评估。尽管没有显示出来，但结合利用计算机程序的最大似然拟合可以让你轻松的完成平行线假设的统计检验。

关于数据和拟合关系的曲线图绘制在图2右侧。在这种情况下，寿命与电源能量有关。当在日志记录表上绘制时，关系图是一条直线。例如，在中间的直线代表拟合关系的中位数（第五十百分位）。你可以利用这张图去推断在任何所需压力值下的寿命中点。为了保证关系图的完整性，还绘制了第10位和第90位的百分比。

其他的相关曲线可以很容易的绘制出来，例如Arrhenius这类寿命与绝对温度呈倒数的关系。像温度和电压这类模型，因为它们的寿命不仅仅只与某一个应力向量有关，所以这类模型也比较容易绘制。同时，上述这些模型也可以很容易与最大似然法拟合。

4.可修复系统

可修复系统就是在故障发生的状态下能还原和修复的系统。在可修复系统中，发生故障之间的时间不一定都具有相同的统计特性。例如，平均故障间隔时间可能会随时间变化。作为一个可修系统应用广泛的时代，它累计需要维修的次数和修理费用。修复数据分析不同于平常的数据分析，在修复数据分析中元件只能发生一次