CN115600317B - 一种载人航天器密封舱气体泄漏失效评估方法及系统 - Google Patents

一种载人航天器密封舱气体泄漏失效评估方法及系统 Download PDF

Info

Publication number
CN115600317B
CN115600317B CN202211267008.5A CN202211267008A CN115600317B CN 115600317 B CN115600317 B CN 115600317B CN 202211267008 A CN202211267008 A CN 202211267008A CN 115600317 B CN115600317 B CN 115600317B
Authority
CN
China
Prior art keywords
sealed cabin
mmod
perforation
aperture
gas
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202211267008.5A
Other languages
English (en)
Other versions
CN115600317A (zh
Inventor
武江凯
迟润强
韩增尧
庞宝君
胡迪奇
刘瑜妍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202211267008.5A priority Critical patent/CN115600317B/zh
Publication of CN115600317A publication Critical patent/CN115600317A/zh
Application granted granted Critical
Publication of CN115600317B publication Critical patent/CN115600317B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

本发明的一种载人航天器密封舱气体泄漏失效评估方法及系统，涉及一种风险评估方法和系统，目的是为了克服现有工程任务中对载人航天器失效的评估方法不够精准以及效率低的问题，方法具体包括步骤如下：步骤一、预测并计算孔径Dh；步骤二、计算密封舱临界穿孔孔径DhE；步骤三、判断孔径Dh是否大于等于密封舱临界穿孔孔径DhE；是则密封舱发生气体泄漏失效；且将失效次数NF累加1；步骤四、检测是否有Ni个MMOD穿孔完成步骤一～步骤三的计算；是则执行步骤五；否则返回步骤一，更换MMOD进行穿孔计算；步骤五、计算载人航天器密封舱气体泄漏失效概率PF＝NF/NS。

Description

一种载人航天器密封舱气体泄漏失效评估方法及系统
技术领域
本发明涉及一种风险评估方法和系统，具体涉及一种在MMOD环境下，载人航天器受到撞击失效的评估方法和系统。
背景技术
密封舱环境是载人航天活动中，为航天员提供必要的压力、温湿度以及宇宙射线防护等保障条件，是航天员工作和生活的重要场所。微流星体和空间碎片(Micro-Meteoroid&Orbit debris，MMOD)超高速撞击是大型载人航天器在轨长期飞行过程中面临的主要威胁，甚至可能导致载人航天器系统失效或航天员伤亡的灾难性失效。
舱内气体泄漏，造成航天员缺氧、失压，是密封舱穿孔后载人航天器失效、航天员伤亡的最主要的失效模式。当前工程任务中以“压力舱击穿”准则进行系统失效风险评估，即认为密封舱结构一旦被MMOD击穿载人航天器即失效，评估方法太过保守，评估结果也不够精准。另外，由于MMOD数量巨大，针对每一个MMOD撞击粒子分别开展穿孔孔径及TCE后总压、氧分压计算量大，评估效率低。
发明内容
本发明的目的是为了克服现有工程任务中对载人航天器失效的评估方法不够精准以及效率低的问题，提供了一种载人航天器密封舱气体泄漏失效评估方法及系统。
本发明提供一种载人航天器密封舱气体泄漏失效评估方法，具体包括步骤如下：
步骤一、预测并计算载人航天器的密封舱上空间碎片和微流星体MMOD穿孔的孔径Dh；
步骤二、计算不同MMOD穿孔的孔径Dh下，密封舱内氧分压和总压的变化规律；并根据临界逃生时间TCE以及保证航天员安全的总压PTE、氧分压
Figure BDA0003893478590000011
要求，计算密封舱临界穿孔孔径DhE；
步骤三、针对每个撞击碎片进行孔径Dh的计算，并判断孔径Dh是否大于等于密封舱临界穿孔孔径DhE；
是则判断在MMOD穿孔的孔径Dh下，密封舱发生气体泄漏失效；且将失效次数NF累加1；NF是初始值为0的自然数；
步骤四、检测是否有Ni个MMOD穿孔完成步骤一～步骤三的计算；Ni为与密封舱撞击的MMOD数量；
是则执行步骤五；否则返回步骤一，更换MMOD进行穿孔计算；
步骤五、计算载人航天器密封舱气体泄漏失效概率PF＝NF/NS，用以对载人航天器密封舱气体泄漏失效进行评估；
其中，NS为MMOD总数，且NS≥Ni。
本发明还提供一种载人航天器密封舱气体泄漏失效评估系统，包括：
MMOD穿孔计算模块，用于预测并计算载人航天器的密封舱上空间碎片和微流星体MMOD穿孔的孔径Dh，并发送至临界孔径计算模块和失效判断模块；
临界孔径计算模块，用于计算不同MMOD穿孔的孔径Dh下，密封舱内氧分压和总压的变化规律；并根据临界逃生时间TCE以及保证航天员安全的总压PTE、氧分压
Figure BDA0003893478590000021
要求，计算密封舱临界穿孔孔径DhE，并发送至失效判断模块；
失效判断模块，用于判断Dh是否大于等于DhE；
是则判断在MMOD穿孔的孔径Dh下，密封舱发生气体泄漏失效；且将失效次数NF累加1，并发送至评估概率计算模块；NF是初始值为0的自然数；并在判断结束后发送当前MMOD穿孔计算完成指令至循环检测模块；
循环检测模块，用于检测是否有Ni个MMOD穿孔完成步骤一～步骤三的计算；Ni为与密封舱撞击的MMOD数量；
是则将最终的失效次数NF发送至评估概率计算模块；否则更换MMOD穿孔，并发送至MMOD穿孔计算模块；
评估概率计算模块，用于计算载人航天器密封舱气体泄漏失效概率PF＝NF/NS，用以对载人航天器密封舱气体泄漏失效进行评估；
其中，NS为MMOD总数，且NS≥Ni。
本发明的有益效果是：
本文提出一种MMOD环境下大型载人航天器密封舱气体泄漏失效高效评估方法。
本发明的目的在于提供一种MMOD环境下大型载人航天器密封舱穿孔气体泄漏失效风险评估方法，实现对MMOD环境下密封舱穿孔失效风险进行有效、定量评估。通过求解满足航天员应急逃生的临界穿孔直径，将航天员缺氧失效判据从到达临界逃生时间时刻的总压、氧分压是否满足预设条件要求转化为MMOD撞击下穿孔孔径是否小于临界穿孔孔径，实现了在给定密封舱防护结构、舱内体积、初始总压和氧分压、初始温度、补气速率等密封舱参数下航天员缺氧和失压等2类失效的快速评估，解决大型载人航天器长期在轨飞行期间MMOD撞击下密封舱结构穿孔失效风险评估难题，为载人航天器防护结构设计、任务期间风险评估提供支撑。
附图说明
图1为本发明的MMOD穿孔的孔径的原理示意图；
图2为本发明的MMOD穿孔后密封舱内总压变化曲线图；
图3为本发明的MMOD穿孔后密封舱内氧分压变化曲线图；
图4为本发明的一种载人航天器密封舱气体泄漏失效评估方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
下面结合附图和具体实施例对本发明作进一步说明，但不作为本发明的限定。
具体实施方式一，本实施方式的一种载人航天器密封舱气体泄漏失效评估方法，具体包括步骤如下：
步骤一、预测并计算载人航天器的密封舱上空间碎片和微流星体MMOD穿孔的孔径Dh；
步骤二、计算不同MMOD穿孔的孔径Dh下，密封舱内氧分压和总压的变化规律；并根据临界逃生时间TCE以及保证航天员安全的总压PTE、氧分压
Figure BDA0003893478590000031
要求，计算密封舱临界穿孔孔径DhE；
步骤三、针对每个撞击碎片进行孔径Dh的计算，并判断孔径Dh是否大于等于密封舱临界穿孔孔径DhE；
是则判断在所述MMOD穿孔的孔径Dh下，密封舱发生气体泄漏失效；且将失效次数NF累加1；NF是初始值为0的自然数；
步骤四、检测是否有Ni个MMOD穿孔完成步骤一～步骤三的计算；Ni为与密封舱撞击的MMOD数量；
是则执行步骤五；否则返回步骤一，更换MMOD进行穿孔计算；
步骤五、计算载人航天器密封舱气体泄漏失效概率PF＝NF/NS，用以对载人航天器密封舱气体泄漏失效进行评估；
其中，NS为MMOD总数，且NS≥Ni。
最佳实施例，本实施例是对实施方式一的进一步说明，本实施例中，步骤二的具体步骤如下：
步骤二一、计算得到密封舱的MMOD穿孔处喷气速度vl和泄漏气体质量流量
Figure BDA0003893478590000041
Figure BDA0003893478590000042
Figure BDA0003893478590000043
其中，vl、Pl、ρl分别表示MMOD穿孔处气体流动速度、压力和密度，
Figure BDA0003893478590000044
表示MMOD穿孔处气体质量流量，γ表示绝热膨胀系数；
步骤二二、计算得到MMOD穿孔处气体压力Pl和密度ρl：
Figure BDA0003893478590000045
Figure BDA0003893478590000046
步骤二三、将公式(3)和公式(4)代入到公式(2)得到MMOD穿孔处气体质量流量为：
Figure BDA0003893478590000047
其中，Pi和ρi分别表示密封舱泄漏过程中第i时刻密封舱内气体压力和密度；则经过第i时刻到第i+1时刻的时间Δt，密封舱内气体总质量为
Figure BDA0003893478590000048
第i时刻舱内气体密度/>
Figure BDA0003893478590000049
其中，Vcabin为密封舱体积；
根据绝热等熵膨胀假设得到：
Figure BDA00038934785900000410
其中，Pi0和ρi0表示密封舱内气体初始压力和密度；并且Pi0、ρi0和Pi、ρi都满足理想气体的热状态方程；
步骤二四、根据公式(1)～(6)和热状态方程，得到密封舱Dh穿孔直径下，舱内氧分压和总压随时间变化规律，以及给定临界逃生时间TCE条件下密封舱临界穿孔直径DhE。
最佳实施例，本实施例是对实施方式一的进一步说明，本实施例中，步骤一的具体方法如下：
步骤一一、选取MMOD环境模型，针对载人航天器运行轨道、飞行姿态和飞行时间，获取各个撞击方向的MMOD通量数据，得到MMOD总数Ns；
步骤一二、针对每个MMOD与载人航天器的密封舱结构面元进行相交计算，记录撞击到密封舱的MMOD数量Ni；
步骤一三、针对撞击到密封舱的MMOD，计算MMOD穿孔的孔径Dh。
最佳实施例，本实施例是对实施方式一的进一步说明，本实施例中，步骤二中还包括：
判断是密封舱内的氧分压情况是否为氧分压紧急情况，
是则通过氧气瓶以设定的氧气补气速率对密封舱进行氧气补充，并根据氧气补气速率计算考虑补充氧气情况下，密封舱内氧分压的变化规律；
所述氧分压紧急情况为当氧分压下降到设定氧分压阈值以下的情况。
最佳实施例，本实施例是对实施方式一的进一步说明，本实施例中，步骤二中还包括：
判断是密封舱内的总压情况是否为总压紧急情况，
是则通过氮气瓶以设定的氮气补气速率对密封舱进行氮气补充，并根据氮气补气速率算考虑补充氮气情况下，密封舱内总压的变化规律；
所述总压紧急情况为当总压下降到设定总压阈值以下的情况。
具体实施方式二，本实施方式的一种载人航天器密封舱气体泄漏失效评估系统，包括：
MMOD穿孔计算模块，用于预测并计算载人航天器的密封舱上空间碎片和微流星体MMOD穿孔的孔径Dh，并发送至临界孔径计算模块和失效判断模块；
临界孔径计算模块，用于计算不同MMOD穿孔的孔径Dh下，密封舱内氧分压和总压的变化规律；并根据临界逃生时间TCE以及保证航天员安全的总压PTE、氧分压
Figure BDA0003893478590000051
要求，计算密封舱临界穿孔孔径DhE，并发送至失效判断模块；
失效判断模块，用于针对每个撞击碎片进行孔径Dh的计算，并判断孔径Dh是否大于等于密封舱临界穿孔孔径DhE；
是则判断在所述MMOD穿孔的孔径Dh下，密封舱发生气体泄漏失效；且将失效次数NF累加1，并发送至评估概率计算模块；NF是初始值为0的自然数；并在判断结束后发送当前MMOD穿孔计算完成指令至循环检测模块；
循环检测模块，用于检测是否有Ni个MMOD穿孔完成步骤一～步骤三的计算；Ni为与密封舱撞击的MMOD数量；
是则将最终的失效次数NF发送至评估概率计算模块；否则更换MMOD穿孔，并发送至MMOD穿孔计算模块；
评估概率计算模块，用于计算载人航天器密封舱气体泄漏失效概率PF＝NF/NS，用以对载人航天器密封舱气体泄漏失效进行评估；
其中，NS为MMOD总数，且NS≥Ni。
最佳实施例，本实施例是对实施方式二的进一步说明，本实施例中，临界孔径计算模块中采用如下步骤得到MMOD穿孔的孔径Dh下，密封舱内氧分压和总压的变化规律，以及给定临界逃生时间TCE条件下密封舱临界穿孔孔径DhE：
步骤二一、计算得到密封舱的MMOD穿孔处喷气速度vl和泄漏气体质量流量
Figure BDA0003893478590000061
Figure BDA0003893478590000062
Figure BDA0003893478590000063
其中，vl、Pl、ρl分别表示MMOD穿孔处气体流动速度、压力和密度，
Figure BDA0003893478590000064
表示MMOD穿孔处气体质量流量，γ表示绝热膨胀系数；
步骤二二、计算得到MMOD穿孔处气体压力Pl和密度ρl：
Figure BDA0003893478590000065
Figure BDA0003893478590000066
步骤二三、将公式(3)和公式(4)代入到公式(2)得到MMOD穿孔处气体质量流量为：
Figure BDA0003893478590000067
其中，Pi和ρi分别表示密封舱泄漏过程中第i时刻密封舱内气体压力和密度；则经过第i时刻到第i+1时刻的时间Δt，密封舱内气体总质量为
Figure BDA0003893478590000068
第i时刻舱内气体密度/>
Figure BDA0003893478590000069
其中，Vcabin为密封舱体积；
根据绝热等熵膨胀假设得到；
Figure BDA00038934785900000610
其中，Pi0和ρi0表示密封舱内气体初始压力和密度；并且Pi0、ρi0和Pi、ρi都满足理想气体的热状态方程；
步骤二四、根据公式(1)～(6)和热状态方程，得到密封舱Dh穿孔直径下，舱内氧分压和总压随时间变化规律，并根据临界逃生时间TCE以及保证航天员安全的总压PTE、氧分压
Figure BDA0003893478590000071
要求，计算密封舱临界穿孔直径DhE。
最佳实施例，本实施例是对实施方式二的进一步说明，本实施例中，MMOD穿孔计算模块包括：
MMOD总数计算模块，用于选取MMOD环境模型，针对载人航天器运行轨道、飞行姿态和飞行时间，获取各个撞击方向的MMOD通量数据，得到MMOD总数Ns，并发送至撞击MMOD数量计算模块；
撞击MMOD数量计算模块，用于针对每个MMOD与载人航天器的密封舱结构面元进行相交计算，记录撞击到密封舱的MMOD数量Ni，并发送至孔径计算模块；
孔径计算模块，用于针对撞击到密封舱的MMOD，计算MMOD穿孔的孔径Dh。
最佳实施例，本实施例是对实施方式二的进一步说明，本实施例中，临界孔径计算模块还包括：
氧分压情况判断计算模块，用于判断密封舱内的氧分压情况是否为氧分压紧急情况，
是则通过氧气瓶以设定的氧气补气速率对密封舱进行氧气补充，并根据氧气补气速率计算考虑补充氧气情况下，密封舱内氧分压的变化规律；
所述氧分压紧急情况为当氧分压下降到设定氧分压阈值以下的情况。
最佳实施例，本实施例是对实施方式二的进一步说明，本实施例中，临界孔径计算模块还包括：
总压情况判断计算模块，用于判断是密封舱内的总压情况是否为总压紧急情况，
是则通过氮气瓶以设定的氮气补气速率对密封舱进行氮气补充，并根据氮气补气速率算考虑补充氮气情况下，密封舱内总压的变化规律；
所述总压紧急情况为当总压下降到设定总压阈值以下的情况。
具体地，本发明提供一种大型载人航天器密封舱穿孔气体泄漏失效风险评估方法，通过计算MMOD环境下密封舱穿孔孔径并与满足航天员临界逃生时间的临界穿孔直径比对，可以实现对大型载人航天器长期在轨期间，MMOD环境下密封舱结构穿孔模式下，航天员安全性进行风险评估，支持防护结构设计和任务风险评估实施。包括：
1)选取空间碎片和微流星体(Micro-Meteoroid&Orbitdebris，以下简称MMOD)环境模型，针对大型载人航天器运行轨道、飞行姿态和飞行时间，获取各个撞击方向的通量数据，记录MMOD总数NS；
2)采用面元法或射击线法，针对每个MMOD与载人航天器密封舱结构面元进行相交计算，记录撞击到舱体的MMOD数量Ni；
3)针对撞击到舱体的MMOD，采用载人航天器密封舱结构穿孔孔径经验公式(可以使用W-S穿孔方程)，计算穿孔孔径大小；
4)失效判据
(1)定义DhE为引起航天员缺氧失效的临界穿孔尺寸，当密封舱穿孔直径Dh小于DhE情况下，到达临界逃生时间TCE时，密封舱内总压和氧分压分别大于PTE、
Figure BDA0003893478590000081
航天员不会产生缺氧失效；当Dh大于DhE情况下，TCE时刻密封舱内总压小于PTE或氧分压小于/>
Figure BDA0003893478590000082
则航天员产生缺氧失效。因此，通过求解满足航天员应急逃生的临界穿孔直径DhE，将航天员缺氧失效判据从TCE时刻的总压、氧分压是否满足PTE、/>
Figure BDA0003893478590000083
要求转化为MMOD撞击下穿孔尺寸是否小于DhE。
5)临界穿孔直径计算
假设MMOD撞击作用下，密封舱舱壁被击穿后，舱壁产生孔径为Dh的通孔，如图1所示。
为简化气体泄漏分析模型，假设密封舱内空气为理想气体，气体泄漏速度垂直于该通孔切面，且沿泄漏方向气体的温度、速度、气压及密度相同。此外，考虑舱内空气初始压力和流速较高，泄漏时间较短，因此忽略泄漏气体与舱体之间热交换，将气体泄漏过程近似为绝热等熵流动，采用一维拉瓦尔喷管模型描述舱壁通孔处的气体泄漏过程。
由于穿孔处气体流动参数取决于孔外压力Po(对于运行在轨道上大型载人航天器Po＝0Pa)和舱内环境气体压力Pi之比，当Po/Pi小于0.5283时，气体泄漏速度达到当地音速(即马赫数Ma＝1)，根据气体连续性方程，得到密封舱穿孔处喷气速度和泄漏气体质量流量如下：
Figure BDA0003893478590000084
Figure BDA0003893478590000085
其中vl、Pl、ρl分别表示穿孔处气体流动速度、压力和密度，
Figure BDA0003893478590000086
表示穿孔处气体质量流量，γ表示绝热膨胀系数(对于空气，γ＝1.4)。
假设在气体泄漏过程中，航天员耗氧及舱内正常补气速率恒定，忽略密封舱内气体流动和气体的正常流出，根据理想气体的绝热等熵膨胀假设，得到穿孔处气体压力和密度如式(3)和式(4)所示：
Figure BDA0003893478590000091
Figure BDA0003893478590000092
代入到式(2)得到穿孔处气体质量流量为：
Figure BDA0003893478590000093
Pi和ρi分别表示密封舱泄漏过程中第i时刻舱内气体压力和密度，
则经过时间Δt密封舱内气体总质量
Figure BDA0003893478590000094
第i时刻舱内气体密度
Figure BDA0003893478590000095
其中Vcabin为密封舱体积。
根据绝热等熵膨胀假设，由式(6)计算得到。
Figure BDA0003893478590000096
Pi0和ρi0表示密封舱内气体初始压力和密度，在密封舱体积Vcabin一定情况下，Pi0、ρi0和Pi、ρi都满足理想气体的热状态方程。根据以上式(1)～(6)，则得到舱内压力变化规律以及给定TCE条件下密封舱临界穿孔直径DhE。
例如，密封舱体积为36m3，初始温度为21℃，初始总压为97.0kPa，初始氧分压为20.6kPa，舱外空间环境气压为1.01E-5Pa，紧急情况下氧气瓶和氮气瓶补气速率分别为12g/s和18g/s，3名航天员平均耗氧速率为0.01g/s/人，当氧分压下降到20.0kPa时，供氧组件(氧气瓶)开启，当总压下降到87.0kPa时，供氮组件(氮气瓶)开启，满足航天员紧急在轨堵漏或者撤离总压PTE为70kPa，氧分压为
Figure BDA0003893478590000097
为14kPa，得到不同穿孔孔径下，密封舱内氧分压和总压变化如图2和图3所示，设定保证密封舱内航天员紧急堵漏或撤离时间1200s，则利用式(3)～(6)及给定的补气速率(包括氧气补气速率和氮气补气速率)下，得到密封舱结构不发生灾难性失效的临界穿孔孔径为1.54cm，即当MMOD穿孔孔径大于1.54cm时，航天员发生缺氧失效。
6)重复进行第2)、3)、4)，完成所有与舱体结构撞击的MMOD(共Ni个)的舱体穿孔直径Dh计算，与6)临界穿孔直径DhE比较，若穿孔直径Dh≥DhE，则系统失效事件次数NF＝NF+1。
7)所有撞击事件中系统失效事件次数NF，系统失效概率为PF＝NF/NS。
并且，在开展MMOD环境下大型载人航天器密封舱结构穿孔失效气体泄漏风险评估前，应首先完成以下数据准备：
(1)选取空间碎片和微流星体(Micro-Meteoroid&Orbit debris，以下简称MMOD)环境模型；
(2)载人航天器密封舱结构(含防护结构)穿孔经验方程，推荐使用W-S方程；
(3)根据航天员在轨安全总压PTE和安全氧分压
Figure BDA0003893478590000101
和临界逃生时间TCE要求，求解密封舱结构临界穿孔直径DhE；可以事先求得不同MMOD穿孔孔径下，密封舱内的总压和氧分压随时间的变化规律。然后再通过拟合的方式，得到一条在TCE时刻对应的总压和氧分压均保持在安全总压PTE和安全氧分压/>
Figure BDA0003893478590000102
以上(或重合)的曲线，根据曲线得到所对应的穿孔孔径，即为临界穿孔孔径DhE。
2)针对大型载人航天器运行轨道、飞行姿态和飞行时间，获取各个撞击方向的通量数据，记录MMOD总数NS；
3)采用面元法或射击线法，针对每个MMOD与载人航天器密封舱结构面元进行相交计算，记录撞击到舱体的MMOD数量Ni；
4)针对撞击到舱体的MMOD，采用载人航天器密封舱结构穿孔孔径经验公式(推荐使用W-S穿孔方程)，计算穿孔孔径大小Dh；
5)若Dh≥DhE，则本次撞击穿孔事件引起系统失效，系统失效次数NF加1，否则，系统安全无影响；
10)对所有射击线依次进行计算，得到引起系统失效总撞击次数NF，系统失效概率为PF＝NF/NS。
综上，通过本发明，可以解决大型载人航天器长期在轨飞行期间，MMOD撞击下密封舱结构穿孔失效风险评估难题，为航天器防护设计以及在轨任务规划提供支撑，并可推广应用于其他影响航天器安全的空间环境或空间主动威胁源作用下，航天器、航天器易损性评估。
虽然在本文中参照了特定的实施方式来描述本发明，但是应该理解的是，这些实施例仅仅是本发明的原理和应用的示例。因此应该理解的是，可以对示例性的实施例进行许多修改，并且可以设计出其他的布置，只要不偏离所附权利要求所限定的本发明的精神和范围。应该理解的是，可以通过不同于原始权利要求所描述的方式来结合不同的从属权利要求和本文中的特征。还可以理解的是，结合单独实施例所描述的特征可以使用在其他实施例中。

