主管单位:中华人民共和国工业和信息化部
主办单位:西北工业大学  中国航空学会
地       址:西北工业大学友谊校区航空楼
网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

飞机空调系统动态安全分析  PDF

  • 孙晨赫 1
  • 曹湘玲 2
1. 上海民航职业技术学院 民航维修学院, 上海 200237; 2. 上海民航职业技术学院 民航运输学院, 上海 200237

中图分类号: V245.3

最近更新:2023-12-14

DOI:10.16615/j.cnki.1674-8190.2023.06.16

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

目前,飞机空调系统的安全性分析多基于固定故障状态,未考虑故障动态发展与变化,且仅能进行定性分析,因此有必要对飞机空调系统开展精准安全性分析。基于空调系统FMEA中的根故障参数,借助Markov时间序列求解故障横向转移路径,利用联合概率分布及Bayes算法求解故障等级与概率,实现横向定性又定量的安全性分析;运用Monte Carlo仿真求解故障纵向升级路径、故障等级与概率,实现纵向定性又定量的安全性分析。结果表明:动态定性又定量的安全性分析可以精准地判定飞机空调系统的安全性、故障状态下的放行结论及放行要求,具有较好的工程应用价值。

0 引 言

飞机系统安全性的判定是飞机设计、制造、运营等工作开展的基础。基于故障、失效的系统安全性分析,可用于指导运营人在符合适航标准和运行要求的基础上,在某些设备缺失或者不工作的情况下安全放行飞机,提升航空器的利用率。要确定哪些失效、故障的系统或部件满足安全标准,需要在对飞机各个系统和部件进行故障分析、故障预测的基础上,开展安全性判定,并确定其最终故障是否影响飞机的正常放行。飞机空调系统的安全状态对飞机的安全运行有着重要影响,空调系统故障可能导致座舱失去新鲜空气、座舱失压、超压并引发乘员身体不适,甚至触发紧急下降等严重危害。但并非所有的空调系统故障都会导致严重的安全危害,在发生某些故障时,飞机在满足一定放行条件的情况下,依然可以安全运行。

系统安全性可通过故障危害程度、故障可能性及故障发展路径等关键参数进行判定。K.Jenab

1提出了一种基于流程图的方法来计算飞机空调系统失效时间,用于制定故障限制和维护策略;M.Gerdes2提出了一种基于实测数据的飞机空调系统状态预测方法,通过迭代算法,逐点进行预测;王辉3利用时间序列算法,对飞机空调系统故障进行预测研究;方崴源4针对某型飞机空调系统故障,用模拟工作情况的实验排故方法分析解决常见故障;石旭东5基于AMESim-Simulink联合仿真平台,模拟飞行过程中空调系统组件性能动态变化过程,并对典型故障(预冷器泄漏、预冷器空气活门卡死、压力传感器冲击、冲压空气进气作动筒卡死等)进行模拟,分析了空调组件性能的变化过程;Chen Lu6提出了一种基于自适应阈值和参数提取(ATPE)的飞机空调系统故障诊断方法,实现了对飞机空调系统控制元件的故障检测和隔离。

