针对多失效模式关联系统的安全性概率风险评价方法技术方案

本发明专利技术属于系统安全性技术研究领域，涉及一种利用集合理论针对多失效模式关联系统的安全性概率风险评价方法。本发明专利技术首先针对复杂系统多失效模式关联特点，利用集理论将系统的失效模式通过多个时变极限状态函数对应的失效域来表示，解决了多失效模式下系统安全性建模难题；再针对失效模式关联系统退化特性，考虑时间因素，得到多模式多失效模式下系统安全性模型随时间变化情况；在得到多失效模式对应的累积失效域基础上，利用蒙特卡洛方法进行安全性风险概率分析评价，解决了多失效模式相关且具有时变特性的系统安全性定量评价难题。

【技术实现步骤摘要】
针对多失效模式关联系统的安全性概率风险评价方法
本专利技术属于系统安全性技术研究领域，涉及一种利用集合理论针对多失效模式关联系统的安全性概率风险评价方法。
技术介绍
随着航空武器装备研制水平提高，现代航空装备系统具有两个特点：一个特点是系统与飞机的多个功能相关或对多个功能有影响，另一个特点是系统的状态复杂，不能通过试验独立证明其安全性，并且在没有分析工具的帮助下难于理解其逻辑性，因此复杂系统安全性分析成为目前研究难点。国外相关系统安全性标准和规范对复杂系统进行了具体定义，同时要求考虑多重失效对系统安全性的影响。ARP4754及ARP4761将复杂系统定义为：系统或者产品的一种使其运行难于理解的属性，系统的这种属性的增高通常是由复杂的部件、多重相互关系等因素导致。ARP4761中通过共模分析(CMA,CommonModeFailure)校验故障树分析中“与门”事件的独立性，FAR25.1309要求通过分析手段来验证飞机系统安全性定量要求(包括灾难的、危险的、主要的和次要的)，其中包括对多重失效的分析与评估。因此，对复杂系统而言，由于影响系统安全性的多个失效状态之间具有相关性，“相关”是其失效的普遍特征，忽略系统失效的相关性，简单地在系统各部分失效状态相互独立的假设下进行系统安全性分析与评估，常常会导致过大的误差。在相关标准和规范明确要求在飞机系统安全性评估过程中对共因失效进行定性和定量分析，并且要验证失效状态之间的独立性，综上所述，如何针对多失效模式关联的复杂系统，进行安全性概率风险评价是目前研究难点。
技术实现思路
本专利技术的目的是提出针对多失效模式关联系统的安全性概率风险评价方法，该方法突破传统的系统安全性评估中失效模式独立性假设，考虑的复杂系统失效模式的相关性和时变性，构建多模式关联系统安全性分析模型，并进行定量的安全性概率风险评价，解决传统失效模式独立的假设下系统安全性分析与评估中误差过大难题。本专利技术的技术解决方案是，步骤一：识别影响模式关联系统安全性的关联失效模式，步骤二：定义模式关联系统安全性概率风险评价模型，假设模式关联系统存在m个失效模式，则对应m个时变极限状态函数，第i个时变极限状态函数定义为：i＝1,2,…,m，t∈[0,T]。其中和Di(.)为第i个极限状态函数的广义强度和广义应力，ω表示样本空间；步骤三：定义步骤二中随机向量X(t,ω)随机特性，包括分布类型、均值和标准差；步骤四：确定安全性概率风险评价模型中的失效域和安全域；根据步骤二定义的时变极限状态函数，当表示安全域，表示极限状态，G(t,X(t,ω))＜0表示失效域；步骤五：定义单一失效事件集合和危险概率；使用期，模式关联系统单个失效表示为：累积失效概率为在时刻ti的瞬时失效概率定义为pf,i(ti)＝P{G(ti,X(ti,ω))≤0}，其中，P(.)表示概率函数，下标f表示失效；步骤六：连续时间离散化；将连续的时间区间[0,T]离散化，设时间步长为h，则离散的时间区间数量为L＝T/h，第l个时间步为tl＝l×h，其中l＝0,1,…,L,T为系统使用周期。步骤七：定义第i个失效模式对应的失效域；设模式关联系统第i个时变极限状态函数为Gi(t,X(t,ω))，则系统在tl时刻的第i个失效模式对应的失效域定义为El,i＝{x∈X,Gi(tl,x(tl,ω))≤0}；步骤八：建立单个并联子系统对应的瞬时失效域；假设系统由并联和串联子系统构成，则由Nj个零部件组成的第j个并联子系统在tl时刻的瞬时失效域为：步骤九：建立单个串联子系统对应的瞬时失效域；同理步骤八，由Sj个零部件组成的第j个串联子系统在tl时刻的瞬时失效域为步骤十：建立若干个子系统构成整个失效模式关联系统瞬时失效域根据模式关联系统的构成方式，采取以下两种方法之一：1）建立若干个并联子系统以串联的形式构成整个系统瞬时失效域；由n个并联子系统以串联的形式构成的系统，在tl时刻的瞬时失效域为：2）建立若干个串联子系统以并联的形式构成整个系统瞬时失效域。由s个串联子系统以并联的形式构成的系统，在tl时刻的瞬时失效域为：步骤十一：确定系统的累积失效域，通过步骤十确定了系统在tl时刻的瞬时失效域，在此基础上，得到系统在[0,tl]内的累积失效域：步骤十二：确定系统离散时间内增加的失效域；系统在离散时间区间[tl,tl+h]内增加的失效域表示为：通过布尔法则进一步简化为步骤十三：通过离散时间内的增长失效域表示系统累积失效域；由步骤十一和十二，系统的累积失效域通过各个离散时间区间的增长失效域来表示，即步骤十四：确定在时域内多失效模式下的危险概率对多失效模式的系统在使用期[0,T]内危险概率表示为：步骤十五：定义子系统瞬时失效和安全状态函数；子系统在tl时刻的失效和安全状态定义为步骤十六：在瞬时时刻进行蒙特卡洛抽样；利用蒙特卡洛抽取N个随机样本，若用Xk表示第k组随机样本，则在tl时刻系统失效的次数统计为：步骤十七：在时间区间内进行蒙特卡洛抽样；在[tl,tl+1]区间内系统失效次数统计为：其中步骤十八：利用蒙特卡洛计算离散事件区间内的累积失效概率。在离散时间区间[tl,tl+1]增加失效域对应的失效概率表达为：步骤十九：将步骤十八的失效概率式子带入步骤十四，即可得到系统在使用期内累积失效概率，从而完成对多失效模式关联复杂系统失效概率仿真分析。本专利技术的优点和积极效果在于：(1)针对复杂系统多失效模式关联特点，利用集理论将系统的失效模式通过多个时变极限状态函数对应的失效域来表示，解决了多失效模式下系统安全性建模难题；(2)针对失效模式关联系统退化特性，考虑时间因素，得到多模式多失效模式下系统安全性模型随时间变化情况；(3)在得到多失效模式对应的累积失效域基础上，利用蒙特卡洛方法进行安全性风险概率分析评价，解决了多失效模式相关且具有时变特性的系统安全性定量评价难题。附图说明图1为本专利技术的多失效模式关联系统安全性概率风险评价方法的整体步骤流程图；图2简单并联系统失效域随时间变化情况，图3两个并联子系统以串联形式构成的整个系统失效域随时间变化情况，图4某型无人机机翼前梁、中梁和后梁结构以及应力云图，图5多失效模式关联系统安全性概率风险评价结果。具体实施方式下面结合附图和实施例子对本专利技术做进一步说明，整体实施流程如图1所示，具体实施步骤如下所述：步骤一：识别影响系统安全性的关联失效模式。通过系统安全性分析，如FTA(FaultTreeAnalysis，故障树分析)，识别FTA底事件“与门”中关联失效模式，并通过FTA将系统安全性概率要求分解至底事件。步骤二：定义系统安全性概率风险评价模型。本专利技术中通过时变极限状态函数Gi(t,X(t,ω))构建复杂系统安全性概率风险评价模型，对于多功能的复杂系统，假设系统存在m个失效模式，则对应m个时变极限状态函数，第i个时变极限状态函数定义为：Gi(t,X(t,ω))＝Di-Di(X(t,ω))，i＝1,2,…,m，t∈[0,T]。其中和Di(.)为第i个极限状态函数的广义强度和广义应力，X(t,ω)为随机向量，ω表示样本空间。步骤三：安全性概率风险评价模型中参数随机特性。在安全性概率风险评价模型中，定义p+q维随机向量X(t,ω)：X(t,ω)＝[X1本文档来自技高网...

【技术保护点】
针对多失效模式关联系统的安全性概率风险评价方法，其特征是，步骤一：识别影响模式关联系统安全性的关联失效模式，步骤二：定义模式关联系统安全性概率风险评价模型，假设模式关联系统存在m个失效模式，则对应m个时变极限状态函数，第i个时变极限状态函数定义为：i＝1,2,…,m，t∈[0,T]。其中和Di(.)为第i个极限状态函数的广义强度和广义应力，ω表示样本空间；步骤三：定义步骤二中随机向量X(t,ω)随机特性，包括分布类型、均值和标准差；步骤四：确定安全性概率风险评价模型中的失效域和安全域；根据步骤二定义的时变极限状态函数，当表示安全域，表示极限状态，G(t,X(t,ω))＜0表示失效域；步骤五：定义单一失效事件集合和危险概率；使用期，模式关联系统单个失效表示为：E={G(τ,X(τ,ω))≤0,∃τ∈[0,T]}累积失效概率为pf,c(0,T)=P(E)=P{G(τ,X(τ,ω))≤0,∃τ∈[0,T]};在时刻ti的瞬时失效概率定义为pf,i(ti)＝P{G(ti,X(ti,ω))≤0}，其中，P(.)表示概率函数，下标f表示失效；步骤六：连续时间离散化；将连续的时间区间[0,T]离散化，设时间步长为h，则离散的时间区间数量为L＝T/h，第l个时间步为tl＝l×h，其中l＝0,1,…,L,T为系统使用周期；步骤七：定义第i个失效模式对应的失效域；设模式关联系统第i个时变极限状态函数为Gi(t,X(t,ω))，则系统在tl时刻的第i个失效模式对应的失效域定义为El,i＝{x∈X,Gi(tl,x(tl,ω))≤0}；步骤八：建立单个并联子系统对应的瞬时失效域；假设系统由并联和串联子系统构成，则由Nj个零部件组成的第j个并联子系统在tl时刻的瞬时失效域为：步骤九：建立单个串联子系统对应的瞬时失效域；同理步骤八，由Sj个零部件组成的第j个串联子系统在tl时刻的瞬时失效域为Ll,j=∪i=1sjEl,i;步骤十：建立若干个子系统构成整个失效模式关联系统瞬时失效域根据模式关联系统的构成方式，采取以下两种方法之一：1）建立若干个并联子系统以串联的形式构成整个系统瞬时失效域；由n个并联子系统以串联的形式构成的系统，在tl时刻的瞬时失效域为：El=Ll,1∪Ll,2∪···∪Ll,n=∪j=1nLl,j---(1)2）建立若干个串联子系统以并联的形式构成整个系统瞬时失效域。由s个串联子系统以并联的形式构成的系统，在tl时刻的瞬时失效域为：El=Ml,1∩Ml,2∩···∩Ml,n=∩j=1sMl,j---(2)步骤十一：确定系统的累积失效域，通过步骤十确定了系统在tl时刻的瞬时失效域，在此基础上，得到系统在[0,tl]内的累积失效域：Al=E0∪E1∪···∪El=∪q=0lEq;步骤十二：确定系统离散时间内增加的失效域；系统在离散时间区间[tl,tl+h]内增加的失效域表示为：Bl=Al+1∩A‾l=(∪q=0l+1Eq)∩(∩q=0lE‾q)---(3)通过布尔法则进一步简化为Bl=El+1∩q=0lE‾q=El+1∩A‾l;步骤十三：通过离散时间内的增长失效域表示系统累积失效域；由步骤十一和十二，系统的累积失效域通过各个离散时间区间的增长失效域来表示，即Al+1=∪q=0l+1Eq=E0∪(E1∩E‾0)∪(E1∩(E‾1∩E‾0))∪···∪(El+1∩A‾l)=E0∪B0∪B1∪···∪Bl步骤十四：确定在时域内多失效模式下的危险概率对多失效模式的系统在使用期[0,T]内危险概率表示为：pf,c(0,T)=pf,c(AL=&cu...

【技术特征摘要】
1.针对多失效模式关联系统的安全性概率风险评价方法，其特征是，步骤一：识别影响模式关联系统安全性的关联失效模式，步骤二：定义模式关联系统安全性概率风险评价模型，假设模式关联系统存在m个失效模式，则对应m个时变极限状态函数，第i个时变极限状态函数定义为：其中和Di(.)为第i个极限状态函数的广义强度和广义应力，ω表示样本空间；步骤三：定义步骤二中随机向量X(t,ω)随机特性，包括分布类型、均值和标准差；步骤四：确定安全性概率风险评价模型中的失效域和安全域；根据步骤二定义的时变极限状态函数，当G(t,X(t,ω))＞0表示安全域，G(t,X(t,ω))＝0表示极限状态，G(t,X(t,ω))＜0表示失效域；步骤五：定义单一失效事件集合和危险概率；使用期[0,T]内，模式关联系统单个失效表示为：累积失效概率为在时刻ti的瞬时失效概率定义为pf,i(ti)＝P{G(ti,X(ti,ω))≤0}，其中，P(.)表示概率函数，下标f表示失效；步骤六：连续时间离散化；将连续的时间区间[0,T]离散化，设时间步长为h，则离散的时间区间数量为L＝T/h，第l个时间步为tl＝l×h，其中l＝0,1,…,L,T为系统使用周期；步骤七：定义第i个失效模式对应的失效域；设模式关联系统第i个时变极限状态函数为Gi(t,X(t,ω))，则系统在tl时刻的第i个失效模式对应的失效域定义为El,i＝{x∈X,Gi(tl,x(tl,ω))≤0}；步骤八：建立单个并联子系统对应的瞬时失效域；假设系统由并联和串联子系统构成，则由Nj个零部件组成的第j个并联子系统在tl时刻的瞬时失效域为：步骤九：建立单个串联子系统对应的瞬时失效域；同理步骤八，由Sj个零部件组成的第j个串联子系统在tl时刻的瞬时失效域为步骤十：建立若干个子系统构成整个失效模式关联系统瞬时失效域根据模式关联系统的构成方式，采取以下两种方法之一：1)建立若干个并联子系统以串联的形式构成整个系统瞬时失效域；由n个并联子系统以串联的形式构成的系统，在tl时刻的瞬时失效域为：

【专利技术属性】
技术研发人员：宫綦，李岩，苏多，刘欣，李晓勋，冯畅，
申请(专利权)人：中国航空综合技术研究所，
类型：发明
国别省市：