CN104898039A - 基于故障传播概率模型的故障模式优选方法 - Google Patents
基于故障传播概率模型的故障模式优选方法 Download PDFInfo
- Publication number
- CN104898039A CN104898039A CN201510280067.XA CN201510280067A CN104898039A CN 104898039 A CN104898039 A CN 104898039A CN 201510280067 A CN201510280067 A CN 201510280067A CN 104898039 A CN104898039 A CN 104898039A
- Authority
- CN
- China
- Prior art keywords
- fault
- fault mode
- mode
- pervasion
- probability
- 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
Links
Landscapes
- Test And Diagnosis Of Digital Computers (AREA)
Abstract
本发明公开了一种基于故障传播概率模型的故障模式优选方法,首先对系统各模块的故障样本量进行分配,然后对于各个模块建立模块的故障传播概率模型,得到各故障模式之间的故障传播概率,然后根据故障传播概率和故障模式所在元件节点的出度计算得到故障模式间的故障扩散强度,根据故障扩散强度搜索出每个故障模式作为原发故障时的故障扩散概率最大路径,将该路径中包括的后继故障模式数量作为故障扩散数量,筛选出故障扩散数量大于预定阈值的故障模式,优先进行选择;最后将各个模块的选定故障模式集合得到系统的选定故障模式集合。本发明基于故障扩散数量进行故障模式的优选,对传播型故障进行重点抽样,提高测试性验证方案可靠性。
Description
技术领域
本发明属于测试性验证技术领域,更为具体地讲,涉及一种基于故障传播概率模型的故障模式优选方法。
背景技术
测试性验证是依据故障检测、隔离和指示结果进行的,若仅对产品自然发生的故障进行检测,这是达不到设计规定的样本量要求的,因此需要人为的引入故障即故障模拟。在开展试验验证之前,需要确定验证的总的工作流程以确定试验样本量,指示故障注入及数据记录与评估。
进行测试性验证,首先必须选取试验的样本量,样本量的选取是随机抽取。如果选取n个样本进行试验,试验失败F次。规定一个正整数C,如果F<C则认为合格,判定接收;如果F≥C则认为不合格,判定拒收。确定抽样方案就是同时确定n和C的值。
在这种抽样试验中,设成功的概率记为q,则在n次试验中出现F次失败的概率为:
接收的概率为n个样本中失败数不超C的概率总和:
在承制方和使用方协商确定FDR(Fault Detection Rate,故障检测率)的设计值q0,FDR的最低可接受值q1,生产方风险α,使用方风险β的基础上,就可以用下面的联立不等式组确定抽样方案,即确定n和C值。
由于实际应用中n和C必须为正整数,正好满足以上方程组的正整数不一定存在,所以一般应该求取满足方程最小的n和C。
由于系统采用分层设计的模式,所以当测试性验证方案的抽样方案确定之后,需要将总的样本量按比例分配到每个模块,然后在每个模块中抽取故障模式组成注入系统的故障模式集。之所以要把样本量分配到各模块中去,是因为复杂的系统的故障模式十分庞大,设故障样本总量为M个故障模式,若从总M个故障模式中直接抽取n个故障模式将会造成故障模式分配不均匀,样本结构不合理等特点,因此通常采用分层分配的方式。分层分配故障样本量的公式为:
式中,nr表示系统中第r个模块分配的故障样本数量,r的取值范围为1≤r≤R,R表示系统中的模块数量,Wr表示系统中第r个模块的分配权重值;Qr表示第r个模块中故障模式的数量,用来表征模块的复杂性;λr表示第r个模块的故障率;Tr表示第r个模块的工作时间系数,等于该单元工作时间与全程工作时间之比。
根据以上步骤完成样本总量的分层分配后,需要对各层模块分别抽取故障模式组成故障模式集来。在目前的板级电路故障诊断中,测试性验证中采用简单随机抽样来抽取故障样本集。
简单随机抽样是在每层内开展抽样,直至抽取的故障模式数量等于被分配的故障模式数量。随机抽样往往导致传播型故障模式没有被抽取到。传播型故障是指一种故障模式的故障率低,但一旦发生则会传播扩散到其它元件,造成影响的规模较大。若对传播型故障的测试性设计不完备,一旦发生而没有正确检测和隔离,给使用方造成的风险将与常见故障风nr总和相提并论。因此在抽取故障模式时应当对传播型故障做重要抽样,以降低使用方风险。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于故障传播概率模型的故障模式优选方法,采用故障传播概率模型分析得到每个故障模式的故障扩散数量,基于故障扩散数量进行故障模式的优选,对传播型故障进行重点抽样,提高测试性验证方案可靠性。
为实现上述发明目的,本发明基于故障传播概率模型的故障模式优选方法,包括以下步骤:
S1:根据测试要求向系统中的各个模块分配故障模式抽取数量nr,nr表示系统中第r个模块分配的故障样本数量,r的取值范围为1≤r≤R,R表示系统中的模块数量;
S2:对系统中的每个模块分别选择故障模式,选择方法包括以下步骤:
S2.1:建立模块的故障传播概率模型,得到模块中每个故障模式fi和其一步传播后继故障模式fj之间的故障传播概率Pij,故障模式序号i,j的取值范围为1≤i,j≤Qr,Qr表示系统模块中的故障模式数量;
S2.2:计算故障模式i和与其一步传播后继故障模式j之间的故障扩散强度Iij:
其中,wi表示故障模式i的跨模块传播系数,wp表示预设的故障传播概率的权值,wd表示预设的节点出度的权值,dj表示故障j所在元件节点的出度,Fi表示故障模式i的一步传播后继故障模式集合;
对故障扩散强度进行归一化处理,归一化故障扩散强度I′ij的计算公式为:
S2.3:根据步骤S2.2得到的故障扩散强度计算元件间的故障扩散概率,计算公式的原则为:故障扩散强度越大,故障扩散概率越大,搜索每个故障模式i作为原发故障时的故障扩散概率最大路径,将该路径中包括的后继故障模式数量作为故障扩散数量Si;
S2.4:抽取故障扩散数量Si大于预设阈值δ的故障,组成集合φ,剩下的组成集合
S2.5:如果集合φ中的元素数量等于nr,集合φ作为选定故障模式集合Vr,本模块故障模式选择结束;
如果集合φ中的元素数量小于nr,从集合中抽取nr-|φ|个故障模式,|φ|表示集合φ中元素数量,这nr-|φ|个故障模式和集合φ中的故障模式共同构成选定故障模式集合Vr,本模块故障模式选择结束;
如果集合φ中的元素数量大于nr,从集合φ中抽取nr个故障模式,构成选定故障模式集合Vr,本模块故障模式选择结束;
S3:将系统中各个模块的选定故障模式集合Vr合并,得到系统的选定故障模式集合V=V1∪V2∪…∪VR。
本发明基于故障传播概率模型的故障模式优选方法,首先对系统各模块的故障样本量进行分配,然后对于各个模块建立模块的故障传播概率模型,得到各故障模式之间的故障传播概率,然后根据故障传播概率和故障模式所在元件节点的出度计算得到故障模式间的故障扩散强度,根据故障扩散强度搜索出每个故障模式作为原发故障时的故障扩散概率最大路径,将该路径中包括的后继故障模式数量作为故障扩散数量,筛选出故障扩散数量大于预定阈值的故障模式,优先进行选择;最后将各个模块的选定故障模式集合得到系统的选定故障模式集合。
本发明采用故障传播概率模型分析得到每个故障模式的故障扩散数量,基于故障扩散数量进行故障模式的优选,对传播型故障进行重点抽样,有效降低了因传播型故障没有被检验而带来的风险,从而提高测试性验证方案可靠性。
附图说明
图1是故障模糊状态和其三角隶属度函数示意图;
图2是本发明故障传播概率模型的故障模式优选方法的具体实施方式流程图;
图3是故障传播过程分析示意图;
图4是致冷组件的故障传播概率模型图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
为了更好地对本发明基于故障传播概率模型的故障模式优选方法的技术内容进行说明,首先对本发明采用的故障传播概率模型进行说明:
故障传播概率模型是在多信号流图模型的基础上提出来的,将系统当中的故障间的影响度和故障状态信息加入信号流图模型即可建立此模型。故障传播概率模型主要是用于推理故障模式的扩散,即当系统中发生了一个故障,利用此模型推理出故障会传播导致哪些器件也发生故障。
故障传播概率模型主要构成因素可以描述如下:
(1)系统模块中P个元件的集合M={m1,m2,…,mP},M表示系统模块中独立的具有相对完整功能的最小功能模块集合,其元素可以用来标识故障发生的位置,即对应的元件mp,1≤p≤P。
(2)模块的故障模式集F={f1,f2,…,fQ},表示系统模块中已知的具体故障模式,一个元件可能包含有多种故障模式。每个故障模式fi(1≤i≤Q)均有一个表示其所处状态,也就是严重程度的属性F_State(fi),采用模糊集合FS={正常,轻微,一般,严重,致命}来表征,模糊集FS的各个等级是通过一定的隶属函数来确定的。
故障状态F_State(fi)由模糊集合来表示,而其中具体的故障状态需要根据故障的信号或征兆来确定。该故障状态是一个模糊的概念,一般不能由一个确定的度量值来刻画,而通常处于某个参数的区间范围,如正常、严重等概念,在工程上一般都采用某个性能参数的模糊区间来进行描述。因此,本发明引入模糊集的理论来描述故障的状态,下面首先给出模糊集的定义。
对于论域U内任意一个元素x与一个集合A的关系只能是x∈A或两种。所以集合可通过特征函数来描述,并且每个集合都有一个特征函数,如:
函数称为集合A的特征函数,它刻画了对集合A的隶属情况,也称为集合A的隶属函数,而模糊集合就是根据隶属函数的概念定义出来的。
定义1:论域U={x}上的模糊子集A由隶属函数来表征,在闭区间[0,1]上取值,取值的大小反映了x对模糊集A的隶属程度。
对于本发明所针对的故障状态而言,论域U={x}就是对故障信号或征兆的度量值,对于模糊集合中的每一个状态都会有一个隶属度函数,故障状态的模糊集一般可用三角形隶属度函数进行恰当描述和表征。本发明采用五种故障状态FS={正常,轻微,一般,严重,致命}。图1是故障模糊状态和其三角隶属度函数示意图。
(3)连接各元件节点的有向边的集合E={Ep,p′},其中Ep,p′表示由元件p指向元件p′的有向边,其中1≤p′≤M,p′≠p。这里有向边有两层含义:一是表示模块间的物理连接或逻辑耦联;二是表示模块间故障的传播关系。
(4)由元件集M、故障模式集F和有向边集合E构成有向图。
基于以上原理,本发明提出了基于故障传播概率模型的故障模式优选方法。图2是本发明故障传播概率模型的故障模式优选方法的具体实施方式流程图。如图2所示,本发明基于故障传播概率模型的故障模式优选方法包括以下步骤:
S201:分配故障样本量:
首先根据测试要求向系统中的各个模块分配故障模式抽取数量nr,nr表示系统中第r个模块分配的故障样本数量,r的取值范围为1≤r≤R,R表示系统中的模块数量。
S202:令模块序号r=1;
S203:得到模块r中故障模式间的故障传播概率:
建立系统模块的故障传播概率模型,得到系统模块中每个故障模式fi和其一步传播后继故障模式fj之间的故障传播概率Pij,故障模式序号i,j的取值范围为1≤i,j≤Qr,Qr表示系统模块中的故障模式数量。
故障传播概率Pij表示故障模式i在前一节点发生后引起相邻后续节点产生故障模式j的概率大小。其概率的大小主要是由两个故障模式所在元件之间的物理链接关系决定。通过分析系统层次结构和物理连接机理以及结合故障领域专家分析系统和故障状态给出故障概率可以得出系统节点间故障传播概率的大小。系统原发故障的故障状态的强弱会对故障传播概率产生一定的影响,因此可以根据故障状态及各故障状态下对后继故障模式的影响度得到故障传播概率。例如系统中元件节点存在一个故障模式fi,由该故障模式诱发后一元件节点上发生故障模式fj的影响度分别为H=h1,h2,h3,h4,h5,hk表示故障模式fi处于第k个故障状态下对故障模式fj的影响度,本发明中故障模式有5种故障状态,因此也就有5个影响度,影响度是通过分析系统模块中各元件和故障的关系来得到的。那么故障模式fi到故障模式fj的故障传播概率Pij的大小可以表示为:
其中,uk(x)表示故障状态模糊集的隶属度函数,x表示各故障状态相应的征兆信号。
S204:计算模块r中故障模式间的故障扩散强度:
计算故障模式i和与其一步传播后继故障模式j之间的故障扩散强度Iij:
其中,wi表示故障模式i的跨模块传播系数,wp表示预设的故障传播概率的权值,wd表示预设的节点出度的权值,dj表示故障j所在元件节点的出度,Fi表示故障模式i的一步传播后继故障模式集合。
当电路系统中的某个元件节点发生故障时,由于元件节点之间的连接关系,此故障会向与此元件节点相连接的后续元件节点扩散,导致后续元件节点产生故障。而故障扩散的方向主要受两个方面的因素影响:(1)故障会选择故障传播概率大的方向的节点进行传播;(2)节点的长程连接以及节点的度数会对传播方向造成影响,因此在分析故障传播方向时需要综合考虑这两方面因素。所以本发明的故障扩散强度Iij是基于故障传播概率和故障所在元件节点的出度来计算的。故障扩散强度越大,表示故障越有可能在此边中发生,也最有可能传播导致后续故障的发生。
在系统模块中,故障之间的传播是通过元件节点的连接逐步传播的,因此本发明计算相邻故障模式间(即所在元件节点具有直接连接关系的故障模式)的故障扩散强度,用于进行故障传播分析。
为了便于后一步分析判断,对故障扩散强度进行归一化,归一化故障扩散强度I′ij的计算公式为:
通常为了方便计算机进行运算会采用邻接矩阵来表示故障传播概率模型。以Qr个元件的模块r为例,邻接矩阵Gr的大小为Qr×Qr,其中每个元素gij=I′ij。
S205:搜索故障扩散概率最大路径:
根据步骤S204得到的故障扩散强度计算元件间的故障扩散概率,计算公式的原则为:故障扩散强度越大,故障扩散概率越大,搜索每个故障模式i作为原发故障时的故障扩散概率最大路径,将该路径中包括的后继故障模式数量作为故障扩散数量Si。
在故障扩散过程中,会优先选择故障扩散强度大的路径进行传播,通过故障模式间的故障扩散强度,可以分析得到故障扩散最大概率传播路径,这条路径也是实际电路中最有可能发生故障传播的路径。图3是故障传播过程分析示意图。如图3所示,图中元件节点之间的有向边代表故障传播方向,有向边上所标识的数字为故障扩散强度。当元件节点m2上发生故障f2时,通过故障扩散强度可知故障最大概率传播路径为f2→f4→f5,即图3中的黑色粗线条所示,因此故障模式f2的后继故障模式有2个,即f4和f5,即故障扩散数量为2。
在大规模集成电路中,器件规模非常大以及电路链接错综复杂,用人工的方法进行故障原理推理显然是不合理,因此有必要借用计算机求出系统模块电路中故障最大概率传播路径。求解最佳路径的方法有很多种,迪杰斯特拉算法、福特算法、弗洛伊德算法、但茨希算法等。不同的算法有不同的故障扩散概率计算公式。其实就是相当于把故障扩散强度作为路径的权值,然后进行搜索。本实施例采用选用鲁棒性强、适合并行运算的蚁群算法来搜寻电路的故障传播最大概率路径。
蚁群算法是由意大利学者M.Dorigo于1991年在其博士论文中提出,并用此方法成功的解决了旅行商(TSP)问题。它是一种用来在图中寻找最优化路径的机率型算法,灵感来源于蚂蚁在寻找食物过程中发现最短路径的行为,其基本思想是,当蚂蚁沿着一条路到达终点以后会马上返回,这样路径短的蚂蚁来回时间短,在路途中洒下的信息素越多,就会有更多蚂蚁被吸引过来,形成正反馈从而产生更多的信息素,最后越来越多的蚂蚁聚集到最短路径上来,从而近似的找到了最短路径。
在系统模块电路中求解故障扩散概率最大的路径的数学模型如下:
其中,dk表示原发故障进行最大概率后继故障扩散所经过的N步路径中的第k步路径,qij k表示第k步路径的故障扩散概率,Fk-1表示第k-1步路径目的元件节点,Fk表示第k步路径目的元件节点。搜索的结束条件为整条路径的扩散概率小于10-8。
蚂蚁根据信息浓度寻找下一个路径,在t时刻,蚁群中的蚂蚁l从故障模式i选择后继故障模式作为行走路径的概率,也就是故障扩散概率ql ij(t)为:
其中,ηij(t)表示启发信息,本发明中ηij(t)=I′ij,α表示信息启发式因子,β表示期望启发式因子,分别表示控制信息与启发信息对蚂蚁选择路径的概率影响参数,Nl i是故障模式i的一步后继故障模式集合,τij(t)表示信息激素,其计算公式为:
τij(t)=(1-ρ)τij(t-1)+Δτij
其中,ρ表示信息激素的挥发程度,τij(t)的初始值为τij(0)根据需要设置,Δτij表示信息素总增量,其中Na表示蚂蚁的数量,表示蚂蚁l的信息素。利用蚁群算法,就可以搜索每个故障模式i作为原发故障时的故障扩散概率最大路径。蚁群算法的详细求解过程可以参见文献“李士勇.蚁群算法及其应用[M]哈尔滨:哈尔滨工业大学出版社,2004.”。
S206:故障模式分类:
对于Q个故障模式fi,将其中故障扩散数量Si大于预设阈值δ的故障模式抽取出来,组成集合φ,剩下的组成集合阈值δ是根据实际情况设置的。
S207:判断集合φ中的元素数量|φ|是否等于nr,如果是,进入步骤S208,否则进入步骤S209。
S208:将集合φ作为选定故障模式集合Vr,即Vr=φ,本模块故障模式选择结束,进入步骤S212。
S209:判断集合φ中的元素数量是否小于nr,如果是,进入步骤S210,否则进入步骤S211。
S210:从集合中抽取nr-|φ|个故障模式,其集合记为|φ|表示集合φ中元素数量,这nr-|φ|个故障模式和集合φ中的故障模式共同构成选定故障模式集合Vr,即本模块故障模式选择结束,进入步骤S212。
本步骤中的故障抽取可以直接采用随机抽取,也可以根据其他参数进行抽取,根据实际情况而定。
S211:从集合φ中抽取nr个故障模式,其集合记为φ*,构成选定故障模式集合Vr,即Vr=φ*,本模块故障模式选择结束,进入步骤S212。
同样地,本步骤中的故障抽取可以直接采用随机抽取,也可以根据其他参数进行抽取,根据实际情况而定。
S212:判断是否r<R,如果是,进入步骤S213,否则进入步骤S214。
S213:令r=r+1,返回步骤S203。
S214:将系统中各个模块的选定故障模式集合Vr合并,得到系统的选定故障模式集合V=V1∪V2∪…∪VR。
为了说明本发明的技术效果,采用一个某型空空导弹进行实验验证。该空空导弹由致冷组件、防振器、整流器、剪切支架、锁制器、电路盒6个模块组成,由于模块较多,各模块的处理方法相类似,所以主要对致冷组件建立故障概率传播信号流图模型,进行故障推理,完成测试性验证方案实施环节。表1是致冷组件中各元件的故障模式信息。
所属部件 | 模式编号 | 故障模式 |
M1 | f1 | 密封圈A失效 |
M2 | f2 | 单向阀失效 |
M3 | f3 | 导管零部件失效 |
M4 | f4 | 排气嘴失效 |
M5 | f5 | 密封圈B失效 |
M6 | f6 | 转接头材料缺陷 |
M7 | f7 | 电磁开关部件失效 |
M8 | f8 | 充气阀漏气 |
M9 | f9 | 致冷导管护口变形 |
M10 | f10 | 压力表失效 |
M11 | f11 | 电磁开关漏气 |
M12 | f12 | 致冷组件失效 |
表1
根据合同约定的极限质量水平q1、设计值q0、承制方风险(α)、使用方风险(β),由抽样方案(n,C)与承制方风险(α)、使用方风险(β)的关系建模公式得到抽样方案为(50,5)。验证流程需要把总的样本量按比例分层分配到系统中的6个模块中去,已知致冷组件在系统中所占的权重w1=0.121,可知故障样本量分配到致冷组件的量为n×w1≈6。根据表1可知致冷组件的故障模式共计12个,因此需要在12个故障模式中抽取6个故障模式组成故障模式集。
根据系统的电路链接方式,可以建立致冷组件的故障概率传播信号流图模型。图4是致冷组件的故障概率传播信号流图。如图4所示,致冷组件的元件集为M={m1,m2,m3,m4,m5,m6,m7,m8,m9,m10,m11,m12},故障模式集为F={f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12},每个元件只有一个故障模式。
首先计算得到各故障间的故障传播概率,然后根据故障扩散强度的运算公式,计算得到致冷组件各个相邻故障模式间(即各个有向边之间)的故障扩散强度:
I1,8=0.71,I2,8=0.29,I3,9=1.0,I4,10=0.50,I5,10=0.25,I6,10=0.167,I7,11=0 I8,12=1.0,I9,12=0.226,I10,12=0,I11,12=0.323,
归一化后的故障扩散强度为:
I′1,8=0.16,I′2,8=0.06,I′3,9=0.22,I′4,10=0.11,I′5,10=0.05,I′6,10=0.03,I′7,11=0 I′8,12=0.22,I′9,12=0.05,I′10,12=0,I′11,12=0.07,
以故障f1为例进行故障传播推理,采用蚁群算法得到故障扩散概率最大路径,其路径为f1→f8→f12。可见,故障f1的后继故障模式有2个,即故障扩散数量为2。用同样的方法推理其它模块求出求最大故障概率传播路径,根据路径求出故障扩散数量。然后筛选出故障扩散数量大于1的故障模式,优先进行选择。表2是本发明与传统简单随机抽样的故障模式抽取结果。
序号 | 故障模式 | 故障扩散强度 | 本发明 | 简单随机抽样 |
M1 | 密封圈A失效 | 2 | √ | √ |
M2 | 单向阀失效 | 2 | √ | |
M3 | 导管零部件失效 | 2 | √ | √ |
M4 | 排气嘴失效 | 1 | ||
M5 | 密封圈B失效 | √ | ||
M6 | 转接头材料缺陷 | 2 | √ | √ |
M7 | 电磁开关部件失效 | |||
M8 | 充气阀漏气 | 1 | √ | √ |
M9 | 致冷导管护口变形 | 1 | √ | |
M10 | 压力表失效 | 1 | ||
M11 | 电磁开关漏气 | 1 | √ | |
M12 | 致冷组件失效 | 1 |
表2
从表2可以看出,传统的简单随机抽样是随机抽取6个故障模式,并没有考虑传播型故障模式。而本发明先筛选故障扩散强度大于1的故障模式,优先进行选择,从表2可知故障扩散强度大于1的故障模式有4个,因此将这4个故障模式作为选择的故障模式,然后在剩下的8个故障模式中随机选取2个故障模式,由这6个故障模式构成选定故障模式。通过比较两种方法可知,简单随机抽样经常会导致传播型故障模式没有被抽取到,而本发明利用故障传播概率模型推理故障传播的后继模式集,根据后继模式集的个数提出基于故障传播概率模型的优化抽样方法。此方法对传播型故障进行重点抽样,有效降低了因传播型故障没有被检验给使用方带来的风险。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (4)
1.一种基于故障传播概率模型的故障模式优选方法,其特征在于,包括以下步骤:
S1:根据测试要求向系统中的各个模块分配故障模式抽取数量nr,nr表示系统中第r个模块分配的故障样本数量,r的取值范围为1≤r≤R,R表示系统中的模块数量;
S2:对系统中的每个模块分别选择故障模式,选择方法包括以下步骤:
S2.1:建立模块的故障传播概率模型,得到模块中每个故障模式fi和其一步传播后继故障模式fj之间的故障传播概率Pij,故障模式序号i,j的取值范围为1≤i,j≤Qr,Qr表示系统模块中的故障模式数量;
S2.2:计算故障模式i和与其一步传播后继故障模式j之间的故障扩散强度Iij:
其中,wi表示故障模式i的跨模块传播系数,wp表示预设的故障传播概率的权值,wd表示预设的节点出度的权值,dj表示故障j所在元件节点的出度,Fi表示故障模式i的一步传播后继故障模式集合;
对故障扩散强度进行归一化处理,归一化故障扩散强度I′ij的计算公式为:
S2.3:根据步骤S2.2得到的故障扩散强度计算元件间的故障扩散概率,计算公式的原则为:故障扩散强度越大,故障扩散概率越大,搜索每个故障模式i作为原发故障时的故障扩散概率最大路径,将该路径中包括的后继故障模式数量作为故障扩散数量Si;
S2.4:抽取故障扩散数量Si大于预设阈值δ的故障,组成集合φ,剩下的组成集合
S2.5:如果集合φ中的元素数量等于nr,集合φ作为选定故障模式集合Vr,本模块故障模式选择结束;
如果集合φ中的元素数量小于nr,从集合中抽取nr-|φ|个故障模式,|φ|表示集合φ中元素数量,这n-|φ|个故障模式和集合φ中的故障模式共同构成选定故障模式集合Vr,本模拟故障模式选择结束;
如果集合φ中的元素数量大于nr,从集合φ中抽取n个故障模式,构成选定故障模式集合Vr,本模拟故障模式选择结束;
S3:将系统中各个模块的选定故障模式集合Vr合并,得到系统的选定故障模式集合V=V1∪V2∪…∪VR。
2.根据权利要求1所述的故障优选模式优选方法,其特征在于,所述步骤S1中故障传播概率Pij根据以下公式得到:
其中,uk(x)表示故障状态模糊集的隶属度函数,x表示各故障状态相应的征兆信号,故障状态模糊集为FS={正常,轻微,一般,严重,致命},hk表示故障模式fi处于第k个故障状态下对故障模式fj的影响度。
3.根据权利要求1所述的故障模式优选方法,其特征在于,所述步骤S3中故障扩散概率最大路径采用蚁群算法搜索得到,故障扩散强度作为启发信息来计算故障扩散概率。
4.根据权利要求1所述的故障模式优选方法,其特征在于,所述步骤S2.5中的故障模式的抽取为随机抽取。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510280067.XA CN104898039B (zh) | 2015-05-27 | 2015-05-27 | 基于故障传播概率模型的故障模式优选方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510280067.XA CN104898039B (zh) | 2015-05-27 | 2015-05-27 | 基于故障传播概率模型的故障模式优选方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104898039A true CN104898039A (zh) | 2015-09-09 |
CN104898039B CN104898039B (zh) | 2017-09-12 |
Family
ID=54030811
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510280067.XA Expired - Fee Related CN104898039B (zh) | 2015-05-27 | 2015-05-27 | 基于故障传播概率模型的故障模式优选方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104898039B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107192915A (zh) * | 2017-05-02 | 2017-09-22 | 国家电网公司 | 一种基于图论法的变压器故障诊断方法 |
CN107490758A (zh) * | 2017-07-14 | 2017-12-19 | 电子科技大学 | 基于故障传播的模块化bp神经网络电路故障诊断方法 |
CN112100456A (zh) * | 2020-09-16 | 2020-12-18 | 广东电网有限责任公司电力科学研究院 | 一次设备共性缺陷或故障的判断方法、装置及终端设备 |
CN114063588A (zh) * | 2020-07-29 | 2022-02-18 | 中车株洲电力机车研究所有限公司 | 一种传动控制单元测试速度的选取方法、装置及设备 |
CN114205289A (zh) * | 2021-10-22 | 2022-03-18 | 清华大学 | 一种分段路由网络中考虑故障的流量工程计算方法和装置 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090083576A1 (en) * | 2007-09-20 | 2009-03-26 | Olga Alexandrovna Vlassova | Fault tree map generation |
EP2110778A2 (en) * | 2008-04-14 | 2009-10-21 | Honeywell International Inc. | Determining corrective actions using a geometrically-based determination of sufficient confidence |
CN101718834A (zh) * | 2009-12-01 | 2010-06-02 | 湖南大学 | 一种模拟电路故障传播特性的分析方法 |
CN102662389A (zh) * | 2012-03-21 | 2012-09-12 | 中国人民解放军国防科学技术大学 | 一种基于故障特性分析的监测参数选择方法 |
CN103064340A (zh) * | 2011-10-21 | 2013-04-24 | 沈阳高精数控技术有限公司 | 一种面向数控机床的故障预测方法 |
CN103298009A (zh) * | 2012-12-31 | 2013-09-11 | 南京理工大学 | 复杂无线网络控制系统的建模及故障传播方法 |
CN103399979A (zh) * | 2013-07-04 | 2013-11-20 | 电子科技大学 | 一种板级电路测试性模型自动生成方法 |
CN103870659A (zh) * | 2014-03-28 | 2014-06-18 | 吉林大学 | 一种数控机床故障分析方法 |
JP2015005077A (ja) * | 2013-06-20 | 2015-01-08 | 株式会社日立製作所 | 障害情報管理方法、障害情報管理装置及びプログラム |
-
2015
- 2015-05-27 CN CN201510280067.XA patent/CN104898039B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090083576A1 (en) * | 2007-09-20 | 2009-03-26 | Olga Alexandrovna Vlassova | Fault tree map generation |
EP2110778A2 (en) * | 2008-04-14 | 2009-10-21 | Honeywell International Inc. | Determining corrective actions using a geometrically-based determination of sufficient confidence |
CN101718834A (zh) * | 2009-12-01 | 2010-06-02 | 湖南大学 | 一种模拟电路故障传播特性的分析方法 |
CN103064340A (zh) * | 2011-10-21 | 2013-04-24 | 沈阳高精数控技术有限公司 | 一种面向数控机床的故障预测方法 |
CN102662389A (zh) * | 2012-03-21 | 2012-09-12 | 中国人民解放军国防科学技术大学 | 一种基于故障特性分析的监测参数选择方法 |
CN103298009A (zh) * | 2012-12-31 | 2013-09-11 | 南京理工大学 | 复杂无线网络控制系统的建模及故障传播方法 |
JP2015005077A (ja) * | 2013-06-20 | 2015-01-08 | 株式会社日立製作所 | 障害情報管理方法、障害情報管理装置及びプログラム |
CN103399979A (zh) * | 2013-07-04 | 2013-11-20 | 电子科技大学 | 一种板级电路测试性模型自动生成方法 |
CN103870659A (zh) * | 2014-03-28 | 2014-06-18 | 吉林大学 | 一种数控机床故障分析方法 |
Non-Patent Citations (2)
Title |
---|
杨成林 等: "改进的故障分辨率计算公式和序替代方法", 《电子测量与仪器学报》 * |
高鑫宇 等: "基于模糊概率多信号流图的故障传播模型研究", 《测试技术学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107192915A (zh) * | 2017-05-02 | 2017-09-22 | 国家电网公司 | 一种基于图论法的变压器故障诊断方法 |
CN107490758A (zh) * | 2017-07-14 | 2017-12-19 | 电子科技大学 | 基于故障传播的模块化bp神经网络电路故障诊断方法 |
CN107490758B (zh) * | 2017-07-14 | 2019-08-13 | 电子科技大学 | 基于故障传播的模块化bp神经网络电路故障诊断方法 |
CN114063588A (zh) * | 2020-07-29 | 2022-02-18 | 中车株洲电力机车研究所有限公司 | 一种传动控制单元测试速度的选取方法、装置及设备 |
CN114063588B (zh) * | 2020-07-29 | 2023-10-31 | 中车株洲电力机车研究所有限公司 | 一种传动控制单元测试速度的选取方法、装置及设备 |
CN112100456A (zh) * | 2020-09-16 | 2020-12-18 | 广东电网有限责任公司电力科学研究院 | 一次设备共性缺陷或故障的判断方法、装置及终端设备 |
CN114205289A (zh) * | 2021-10-22 | 2022-03-18 | 清华大学 | 一种分段路由网络中考虑故障的流量工程计算方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104898039B (zh) | 2017-09-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104898039A (zh) | 基于故障传播概率模型的故障模式优选方法 | |
CN103870659B (zh) | 一种数控机床故障分析方法 | |
CN106341414A (zh) | 一种基于贝叶斯网络的多步攻击安全态势评估方法 | |
CN106021062B (zh) | 关联故障的预测方法和系统 | |
Niu et al. | A new efficient algorithm for finding all d-minimal cuts in multi-state networks | |
CN104298843A (zh) | 一种基于着色随机Petri网的复杂机构动态级联可靠性建模方法 | |
CN105868116A (zh) | 基于语义变异算子的测试用例生成和优化方法 | |
CN104125112B (zh) | 基于物理‑信息模糊推理的智能电网攻击检测方法 | |
CN108732528A (zh) | 一种基于深度置信网络的数字化电能表故障诊断方法 | |
CN104899327A (zh) | 一种无类别标签的时间序列异常检测方法 | |
CN103840967A (zh) | 一种电力通信网中故障定位的方法 | |
CN104466959A (zh) | 电力系统关键线路辨识方法和系统 | |
CN103778467A (zh) | 一种选择电力系统暂态稳定评估输入特征量的方法 | |
CN104361169A (zh) | 一种基于分解法建模的可靠性监测方法 | |
CN105117772A (zh) | 一种多状态系统可靠性模型的参数估计方法 | |
CN104090227A (zh) | 一种模拟电路故障诊断中的测点选择方法 | |
CN105159827A (zh) | 一种面向gui软件的可靠性加速测试方法 | |
CN104539601A (zh) | 动态网络攻击过程可靠性分析方法及系统 | |
CN105068928A (zh) | 一种基于复杂网络理论的软件测试用例生成方法 | |
CN104486222B (zh) | 基于蚁群优化算法的小时延缺陷测试关键路径选择方法 | |
CN103065200A (zh) | 动态社区发现及追踪系统及其方法 | |
CN104597375A (zh) | 一种电网故障诊断方法 | |
CN105786763A (zh) | 一种设备集成系统网络的故障传播路径的生成方法 | |
CN103400213A (zh) | 一种基于lda与pca的骨干网架生存性评估方法 | |
CN108932197A (zh) | 基于参数Bootstrap重抽样的软件失效时间预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20170912 Termination date: 20200527 |