上述研究主要借助了故障树分析(FTA

7-8、功能危险性评估(FHA9、故障模式与影响分析(FMEA10-11、系统安全性评估(SSA12等评估方法,分别给出了故障逻辑关系、功能危害等级、故障模式与故障影响以及安全性结论。不足的是,上述分析方法都是基于某一确定故障,对动态转移和升级的故障无法做出精准判定,而实际运行时,系统与部件的故障状态一直在动态变化。故障的动态变化主要有两个方向:故障转移和故障升级。同时,FMEA、FTA及SSA都是定性的安全性分析方法,仅能判断故障发生等级,无法分析故障发生概率。而在系统实际运行过程中,低概率的严重故障和大概率的轻微故障都足以触发整个系统的故障,因此需要对系统故障进行既定性又定量的安全性分析。

目前对飞机空调系统的安全性分析,主要基于确定的故障、确定的故障升级路径,并未考虑故障转移、故障发展路径的可能变化;同时对故障升级带来的安全性影响评判上仅进行了定性分析,未实现定量分析,即只考虑了故障等级而未对故障的可能性进行准确界定。故障对系统的安全性影响由故障危害等级、故障概率共同决定,因此需要进行定性又定量的飞机空调系统动态安全性判定。

根据空调系统各部件根故障动态变化路径,把空调系统安全性分析分为2个层次:level-1层的故障横向转移安全性分析,level-2层的故障纵向升级安全性分析。本文通过Markov算

13-14分析转移故障的概率及危害等级,使用Monte Carlo15-17仿真分析升级故障的概率及危害等级,给出空调系统在各类故障、失效状态下飞机的放行结果。

1 安全性适航要求

根据飞机设计及适航要求,飞机系统及部件对飞机安全性的影响需要按照故障及危害等级不同而满足一定的适航安全标准参数ϕ(i),这个条件也称为飞机安全性门限

18,飞机设计中所有系统以及部件都要满足如下门限值:

ϕ(i)1×10-9             (类故)1×10-7             (类故)1×10-5            (类故)1×10-3             (类故)no limitation     (类故) (1)

式中:ϕ(i)为满足飞机安全性的最低要求标准。

各飞机设备和系统必须明确证明能满足ϕ(i),并尽可能达到更高的安全性参数标准。

2 转移故障安全性

飞机的故障状况随着时间在动态地变化,飞机系统和部件在某一时刻的状态仅仅取决于前一时刻的状态,这种性质称为马尔可夫特性。要精准地分析飞机转移到的另一个根故障的安全性,必须以动态变化的时间序列(即马尔可夫序列)下的根故障作为分析对象。以流量控制活门(FCV)为例,对于任意非负整数k组成的马尔可夫时间序列tk,k0,在k时刻FCV的故障状态为si,故障si发生的概率为pi(tk)i=1,2,…,M),其中M是部件FCV所有的故障状态数;在k+1时刻FCV的故障状态为sj,故障sj发生的概率为pj(tk+1)。从根故障si转移到根故障sj的故障概率为pij(k)。根据马尔可夫性质可以得到:

p(ijk)=p(tk+1=sj|tk=si)0 (2)

式中:pij(k)为马尔可夫链的一步转移概率,也即FCV从si状态转移到sj的概率。

根据FMEA分析数据,可得FCV根故障如表1所示。

表1  FCV的根故障
Table 1  FCV root fault

故障

编码

故障状

态编号

故障状态

故障

等级

故障

提示

故障是

否升级

1 s1 失效在关位-A
2 s2 失效在开位-B
3 s3 内部泄漏-C
4 s4 外部泄漏-D
5 s5 丧失控制功能-E
6 s6 示值偏高-F
7 s7 示值偏低-G

不同时刻,不同初始根故障,下一个时刻转移到的另一个根故障及概率各不相同。综合FCV所有可能触发的根故障,可以得到FCV从tk时刻到tk+1时刻根故障横向转移矩阵,如表2所示。

表2  FCV根故障转移矩阵
Table 2  FCV root fault state transition matrix
列号行号
12345678
1 p11 p12 p13 p14 p15 p16 p17 j=17p1j
2 p21 p22 p23 p24 p25 p26 p27 j=17p2j
3 p31 p32 p33 p34 p35 p36 p37 j=17p3j
4 p41 p42 p43 p44 p45 p46 p47 j=17p4j
5 p51 p52 p53 p54 p55 p56 p57 j=17p5j
6 p61 p62 p63 p64 p65 p66 p67 j=17p6j
7 p71 p72 p73 p74 p75 p76 p77 j=17p7j

根据FMEA分析报告数据,FCV根故障概率pi={p1,p2,p3,p4,p5,p6,p7}=           10.03×10-6,9.78×10-6,19.36×10-6,

0.52×10-6,24.17×10-6,0.24×10-6,
29.76×10-6

以FCV的第4个根故障为例,FCV 外部泄漏的最大故障转移概率为:max p41,p42,p43,p45,p46,p47。根据式(2)可得,在任意时刻p41=p(tk+1=s1|tk=s4)0p41仅受故障状态s1s4影响,通过贝叶斯(Bayes

19-20算法可得:

p41=p(s1|s4)=p(s4|s1)×p(s4)p(s1) (3)

设有N条故障路径会触发FCV第1个故障,p41仅是其中一条,则可得:

p41<j=1Npj1=p(s1) (4)

表明在出现故障4的基础上再出现故障1的概率明显小于故障1的触发概率。把FCV根故障概率参数带入式(4),可得p41<10-5。综上同理可得:

max {p41,p42,p43,p45,p46,p47}<10-5 (5)

故障4转移到故障1的故障分布,是两个随机变量p(s4)p(s1)组成的随机向量的概率分布,属于联合概率分布,可表示为p(s4,s1),根据联合概率与贝叶斯条件概率关系,可得:

p(s4,s1)=p(s1|s4)×p(s4) (6)

式(4)带入式(6),可得:

p(s4,s1)<p(s1)×p(s4)=5×10-12 (7)

式(7)带入式(1),可得部件FCV的4号根故障转移到1号根故障满足安全性标准要求。 重复上述步骤,可得FCV根故障之间状态转移都满足ϕ(i)限定。

重复以上算法计算FCV的1~3号、5~7号根故障的横向转移故障发生概率:

pij=0,  pi×pj            (ij)1-pi×pj,  1    (i=j) (8)

综上,不考虑故障升级的情况下,FCV根故障的横向转移故障概率满足适航安全门限ϕ(i)要求。

FCV不同根故障间的故障转移属于同一个部件之间的故障状态动态变化,在某些情况下,不同的部件之间也存在故障转移的可能。故障转移的本

21是物质、能量、信息之间的交换,具体的转移过程如图1所示。

图1  故障转移的过程

Fig.1  The process of failover

故障要发生转移,在物质、能量、信息这3个要素中至少要有1种存在交换。对于典型的机电类设备——飞机空调系统部件,不存在物质交换基础和信息传输电路,因此不同部件间的故障转移仅由能量传输导致。能量传输有传导、对流和辐射3种模式,空调系统部件不具备对流和辐射的条件,只能通过传导的形式传输能量。能量传导是一种通过物理接触,实现能量交换的现象,空调系统部件互相独立,无物理接触。据此可建立空调部件间故障转移判断矩阵,如表3所示。

表3  部件间故障转移判断矩阵
Table 3  Judgment matrix of inter component fail-over
部件编号部件名称传输路径是否转移
物质信息能量
1 流量控制活门(FCV)
2 应急冲压空气作动器(ERAIA)
3 地面接头(LPGC)
4 货舱关断活门(BBSOV)
5 再循环风扇(RFAN)
6 电子设备风扇(AVS FAN)
7 电子设备通风活门(AVV)
8 配平空气活门(TAV)
9 配平空气压力调节活门(TAPRV)
10 飞行员脚部加热器

基于马尔可夫时间序列,对空调系统的其余部件故障,按照上述分析方法,逐项进行分析,并带入式(1)进行验证,可得到空调系统level-1层横向转移故障安全性分析结论,如表4所示。

表4  空调系统横向故障转移安全分析结论
Table 4  Conclusion of fail-over safety analysis of air conditioning system

部件

编号

部件名称

故障

等级

故障

提示

故障转移

安全性

故障是

否升级

1 流量控制活门(FCV)
2 应急冲压空气作动器(ERAIA)
3 地面接头(LPGC)
4 货舱关断活门(BBSOV)
5 再循环风扇(RFAN)
6 电子设备风扇(AVS FAN)
7 电子设备通风活门(AVV)
8 配平空气活门(TAV)
9 配平空气压力调节活门(TAPRV)
10 飞行员脚部加热器

3 升级故障安全性

对于发生故障升级的根故障而言,通过马尔可夫分析可以得到其转移到另一个根故障的概率,但无法描述根故障与下一个升级故障之间的逻辑关系以及故障转移路径,故无法确定当前根故障对飞机造成的最终危害。为了确定根故障对飞机的最终危害,通常通过遍历故障树分析

22进行分析,但是遍历分析都是基于某一特定故障或某条特定故障升级路径,当故障动态变化,故障升级路径交叉多变时,常规FTA故障分析法难以为继。

同时,常规FTA方法分析升级故障时,无法得到升级故障准确的触发概率,因而无法对升级故障做出准确的安全性判定。为了得到准确的升级故障概率,需要对全部根故障在其故障发生周期内的工作状态进行分析验证,分析验证的次数越多,验证结果越接近真实值。传统的数理方法无法解析这么多次且如此长时间序列的验证问题,为此,应用Monte Carlo仿真分析法对空调系统部件的升级故障进行全时段、全故障状态仿真分析。飞机空调系统具有通风、温控和流量调节功能。空调系统部件与功能之间的逻辑关系如图2所示。

图2  空调系统部件与功能之间的逻辑关系

Fig.2  Logical relationship between components and functions of an air conditioning system

3.1 Monte Carlo分析算法

设部件的失效分布函

23-24fij(t)λij 为失效分布的参数,i为故障部件编号,i1,2,3,4,5,6,7,8,9j为部件i的故障序号,根据FMEA分析数据,部件失效分布参数 如表5所示。

表5  空调系统部件失效分布参数表
Table 5  Component failure distribution parameters of the air conditioner
部件编号部件失效概率密度分布fij(t)fij(t)的分布参数(1/λ
1 FCV Index 0.99, 1.02, 0.52, 19.23, 0.41, 41.67, 0.34×105 h
2 ERAIA Index 3.33, 3.33, 3.33, 3.33×105 h
3 LPGC Index 10, 2, 0.14×106 h
4 BBSOV Index 1.67, 1.67, 25.33, 25.33, 63.25×105 h
5 RFAN Index 0.31, 1.13, 80.16, 77.32, 93.25×105 h
6 AVS FAN Index 0.43, 1.27, 90.31, 80.32, 106.47, 68.05, 76.32}×105 h
7 AVV Index 0.87, 32.67, 53.62, 1.43,  26.56, 41.29×105 h
8 TAV Index 0.53, 143.12, 38.76, 1.86, 37.26, 67.51×105 h
9 TAPRV Index 0.87, 4.37, 1.21, 5.32, 1.46, 0.92, 5.76, 0.34}×105 h

仿真算法及逻辑:

1) 为了模拟飞机生命周期的故障状

25,设定最大工作时间tmax=106 h,仿真的推进步长为

3 h,时间间隔为1 h。

2) 设定的仿真样本个

26Ns=2.5×105

3) 仿真运算逻辑:根据FMEA中的根故障状态,基于FTA升级路径,分别建立部件1~9仿真逻辑

