第9章可靠性数据检验与分布参数估计ppt课件.ppt

《第9章可靠性数据检验与分布参数估计ppt课件.ppt》由会员分享，可在线阅读，更多相关《第9章可靠性数据检验与分布参数估计ppt课件.ppt（50页珍藏版）》请在三一办公上搜索。

1、第9章数据检验与分布参数估计,9.1 可靠性试验与可靠性数据9.1.1可靠性试验可靠性分析需要可靠性数据。可靠性数据可以通过记录现场失效获得，也可以通过可靠性试验获得。可靠性试验是指为了获得可靠性数据、验证及提高产品的可靠性水平而进行的试验。,现场试验反映的是产品在使用条件下实际失效情况，最能说明产品的可靠性水平。实验室试验是在实验室内模拟现场情况进行的试验。在实验室中也可以进行影响寿命的单项因素试验，还可以进行加速试验。,第9章可靠性数据检验与分布参数估计,9.1.1可靠性试验无论通过那种途径获取可靠性数据，要等到全部观测对象失效，都需要很长的时间。根据统计学理论，在只有部分试样失效时，也能

2、获得有关的可靠性指标。因而，视可靠性试验的目的不同，有些试验可以在进行到一定时间即可停止。按试验截止情况，可分为定数截尾和定时截尾试验两种。定数截尾试验是试验到预定的失效样品数时停止试验。定时截尾试验是试验到规定的时间时停止试验。根据试验中试样失效后是否用新试样替换，可分为有替换试验和无替换试验。,截尾数据以可划分为：,型截尾（定时截尾）：在一项试验中，对多件样品同时进行试验，出现失效后也可以补充新样品，试验在一段预定时间后结束。例如，在试验台上同时测试100个轴承，并且在1000小时后结束测试，而不管发生失效的轴承数量有多少，这样的试验就是定时截尾试验。在这个例子当中，试验时间是固定的，而失

3、效数据数量为随机变量。II型截尾（定数截尾）：试验中最初有多个样品，试验到指定数量的样品失效时结束。例如，在在试验台上同时对100个轴承进行试验，并且在出现30个失效后终止试验。这样的试验就是定数截尾试验。在这个例子中，失效样品的数量是事先指定的，而试验时间是随机变量。,以上属于右截尾的情况，即如果一个样品失效的时间未知的原因是，该样品在试验结束时还没有失效，那么就产生了右侧被截尾的数据。如果一个样品已经失效，但不能确定准确失效时间的原因是不知道该样品是何时投入试验的，那么就产生了左截尾数据。单一截尾与多重截尾：单一截尾是只有一个截尾点。如果100个传感器置于一个检测台上，检测在1000小时后

4、完成，在1000小时处有一个截尾点。如果在测试1000小时后有20个没有失效，测试1200小时后还有15个未失，这时就有2个截尾点，结果数据是多重截尾的。如果确切的失效时间未知，但是在一个时间区间内失效的次数被记录到了，那么这就是区间数据或成组数据。,四种组合形式：,有替换定时截尾寿命试验；有替换定数截尾寿命试验；无替换定时截尾寿命试验；无替换定数截尾寿命试验。,9.1.2 可靠性试验数据,可靠性数据多种多样，很重要的一类是截尾数据。通过试验或现场记录获得的数据可能既包括失效数据，也包括不失效数据。以下表种的数据为例，在试验台上进行八个产品的试验；其中的三个失效了，五个到试验结束时还没有失效。

5、样品序号失效时间1 35245355460560660760860,很明显，只用三个失效样本平均值及标准偏差不能用来估算这个案例中的分布参数。这三个样本平均值为（35+45+55）/3=45。存活下来的五个数据中每个对应的失效时间都大于60，因此真实的样本均值要比45大的多。因此，需要处理截尾试验数据的方法。,9.1.3截尾数据统计推断,假设观测母体包括n个样本，在t0=0时投入运行。运行过程中，有些样本发生了失效，也会有样本退出运行。记录的数据包括每个失效产品的失效时间和各失效发生前仍在运行的样本数量。令发生失效的时间为t1,t2,tk按从不到大的顺序排列，ni为时刻ti之前仍在运行的样本数

