发明名称 作战飞机系统易损性指标分配方法
摘要 一种作战飞机系统易损性指标分配方法,通过建立某型战斗机易损性模型,在确定威胁打击下,确定易损性初始参数,然后确定飞机遭受打击的次数以及每个部件在遭受打击后的杀伤概率和生存概率,接着确定飞机独立存在的状态和各状态被打击后的概率,将上述所得数据代入综合重要度和Griffith重要度计算公式,得到各个部件的综合重要度和Griffith重要度,最后利用易损性指标分配法对全机易损性进行分配,可得到各个部件的易损性。得到每个部件分配的易损性后,针对易损性大的部件进行易损性减缩,降低其易损性,从而提高飞机的生存力。
申请公布号 CN104820778A 申请公布日期 2015.08.05
申请号 CN201510205560.5 申请日期 2015.04.27
申请人 西北工业大学 发明人 裴扬;赵倩;田晨;谢芳敏;宋笔锋
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 西北工业大学专利中心 61204 代理人 慕安荣
主权项 一种作战飞机系统易损性指标分配方法,其特征在于,具体过程是:步骤1:建立飞机易损性模型:确定飞机的致命性部件;所述致命性部件是飞行员B、尾部燃油箱C、左发动机D、右发动机F、左液压油箱G、右液压油箱H、左发电机I和右发电机J;其中,飞行员B、尾部燃油箱C、左发动机D和右发动机F均为非余度部件;左液压油箱G和右液压油箱H互为余度部件,左发电机I和右发电机J互为余度部件;所述致命性部件的最小割集为:{B、C、D、F、{G、H}、{I、J}};根据最小割集绘制飞机杀伤树;步骤2:确定战斗机被攻击的方向:确定飞机被攻击方向的方位角A为90°,俯仰角E为0°;步骤3:定义易损性初始参数;所述的易损性初始参数包括部件i的暴露面积<img file="FDA0000706461950000011.GIF" wi="101" he="80" />威胁打击部件时部件i的杀伤概率<img file="FDA0000706461950000012.GIF" wi="82" he="81" />和威胁打击飞机时部件i的杀伤概率<img file="FDA0000706461950000013.GIF" wi="131" he="82" />通过得到的部件i的暴露面积<img file="FDA0000706461950000014.GIF" wi="68" he="82" />和威胁打击部件时部件i的杀伤概率<img file="FDA0000706461950000015.GIF" wi="116" he="77" />确定威胁打击飞机时部件i的杀伤概率<img file="FDA0000706461950000016.GIF" wi="130" he="78" />将暴露面积<img file="FDA00007064619500000120.GIF" wi="100" he="70" />威胁打击部件时部件i的杀伤概率<img file="FDA0000706461950000017.GIF" wi="87" he="76" />带入A<sub>Vi</sub>=A<sub>Pi</sub>P<sub>k/hi</sub>,得到每个部件的易损面积<img file="FDA0000706461950000018.GIF" wi="97" he="80" />将得到的易损面积<img file="FDA0000706461950000019.GIF" wi="63" he="78" />和飞机的暴露面积A<sub>P</sub>带入<img file="FDA00007064619500000110.GIF" wi="322" he="84" />得到威胁打击飞机时部件i的杀伤概率<img file="FDA00007064619500000111.GIF" wi="125" he="77" />所述飞机暴露面积A<sub>P</sub>中A表示面积,下标P表示暴露在外;飞机易损面积A<sub>V</sub>中A表示面积,下标V表示易损性;部件i的暴露面积<img file="FDA00007064619500000112.GIF" wi="68" he="78" />中A表示面积,下标P表示暴露在外,i表示部件编号;威胁打击部件时部件i的杀伤概率<img file="FDA00007064619500000113.GIF" wi="84" he="80" />中P表示概率,下标k表示杀伤,h表示打击在部件上,i表示部件编号;威胁打击飞机时部件i的杀伤概率<img file="FDA00007064619500000114.GIF" wi="98" he="85" />中P表示概率,下标k表示杀伤,H表示打击在飞机上,i表示部件编号;基于暴露面积上的期望打击次数<img file="FDA00007064619500000115.GIF" wi="80" he="84" />中E是所述的期望打击次数,下标A表示面积,P表示暴露在外;步骤4:确定每个部件遭受<img file="FDA00007064619500000116.GIF" wi="80" he="79" />次打击后的杀伤概率和生存概率:对于飞机遭受<img file="FDA00007064619500000117.GIF" wi="79" he="77" />次随机打击,利用二项式方法确定飞机遭遇<img file="FDA00007064619500000118.GIF" wi="76" he="77" />次打击后部件i的生存概率为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>P</mi><mrow><mi>s</mi><mo>/</mo><msub><mi>H</mi><mi>i</mi></msub></mrow><msub><mi>E</mi><msub><mi>A</mi><mi>P</mi></msub></msub></msubsup><mo>=</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>/</mo><msub><mi>H</mi><mi>i</mi></msub></mrow></msub><mo>)</mo></mrow><msub><mi>E</mi><msub><mi>A</mi><mi>P</mi></msub></msub></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00007064619500000119.GIF" wi="1460" he="87" /></maths>公式(1)中,P表示概率,s表示生存,<img file="FDA0000706461950000021.GIF" wi="80" he="83" />表示基于暴露面积上的期望打击次数,威胁打击飞机时部件i的杀伤概率<img file="FDA0000706461950000022.GIF" wi="110" he="77" />中P表示概率,下标中:k表示杀伤,H表示打击在飞机上,i表示部件编号;飞机遭受<img file="FDA0000706461950000023.GIF" wi="79" he="77" />次打击后部件i的杀伤概率为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>P</mi><mrow><mi>k</mi><mo>/</mo><msub><mi>H</mi><mi>i</mi></msub></mrow><msub><mi>E</mi><msub><mi>A</mi><mi>P</mi></msub></msub></msubsup><mo>=</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>/</mo><msub><mi>H</mi><mi>i</mi></msub></mrow></msub><mo>)</mo></mrow><msub><mi>E</mi><msub><mi>A</mi><mi>P</mi></msub></msub></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000706461950000024.GIF" wi="1442" he="95" /></maths>公式(2)中,P表示概率,k表示杀伤,<img file="FDA0000706461950000025.GIF" wi="85" he="76" />表示基于暴露面积上的期望打击次数,威胁打击飞机时部件i的杀伤概率<img file="FDA0000706461950000026.GIF" wi="106" he="80" />中P表示概率,下标中:k表示杀伤,H表示打击在飞机上,i表示部件编号;将表1中相应数据带入公式(1)和公式(2)中,得到每个部件遭受<img file="FDA0000706461950000027.GIF" wi="84" he="78" />次打击后的杀伤概率和生存概率;步骤5:确定飞机的独立存在状态j及飞机受到<img file="FDA0000706461950000028.GIF" wi="72" he="77" />次打击后所述各独立存在状态发生的概率:利用排列组合法确定飞机所有的杀伤状态;分别把导致飞机杀伤的状态合并,得到所述各飞机独立存在的状态在<img file="FDA0000706461950000029.GIF" wi="70" he="80" />次打击后的概率值;步骤6:给定各独立状态的系统性能a<sub>j</sub>,其中j=0,1,2,…,q,分别代表了飞机独立存在的各种状态;步骤7:确定飞机各部件的重要度VI:所述的飞机各部件的重要度包括综合重要度VI<sup>IIM</sup>和Griffith重要度VI<sup>G</sup>;所述确定综合重要度和Griffith重要度是通过得到的各部件的杀伤概率,采用公式(3)和公式(4)确定;飞机在作战前,各部件i的可用度P<sub>i1</sub>=1;部件的状态空间为{0,1},其中0表示杀伤状态,1表示生存状态;飞机状态空间为{0,1,2,...,M},部件i的综合重要度通过公式(3)确定:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>VI</mi><mn>1</mn><mi>IIM</mi></msubsup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>P</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>&CenterDot;</mo><msubsup><mi>&lambda;</mi><mn>1,0</mn><mi>i</mi></msubsup><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>a</mi><mi>j</mi></msub><mo>[</mo><msub><mi>Pr</mi><msub><mi>E</mi><msub><mi>A</mi><mi>P</mi></msub></msub></msub><mrow><mo>(</mo><mi>&Phi;</mi><mrow><mo>(</mo><msub><mn>1</mn><mi>i</mi></msub><mo>,</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Pr</mi><msub><mi>E</mi><msub><mi>A</mi><mi>P</mi></msub></msub></msub><mrow><mo>(</mo><mi>&Phi;</mi><mrow><mo>(</mo><msub><mn>0</mn><mi>i</mi></msub><mo>,</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00007064619500000210.GIF" wi="1712" he="156" /></maths>部件i的Griffith重要度通过公式(4)确定:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>VI</mi><mn>1</mn><mi>G</mi></msubsup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>a</mi><mi>j</mi></msub><mo>[</mo><msub><mi>Pr</mi><msub><mi>E</mi><msub><mi>A</mi><mi>P</mi></msub></msub></msub><mrow><mo>(</mo><mi>&Phi;</mi><mrow><mo>(</mo><msub><mn>1</mn><mi>i</mi></msub><mo>,</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Pr</mi><msub><mi>E</mi><msub><mi>A</mi><mi>P</mi></msub></msub></msub><mrow><mo>(</mo><mi>&Phi;</mi><mrow><mo>(</mo><msub><mn>0</mn><mi>i</mi></msub><mo>,</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00007064619500000211.GIF" wi="1668" he="142" /></maths>公式(3)中,<img file="FDA0000706461950000031.GIF" wi="77" he="81" />为部件i遭受单次打击从完好状态劣化到杀伤状态的转化率,该部分的转化是由于遭受威胁的打击引起的,因此<img file="FDA0000706461950000032.GIF" wi="242" he="88" />所述的完好状态和杀伤状态分别用1和0表示;公式(3)和式(4)中:<img file="FDA0000706461950000033.GIF" wi="406" he="91" />表示当部件i处于完好状态时,飞机为第j种状态的概率;<img file="FDA0000706461950000034.GIF" wi="418" he="99" />表示部件i杀伤状态时,飞机为状态j的概率;其中Pr表示概率,<img file="FDA0000706461950000035.GIF" wi="87" he="77" />表示基于暴露面积上的期望打击次数,1表示完好,0表示杀伤,X表示部件向量,Φ(1<sub>i</sub>,X)表示飞机结构函数,j表示飞机的第j个独立状态;a<sub>j</sub>表示系统系能;M表示飞机独立状态总数;VI<sup>IIM</sup>表示综合重要度;VI<sup>G</sup>表示Griffith重要度;步骤8:易损性指标分配:采用比例分配法对综合重要度和Griffith重要度进行易损性指标的分配;飞机系统各部件的综合重要度为VI<sup>IIM</sup>(i),Griffith重要度为VI<sup>G</sup>(i);系统的易损性为A<sub>V</sub>,依据系统性能a<sub>j</sub>对飞机系统各部件的重要度VI(i)分配时应赋予部件i的加权因子为ω<sub>i</sub>,通过公式(5)确定该加权因子:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><mi>VI</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>VI</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></mfrac><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000706461950000036.GIF" wi="1386" he="218" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&omega;</mi><mi>i</mi></msub><mo>=</mo><mn>1</mn><mo>;</mo></mrow>]]></math><img file="FDA0000706461950000037.GIF" wi="217" he="132" /></maths>公式(5)中,ω<sub>i</sub>是部件i的加权因子,下标i的取值范围都为1,…,n;通过公式(6)确定部件i分配的易损性指标为:A<sub>Vi</sub>=ω<sub>i</sub>·A<sub>V</sub>                  (6)公式(6)中,A<sub>V</sub>为飞机的易损性,A<sub>Vi</sub>为各个部件所分配得到的易损性;至此,完成了作战飞机系统易损性指标的分配。
地址 710072 陕西省西安市友谊西路127号