27与根故障代码,其中底层是根故障,终点是根故障的下一个升级故障。

4) 故障结构函数为φx(t),其中x是部件的故障状态,若部件的故障能够向上一个层级升级,则此时y=φx(t)=1,代表系统故障,一次仿真结束,对应的部件故障事件计数加1;若部件故障无升级路径,则退出本次仿真,并开始下一次故障仿真。

5) 在第N次仿真运行中,第i个部件故障zij的失效时间tijN由失效分布密度函数fij(t)抽样产生,因而在第N次仿真中9个部件的失效时间抽样值为t1jN, t2jN, t3jN, t4jN, t5jN, t6jN, t7jN, t8jN, t9jN

6) 在每次仿真运行中,按照仿真逻辑树逻辑,利用结构函数来判断升级故障发生时间tijN

7) 重复上述过程直到仿真运行Ns次。

8) 借助Python语言,其转化逻辑如图3所示。

图3  Python转化逻辑

Fig.3  Conversion logical process by Python

3.2 仿真故障树

根故障代码及状态如表6所示,仿真逻辑树如图4所示。

表6  根故障代码及状态
Table 6  Root fault codes and status
故障代码故障状态故障代码故障状态故障代码故障状态
FCV-1 失效在开位 BBSOV-4 内部泄漏 AVV-1 失效在开位
FCV-2 失效在关位 BBSOV-5 位置信息失效 AVV-2 失效在关位
FCV-3 内部泄漏 RFAN-1 未检测到的失效 AVV-3 显著的外部泄漏
FCV-4 示值偏高 RFAN-2 鼓风性能下降 AVV-4 内部泄漏
FCV-5 外部泄漏 RFAN-3 单向活门失效在开位 AVV-5 位置信息失效,全关位置微调失效
FCV-6 丧失控制能力 RFAN-4 电子保护失效 AVV-6 位置信息失效,未全关位置微调失效
FCV-7 示值偏低 RFAN-5 单向活门失效在关位 TAV-1 显著外部泄漏
ERAIA-1 失效在关位 RFAN-6 鼓风失效 TAV-2 失效在关位
ERAIA-2 失效在开位 AVSFAN-1 鼓风性能下降 TAV-3 失效在开位
ERAIA-3 位置指示失效 AVSFAN-2 单向活门失效 TAV-4 内部泄漏
LGPC-1 显著泄漏 AVSFAN-3 电子保护失效 TAV-5 位置信息失效
LGPC-2 失效在关位 AVSFAN-4 电子保护下降 TAPRV-1 调压偏低
LGPC-3 轻微泄漏 AVSFAN-5 过热保护丧失 TAPRV-2 调压偏高
BBSOV-1 失效在开位 AVSFAN-6 机械保护丧失 TAPRV-3 显著的外部泄漏
BBSOV-2 失效在关位 AVSFAN-7 鼓风失效 TAPRV-4 失效在开位
BBSOV-3 显著外部泄漏 AVSFAN-8 排风失效 TAPRV-5 失效在关位