Claims (8)

1.一种载人航天器密封舱气体泄漏失效评估方法，其特征在于，具体包括步骤如下：
步骤一、预测并计算载人航天器的密封舱上空间碎片和微流星体MMOD穿孔的孔径Dh；
步骤二、计算不同MMOD穿孔的孔径Dh下，密封舱内氧分压和总压的变化规律；并根据临界逃生时间TCE以及保证航天员安全的总压PTE、氧分压
Figure FDA0004223594380000011
要求，计算密封舱临界穿孔孔径DhE；
步骤三、针对每个撞击碎片进行孔径Dh的计算，并判断孔径Dh是否大于等于密封舱临界穿孔孔径DhE；
是则判断在所述MMOD穿孔的孔径Dh下，密封舱发生气体泄漏失效；且将失效次数NF累加1；NF是初始值为0的自然数；
步骤四、检测是否有Ni个MMOD穿孔完成步骤一～步骤三的计算；Ni为与密封舱撞击的MMOD数量；
是则执行步骤五；否则返回步骤一，更换MMOD进行穿孔计算；
步骤五、计算载人航天器密封舱气体泄漏失效概率PF＝NF/NS，用以对载人航天器密封舱气体泄漏失效进行评估；
其中，NS为MMOD总数，且NS≥Ni；
步骤二的具体步骤如下：
步骤二一、计算得到密封舱的MMOD穿孔处喷气速度vl和泄漏气体质量流量
Figure FDA0004223594380000012
Figure FDA0004223594380000013
Figure FDA0004223594380000014
其中，vl、Pl、ρl分别表示MMOD穿孔处气体流动速度、压力和密度，
Figure FDA0004223594380000015
表示MMOD穿孔处气体质量流量，γ表示绝热膨胀系数；
步骤二二、计算得到MMOD穿孔处气体压力Pl和密度ρl：
Figure FDA0004223594380000016
Figure FDA0004223594380000017
步骤二三、将公式(3)和公式(4)代入到公式(2)得到MMOD穿孔处气体质量流量为：
Figure FDA0004223594380000021
其中，Pi和ρi分别表示密封舱泄漏过程中第i时刻密封舱内气体压力和密度；则经过第i时刻到第i+1时刻的时间Δt，密封舱内气体总质量为
Figure FDA0004223594380000022
第i时刻舱内气体密度/>
Figure FDA0004223594380000023
其中，Vcabin为密封舱体积；
并且根据绝热等熵膨胀假设得到；
Figure FDA0004223594380000024
其中，Pi0和ρi0表示密封舱内气体初始压力和密度；并且Pi0、ρi0和Pi、ρi都满足理想气体的热状态方程；
步骤二四、根据公式(1)～(6)和热状态方程，得到密封舱Dh穿孔直径下，舱内氧分压和总压随时间变化规律，以及给定临界逃生时间TCE条件下密封舱临界穿孔直径DhE；
事先求得不同MMOD穿孔孔径下，密封舱内的总压和氧分压随时间的变化规律；然后再通过拟合的方式，得到一条在TCE时刻对应的总压和氧分压均保持在安全总压PTE和安全氧分压
Figure FDA0004223594380000025
以上或重合的曲线，根据曲线得到所对应的穿孔孔径，即为临界穿孔孔径DhE。
2.根据权利要求1所述的一种载人航天器密封舱气体泄漏失效评估方法，其特征在于，步骤一的具体方法如下：
步骤一一、选取MMOD环境模型，针对载人航天器运行轨道、飞行姿态和飞行时间，获取各个撞击方向的MMOD通量数据，得到MMOD总数Ns；
步骤一二、针对每个MMOD与载人航天器的密封舱结构面元进行相交计算，记录撞击到密封舱的MMOD数量Ni；
步骤一三、针对撞击到密封舱的MMOD，计算MMOD穿孔的孔径Dh。
3.根据权利要求2所述的一种载人航天器密封舱气体泄漏失效评估方法，其特征在于，步骤二中还包括：
判断密封舱内的氧分压情况是否为氧分压紧急情况，
是则通过氧气瓶以设定的氧气补气速率对密封舱进行氧气补充，并根据氧气补气速率计算考虑补充氧气情况下，密封舱内氧分压的变化规律；
所述氧分压紧急情况为当氧分压下降到设定氧分压阈值以下的情况。
4.根据权利要求3所述的一种载人航天器密封舱气体泄漏失效评估方法，其特征在于，步骤二中还包括：
判断密封舱内的总压情况是否为总压紧急情况，
是则通过氮气瓶以设定的氮气补气速率对密封舱进行氮气补充，并根据氮气补气速率算考虑补充氮气情况下，密封舱内总压的变化规律；
所述总压紧急情况为当总压下降到设定总压阈值以下的情况。
5.一种载人航天器密封舱气体泄漏失效评估系统，其特征在于，包括：
MMOD穿孔计算模块，用于预测并计算载人航天器的密封舱上空间碎片和微流星体MMOD穿孔的孔径Dh，并发送至临界孔径计算模块和失效判断模块；
临界孔径计算模块，用于计算不同MMOD穿孔的孔径Dh下，密封舱内氧分压和总压的变化规律；并根据临界逃生时间TCE以及保证航天员安全的总压PTE、氧分压
Figure FDA0004223594380000031
要求，计算密封舱临界穿孔孔径DhE，并发送至失效判断模块；
失效判断模块，用于针对每个撞击碎片进行孔径Dh的计算，并判断孔径Dh是否大于等于密封舱临界穿孔孔径DhE；
是则判断在所述MMOD穿孔的孔径Dh下，密封舱发生气体泄漏失效；且将失效次数NF累加1，并发送至评估概率计算模块；NF是初始值为0的自然数；并在判断结束后发送当前MMOD穿孔计算完成指令至循环检测模块；
循环检测模块，用于检测是否有Ni个MMOD穿孔完成步骤一～步骤三的计算；Ni为与密封舱撞击的MMOD数量；
是则将最终的失效次数NF发送至评估概率计算模块；否则更换MMOD穿孔，并发送至MMOD穿孔计算模块；
评估概率计算模块，用于计算载人航天器密封舱气体泄漏失效概率PF＝NF/NS，用以对载人航天器密封舱气体泄漏失效进行评估；
其中，NS为MMOD总数，且NS≥Ni；
临界孔径计算模块中采用如下步骤得到MMOD穿孔的孔径Dh下，密封舱内氧分压和总压的变化规律，以及给定临界逃生时间TCE条件下密封舱临界穿孔孔径DhE：
步骤二一、计算得到密封舱的MMOD穿孔处喷气速度vl和泄漏气体质量流量
Figure FDA0004223594380000032
Figure FDA0004223594380000033
Figure FDA0004223594380000034
其中，vl、Pl、ρl分别表示MMOD穿孔处气体流动速度、压力和密度，
Figure FDA0004223594380000041
表示MMOD穿孔处气体质量流量，γ表示绝热膨胀系数；
步骤二二、计算得到MMOD穿孔处气体压力Pl和密度ρl：
Figure FDA0004223594380000042
Figure FDA0004223594380000043
步骤二三、将公式(3)和公式(4)代入到公式(2)得到MMOD穿孔处气体质量流量为：
Figure FDA0004223594380000044
其中，Pi和ρi分别表示密封舱泄漏过程中第i时刻密封舱内气体压力和密度；则经过第i时刻到第i+1时刻的时间Δt，密封舱内气体总质量为
Figure FDA0004223594380000045
第i时刻舱内气体密度/>
Figure FDA0004223594380000046
其中，Vcabin为密封舱体积；
根据绝热等熵膨胀假设得到；
Figure FDA0004223594380000047
其中，Pi0和ρi0表示密封舱内气体初始压力和密度；并且Pi0、ρi0和Pi、ρi都满足理想气体的热状态方程；
步骤二四、根据公式(1)～(6)和热状态方程，得到密封舱Dh穿孔直径下，舱内氧分压和总压随时间变化规律，以及给定临界逃生时间TCE条件下密封舱临界穿孔直径DhE；
事先求得不同MMOD穿孔孔径下，密封舱内的总压和氧分压随时间的变化规律；然后再通过拟合的方式，得到一条在TCE时刻对应的总压和氧分压均保持在安全总压PTE和安全氧分压
Figure FDA0004223594380000048
以上或重合的曲线，根据曲线得到所对应的穿孔孔径，即为临界穿孔孔径DhE。
6.根据权利要求5所述的一种载人航天器密封舱气体泄漏失效评估系统，其特征在于，MMOD穿孔计算模块包括：
MMOD总数计算模块，用于选取MMOD环境模型，针对载人航天器运行轨道、飞行姿态和飞行时间，获取各个撞击方向的MMOD通量数据，得到MMOD总数Ns，并发送至撞击MMOD数量计算模块；
撞击MMOD数量计算模块，用于针对每个MMOD与载人航天器的密封舱结构面元进行相交计算，记录撞击到密封舱的MMOD数量Ni，并发送至孔径计算模块；
孔径计算模块，用于针对撞击到密封舱的MMOD，计算MMOD穿孔的孔径Dh。
7.根据权利要求6所述的一种载人航天器密封舱气体泄漏失效评估系统，其特征在于，临界孔径计算模块还包括：
氧分压情况判断计算模块，用于判断密封舱内的氧分压情况是否为氧分压紧急情况，
是则通过氧气瓶以设定的氧气补气速率对密封舱进行氧气补充，并根据氧气补气速率计算考虑补充氧气情况下，密封舱内氧分压的变化规律；
所述氧分压紧急情况为当氧分压下降到设定氧分压阈值以下的情况。
8.根据权利要求7所述的一种载人航天器密封舱气体泄漏失效评估系统，其特征在于，临界孔径计算模块还包括：
总压情况判断计算模块，用于判断密封舱内的总压情况是否为总压紧急情况，
是则通过氮气瓶以设定的氮气补气速率对密封舱进行氮气补充，并根据氮气补气速率算考虑补充氮气情况下，密封舱内总压的变化规律；
所述总压紧急情况为当总压下降到设定总压阈值以下的情况。
CN202211267008.5A 2022-10-17 2022-10-17 一种载人航天器密封舱气体泄漏失效评估方法及系统 Active CN115600317B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211267008.5A CN115600317B (zh) 2022-10-17 2022-10-17 一种载人航天器密封舱气体泄漏失效评估方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211267008.5A CN115600317B (zh) 2022-10-17 2022-10-17 一种载人航天器密封舱气体泄漏失效评估方法及系统

