故障的统计分析典型的故障率分布曲线

1、题目：故障的统计分析与典型的故障率分布曲线学号：120606325 姓名：王逢雨摘 要机械故障诊断是一门起源于20世纪60年代的新兴学科，其突出特点 是理论研究与工程实际应用紧密结合。 该学科经过半个世纪的发展逐渐成熟， 在 信号获取与传感技术、故障机理与征兆联系、信号处理与诊断方法、智能决策与 诊断系统等方面形成较完善的理论体系，涌现了如全息谱诊断、小波有限元裂纹 动态定量诊断等原创性理论成果，在机械、冶金、石化、能源和航空等行业取得 了大量卓有成效的工程应用。统计分析工作是机械故障诊断中的核心环节，统计分析工作的质量和水平将会对机械设备的检修工作产生重要影响，关系到机械设备的安全与可靠运行

2、。本文在对机械故障的特性等问题进行阐述的基础上，重点就机械故障统计分析工作中数据的收集和统计分析的方法进行重点探讨，希望对提高机械故障的管理水平能够有所帮助。关键词机械故障；统计分析；数据收集；方法一、统计分析工作中机械故障的特性机械设备在使用过程中，由于会受荷载应力等环境因素的影响， 随着机械设备部 件之间磨损的不断增加，结构参数与随之变化，进而会对机械功能的输出参数产 生影响，甚至使其偏离正常值，直至产生机械故障。概括说来，主要有以下几方 面的特性。（一）耗损性在机械设备运行过程中，不断发生着质量与能量的变化，导致设备的磨损、疲劳、 腐蚀与老化等，这是不可避免的，随着机械设备使用时间延长，

3、故障发生的概率 也在不断增加，即使可以采取一定的维修措施，但是由于机械故障的耗损性，不 可能恢复到原先的状态，在经过统计分析工作后，必要时需要对设备进行报废。（二）渐损性机械故障的发生大多是长期运行的老化或疲劳引起的，所以具有渐损性，而且与设备的运行时间有一定的关系，所以做好机械设备的统计分析工作是很有必要 的，当掌握了设备故障的渐损规律后， 可以通过事前监控或测试等手段， 有效预 防机械故障的发生。（三）随机性虽然有的机械故障具有一定的规律性， 但这并不是绝对的，因为机械故障的发生 还会受到使用环境、制造技术、设备材料、操作方式等多种因素的影响，因此故 障的发生会具有一定的分散性和随机性，这

4、在一定程度上增肌了机械设备预防维 修与统计分析工作的难度。（四）多样性随着科学技术的发展与应用，机械设备的工作原理日趋复杂，零部件的数量在不 多增多，这就使得机械故障机理发生的形式日趋多样化。 机械故障的发生不仅存 在多种形式，而且分布模型及在各级的影响程度也不同， 在统计分析工作中需要 引起足够的重视。二、机械故障管理中统计数据的收集在对机械故障的统计分析工作中，数据的收集是最基础的环节，因此必须保障数 据收集的及时性、准确性和规范性，这样才能为接下来的数据分析工作奠定良好 的基础。（一）做好日常检点数据的收集机械设备的操作人员和统计人员要重视对日常点检记录数据的收集，只有这样才能保证统计数

5、据收集的全面性。对此，相关人员要严格按照点检表对设备进行检 查和记录，对于发生的问题或故障要在第一时间记录在 设备故障报修单上, 并及时上交机械故障管理部门。（二）安装调阶段相关数据的统计与处理对于机械设备在安装调试阶段发生的故障， 统计人员要引起高度重视，并详细记 录在案，以备后期的参考与分析，对于搜集到的同行业相同设备的故障数据， 在 确保其真实性的基础上，也可以将其纳入到统计范围之内。（三）确保采集与整理的规范性为了保证故障数据分析的准确性和规范性， 统计人员及设备检修人员在日常工作 中必须详细、规范地填写 设备故障分析报告、设备故障处理单 等资料，对机械故障发生的部位、原因、时间、表现