(a)  流量控制活门故障逻辑

(b)  应急冲压空气作动器故障逻辑

(c)  地面接头故障逻辑

(d)  货舱关断活门故障逻辑

(e)  再循环风扇故障逻辑

(f)  电子设备风扇故障逻辑

(g)  电子设备通风活门故障逻辑

(h)  配平空气活门故障逻辑

(i)  配平空气压力调节活门故障逻辑

图4  仿真逻辑树

Fig.4  Simulation logic tree

3.3 结果分析

表5中的各部件根故障参数输入训练好的基于Python环境的Monte Carlo仿真环境进行故障分析,运行程序后,得到空调系统1~9号部件的升级故障概率pi(t),如表7所示。带入式(1)可知其满足适航安全标准门限ϕ(i),表明空调系统1~9号部件故障的危害程度及触发频次处于可接受范围。

表7  空调系统部件最终故障
Table 7  Final failure of air conditioning system components
部件最终升级故障最终升级故障发生概率pi(t)等级
FCV 丧失单侧空调系统气源 4.32×10-5
ERAIA 丧失应急通风功能 1.12×10-5
LPGC 丧失单侧空调系统气源 21.83×10-6
BBSOV 丧失货舱通风和货舱隔离功能 5.21×10-5
RFAN 丧失空调系统冗余功能 31.87×10-5
AVS FAN 电子设备舱通风冗余丧失 56.13×10-5
AVV 丧失电子设备通风功能的冗余 71.02×10-6
TAV 温度控制系统失效 26.39×10-6
TAPRV 丧失单侧空调系统气源 13.24×10-6
配平空气失效 41.57×10-6

