钟形返回舱空气动力特性数值模拟

文章基于Euler方程及N_S方程的数值求解方法 ,对返回舱亚、跨声速和高超声速的流场及气动特性进行了数值模拟 ,其中Euler方程数值求解采用二阶Godunov有限体积法 ;N_S方程数值求解采用二阶Harten_Yee格式的差分法。得到同实验值一致的物面压力、气动力系数和在不同速度范围出现的激波、流动分离及旋涡等流场特征。通过完全气体、平衡气体和非平衡气体的流场数值模拟结果分析比较 ,得出真实气体效应对返回舱气动力特性影响较小这一结论。计算结果表明数值模拟方法是预测返回舱气动特性的有效手段     

钟形返回舱空气动力特性数值模拟

文章基于Euler方程及N－S方程的数值求解方法，对回舱亚、跨声速和高超声速的流场及气动特性进行了数值模拟，其中Euler方程数值求解采用二阶Godunov有限体积法；N－S方程数值求解采用二阶Harten-Yee格式的差分法，得到同实验值一致的物面压力、气动力系数和在不同速度范围出现的激波、流动分离及旋涡等流场特征。通过完全气体、平衡气体和非平衡气体的流场数值模拟结果分析比较，得出真实气体效应对返回舱气动力特性影响较小这一结论。计算结果表明数值模拟方法是预测返回舱气动特性的有效手段。     

高超声速导弹多场耦合仿真

为准确预测高超声速导弹气动力与气动热环境,采用分区求解方法,实现气动力/热/结构多场耦合计算,并通过与试验值对比确立了耦合仿真方法的有效性。建立大长细比高超声速导弹仿真计算模型,对导弹弹体结构变形、温度场和压力场进行了数值模拟。仿真结果表明：弹体结构随着攻角增加发生轴向拉伸和横向弯曲;拉伸变形主要由气动热引起,弯曲变形主要由气动力引起;导弹头部气动热与气动力载荷大;耦合效应随着导弹攻角增加更加明显,造成导弹气动热力学环境更加严酷。     

气动热作用下的充气式减速器性能研究

为了解高超声速再入时气动热载荷对充气式减速器柔性结构的影响,文章基于松散耦合方法开展了极端热载荷工况下的耦合数值研究。文章首先建立了流固耦合和热固耦合两种模型,分别对比研究了气动力和气动热两种气动载荷对蒙皮结构的影响。结果表明,气动热对结构的影响远大于气动力,在高超声速再入时应重点考虑。之后研究了气动热载荷下充气式减速器防热层各功能层温度分布,结果表明,绝热层隔热效果最为显著,绝热层导热系数增大一倍,内部最高温度升高21.7%,热变形最大值升高10.7%。上述成果为充气式减速器的设计提供了一定的理论依据。     

飞行器结构部件导热/辐射耦合传热特性预测方法

采用数值方法求解气动热环境,基于蒙特卡罗方法求解内部空腔辐射换热及有限元方法求解三维导热,建立了沿弹道求解导热/辐射耦合的热响应预测方法,计算并获得了沿给定飞行弹道条件下的考虑内部辐射和不考虑内部辐射时的某高超声速飞行器结构部件的热响应特性。研究表明,所发展的耦合计算方法具有较高的精度和较好的工程适用性;考虑内部辐射时,结构部件局部最高温度明显低于不考虑辐射时,最高相差400K以上,且温度分布趋于均匀,温度梯度减小。相关研究对高超声速飞行器防热结构设计与优化具有重要参考意义。     

高超声速飞行器主动气膜冷却热防护数值仿真研究

为研究高超声速飞行器头部全覆盖保护情况下的气动热分布特征,文章提出一种微孔射流的主动气膜热防护方案,并对射流微孔分布进行优化;通过数值求解N-S方程,得到高超声速飞行器头部的驻点压力及表面、附近温度分布。研究结果表明,在主动气膜冷却热防护下,高超声速飞行器壁面温度可以降到1000 K以下。该方案可为未来高超声速飞行器的外壳设计提供参考。     

基于面元法的栅格翼翼身组合体高超声速气动特性研究