6、以及后期的处理与改进等详细 记录在案。三、机械故障管理中统计数据的分析方法（一）焦点分析法焦点分析法是一种最直接、最简单的方式，是以机械故障问题点为中心的分析方 法，其分析结果简单明了，实用性比较强。首先，我国要根据需要把一个圆分成 等分成若干块，每一块分别代表着生产线机械设备有标准化问题点的一部分，分别记录着该部分发生故障的次数，然后用有量线段进行表示，最后将这些点进行 连接，所形成的多边形就是带有评价性质的焦点分析图。（二）直方图对比分析法该方法要求预先对计划指标数值进行设定，然后按照机械故障发生的实际录入实 绩值，然后将实绩值与计划值进行对比， 看其差距之间的大小，并参照历史实绩 值进行

7、分析，这可以反映出机械故障发生时计划值与实绩值的科学性，以及设备故障发生概率的大小，便于及时采取相应的检修措施。（三）排列图分析法排列图分析法也被称为帕洛特图法、 主次因素分析法，它是找出造成设备故障并 进行分析的一种简便有效的图表分析方法。排列图是根据“关键的少数和次要的 多数”的原理而制作的。即对影响机械设备故障的因素按照影响程度的大小用直 方图进行排列，找出最主要的因素，其结构包括一个横坐标和两个纵坐标， 若干 个直方形和一条折线构成，通常将影响因素分为三类：A类因素（占比80%以下）、B类因素 （占比80%90% 、C类因素 （占比90%以上）。其中， A类因素为主要因素，也是设备故障

8、管理中需要重点解决的因素。四、典型的故障率分布曲线现代的设备管理中，典型的故障率分布曲线一一浴盆曲线仍然占有很重要的 地位。很多故障的分析都是基于浴盆曲线发展的。无故障工作期就是在浴盆曲线上发展而来的。 与传统可靠性指标中假设产品的随机失效不可避免不 同，无故障工作期（FFOP） 内产品不会发生任何故障(即零故障)。首先阐述了 FFOP的概念内涵、 与平均故 障间隔时间(MTBF) 的区别和联系， 提出了一种 FFOP的预计方法。 该方法假设产品的故障率函数具有浴盆曲线特征、故障发生过程为泊松过程、产品具有固定的免维修工作期。然后以一种改进的We 1 b u 1 1分布函数描述具有浴盆曲线函数

9、特征产品的故障率。 基于泊松过程理论， 给出了 F FOP 的预计算法、流程和仿真验证手段。最后以某型无人机舵机为案例对研究方法的可用性进行了验证。结果表明：F FOP 与免维修工 作期(MFOP)、置信度水平密切相关，及时维修的产品能够保证较长的 F FOPo在工程应用时， F FOP的确定应综合考虑运行维护费用进行权衡。无故障工作期(Fai lureFreeOp eratin gPerio d , FFOP)定义为产品不会 发生任何故障(即零故障)的时 问。对于符合设计要求、 质量合格的产品， 往往都要求其具有一定的无故障工 作期，尤其是具有高可靠性/安全性需求的系统，如武器装备、 核能系

