CN108305014A - 一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法 - Google Patents

一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法 Download PDF

Info

Publication number
CN108305014A
CN108305014A CN201810155798.5A CN201810155798A CN108305014A CN 108305014 A CN108305014 A CN 108305014A CN 201810155798 A CN201810155798 A CN 201810155798A CN 108305014 A CN108305014 A CN 108305014A
Authority
CN
China
Prior art keywords
failure mode
coarse
risk factors
weight
failure
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810155798.5A
Other languages
English (en)
Other versions
CN108305014B (zh
Inventor
孙岗
高建民
王宁华
张鹏飞
崔博源
陈允
王荣喜
刘贞瑶
陈轩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Xian Jiaotong University
Maintenance Branch of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Xian Jiaotong University
Maintenance Branch of State Grid Jiangsu Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, Xian Jiaotong University, Maintenance Branch of State Grid Jiangsu Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201810155798.5A priority Critical patent/CN108305014B/zh
Publication of CN108305014A publication Critical patent/CN108305014A/zh
Application granted granted Critical
Publication of CN108305014B publication Critical patent/CN108305014B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0635Risk analysis of enterprise or organisation activities
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Operations Research (AREA)
  • Marketing (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种基于可靠性屋和粗糙理想点发的故障模式与影响分析方法,以可靠性屋的矩阵结构为基础,通过专家打分的方式,获得失效模式的评分信息;采用粗糙数融合不同专家给出的评分信息,建立失效模式粗糙数评分矩阵;通过分析失效模式的传播机理,建立失效模式传播链,以失效模式的传播方向和强度作为输入,重新构建失效模式粗糙数评分矩阵,在充分考虑风险因素权重差异的情况下,引入理想点(VIKOR)法对失效模式的风险等级进行排序;相比与传统的FMEA方法及其现有的改进发放,本发明能够有效地处理专家评分信息的主观性与不确定性问题,充分考虑了失效模式间的传播影响关系,使得出的分析结果更容易被决策者接受。

Description

一种基于可靠性屋和粗糙理想点法的失效模式与影响分析 方法
技术领域
本发明属于失效模式与影响分析(FMEA)领域,涉及一种设备失效模式风险优先级辨识与分析方法,具体涉及一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法。
背景技术
失效模式和影响分析(FMEA)是一种广泛使用的工程技术,用于在产品设计和制造过程中,通过识别和分析系统及工艺过程中各种潜在失效模式、失效原因,确定其对系统所产生的影响,从而识别系统中的薄弱环节和关键项目,为制定改进控制措施提供依据。传统的FMEA采用风险优先数(Risk Priority Number,RPN)来表示失效模式对系统造成的影响,即风险等级的大小,从而有针对性地提出预防或改进措施,以保证系统的可靠性。一般来说,通常选取失效模式的发生率(Occurrence,0)、严酷度(Severity,S)和检测难度(Detection,D)来对失效模式进行评估,每个风险因素(O、S和D)有10个等级,RPN值由3个风险因素的等级数相乘确定,RPN值越大,说明该故障模式的风险等级越高,对系统的危害性越大。但传统的FMEA方法在实践中存在以下问题:
1)传统FMEA没有考虑到发生率(O)、严酷度(S)和检测难度(D)三者之间的相对重要性,认为三者的重要性相同。然而在实际中,对于不可修复系统来说,发生率D是关键影响因素,其权重不能与严重程度S和检测难易程度D等同。
2)传统FMEA是通过RPN(风险优先数)的大小来评估各故障的风险等级,由于不同的O、S、D相乘可以得到相同的RPN值,但是各种故障对于系统的风险作用却不相同。这就使得仅仅以RPN来判断各种失效模式的风险顺序在实际应用中可能会造成资源和时间的浪费,或者在某些情况下,一些高风险的失效模式反而被忽视掉。
3)PPN值的计算方法存在很大的争议,因为其对风险因素的变化非常敏感,可能一个风险因素微小的变动会导致RPN值截然不同。比如2个故障模式的O、S和D等级数分别为1、9、9和2、9、9,则它们的RPN值分别为81和162,相差2倍,但由于是发生概率的微小不同,可能对系统的危害相差并不是很大。
4)由于专家的经验和知识通常是以主观或定性的语言描述出来,而传统FMEA又很难对这些语言变量做出准确判断。
5)各专家对于同一故障给出的评价信息具有多样性和不确定性,传统的FMEA只是简单的进行加权平均,未能有效的提取各专家的评价信息。
6)传统的FMEA没有考虑到失效模式间的传播影响关系,在实际应用中,由于零部件间复杂的连接及功能偶合关系,一个零部件故障可能会引发其它零部件发生故障,即故障模式间具有传播效应。因此,在进行故障模式的风险等级评估时有必要考虑故障间的相关性,通过对故障传播链中某一故障引发后续故障的复合影响进行综合评估。
为了克服传统FMEA在实际应用中的不足,近年来,许多改进的FMEA方法相继被提出,Bowles and Peláez将模糊集理论引入到FMEA以解决传统FMEA对专家给出的语言变量很难做出准确判断的问题。Hua等针对传统FMEA中各风险因素等权重的缺陷,采用模糊层次分析法计算各风险因素的相对重要度。Hu-Chen Liu等提出了一种组合权重与模糊理想点法相结合的FMEA方法,通过引入模糊层次分析法和信息熵来共同计算各风险因素的权重。为解决传统FMEA中风险优先数计算方法缺陷造成的问题,一些学者将失效模式的风险优先级排序问题看成是多属性决策问题(multiple criteria decision-making:MCDM)。Braglia等将多属性决策问题的常用方法-优劣解距离法(TOPSIS)引入到FMEA中,提出了模糊优劣解距离法(Fuzzy TOPSIS)来对失效模式的风险优先级进行排序。杜晗恒等将模糊证据推理和TOPSIS相结合,提出了模糊TOPSIS的FMEA方法以解决影响因素难以评价的问题,该方法采用模糊置信结构对影响因素的评价进行表示,并通过去模糊化和加权平均建立明确置信矩阵。同样作为多属性决策问题的常用方法,理想点法(VIKOR)也被应用于FMEA中的失效模式风险优先级排序。Mohsen等提出了一种基于失效模式评价信息熵的VIKOR方法来对失效模式进行排序,该方法通过使用信息熵获取影响因素的客观权重,通过使用VIKOR方法获取失效模式的团体效用最大化排序、个体遗憾最小化排序以及折中排序。另外,为解决传统FMEA团队成员评估意见的多样性问题,有研究者将模糊证据推理方法、D-S证据理论和专家系统用于FMEA实施过程。Yang等将D-S证据理论应用于FMEA,以分析某型飞机发动机转子叶片在不同来源、不同评价信息下各失效模式的风险优先级。Su等为解决不同来源的评价信息间的冲突问题,提出了一种改进的D-S证据理论应用于FMEA,该方法假定评价信息服从高斯分布,通过使用不确定性推理的方法来重构评价信息的基本概率分配。
从上面的分析中可以看出,现有的FMEA改良方法大多通过将因子权重和专家模糊度引入到FMEA中以解决风险因素等权重和评价信息难以获取的问题,同时引入多属性决策方法来弥补RPN的计算缺陷,虽然在一定程度上克服了传统FMEA的不足,但是有些问题依然存在。比如专家评价信息的主观性与不确定性问题、未考虑失效模式间相关性的问题等。
发明内容
为了解决现有技术中的问题,本发明公开了一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法,有效解决了FMEA中专家评价信息主观性与不确定性以及未考虑失效模式间相关性的问题,提高了失效模式风险等级分析的准确性,能够更为有效地识别系统中的薄弱环节和关键项目,进而指导设计人员提出更为合理的预防或改进措施,以保证系统的可靠性。
为了实现以上目的,本发明所采用的技术方案为,一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法,包括以下步骤:
1)以可靠性屋(HoR)为理论基础,建立各专家的失效模式评价可靠性屋,其中包括严酷度评价指标向量、零部件失效模式向量、失效模式与严酷度评价指标关系矩阵、失效模式自相关矩阵、失效模式发生概率向量、失效模式检测难易程度向量和失效模式严酷度计算结果向量;建立风险因素的权重评价矩阵;
2)根据步骤1)获得的各专家的失效模式评价可靠性屋和风险因素权重评价矩阵,采用粗糙数方法融合各专家给出的关于失效模式的发生率(O)、严酷度(S)和检测难度(D)的评价信息以及各风险因素的权重评价信息,获得各失效模式的发生概率粗糙数、严酷度粗糙数和检测难易程度粗糙数以及各风险因素的权重粗糙数;
3)失效模式风险优先级分析:
3.1)根据步骤2)获得的各失效模式的发生概率粗糙数、严酷度粗糙数、检测难易程度粗糙数以及各风险因素的权重粗糙数,构建失效模式的粗糙数评分矩阵FM和风险因素的权重粗糙数矩阵RF:
3.2)根据步骤3.1)获得的风险因素的权重粗糙数矩阵RF计算各风险因素的权重wj
其中:wj是第j个风险因素的权重,λ是折扣因子,表示粗糙数区间的大小(即不确定度)对风险因素权重的影响程度0≤λ≤1,λ值越大,其对风险因素的权重影响越大;
3.3)根据步骤1)中给出的失效模式自相关矩阵获得失效模式间的传播方向及强度,通过考虑失效模式间的传播关系,重构失效模式的粗糙数评价矩阵FM:
其中:ait为失效模式i到失效模式t的传播强度;ati为失效模式t到失效模式i的传播强度;
重构后的失效模式粗糙数评分矩阵为:
3.4)根据步骤3.3)中获得的失效模式粗糙数评价矩阵确定各风险因素下的正理想点与负理想点(理想解是指各失效模式在各评价准则中的最佳值,而负理想解指各失效模式在各评价准则中的最差值):
3.5)计算各失效模式到各风险因素正理想点的距离dij
3.6)根据步骤3.5)中各失效模式到各风险因素正理想点的距离dij和风险因素的权重wj,计算各失效模式的团体效用最大化排序Si、个体遗憾最小化排序Ri和折中排序Qi
其中:S*=minSi,S-=maxSi,R*=minRi,R-=maxRi;v为“大多数准则”策略的权重或最大群体效用权重;v>0.5时,表示根据大多数人的意见制定决策;v=0.5时,表示根据赞同的情况制定决策;v<0.5时,表示根据拒绝的情况制定决策;
3.7)对失效模式风险等级进行排序
当下面两个条件成立时,可根据Qi的大小对失效模式进行升序排序,Qi的值越小,意味着失效模式的风险等级越高,假设fm(1)是按Qi值升序排序第一的失效模式;
条件1:可接受优势
Q(fm(2))-Q(fm(1))≥1/(m-1)
其中:Q(fm(2))为排序第二的失效模式对应的值,m为失效模式的个数;
条件2:决策过程中可接受的稳定性
根据Qi值排序第一的失效模式fm(1)必须是Si或Ri排序第一的失效模式;
如果以上条件中有一个不满足,将得到一个折衷解;
c)如果不满足条件2,则fm(1)和fm(2)均为风险等级最高的失效模式;
d)如果不满足条件1,通过
Q(fm(M))-Q(fm(1))<1/(m-1)
得到最大的M,失效模式fm(1),fm(2),…,fm(M)都贴近风险等级最高的失效模式。
进一步的,所述步骤1)包括以下步骤:
1.1)假设有K个专家组成的FMEA团队,对某一系统或设备辨识出了m个失效模式,现需对辨识出的m个失效模式的风险等级进行分析,在失效模式的评分过程中,每个专家基于可靠性屋的矩阵结构给出失效模式的评分可靠性屋;
1.2)根据失效模式的评分可靠性屋,获得失效模式评分矩阵:
同理,风险因素的权重矩阵为:
其中:FMk为第k个专家给出的失效模式评分矩阵;为第k个专家给出的第i个失效模式的第j个风险因素的评分值;K为FMEA团队中专家的个数;RFk为专家k给出的风险因素的权重向量,为专家k给出的第j个风险因素的权重值。
进一步的,所述步骤2)包括以下步骤:
2.1)将步骤1.2)中获得的失效模式评分矩阵FMk(k=1,2,…,K)中的元素以及风险因素权重矩阵RFk中的元素转换为粗糙数;
其中:为粗糙数的下边界,为粗糙数的上边界;为粗糙数的下边界,为粗糙数的上边界;
则关于风险因素j的失效模式i的评分粗糙数序列及第j个风险因素的权重粗糙数序列为:
2.2)根据粗糙数的运算规则,将获得的评分粗糙数序列和风险因素权重粗糙数序列进行融合:
其中:为融合后的失效模式评分粗糙数,为融合后的风险因素权重粗糙数。进一步的,所述步骤2.1)包括:
2.1.1)确定失效模式评分值的下近似、上近似:
同理,风险因素权重的下近似、上近似为:
2.1.2)确定失效模式评分值和风险因素权重的下边界、上边界:
其中:分别为的下近似和上近似中包含的元素的个数;分别为的下近似和上近似中包含的元素的个数;
2.1.3)失效模式评分值和风险因素权重的粗糙数为:
进一步的,所述步骤3.2)中,λ=0.5。
进一步的,所述步骤3.3)中,ait不等于ati
进一步的,所述步骤3.6)中,v=0.5。
针对传统FMEA方法的不足,本发明提出了一种基于可靠性屋和粗糙理想点法的FMEA分析方法,以可靠性屋的矩阵结构为基础,结合专家知识和经验,建立评价失效模式风险等级的评分可靠性屋,通过专家打分的方式,获得失效模式与严酷度评价指标关系矩阵、失效模式自相关矩阵、失效模式发生率向量、失效模式检测难度向量和失效模式严酷度计算结果向量,并构建各专家的失效模式评分矩阵;采用粗糙数融合各专家给出的评分信息,建立失效模式评分粗糙数矩阵,以消除专家评分信息中的主观性与不确定性;融合过程不需要依赖任何的先验知识,只依赖于专家给出的原始评分信息;采用粗糙数融合各专家给出的风险因素评分信息,建立风险因素权重粗糙数矩阵,并计算各风险因素的权重;通过分析失效模式的失效机理和零部件的功能偶合关系,建立失效模式传播链,在充分考虑失效模式传播影响关系的情况下,以失效模式传播链的传播方向和强度为输入,重新构建失效模式的粗糙数评分矩阵;以获得的计及失效模式传播影响的粗糙数评分矩阵作为最终的决策矩阵,引入理想点法(VIKOR)对失效模式的风险等级进行排序;基于VIKOR方法,本发明能在属性冲突的条件下,对失效模式的风险等级进行排序,以最接近理想解的思想给出失效模式的排序指标,在分析失效模式的风险等级时最大化群体效用,同时最小化个体遗憾,使得出的分析结果更容易被决策者接受。
附图说明
图1为本发明的方法流程图;
图2为本发明的实施例中床头箱传动系统故障树;
图3为本发明的实施例中不同零部件间故障的传播关系图;
图4为本发明方法中失效模式评分可靠性屋示意图。
具体实施方式
下面结合具体的实施例和说明书附图对本发明作进一步的解释说明。
参见图1,本发明具体包括以下步骤:
1)建立FMEA团队,并确定评价失效模式风险等级的风险因素项;
2)辨识系统或设备中潜在的失效模式(以某型号数控卧式镗车床床头箱传动系统为例):
3)基于可靠性屋的矩阵结构,对辨识出的失效模式进行评分,包含以下步骤:
3.1)确定失效模式间的传播影响关系。为获取可靠性屋中失效模式自相关矩阵中级联故障传播关系数据,采用演绎法建立床头箱传动系统故障树,故障树顶事件为传动系统故障,次级事件进一步细分为部件的故障现象,底事件为各部件具体的失效模式,参见图2;在建立的故障树中,不同零部件之间有直接相关关系的失效模式用虚线圈出,失效模式传播关系参见图3。在图3中可以观察到主要存在以下几条故障传播链:
(1)轴承同轴度差(D8)→主轴同轴度差(B02);
(2)轴承预紧力过大(B05)→齿轮啮合精度下降(D5);
(3)轴承磨损(B09)→齿轮啮合精度下降(D5);
(4)轴承同轴度差(D8)→齿轮啮合精度下降(D5);
(5)螺钉、螺母松动(B01)→接近开关未压合(B12)→主轴不转动(D3);
(6)液压油管泄漏(B14)→液压变速机构故障(D7)→主轴不转动(D3);
(7)轴断裂(B07)→主轴不转动(D3)。
3.1)参见图4,对失效模式进行评分,对于每一个专家,给出失效模式的评分可靠性屋,包括失效模式与严酷度评价指标关系矩阵、潜在失效模式自相关矩阵、失效模式发生概率向量、失效模式检测难易程度向量和失效模式严酷度计算结果向量;
3.1.1)对于失效模式严酷度计算结果向量,其计算公式为:
其中:Si为第i个失效模式的严酷度值;wj为第j个严酷度评价指标的权重;cij为第i个失效模式与第j个严酷度评价指标的相关系数;n为细分的严酷度评价指标的个数(即严酷度的子项)。
3.1.2)根据建立的失效模式评分可靠性屋,获得失效模式评分矩阵和风险因素权重矩阵:
4)采用粗糙数理论融合各专家给出的评分信息,包含以下步骤:4.1)确定失效模式评分值的下近似和上近似:
4.2)确定失效模式评分值和风险因素权重的下边界、上边界:
4.3)获得失效模式评分值和风险因素权重的粗糙数:
5)根据获得的风险因素的权重粗糙数计算各风险因素的权重wj
6)根据获得的各失效模式的评分粗糙数和风险因素的权重,建立失效模式评分矩阵和风险因素权重矩阵:
RF=[0.38 0.34 0.28]
7)根据失效模式间的传播影响关系,重构失效模式评分矩阵:
参照图3,以失效模式传播链B01→B12→D3为例,在考虑失效模式间的传播影响关系后,失效模式B01的严酷度粗糙数和失效模式B12的发生率粗糙数重构过程如下:
其中:αB01,D3=αB01,B12×αB12,D3为失效模式B01到失效模式D3的传播强度。
重构后的失效模式评分矩阵如下:
8)确定各风险因素下的理想点与负理想点
9)计算各失效模式到各风险因素理想点的距离dij
10)根据失效模式到理想点的距离dij和风险因素的权重wj,计算各失效模式的团体效用最大化排序Si、个体遗憾最小化排序Ri和折中排序Qi
11)对失效模式风险等级进行排序,根据以下根据排序规则:
当下面两个条件成立时,可根据Qi的大小对失效模式进行排序,Qi的值越小,意味着失效模式的风险等级越高。假设fm(1)是按Qi值排序第一的失效模式。
条件1:可接受优势。
Q(fm(2))-Q(fm(1))≥1/(m-1)
其中:Q(fm(2))为Qi值排序第二的失效模式对应的值,m为失效模式的个数。
条件2:决策过程中可接受的稳定性。
根据Qi值排序第一的失效模式fm(1)必须是Si或Ri排序第一的失效模式。
如果以上条件中有一个不满足,将得到一个折衷解。
a)如果不满足条件2,则fm(1)和fm(2)均为风险等级最高的失效模式。
b)如果不满足条件1,通过Q(fm(M))-Q(fm(1))<1/(m-1)得到最大的M,失效模式fm(1),fm(2),…,fm(M)都贴近风险等级最高的失效模式。
得到风险等级最高的失效模式为D5(齿轮啮合精度下降)。
本发明针对传统FMEA的不足,发明了一种改进的FMEA分析方法,通过引入可靠性屋(HoR)、粗糙集理论和理想点法(VIKOR)对失效模式的风险等级进行分析;首先基于可靠性屋的矩阵结果,结合专家知识和经验建立评价失效模式风险等级的评分可靠性屋,获得失效模式与严酷度评价指标关系矩阵、失效模式自相关矩阵、失效模式发生概率向量、失效模式检测难易程度向量和失效模式严酷度计算结果向量。根据获得的评分可靠性屋,采用粗糙数融合各专家给出的评分信息,建立失效模式粗糙数评分矩阵。在融合过程中,基于粗糙数灵活的处理机制,本发明能够有效地消除专家评分信息的主观性与不确定性,而且不需要依赖任何的先验知识,只依赖于专家给出的原始评价信息;在分析失效模式的风险等级时,本发明综合考虑了风险因素权重的差异性以及失效模式间的传播影响关系,通过分析失效模式的失效机理和零部件间的拓扑耦合关系,建立失效模式传播链,根据失效模式的传播方向和强度重新构建失效模式评分矩阵。最后,引入了理想点法(VIKOR)对失效模式风险等级进行排序,其能在属性冲突的条件下,对失效模式的风险等级进行排序,以最接近理想解的思想给出失效模式的排序指标,在分析失效模式的风险等级时最大化群体效用,同时最小化个体遗憾,使得出的分析结果更容易被决策者接受。