虽然空调系统1~9号部件发生故障后,还能满足适航安全要求,但是会对飞机的正常运行造成影响,需要增加额外的运行限制,是否放行飞机还需按照运行要求进行综合考虑。例如,1套FCV失效会导致对应侧空调组件(PACK)失效,进而导致飞机的理论升限降至31 000 ft(1 ft=0.304 8 m)以下。此外,相同部件,故障数量不同也会产生不同的影响。例如,单TAPRV失效时,空调系统会丧失单侧气源,温度控制功能降级,导致系统冗余度降低;而双TAPRV失效则会导致整个空调系统失效。 最后还需考虑故障对飞机特定状态,特别是应急状态的运行影响。因此结合CCAR-121部运行规

28及实际运行要求,可得空调系统故障放行结论及放行要求,如表8所示。在飞机维护中,若发现表8中不可放行的故障,需要立即安排维修;可放行的故障在一定的期29内可以延后修理。

表8  空调系统故障放行结论及放行条件
Table 8  Conclusions and requirements for air conditioning system failure release
部件安装数量失效数量是否存在特殊运行要求是否放行放行条件/拒绝放行原因
FCV 2 1 运行高度小于31 000 ft
2 座舱无法增压
ERAIA 1 1 影响应急通风
LPGC 1 1 不可使用地面空调车
BBSOV 2 1 货舱不可装载活物
2 货舱不可装载活物
RFAN 2 1 空调系统丧失冗余度
2 空调系统丧失冗余度
AVS FAN 2 1