10、统、 载人航空航天器、 高速列车等。作为耐久性度量指标，FFOP的长短与维修费用、 保障费用紧密相关。准确预计FFOP,结合 合理的维修策略，能够实现对产品的充分使用，降低运行成本。F FOP概念最早在美国空军颁布的军用规范MI L A87244航空电子设备完整性大纲要求 中提出1,其中FFOP作为耐久性 参数，对传统的可靠性参数进行了补充， 并指导设计和生产。后来美国又 颁布了一系列规范和指南， 都对装备的FFOP指标有了明确的要求2 3。在 1 9 9 6 年英国国防部( Ministr yofDe fence, M O D) 提出免维修工作期( Ma i ntenanceF reeOp

11、-era tin gPeriod, MFOP) 的概念以后,FFOP就通常与 MFOP结合度量产品的耐久性。MF OP概念比FFOP严 格，在MFOP内，产品不允许出现任何影响性能和任务的 失效事件；而FFOP内不允许故障但允许维修活动，F F OP是一系列免维修工作期的集合。文分析了英国国防 部为新一代战机提出的MFOP概念，与平均故障间隔时 问(MeanTimeBetweenFailures, MTBF)进行了 对比，并分别研究了基于任务可靠度和更新理论的MF O P预计方法，于英国的超高可靠飞行器(UltraReli ableAir craft, URA) 和未来攻击飞行器(Fu tur

12、eOf fen s iveAi rcraft, FOA) 项目。当前国 内外的研究大多集中在对无故障工作期/免维修工作期 (F /M-FOP)概 念的阐述以及适场合分析等方面17 11,证明了基于F/ MF O P维修策略的有效性。文献12文献14假设产品故障为有限时间区间内的离散事件， 基于统计方法估计了产品存在某固定长度MFO P的概率。文献15以典型机电产品为案例， 研究故障事件为齐次泊 松过程情况下 FFOP的评估方法，并对结果进行了合理性分析。文献1 6 和文献17基于Petri网络，使用仿真方法分析了周 定MFOP系统的可靠度。以上研究集中在MFOP预计方法方 面，没有考虑维修策略

13、对 FFOP的影响。然而， 为促进基于FFO P维修策略的应用，需要进一步研究 F F O P的预计方法与模型。 在很多情况下，产品(系统)的F /MFOP大多由运行过程中 随机故障事件之间的相对位置(时间、空间距离) 决定，相对 位置的远近直接影响产品的 FFOP。 以图1所示的时间(空间) 区间0, L为例，假设系统是一个客户服务系统， 为一个客户服务的 免维修周期为s。如果两个或者更多的客户集中在s内出现， 如图1 ( a)所 示，则系统会出现过载(故障)，此系统的FFOP为s的概 率就是Pnt, t + s 01, nt, t + s为t, t + s区 间内的客户数量。类似的方法也可