Claims (7)

1.一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法,其特征在于,包括以下步骤:
1)以可靠性屋为理论基础,建立各专家的失效模式评价可靠性屋,其中包括严酷度评价指标向量、零部件失效模式向量、失效模式与严酷度评价指标关系矩阵、失效模式自相关矩阵、失效模式发生概率向量、失效模式检测难易程度向量和失效模式严酷度计算结果向量;建立风险因素的权重评价矩阵;
2)根据步骤1)获得的各专家的失效模式评价可靠性屋和风险因素权重评价矩阵,采用粗糙数方法融合各专家给出的关于失效模式的发生率(O)、严酷度(S)和检测难度(D)的评价信息以及各风险因素的权重评价信息,获得各失效模式的发生概率粗糙数、严酷度粗糙数和检测难易程度粗糙数以及各风险因素的权重粗糙数;
3)失效模式风险优先级分析:
3.1)根据步骤2)获得的各失效模式的发生概率粗糙数、严酷度粗糙数、检测难易程度粗糙数以及各风险因素的权重粗糙数,构建失效模式的粗糙数评分矩阵FM和风险因素的权重粗糙数矩阵RF:
3.2)根据步骤3.1)获得的风险因素的权重粗糙数矩阵RF计算各风险因素的权重wj
其中:wj是第j个风险因素的权重,λ是折扣因子,表示粗糙数区间的大小对风险因素权重的影响程度,0≤λ≤1,λ值越大,其对风险因素的权重影响越大;
3.3)根据步骤1)中给出的失效模式自相关矩阵获得失效模式间的传播方向及强度,通过考虑失效模式间的传播关系,重构失效模式的粗糙数评价矩阵FM:
其中:ait为失效模式i到失效模式t的传播强度;ati为失效模式t到失效模式i的传播强度;
重构后的失效模式粗糙数评分矩阵为:
3.4)根据步骤3.3)中获得的失效模式粗糙数评价矩阵确定各风险因素下的正理想点与负理想点
3.5)计算各失效模式到各风险因素正理想点的距离dij
3.6)根据步骤3.5)中各失效模式到各风险因素正理想点的距离dij和风险因素的权重wj,计算各失效模式的团体效用最大化排序Si、个体遗憾最小化排序Ri和折中排序Qi
其中:S*=min Si,S-=max Si,R*=min Ri,R-=max Ri;v为“大多数准则”策略的权重或最大群体效用权重;v>0.5时,表示根据大多数人的意见制定决策;v=0.5时,表示根据赞同的情况制定决策;v<0.5时,表示根据拒绝的情况制定决策;
3.7)对失效模式风险等级进行排序
当下面两个条件成立时,可根据Qi的大小对失效模式进行升序排序,Qi的值越小,意味着失效模式的风险等级越高,假设fm(1)是按Qi值升序排序第一的失效模式;
条件1:可接受优势
Q(fm(2))-Q(fm(1))≥1/(m-1)
其中:Q(fm(2))为排序第二的失效模式对应的值,m为失效模式的个数;
条件2:决策过程中可接受的稳定性
根据Qi值排序第一的失效模式fm(1)必须是Si或Ri排序第一的失效模式;
如果以上条件中有一个不满足,将得到一个折衷解;
a)如果不满足条件2,则fm(1)和fm(2)均为风险等级最高的失效模式;
b)如果不满足条件1,通过
Q(fm(M))-Q(fm(1))<1/(m-1)
得到最大的M,失效模式fm(1),fm(2),…,fm(M)都贴近风险等级最高的失效模式。
2.根据权利要求1所述的一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法,其特征在于,所述步骤1)包括以下步骤:
1.1)假设有K个专家组成的FMEA团队,对某一系统或设备辨识出了m个失效模式,现需对辨识出的m个失效模式的风险等级进行分析,在失效模式的评分过程中,每个专家基于可靠性屋的矩阵结构给出失效模式的评分可靠性屋;
1.2)根据失效模式的评分可靠性屋,获得失效模式评分矩阵:
同理,风险因素的权重矩阵为:
其中:FMk为第k个专家给出的失效模式评分矩阵;为第k个专家给出的第i个失效模式的第j个风险因素的评分值;K为FMEA团队中专家的个数;RFk为专家k给出的风险因素的权重向量,为专家k给出的第j个风险因素的权重值。
3.根据权利要求2所述的一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法,其特征在于,所述步骤2)包括以下步骤:
2.1)将步骤1.2)中获得的失效模式评分矩阵FMk(k=1,2,…,K)中的元素以及风险因素权重矩阵RFk中的元素转换为粗糙数;
其中:为粗糙数的下边界,为粗糙数的上边界;为粗糙数的下边界,为粗糙数的上边界;
则关于风险因素j的失效模式i的评分粗糙数序列及第j个风险因素的权重粗糙数序列为:
2.2)根据粗糙数的运算规则,将获得的评分粗糙数序列和风险因素权重粗糙数序列进行融合:
其中:为融合后的失效模式评分粗糙数,为融合后的风险因素权重粗糙数。
4.根据权利要求3所述的一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法,其特征在于,所述步骤2.1)包括:
2.1.1)确定失效模式评分值的下近似、上近似:
同理,风险因素权重的下近似、上近似为:
2.1.2)确定失效模式评分值和风险因素权重的下边界、上边界:
其中:分别为的下近似和上近似中包含的元素的个数;分别为的下近似和上近似中包含的元素的个数;
2.1.3)失效模式评分值和风险因素权重的粗糙数为:
5.根据权利要求1所述的一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法,其特征在于,所述步骤3.2)中,λ=0.5。
6.根据权利要求1所述的一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法,其特征在于,所述步骤3.3)中,ait不等于ati
7.根据权利要求1所述的一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法,其特征在于,所述步骤3.6)中,v=0.5。
CN201810155798.5A 2018-02-23 2018-02-23 一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法 Expired - Fee Related CN108305014B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810155798.5A CN108305014B (zh) 2018-02-23 2018-02-23 一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810155798.5A CN108305014B (zh) 2018-02-23 2018-02-23 一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法