栅格翼是由外部框架和内部栅格布置形成的空间多升力面系统,其在高超声速状态具有突出优点。文中结合高超声速流动特点,考虑头部激波及弹身上洗干扰,基于面元法建立栅格翼翼身组合体高超声速气动特性计算模型,根据研究对象特点自动划分面元,进而完成典型栅格翼翼身组合体高超声速气动特性计算,并与数值计算结果进行对比分析。结果表明,在研究范围内,栅格翼翼身组合体高超声速气动特性面元法计算结果与数值计算结果吻合较好,前者较后者相对误差绝对值随攻角增加而增大。     

吸气式高超声速飞行器动力学特性分析

吸气式高超声速飞行器是下一代单级入轨和高超声速巡航研究的重点飞行器。其机身/发动机组合对飞行器动力学影响不同于常规飞行器,不能通过常规的气动分析方法对整体飞行器进行气动力计算分析。提出了利用激波膨胀波理论计算飞行器上表面的气动力,利用二维流场中的牛顿激波理论计算飞行器前体下表面的气动力,采用一维定常流模拟发动机内部工作流场,膨胀波理论计算飞行器后体气动力的吸气式高超声速飞行器整体气动分析方法。仿真分析验证了本方法的可行性和正确性。     

火星探测器表面材料催化特性数值模拟研究

针对火星探测器高超声速进入过程中的表面材料催化作用及其对气动热环境影响这一问题,建立了变壁面温度的火星大气表面材料催化作用模型,并基于火星大气物理化学模型和求解三维热化学非平衡N-S方程的数值方法,对典型火星探测器防热大底进行了数值模拟,获得了不同催化特性下的高超声速非平衡流场和气动热数据,分析了表面材料催化特性对气动热环境影响的规律性。研究结果表明：表面催化特性对壁面附近组分分布影响很大,催化反应进程主要受O原子浓度限制;有限催化热流随催化效率增大而增大,完全催化峰值热流比催化效率为1的有限催化峰值热流高25%~64%;表面温度随催化特性的变化规律与热流变化规律类似。有限催化模型能根据表面材料的催化特性精细化预测表面热流和温度,为防热设计提供更精确合理的参考标准。     

壁面催化条件对热环境预测的影响

基于高超声速飞行器热流精确预测的需要,文章采用数值模拟的方法研究了高超声速热化学非平衡流动中壁面催化条件对热环境的影响.针对ram-cⅡ气动外形,使用ausm+up格式耦合lu-sgs方法求解带化学反应源项的多组分ns方程.得到如下结论:有限催化壁结果介于完全催化壁以及非催化壁之间:壁面催化特性的差异主要体现在热流的组分扩散项部分,因此可以通过评估完全催化壁时组分扩散项所占比率大小来预测壁面热流准确值.     

用样条有限点法分析梁和薄板的弹塑性弯曲

本文首次用样条有限点法(SFPM),分析了梁和薄板的弹塑性弯曲,并给出了计算实例。通过计算表明:SFPM具有方法简单、节约内存、计算精度高等优点。对于规则区域的特殊问题,SFPM做为有限元法的一种补充,是一种有效的数值分析方法。     

高超声速飞行器材料与结构气动热环境模拟方法及试验研究

文章介绍了自行研制的石英灯红外辐射式气动加热试验模拟系统以及使用该系统对高超声速飞行器材料与结构进行的高温热评价试验。本热试验系统可实现升温速率高至200 ℃/s的非线性热冲击过程的动态模拟;能够生成1.8 MW/m2热流密度的瞬态非线性热试验模拟环境;能将试验环境温度提高到1 500 ℃。在该热试验系统上完成了如下试验研究: 1）金属蜂窝板结构在高温950 ℃非线性热环境下的隔热性能评价试验和数值模拟;2）对SiC/SiC复合材料试件在1 300～1 500 ℃下的隔热性能评价试验;3）采用轴向非分段加热试验方式对圆柱型壳体结构（长2.1 m）内壁进行高温热环境试验。本试验系统在可控的非线性温升速率、高温高热流密度变化过程的动态模拟、热试验环境模拟的准确性以及非接触式全场高温变形测量等方面的研究成果达到了国际先进水平。     

进入火星大气的高温真实气体效应与气动加热研究