另外一套工作正常

AVV不能坏

2 电子舱过热
AVV 1 1 影响应急条件下电子舱通风
TAV 3 1 丧失温度控制
2 丧失温度控制
3 丧失温度控制
TAPRV 2 1 运行高度小于31 000 ft
2 座舱无法增压

4 结 论

1) 本文所运用的空调系统动态安全性分析方法,不基于固定故障,充分考虑了故障的动态发展与变化中的转移、升级两个维度,更加符合实际故障状态。

2) 实现了定性又定量的空调系统安全性判定,既判定了故障的危害等级,又求解了故障发生的概率,提升了安全性分析的精准度。

3) 通过动态安全性分析,得到了空调系统故障条件下飞机的放行结论及放行条件。飞机放行与故障数量、特殊的运行要求相关。不可放行的故障需要立即维修,可放行的故障需按照规定期限完成维修。

参考文献

1

JENAB KRASHIDI K. Operational reliability assessment of an aircraft environmental control system[J]. Reliability Engineering and System Safety2008942): 456-462. [百度学术] 

2

GERDES M. Decision trees and genetic algorithms for condition monitoring forecasting of aircraft air conditioning[J]. Expert Systems with Applications20134012): 5021-5026. [百度学术] 

3

王辉齐亚东. 飞机空调系统频发故障模型构建及预测[J]. 航空维修与工程20173): 39-43. [百度学术] 

WANG HuiQI Yadong. Modeling and forecasting frequent fault of aircraft air conditioning system[J]. Aviation Maintenance and Engineering20173): 39-43.(in Chinese) [百度学术] 

4

方崴源谌凯欣李伟. 某型飞机空调系统故障的实验排除分析[J]. 装备制造技术20155): 146-147. [百度学术] 

FANG WeiyuanCHEN KaixinLI Weiet al. Experimental analysis of air conditioning system failure of an aircraft[J]. Equipment Manufacturing Technology20155): 146-147.(in Chinese) [百度学术] 

5

石旭东蒋贵嘉张宇. 基于联合仿真的飞机空调系统故障影响[J]. 航空学报2020418): 295-303. [百度学术] 

SHI XudongJIANG GuijiaZHANG Yuet al. Fault influence of aircraft air conditioning system based on joint simulation[J]. Acta Aeronautica et Astronautica Sinica2020418): 295-303.(in Chinese) [百度学术] 

