发明名称 基于风险与成本分析的飞机结构检查任务规划方法
摘要 本发明公开基于风险与成本分析的飞机结构检查任务规划方法,检查方式包括传统的结构无损检测和结构健康监测。首先基于概率损伤容限分析原理,建立飞机结构元件单次飞行失效概率的计算模型,并在此基础上建立机队某一时间段内的累积风险计算模型;然后在现行的飞机结构计划维修模式下,建立飞机结构件寿命周期内维修成本计算模型,包括结构检查成本、修理成本以及虚警成本的计算模型;最后,以单次飞行失效风险和机队累积风险作为结构维修策略的风险评价指标,在满足安全风险要求的前提下,依据结构全寿命周期维修成本的分析,可确定最佳的结构检查策略。可以在权衡结构失效风险和寿命周期维护成本的基础上确定最优的结构检查间隔和修理阈值。
申请公布号 CN106295897A 申请公布日期 2017.01.04
申请号 CN201610671617.5 申请日期 2016.08.15
申请人 南京航空航天大学 发明人 孙见忠;陈丹;李超役;闫洪胜;左洪福
分类号 G06Q10/04(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 张耀文
主权项 基于风险与成本分析的飞机结构检查任务规划方法,其特征在于,包括如下步骤:第一步:基于概率损伤容限分析原理,建立飞机结构件单次飞行失效概率计算模型,进而建立机队累积风险计算模型;第二步:基于现行的飞机结构计划维修模式,建立飞机结构寿命周期内维修成本计算模型,包括结构检查成本、修理成本以及虚警事件引发的结构检查成本的计算模型;所述结构检查成本,具体如下:单次检查成本用工时费乘以人工时,因此对于一个给定的维修策略在寿命周期内的m次结构检查的总成本如下:<maths num="0001"><math><![CDATA[<mrow><msub><mi>C</mi><mi>I</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msubsup><mi>c</mi><mi>I</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup></mrow>]]></math><img file="FDA0001078322220000011.GIF" wi="197" he="103" /></maths>其中,式中<img file="FDA0001078322220000012.GIF" wi="59" he="54" />第j次检查成本;所述修理成本,具体如下:引进裂纹修复概率PCR,定义如下:<maths num="0002"><math><![CDATA[<mrow><mi>P</mi><mi>C</mi><mi>R</mi><mo>=</mo><mi>E</mi><mrow><mo>(</mo><mi>P</mi><mi>O</mi><mi>R</mi><mo>(</mo><mi>a</mi><mo>)</mo><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>&infin;</mi></msubsup><mi>P</mi><mi>O</mi><mi>R</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mi>f</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mi>d</mi><mi>a</mi></mrow>]]></math><img file="FDA0001078322220000013.GIF" wi="677" he="76" /></maths>式中:f(a)表示检查时裂纹尺寸的概率密度函数,POR(a)表示当裂纹尺寸为a的时候裂纹修理的概率;POR(a)定义如下:<maths num="0003"><math><![CDATA[<mrow><mi>P</mi><mi>O</mi><mi>R</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><mi>&phi;</mi><mrow><mo>(</mo><mfrac><mrow><mi>l</mi><mi>n</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>-</mo><mi>&mu;</mi></mrow><mi>&sigma;</mi></mfrac><mo>)</mo></mrow><mo>=</mo><mi>&phi;</mi><mrow><mo>(</mo><mfrac><mrow><mi>l</mi><mi>n</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>-</mo><mi>l</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>a</mi><mrow><mi>r</mi><mi>e</mi><mi>s</mi></mrow></msub><mo>)</mo></mrow></mrow><mi>&sigma;</mi></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001078322220000014.GIF" wi="790" he="121" /></maths>式中:σ是根据无损检测系统可靠性评估手册1823A推荐的方法来分析NDE或SHM系统实验数据获得的,μ取修理决策阈值的对数值,即log(a<sub>res</sub>);根据修理决策阈值的定义,当裂纹尺寸等于或超过阈值时进行修理,否则就不进行修理;所述虚警事件引发的结构检查成本,具体如下:定义结构修理虚警事件为NDE或SHM系统检测到损伤尺寸已经达到了修理决策阈值,必须停机修理,但是实际上裂纹尺寸并没有超过阈值;对于一个给定的维修策略在寿命周期内的m次检查计划中,由于虚警引起的额外维修成本表示如下:<maths num="0004"><math><![CDATA[<mrow><msub><mi>C</mi><mrow><mi>F</mi><mi>C</mi><mi>R</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msubsup><mi>c</mi><mrow><mi>F</mi><mi>C</mi><mi>R</mi></mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>&CenterDot;</mo><msubsup><mi>p</mi><mrow><mi>F</mi><mi>C</mi><mi>R</mi></mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup></mrow>]]></math><img file="FDA0001078322220000015.GIF" wi="342" he="103" /></maths>式中:<img file="FDA0001078322220000016.GIF" wi="67" he="55" />表示第j次检查时由于虚警事件引起的额外的检查和维修成本,<img file="FDA0001078322220000017.GIF" wi="74" he="55" />表示第j次检查时修理虚警事件发生的概率;第三步:在同时满足飞机结构单次飞行失效风险及机队累积失效风险的前提下,选择无损检测设备NDE或在线结构健康监测系统SHM,优化结构检查间隔和修理决策阈值,使飞机结构寿命周期内维修成本最低;<maths num="0005"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>min</mi><mo>:</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>C</mi><mrow><mo>(</mo><mi>I</mi><mo>,</mo><msub><mi>a</mi><mrow><mi>r</mi><mi>e</mi><mi>s</mi><mo>,</mo><mi>d</mi><mi>e</mi><mi>p</mi><mi>o</mi><mi>t</mi></mrow></msub><mo>,</mo><msub><mi>a</mi><mrow><mi>r</mi><mi>e</mi><mi>s</mi><mo>,</mo><mi>f</mi><mi>i</mi><mi>e</mi><mi>l</mi><mi>d</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>min</mi><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><msubsup><mi>c</mi><mi>I</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>+</mo><msubsup><mi>c</mi><mi>R</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>&CenterDot;</mo><msubsup><mi>p</mi><mi>R</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>+</mo><msubsup><mi>c</mi><mrow><mi>F</mi><mi>C</mi><mi>R</mi></mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>&CenterDot;</mo><msubsup><mi>p</mi><mrow><mi>F</mi><mi>C</mi><mi>R</mi></mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>max</mi><mrow><mo>(</mo><msub><mi>SFPOF</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>&lt;</mo><msub><mi>SFPOF</mi><mrow><mi>t</mi><mi>h</mi><mi>r</mi><mi>e</mi><mi>s</mi><mi>h</mi><mi>o</mi><mi>l</mi><mi>d</mi></mrow></msub><mo>,</mo><mi>t</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>m</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>R</mi><mi>m</mi></msub><mo>&lt;</mo><msub><mi>R</mi><mrow><mi>t</mi><mi>h</mi><mi>r</mi><mi>e</mi><mi>s</mi><mi>h</mi><mi>o</mi><mi>l</mi><mi>d</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0001078322220000021.GIF" wi="1022" he="439" /></maths>其中:C(I,a<sub>res,depot</sub>,a<sub>res,field</sub>)表示结构寿命周期内维修成本;I表示检查间隔;a<sub>res,depot</sub>表示基地结构修理阈值;a<sub>res,field</sub>表示航线结构修理阈值;s.t.表示约束条件;max(SFPOF<sub>t</sub>)表示结构寿命周期内SFPOF最大值;SFPOF<sub>threshold</sub>表示能够接受的SFPOF阈值;R<sub>m</sub>表示未来m个飞行循环内机队累积风险;R<sub>threshold</sub>表示能够接受的机队累积风险阈值。
地址 210016 江苏省南京市秦淮区御道街29号