Publications (2)

Publication Number Publication Date
CN108305014A true CN108305014A (zh) 2018-07-20
CN108305014B CN108305014B (zh) 2021-12-03

Family

ID=62848807

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810155798.5A Expired - Fee Related CN108305014B (zh) 2018-02-23 2018-02-23 一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法

Country Status (1)

Country Link
CN (1) CN108305014B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109657260A (zh) * 2018-09-19 2019-04-19 北京航空航天大学 一种考虑失效相关性的涡轮转子系统可靠性分配方法
CN110222375A (zh) * 2019-05-13 2019-09-10 北京航空航天大学 一种舰载机着舰过程的安全监测方法
CN112395766A (zh) * 2020-11-24 2021-02-23 华能新能源股份有限公司 一种考虑故障传播影响的风电机组故障模式与影响分析方法
CN112685825A (zh) * 2021-01-22 2021-04-20 西安航空职业技术学院 一种逐步等效平面法的优化方法
CN113688527A (zh) * 2021-08-27 2021-11-23 沈阳工业大学 一种考虑指标关联性与评价模糊性的激光增材制造装备工效学优化设计效果综合评价方法
CN115422620A (zh) * 2022-11-08 2022-12-02 北自所(北京)科技发展股份有限公司 基于数字孪生的物流装备再设计方法、装置及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103186708A (zh) * 2011-12-31 2013-07-03 电子科技大学 一种采用两个rpn数的故障模式影响及危害性分析方法
CN103279640A (zh) * 2013-04-23 2013-09-04 中国电子科技集团公司第三十六研究所 基于粗糙集理论和可靠性预计的星载电子系统fmeca方法
CN107133727A (zh) * 2017-04-22 2017-09-05 南昌航空大学 基于鉴别力系数的失效影响评估维度优化方法
CN107239626A (zh) * 2017-06-13 2017-10-10 南京理工大学 一种基于焊接机器人的模糊fmeca分析方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103186708A (zh) * 2011-12-31 2013-07-03 电子科技大学 一种采用两个rpn数的故障模式影响及危害性分析方法
CN103279640A (zh) * 2013-04-23 2013-09-04 中国电子科技集团公司第三十六研究所 基于粗糙集理论和可靠性预计的星载电子系统fmeca方法
CN107133727A (zh) * 2017-04-22 2017-09-05 南昌航空大学 基于鉴别力系数的失效影响评估维度优化方法
CN107239626A (zh) * 2017-06-13 2017-10-10 南京理工大学 一种基于焊接机器人的模糊fmeca分析方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
GUO JIAGANG,ZHANG WEI: "Selection of Suppliers Based on Rough Set Theory and VIKOR Algorithm", 《2008 INTERNATIONAL SYMPOSIUM ON INTELLIGENT INFORMATION TECHNOLOGY APPLICATION WORKSHOPS》 *
伍晓榕: "零件工艺失效模式关联建模与工艺方案绿色性分析关键技术研究", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 *
石洁 等: "基于FMECA/HoR的气动薄膜调节阀可靠性分析", 《轻工机械》 *
程伟平: "基于QFD汽车离合器设计方法及评价决策研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109657260A (zh) * 2018-09-19 2019-04-19 北京航空航天大学 一种考虑失效相关性的涡轮转子系统可靠性分配方法
CN110222375A (zh) * 2019-05-13 2019-09-10 北京航空航天大学 一种舰载机着舰过程的安全监测方法
CN112395766A (zh) * 2020-11-24 2021-02-23 华能新能源股份有限公司 一种考虑故障传播影响的风电机组故障模式与影响分析方法
CN112395766B (zh) * 2020-11-24 2023-10-03 华能新能源股份有限公司 一种考虑故障传播影响的风电机组故障模式与影响分析方法
CN112685825A (zh) * 2021-01-22 2021-04-20 西安航空职业技术学院 一种逐步等效平面法的优化方法
CN112685825B (zh) * 2021-01-22 2024-06-11 西安航空职业技术学院 一种逐步等效平面法的优化方法
CN113688527A (zh) * 2021-08-27 2021-11-23 沈阳工业大学 一种考虑指标关联性与评价模糊性的激光增材制造装备工效学优化设计效果综合评价方法
CN115422620A (zh) * 2022-11-08 2022-12-02 北自所(北京)科技发展股份有限公司 基于数字孪生的物流装备再设计方法、装置及存储介质