6、，wi为在(ti-1)ti时间段内退出运行（并未发生失效）的样本数。显然，n1=n-w1,n2=n1-1-w2，以此类推。,要根据这样的观测数据估计样本的可靠度，有如下的Kaplan-Meier公式：(9-1),9.1.4参数估计1.矩估计,一阶原点矩等于分布的平均值：(9-2)对均值的二阶矩等于分布的方差：(9-3)其中是分布均值。对于正态分布等，这些矩提供了参数的直接估计量。对于威布尔、对数正态和伽玛分布等，分布参数需要通过令样本矩等于理论矩来解出分布参数。所需矩的数目取决于被估计的参数的数目。,2.极大似然估计,极大似然估计法是一种应用广泛的估计参数的方法。其基本原理是，通过使已知样本出

7、现的概率最大化的方法来确定参数值。令 是服从概率密度函数 的独立随机变量，其中 是唯一的分布参数。那么(9-4)就是随机变量的联合分布或似然函数。分布参数的极大似然估计值，是使似然函数取最大值的。,为了便于计算，通常先将似然函数取自然对数，即(9-5)解方程(9-6)即可求得分布参数 的估计值。,对于概率密度函数为f(x,)，累积分布函数为F(x,)的随机变量，如果只知道其样本中的前k个顺序统计量x(1),x(2),x(k)，而关于其它n-k个样本的已知信息只是大于x(k)，这种情况下的似然函数形式为(9-7),如果只知道第i个样本在运行时间段0-xi（xi为已知量）内是否失效，不知道准确的失

8、效时间，则有如下似然函数(9-8)式中，k为失效样本数，n-k为未失效样本数。,9.1.5 分布类型假设检验,分布类型的判断有理论法和统计法两种。理论法是根据失效机理制定的数学模型或根据某种分布的性质推导出来的。例如，失效率为常数的寿命分布为指数分布；失效由“最弱”环节决定的寿命分布为威布尔分析或极值分布；受很多独立随机因素的影响，且没有一个因素起主导作用，这种分布为正态分布等。统计法是根据大量试验数据经统计得出的。很多同类性能在以往大量试验的基础上已经验证了其分布。例如，几何尺寸、材料性能、硬度等多服从正态分布；零件疲劳寿命则服从对数正态分布或威布尔分布等。在使用统计法时，对分布未知的情况下

9、应做大样本的试验，以判定其分布类型；对已有经验参考的情况则可做较小样本的试验，假设其分布类型再进行相应的拟合性检验。下面介绍通用的2检验法和K-S检验法。,9.2 检验法,检验法适用于大样本试验数据。其基本思想是，将随机试验的全部可能结果划分为 k 个互不相容的事件 A1，A2，Ak，在假设成立的条件下计算 P(Ai)=pi，i=1,2,k。在 n 次试验中，事件 Ai 出现的频率 ni/n 与 pi 常有差异。但由大数定律可知，如果试验次数很多，在假设成立的条件下，ni/n-pi的值应该很小。由此，检验法首先计算理论频数与实际频数间的差异，将统计量的观测值与临界值比较。满足条件则接受原假设；

10、否则拒绝原假设。公式如下：（9-9）,式中:n-样本数；k-分组数，按样本大小取；-第组的实际频数，；-第组的理论频数（概率）；-未知参数的数目；-显著性水平；-临界值，可由分布表查出。,例9-1 某220件产品的失效时间记录如下所示，试检验该产品的寿命是否服从指数分布。,产品失效时间的数据记录时间t：0-100-200-300-400-500-600-700-800-900失效数i：39 5035322818124 2解 假设该产品的寿命服从指数分布，参数未知。取每组的中值作为该组时间的表征值 ti,平均寿命的点估计为失效率的点估计为 假设H0：,为了使用检验法，首先对数据进行分组。由于每组