针对火星和地球大气分子热力学和化学行为的差异性,采用理论分析和数值模拟两种手段,研究探测器进入过程高超声速流动的分子振动激发、离解反应及热力学和化学非平衡等真实气体效应,获得不同气体模型条件下的高超声速气动加热规律,探究引起地火差异的根本原因。分析认为,探测器进入火星大气层的稀薄气体效应明显;激波层内发生CO 2气体为主的大规模离解,在极高温环境下O 2和CO也将离解;沿进入轨道的高超声速流动基本处于化学非平衡但热力学平衡状态;激波层内能量储存和分配模式因分子振动激发和化学反应而改变,分子振动激发会增强气动加热量,但均介于化学反应模型的完全非催化和完全催化壁结果之间;相同来流条件下CO 2介质高超声速气动加热强于空气介质,但真实的火星进入热载荷因大气稀薄而弱于地球再入环境。相关研究为我国未来火星探测器热防护系统设计提供技术支持。     

对接过程中捕获锁质量和锁合面间隙影响分析

采用弹性接触动力学有限元和多刚体动力学方法分别对航天器对接接触碰撞过程和对接机构捕获、缓冲与校正过程进行了数值模拟,对比分析了捕获锁质量对捕获过程的影响及捕获锁锁合面间隙对捕获后对接过程的影响。从仿真结果可以看出,附加质量有利于主、被动对接环之间形成接触滑移,捕获后锁合面间隙引起两飞行器间相对速度震荡,并延迟了碰撞力和缓冲时间。     

热防护系统分区协调耦合推进方法

提出一种适用于热防护系统(TPS)热控性能研究的分区协调耦合推进方法,其中采用有限体积法(FVM)进行气动热分析,FVM空间离散采用NND格式,而结构传热采用有限元法(FEM)进行分析,且在耦合面采用基于控制面的双向映射插值方法进行数据传递。进行了圆管算例分析,2 s时刻驻点处温度计算值与试验值相对误差为4.95%。研究了空天飞行器头锥TPS的热控性能,非耦合方法获得的防热瓦和应变隔离垫(SIP)最高温度分别比耦合结果高114.4 K和32.6 K,这是由于非耦合方法未考虑壁面温度升高对气动热的反馈作用,而耦合方法充分考虑了此影响。采用高热辐射率的涂层、低导热系数和较厚的防热瓦能有效提高热防护系统的隔热性能和降低主动冷却系统的功率和重量,而防热瓦最高温度对其导热系数和厚度不敏感。     

真空热试验中红外笼对伞状天线非稳态加热过程数值仿真

文章针对真空热试验红外加热装置计算机辅助设计中的非稳态过程模拟问题,以一种较为复杂的伞状天线结构为例,以红外笼作为加热装置,对真空热试验环境下的加热过程进行数值仿真研究。采用红外笼的等效加热面模型简化仿真模型,使计算负荷可以达到工程应用的水平,同时满足模拟精度的要求。数值仿真与试验数据对比证明,所采用的计算方法可以为红外笼等红外加热装置的设计提供准确的非稳态加热过程模拟。     

考虑温度大范围变化的瞬态热固耦合方法研究

针对传统瞬态热固耦合方法不能准确反映温度大范围变化时传热与结构变形之间的耦合效应,基于Galerkin与Newmark算法,建立了充分考虑温度大范围变化对耦合项影响的瞬态热固耦合有限元计算方法。通过求解铝质薄板在强热流作用下的热响应问题,校验了该方法的正确性,并利用该方法对某型高超声速飞行器进气道热防护面板进行了瞬态热固耦合分析。利用理论公式证明并通过分析算例得出结论：温度剧烈变化对热固耦合问题的计算结果影响十分显著;考虑温度剧烈变化的耦合项对结构的温度场影响较小,但会使温度变化率呈现震荡状态,耦合项对结构的位移、变形速度及变形加速度有显著影响,对结构的振动起到阻尼作用。     

基于拟合模态法的大范围运动柔性曲线梁动力学建模方法研究

大范围运动柔性曲线梁的动力学建模主要采用有限元法，这是由于很难得到曲线梁模态的解析模态表达式，求解的复杂程度随着单元的增加而增加。本文提出一种拟合模态法进行动力学建模，把大型有限元分析软件计算得出的模态拟合成为多项式的形式，得到曲线梁模态的近似解析表达式，从而大大简化了建模与计算的过程。以某型号卫星抛物面天线为例，验证了方法的可行性与准确型。