• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 11
    • 下载费用:30 金币  

    重庆时时彩中奖条件: 一种大气中子对航空器件单粒子效应诱发概率的评估方法.pdf

    摘要
    申请专利号:

    重庆时时彩单双窍门 www.4mum.com.cn CN201610978929.0

    申请日:

    2016.11.08

    公开号:

    CN106570645A

    公开日:

    2017.04.19

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06Q 10/06申请日:20161108|||公开
    IPC分类号: G06Q10/06(2012.01)I; G06Q50/26(2012.01)I; G06Q50/30(2012.01)I 主分类号: G06Q10/06
    申请人: 北京天工科仪空间技术有限公司
    发明人: 杨垂柏; 冯宇波; 王世金; 王月
    地址: 100193 北京市海淀区天秀路10号中国农大国际创业园3号楼4020室
    优先权:
    专利代理机构: 代理人:
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201610978929.0

    授权公告号:

    |||

    法律状态公告日:

    2017.05.17|||2017.04.19

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明提出了一种用于评估大气中子诱发航空电子器件发生单粒子事件效应的概率的方法,该方法利用航空飞行器的飞行航迹数据作为持续动态输入,包括时间、高度、纬度、经度,利用航空航迹计算出大气中子通量,利用粒子输运工具计算出不同能量中子在卫星结构材料中的射程,并结合需要评估的空间器件的相关参数,在获得目标位置处的大气中子通量后,结合不同能量下航空飞行器电子器件发生单粒子事件效应的概率之间的关系,最后获得在目标位置处在航空飞行器电子器件中发生单粒子事件效应的概率。

    权利要求书

    1.一种航空飞行器的器件单粒子效应发生概率的评估方法,其特征在于包括:
    步骤1:获得航空飞行器在飞行过程所处太阳活动周期内所处阶段参数;
    步骤2:获得航空器在飞行过程的航迹点的高度、纬度、经度数值序列;
    步骤3:设置大气中子的能量范围和能量点;
    步骤4:依据飞行时间、航迹点高度、纬度、经度数值序列及中子能谱设置,利用大气中
    子计算模型计算出不同航迹点处的大气中子密度分布、能谱分布;
    步骤4.1:依据航迹点纬度、经度数值序列,计算出地磁纬度、经度数值序列;
    Latm=90-arcos[sin(Latg)*sin(lat0)+cos(Latg)*cos(Lat0)
    *cos(long-lon0)]
    其中:Latm为航迹点地磁纬度,Latg为航迹点地理纬度,Long为航迹点地理经度,Lat0为
    地磁北极地理纬度,Lon0为地磁北极地理经度;
    步骤4.2:依据航迹点高度、地磁纬度、地磁经度数值序列,计算出航迹点处的磁场刚度
    值序列,利用如下公式:
    <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>=</mo> <mn>14.9</mn> <mo>*</mo> <msup> <mi>cos</mi> <mn>4</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>Lat</mi> <mi>m</mi> </msub> <mo>)</mo> </mrow> <mo>/</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>x</mi> <mo>+</mo> <mi>R</mi> </mrow> <mi>R</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow>
    其中:Rc为磁刚度,x为航迹高度值,R为地球半径值;
    步骤4.3:依据步骤1、2的飞行时间、航迹纬度值及磁刚度序列,计算出所有航迹点大气
    中子的通量值,利用如下公式:
    λ=165+2*Rc
    Xm=50+ln{2000+exp[-x*(Cs-100)]}
    <mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&Phi;</mi> <mrow> <mn>1</mn> <mo>-</mo> <mn>10</mn> </mrow> </msub> <mrow> <mo>(</mo> <mn>250</mn> <mo>,</mo> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mn>0.17</mn> <mo>+</mo> <mo>&lsqb;</mo> <mn>0.787</mn> <mo>+</mo> <mn>0.035</mn> <mo>*</mo> <mrow> <mo>(</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>-</mo> <mn>100</mn> <mo>)</mo> </mrow> <mo>&rsqb;</mo> <mo>*</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msup> <msub> <mi>R</mi> <mi>c</mi> </msub> <mn>2</mn> </msup> <mo>/</mo> <mn>25</mn> </mrow> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mo>{</mo> <mo>-</mo> <mn>0.107</mn> <mo>-</mo> <mn>0.0265</mn> <mo>*</mo> <mrow> <mo>(</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>-</mo> <mn>100</mn> <mo>)</mo> </mrow> <mo>+</mo> <mn>0.612</mn> <mo>*</mo> <mi>exp</mi> <mo>&lsqb;</mo> <mfrac> <mrow> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>-</mo> <mn>100</mn> </mrow> <mn>3.73</mn> </mfrac> <mo>&rsqb;</mo> <mo>}</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>*</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msup> <msub> <mi>R</mi> <mi>c</mi> </msub> <mn>2</mn> </msup> <mo>/</mo> <mn>139.2</mn> </mrow> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced>
    <mrow> <mi>f</mi> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mfrac> <mn>250</mn> <mi>&lambda;</mi> </mfrac> <mo>)</mo> </mrow> <mo>*</mo> <msub> <mi>&Phi;</mi> <mrow> <mn>1</mn> <mo>-</mo> <mn>10</mn> </mrow> </msub> <mrow> <mo>(</mo> <mn>250</mn> <mo>,</mo> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> </mrow>
    <mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&Phi;</mi> <mrow> <mn>1</mn> <mo>-</mo> <mn>10</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>m</mi> </msub> <mo>,</mo> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>R</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mn>0.23</mn> <mo>+</mo> <mo>&lsqb;</mo> <mn>1.1</mn> <mo>+</mo> <mn>0.0167</mn> <mo>*</mo> <mrow> <mo>(</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>-</mo> <mn>100</mn> <mo>)</mo> </mrow> <mo>&rsqb;</mo> <mo>*</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <msup> <msub> <mi>R</mi> <mi>c</mi> </msub> <mn>2</mn> </msup> </mrow> <mn>81</mn> </mfrac> </mrow> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mo>{</mo> <mn>0.991</mn> <mo>+</mo> <mn>0.051</mn> <mo>*</mo> <mrow> <mo>(</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>-</mo> <mn>100</mn> <mo>)</mo> </mrow> <mo>+</mo> <mn>0.4</mn> <mo>*</mo> <mi>exp</mi> <mo>&lsqb;</mo> <mrow> <mo>(</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>-</mo> <mn>100</mn> <mo>)</mo> </mrow> <mo>/</mo> <mn>3.73</mn> <mo>&rsqb;</mo> <mo>}</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>*</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msup> <msub> <mi>R</mi> <mi>c</mi> </msub> <mn>2</mn> </msup> <mo>/</mo> <mn>12.96</mn> </mrow> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced>
    <mrow> <mi>&Lambda;</mi> <mo>=</mo> <mo>&lsqb;</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>&Phi;</mi> <mrow> <mn>1</mn> <mo>-</mo> <mn>10</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>m</mi> </msub> <mo>,</mo> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>*</mo> <msup> <mi>e</mi> <mfrac> <msub> <mi>X</mi> <mi>m</mi> </msub> <mi>&lambda;</mi> </mfrac> </msup> <mo>/</mo> <mi>f</mi> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>&rsqb;</mo> </mrow>
    <mrow> <mi>F</mi> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Lambda;</mi> <mi>&lambda;</mi> </mfrac> <mo>)</mo> </mrow> <mo>*</mo> <mi>f</mi> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>*</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>X</mi> <mi>m</mi> </msub> <mi>&Lambda;</mi> </mfrac> <mo>-</mo> <msub> <mi>X</mi> <mi>m</mi> </msub> <mo>/</mo> <mi>&lambda;</mi> <mo>)</mo> </mrow> </mrow>
    <mrow> <msub> <mi>&Phi;</mi> <mrow> <mn>1</mn> <mo>-</mo> <mn>10</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>f</mi> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>*</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mi>X</mi> <mi>&lambda;</mi> </mfrac> </mrow> </msup> <mo>-</mo> <mi>F</mi> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>*</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>X</mi> <mo>/</mo> <mi>&Lambda;</mi> </mrow> </msup> </mrow>
    其中:x表示航迹高度值(单位:g/cm2),Cs为太阳活动参数,表征太阳活动时节,λ、Xm及
    Λ是计算过程的参数;
    步骤4.4:依据步骤3的中子能量范围和能量点和步骤4.3的通量计算结果,计算出所有
    航迹点随航迹点不同能量的大气中子通量值,利用如下公式计算航迹点的大气中子通量:
    <mrow> <mfrac> <mrow> <mi>d</mi> <mi>N</mi> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mi>n</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>dE</mi> <mi>n</mi> </msub> </mrow> </mfrac> <mo>=</mo> <mn>0.3459</mn> <mo>*</mo> <msup> <mn>10.0</mn> <mrow> <mo>-</mo> <mn>0.9219</mn> </mrow> </msup> <mo>*</mo> <mi>exp</mi> <mo>&lsqb;</mo> <mo>-</mo> <mn>0.01522</mn> <mo>*</mo> <msup> <mrow> <mo>(</mo> <mi>l</mi> <mi>n</mi> <mo>(</mo> <msub> <mi>E</mi> <mi>n</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&rsqb;</mo> </mrow>
    其中:N(En)表示大气中子通量,En表示大气中子能量点;
    步骤5:获取待评估的空间器件的参数,包括σsat、LETth、W、S、t,空间器件参数由地面试
    验或者数值仿真获得,
    其中:σsat为重离子造成的器件单粒子事件效应发生概率的饱和截面,
    LETth、W、S为概率分布函数Weibull的系数值,
    t为电子器件的灵敏区厚度;
    步骤6:依据步骤3给出的中子能谱,利用蒙特卡洛计算工具计算给出不同能量中子的
    反冲核概率BGR(E,Eri)值,
    其中:E表示中子能量,
    Eri表示粒子在电子器件灵敏区内的沉积能量,
    BGR(E,Eri)表示能量为E的中子产生能量大于Eri的反冲核的概率;
    步骤7:依据步骤4获得参数,利用如下公式计算单粒子事件发生截面σsee(i),
    <mrow> <msub> <mi>&sigma;</mi> <mrow> <mi>s</mi> <mi>e</mi> <mi>e</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&sigma;</mi> <mrow> <mi>s</mi> <mi>a</mi> <mi>t</mi> </mrow> </msub> <mo>&CenterDot;</mo> <mo>{</mo> <mn>1</mn> <mo>-</mo> <mi>exp</mi> <mo>&lsqb;</mo> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>LET</mi> <mrow> <mi>i</mi> <mi>o</mi> <mi>n</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>LET</mi> <mrow> <mi>t</mi> <mi>h</mi> </mrow> </msub> </mrow> <mi>W</mi> </mfrac> <mo>)</mo> </mrow> <mi>s</mi> </msup> <mo>&rsqb;</mo> <mo>}</mo> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msub> <mi>LET</mi> <mrow> <mi>i</mi> <mi>o</mi> <mi>n</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>&GreaterEqual;</mo> <msub> <mi>LET</mi> <mrow> <mi>t</mi> <mi>h</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>0</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msub> <mi>LET</mi> <mrow> <mi>i</mi> <mi>o</mi> <mi>n</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>&lt;</mo> <msub> <mi>LET</mi> <mrow> <mi>t</mi> <mi>h</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
    其中,LETion(i)表示重离子的LET值;
    步骤8:依据步骤6获得参数,利用如下公式计算单粒子事件发生概率,
    <mrow> <msub> <mi>f</mi> <mrow> <mi>s</mi> <mi>e</mi> <mi>e</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>C</mi> <mo>*</mo> <munder> <mo>&Sigma;</mo> <mi>i</mi> </munder> <msub> <mi>&Delta;V</mi> <mi>i</mi> </msub> <mo>&Integral;</mo> <mi>B</mi> <mi>G</mi> <mi>R</mi> <mrow> <mo>(</mo> <mi>E</mi> <mo>,</mo> <mi>E</mi> <mi>r</mi> <mi>i</mi> <mo>)</mo> </mrow> <mo>*</mo> <mrow> <mo>(</mo> <mi>d</mi> <mi>J</mi> <mo>/</mo> <mi>d</mi> <mi>E</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>E</mi> </mrow>
    其中,C表示吸收系数,对于绝对大多数电子器件C取0.5时于结果符合实验验证,ΔVi参
    数表示为tΔσi=t(σi-σi-1)。
    2.根据权利要求1所述的用于评估航空电子器件发生单粒子事件效应的概率的方法,
    其特征在于,还包括:
    步骤9:采用Python或Matlab软件将计算结果可视化。
    3.根据权利要求1所述的用于评估航空电子器件发生单粒子事件效应的概率的方法,
    其特征在于,步骤1的太阳活动参数Cs的范围为从65到102,在11年的太阳活动周期内,活动
    谷年的参数值为102,活动峰年的参数值为65。
    4.根据权利要求1或2所述的用于评估航空电子器件发生单粒子事件效应的概率的方
    法,其特征在于,所述航空飞行器的航迹取样点时间间隔不小于1个小时,否则将极大降低
    准确度。
    5.根据权利要求1或2所述的用于评估航空电子器件发生单粒子事件效应的概率的方
    法,其特征在于,所述输入的大气中能量范围,其能段低端小于等于10MeV、高端大于等于
    10MeV。
    6.根据权利要求1或2所述的用于评估航空电子器件发生单粒子事件效应的概率的方
    法,其特征在于,步骤4中的粒子输运仿真工具包括Geamt3、Geant4及MCNP系列。
    7.根据权利要求1或2所述的用于评估航空电子器件发生单粒子事件效应的概率的方
    法,其特征在于,所述航空飞行器结构物质包括金属铝、铝镁合金、钢制及复合材料,所述电
    子器件的灵敏物质包括硅、锗及砷化镓。
    8.根据权利要求1或2所述的用于评估航空电子器件发生单粒子事件效应的概率的方
    法,其特征在于,所述航空飞行器包括各类民用飞机、军用飞机、气球、导弹及高超音速飞行
    器。
    9.根据权利要求1或2所述的用于评估航空电子器件发生单粒子事件效应的概率的方
    法,其特征在于,所述航空飞行器的飞行高度不高于100km,不低于0km,所述航空飞行器的
    飞行纬度为-90度到90度以内。

    说明书

    一种大气中子对航空器件单粒子效应诱发概率的评估方法

    技术领域

    本发明涉及航空飞行过程的太空辐射环境效应危害防护管理领域,尤其涉及一种
    用于评估太空宇宙线与高层大气发生碰撞后的大气中子对航空电子器件造成的发生单粒
    子事件效应的概率的方法。

    背景技术

    太空宇宙线的高能电子、高能质子及重离子等要素都会对地球上空运行的航行器
    造成辐射效应危害,高能质子和重离子会造成被称为单粒子事件效应(Single Event
    Effect-SEE)的危害,在航空飞行器上也存在这类效应,但却是由于宇宙线于大气碰撞后的
    中子造成的。

    宇宙线粒子在穿越物质的过程中会由于受到原子核和核外电子形成的阻力而损
    失掉能量,而损失掉的能量将会转移给阻滞其的物质,单位内转移给阻滞物质的能量被称
    之为线性能量转移(Linear Energy Transfer-LET),当LET值超过一定值就可能导致空间
    器件发生单粒子事件效应。根据在器件内的危害表现单粒子事件效应分为单粒子翻转、单
    粒子锁闭、单粒子瞬态等。

    如果飞行器内的电子器件发生单粒子事件效应,轻则造成数字逻辑位的翻转,诸
    如由“0”变“1”或由“1”变“0”,造成干扰,重则造成空间器件烧毁等,导致单机故障、甚至整
    星故障,前者被称之为单粒子翻转效应(Single Event Upset-SEU),后者被称之为单粒子
    烧蚀效应(Single Event Burnout-SEB)。

    由于航空电子器件单粒子事件效应会对航空飞行器造成干扰危害,特别是随着高
    精度、长航程、长航时的无人机日益广泛的应用,因此,在航空器研制阶段、操控维护阶段及
    事后异常诊断阶段,开展对电子器件发生单粒子事件效应的概率的评估是一种降低由于单
    粒子事件效应造成的危害的重要手段。

    发明内容

    本发明的目的在于,为了避免航空飞行器在飞行执勤阶段的电子器件由于大气中
    子造成的单粒子事件效应的危害而评估该事件效应发生的概率,本发明提供一种用于评估
    航空飞行器的电子器件发生单粒子事件效应的概率的方法。

    为实现上述目的,本发明提出了一种用于评估大气中子诱发航空电子器件发生单
    粒子事件效应的概率的方法,该方法利用航空器的飞行航迹数据作为持续动态输入,包括
    时间、高度、纬度、经度,利用航空航迹计算出大气中子通量,并结合需要评估的电子器件的
    相关参数,在获得目标位置处的大气中子通量后,结合不同能量下航空器电子器件发生单
    粒子事件效应的概率之间的关系,最后获得在目标位置处在航空器电子器件中发生单粒子
    事件效应的概率。

    具体步骤如下:

    步骤1:获得航空器在飞行过程所处太阳活动周期内所处阶段参数;

    步骤2:获得航空器在飞行过程的航迹点的高度、纬度、经度数值序列;

    步骤3:设置大气中子的能量范围和能量点;

    步骤4:依据飞行时间、航迹点高度、纬度、经度数值序列及中子能谱设置,利用大
    气中子计算模型计算出不同航迹点处的大气中子密度分布、能谱分布;

    步骤4.1:依据航迹点纬度、经度数值序列,计算出地磁纬度、经度数值序列;

    Latm=90-arcos[sin(Latg)*sin(lat0)+cos(Latg)*cos(Lat0)*cos(long-lon0)]

    其中:Latm为航迹点地磁纬度,Latg为航迹点地理纬度,Long为航迹点地理经度,
    Lat0为地磁北极地理纬度,Lon0为地磁北极地理经度;

    步骤4.2:依据航迹点高度、地磁纬度、地磁经度数值序列,计算出航迹点处的磁场
    刚度值序列,利用如下公式:


    其中:Rc为磁刚度,x为航迹高度值,R为地球半径值;

    步骤4.3:依据步骤1、2的飞行时间、航迹纬度值及磁刚度序列,计算出所有航迹点
    大气中子的通量值,利用如下公式:

    λ=165+2*Rc

    Xm=50+ln{2000+exp[-x*(Cs-100)]}







    其中:x表示航迹高度值(单位:g/cm2),Cs为太阳活动参数,表征太阳活动时节,λ、
    Xm及Λ是计算过程的参数;

    步骤4.4:依据步骤3的中子能量范围和能量点和步骤4.3的通量计算结果,计算出
    所有航迹点随航迹点不同能量的大气中子通量值,利用如下公式计算航迹点的大气中子通
    量:


    其中:N(En)表示大气中子通量,En表示大气中子能量点;

    步骤5:获取待评估的空间器件的参数,包括σsat、LETth、W、S、t,空间器件参数由地
    面试验或者数值仿真获得,

    其中:σsat为重离子造成的器件单粒子事件效应发生概率的饱和截面,

    LETth、W、S为概率分布函数Weibull的系数值,

    t为电子器件的灵敏区厚度;

    步骤6:依据步骤3给出的中子能谱,利用蒙特卡洛计算工具计算给出不同能量中
    子的反冲核概率BGR(E,Eri)值,

    其中:E表示中子能量,

    Eri表示粒子在电子器件灵敏区内的沉积能量,

    BGR(E,Eri)表示能量为E的中子产生能量大于Eri的反冲核的概率;

    步骤7:依据步骤4获得参数,利用如下公式计算单粒子事件发生截面σsee(i),


    其中,LETion(i)表示重离子的LET值;

    步骤8:依据步骤6获得参数,利用如下公式计算单粒子事件发生概率,


    其中,C表示吸收系数,对于绝对大多数电子器件C取0.5时于结果符合实验验证,
    ΔVi参数表示为tΔσi=t(σi-σi-1)。

    较好地,本发明的方法还包括步骤9:采用Python或Matlab软件将计算结果可视
    化。

    根据本发明的一个实施例,所述航空飞行器的航迹取样点时间间隔不小于1个小
    时,否则将极大降低准确度。

    作为上述技术方案的一个实施例,所述输入的大气中能量范围,其能段低端小于
    等于10MeV、高端大于等于10MeV。

    根据本发明的一个实施例,第一步骤中的粒子输运仿真工具包括Geant4、MCNP及
    SRIM。

    根据本发明的一个实施例,所述航空飞行器结构物质包括金属铝、铝镁合金、钢制
    及复合材料,所述电子器件的灵敏物质包括硅、锗及砷化镓,所述航空器包括各类民用飞
    机、军用飞机、气球、导弹及高超音速飞行器。

    本发明的评估方法适合于地球近地范围的航空飞行器,所述航空器的飞行高度不
    高于100km,不低于0km,所述航空飞行器的飞行纬度为-90度到90度以内。

    本发明的优点在于:针对研究大气中子对航空飞行器电子器件发生单粒子事件效
    应的危害分布和防护策略、方法技术的需求,利用本发明的电子器件发生单粒子事件效应
    的概率的评估方法进行计算和分析,可以根据实测航空飞行器的航迹数据对器件发生单粒
    子事件效应的概率进行评估,从而便于在航空飞行器工程设计、故障诊断、飞行管理等工程
    阶段应用。

    附图说明

    图1为本发明的评估航空飞行电子器件发生单粒子事件效应的概率的方法的流程
    图;

    具体实施方式

    下面结合附图和优选实施例对本发明的一种用于评估航空飞行器电子器件发生
    单粒子事件效应的概率的方法进行详细说明,此实施例中以航空飞行器的结构物质铝(Al)
    来作为范例进行说明

    首先,步骤1:获得航空器在飞行过程的2012年所处太阳活动周期内所处阶段太阳
    活动参数为65;

    步骤2:获得航空器在飞行过程的航迹点的高度、纬度、经度数值序列,每5分钟选
    取一个航迹点;

    步骤3:设置大气中子的能量范围为1MeV到10GeV和能量点包括100个对数划分;

    步骤4:依据飞行时间、航迹点高度、纬度、经度数值序列及中子能谱设置,利用大
    气中子计算模型计算出不同航迹点处的大气中子密度分布、能谱分布;

    步骤4.1:依据航迹点纬度、经度数值序列,计算出地磁纬度、经度数值序列;

    Latm=90-arcos[sin(Latg)*sin(lat0)+cos(Latg)*cos(Lat0)*cos(long-lon0)]

    其中:Latm为航迹点地磁纬度,Latg为航迹点地理纬度,Long为航迹点地理经度,
    Lat0为地磁北极地理纬度,Lon0为地磁北极地理经度;

    步骤4.2:依据航迹点高度、地磁纬度、地磁经度数值序列,计算出航迹点处的磁场
    刚度值序列,利用如下公式:


    其中:Rc为磁刚度,x为航迹高度值,R为地球半径值;

    步骤4.3:依据步骤1、2的飞行时间、航迹纬度值及磁刚度序列,计算出所有航迹点
    大气中子的通量值,利用如下公式:

    λ=165+2*Rc

    Xm=50+ln{2000+exp[-x*(Cs-100)]}







    其中:x表示航迹高度值(单位:g/cm2),Cs为太阳活动参数,表征太阳活动时节,λ、
    Xm及Λ是计算过程的参数;

    步骤4.4:依据步骤3的中子能量范围和能量点和步骤4.3的通量计算结果,计算出
    所有航迹点随航迹点不同能量的大气中子通量值,利用如下公式计算航迹点的大气中子通
    量:


    其中:N(En)表示大气中子通量,En表示大气中子能量点;

    步骤5:获取待评估的空间器件的参数,包括σsat、LETth、W、S、t,空间器件参数由地
    面试验或者数值仿真获得,

    其中:σsat为重离子造成的器件单粒子事件效应发生概率的饱和截面,

    LETth、W、S为概率分布函数Weibull的系数值,

    t为电子器件的灵敏区厚度;

    步骤6:依据步骤3给出的中子能谱,利用蒙特卡洛计算工具计算给出不同能量中
    子的反冲核概率BGR(E,Eri)值,

    其中:E表示中子能量,

    Eri表示粒子在电子器件灵敏区内的沉积能量,

    BGR(E,Eri)表示能量为E的中子产生能量大于Eri的反冲核的概率;

    步骤7:依据步骤5获得参数,利用如下公式计算单粒子事件发生截面σsee(i),


    其中,LETion(i)表示重离子的LET值;

    步骤8:依据步骤6获得参数,利用如下公式计算单粒子事件发生概率,


    其中,C表示吸收系数,对于绝对大多数电子器件C取0.5时于结果符合实验验证,
    ΔVi参数表示为tΔσi=t(σi-σi-1)。

    较好地,本发明的实施案例还包括步骤9:采用Python或Matlab软件将计算结果可
    视化。

    本发明的评估方法适合于航空飞行器内各类电子半导体器件的单粒子事件效应
    发生概率的评估。此外,本发明的评估方法适合于各类近地航行器,包括各类民用飞机、军
    用飞机、气球、导弹及高超音速飞行器。

    最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参
    照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方
    案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明
    的权利要求范围当中。

    关 键 词:
    一种 大气 中子 航空 器件 粒子 效应 诱发 概率 评估 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    关于本文
    本文标题:一种大气中子对航空器件单粒子效应诱发概率的评估方法.pdf
    链接地址://www.4mum.com.cn/p-6092714.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03