11、中实际频数不宜少于5个，故将前7段时间各作为一组，最后两段时间合为一组。总计组数k=8，在7-14范围内。具体计算列表如下。例9-1计算列表组号 ii npi=220pi i npi 1 39 0.2827 62.194-23.194 537.962 8.650 2 50 0.2055 45.210 4.790 22.944 0.507 3 35 0.1461 32.140 2.860 8.180 0.254 4 32 5 28 6 18 7 12 8 6,统计量为取显著性水平，由，查表，得由于，故拒绝原假设，既不能认为该产品的寿命服从指数分布。,9.3 K-S检验法,K-S 检验法（亦称 D

12、 检验法）也是比较样本经验分布函数Fn(x)与总体分布函数 F0(x)。但不是在划分的区间上考虑样本经验分布函数 Fn(x)与假设的总体分布函数F0(x)之间的偏差，而是在每一点上考虑它们之间的偏差。K-S检验法比 2 检验法精确。但要求所检验的分布中不含未知参数，且要假定总体分布函数为连续函数。具体做法是，先将n个试验数据由小到大的次序排列。根据假设的分布，计算每个数据对应的分布函数 F0(xi)，并将其与经验分布函数 Fn(xi）比较。其中，差的最大绝对值就是检验统计量 Dn 的观测值。将 Dn 与临界值 Dn,比较。满足条件，接受原假设；否则拒绝原假设。条件如下：,（9-10）,式中 F

13、0(x)原假设的分布函数；Fn(x)经验分布函数。具体计算公式如下（9-11）（9-12）Dn,a临界值，由K-S检验临界值表查得。,例9-2 测得某合金9个试件的强度极限数据453，436，429，419，405，416，432，423，440 MPa。检验该合金的强度极限是否服从均值=428 MPa，标准差=15 MPa的正态分布。解 用X 表示该合金的强度极限，将试验数据按由小到大次序排列。假设的 X 分布函数为式中的可由正态分布表查得。具体计算过程如下：,序号i xi(i-1)/n i/n di14050.0630 0.000 0.111 0.04812 416 0.2119 0.11

14、1 0.222 0.10093 419 0.2743 0.222 0.333 0.0587 由上述中计算结果知，Dn的观测值按式(9-2)取显著性水平=0.10，由表K-S检验表查得。由于，故接受原假设，即认为该合金的强度极限服从=428MPa，=15 MPa的正态分布。,9.4 回归分析检验法,对于给定的n个数据对（xi,yi），按最小二乘原理拟合出一条直线（9-13）称为回归直线。斜率 称为回归系数，截距 为常数项。计算公式如下：(9-14)y与x是否具有线性相关的关系，可用相关系数检验法进行检验。相关系数定义为(9-15),若，则认为y与x具有线性相关的关系。是显著性水平为 时的相关系数

15、临界值，由相关系数表查得（自由度）。表-相关系数临界值表 0.10 0.05 0.02 0.01 0.001 5 0.6694 0.7545 0.8329 0.8745 0.9507,对于可靠性工程中常用的概率分布，其分布函数与自变量之间一般在直角坐标上并不成线性关系。为了应用线性回归检验方法，需要先进行适当的变换，使变换量与自变量之间呈线性关系（常用概率分布的变换关系列于表中-）。表-几种常用概率分布的变换关系名 称 分布函数 y x B A指数分布 t 0,在使用回归分析法时，首先将试验获得的 n 个数据按由小到大的次序排列，即t1t2tn，取中位秩作为各数据点对应的分布函数，即（9-16

16、）假设一种分布，变换后即可用上述各公式进行计算。若相关系数检验通过则接受原假设。估计得B、A后再估计原分布函数的参数。,例9-3 在某应力水平下，某合金材料10个试件的疲劳寿命分别为211，229，272，276，295，303，332，354，382，409千次。试进行分布类型的判断并进行参数估计。解 一般金属疲劳寿命服从对数正态分布，故先假设该材料疲劳寿命服从对数正态分布。进行变换，计算结果如下：序号Ni12115.3520.067-1.50028.6422.250-8.028,由式（9-15），算出相关系数为 由表查得，当，。由于，故接受疲劳寿命服从对数正态分布的假设。,9.5 指数分布

17、假设检验与参数估计,指数分布的概率密度函数为（9-17）,9.5.1 拟合性检验,检验统计量为（9-18）式中-总累积试验时间；Tk-第（=1，2，）次失效时的累积试验时间；t0指定的定时截尾时间；tr指定的定数截尾时间。,上述检验统计量服从自由度为2d的2分布。给定显著度水平，接受指数分布的假设的条件如下：（9-19）式中，表示自由度为d 的 分布的分位数，查表确定。总累积的试验时间 是所有投入试验的试样（包括失效的、中止的、截尾未失效的）试验到规定时间的试验时间总和。当开始投入 n 个试样同时试验，试验中有 b 个中止，中止时间为（j=1，2，）；有 r 个失效，失效时间为i（i=1，2，

18、r）。规定试验到 t0 停止试验，则试验总累积时间计算方法如下。,无替换试验：（9-20）有替换试验：（9-21）t0定时截尾规定的截尾时间或定数截尾试验中第 r 个失效时间 tr；b中途中止试验的试样个数。,第 k 次失效时的累积试验时间 Tk的计算公式为无替换试验：（9-22）有替换试验：（9-23）式中 tk第 k（k=1，2，r）个试样失效时间；bk第 k 个试样失效前中止试验的试样个数。,例9-4 取10个产品样本进行寿命试验，失效5个时停止试验。这个样本的寿命分别为：76，143，152，275，326 h。检验该产品寿命是否服从指数分布。解 假设该产品的寿命服从指数分布。这是无替

19、换定数截尾寿命试验。由式（9-19）计算总累积试验时间由式（9-17）计算第k次失效时的累积试验时间,本例，由式（9-13），取显著性水平=0.10，由表查得，可见，检验统计量满足条件，故接受原假设，即认为该产品的寿命服从指数分布。,9.5.2 参数估计,指数分布进行参数估计的最简单的方法是极大似然法。极大似然法提供一个无偏估计，但不表明相符合的程度。指数分布概率密度函数是：（9-24）对于参数的极大似然估计是：（9-25）其中，xi第i个数据，它可能是失效数据，或截尾数据，n包括截尾和未截尾的数据总数，r失效数据个数。此估计是无偏估计和最小方差估计。,例 9-5 若7个弹簧失效时经历的载荷循

20、环数分别为30183，14871，35031，76321，43891，31650和12310。假设寿命循环数服从指数分布，估计平均寿命和平均失效率。解：平均寿命：失效率：,例 9-6 如果上例中的数据代表7根弹簧的失效循环数，但另有10根弹簧测试了80000个载荷循环而未失效的。估计平均失效时间和平均失效率。解：平均失效时间：次循环平均失效率是：,9.6 正态分布统计检验,正态分布的概率密度为（9-26）9.6.1 拟合性检验样本容量不大、分布参数、未知时，用两参数的点估计时，假设（9-27）,与K-S检验法类似，满足下列条件则接受原假设：（9-28）式中 经验分布函数；（9-29）临界值，查表。,例9-10 对某钢材进行静强度试验，9个试件的强度极限按由小到大次序分别为625，650，656，659，661，662，663，668，672 N/mm2。检验该钢材强度极限是否服从正态分布。解：假设该钢材的强度极限服从正态分布。由于分布参数未知，故先进行估计：,假设：计算列表如下,由上表中的计算结果，取显著性水平=0.10，由表12-5查得：接受原假设，即认为该钢材的强度极限服从正态分布。,END,