一种航空器系统定量相似性分析方法、装置及介质

1.本发明涉及相似性分析技术领域，具体的说，涉及一种航空器系统定量相似性分析方法、装置及介质。


背景技术：

2.空客a330、a380、a350、a320neo飞机属于新一代民用航空器，它们的维修大纲、重要维修项目(mrbr)、维修时间间隔及主最低设备清单(mmel)的制定都是依据前一代成熟的相似机型方案借鉴而来。相同对于波音b777、b787、b787-8/9飞机也是新一代依据老机型的大纲清单等进行相似性参考。目前对于mmel的制定主要由飞机制造厂商制定，飞机制造商在民航局foeb审批以前提交的为主最低设备清单建议稿(pmmel)，待中国民用航空局审批后以中国民用航空局的名义发布为mmel，作为航空公司制定相关机型最低设备清单(mel)的主要依据文件，在不影响飞机安全的基础下，通过高质量的mel手册提高飞机的签派率和经济性。
3.通过分析pmmel制定需求来源的三个方面，即用户需求分析、相似机型对比及安全性分析。理清用户需求内容，建立规范的用户需求数据标准，确定历史用户案例库相似pmmel项目提取方法，相似机型的相似pmmel项目提取方法以及安全分析备选pmmel项目提取方法，验证目标机型pmmel项目清单对局方要求的满足性，形成一套落实pmmel利益攸关方的需求，具有继承性、系统性和通用性的pmmel需求工程分析方法和分析程序，并开发pmmel需求分析和编制原型系统，为国产民机型号pmmel制定提供技术支撑。
4.现有技术中,对系统相似判断方法现有三种方法。
5.方法一：首先分析相同系统间相似机理的内涵及其数学表达，根据数学表达进行相似对比。研究分析相似系统的一般性，找到存在的主要因素，再对相似系统进行并列、细化分解，发现相似系统之间的相似致因。将系统并列，对比分析确定相似元及其相似特征值，然后进行相似度计算。相似度的大小反映系统之间的相似程度，是对应相似元的数量及每个相似元对相似度影响权重的函数。权重w的客观性也影响着相似度的准确性。系统的相似度作为整体性判断的度量指标，可以通过各子系统、各相似特征的相似度，为从系统内部剖析对比提供可能。
6.方法二：从系统论的角度出发，对系统的“功能-结构”进行相似分析。将系统分为第k层次、第m层次、第m-1层次。第k层次是以子系统进行划分，第m层次直接划分为元素也成为“部件”层次，第m-1层次由多个元素通过一定的作用关系连接而成并且此时系统不可再分。分别以这三个层次通过数理方程表达式、结构功能图对功能相似和结构相似及两者间的固有关系相似进行分析论述，总结出系统相似性。
7.方法三：对于研究相似系统失效的相似相关性，从系统的可靠性和重要性进行分析对比，通过系统组件可靠性和重要性的相似来判断整体相似。对系统发生故障的数据进行统计，对数据处理后基于无模型抽样进行可靠性分析，建立可靠性框图并使用无模型抽样算法扩容各部件目标样本，做出失效曲线后看是否重合度较高。然后建立故障树使用sdp
方法进行不确定性变量的重要度分析，确定重要度指标后再将两组数据，利用蒙特卡罗法和概率密度演化法得到图形的分布状况，对比后得出系统相似相关性。
8.现有的飞机系统相似性判断主要依靠工程师经验和主观定性判断，缺乏严谨的分析流程和数学模型计算。现有的解决方案的很难评估高度复杂的机电系统，相似性判断主要停留在系统组成架构和部件件号，无法对系统功能及内在控制信号流进行评估。若飞机系统控制交联，功能高度耦合的相似性评估结果的可信度就非常低，工程应用安全风险较大。


技术实现要素：