6

CHEN LuCHENG YujieLIU Hongmeiet al. An approach to fault detection and isolation for control components in the aircraft environmental control system[J]. Proceedings of the Institution of Mechanical Engineering20142287): 1202-1214. [百度学术] 

7

陈洪转赵爱佳李腾蛟. 基于故障树的复杂装备模糊贝叶斯网络推理故障诊断[J]. 系统工程与电子技术2021435): 1250-1251. [百度学术] 

CHEN HongzhuanZHAO AijiaLI Tengjiaoet al. Fuzzy Bayesian network inference for fault diagnosis of complex equipment based on fault tree[J]. System Engineering and Electronics2021435): 1250-1251.(in Chinese) [百度学术] 

8

李小彪马庆岩. 航空发动机高压转子转速摆动故障分析与排除[J]. 航空发动机20214710): 63-64. [百度学术] 

LI XiaobiaoMA Qingyan. Analysis and elimination of high pressure rotor speed swing fault of aero-engine[J]. Aeroengine20214710): 63-64.(in Chinese) [百度学术] 

9

肖女娥闫芳王鹏. 基于安全论证的民机机载系统安全性评估[J]. 中国安全科学学报20192912): 73-75. [百度学术] 

XIAO Nyu’eYAN FangWANG Peng. Safety assessment of civil aircraft airborne system based on safety demonstration[J]. China Safety Science Journal20192912): 73-75.(in Chinese) [百度学术] 

10

李冰月孙建红. 基于FMEA的飞机空调系统故障诊断与仿真[J]. 振动、测试与诊断2017373): 588-590. [百度学术] 

LI BingyueSUN Jianhong. Fault diagnosis and simulation of aircraft air conditioning system based on FMEA[J]. Vibration, Testing and Diagnosis2017373): 588-590.(in Chinese) [百度学术] 

11

李贺. FMEA在民用飞机研制中的应用研究[J]. 质量与可靠性20211): 11-15. [百度学术] 

LI He. Application of FMEA in civil aircraft development[J]. Quality and Reliability20211): 11-15.(in Chinese) [百度学术] 

12

王祥. 面向持续适航的民机系统安全性分析技术研究[D]. 南京南京航空航天大学2017. [百度学术] 

WANG Xiang. Research on safety analysis technology of civil aircraft system for continuous airworthiness[D]. NanjingNanjing University of Aeronautics and Astronautics2017.(in Chinese) [百度学术] 

13

陶长琪. 决策理论与方法[M]. 北京中国人民大学出版社2010160-164. [百度学术] 

TAO Changqi. Decision theory and method[M]. BeijingChina Renmin University Press2010160-164.(in Chinese) [百度学术] 

14

陈雯吕王勇李思奇. 两种马尔可夫链状态转移概率矩阵的估计与比较[J]. 重庆理工大学学报(自然科学版)2021358): 218-220. [百度学术] 

CHEN WenLYU WangyongLI Siqiet al. Estimation and comparison of state transition probability matrices for two Markov chains[J]. Journal of Chongqing University of Technology (Natural Science)2021358): 218-220.(in Chinese) [百度学术] 

15

FAN YongleiCHEN Guoguang. Application of Monte Carlo method in rudder control precision[J]. Journal of Measurement Science and Instrumentation201564): 378-379. [百度学术] 

16

ESTÉBANEZ L RSHEHAB ESYDOR Pet al. An integrated aerospace requirement setting and risk analysis tool for life cycle cost reduction and system design improvement[J]. Procedia CIRP (The International Academy for Production Engineering)20153878-83. [百度学术] 

17

WENSKY TWINKLER LFRIEDRICHS J. Environmental influences on engine performance degradation: GT2010-22748[R]. UKTurbo Expo, Power for Land, Sea, and Air2010. [百度学术] 

18

