可靠性试验及数据处理 (3).ppt

可靠性试验及数据处理方法,东北大学孙志礼,可靠性试验及数据处理方法,可靠性试验及分类分布类型的假设检验指数分布的分析法正态及对数正态分布的分析法威布尔分布的分析法,一、可靠性试验及分类,可靠性试验就是为了提高和证实产品的可靠性水平而进行的各种试验的总称,可靠性试验是为了获得统计数据，所用时间较长、所花的费用较大。但从提高和保证产品的质量角度来讲是值得的，费效比是较高的。寿命试验是可靠性试验的重要组成部分，是评价、分析产品寿命可靠性特征量所进行的试验。下面列出几种寿命试验的分类,按试验场所分：现场试验和实验室试验两种。现场试验是产品在使用条件下观测到的寿命数据。最能说明产品的可靠性水平，是最终的客观标准。，收集现场数据重要。但会遇到很多困难，需要的时间较长、工作情况难以一致，要有详细的产品使用记录，很难获得比较准确的数据实验室试验是模拟现场情况的试验。它将现场重要的应力条件搬到实验室，并加以人工控制。还可设法缩短试验时间以加速取得试验结果,一、可靠性试验及分类,按试验截止情况分：分为全数试验和截尾试验两种全数试验是当试样全部失效才停止的试验这种试验方式可获得较完整的试验数据，统计分析结果也较好。但这种试验所需时间较长，有时甚至难以实现,一、可靠性试验及分类,截尾试验又可分为定数和定时截尾试验两种定数截尾试验就是试验到规定的失效数即停止的试验定时截尾试验就是试验到规定的时间，此时不管试样失效多少都停止的试验根据试验中试样失效后是否用新试样替换继续试验，还可分为有替换和无替换两种,一、可靠性试验及分类,一般可归纳为如下四种试验：有替换定时截尾寿命试验；有替换定数截尾寿命试验；无替换定时截尾寿命试验；无替换定数截尾寿命试验。全数寿命试验也可看成是截尾数是n的无替换定数截尾寿命试验。此外，尚有分组最小寿命试验、序贯寿命试验、有中止的寿命试验等,一、可靠性试验及分类,二、分布类型的假设检验,分布类型的判断有理论法和统计法两种理论法是根据失效机理制定的数学模型或根据某种分布的性质推导出来的例如，失效率为常数的寿命分布为指数分布；失效由“最弱”环节决定的寿命分布为极值分布；受很多独立随机因素和的影响，且没有一个因素起主导作用，这种分布为正态分布等。统计法是根据大量试验数据经统计求得的。很多同类性能在以往大量试验的基础上已经验证了其分布例如，几何尺寸、材料性能、硬度等多服从正态分布；金属的疲劳寿命则服从对数正态分布或威布尔分布等,下面仅介绍统计法在使用统计法时：对分布不明的情况应做大样本的试验以判定其分布类型；对已有经验参考的情况则可做较小样本的试验，假设其分布类型再进行相应的拟合性检验,下面给出通用的,检验法和K-S检验法。,二、分布类型的假设检验,检验法一般只用于大样本,计算理论频数与实际频数间的差异，将检验统计量,的观测值与临界值,满足下列条件，接受原假设；否则，拒绝原假设,比较,二、分布类型的假设检验,式中n样本大小k分组数，按样本大小宜取第组的实际频数，第组的理论频数（概率）未知参数的数目显著性水平,临界值，查表,二、分布类型的假设检验,例1：220个某产品的失效时间记录列于表中。试检验该产品的寿命是否服从指数分布。某产品失效时间的数据记录,二、分布类型的假设检验,解假设该产品的寿命服从指数分布，参数未知。取组中值作为该组时间的代表值ti,则的点估计,h,1/h,二、分布类型的假设检验,假设H0：,为了使用,首先按规定分组。由于每组中实际频数不宜少于5，故将前7段时间各作为一组，最后两段时间合为一组。总计组数k=8，正好在714范围内,检验法,二、分布类型的假设检验,表2-2例2-1的列表计算,二、分布类型的假设检验,二、分布类型的假设检验,取显著性水平,，由,查表,由于,故拒绝原假设，既不能认为该产品的寿命服从指数分布。,K-S检验法（亦称d检验法）适用于小样本的情况,K-S检验法要求所检验的分布中不含未知参数。当指定分布中含有未知参数时，对某些分布应该用专门的临界值表,二、分布类型的假设检验,K-S检验法是将n个试验数据由小到大的次序排列。根据假设的分布，计算每个数据对应的F0(xi)，将其与经验分布函数Fn(xi）相比较。其中，差的最大绝对值就是检验统计量Dn的观测值。将Dn与临界值Dn,比较。满足下列条件，接受原假设；否则，拒绝原假设,F0（x）原假设的分布函数；Fn（x）经验分布函数,二、分布类型的假设检验,临界值，查表,二、分布类型的假设检验,例2某合金9个试件测得的强度极限为453，436，429，419，405，416，432，423，440N/mm2。检验该合金的强度极限是否服从均值=28N/mm2，标准差=15N/mm2的正态分布。,二、分布类型的假设检验,解令该合金的强度极限b=X，将数据按由小到大次序排列。假设X服从正态分布，分布函数,式中的,查表。计算结果见表,由上述中计算结果知，Dn的观测值按前式求,二、分布类型的假设检验,取显著性水平=0.10，由表2-3查得,由于,故接受原假设，即认为该合金的强度极限服从=28N/mm2，=15N/mm2的正态分布,二、分布类型的假设检验,二、分布类型的假设检验,K-S检验临界值表,n,回归分析法：回归分析法就是图解分析法的解析。在直角坐标纸上描得几个试验点（x1，y1），（x2，y2），（n，yn），如图2-1所示。按最小二乘原理确定直线,二、分布类型的假设检验,它反映出试验点散布状态的一条最佳直线，称为回归直线。斜率,称为回归系数，截距,为常数项。,二、分布类型的假设检验,各试验点是否在一直线上，即是否具有线性相关的关系，可用相关系数检验法进行检验。相关系数,二、分布类型的假设检验,当,，则认为具有线性相关的关系。,是显著性水平为,时的相关系数起码值，查表。查表时取自由度,对于可靠性分析中常用的概率分布，其分布函数与自变量之间一般在直角坐标上并不成线性关系，因此应先进行适当的变换。几种常用概率分布的变换关系列于下表中,二、分布类型的假设检验,相关系数起码值,分布函数,二、分布类型的假设检验,在使用回归分析法时，首先将试验获得的n个数据按由小到达的次序排列，即t1t2tn，取中位秩作为各试验点相应的分布函数，即,假设一种分布，按上表进行变换后即可用前式进行计算。若相关系数检验通过则接受原假设。估计得B、A后再按上表关系估计原分布函数的参数,二、分布类型的假设检验,例3某合金材料在某应力水平做疲劳寿命试验，10个试件的疲劳寿命分别为211，229，272，276，295，303，332，354，382，409千次。试进行分布类型的判断并进行参数估计。解一般金属疲劳寿命多较好地服从对数正态分布，故先假设该材料疲劳寿命服从对数正态分布。按表,二、分布类型的假设检验,二、分布类型的假设检验,计算结果见表,相关系数检验：,二、分布类型的假设检验,由表，当,查得,，故接受疲劳寿命服从对数正态分布的假设,二、分布类型的假设检验,估计分布参数和，由表2-6知：,二、分布类型的假设检验,三、指数分布的分析法,拟合性检验,这种检验法适用于截尾试验、全数试验和有中止的试验,计算检验统计量,满足下列条件则接受指数分布的假设，否则拒绝指数分布的假设,三、指数分布的分析法,总累积试验时间；,第（=1，2，）次失效时的累积试验时间,t0指定的定时截尾时间tr指定的定数截尾时间,显著性水平；,自由度为d的,分布的分位数，查表,三、指数分布的分析法,总累积的试验时间,是所有投入试验的试样（包括失效的、中止的、截尾未失效的）试验到规定时间的试验时间总和。当开始投入n个试样同时试验，试验中有b个中止，中止时间为,（j=1，2，）；有r个失效，失效时间为i（i=1，2，r）。规定试验到t0停止试验，则试验总累积时间,三、指数分布的分析法,无替换：,有替换：,第k次失效时的累积试验时间,无替换：,有替换：,第k（k=1，2，r）个失效的时间,三、指数分布的分析法,例4抽取某产品10个进行寿命试验，失效5个即停止试验。试验结果为：76，143，152，275，326h。检验该产品寿命是否服从指数分布。,解假设该产品的寿命服从指数分布。这是无替换定数截尾、无中止的寿命试验。总累积试验时间,三、指数分布的分析法,第k次失效时的累积试验时间,三、指数分布的分析法,本例,取显著性水平=0.10，由表查得,满足,故接受原假设，即认为该产品的寿命服从指数分布。,三、指数分布的分析法,参数估计和可靠度估计,全数试验,截尾试验,式中,总累积试验时间r观测失效数，当r=0进行的点估计时，建议取,。,三、指数分布的分析法,区间估计,单侧置信下限,双侧置信下限,双侧置信上限,三、指数分布的分析法,参数的区间估计,先由上表计算平均寿命的置信限,的置信下限,的置信上限,可靠度的点估计,可靠度的置信下限,三、指数分布的分析法,例5某产品寿命服从指数分布，抽取11个进行寿命试验。在试验到500h式中止1个，900h时失效1个，其它试样达到1000h均未失效即停止试验。求平均寿命、失效率及工作到100h时可靠度的点估计。若要求置信水平,求平均寿命的单侧置信下限、失效率的单侧置信上限及工作到100h时可靠度的单测置信下限。,三、指数分布的分析法,解这是n=11，失效数r=1，中止数b=2，截尾时间t0=1000h的无替换定时截尾寿命试验,总累积试验时间,平均寿命的点估计,失效率的点估计,三、指数分布的分析法,t=100h时可靠度的点估计,由表查得,故平均寿命的单侧置信下限,失效率的单侧置信上限,t=100h时可靠度的单侧置信下限,四、正态及对数正态分布的分析法,拟合性检验,假设,与K-S检验法类似，满足下列条件则接受原假设，否则拒绝原假设,经验分布函数,四、正态及对数正态分布的分析法,临界值，查表,四、正态及对数正态分布的分析法,例6对某钢材进行静强度试验，9个试件的强度极限按由小到大次序分别为625，650，656，659，661，662，663，668，672N/mm2。检验该钢材强度极限是否服从正态分布。,四、正态及对数正态分布的分析法,解假设该钢材的强度极限服从正态分布。由于分布参数未知，故先进行估计,=13.69,四、正态及对数正态分布的分析法,假设：,列表计算,由计算结果知,取显著性水平=0.10,因为,故接受原假设，即认为该钢材的强度极限服从正态分布,四、正态及对数正态分布的分析法,完全样本的参数估计,点估计,区间估计,显著性水平，1为置信水平；n样本大小；自由度。当标准差已知，sx用代替，则=；当标准差为点估计sx，则=n1；,四、正态及对数正态分布的分析法,自由度。当均值已知时，=n；当未知时，=n1。,若随机变量y服从对数正态分布，即YLn（，2），则lnYN（，2），故取,四、正态及对数正态分布的分析法,截尾寿命试验的参数估计,极大似然估计,失效时间,四、正态及对数正态分布的分析法,式中,查表。,四、正态及对数正态分布的分析法,最佳线性无偏估计,式中n样本大小；r失效数；j寿命由小到大排列的次序；xj第j个寿命值；,的最佳线性无偏估计系数,的最佳线性无偏估计系数,四、正态及对数正态分布的分析法,简单线性无偏估计,式中nkr,n系数查表,E（Yr,n）系数，查表,其余同前,四、正态及对数正态分布的分析法,可靠寿命和可靠度的估计,指定可靠度R，可靠寿命的点估计,指定寿命x的可靠度R的点估计,四、正态及对数正态分布的分析法,指定可靠度RL的可靠寿命置信下限,zR单侧置信限系数，按指定的可靠度和置信水平由下式求取,z按指定的置信水平查表,四、正态及对数正态分布的分析法,指定寿命xL的可靠度置信下限,式中,应该指出，这里的x并不限于寿命，也可以是服从正态分布的其它特性值，例如材料的机械强度等,四、正态及对数正态分布的分析法,极大似然估计用表（正态及对数正态分布）,最佳线性无偏估计用表（正态及对数正态分布）,四、正态及对数正态分布的分析法,简单线性无偏估计用表（正态及对数正态分布）,四、正态及对数正态分布的分析法,zR与可靠度R的关系（正态分布）,不同置信水平时的z,四、正态及对数正态分布的分析法,例某钢材的强度极限服从正态分布，11个试件测得的强度极限为608，622，630，638，642，648，652，660，673，688N/mm2。求均值和标准差的点估计和置信水平=80的双侧置信限，失效概率为0.10时强度极限的点估计和置信水平为90的单侧置信下限。,四、正态及对数正态分布的分析法,解本例为完全样本试验。均值的点估计，按式,标准差的点估计，按式,四、正态及对数正态分布的分析法,均值的双侧置信限，按式并由表查得,故得,四、正态及对数正态分布的分析法,标准差的双侧置信限，按式并由表查得,故得,四、正态及对数正态分布的分析法,失效概率F=0.10的强度极限，可借用可靠寿命与可靠度的关系式。这时相当于R=1-F=0.90时的强度极限。点估计可用式来求，并由R=0.90查表2-15得zR=1.28，故,四、正态及对数正态分布的分析法,失效概率F=0.10，置信水平=90强度极限的单侧置信下限可用式来求，并由R=1-F=0.90查表得zR=2.01129，故,五、威布尔分布的分析法,拟合性检验,样本大小为n，截尾寿命试验得t1t2tr。检验统计量观测值,式中,即取括号内整数部分,E(zi)查表中的E(zi，n),五、威布尔分布的分析法,满足下式条件则接受两参数威布尔分布的假设，否则拒绝两参数威布尔分布的假设,式中显著性水平,自由度为1=2(r-r0-1)，2=2r0的F分布的分位数，查表,五、威布尔分布的分析法,参数估计,矩法估计,样本均值,样本标准差,样本偏态系数,五、威布尔分布的分析法,根据kk值由表查得形状参数k的点估计,同时可查得系数ka、kb，故,尺度参数b的点估计,位置参数a的点估计,五、威布尔分布的分析法,若位置参数a已知（例如两参数威布尔分布a=0），则求,根据kc值由表查得形状参数k的点估计,同时查得kb，再由式求尺度参数b的点估计,五、威布尔分布的分析法,极大似然估计,五、威布尔分布的分析法,最佳线性无偏估计和简单线性无偏估计,先做如下变换：,原分布函数,令：,即,即,即,五、威布尔分布的分析法,则原分布函数变为,式中，参数和的最佳线性无偏估计,五、威布尔分布的分析法,当样本容量较大时，和用简单线性无偏估计。和的简单线性无偏估计,s，nkr,n，E(zr,n)查表,五、威布尔分布的分析法,形状参数k的点估计按式求得后，为得到k的无偏估计再加以修正得,式中gr,n修偏系数，查表,尺度参数的点估计,五、威布尔分布的分析法,可靠度和可靠寿命的估计,指定寿命t时可靠度R的点估计,指定可靠度R时可靠