9.本发明针对现有技术存在的不足，提供了一种航空器系统定量相似性分析方法、装置及介质，首先通过对相似机型的结构进行分析，选取相同类似的系统结构建立故障树后进行结构重要度计算。确定底事件发生的概率后，再进行定量分析，从故障发生率数据入手进行概率重要度及关键重要度计算。对得到的重要度数据运用数理统计方法构建三个参数的重要矩阵，利用余弦相似度函数算法，输入相似系统的矩阵值后输出余弦值来衡量特征之间的接近程度，得到是否相似相关的结论。
10.本发明的具体技术方案如下：
11.根据本发明的第一技术方案，提供了一种航空器系统定量相似性分析方法，包括以下步骤：基于两种不同型号的航空器的功能系统，分别建立故障树；通过如下方法分别计算两种不同型号的航空器的功能系统的结构重要度、概率重要度以及关键重要度：通过分析故障树的结构上各个事件的重要程度计算结构重要度；通过计算故障树的顶事件的发生概率对底事件发生概率的偏导数，得到概率重要度；通过计算故障树的底事件故障概率的变化率及其引起的顶事件发生概率的变化率之比，得到关键重要度；对所述两种不同型号的航空器的功能系统的概率重要度和关键重要度分别进行余弦相似度计算，得到两个相似度值；基于所述两个相似度值，来确定两种不同型号的航空器的相似度。
12.作为优选的技术方案，所述通过分析故障树的结构上各个事件的重要程度计算结构重要度，包括：
13.通过如下公式计算所述结构重要度：
[0014][0015]
式中，式中，为系统的结构重要度，仅取决于第i个部件在系统中所处的位置，与顶事件的发生概率无关，[φ(1i,x)-φ(0i,x)]为系统中第i个部件由正常状态(0)转为故障状态(1)，其他部件状态不变时，系统的结构函数发生变化。
[0016]
作为优选的技术方案，所述通过计算故障树的顶事件的发生概率对底事件发生概率的偏导数，得到概率重要度，包括：
[0017]
所述概率重要度通过如下公式计算：
[0018][0019]
式中，i
p(i)
为第i个基本事件的概率重要度，p(xi)为第i个事件发生的概率，p(t)为某个事件在t时刻发生的概率。
[0020]
作为优选的技术方案，所述通过计算故障树的底事件故障概率的变化率与它引起的顶事件发生概率的变化率之比，得到关键重要度，包括：
[0021]
所述关键重要度通过如下公式计算：
[0022][0023]
其中，i
c(i)
为第i个基本事件的关键重要度，p(xi)为第i个事件发生的概率，p(t)为某个事件在t时刻发生的概率。
[0024]
作为优选的技术方案，对所述两种不同型号航空器的功能系统的概率重要度和关键重要度分别进行余弦相似度计算，得到两个相似度值，包括：
[0025]
将所述两种不同型号的航空器的功能系统的概率重要度分别转换成两个向量，或者将所述两种不同型号的航空器的功能系统的关键重要度分别转换成两个向量；
[0026]
根据如下公式计算得到相似度值：
[0027][0028]
其中，cos(θ)表示相似度值，向量a表示其中一种型号的航空器的功能系统的概率重要度或关键重要度，向量b表示另一种型号的航空器的功能系统的概率重要度或关键重要度。
[0029]
作为优选的技术方案，基于所述两个相似度值，来确定两种不同型号的航空器的相似度，包括：
[0030]
所述相似度是余弦值或弧度或角度；当余弦值越接近1，表明两种不同型号的航空器的相似度越高；当弧度和角度的值越小，表明两种不同型号的航空器的相似度越高。
[0031]
根据本发明的第二技术方案，提供了一种航空器系统定量相似性分析装置，包括处理器，所述处理器配置为：基于两种不同型号的航空器的功能系统，分别建立故障树；通过如下方法分别计算两种不同型号的航空器的功能系统的结构重要度、概率重要度以及关键重要度：通过分析故障树的结构上各个事件的重要程度计算结构重要度；通过计算故障树的顶事件的发生概率对底事件发生概率的偏导数，得到概率重要度；通过计算故障树的底事件故障概率的变化率及其引起的顶事件发生概率的变化率之比，得到关键重要度；对所述两种不同型号的航空器的功能系统的概率重要度和关键重要度分别进行余弦相似度计算，得到两个相似度值；基于所述两个相似度值，来确定两种不同型号的航空器的相似度。
[0032]
作为优选的技术方案，所述处理器还配置为：
[0033]
将所述两种不同型号的航空器的功能系统的概率重要度分别转换成两个向量，或者将所述两种不同型号的航空器的功能系统的关键重要度分别转换成两个向量；
[0034]
根据如下公式计算得到相似度值：
[0035][0036]
其中，cos(θ)表示相似度值，向量a表示其中一种型号的航空器的功能系统的概率
重要度或关键重要度，向量b表示另一种型号的航空器的功能系统的概率重要度或关键重要度。
[0037]
作为优选的技术方案，所述处理器还配置为：
[0038]
所述相似度是余弦值或弧度或角度；当余弦值越接近1，表明两种不同型号的航空器的相似度越高；当弧度和角度的值越小，表明两种不同型号的航空器的相似度越高。
[0039]
作为优选的技术方案，所述处理器还配置为：
[0040]
通过如下公式计算所述结构重要度：
[0041][0042]
式中，式中，为系统的结构重要度，仅取决于第i个部件在系统中所处的位置，与顶事件的发生概率无关，[φ(1i,x)-φ(0i,x)]为系统中第i个部件由正常状态(0)转为故障状态(1)，其他部件状态不变时，系统的结构函数发生变化。
[0043]
作为优选的技术方案，所述处理器还配置为：所述概率重要度通过如下公式计算：
[0044][0045]
式中，i
p(i)
为第i个基本事件的概率重要度，p(xi)为第i个事件发生的概率，p(t)为某个事件在t时刻发生的概率。
[0046]
作为优选的技术方案，所述处理器还配置为：所述关键重要度通过如下公式计算：
[0047][0048]
其中，i
c(i)
为第i个基本事件的关键重要度，p(xi)为第i个事件发生的概率，p(t)为某个事件在t时刻发生的概率。
[0049]
根据本发明的第三技术方案，提供了一种计算机可读存储介质，其上存储有计算机可读指令，当所述计算机可读指令被计算机的处理器执行时，使计算机执行根据本发明实施例中的任一所述的方法。
[0050]
有益效果在于：
[0051]
根据本发明实施例的航空器系统定量相似性分析方法、装置及介质，其对国产大飞机(如c919型号飞机)pmmel手册制定提供一种全新的想法，通过对机型相似性分析方法的研究，为国产大飞机运行技术手册的制定提供科学的技术方法。pmmel手册项目清单能解决因为维修时间过长而导致航班运行延误的状况，合理运用可有效提高飞机的利用率和航班正点率，降低了运营成本，完整的持续适航文件体系是国产大飞机最终取得国际适航证的必要基础。
附图说明
[0052]
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍。在所有附图中，类似的元件或部分一般由类似的附图标记标识。附图中，各元件或部分并不一定按照实际的比例绘制。
[0053]
图1是液压系统组成图。
[0054]
图2是液压作动系统组成图。
[0055]
图3是正常刹车系统液压故障树。
[0056]
图4是襟翼液压系统故障树。
[0057]
图5是副翼液压系统故障树。
[0058]
图6是方向舵液压系统故障树。
[0059]
图7是扰流板(1、5)液压系统故障树。
[0060]
图8是一号发动机反推液压系统故障树。
[0061]
图9是正常刹车液压系统故障树。
[0062]
图10是襟翼液压系统故障树。
[0063]
图11是副翼液压系统故障树。
[0064]
图12是方向舵液压系统故障树。
[0065]
图13是右发动机反推液压系统故障树。
[0066]
图14是备用刹车液压系统故障树。
[0067]
图15是方向舵液压系统故障树。
[0068]
图16是扰流板(2、4)液压系统故障树。
[0069]
图17是二号发动机反推液压系统故障树。
[0070]
图18是前轮转弯液压系统故障树。
[0071]
图19是刹车液压系统故障树。
[0072]
图20是扰流板液压系统故障树。
[0073]
图21是左发动机反推液压系统故障树。
[0074]
图22是扰流板液压系统故障树。
[0075]
图23是襟缝翼液压系统故障树。
[0076]
图24是副翼液压系统故障树。
[0077]
图25是方向舵液压系统故障树。
[0078]
图26是扰流板(3)液压系统故障树。
[0079]
图27是襟缝翼液压系统故障树。
[0080]
图28是备用方向舵液压系统故障树。
[0081]
图29是发动机反推液压系统故障树。
具体实施方式
[0082]
下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0083]
现在结合说明书附图对本发明做进一步的说明。
[0084]
本发明实施例提供了一种航空器系统定量相似性分析方法。该方法基于两种不同型号的航空器的功能系统，分别建立故障树；通过如下方法分别计算两种不同型号的航空器的功能系统的结构重要度、概率重要度以及关键重要度：通过分析故障树的结构上各个事件的重要程度计算结构重要度；通过计算故障树的顶事件的发生概率对底事件发生概率的偏导数，得到概率重要度；通过计算故障树的底事件故障概率的变化率及其引起的顶事
件发生概率的变化率之比，得到关键重要度；对所述两种不同型号的航空器的功能系统的概率重要度和关键重要度分别进行余弦相似度计算，得到两个相似度值；基于所述两个相似度值，来确定两种不同型号的航空器的相似度。
[0085]
其中，故障树分析是通过将单调联系的故障树转化为只含与、或门和底事件的故障树，按故障树的逻辑结构根据结构函数向上进行运算，从而求出顶事件的发生概率，即顶事件的不可靠度
[19]
。设底事件为x1，x2，...，xn，则底事件发生故障的概率为f1(t)，f2(t)，...，fn(t)。
[0086]
逻辑门构造函数如下:
[0087]
与门结构函数表达式：
[0088][0089]
或门结构函数表达式：
[0090][0091]
下面本发明实施例将详细阐述各重要度的计算方法及其意义。
[0092]
结构重要度是从故障树的结构上分析各个事件的重要程度，是故障树定性分析的一部分。可以根据结构重要度系数大小精确计算排列各基本事件的重要顺序，或者根据最小割集(径集)近似判断结构重要度顺序
[
。
[0093][0094]
式中：式中：为系统的结构重要度，仅取决于第i个部件在系统中所处的位置，与顶事件的发生概率无关。[φ(1i,x)-φ(0i,x)]为系统中第i个部件由正常状态(0)转为故障状态(1)，其他部件状态不变时，系统的结构函数发生变化。
[0095]
概率重要度是表示底事件发生概率变化会给顶事件发生概率以多大的影响，用顶事件的发生概率对底事件发生概率的偏导数来表示
[21]
。概率重要度计算公式：
[0096][0097]
式中i
p(i)
为第i个基本事件的概率重要度，p(xi)为第i个事件发生的概率。
[0098]
关键重要度指的是底事件故障概率的变化率与它引起的顶事件发生概率的变化率之比
[21]
。用于判断底事件发生概率对于顶事件的关键程度，表达式为：
[0099][0100]
其中i
c(i)
为第i个基本事件的关键重要度。
[0101]
随着设备性能不断提升，系统结构复杂程度不断增加其所面临的寿命与可靠性、风险、可用性以及维修故障等问题日益突出。可靠性作为保障重大装备自主创新能力和核心竞争力的共性关键技术与基础性问题，已经引起了深度关注。重要性测度是近年来可靠性领域发展起来的新兴研究方向。各种重要性测度之间，以及它们与可靠性、风险及优化的紧密联系是一个不可缺失的环节，特别是在复杂系统的可靠性与风险评估方面，对建立有效的定性、定量分析方法有重大的指导意义和参考价值。
[0102]
在直角坐标系中，向量表示三角形的余弦函数，向量a(x1,y1)，向量b(x2,y2)。向量a和向量b在直角坐标中的长度为向量a和b的距离用向量c表示，那么向量c在直角坐标系中的长度为将a,b,c代入三角函数公式后得到二维空间中余弦函数的表达式：
[0103][0104]
那么多维空间余弦函数的表达式为：
[0105][0106]
余弦相似度用于计算个体间的相似度，距离越大相似度越小，相反如果距离越小那么相似度就越大。matlab中余弦相似度算法的原理是通过一个向量空间中两个向量夹角间的余弦值作为衡量两个个体或系统之间差异的大小，通过编程对向量间角度的计算得到余弦值。并且余弦相似度在多维空间的计算较也为方便，可以把矩阵转化成向量的形式进行计算，用于判断两个矩阵的相似相关性。当余弦值接近1时夹角趋于0，表明两个向量越相似；余弦值如果接近于0，夹角趋于90度，则表明两个向量相似性越弱。同时还可以得出两个向量之间的弧度和角度，当弧度和角度的值越小，表明两个向量越相似。
[0107]
在一些实施例中，也可以对得到的重要度数据运用数理统计方法做成表形成矩阵，利用余弦相似度函数算法，输入相似系统的矩阵值后输出余弦值来衡量特征之间的接近程度，得到是否相似相关的结论。
[0108]
本发明实施例还提供了一种航空器系统定量相似性分析装置，包括处理器，所述处理器配置为：基于两种不同型号的航空器的功能系统，分别建立故障树；通过如下方法分别计算两种不同型号的航空器的功能系统的结构重要度、概率重要度以及关键重要度：通过分析故障树的结构上各个事件的重要程度计算结构重要度；通过计算故障树的顶事件的发生概率对底事件发生概率的偏导数，得到概率重要度；通过计算故障树的底事件故障概率的变化率及其引起的顶事件发生概率的变化率之比，得到关键重要度；对所述两种不同型号的航空器的功能系统的概率重要度和关键重要度分别进行余弦相似度计算，得到两个相似度值；基于所述两个相似度值，来确定两种不同型号的航空器的相似度。
[0109]
需要注意的是，处理器可以是包括一个以上通用处理设备的处理设备，诸如微处理器、中央处理单元(cpu)、图形处理单元(gpu)等。更具体地，处理器可以是复杂指令集计算(cisc)微处理器、精简指令集计算(risc)微处理器、超长指令字(vliw)微处理器、运行其他指令集的处理器或运行指令集的组合的处理器。处理器还可以是一个以上专用处理设备，诸如专用集成电路(asic)、现场可编程门阵列(fpga)、数字信号处理器(dsp)、片上系统(soc)等。
[0110]
本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机可读指令，当所述计算机可读指令被计算机的处理器执行时，使计算机本发明实施例中的任一所述的方法。
[0111]
下面本发明实施例将结合具体的实施案例来阐述本发明实施例的原理及其有益
效果。
[0112]
案例一：相似机型液压系统分析
[0113]
液压系统是指从发动机、电机、apu、epu、rat等获得电能后转化为液压能，在对液压能进行分配和调节，传输到飞机各类液压作动装置上的系统
[8]
。现代飞机多采用两套(或多套)相互独立的液压系统，一套称为主系统，主要利用液压进行作动转换或伺服操纵起落架、扰流板、襟翼等的收放，并向液压操纵系统助力器作动筒供压；另一套称为助力系统，主要用于操纵飞机的舵面。另外还设了一套应急系统，应急系统的液压由单独的电瓶驱动电动液压泵或靠冲压空气涡轮驱动的液压泵提供。
[0114]
液压系统对保证飞行安全，完成设计性能正常飞行起到了至关重要的作用。由于液压系统是直接影响飞机安全性的关键系统，不允许因为单个系统的失效而可能引发机毁人亡的事故。飞机液压系统被设计成多套相互独立的冗余系统，在每个单套的液压系统中又有多个液压泵以提高整个液压系统的可靠度。目前液压系统仍是飞机飞行操纵系统、刹车系统、起落架收放系统的主要动力来源。对于全电/多电飞机概念的兴起，在过去的40多年中由于发展不够成熟，使得电传的电作动系统相比于液压作动稍处于劣势，如表1所示：
[0115]
表1电作动与液压作动
[0116][0117]
通过比较后可以得到液压操纵部件由于工作刚度大、快速、执行机构有自润滑特性且故障回中等优势，在如今仍是飞控作动系统首选的。随着民航飞机的不断发展，对安全性和功重比的追求越来越高。为了实现飞机更好的性能和更高的安全性，模拟和数字技术在液压系统监控及执行控制功能的运用更加广泛。液压系统故障预测和智能变压力系统及健康管理等新技术也在不断地研究完善中，用于克服液压系统存在的某些常规缺点，提高液压系统的可靠性。
[0118]
空客相似机型：a319、a320、a321、a330，根据故障保留清单数据统计如表2所示：
[0119]
表2空客保留故障清单数据表
[0120][0121]
波音相似机型：b737、b757、b787，根据故障保留清单数据统计如表3所示：
[0122]
表3波音保留故障清单数据表
[0123][0124]
相似机型保留故障清单中进行保留的故障多数都发生在设施/设备及灯光系统，是由于被经常使用而导致容易损坏，且这些部件对飞机的适航并不存在较大的影响。空调、通讯、起落架、apu等系统的故障保留大多数原因也都与气源、电源、液压系统有直接关系。从故障保留清单中可知液压系统具有较高的可靠性进行保留的故障占比小，但是在多数情况下如果液压系统出现故障都会直接影响飞机运行的安全性。同时在飞机的整体结构中涉及液压的关键系统较多重要性也较大，所以液压系统的故障保留要求较为复杂，下面对液压系统进行深度的分析，如图1所示：
[0125]
液压作动系统如图2所示：
[0126]
对相似机型作动系统的液压部件组成进行研究对比分析得到下表如表4所示：
[0127]
表4液压作动系统部件组成表
[0128]
系统名称液压相关组成部件起落架电磁阀、电液换向阀、溢流阀、单向阀、优先阀等扰流板电磁阀、电液换向阀、溢流阀、节流阀、多功能扰流板液压缸等副翼电磁阀、电液换向阀、溢流阀、副翼液压缸、比例调速阀等襟翼电磁阀、电液换向阀、溢流阀、襟翼液压缸、比例调速阀等缝翼电磁阀、电液换向阀、溢流阀、缝翼液压缸、比例调速阀等升降舵电磁阀、电液换向阀、溢流阀、升降舵液压缸、比例调速阀等方向舵电磁阀、电液换向阀、溢流阀、方向舵液压缸、比例调速阀等刹车系统电磁阀、电液换向阀、溢流阀、单向阀、刹车液压缸等
[0129]
将波音和空客的相似机型结构系统进行详细的分析对比后发现相似机型作动系统中各系统的液压组成都较为相似，故进行具体详细分析。
[0130]
4.1.2主要故障模式及安全性影响分析
[0131]
相似机型系统的故障模式与所需完成的功能直接相关。对于系统的安全性分析是在确定主要的故障模式、故障等级及对应的安全性展开分析。通过系统结构的指导性文件进行系统的功能危险性分析(fha)，在分析系统工作原理及功能的基础上，结合fha分析得到的系统灾难性、危险性及影响较大的故障模式，确定为故障树的顶事件。通过对每个底事件重要度的判断，根据所得到的历史数据进行故障发生概率的估计并对发生故障后的影响进行判断。以液压系统为例进行分析，如表5所示；
[0132]
表5液压系统故障事件fha表
[0133]
故障编号故障事件名称故障等级指标要求h1正常刹车系统失效i级灾难的10-9
/飞行小时h2襟翼收放不正常ii级危险的10-7
/飞行小时h3副翼收放不正常ii级危险的10-7
/飞行小时h4方向舵无法偏转ii级危险的10-7
/飞行小时h5扰流板收放不正常ii级危险的10-7
/飞行小时h6发动机反推无法打开i级灾难的10-9
/飞行小时h7备用刹车系统失效ii级危险的10-7
/飞行小时h8前轮转弯失效iii级主要的10-5
/飞行小时
[0134]
建立相似机型液压系统故障树进行结构重要度的详细分析计算。
[0135]
空客a320相似机型绿液压系统失效故障树如图3至图8所示：
[0136]
在图3中，各参数含义如下，t1:正常刹车系统无压力或压力低；x1:正常刹车活门失效；x2:电磁阀卡滞；x3:电液换向阀失效；x4:流量限制器故障；x5:刹车作动筒卡阻；x6:传感器故障；x7:连接管路渗漏。
[0137]
由上图的故障树可计算出最小割集，最小割集代表系统的危险性，每个最小割集都是直接导致顶事件故障发生的集合。最小割集越多说明顶事件发生事故的危险性越大。故选取最小割集进行结构重要度系数的研究计算。
[0138]
t1＝x1+x2+x3+x4+x5+x6+x7[0139]
通过上式计算得到最小割集。即：p1＝{x1,x2,x3,x4,x5,x6,x7}。
[0140]
对于仅在同一个最小割集中出现的所有基本事件，且在其他最小割集中不再出现，则此割集中的基本事件结构重要度相等，故对最小路集p1的基本事件进行结构重要度系数的排序，可得到：
[0141]iφ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)
[0142]
结构重要度是原件关于系统临界的状态向量x的数量与所有2n个状态向量数值的比率，即原件的临界割集(路集)向量对2n个状态向量数值的比率。根据以下公式进行基本事件结构重要度近似值的计算：
[0143][0144]
式中：
[0145]
基本事件x1、x2、x3、x4、x5、x6、x7同处在p1中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝1/26≈0.0156。
[0146]
在图4中，各参数含义如下，t2:襟翼收放不正常；x1:液压泵故障；x2:连接管路渗漏；x3:节流阀卡阻；x4:传感器故障；x5:作动筒卡阻。
[0147]
同上故障树计算：
[0148]
t2＝x1+x2+x3+x4+x5[0149]
通过上式计算得到最小割集。即：p2＝{x1,x2,x3,x4,x5}。
[0150]
基本事件x1、x2、x3、x4、x5同处在p2中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝1/24＝0.0625。
[0151]
在图5中，各参数含义如下，t3:副翼收放不正常；x1:液压泵故障；x2:连接管路渗漏x3:伺服阀卡阻；x4:蓄压器故障；x5:作动筒卡阻；x6:传感器故障；x7:电磁阀卡滞。
[0152]
同上故障树计算：
[0153]
t3＝x1+x2+x3+x4+x5+x6+x7[0154]
通过上式计算得到最小割集。即：p3＝{x1,x2,x3,x4,x5,x6,x7}。
[0155]
基本事件x1、x2、x3、x4、x5、x6、x7同处在p3中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝1/26≈0.0156。
[0156]
在图6中，各参数含义如下，t4:方向舵无法偏转；x1:电磁阀卡滞；x2:连接管路渗漏x3:释压阀卡在开位；x4:伺服阀卡阻；x5:作动筒卡阻；x6:传感器故障。
[0157]
同上故障树计算：
[0158]
t4＝x1+x2+x3+x4+x5+x6[0159]
通过上式计算得到最小割集。即：p4＝{x1,x2,x3,x4,x5,x6}。
[0160]
基本事件x1、x2、x3、x4、x5、x6同处在p4中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝1/25≈0.0313。
[0161]
在图7中，各参数含义如下，t5:扰流板(1、5)收放不正常；x1:关断阀卡阻x2:连接管路渗漏；x3:电磁阀卡滞；x4:传感器故障；x5:作动筒卡阻。
[0162]
同上故障树计算：
[0163]
t5＝x1+x2+x3+x4+x5[0164]
通过上式计算得到最小割集。即：p5＝{x1,x2,x3,x4,x5}。
[0165]
基本事件x1、x2、x3、x4、x5同处在p5中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝1/24＝0.0625。
[0166]
在图8中，各参数含义如下，t6:一号发动机反推无法打开；x1:液压泵故障x2:连接管路渗漏；x3:释压阀卡在开位；x4:电磁阀卡滞；x5:作动筒卡阻；x6:传感器故障x7:关断阀卡阻。
[0167]
同上故障树计算：
[0168]
t6＝x1+x2+x3+x4+x5+x6+x7[0169]
通过上式计算得到最小割集。即：p6＝{x1,x2,x3,x4,x5,x6,x7}。
[0170]
基本事件x1、x2、x3、x4、x5、x6、x7同处在p6中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝1/26≈0.0156。
[0171]
波音b737ng相似机型b液压系统失效故障树如图9至图13所示：
[0172]
在图9中，各参数含义如下，t7:正常刹车系统无压力或压力低；x1:正常刹车活门失效；x2:电磁阀卡滞；x3:电液换向阀失效；x4:流量限制器故障；x5:刹车作动筒卡阻；x6:传感器故障；x7:连接管路渗漏；x8:释压阀卡在开位。
[0173]
同上故障树计算：
[0174]
t7＝x1+x2+x3+x4+x5+x6+x7+x8[0175]
通过上式计算得到最小割集。即：p7＝{x1,x2,x3,x4,x5,x6,x7,x8}。
[0176]
基本事件x1、x2、x3、x4、x5、x6、x7、x8同处在p7中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝i
φ
(x8)＝1/27≈0.0078。
[0177]
在图10中，各参数含义如下，t8:襟翼收放不正常；x1:拉力管、扭力管卡阻x2:电磁
阀卡滞；x3:传感器故障；x4:液压马达故障；x5:连接管路渗漏x6:作动筒卡阻；x7:流量限制器故障；x8:液压泵故障；x9:电液换向阀失效。
[0178]
同上故障树计算：
[0179]
t8＝x1+x2+x3+x4+x5+x6+x7+x8+x9[0180]
通过上式计算得到最小割集。即：p8＝{x1,x2,x3,x4,x5,x6,x7,x8,x9}。
[0181]
基本事件x1、x2、x3、x4、x5、x6、x7、x8同处在p8中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝i
φ
(x8)＝i
φ
(x9)＝1/28≈0.0039。
[0182]
在图11中，各参数含义如下，t9:副翼收放不正常；x1:单向阀故障；x2:电磁阀卡滞x3:电液换向阀失效；x4:流量限制器故障；x5:作动筒卡阻x6:传感器故障x7:连接管路渗漏；x8:释压阀卡在开位。
[0183]
同上故障树计算：
[0184]
t9＝x1+x2+x3+x4+x5+x6+x7+x8[0185]
通过上式计算得到最小割集。即：p9＝{x1,x2,x3,x4,x5,x6,x7,x8}。
[0186]
基本事件x1、x2、x3、x4、x5、x6、x7、x8同处在p9中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝i
φ
(x8)＝1/27≈0.0078。
[0187]
在图12中，各参数含义如下，t
10
:方向舵无法偏转；x1:单向阀故障；x2:电磁阀卡滞x3:释压阀卡在开位；x4:伺服阀故障；x5:关断阀卡阻x6:作动筒卡阻x7:传感器故障；x8:连接管路渗漏。
[0188]
同上故障树计算：
[0189]
t
10
＝x1+x2+x3+x4+x5+x6+x7+x8[0190]
通过上式计算得到最小割集。即：p
10
＝{x1,x2,x3,x4,x5,x6,x7,x8}。
[0191]
基本事件x1、x2、x3、x4、x5、x6、x7、x8同处在p
10
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝i
φ
(x8)＝1/27≈0.0078。
[0192]
在图13中，各参数含义如下，t
11
:右发动机反推无法打开；x1:液压泵故障x2:连接管路渗漏；x3:控制阀故障；x4:电磁阀卡滞；x5:作动筒卡阻；x6:传感器故障。
[0193]
同上故障树计算：
[0194]
t
11
＝x1+x2+x3+x4+x5+x6[0195]
通过上式计算得到最小割集。即：p
11
＝{x1,x2,x3,x4,x5,x6}。
[0196]
基本事件x1、x2、x3、x4、x5、x6同处在p
11
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝1/25≈0.0313。
[0197]
空客a320相似机型黄液压系统失效故障树如图14至图18所示：
[0198]
在图14中，各参数含义如下，t
12
:备用刹车无法保持压力；m1:蓄压器压力低；x1:单向阀故障；x2:释压阀卡在开位；x3:流量限制器故障；x4:电磁阀卡滞；x5:传感器故障x6:连接管路渗漏；x7:作动筒卡阻。
[0199]
同上故障树计算：
[0200]
t
12
＝x1+x2+x3+x4+x5+x6+x7[0201]
通过上式计算得到最小割集。即：p
12
＝{x1,x2,x3,x4,x5,x6,x7}。
[0202]
基本事件x1、x2、x3、x4、x5、x6同处在p
12
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝1/26≈0.0156。
[0203]
在图15中，各参数含义如下，t
13
:方向舵无法偏转；x1:作动弹簧故障x2:连接管路渗漏；x3:释压阀卡在开位；x4:控制阀卡阻；x5:作动筒卡阻x6:传感器故障。
[0204]
同上故障树计算：
[0205]
t
13
＝x1+x2+x3+x4+x5+x6[0206]
通过上式计算得到最小割集。即：p
13
＝{x1,x2,x3,x4,x5,x6}。
[0207]
基本事件x1、x2、x3、x4、x5、x6同处在p
13
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝1/25≈0.0313。
[0208]
在图16中，各参数含义如下，t
14
:扰流板(2、4)收放不正常；x1:作动弹簧故障x2:连接管路渗漏；x3:关断阀卡阻；x4:电磁阀卡滞；x5:作动筒卡阻；x6:传感器故障。
[0209]
同上故障树计算：
[0210]
t
14
＝x1+x2+x3+x4+x5+x6[0211]
通过上式计算得到最小割集。即：p
14
＝{x1,x2,x3,x4,x5,x6}。
[0212]
基本事件x1、x2、x3、x4、x5、x6同处在p
14
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝1/25≈0.0313。
[0213]
在图17中，各参数含义如下，t
15
:二号发动机反推无法打开；x1:液压泵故障x2:连接管路渗漏；x3:释压阀卡在开位；x4:电磁阀卡滞；x5:作动筒卡阻；x6:传感器故障x7:关断阀卡阻。
[0214]
同上故障树计算：
[0215]
t
15
＝x1+x2+x3+x4+x5+x6+x7[0216]
通过上式计算得到最小割集。即：p
15
＝{x1,x2,x3,x4,x5,x6,x7}。
[0217]
基本事件x1、x2、x3、x4、x5、x6、x7同处在p
15
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝1/26≈0.0156。
[0218]
在图18中，各参数含义如下，t
16
:前轮转弯失效；x1:单向阀故障；x2:电磁阀卡滞；x3:伺服阀故障；x4:流量限制器故障；x5:作动筒卡阻；x6:传感器故障x7:连接管路渗漏；x8:蓄压器故障。
[0219]
同上故障树计算：
[0220]
t
16
＝x1+x2+x3+x4+x5+x6+x7+x8[0221]
通过上式计算得到最小割集。即：p
16
＝{x1,x2,x3,x4,x5,x6,x7,x8}。
[0222]
基本事件x1、x2、x3、x4、x5、x6、x7、x8同处在p
16
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝i
φ
(x8)＝1/27≈0.0078。
[0223]
波音b737ng相似机型a液压系统失效故障树如图19至图22所示：
[0224]
在图19中，各参数含义如下，t
17
:备用刹车系统压力无法保持；x1:正常刹车活门失效x2:电磁阀卡滞；x3:电液换向阀失效；x4:流量限制器故障；x5:刹车作动筒卡阻；x6:传感器故障；x7:连接管路渗漏；x8:释压阀卡在开位。
[0225]
同上故障树计算：
[0226]
t
17
＝x1+x2+x3+x4+x5+x6+x7+x8[0227]
通过上式计算得到最小割集。即：p
17
＝{x1,x2,x3,x4,x5,x6,x7,x8}。
[0228]
基本事件x1、x2、x3、x4、x5、x6、x7、x8同处在p
17
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝i
φ
(x8)＝1/27≈0.0078。
[0229]
在图20中，各参数含义如下，t
18
:扰流板收放不正常；x1:释压阀卡在开位x2:连接管路渗漏；x3:控制阀卡阻；x4:电磁阀卡滞；x5:作动筒卡阻；x6:传感器故障。
[0230]
同上故障树计算：
[0231]
t
18
＝x1+x2+x3+x4+x5+x6[0232]
通过上式计算得到最小割集。即：p
18
＝{x1,x2,x3,x4,x5,x6}。
[0233]
基本事件x1、x2、x3、x4、x5、x6同处在p
18
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝1/25≈0.0313。
[0234]
在图21中，各参数含义如下，t
19
:左发动机反推无法打开；x1:液压泵故障x2:连接管路渗漏；x3:控制活门故障；x4:电磁阀卡滞；x5:作动筒卡阻；x6:传感器故障。
[0235]
同上故障树计算：
[0236]
t
19
＝x1+x2+x3+x4+x5+x6[0237]
通过上式计算得到最小割集。即：p
19
＝{x1,x2,x3,x4,x5,x6}。
[0238]
基本事件x1、x2、x3、x4、x5、x6同处在p
19
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝1/25≈0.0313。
[0239]
在图22中，各参数含义如下，t
20
:前轮转弯失效；x1:单向阀故障；x2:连接管路渗漏；x3:流量限制器故障；x4:传感器故障；x5:作动筒卡阻。
[0240]
同上故障树计算：
[0241]
t
20
＝x1+x2+x3+x4+x5[0242]
通过上式计算得到最小割集。即：p
20
＝{x1,x2,x3,x4,x5}。
[0243]
基本事件x1、x2、x3、x4、x5同处在p
20
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝1/24＝0.0625。
[0244]
空客a320相似机型蓝液压系统失效故障树如图4-23至图4-26所示：
[0245]
在图23中，各参数含义如下，t
21
:襟缝翼收放不正常；x1:释压阀卡在开位；x2:连接管路渗漏；x3:控制活门故障；x4:电磁阀卡滞；x5:作动筒卡阻x6:传感器故障。
[0246]
同上故障树计算：
[0247]
t
21
＝x1+x2+x3+x4+x5+x6[0248]
通过上式计算得到最小割集。即：p
21
＝{x1,x2,x3,x4,x5,x6}。
[0249]
基本事件x1、x2、x3、x4、x5、x6同处在p
21
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝1/25≈0.0313。
[0250]
在图24中，各参数含义如下，t
22
:副翼收放不正常；x1:液压泵故障x2:连接管路渗漏；x3:伺服阀卡阻；x4:蓄压器故障；x5:作动筒卡阻x6:传感器故障x7:电磁阀卡滞。
[0251]
同上故障树计算：
[0252]
t
22
＝x1+x2+x3+x4+x5+x6+x7[0253]
通过上式计算得到最小割集。即：p
22
＝{x1,x2,x3,x4,x5,x6,x7}。
[0254]
基本事件x1、x2、x3、x4、x5、x6、x7同处在p
22
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝1/26≈0.0156。
[0255]
在图25中，各参数含义如下，t
23
:方向舵无法偏转；x1:作动弹簧故障x2:连接管路渗漏；x3:释压阀卡在开位；x4:控制阀卡阻；x5:作动筒卡阻；x6:传感器故障。
[0256]
同上故障树计算：
[0257]
t
23
＝x1+x2+x3+x4+x5+x6[0258]
通过上式计算得到最小割集。即：p
23
＝{x1,x2,x3,x4,x5,x6}。
[0259]
基本事件x1、x2、x3、x4、x5、x6同处在p
23
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝1/25≈0.0313。
[0260]
在图26中，各参数含义如下，t
24
:扰流板(3)收放不正常；x1:作动弹簧故障x2:连接管路渗漏；x3:关断阀卡阻；x4:电磁阀卡滞；x5:作动筒卡阻；x6:传感器故障。
[0261]
同上故障树计算：
[0262]
t
24
＝x1+x2+x3+x4+x5+x6[0263]
通过上式计算得到最小割集。即：p
24
＝{x1,x2,x3,x4,x5,x6}。
[0264]
基本事件x1、x2、x3、x4、x5、x6同处在p
24
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝1/25≈0.0313。
[0265]
波音b737ng相似机型备用液压系统失效故障树如图4-27至图4-29所示：
[0266]
在图27中，各参数含义如下，t
25
:襟缝翼收放不正常；x1:单向阀故障；x2:连接管路渗漏；x3:释压阀卡在开位；x4:关断阀卡阻；x5:作动筒卡阻；x6:传感器故障x7:流量限制器故障。
[0267]
同上故障树计算：
[0268]
t
25
＝x1+x2+x3+x4+x5+x6+x7[0269]
通过上式计算得到最小割集。即：p
25
＝{x1,x2,x3,x4,x5,x6,x7}。
[0270]
基本事件x1、x2、x3、x4、x5、x6、x7同处在p
25
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝1/26≈0.0156。
[0271]
在图28中，各参数含义如下，t
26
:备用方向舵无法偏转；x1:单向阀故障；x2:连接管路渗x3:控制阀故障；x4:关断阀卡阻；x5:作动筒卡阻x6:传感器故障；x7:电磁阀卡滞。
[0272]
同上故障树计算：
[0273]
t
26
＝x1+x2+x3+x4+x5+x6+x7[0274]
通过上式计算得到最小割集。即：p
26
＝{x1,x2,x3,x4,x5,x6,x7}。
[0275]
基本事件x1、x2、x3、x4、x5、x6、x7同处在p
26
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝i
φ
(x7)＝1/26≈0.0156。
[0276]
在图29中，各参数含义如下，t
27
:发动机反推无法打开；x1:液压泵故障x2:连接管路渗漏；x3:控制阀故障；x4:电磁阀卡滞；x5:作动筒卡阻；x6:传感器故障。
[0277]
同上故障树计算：
[0278]
t
27
＝x1+x2+x3+x4+x5+x6[0279]
通过上式计算得到最小割集。即：p
27
＝{x1,x2,x3,x4,x5,x6}。
[0280]
基本事件x1、x2、x3、x4、x5、x6同处在p
27
中，则i
φ
(x1)＝i
φ
(x2)＝i
φ
(x3)＝i
φ
(x4)＝i
φ
(x5)＝i
φ
(x6)＝1/25≈0.0313。
[0281]
故障底事件的故障数据如表6所示：
[0282]
表6相似机型液压刹车系统部件主要故障数据表
[0283][0284][0285]
故障树顶事件失效概率：
[0286]
p(t1)＝1-(1-p(y8))(1-p(y7))(1-p(y
13
))(1-p(y3))(1-p(y4))(1-p(y5))(1-p(y6))
[0287]
＝1-(1-0.04596)(1-0.02945)(1-0.02186)(1-0.01398)(1-0.00938)(1-0.07512)(1-0.01732)
[0288]
＝1-0.95404
×
0.97055
×
0.97814
×
0.98602
×
0.99062
×
0.92488
×
0.98268
[0289]
≈0.1960
[0290]
p(t2)＝1-(1-p(y1))(1-p(y6))(1-p(y2))(1-p(y5))(1-p(y4))
[0291]
＝1-(1-0.03487)(1-0.01732)(1-0.00926)(1-0.07512)(1-0.00938)
[0292]
＝1-0.96513
×
0.98268
×
0.99074
×
0.92488
×
0.99062
[0293]
≈0.1391
[0294]
p(t3)＝1-(1-p(y1))(1-p(y6))(1-p(y9))(1-p(y
10
))(1-p(y4))(1-p(y5))(1-p(y7))
[0295]
＝1-(1-0.03487)(1-0.01732)(1-0.00958)(1-0.01304)(1-0.00938)(1-0.07512)(1-0.02945)
[0296]
＝1-0.96513
×
0.98268
×
0.99042
×
0.98696
×
0.99062
×
0.92488
×
0.97055
[0297]
≈0.1756
[0298]
p(t4)＝1-(1-p(y7))(1-p(y6))(1-p(y
11
))(1-p(y9))(1-p(y4))(1-p(y5))
[0299]
＝1-(1-0.02945)(1-0.01732)(1-0.01064)(1-0.00958)(1-0.00938)(1-0.07512)
[0300]
＝1-0.97055
×
0.98268
×
0.98936
×
0.99042
×
0.99062
×
0.92488
[0301]
≈0.1438
[0302]
p(t5)＝1-(1-p(y
12
))(1-p(y6))(1-p(y7))(1-p(y5))(1-p(y4))
[0303]
＝1-(1-0.00916)(1-0.01732)(1-0.02945)(1-0.07512)(1-0.00938)
[0304]
＝1-0.99084
×
0.98268
×
0.97055
×
0.92488
×
0.99062
[0305]
≈0.1342
[0306]
p(t6)＝1-(1-p(y1))(1-p(y6))(1-p(y
11
))(1-p(y7))(1-p(y4))(1-p(y5))(1-p(y
12
))
[0307]
＝1-(1-0.03487)(1-0.01732)(1-0.01064)(1-0.02945)(1-0.00938)(1-0.07512)(1-0.00916)
[0308]
＝1-0.96513
×
0.98268
×
0.98936
×
0.97055
×
0.99062
×
0.92488
×
0.99084
[0309]
≈0.1733
[0310]
p(t7)＝1-(1-p(y8))(1-p(y7))(1-p(y1))(1-p(y3))(1-p(y4))(1-p(y5))(1-p(y6))(1-p(y
11
))
[0311]
＝1-0.96288
×
0.97222
×
0.97124
×
0.98572
×
0.99028
×
0.92428
×
0.98268
×
0.98936
[0312]
≈0.2151
[0313]
p(t8)＝1-(1-p(y
14
))(1-p(y7))(1-p(y5))(1-p(y
15
))(1-p(y6))(1-p(y4))(1-p(y3))(1-p(y1))(1-p(y
13
))
[0314]
＝1-0.97632
×
0.97222
×
0.94283
×
0.98047
×
0.93084
×
0.99028
×
0.98572
×
0.97124
×
0.9811
[0315]
≈0.2359
[0316]
p(t9)＝1-(1-p(y
16
))(1-p(y7))(1-p(y
13
))(1-p(y3))(1-p(y4))(1-p(y5))(1-p(y6))(1-p(y
11
))
[0317]
＝1-0.97129
×
0.97222
×
0.9811
×
0.98572
×
0.99028
×
0.94283
×
0.93084
×
0.9874
[0318]
≈0.2163
[0319]
p(t
10
)＝1-(1-p(y
16
))(1-p(y7))(1-p(y
11
))(1-p(y9))(1-p(y
12
))(1-p(y4))(1-p(y5))(1-p(y6))
[0320]
＝1-0.9729
×
0.97222
×
0.9874
×
0.9968
×
0.99107
×
0.99028
×
0.94283
×
0.93084
[0321]
≈0.1981
[0322]
p(t
11
)＝1-(1-p(y1))(1-p(y6))(1-p(y
17
))(1-p(y7))(1-p(y4))(1-p(y5))
[0323]
＝1-(1-0.02876)(1-0.06916)(1-0.01046)(1-0.02778)(1-0.00972)(1-0.05717)
[0324]
＝1-0.97124
×
0.973084
×
0.98954
×
0.97222
×
0.99028
×
0.94283
[0325]
≈0.1511
[0326]
p(t
12
)＝1-(1-p(y
16
))(1-p(y
11
))(1-p(y3))(1-p(y7))(1-p(y5))(1-p(y6))(1-p(y4))
[0327]
＝1-(1-0.02927)(1-0.01064)(1-0.01398)(1-0.02945)(1-0.07512)(1-0.01732)(1-0.00938)
[0328]
＝1-0.97073
×
0.98936
×
0.98602
×
0.97055
×
0.92488
×
0.98268
×
0.99062
[0329]
≈0.1770
[0330]
p(t
13
)＝1-(1-p(y
18
))(1-p(y6))(1-p(y
11
))(1-p(y
17
))(1-p(y4))(1-p(y5))
[0331]
＝1-(1-0.00764)(1-0.01732)(1-0.01064)(1-0.00956)(1-0.00938)(1-0.07512)
[0332]
＝1-0.97055
×
0.98268
×
0.98936
×
0.99042
×
0.99062
×
0.92488
[0333]
≈0.1438
[0334]
p(t
14
)＝1-(1-p(y
18
))(1-p(y6))(1-p(y
12
))(1-p(y7))(1-p(y4))(1-p(y5))
[0335]
＝1-(1-0.00764)(1-0.01732)(1-0.00916)(1-0.02945)(1-0.00938)(1-0.07512)
[0336]
＝1-0.97055
×
0.98268
×
0.99084
×
0.97055
×
0.99062
×
0.92488
[0337]
≈0.1597
[0338]
p(t
15
)＝1-(1-p(y1))(1-p(y6))(1-p(y
11
))(1-p(y7))(1-p(y4))(1-p(y5))(1-p(y
12
))
[0339]
＝1-(1-0.03487)(1-0.01732)(1-0.01064)(1-0.02945)(1-0.00938)(1-0.07512)(1-0.00916)
[0340]
＝1-0.96513
×
0.98268
×
0.98936
×
0.97055
×
0.99062
×
0.92488
×
0.99084
[0341]
≈0.1733
[0342]
p(t
16
)＝1-(1-p(y
16
))(1-p(y7))(1-p(y9))(1-p(y3))(1-p(y4))(1-p(y5))(1-p(y6))(1-p(y
10
))
[0343]
＝1-0.97073
×
0.97055
×
0.99042
×
0.98602
×
0.99062
×
0.92488
×
0.98268
×
0.98696
[0344]
≈0.1824
[0345]
p(t
17
)＝1-(1-p(y8))(1-p(y7))(1-p(y
13
))(1-p(y3))(1-p(y4))(1-p(y5))(1-p(y6))(1-p(y
11
))
[0346]
＝1-0.96288
×
0.97222
×
0.9811
×
0.98572
×
0.99028
×
0.94283
×
0.98696
×
0.9874
[0347]
≈0.1763
[0348]
p(t
18
)＝1-(1-p(y
11
))(1-p(y6))(1-p(y
17
))(1-p(y7))(1-p(y4))(1-p(y5))
[0349]
＝1-(1-0.0126)(1-0.06916)(1-0.01064)(1-0.02778)(1-0.00972)(1-0.05717)
[0350]
＝1-0.9874
×
0.93084
×
0.98954
×
0.97222
×
0.99028
×
0.94283
[0351]
≈0.1744
[0352]
p(t
19
)＝1-(1-p(y1))(1-p(y6))(1-p(y
17
))(1-p(y7))(1-p(y4))(1-p(y5))
[0353]
＝1-(1-0.02876)(1-0.06916)(1-0.01064)(1-0.02778)(1-0.00972)(1-0.05717)
[0354]
＝1-0.97124
×
0.93084
×
0.98954
×
0.97222
×
0.99028
×
0.94283
[0355]
≈0.1879
[0356]
p(t
20
)＝1-(1-p(y
16
))(1-p(y6))(1-p(y3))(1-p(y4))(1-p(y5))
[0357]
＝1-(1-0.02871)(1-0.06916)(1-0.01428)(1-0.00972)(1-0.05717)
[0358]
＝1-0.97129
×
0.93084
×
0.98572
×
0.99028
×
0.94283
[0359]
≈0.1679
[0360]
p(t
21
)＝1-(1-p(y
11
))(1-p(y6))(1-p(y
17
))(1-p(y7))(1-p(y4))(1-p(y5))
[0361]
＝1-(1-0.01064)(1-0.01732)(1-0.00956)(1-0.02945)(1-0.00938)(1-0.07512)
[0362]
＝1-0.98936
×
0.98268
×
0.99044
×
0.97055
×
0.99062
×
0.92488
[0363]
≈0.1437
[0364]
p(t
22
)＝1-(1-p(y1))(1-p(y6))(1-p(y9))(1-p(y7))(1-p(y4))(1-p(y5))(1-p(y
10
))
[0365]
＝1-(1-0.03487)(1-0.01732)(1-0.00958)(1-0.02945)(1-0.00938)(1-0.07512)(1-0.01304)
[0366]
＝1-0.96513
×
0.98268
×
0.99042
×
0.97055
×
0.99062
×
0.92488
×
0.98696
[0367]
≈0.1756
[0368]
p(t
23
)＝1-(1-p(y
18
))(1-p(y6))(1-p(y
11
))(1-p(y
17
))(1-p(y4))(1-p(y5))
[0369]
＝1-(1-0.00764)(1-0.01732)(1-0.01064)(1-0.00956)(1-0.00938)(1-0.07512)
[0370]
＝1-0.99236
×
0.98268
×
0.98936
×
0.99044
×
0.99062
×
0.92488
[0371]
≈0.1245
[0372]
p(t
24
)＝1-(1-p(y
18
))(1-p(y6))(1-p(y
12
))(1-p(y7))(1-p(y4))(1-p(y5))
[0373]
＝1-(1-0.00764)(1-0.01732)(1-0.00916)(1-0.02945)(1-0.00938)(1-0.07512)
[0374]
＝1-0.99236
×
0.98268
×
0.99084
×
0.97055
×
0.99062
×
0.92488
[0375]
≈0.1408
[0376]
p(t
25
)＝1-(1-p(y
16
))(1-p(y6))(1-p(y
11
))(1-p(y
12
))(1-p(y4))(1-p(y5))(1-p(y3))
[0377]
＝1-(1-0.02871)(1-0.06916)(1-0.0126)(1-0.00893)(1-0.00972)(1-0.05717)(1-0.01428)
[0378]
＝1-0.97129
×
0.93084
×
0.9874
×
0.99107
×
0.99028
×
0.94283
×
0.98572
[0379]
≈0.1857
[0380]
p(t
26
)＝1-(1-p(y
16
))(1-p(y6))(1-p(y
17
))(1-p(y
12
))(1-p(y4))(1-p(y5))(1-p(y7))
[0381]
＝1-(1-0.02871)(1-0.06916)(1-0.01046)(1-0.00893)(1-0.00972)(1-0.05717)(1-0.02778)
[0382]
＝1-0.97129
×
0.93084
×
0.98954
×
0.99107
×
0.99028
×
0.94283
×
0.97222
[0383]
≈0.1951
[0384]
p(t
27
)＝1-(1-p(y1))(1-p(y6))(1-p(y
17
))(1-p(y4))(1-p(y5))(1-p(y7))
[0385]
＝1-(1-0.02876)(1-0.06916)(1-0.01046)(1-0.00972)(1-0.05717)(1-0.02778)
[0386]
＝1-0.97124
×
0.93084
×
0.98954
×
0.99028
×
0.94283
×
0.97222
[0387]
≈0.1879
[0388]
故障树底事件概率重要度：
[0389]
概率重要度是表示底事件xi发生概率变化引起顶事件t发生概率变化的变化程度，用顶事件p(t)的发生概率对底事件xi发生概率p(xi)的偏导数来表示。
[0390]
概率重要度计算公式：
[0391][0392]
式中i
p(i)
为第i个基本事件的概率重要度，p(xi)为第i个事件发生的概率。
[0393]
根据公式3-4可求得t1中部件概率重要度
[0394][0395][0396][0397][0398][0399][0400][0401]
根据公式3.4可求得t2中部件概率重要度
[0402][0403][0404][0405][0406][0407]
根据公式3.4可求得t3中部件概率重要度
[0408][0409][0410][0411]
[0412][0413][0414][0415]
根据公式3.4可求得t4中部件概率重要度
[0416][0417][0418][0419][0420][0421]
根据公式3.4可求得t5中部件概率重要度
[0422][0423][0424][0425][0426][0427]
根据公式3.4可求得t6中部件概率重要度
[0428][0429][0430][0431][0432][0433]
[0434][0435]
根据公式3.4可求得t7中部件概率重要度
[0436][0437][0438][0439][0440][0441][0442][0443][0444]
根据公式3.4可求得t8中部件概率重要度
[0445][0446][0447][0448][0449][0450][0451][0452][0453][0454]
根据公式3.4可求得t9中部件概率重要度
[0455][0456][0457]
[0458][0459][0460][0461][0462][0463]
根据公式3.4可求得t
10
中部件概率重要度
[0464][0465][0466][0467][0468][0469][0470][0471][0472]
根据公式3.4可求得t
11
中部件概率重要度
[0473][0474][0475][0476][0477][0478][0479]
根据公式3.4可求得t
12
中部件概率重要度
[0480][0481]
[0482][0483][0484][0485][0486][0487]
根据公式3.4可求得t
13
中部件概率重要度
[0488][0489][0490][0491][0492][0493][0494]
根据公式3.4可求得t
14
中部件概率重要度
[0495][0496][0497][0498][0499][0500][0501]
根据公式3.4可求得t
15
中部件概率重要度
[0502][0503][0504]
[0505][0506][0507][0508][0509]
根据公式3.4可求得t
16
中部件概率重要度
[0510][0511][0512][0513][0514][0515][0516][0517][0518]
根据公式3.4可求得t
17
中部件概率重要度
[0519][0520][0521][0522][0523][0524][0525][0526][0527]
根据公式3.4可求得t
18
中部件概率重要度
[0528]
[0529][0530][0531][0532][0533][0534]
根据公式3.4可求得t
19
中部件概率重要度
[0535][0536][0537][0538][0539][0540][0541]
根据公式3.4可求得t
20
中部件概率重要度
[0542][0543][0544][0545][0546][0547]
根据公式3.4可求得t
21
中部件概率重要度
[0548][0549][0550][0551]
[0552][0553][0554]
根据公式3.4可求得t
22
中部件概率重要度
[0555][0556][0557][0558][0559][0560][0561][0562]
根据公式3.4可求得t
23
中部件概率重要度
[0563][0564][0565][0566][0567][0568][0569]
根据公式3.4可求得t
24
中部件概率重要度
[0570][0571][0572][0573][0574]
[0575][0576]
根据公式3.4可求得t
25
中部件概率重要度
[0577][0578][0579][0580][0581][0582][0583][0584]
根据公式3.4可求得t
26
中部件概率重要度
[0585][0586][0587][0588][0589][0590][0591][0592]
根据公式3.4可求得t
27
中部件概率重要度
[0593][0594][0595][0596][0597]
[0598][0599]
关键重要度指的是底事件xi故障概率的变化率与它引起的顶事件t发生概率p(t)的变化率之比。表达式为：
[0600][0601]
其中i
c(i)
为第i个基本事件的关键重要度。
[0602]
由公式3.5可求得t1中
[0603]
同理可得到下列底事件的关键重要度：
[0604]ic(7)
＝0.1245 i
c(13)
＝0.0917 i
c(3)
＝0.0520 i
c(4)
＝0.0388 i
c(5)
＝0.3332 i
c(6)
＝0.0723
[0605]
由公式3.5可求得t2中
[0606]
同理可得到下列底事件的关键重要度：
[0607]ic(6)
＝0.1091 i
c(2)
＝0.0578 i
c(5)
＝0.5027 i
c(4)
＝0.0586
[0608]
由公式3.5可求得t3中
[0609]
同理可得到下列底事件的关键重要度：
[0610]ic(6)
＝0.0827 i
c(9)
＝0.0454 i
c(10)
＝0.0620 i
c(4)
＝0.0445 i
c(5)
＝0.3813 i
c(7)
＝0.1425
[0611]
由公式3.5可求得t4中
[0612]
同理可得到下列底事件的关键重要度：
[0613]ic(6)
＝0.1049 i
c(11)
＝0.0640 i
c(9)
＝0.0576 i
c(4)
＝0.0543 i
c(5)
＝0.4836
[0614]
由公式3.5可求得t5中
[0615]
同理可得到下列底事件的关键重要度：
[0616]ic(6)
＝0.1137 i
c(7)
＝0.1958 i
c(5)
＝0.5240 i
c(4)
＝0.0608
[0617]
由公式3.5可求得t6中
[0618]
同理可得到下列底事件的关键重要度：
[0619]ic(6)
＝0.0841i
c(11)
＝0.1054 i
c(7)
＝0.1448 i
c(4)
＝0.0452 i
c(5)
＝0.3874 i
c(12)
＝0.0441
[0620]
由公式3.5可求得t7中
[0621]
同理可得到下列底事件的关键重要度：
[0622]ic(7)
＝0.1043 i
c(13)
＝0.0703 i
c(3)
＝0.0529 i
c(4)
＝0.0358 i
c(5)
＝0.2213
[0623]ic(6)
＝0.2711
[0624]ic(11)
＝0.0096
[0625]
由公式3.5可求得t8中
[0626]
同理可得到下列底事件的关键重要度：
[0627]ic(7)
＝0.0925 i
c(5)
＝0.1964 i
c(15)
＝0.0645 i
c(6)
＝0.2407 i
c(4)
＝0.0318 i
c(3)
＝0.0480 i
c(1)
＝0.0959 i
c(13)
＝0.0624
[0628]
由公式3.5可求得t9中
[0629]
同理可得到下列底事件的关键重要度：
[0630]ic(7)
＝0.1035 i
c(13)
＝0.0700 i
c(3)
＝0.0525 i
c(4)
＝0.0351 i
c(5)
＝0.219 7i
c(6)
＝0.2692 i
c(11)
＝0.0462
[0631]
由公式3.5可求得t
10
中
[0632]
同理可得到下列底事件的关键重要度：
[0633]ic(7)
＝0.1164 i
c(11)
＝0.0518 i
c(9)
＝0.0382 i
c(12)
＝0.0366 i
c(4)
＝0.0398 i
c(5)
＝0.2489 i
c(6)
＝0.3059
[0634]
由公式3.5可求得t
11
中
[0635]
同理可得到下列底事件的关键重要度：
[0636]ic(6)
＝0.4174 i
c(17)
＝0.0594 i
c(7)
＝0.1605 i
c(4)
＝0.0551 i
c(5)
＝0.3406
[0637]
由公式3.5可求得t
12
中
[0638]
同理可得到下列底事件的关键重要度：
[0639]ic(11)
＝0.0522 i
c(3)
＝0.0688 i
c(7)
＝0.1478 i
c(5)
＝0.3996 i
c(6)
＝0.0856 i
c(4)
＝0.0459
[0640]
由公式3.5可求得t
13
中
[0641]
同理可得到下列底事件的关键重要度：
[0642]ic(6)
＝0.1049 i
c(11)
＝0.0640 i
c(17)
＝0.0575 i
c(4)
＝0.0564 i
c(5)
＝0.4836
[0643]
由公式3.5可求得t
14
中
[0644]
同理可得到下列底事件的关键重要度：
[0645]ic(6)
＝0.0927 i
c(12)
＝0.0486 i
c(7)
＝0.1578 i
c(4)
＝0.0498 i
c(5)
＝0.4274
[0646]
由公式3.5可求得t
15
中
[0647]
同理可得到下列底事件的关键重要度：
[0648]ic(6)
＝0.0841 i
c(11)
＝0.0513 i
c(7)
＝0.1448 i
c(4)
＝0.0452 i
c(5)
＝0.3874 i
c(12)
＝0.0441
[0649]
由公式3.5可求得t
16
中
[0650]
同理可得到下列底事件的关键重要度：
[0651]ic(7)
＝0.1360 i
c(9)
＝0.0434 i
c(3)
＝0.0436 i
c(4)
＝0.0424 i
c(5)
＝0.3641 i
c(6)
＝0.0790
[0652]ic(10)
＝0.0592
[0653]
由公式3.5可求得t
17
中
[0654]
同理可得到下列底事件的关键重要度：
[0655]ic(7)
＝0.3475 i
c(13)
＝0.0900 i
c(3)
＝0.0677 i
c(4)
＝0.0459 i
c(5)
＝0.2832 i
c(6)
＝0.3471
[0656]ic(11)
＝0.0596
[0657]
由公式3.5可求得t
18
中
[0658]
同理可得到下列底事件的关键重要度：
[0659]ic(6)
＝0.3513 i
c(17)
＝0.0500 i
c(7)
＝0.1353 i
c(4)
＝0.0465 i c(5)
＝0.2871
[0660]
由公式3.5可求得t
19
中
[0661]
同理可得到下列底事件的关键重要度：
[0662]ic(6)
＝0.3211 i
c(17)
＝0.0457 i
c(7)
＝0.1235 i
c(4)
＝0.0424 i
c(5)
＝0.2621
[0663]
由公式3.5可求得t
20
中
[0664]
同理可得到下列底事件的关键重要度：
[0665]ic(6)
＝0.3682 i
c(3)
＝0.0718 i
c(4)
＝0.0486 i
c(5)
＝0.3636
[0666]
由公式3.5可求得t
21
中
[0667]
同理可得到下列底事件的关键重要度：
[0668]ic(6)
＝0.1050 i
c(17)
＝0.0575 i
c(7)
＝0.1808 i
c(4)
＝0.0564 i
c(5)
＝0.4840
[0669]
由公式3.5可求得t
22
中
[0670]
同理可得到下列底事件的关键重要度：
[0671]ic(6)
＝0.0827 i
c(9)
＝0.0454 i
c(10)
＝0.0620 i
c(4)
＝0.0445 i
c(5)
＝0.3813 i
c(7)
＝0.1425
[0672]
由公式3.5可求得t
23
中
[0673]
同理可得到下列底事件的关键重要度：
[0674]ic(6)
＝0.1239 i
c(11)
＝0.0756 i
c(17)
＝0.0679 i
c(4)
＝0.0666 i
c(5)
＝0.5712
[0675]
由公式3.5可求得t
24
中
[0676]
同理可得到下列底事件的关键重要度：
[0677]ic(6)
＝0.1075 i
c(12)
＝0.0564 i
c(7)
＝0.1852 i
c(4)
＝0.0578 i
c(5)
＝0.4956
[0678]
由公式3.5可求得t
25
中
[0679]
同理可得到下列底事件的关键重要度：
[0680]ic(6)
＝0.3258 i
c(11)
＝0.0560 i
c(12)
＝0.0395 i
c(4)
＝0.0430 i
c(5)
＝0.2659 i
c(3)
＝0.0635
[0681]
由公式3.5可求得t
26
中
[0682]
同理可得到下列底事件的关键重要度：
[0683]ic(6)
＝0.3065 i
c(17)
＝0.0436 i
c(12)
＝0.0372 i
c(4)
＝0.0405 i
c(5)
＝0.250 2i
c(7)
＝0.1179
[0684]
由公式3.5可求得t
27
中
[0685]
同理可得到下列底事件的关键重要度：
[0686]ic(6)
＝0.3216 i
c(17)
＝0.0457 i
c(7)
＝0.1235 i
c(4)
＝0.0424i c(5)
＝0.2621
[0687]
概率重要度矩阵表如表7至表12所示：
[0688]
表7空客a320相似机型绿液压系统概率重要度表
[0689][0690]
表8波音b737ng相似机型b液压系统概率重要度表
[0691]
[0692][0693]
表9空客a320相似机型黄液压系统概率重要度表
[0694][0695]
表10波音b737ng相似机型a液压系统概率重要度表
[0696]
[0697][0698]
表11空客a320相似机型蓝液压系统概率重要度表
[0699][0700]
表12波音b737ng相似机型备用液压系统概率重要度表
[0701][0702]
关键重要度矩阵表如表13至表18所示：
[0703]
表13空客a320相似机型绿液压系统关键重要度表
[0704]
[0705][0706]
表14波音b737ng相似机型b液压系统关键重要度表
[0707][0708]
表15空客a320相似机型黄液压系统关键重要度表
[0709]
[0710][0711]
表16波音b737ng相似机型a液压系统关键重要度表
[0712][0713]
表17空客a320相似机型蓝液压系统关键重要度表
[0714][0715][0716]
表18波音b737ng相似机型备用液压系统关键重要度表
[0717][0718]
将上述表格所得到的数据编写成矩阵，运用matlab进行余弦相似度的编程计算后，可得到六组余弦相似度数据，结果为(cos1为余弦值，cos2为弧度值，v为角度值)：
[0719]
(概率重要度矩阵)
[0720]
％％1空客a320相似机型绿液压系统与波音b737ng相似机型b液压系统
[0721]
cos1＝0.9994 cos2＝0.0340 v＝1.9507；
[0722]
％％2空客a320相似机型黄液压系统与波音b737ng相似机型a液压系统
[0723]
cos1＝0.9841 cos2＝0.1785 v＝10.2262；
[0724]
％％3空客a320相似机型蓝液压系统与波音b737ng相似机型备用液压系统
[0725]
cos1＝0.9961 cos2＝0.0889 v＝5.0921。
[0726]
(关键重要度矩阵)
[0727]
％％4空客a320相似机型绿液压系统与波音b737ng相似机型b液压系统
[0728]
cos1＝0.9998 cos2＝0.0189 v＝1.0809；
[0729]
％％5空客a320相似机型黄液压系统与波音b737ng相似机型a液压系统
[0730]
cos1＝1.0000 cos2＝0.0000 v＝0.0000；
[0731]
％％6空客a320相似机型蓝液压系统与波音b737ng相似机型备用液压系统
[0732]
cos1＝0.9999 cos2＝0.0112 v＝0.6437。
[0733]
由于余弦相值都近似趋近于1，弧度和角度值都较小趋近于0，可得到空客a320相似机型与波音b737ng相似机型液压系统的相似相关性较强。根据由系统可类比到整体机型相似理论，当两类相似机型的系统都存在近似相似相关时，就可得到机型相似相关的结论。
对于空客a320相似机型与波音b737ng相似机型的法规手册及保留故障清单上还是存在差别，但是对于新机型法规性文件基础的的制定仍是建立在根据已有成熟相似机型发动机的进行修改完善后再进行单独化设计，不仅能提高可靠性还节省了在实验上时间和金钱的花费。
[0734]
以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围，其均应涵盖在本发明的权利要求和说明书的范围当中。

技术特征：
1.一种航空器系统定量相似性分析方法，其特征在于，包括以下步骤：基于两种不同型号的航空器的功能系统，分别建立故障树；通过如下方法分别计算两种不同型号的航空器的功能系统的结构重要度、概率重要度以及关键重要度：通过分析故障树的结构上各个事件的重要程度计算结构重要度；通过计算故障树的顶事件的发生概率对底事件发生概率的偏导数，得到概率重要度；通过计算故障树的底事件故障概率的变化率及其引起的顶事件发生概率的变化率之比，得到关键重要度；对所述两种不同型号的航空器的功能系统的概率重要度和关键重要度分别进行余弦相似度计算，得到两个相似度值；基于所述两个相似度值，来确定两种不同型号的航空器的相似度。2.按照权利要求1所述的方法，其特征在于，所述通过分析故障树的结构上各个事件的重要程度计算结构重要度，包括：通过如下公式计算所述结构重要度：式中，式中，为系统的结构重要度，仅取决于第i个部件在系统中所处的位置，与顶事件的发生概率无关，[φ(1
i
,x)-φ(0
i
,x)]为系统中第i个部件由正常状态(0)转为故障状态(1)，其他部件状态不变时，系统的结构函数发生变化。3.按照权利要求2所述的方法，其特征在于，所述通过计算故障树的顶事件的发生概率对底事件发生概率的偏导数，得到概率重要度，包括：所述概率重要度通过如下公式计算：式中，i
p(i)
为第i个基本事件的概率重要度，p(x
i
)为第i个事件发生的概率，p(t)为某个事件在t时刻发生的概率。4.按照权利要求1所述的方法，其特征在于，所述通过计算故障树的底事件故障概率的变化率与它引起的顶事件发生概率的变化率之比，得到关键重要度，包括：所述关键重要度通过如下公式计算：其中，i
c(i)
为第i个基本事件的关键重要度，p(x
i
)为第i个事件发生的概率，p(t)为某个事件在t时刻发生的概率。5.按照权利要求1所述的方法，其特征在于，对所述两种不同型号的航空器的功能系统的概率重要度和关键重要度分别进行余弦相似度计算，得到两个相似度值，包括：将所述两种不同型号的航空器的功能系统的概率重要度分别转换成两个向量，或者将所述两种不同型号的航空器的功能系统的关键重要度分别转换成两个向量；根据如下公式计算得到相似度值：
其中，cos(θ)表示相似度值，向量a表示其中一种型号的航空器的功能系统的概率重要度或关键重要度，向量b表示另一种型号的航空器的功能系统的概率重要度或关键重要度。6.按照权利要求1所述的方法，其特征在于，基于所述两个相似度值，来确定两种不同型号的航空器的相似度，包括：所述相似度是余弦值或弧度或角度；当余弦值越接近1，表明两种不同型号的航空器的相似度越高；当弧度和角度的值越小，表明两种不同型号的航空器的相似度越高。7.一种航空器系统定量相似性分析装置，其特征在于，包括处理器，所述处理器配置为：基于两种不同型号的航空器的功能系统，分别建立故障树；通过如下方法分别计算两种不同型号的航空器的功能系统的结构重要度、概率重要度以及关键重要度：通过分析故障树的结构上各个事件的重要程度计算结构重要度；通过计算故障树的顶事件的发生概率对底事件发生概率的偏导数，得到概率重要度；通过计算故障树的底事件故障概率的变化率及其引起的顶事件发生概率的变化率之比，得到关键重要度；对所述两种不同型号的航空器的功能系统的概率重要度和关键重要度分别进行余弦相似度计算，得到两个相似度值；基于所述两个相似度值，来确定两种不同型号的航空器的相似度。8.按照权利要求7所述的装置，其特征在于，所述处理器还配置为：将所述两种不同型号的航空器的功能系统的概率重要度分别转换成两个向量，或者将所述两种不同型号的航空器的功能系统的关键重要度分别转换成两个向量；根据如下公式计算得到相似度值：其中，cos(θ)表示相似度值，向量a表示其中一种型号的航空器的功能系统的概率重要度或关键重要度，向量b表示另一种型号的航空器的功能系统的概率重要度或关键重要度。9.按照权利要求7所述的装置，其特征在于，所述处理器还配置为：所述相似度是余弦值或弧度或角度；当余弦值越接近1，表明两种不同型号的航空器的相似度越高；当弧度和角度的值越小，表明两种不同型号的航空器的相似度越高。10.一种计算机可读存储介质，其特征在于，其上存储有计算机可读指令，当所述计算机可读指令被计算机的处理器执行时，使计算机执行权利要求1-6中的任一项所述的方法。

技术总结
本发明公开了一种航空器系统定量相似性分析方法、装置及介质，该方法通过对两种型号航空器的功能系统进行相似相关性计算，作为新型航空器、新航空器技术手册等制定上可以根据已有成熟相似机型的系统作为基础进行修改完善的相似性设计的依据。对相似机型的所有系统都进行重要度分析及相似度的相关计算得出定量相似度，有效控制新型航空器的研制风险，降低研制周期和成本。对航空器系统分别进行结构重要度、概率重要度、关键重要度的计算，并且对概率及关键重要度数据矩阵进行余弦相似度计算，最终得出两种机型系统的相似度值。最终得出两种机型系统的相似度值。最终得出两种机型系统的相似度值。


技术研发人员：闫锋 付尧明 林一帆
受保护的技术使用者：中国民用航空飞行学院
技术研发日：2021.12.09
技术公布日：2022/3/8