中国民用航空局. 运输类飞机适航标准: CCAR-25-1309[S]. 北京中国民用航空局2011. [百度学术] 

Civil Aviation Administration of China. Airworthiness standard for transport aircraft: CCAR-25-1309[S]. BeijingCAAC2011.(in Chinese) [百度学术] 

19

YUN TaoliXIAO NingheJIAN Shuai. Risk analysis and maintenance decision making of natural gas pipelines with external corrosion based on Bayesian network[J]. Petroleum Science2022193): 1251-1253. [百度学术] 

20

WANG LihuaCUI YahuiZHANG Fengqiet al. Stochastic speed prediction for connected vehicles using improved Bayesian networks with back propagation[J]. Science China(Technological Sciences)2022657): 1525-1527. [百度学术] 

21

柴毅张可毛永芳. 动态系统运行安全性分析与技术[M]. 北京化学工业出版社201948-50. [百度学术] 

CHAI YiZHANG KeMAO Yongfanget al. Analysis and technology of dynamic system operation safety[M]. BeijingChemical Industry Press201948-50.(in Chinese) [百度学术] 

22

李靖杨帆. 分布式协同多机器人多任务目标遍历路径规划[J]. 天津工业大学学报2020126): 70-74. [百度学术] 

LI JingYANG Fan. Distributed collaborative multi-robot multi-task target traversal path planning[J]. Journal of Tianjin Polytechnic University2020126): 70-74.(in Chinese) [百度学术] 

23

丁勇. 独立同均匀分布随机变量和的性质[J]. 统计与决策20213717): 71-73. [百度学术] 

DING Yong. Properties of sum of independent and uniformly distributed random variables[J]. Statistics & Decision20213717): 71-73.(in Chinese) [百度学术] 

24

张斯龙何世熠陈亮. 基于概率密度函数的加权变换在能谱测量中的应用[J]. 现代应用物理2022131): 45-51. [百度学术] 

ZHANG SilongHE ShiyiCHEN Lianget al. Application of weighting transformation based on probability density function in energy spectrum measurement[J]. Modern Applied Physics2022131): 45-51.(in Chinese) [百度学术] 

25

崔铁军李莎莎. 系统故障演化过程最终事件状态及发生概率研究[J]. 中国安全科学学报2021318): 2-6. [百度学术] 

CUI TiejunLI Shasha. Research on the final event state and occurrence probability of system failure evolution process[J]. China Safety Science Journal2021318): 2-6.(in Chinese) [百度学术] 

26

陈敏捷羊荣金. 液压支架液压系统故障仿真与诊断技术研究[J]. 液压与气动20218): 128-131. [百度学术] 

CHEN MinjieYANG Rongjin. Research on fault simulation and diagnosis technology of hydraulic system of hydraulic support[J]. Hydraulic and pneumatic20218): 128-131.(in Chinese) [百度学术] 

27

张颖董国强. 基于DEMATEL与逻辑树的海底原油管道泄漏风险分析[J]. 安全与环境学报20202010): 1692-1693. [百度学术] 

ZHANG YingDONG Guoqiang. Leakage risk analysis of submarine crude oil pipeline based on DEMATEL and logic tree[J]. Journal of Safety and Environment20202010): 1692-1693.(in Chinese) [百度学术] 

28

中国民用航空局. 大型飞机公共航空运输承运人运行合格审定规则: CCAR-121.329-331[S]. 北京中国民用航空局2017. [百度学术] 

Civil Aviation Administration of China. Operation qualification certification rules for large aircraft public air transport carriers: CCAR-121.329-331[S]. BeijingCAAC2017.(in Chinese) [百度学术] 

29

中国民用航空局. 航空器的运行文件: AC-91-24.6[S]. 北京中国民用航空局2014. [百度学术] 

Civil Aviation Administration of China. Operating documents of aircraft: AC-91-24.6[S]. BeijingCAAC2014.(in Chinese) [百度学术]