Publications (2)

Publication Number Publication Date
CN115600317A CN115600317A (zh) 2023-01-13
CN115600317B true CN115600317B (zh) 2023-06-20

Family

ID=84846435

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211267008.5A Active CN115600317B (zh) 2022-10-17 2022-10-17 一种载人航天器密封舱气体泄漏失效评估方法及系统

Country Status (1)

Country Link
CN (1) CN115600317B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1191588A2 (en) * 2000-09-26 2002-03-27 Yukihiro Hirai A spiral contactor and manufacturing method therefor
RU2213943C2 (ru) * 2001-05-28 2003-10-10 Открытое акционерное общество "Ракетно-космическая корпорация "Энергия" им. С.П.Королева" Способ определения герметичности изолированного объема космического аппарата в условиях космического полета
US9396294B1 (en) * 2013-08-27 2016-07-19 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Method of modeling and simulation of shaped external occulters
CN107657069A (zh) * 2016-07-26 2018-02-02 北京空间技术研制试验中心 载人航天器舱内流场与温度场松耦合计算方法
CA3086756A1 (en) * 2019-08-28 2021-02-28 The Boeing Company Additively manufactured spacecraft panel

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009260176A (ja) * 2008-04-21 2009-11-05 Nec Electronics Corp 半導体装置の信頼性予測方法、及びそのプログラム
CN102507106B (zh) * 2011-11-21 2014-12-31 北京卫星环境工程研究所 一种航天器在轨状态泄漏监测的方法
RU2741401C1 (ru) * 2020-01-29 2021-01-25 Андрей Иванович Шумейко Модуль с многоканальной плазменной двигательной установкой для малого космического аппарата
CN112182757B (zh) * 2020-09-27 2022-11-08 北京空间飞行器总体设计部 空间碎片和微流星体环境下航天器易损性评估方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1191588A2 (en) * 2000-09-26 2002-03-27 Yukihiro Hirai A spiral contactor and manufacturing method therefor
RU2213943C2 (ru) * 2001-05-28 2003-10-10 Открытое акционерное общество "Ракетно-космическая корпорация "Энергия" им. С.П.Королева" Способ определения герметичности изолированного объема космического аппарата в условиях космического полета
US9396294B1 (en) * 2013-08-27 2016-07-19 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Method of modeling and simulation of shaped external occulters
CN107657069A (zh) * 2016-07-26 2018-02-02 北京空间技术研制试验中心 载人航天器舱内流场与温度场松耦合计算方法
CA3086756A1 (en) * 2019-08-28 2021-02-28 The Boeing Company Additively manufactured spacecraft panel