14、以用于分析交通处理系统， 如图1 ( b)所 示，如果一个交通意外的恢复周期为s,在这段周期内出 现的其他意外则会导致拥堵(故障)；如果把事件区间换作一段钢结构(见图1 ( c) 或者电缆(见图1 ( d),也存在一个极限区间 s,在这个区间内应力集中点或缺陷次数要低于某一确定数量， 否则会出现故障。以上案例中，客户出现与事故发生时刻Random demands(cuitumers，4-*0(a) Minimum processing periodfailureslb) Minimum cecoven periodKandatn sire%(cl Maximum acceptable suct

15、ion bel ccn sirensRandom(di Miiximum acccptibli! distaiicj between defects图1随机事件聚集导致系统故障Fi耳 1 System failures due to clustering o： random events泊松过程是描述随机事件发生的基本数学模型之一，实际生活或自然世界中的随机事件，大多可以用泊松过程描述1 8 。对于寿命服从指 数 分布的产品， 故障率是一个常数， 寿命周期内随机故障事件可以用齐次泊 松过程描述。然而，实践证明，大多数产品的故障率随时间变化的曲线是浴盆 曲线19,故障率是时变函数，故障事件需要用

16、非齐次泊 松过程描述。本文首先阐述 FFOP 与MFOP 之 间的区 别与联系，然 后提出一种FFOP预计方法，预计故障率函数为浴盆曲线的产品 的无故障工作期。该方法作了如下假设： 故障事件服从泊松过程； 故障率函数为浴盆曲线；F FOP内允许固定周期的计划维修， 产品 修复如 新；一个MFOP内不允许有任何影响产品正常运行的故障事件，一个维修恢复期(Ma i ntenanceRecove r yPeriod, MRP) 只能处理一次随机故障。在以上假设的基础上， 给出了 FFOP的预计方法、 模型和预计步骤， 并通过某型无人机舵机对所提方法进行了应用验证。1 FFOP概念与内涵在 MIL A

17、87244中， FFOP 被定义为故障概率达到2 % 的时 间。图2描述了概率密度函数(Probabilit y Den sit y Function, PDF)、F F O P这3者之间的区别与联系。图2 寿命分布与FFOPFig. 2 Life distribution and FFOP根据FFOP和MTBF的定义，有/OFFOPf( t) dt = 2%(1)MTBF = /0 8R ( t) d t=/ 0/ootoof(r)drdt(2)式中：f ( t)为故障密度函数； R ( t)为可靠度函数。对于大多数产 品来说，由于不可避免 的随机失效，图2所示时间t 0通常为0,这样 就导

18、致产品的FFOP很短。然而对于具有高可靠性/安全性需求的系统，又需要具有一定长度的FFOP。这个要求既可以通过设计手段降低产 品的故障率实现，对于可修复产品，又可以通过固定周期的 维护，使产品始终工作在比较“新”的状态，进而降低随机故障事件发生的概率来实现。对于可修复的产品，FFOP 与MFOP 密切相关10。如果维护频繁，并且能够保证修复如新的话，F FOP会比维护不力的设备要长。建立FFOP预计模型是预计FFOP的关键步骤。若要使产 品在整个工作周期0, L内无故障运行，则要求在每次故障发生前进行 维护并恢复到完好状态。由于一个维修恢复期只能处理一次随机故障，所以要求维修次数要和随机故障的

19、次数一致，并且在故障事件实际发生之前就已经得到维修并完全修复，即第1次和第1 + 1次维修之间的间隔时间s1,1 +1小于第1次 和第1 + 1次实际故障间隔时间S1 , 1 + 1o若在整个寿命周期0, L内出现k 次故障， 设定免维修工作期 MFOP 1, 1 + 1 = s i , i + 1 , 那么存在长度为L的FFOP的概率PF FOP (故障发生前都能被完全修复以避免故障实际发生的概率)为1 5 PFFOP =P( S 0, 1 0S 0, 1 ns 1, 2 0S 1, 2 n ns k 1, k t)=e 入 t (5)假设随机事件是故障事件，在t时刻，随机故障事件导致的系统

20、不可靠度为F(t)=P(Tt)=le 入 t (6)2. 2浴盆曲线的故障率函数已有的研究成果表明， 基于浴盆曲线的故障密度函数有如下形式20: f(t)=YB(t/a) Blexp(t/a) B+ Ya(lexp(t/a) B)(7)对应的可靠度函数为R(t)=exp(Ya(lexp(t/a) B)(8)故障率函数为入(t)=YB(t/a) p lex p(t/a)0)(9)式中：a、 B、 丫均为分布函数中的参数。绘制故障率函数曲 线，如图3所示0.050.04S03002*020406080Tiine/h图3 故障率函数曲线Fig. 3 Plot of failure rate func

21、tion从图3可以看出， 产品的故障率明显 呈浴 盆曲线特性， 可以描述分布特征 为浴盆曲线的产品故障率。2.3 F F O P的预计步骤研究具有浴盆曲线故障率函数的产 品，与 指数分布不同，其故障率为非常值, 且导致故障发生为非齐次泊松过程， 对比文献15中P F F O P的计算 公式，可以得到(10)式中：r为允许的维护次数。在进行FFOP预计之前，需要根据式（7） 式（9）确定产品的入（t ）。F F O P的预计步骤如图4图4 FFOP预计步骤Fi珞 4 Prediction steps of FFOP首先，设置FFOP为L的置信度PFFOP也和设置初始维护 次数r = 1。按照图4所示的流程， 对维护次数递增， 得到满足式（11）的最大维护次 数r。f*FF()P XP |l - J A(Z)d/A(Z)d/(11)(1-r)F F 0 P 的估计区间为 FFOP rs,（r + l） s （12）对于大多数工程应用，式（12）所描述的FFOP区间已经足够。 更精确的预计结果可以通过在区间内多点取值，由式（11）反复校验的方式获取。基于浴盆曲线故障率函数的FFOP预计方法，能够预计失效过 程为泊松过程， 并且故障率函数服从浴 盆曲线特征情况下的产品 无故障工作期。将FF