Also Published As

Publication number Publication date
CN108305014B (zh) 2021-12-03

Similar Documents

Publication Publication Date Title
CN108305014A (zh) 一种基于可靠性屋和粗糙理想点法的失效模式与影响分析方法
CN105096053B (zh) 一种适用于复杂工艺系统的健康管理决策方法
CN103337043B (zh) 电力通信设备运行状态的预警方法和系统
CN104573740B (zh) 一种基于svm分类模型的设备故障诊断方法
Guo et al. Criticality evaluation of petrochemical equipment based on fuzzy comprehensive evaluation and a BP neural network
CN106447205A (zh) 一种基于层次分析法的配电自动化终端状态评价方法
CN106650122B (zh) 一种设备变工况运行风险评估方法
Jiang et al. Fault diagnosis for rolling bearing using a hybrid hierarchical method based on scale-variable dispersion entropy and parametric t-SNE algorithm
CN108985649A (zh) 基于风电设备异常分析数据的时序模糊层次分析评估方法
CN111797365B (zh) 一种换流变压器温度异常判断方法及系统
CN104123678A (zh) 一种基于状态等级评估模型的电力继电保护状态检修方法
Wang et al. Evaluation of EDM process for green manufacturing
CN106295332A (zh) 基于区间数和理想解的信息安全风险评估方法
CN115187832A (zh) 一种基于深度学习与格拉姆角场图像的能源系统故障诊断方法
CN108009937A (zh) 一种配电主设备健康状态的评估方法
CN109144028A (zh) 一种精馏塔能效退化检测方法
CN115949890A (zh) 一种城市燃气管网泄漏监测分级报警和处置方法
CN112395766B (zh) 一种考虑故障传播影响的风电机组故障模式与影响分析方法
CN108171397A (zh) 一种配电二次设备状态风险评估方法
CN111312414A (zh) 核电厂疲劳监测和寿命评估系统的疲劳评价方法
Liu et al. Risk assessment based on failure mode and effects analysis (FMEA) and WASPAS methods under probabilistic double hierarchy linguistic term sets
Jeni Seles Martina et al. Operations on multi-valued neutrosophic matrices and its application to neutrosophic simplified-TOPSIS method
CN108761263A (zh) 一种基于证据理论的故障诊断系统
CN204808407U (zh) 储气库地面设备的风险评估系统
CN110289112A (zh) 基于层次分析和模糊评价的核动力装置健康状态诊断方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20211203