发明名称 针对多失效模式关联系统的安全性概率风险评价方法
摘要 本发明属于系统安全性技术研究领域,涉及一种利用集合理论针对多失效模式关联系统的安全性概率风险评价方法。本发明首先针对复杂系统多失效模式关联特点,利用集理论将系统的失效模式通过多个时变极限状态函数对应的失效域来表示,解决了多失效模式下系统安全性建模难题;再针对失效模式关联系统退化特性,考虑时间因素,得到多模式多失效模式下系统安全性模型随时间变化情况;在得到多失效模式对应的累积失效域基础上,利用蒙特卡洛方法进行安全性风险概率分析评价,解决了多失效模式相关且具有时变特性的系统安全性定量评价难题。
申请公布号 CN103559422B 申请公布日期 2017.04.19
申请号 CN201310598576.8 申请日期 2013.11.25
申请人 中国航空综合技术研究所 发明人 宫綦;李岩;苏多;刘欣;李晓勋;冯畅
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 中国航空专利中心 11008 代理人 李建英
主权项 针对多失效模式关联系统的安全性概率风险评价方法,其特征是,步骤一:识别影响模式关联系统安全性的关联失效模式,步骤二:定义模式关联系统安全性概率风险评价模型,假设模式关联系统存在m个失效模式,则对应m个时变极限状态函数,第i个时变极限状态函数定义为:<img file="FDA0001133694620000011.GIF" wi="1198" he="63" />其中<img file="FDA0001133694620000012.GIF" wi="51" he="61" />和D<sub>i</sub>(.)为第i个极限状态函数的广义强度和广义应力,ω表示样本空间;步骤三:定义步骤二中随机向量X(t,ω)随机特性,包括分布类型、均值和标准差;步骤四:确定安全性概率风险评价模型中的失效域和安全域;根据步骤二定义的时变极限状态函数,当G(t,X(t,ω))>0表示安全域,G(t,X(t,ω))=0表示极限状态,G(t,X(t,ω))<0表示失效域;步骤五:定义单一失效事件集合和危险概率;使用期[0,T]内,模式关联系统单个失效表示为:<img file="FDA0001133694620000013.GIF" wi="699" he="59" />累积失效概率为<img file="FDA0001133694620000014.GIF" wi="1067" he="63" />在时刻t<sub>i</sub>的瞬时失效概率定义为p<sub>f,i</sub>(t<sub>i</sub>)=P{G(t<sub>i</sub>,X(t<sub>i</sub>,ω))≤0},其中,P(.)表示概率函数,下标f表示失效;步骤六:连续时间离散化;将连续的时间区间[0,T]离散化,设时间步长为h,则离散的时间区间数量为L=T/h,第l个时间步为t<sub>l</sub>=l×h,其中l=0,1,…,L,T为系统使用周期;步骤七:定义第i个失效模式对应的失效域;设模式关联系统第i个时变极限状态函数为G<sub>i</sub>(t,X(t,ω)),则系统在t<sub>l</sub>时刻的第i个失效模式对应的失效域定义为E<sub>l,i</sub>={x∈X,G<sub>i</sub>(t<sub>l</sub>,x(t<sub>l</sub>,ω))≤0};步骤八:建立单个并联子系统对应的瞬时失效域;假设系统由并联和串联子系统构成,则由N<sub>j</sub>个零部件组成的第j个并联子系统在t<sub>l</sub>时刻的瞬时失效域为:<img file="FDA0001133694620000015.GIF" wi="236" he="138" />步骤九:建立单个串联子系统对应的瞬时失效域;同理步骤八,由S<sub>j</sub>个零部件组成的第j个串联子系统在t<sub>l</sub>时刻的瞬时失效域为<img file="FDA0001133694620000016.GIF" wi="267" he="143" />步骤十:建立若干个子系统构成整个失效模式关联系统瞬时失效域根据模式关联系统的构成方式,采取以下两种方法之一:1)建立若干个并联子系统以串联的形式构成整个系统瞬时失效域;由n个并联子系统以串联的形式构成的系统,在t<sub>l</sub>时刻的瞬时失效域为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>E</mi><mi>l</mi></msub><mo>=</mo><msub><mi>L</mi><mrow><mi>l</mi><mo>,</mo><mn>1</mn></mrow></msub><mo>&cup;</mo><msub><mi>L</mi><mrow><mi>l</mi><mo>,</mo><mn>2</mn></mrow></msub><mo>&cup;</mo><mo>...</mo><mo>&cup;</mo><msub><mi>L</mi><mrow><mi>l</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><msubsup><mo>&cup;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>L</mi><mrow><mi>l</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001133694620000021.GIF" wi="1470" he="84" /></maths>2)建立若干个串联子系统以并联的形式构成整个系统瞬时失效域;由s个串联子系统以并联的形式构成的系统,在t<sub>l</sub>时刻的瞬时失效域为:<maths num="0002"><math><![CDATA[<mrow><msub><mi>E</mi><mi>l</mi></msub><mo>=</mo><msub><mi>M</mi><mrow><mi>l</mi><mo>,</mo><mn>1</mn></mrow></msub><mo>&cap;</mo><msub><mi>M</mi><mrow><mi>l</mi><mo>,</mo><mn>2</mn></mrow></msub><mo>&cap;</mo><mo>...</mo><mo>&cap;</mo><msub><mi>M</mi><mrow><mi>l</mi><mo>,</mo><mi>s</mi></mrow></msub><mo>=</mo><msubsup><mo>&cup;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>s</mi></msubsup><msub><mi>M</mi><mrow><mi>l</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001133694620000022.GIF" wi="1507" he="85" /></maths>步骤十一:确定系统的累积失效域,通过步骤十确定了系统在t<sub>l</sub>时刻的瞬时失效域,在此基础上,得到系统在[0,t<sub>l</sub>]内的累积失效域:<img file="FDA0001133694620000023.GIF" wi="628" he="92" />步骤十二:确定系统离散时间内增加的失效域;系统在离散时间区间[t<sub>l</sub>,t<sub>l</sub>+h]内增加的失效域表示为:<maths num="0003"><math><![CDATA[<mrow><msub><mi>B</mi><mi>l</mi></msub><mo>=</mo><msub><mi>A</mi><mrow><mi>l</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&cap;</mo><msub><mover><mi>A</mi><mo>&OverBar;</mo></mover><mi>l</mi></msub><mo>=</mo><mrow><mo>(</mo><msubsup><mo>&cup;</mo><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>l</mi><mo>+</mo><mn>1</mn></mrow></msubsup><msub><mi>E</mi><mi>q</mi></msub><mo>)</mo></mrow><mo>&cap;</mo><mrow><mo>(</mo><msubsup><mo>&cap;</mo><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>l</mi></msubsup><msub><mover><mi>E</mi><mo>&OverBar;</mo></mover><mi>q</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001133694620000024.GIF" wi="1502" he="93" /></maths>通过布尔法则进一步简化为<img file="FDA0001133694620000025.GIF" wi="571" he="114" />步骤十三:通过离散时间内的增长失效域表示系统累积失效域;由步骤十一和十二,系统的累积失效域通过各个离散时间区间的增长失效域来表示,即<img file="FDA0001133694620000026.GIF" wi="1625" he="87" />步骤十四:确定在时域内多失效模式下的危险概率对多失效模式的系统在使用期[0,T]内危险概率表示为:<maths num="0004"><math><![CDATA[<mrow><msub><mi>p</mi><mrow><mi>f</mi><mo>,</mo><mi>c</mi></mrow></msub><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mi>T</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>p</mi><mrow><mi>f</mi><mo>,</mo><mi>c</mi></mrow></msub><mrow><mo>(</mo><msub><mi>A</mi><mi>L</mi></msub><mo>=</mo><msubsup><mo>&cup;</mo><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>L</mi></msubsup><msub><mi>E</mi><mi>q</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>P</mi><mrow><mo>(</mo><msub><mi>E</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>P</mi><mrow><mo>(</mo><msub><mi>B</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>P</mi><mrow><mo>(</mo><msub><mi>B</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mo>...</mo><mo>+</mo><mi>P</mi><mrow><mo>(</mo><msub><mi>B</mi><mi>l</mi></msub><mo>)</mo></mrow><mo>+</mo><mo>...</mo><mi>P</mi><mrow><mo>(</mo><msub><mi>B</mi><mrow><mi>L</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001133694620000027.GIF" wi="1716" he="89" /></maths>步骤十五:定义子系统瞬时失效和安全状态函数;子系统在t<sub>l</sub>时刻的失效和安全状态定义为<img file="FDA0001133694620000028.GIF" wi="660" he="150" />步骤十六:在瞬时时刻进行蒙特卡洛抽样;利用蒙特卡洛抽取N个随机样本,若用X<sub>k</sub>表示第k组随机样本,则在t<sub>l</sub>时刻系统失效的次数统计为:<img file="FDA0001133694620000029.GIF" wi="454" he="126" />其中<img file="FDA00011336946200000210.GIF" wi="1045" he="147" />步骤十七:在时间区间内进行蒙特卡洛抽样;在[t<sub>l</sub>,t<sub>l+1</sub>]区间内系统失效次数统计为:<img file="FDA0001133694620000031.GIF" wi="622" he="133" />其中<img file="FDA0001133694620000032.GIF" wi="1670" he="149" />步骤十八:利用蒙特卡洛计算离散事件区间内的累积失效概率,在离散时间区间[t<sub>l</sub>,t<sub>l+1</sub>]增加失效域对应的失效概率表达为:<maths num="0005"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>B</mi><mi>l</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>N</mi><mi>f</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>l</mi></msub><mo>,</mo><msub><mi>t</mi><mrow><mi>l</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>N</mi><mi>f</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>l</mi></msub><mo>)</mo></mrow></mrow><mi>N</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001133694620000033.GIF" wi="1620" he="125" /></maths>步骤十九:将步骤十八的失效概率式子带入步骤十四,即可得到系统在使用期内累积失效概率,从而完成对多失效模式关联复杂系统失效概率仿真分析。
地址 100028 北京市朝阳区京顺路7号