汽车可靠性理论基础课件

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,汽车可靠性理论基础,汽车可靠性理论基础,内 容 提 纲,概率统计基础知识,可靠性函数,可靠性理论分布,可靠性数据分布,3/106,内 容 提 纲概率统计基础知识3/106,一、概率统计基础知识,一、概率统计基础知识,总体,：,研究对象的全体。,通常指研究对象的某项数量指标。,组成总体的元素称为,个体。,从本质上讲，总体就是所研究的随机变量或随机变量的分布。,1,、总体与样本,5/106,总体：研究对象的全体。从本质上讲，总体就是所研究的随机变量,样本：,来自总体的部分个体,X,1,，,，,X,n,如果满足：,（,1,）,代表性：,X,i,，,i=1,n,与总体同分布,.,（,2,）,独立性：,X,1,，,，,X,n,相互独立；,则称为,容量为,n,的简单随机样本,，简称,样本,。,而称对样本进行一次观察，得到一组 确定的值,x,1,，,，,x,n,为,样本观察值,。,6/106,样本：来自总体的部分个体X1， ， Xn （1）代表性：,总体、样本、样本观察值的关系,总体,样本,样本观察值,？,理论分布,统计是从手中已有的资料,样本观察值，去推断总体的情况,总体分布。样本是联系两者的桥梁。总体分布决定了样本取值的概率规律，也就是样本取到样本观察值的规律，因而可以用样本观察值去推断总体,7/106,总体、样本、样本观察值的关系总体 样本 样本观察值 ？理论分,2,、随机事件及概率,随机现象,在一定条件下，并不总出现相同结果的现象,随机事件,在一定条件下，可能发生也可能不发生的试验结果,随机事件间相互关系,事件包含关系,和事件示意图,8/106,2、随机事件及概率随机现象事件包含关系和事件示意图8/106,积事件示意图,差事件示意图,互不相容事件示意图,对立事件示意图,9/106,积事件示意图差事件示意图互不相容事件示意图对立事件示意图9/,汽车可靠性技术,第六讲,主讲教师：杨志发学时：,32,汽车可靠性技术第六讲主讲教师：杨志发,基本定理,概率加法定理,两事件不相容,两事件相容,概率乘法定理,两事件不独立,两事件独立,概率互补与全概率公式,10/106,互不相容则是说,发生则,必不发生；相互独立，是指一个事件的发生不会使另一个事件的发生增加或减少可能。,基本定理10/106互不相容则是说发生则必不发生；相互独,11/106,11/106,步骤,步骤,步骤,未,12/106,步骤步骤步骤未12/106,3,、统计特征数,平均值,中位数,n,个数据从小到大排列，居于中央位置的数，称为中位数。当,n,为奇数时，中位数 ；当,n,为偶数时，中位数是 和 的均值。,众数,在一批数据中，出现次数最多的一个数。,13/106,3、统计特征数平均值中位数众数13/106,1.,样本极差,表示离散程度的特征数,14/106,1.样本极差表示离散程度的特征数14/106,汽车可靠性技术,第七讲,主讲教师：杨志发学时：,32,汽车可靠性技术第七讲主讲教师：杨志发,二、可靠性函数,二、可靠性函数,2.1,可靠度函数,1,）可靠度函数的定义,可靠度,(Reliable),的定义：产品在规定的条件下，在规定的时间内完成规定功能的概率，称为,可靠度,。这种用概率来度量可靠度的函数，称可靠度函数，简称可靠度，用,R(t),表示。,式中,N,0, t=0,时，在规定条件下进行工作的产品数；,n,(t),在,0,到,t,时刻的工作时间内，产品的累计故障数。,16/106,2.1 可靠度函数1）可靠度函数的定义式中16/106,2,）可靠度函数的含义,产品在规定的时间内是否失效,?,是一种随机现象以用概率来度量产品在某一时刻的可靠性。可靠度是时间,t,的函数，取值在,0,与,1,之间。,T,为,产品固有寿命,?,t,为规定时间,2.1,可靠度函数,17/106,2）可靠度函数的含义T为产品固有寿命?,定义,失效概率，就是产品在规定的条件下，在规定的时间内未完成规定功能的概率，也称,不可靠度,F(t),；通常用,累积故障概率,的,分布函数来表示，这种函数，称不可靠度函数或累积失效概率分布函数，简称,失效概率分布函数,。,含义,2.2,失效概率分布函数,设有,N,件产品，在规定的条件下工作到规定的时间,T,发生故障（失效）的件数为,n(t),仍有（,N,n,）,个产品继续工作，如果产品总数,N,足够大时，则可靠度近似为,则失效概率近似为,18/106,定义2.2 失效概率分布函数设有N件产品，在规定的条件下工作,显然，以下关系成立：,失效概率密度函数,f(t),f(t),的量纲为时间的倒数，是一个广义的时间。,2.2,失效概率分布函数,19/106,显然，以下关系成立：失效概率密度函数f(t)f(t)的量纲为,可靠度函数与失效概率分布函数的性质,与 的性质如下表 所示：,对偶性,非减函数,非增函数,单调性,0,，,1,0,，,1,取值范围,2.2,失效概率分布函数,20/106,可靠度函数与失效概率分布函数的性质 与,2.2,失效概率分布函数,21/106,2.2 失效概率分布函数21/106,故障率,工作到某时刻,尚未发生故障,的产品，在该时刻后单位时间内发生故障的概率，称之为产品的,故障率,。,用数学符号表示为：,式中,故障率；,t,时刻后， 时间内,故障,的产品数；,残存产品数，即到,t,时刻尚,未故障,的产品数。,2.3,失效率函数,22/106,故障率式中 故障率； t 时刻后， 时,可按下式进行工程计算：,式中,时刻后， 时间内故障的产品数；,所取时间间隔；,残存产品数。,对于低故障率的元部件常以 为故障率的单位，称之为菲特（,Fit,）。,2.3,失效率函数,23/106,可按下式进行工程计算： 式中 ,故障率数学定义,故障率含义,产品的寿命为非负连续型随机变量,x,，失效概率分布函数为,F(t),失效概率密度函数为,f(t),时，定义 为失效率函数,2.3,失效率函数,24/106,故障率数学定义故障率含义 产品的寿命为非负连续型随机,由于 ，所以：,故障率与可靠度、故障密度函数的关系,2.3,失效率函数,25/106,故障率与可靠度、故障密度函数的关系2.3 失效率函数25/1,2.3,失效率函数,26/106,2.3 失效率函数26/106,平均故障前时间,(MTTF),2.3,失效率函数,27/106,平均故障前时间(MTTF)2.3 失效率函数27/106,平均故障间隔时间,(MTBF),2.3,失效率函数,28/106,平均故障间隔时间(MTBF)2.3 失效率函数28/106,某微型计算机的,MTBF=10000,小时，是否意味着该计算机每工作,10000,小时才出一次故障？,2.3,失效率函数,29/106,某微型计算机的MTBF=10000小时，是否意味着该计算机每,汽车可靠性技术,第八讲,主讲教师：杨志发学时：,32,汽车可靠性技术第八讲主讲教师：杨志发,寿命特征,2.3,失效率函数,30/106,寿命特征2.3 失效率函数30/106,2.3,失效率函数,31/106,2.3 失效率函数31/106,首翻期、翻修间隔期和使用寿命,(t),t,首次翻修期,规定的故障率,A,B,(=1/MTBF),翻修间隔期,使用寿命,2.3,失效率函数,32/106,首翻期、翻修间隔期和使用寿命(t)t首次翻修期规定的故障率,可靠性参数分类,可靠性参数分为基本可靠性参数和任务可靠性参数,基本可靠性反映了产品对维修人力费用和后勤保障资源的需求。,确定基本可靠性指标时应统计产品的所有寿命单位和,所有的故障,。,任务可靠性是产品在规定的任务剖面中完成规定功能的能力。,确定任务可靠性指标时仅考虑在任务期间那些影响任务完成的故障,(,即,致命性故障,),。,2.3,失效率函数,33/106,可靠性参数分类可靠性参数分为基本可靠性参数和任务可靠性参数2,2.3,失效率函数,可靠性参数的相关性,平均故障间隔时间,(MTBF),与平均故障间隔飞行小时,(MFHBF),任务成功概率与致命故障间的任务时间,MTBF,与故障率,平均维修间隔时间与,MTBF,平均拆卸间隔时间与,MTBF,34/106,2.3 失效率函数可靠性参数的相关性平均故障间隔时间(MTB,思考问题,失效率是概率值么？,故障率有量纲么？,故障率和失效率有什么关系？,2.3,失效率函数,35/106,思考问题失效率是概率值么？2.3 失效率函数35/106,分布状态描述,频数直方图,1,、概念：频数直方图是通过对随机收集的样本数据进行分组整理，并用图形描述总体分布状态 的一种常用工具,2,、,绘制程序,3,、,频数直方图、频率直方图、频率密度直方图和频率密度曲线,4,、,频数直方图的应用,2.4,直方图,36/106,分布状态描述频数直方图1、概念：频数直方图是通过对随机收,2,、绘制程序,例,从一批螺栓中随机抽取,100,件测量其外径数据如下表所示，螺栓外径规格为 ，试绘出频数直方图。,7.938,7.930,7.938,7.914,7.924,7.929,7.928,7.920,7.918,7.923,7.930,7.925,7.930,7.930,7.925,7.918,7.920,7.918,7.928,7.928,7.918,7.913,7.925,7.926,7.928,7.924,7.922,7.923,7.915,7.919,7.925,7.925,7.925,7.925,7.927,7.920,7.922,7.927,7.923,7.925,7.923,7.927,7.927,7.927,7.923,7.922,7.923,7.929,7.931,7.922,7.930,7.920,7.924,7.925,7.929,7.922,7.925,7.930,7.926,7.918,7.920,7.925,7.930,7.926,7.923,7.920,7.929,7.930,7.925,7.922,7.929,7.928,7.930,7.935,7.930,7.939,7.925,7.924,7.930,7.935,7.922,7.918,7.922,7.925,7.925,7.920,7.927,7.922,7.930,7.930,7.925,7.938,7.922,7.915,7.918,7.927,7.935,7.931,7.919,7.922,单位：,mm,2.4,直方图,37/106,2、绘制程序例 从一批螺栓中随机抽取100件测量其外径数据,绘制步骤,1,）收集数据，并找出数据中最大值,x,L,和最小值,x,S,数据个数应,50,，并将数据排成矩阵形式。本例数据个数,n,=100,。最大值,x,L,=7.938,，最小值,x,S,=7.913,。,2,）计算极差,3,）确定分组组数,k,k,值的选择一般参考下表给出的经验数值确定本例选择,k,=10,4,）,确定组距,h,5,）,计算各组的上、下边界值,6,）,计算各组的组中值,x,i,7,）,统计落入各组的数据个数，整理成频数表,8,）,作直方图,n,K,50-100,6-10,100-250,7-12,250,以上,10-20,数据量与分组数对照表,2.4,直方图,38/106,绘制步骤1）收集数据，并找出数据中最大值xL和最小值xS数据,4,）确定组距,h,组距即每个小组的宽度，或组与组之间的间隔,本例中,为分组方便，常在,h,的计算值基础上将其修约为测量单位的整数倍，并作适当调整。 如本例测量单位为,0.001,，将,h,修约为,0.003,。,2.4,直方图,39/106,4）确定组距h组距即每个小组的宽度，或组与组之间的间隔2.4,5,）计算各组的上、下边界值,为了不使数据漏掉，应尽可能使边界值最末一位为测量单位的,1/2,。,当,h,为奇数时， 第一组边界值应为,当,h,为偶数时，可以下式计算第一组边界值,第一组上边界值,=,x,S,测量单位,/2,第一组下边界值,=,上边界值,+,h,一直计算到最末一组将,x,L,包括进去为止。,本例,h,为奇数，故第一组上下边界值为,其余各组的上下边界值为：,某组上边界值,=,上组下边界值,某组下边界值,=,该组上边界值,+,h,本例第二组上下边界值为,7.9175,7.9145,；第三组为,7.9175,7.9205,依次类推，最后 一组为,7.9385,，包括了最大值,7.938(,见频数表,),。 ,7.9145,2.4,直方图,40/106,5）计算各组的上、下边界值为了不使数据漏掉，应尽可能使边界值,6,）计算各组的组中值,x,i,如本例,2.4,直方图,41/106,6）计算各组的组中值xi如本例2.4 直方图41/106,7,）统计落入各组的数据个数，整理成频数表,组号,i,组边界值,组中值,x,i,频数统计,f,i,1,7.9115,7.9145,7.913,2,2,7.9145,7.9175,7.916,2,3,7.9175,7.9205,7.919,正 正 正 一,16,4,7.9205,7.9235,7.922,正 正 正,18,5,7.9235,7.9265,7.925,正 正 正 正,23,6,7.9265,7.9295,7.928,正 正 正,17,7,7.9295,7.9325,7.931,正 正 正,15,8,7.9325,7.9355,7.934,3,9,7.9355,7.9385,7.937,4,2.4,直方图,42/106,7）统计落入各组的数据个数，整理成频数表组号i组边界值组中值,8,）作直方图,以频数为纵坐标，质量特性为横坐标画出坐标系，以一系列直方形画出各组频数，并在图中标出规格界限和数据简 历，组成频数直方图,规格要求,频数,7.937,4,3,15,17,23,18,16,2,25,20,15,10,5,0,2,7.934,7.931,7.928,7.925,7.922,7.919,7.916,7.9,7.913,7.95,n,=100 96.5.7,96.5.15,3,#,件,S,=0.00519,2,号机床,X,=7.92524,x,f,2.4,直方图,43/106,8）作直方图 以频数为纵坐标，质量特性为横坐标画出坐标,3,、频数直方图、频率直方图、频率密度直方图和频率密度曲线,频数直方图,以样本数据表征的质量特性值为横坐标，以频数为纵坐标 作出的描述数据分布规律的图形。,频率直方图,将频数直方图的纵坐标改为频率做出的频率直方图，其形状与 频数直方图应完全一样,频率直方图,若将纵坐标改为频率密度，横坐标不变，直方图的形状也不变。,频率密度曲线,当样本数据的大小,n,，组距,h0,时，直方的数量将趋于；随机变量,(,即质量特征,),在 某 区间,h,的频率密度将趋于概率密度；直方顶端联成的折线将形成一条光滑的曲线,概率密 度曲线,1,）,频数直方图、频率直方图、频率密度直方图和频率密度曲线图形演变,2.4,直方图,44/106,3、频数直方图、频率直方图、频率密度直方图和频率密度曲线 频,2,）区别与联系,频数直方图、频率直方图、频率密度直方图与概率密度曲线，虽然它们的坐标不同，描述 分布状态的方式有的是折线、有的是曲线，但其大致形状是相似的。概率密度曲线表明了总 体的分布状态；而频数直方图等是对总体分布状态的描述,3,）正态分布及其频数直方图的特征,实践和理论证明：当一个连 续型随机变量受到许多相互独立的随机因素的影响时，如果这许多因素的影响虽然有的大一 些，有的小一些，但每一个因素在影响的总和中都不起主导作用时，这个随机变量将服从正 态分布。,许多产品的计量质量指标，如强度、长度、寿命、测量误差等在生产条件稳定、正常的前 提下，均服从正态分布。因此，测量这些指标得到的数据，其频数直方图的形状应具有正态 分布概率密度曲线的特征,为中间高、两边低、左右大致对称的山峰型。,2.4,直方图,45/106,2）区别与联系2.4 直方图45/106,4,、频数直方图的应用,1,）,观察工序状态,（,1,） 原理,大部分计量指标服从正态分布，即在稳定正常生产 状态下得到的数据，其频数 直方图的形状是“中间高、两边低、 左右对称的山峰型”，我们称这种形状的直方图为,正常型直方图,。,当影响产品质量特性的因素中，有的因素在影响的总和中占据了主导地位，成为“异常 因素”时，质量特性的正态分布状态将被打破，频数直方图的形状将出现异常型。此时，现 场人员应根据直方图形状迅速分析判断异常原因，采取措施，使工序恢复正常状态。,2.4,直方图,46/106,4、频数直方图的应用1）观察工序状态2.4 直方图46/10,孤岛型, 在直方图旁边有孤立的小岛出现。其原因是在加工和测量过程中有异常情况出现。如原材 料的突然变化，刃具的严重磨损，测量仪器的系统偏差，不熟练工人的临时替班等。,偏向型, 偏向型也称偏峰型。即直方的高峰偏向一边。这常常是由于某些加工习惯造成的。如加工 孔时，有意识地使孔的尺寸偏下限，其直方图的峰则偏左；当加工孔时，有意识 地使轴的尺寸偏上限，其直方图的峰则偏右。,双峰型, 直方图出现了两个高峰。这往往是由于将不同加工者、不同机床、不同操作方法等加工 的产品混在一起造成的。因此，必须先对数据进行分层，再作频数直方图。,平顶型, 平顶型即直方图的峰顶过宽过平。这往往是由于生产过程中某种因素在缓慢的起作用造成 的。如刃具的磨损、操作者逐渐疲劳使质量特性数据的中心值缓慢的移动造成的。,折齿型, 测量误差太大或分组组数不当都会使直方图出现凸凹不平的折齿形状。,陡壁型, 直方图在某一侧出现了高山上陡壁的形状。这往往是在生产中通过检查，剔除了不合格品后 的数据作出的直方图形状。,2.4,直方图,（,2,）几种常见的,异常型频数直方图,47/106, 孤岛型 在直方图旁边有孤立的小岛出现。其原因,a.,正常型,b.,孤岛型,c.,偏向型,d.,双峰型,e.,平顶型,g.,陡壁型,f.,折齿型,规格范围,2.4,直方图,48/106,a.正常型b.孤岛型c.偏向型d.双峰型e.平顶型g.陡壁型,汽车可靠性技术,第九讲,主讲教师：杨志发学时：,32,汽车可靠性技术第九讲主讲教师：杨志发,2.4,直方图,2,）与规格比较，明确改进方向,（,1,）原理：在直方图上标明规格上限及下限，可直观地将直方图的位置、 分散范围与规格比较，从而分析质量状况，明确改进方向。,（,2,）,与规格比较的几种情况,49/106,2.4 直方图2）与规格比较，明确改进方向 49/106,产品分布范围,规格范围,产品分布范围,规格范围,产品分布范围,规格范围,产品分布范围,规格范围,产品分布范围,规格范围,最理想的直方图,直方图的分布范围仍在规格范围内，但中心偏向一侧。 此时，已存在出现不 合格品的潜在危险，应立即采取措施，将分布中心调至规格中心,直方图的分布范围已充满整个规格界限。此时，存在更多出现不合格品的潜 在危险，必须立即采取措施，减小分散。,直方图的分布范围已超出现规格界限，并已出现一定数量的不合格品。应立 即采取措施，减小分散；对产品 实施全数检查；或适当放宽规格界限，以减小损失,分布非常集中。在此情况下，应充分考虑经济效果，采取适当放宽工艺 条件或加严规格要求等措施。,50/106,产品分布范围规格范围产品分布范围规格范围产品分布范围规格范围,1,、对故障发生规律认识的变化,2.5,失效率曲线,51/106,1、对故障发生规律认识的变化2.5 失效率曲线51/106,2,、故障发生规律的六种模式,六种模式所占的比率（美国联合航空公司统计）,2.5,失效率曲线,52/106,2、故障发生规律的六种模式六种模式所占的比率（美国联合航空公,A:,浴盆曲线，（有明显的耗损故障区，如各种结构零件，轮胎、刹车片、气缸、压气机叶片等）,B:,没有明显的早期故障，在生命期内故障增长缓慢，以耗损故障终结，（有明显的耗损故障区）,C:,在生命期内故障始终缓慢增长，没有明显的耗损故障区转折点（航空涡轮发动机）,D:,刚出厂时故障率较低，然后迅速增长到一个稳定的水平,E:,整个生命周期内，故障率始终保持一个稳定的水平,F:,有早期故障，没有明显的耗损故障（航空电子设备，,software,）,定期大修是否能保持装备的可靠性水平？,事实上定期大修由于引入了新的早期故障而降低装备的可靠性水平,2.5,失效率曲线,53/106,A:浴盆曲线，（有明显的耗损故障区，如各种结构零件，轮胎、刹,3,、航天产品的统计数据,六种模式所占的比率（航天产品的统计数据）,2.5,失效率曲线,54/106,3、航天产品的统计数据六种模式所占的比率（航天产品的统计数据,4,、我国的情况,我国海军、装甲兵、通讯装备的一些统计资料都证明了许多产品都,没有明显的耗损故障区,的结论,?,2.5,失效率曲线,55/106,4、我国的情况 我国海军、装甲兵、通讯装备的一些统计资,人类健康的曲线,5,、产品故障浴盆曲线（一）,2.5,失效率曲线,56/106,人类健康的曲线5、产品故障浴盆曲线（一）2.5 失效率曲线5,5,、产品故障浴盆曲线（二）,浴盆曲线,大多数产品的故障率随时间的变化曲线形似浴盆，称之为,浴盆曲线,。由于产品故障机理的不同，产品的故障率随时间的变化大致可以分为三个阶段,：,2.5,失效率曲线,57/106,5、产品故障浴盆曲线（二） 浴盆曲线2.5 失效率曲线57/,6,、可靠性工程内涵,浴盆曲线,典型的故障率模型，许多产品故障率呈此模型,分为三个阶段：,早期寿命区：故障率随时间而降低*,恒定故障率区：故障率保持不变*,耗损区：故障率随时间而增加*,器件的缺陷，制造的缺陷或储存运输的原因造成的故障,产品进入成熟期，故障多为偶然原因造成,产品进入耗损期，故障多为老化、耗损，物理化学机理,通过浴盆曲线,确定保修期,确定使用寿命,零部件个未能更换时间,2.5,失效率曲线,58/106,6、可靠性工程内涵浴盆曲线器件的缺陷，制造的缺陷或储存运输的,7,、失效率函数曲线,1,）早期,(,初期,),失效期,早期失效期的基本特征：开始失效率较高，随着时间的推移失效率逐渐降低，因而也称递减型失效，记作,DFR,型,(Decreas-ing Failure Rate),。早期失效类型常见于使用产品的初期阶段,.,2.5,失效率曲线,59/106,7、失效率函数曲线1）早期(初期)失效期 早期失效期的,2,）偶然失效期,2.5,失效率曲线,60/106,2）偶然失效期2.5 失效率曲线60/106,3,）耗损失效期,耗损失效期的基本特征：随着时间的增长失效率急剧加大，因而也称递增型失效，记作,IFR(Increasing Failure Rate),型。,这种类型的故障是由于产品的老化、磨损、疲劳引起的，故障集中发生在某个很小的时间区间内，所以又叫做集中故障型。,2.5,失效率曲线,61/106,3）耗损失效期 耗损失效期的基本特征：随着时间的增长,4,）汽车失效率与浴盆曲线,以汽车使用的时间为横坐标，以汽车的失效率为纵坐标制成的曲线，称为失效率曲线，也称寿命曲线。由于寿命曲线的形状象浴盆，通常也称浴盆曲线。,2.5,失效率曲线,62/106,4）汽车失效率与浴盆曲线 以汽车使用的时间为横坐标，以,课堂作业,怎样制作直方图,?,课后作业,一,用浴盆曲线分析，怎样延长汽车的使用寿命？,63/106,课堂作业课后作业 一用浴盆曲线分析，怎样延长汽车的使用寿命,三、可靠性理论分布,三、可靠性理论分布,可靠度函数由不同的故障形式而服从不同的分布规律。机械设备中最常见的失效形式,随机故障服从于,指数分布,规律。,汽车可靠性研究中所用的理论概率分布类型很多，较常用的分布有：,指数分布、威布尔分布、正态分布,、超几何分布、二项分布、泊松分布、对数分布等。此外还有,2,分布、,t,分布、,F,分布等等。,3.1,指数分布,65/106,可靠度函数由不同的故障形式而服从不同的分布规律。机械,对于指数分布来说，故障是随机发生的，因而每个相同时间间隔内出现故障的几律是相等的，即指数分布中的故障率,是一个常数。,指数分布条件下，故障率,与平均故障间隔期,MTBF,互为倒数，即：,指数分布下的,可靠度,函数,:,故障概率密度函数：,f(t)=e,t,累计故障概率密度函数：,F(t)=1,e,t,3.1,指数分布,66/106,对于指数分布来说，故障是随机发生的，因而每个,0.632,F(t)=1,e,-,e,0.368,R(t)=1-,e,e,f(t),f(t)=e,-t,t (a),0,t,F(t),0,t,t (b),R(t),0,t,t,（,c,）,当工作时间,t,等于平均故障间隔期时，此时产品的可靠度为,0.368,，相应的不可靠度为,0.632,。即：,不同的检修周期将使使用中的产品具有不同的可靠度， 检修周期越短，产品的可靠度越高。,3.1,指数分布,67/106,0.632 F(t),例题：故障遵从指数分布的产品有,10,件，其故障前的工作时间分别为,30,、,70,、,110,、,120,、,140,、,250,、,320,、,360,、,400,、,450,小时，试求其平均寿命、故障率、工作,300,小时无故障的可靠度及可靠度为,99%,时的工作时间。,解：平均寿命,故障率,=0.044,工作,300,小时的可靠度,=0.263=26.3%,R=0.99,时的工作时间,3.1,指数分布,68/106,例题：故障遵从指数分布的产品有10件，其故障,汽车可靠性技术,第十讲,主讲教师：杨志发学时：,32,汽车可靠性技术第十讲主讲教师：杨志发,:,随机变量,t,平均值的估计值；,:,随机变量,t,方差的估计值。,正态分布,又称高斯分布。在数理统计的各种分布,中居于首位。设备管理中主要用于描述磨损阶段的,故障特性，轴承、齿轮、密封件等磨损件的寿命分,布大多服从正态分布。,故障概率密度函数：,3.2,正态分布,69/106,: 随机变量 t 平均值的估计值；,一般用,N (,2,),表示正态分布。参数,表示曲线对称轴距纵轴的距离，参数,则决定曲线的形状，数值越大曲线越呈扁平状。如图示：,f(t),0,t,3.2,正态分布,70/106, 一般用N (,2 ),参数,=0,、,=1,的正态分布，即,N,（,0,，,1,）称,为标准正态分布，其故障概率密度为：,标准正态分布,N,（,0,，,1,）随机变量的标准离差,u,定义为：,标准正态分布的可靠度函数：,3.2,正态分布,71/106,参数=0、=,例题：,从某批故障服从正态分布的仪器中抽取,5,台进行寿命试验，各台仪器到发生故障的时间分别为,10.5,、,11,、,11.2,、,12.5,、,12.8,（,10,3,h,），试求该批仪器工作,1210,3,h,时的可靠度。,解：,3.2,正态分布,求正态分布下的可靠度，只需求出参数,、,并据此求出参数,u,即可查表求出,（,u,）值，在此基础上计算可靠度的数值。,72/106,例题：3.2 正态分布 求正态,威布尔分布是对指数分布和正态分布的扩展，,指数分布是威布尔分布的一种特殊情况，正态分布,则可运用威布尔分布作近似计算。,对于疲劳和磨损而失效的系统，威布尔分布是最好的描述。,威布尔,分布的可靠度函数,：,3.3,威布尔分布,73/106,威布尔分布是对指数分布和正态分布的扩展，3.,故障概率密度函数：,故障率函数：,m:,形状参数；,t,0:,尺度参数。,t,m=1,m=4,m=3,m=2,f (t),0,形状参数有改变威布尔曲线形状的作用，如,f(t),曲线。当,m=2,4,时，威布尔分布的故障概率密度曲线形状接近于正态分布（,m=3.3086,时为正态分布）；,m=1,时，曲线呈指数分布。对于威布尔分布，改变,m,即可得到其它故障形式的函数式。,3.3,威布尔分布,74/106,故障概率密度函数： 故障率函数：,只需求出参数,m,、,t,0,即可计算产品的可靠度。实际应用中一般是通过威布尔概率纸加以估算。,（,1,）,威布尔概率纸的结构原理：,将式 两边取自然对数，得：,将其变形为：,两边再取对数：,令 、 、,使其具有直线方程的形式，即：,Y=AX+B,3.3,威布尔分布,75/106,只需求出参数 m、t0 即可,（,2,）,.,威布尔概率纸的结构：,1,） 两组坐标：,X,轴和,Y,轴构成的直角坐标系。,X,轴：概率纸上边线，,lnt,；,Y,轴：概率纸右边线,t,轴和,F(t),轴构成的威布尔坐标系。,t,轴：概率纸下边线，对数坐标刻度；,F(t),轴：概率纸左边线，威布尔刻度。,2,）横主轴：平行,t,轴，过,F(t)=63.2%,和,3,）纵主轴：垂直于横主轴，过,t=1,和,lnt=0,点。,4,）形状参数,m,值估算点：位于,X,轴（,1,）和,Y,轴（,0,）点。,3.3,威布尔分布,76/106,（2）. 威布尔概率纸的结构：3.3 威布尔,（,3,）威布尔概率纸的应用程序,1,）计算累积故障率,从实测数据（或按某时段的统计数）求累积故障概率,F(t),不同时间内的累积故障率为：,式中：,K,i,到时间,t,i,时的故障累计数；,n,设备或零件数。,求出,F(t,i,),后，将不同的,t,值及相应的,F,（,t,）值列表如下：,t,t,1,t,2,t,3,t,i,F(t),F(t,1,) F(t,2,) F(t,3,) F(t,i,),3.3,威布尔分布,78/106,（3）威布尔概率纸的应用程序t t1,2,）描点：,将与时间,t,相应的,F(t),值标在威布尔概率纸上并将各点拟合成为一条直线，即为威布尔直线。,3,）估算各参数值：,形状参数,m,值的估算：,将威布尔直线平移，使之,通过,m,值估算点，从通过,m,点,的直线与纵主轴的交点右引水,平线与,Y,轴相交，交点刻度的,绝对值即,m,值。,纵主轴,威布尔直线,m,值估算点,横主轴,m,值,Y,轴,3.3,威布尔分布,79/106,2）描点：纵主轴威布尔直线m值估算点横主轴m值Y轴,尺度参数 的估算：,从威布尔直线与纵主轴交点右引水平线与,Y,轴相交，将该交点所示刻度的绝对值移至,X,轴刻度上，向下作垂线与,t,轴相交所显示的刻度即,t,0,值，如图所示。,应用威布尔概率纸时，求出形状参数,m,后，可大致判断故障分布类型并可据此运用相应公式求解任意时间的可靠度参数。,a,0,0,a,t,3.3,威布尔分布,80/106,尺度参数 的估算： 应,例题：,某企业有,100,台设备，按,1,10,月份的故障统计，每月故障数和累积故障数如下表所示：试用威布尔概率纸求,m,、,t,0,、,(10),、,R(10),的值。,3.3,威布尔分布,月份,故障数,累积,故障,累 积 故,障 率,F(t),月份,故障数,累积,故障,累 积 故,障 率,F(t),1,3,3,3%,6,2,22,22%,2,3,6,6%,7,8,30,30%,3,4,10,10%,8,2,32,32%,4,5,15,15%,9,8,40,40%,5,5,20,20%,10,2,42,42%,81/106,例题：3.3 威布尔分布月份故障数累积累 积 故月份,汽车可靠性理论基础课件,解：将表中数据描在威布尔概率纸上，将各点,拟合成一条直线，如图所示，可求得：,形状参数：,m=1.35,尺度参数：,t,0,=38.5,3.3,威布尔分布,83/106,解：将表中数据描在威布尔概率纸上，将各点3.3 威布尔,威布尔分布另一种常用定义：,威布尔分布是用三个参数来描述，这三个参数分别是尺度参数,，形状参数,、位置参数,，其概率密度函数为：,(,t,，,0,，,0),3.3,威布尔分布,84/106,威布尔分布另一种常用定义：3.3 威布尔分布84/106,不同,值的威布尔分布 （,=2,，,=0,）, =1/3, =1/2, =2, =1,f(t),t,3.3,威布尔分布,85/106,不同值的威布尔分布 （=2，=0） =1/3 =1,不同,值的威布尔分布 （, =1,，,=0,）, =3, =1/2, =2, =1,f(t),t,3.3,威布尔分布,86/106,不同 值的威布尔分布 （ =1，=0） =3 =, =0, =0.5, = - 0.5, =1,f(t),t,不同,值的威布尔分布 （, =1,，, =2,）,3.3,威布尔分布,87/106, =0 =0.5 = - 0.5 =1f(t)t,四、可靠性数据分析,四、可靠性数据分析,参数估计在统计方法中的地位,统计方法,描述统计,推断统计,参数估计,假设检验,4.1,概述,89/106,参数估计在统计方法中的地位统计方法描述统计推断统计参数估计假,统计推断的过程,样本,总体,样本统计量,例如：样本均值、比例、方差,总体均值、比例、方差,抽样,参数估计,样本统计,4.1,概述,90/106,统计推断的过程样本总体样本统计量总体均值、比例、方差抽样参数,（假定身高服从正态分布 ）,设这,5,个数是,:,1.65 1.67 1.68 1.78 1.69,估计 为,1.68,，,这是,点估计,.,这是,区间估计,.,估计,在区间,1.57, 1.84,内，,例如我们要估计某队男生的平均身高,.,现从该总体选取容量为,5,的样本，我们的任务是要根据选出的样本（,5,个数）求出总体均值 的估计,.,而全部信息就由这,5,个数组成,.,91/106,（假定身高服从正态分布,参数估计的方法,矩估计法,最小二乘法,最大似然法,顺序统计量法,估 计 方 法,点 估 计,区间估计,4.1,概述,92/106,参数估计的方法矩估计法最小二乘法最大似然法顺序统计量法估 计,被估计的总体参数,总体参数,符号表示,用于估计的样本统计量,一个总体,均值,比例,方差,两个总体,均值之差,比例之差,方差比,4.1,概述,93/106,被估计的总体参数总体参数符号表示用于估计的样本统计量一个总体,1.,从总体中抽取一个样本，根据该样本的统计量对总体的未知参数作出一个数值点的估计,例如,:,用样本均值作为总体未知均值的估计值就是一个点估计,2.,点估计没有给出估计值接近总体未知参数程度的信息,3.,点估计的方法有矩估计法、顺序统计量法、最大似然法、最小二乘法等,4.2,点估计,概念,94/106,1.从总体中抽取一个样本，根据该样本的统计量对总体的未知参数,4.,用于估计总体某一参数的随机变量,如样本均值，样本比例、样本中位数等,例如,:,样本均值就是总体均值,的一个估计量,如果样本均值 ,x,= 3,，则,3,就是 的,估计值,4.2,点估计,95/106,4. 用于估计总体某一参数的随机变量4.2 点估计95/10,4.2,点估计,随机抽查,100,个婴儿,得,100,个体重数据,10, 7, 6, 6.5, 5, 5.2,呢,?,据此,我们应如何估计,和,而全部信息就由这,100,个数组成,.,例,已知某地区新生婴儿的体重,（,、,未知）,96/106,4.2 点估计随机抽查100个婴儿 ,得100个体重数据,4.2,点估计,为估计,:,我们需要构造出适当的样本的函数,T,(,X,1,X,2,X,n,),， 每当有了样本，就代入该函数中算出一个值，用来作为 的估计值,.,把样本值代入,T,(,X,1,X,2,X,n,),中，,估计值,.,T,(,X,1,X,2,X,n,),称为参数,的,点估计量,，,得到 的一个,点,97/106,4.2 点估计 为估计 :我们需要构造出适当的样本的,4.2,点估计,问题是：,使用什么样的统计量去估计 ？,可以用样本均值,;,也可以用样本中位数,;,还可以用别的统计量,.,矩估计法,极大似然法,最小二乘法,贝叶斯方法,常用估计量需求方法：,98/106,4.2 点估计问题是：使用什么样的统计量去估计 ？,1.,根据一个样本的观察值给出总体参数的估计范围,2.,给出总体参数落在这一区间的概率,3.,例如,:,总体均值落在,5070,之间，置信度为,95%,样本统计量,(,点估计,),置信区间,置信下限,置信上限,4.3,区间估计,概念,98/106,1. 根据一个样本的观察值给出总体参数的估计范围样本统计量,2,已知,2,未知,均 值,方 差,比 例,置 信 区 间,4.3,区间估计,内容,100/106,2 已知2 未知均 值方 差比 例置 信 区 间4.,x,_,X,X,=,Z,x,95%,的样本,-1.96 ,x,+1.96,x,99%,的样本,- 2.58,x,+ 2.58,x,90%,的样本,-1.65 ,x,+1.65,x,4.3,区间估计,落在总体均值某一区间内的样本,101/106,x_XX = Zx95% 的样本 -1.96,均值的抽样分布,(1 -,) %,区间包含了,%,的区间未包含,1 - a,a/,2,a/,2,4.3,区间估计,置信区间与置信水平,102/106,均值的抽样分布(1 - ) % 区间包含了1 - aa/,无偏性,：,估计量的数学期望等于被估计的总体参数,P,(,X,),X,C,A,无偏,有偏,4.4,估计量的优良准则,无偏性,103/106,无偏性：估计量的数学期望等于被估计的总体参数P( X )XC,A,B,中位数的抽样分布,均值的抽样分布,X,P,(,X,),有效性：,一个方差较小的无偏估计量称为一个更 有效的估计量。如，与其他估计量相比，样本均值是一个更有效的估计量,4.4,估计量的优良准则,有效性,104/106,AB中位数的抽样分布均值的抽样分布XP(X )有效性：一个,一致性：,随着样本容量的增大，估计量越来越接近被估计的总体参数,A,B,较小的样本容量,较大的样本容量,P,(,X,),X,4.4,估计量的优良准则,一致性,105/106,一致性：随着样本容量的增大，估计量越来越接近被估计的总体参数,谢谢大家！,谢谢大家！,