Also Published As

Publication number Publication date
CN115600317A (zh) 2023-01-13

Similar Documents

Publication Publication Date Title
Hughes et al. Inflatable re-entry vehicle experiment (IRVE) design overview
CN115600317B (zh) 一种载人航天器密封舱气体泄漏失效评估方法及系统
Silverman et al. Development of the crew dragon ECLSS
Poteet et al. Preliminary thermal-mechanical sizing of metallic TPS-Process development and sensitivity studies
Hamlin et al. 2009 Space Shuttle Probabilistic Risk Assessment Overview
Ocampo Limitations of spacecraft redundancy: A case study analysis
Tobias et al. Evolution of Hardware and Philosophy of Emergency Response Actions on the International Space Station and Future Spacecrafts
Christiansen MMOD protection and degradation effects for thermal control systems
McNamara Entry Atmospheric Flight Control Authority Impacts on GN&C and Trajectory Performance for Orion Exploration Flight Test 1
Maddock et al. Passive vs. parachute system architecture for robotic sample return vehicles
Handley A Pilot Model for the NASA Simplified Aid for EVA Rescue (SAFER)(Single-Axis Pitch Task)
Hirsch Current Trends on the Applicability of Ground Aerospace Materials Test Data to Space System Environments
Mitcheltree et al. An Earth Entry Vehicle for Returning Samples from Mars
Boyer et al. Space Shuttle Probabilistic Risk Assessment Overview
Klinkrad Meteoroid and debris protection
GAMBLE et al. The development and application of aerodynamic uncertainties in the design of the entry trajectory and flight control system of the SpaceShuttle Orbiter
Edquist Correlation of radiative heating calculations for Venus entry
Skeen Conceptual modeling and analysis of drag-augmented supersonic retropropulsion for application in Mars entry, descent, and landing vehicles
Blue Human Spaceflight Accidents: The US Space Program
Datta The application of NASA risk management to the SOFIA program
Johnston et al. Aerothermodynamics of the ARD-Postflight numerics and shock-tunnel experiments
Williamsen et al. Effect of multi-wall system composition on survivability for spacecraft impacted by orbital debris
Celestian et al. Dynamic response of the space station freedom due to a module perforation by a hypervelocity impact
Perino et al. A new era and a new tradespace: evaluating Earth entry vehicles concepts for a potential 2026 Mars Sample Return
Hanson et al. Launch vehicle failure dynamics and abort triggering analysis

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant