激波与物质界面相互作用的数值模拟

计算激波与物质界面相互作用问题的关键是精确捕获运动物质界面.采用level set捕获物质界面,用改进虚拟流(modified ghost fluid method)方法定义虚拟流节点参数,从而将多介质流问题分解成多个单一介质流问题求解.应用该方法对激波与物质界面相互作用系列问题进行了数值模拟,并给出了物质界面的演化过程.计算结果发现激波强度、物质界面对涡的形式有很大影响,界面两边的密度比决定了界面是否反转.     

一种新的用于MHD方程逆风格式的Jacobian矩阵分裂方法

针对理想MHD方程,提出了一种新的基于MacCormack算法的雅可比矩阵分裂方法,克服了原有方法稳定性差的问题,并成功地应用于理想MHD方程的求解.控制方程在非结构混合网格上进行空间离散,其中对流项采用本文发展的逆风向量分裂格式,并引入了双曲型磁场散度清除技术,时间推进为显式5步龙格-库塔方法.对MHD激波管流动和带均匀磁场干扰的二维高超声速钝头体绕流流场进行了数值模拟,得到了与参考文献相吻合的数值结果,表明本文发展的数值分裂方法可以有效地捕捉MHD流场的流动特征,并且具有比MacCormack方法更高的稳定性和计算精度.     

航天器空间碎片超高速撞击防护的若干问题

空间碎片对航天器的超高速撞击损伤已受到国内外的普遍重视,如何使在轨航天器对空间碎片进行有效防护是航天器长寿命、高可靠安全运行的重要保障。文章概述了空间碎片环境现状、空间碎片超高速撞击危害以及国内外空间碎片防护的研究现状和趋势。重点介绍了航天器常用的Whipple防护结构及其各种衍生结构的防护性能和弹道极限方程（BLE）,评述了这些防护结构防护性能的优缺点。     

基于Kriging模型的高超声速舵面优化设计

采用RANS方程进行了高超声速舵面的气动特性计算与平面形状优化设计研究。为了提高优化效率,采用拉丁超立方抽样试验设计方法得到样本点,并建立Kriging模型替代费时的流动数值模拟。以最大化舵面操纵力矩和最小化平面面积为目标,以阻力不增为约束条件,进行了舵面平面形状优化设计。优化结果表明,优化舵面的平面面积减少了4.78%,力矩绝对值在设计点增加了24.69%,且在较大的攻角范围均有提高。将初始舵面和优化舵面分别装配到方形截面导弹上,验证计算表明,优化舵面的操纵效率更高。     

盒式机翼布局气动特性研究

采用基于S-A一方程湍流模型的雷诺平均法求解N-S方程,研究了盒式机翼布局前后翼几何参数对气动特性的影响规律,为盒式布局飞机气动设计提供参考。结果显示:随着前后翼垂直距离的增大,盒翼升力和升阻比均增大;增大前后翼的纵向距离,盒翼升力增加;盒翼气动性能对两翼纵向距离的敏感程度不及对两翼的垂直距离。     

鸟撞平板试验与鸟体本构参数识别方法

准确可靠的鸟体本构模型与参数是开展结构抗鸟撞分析和设计的基础。为获取鸟体本构参数,在139m/s的速度下,开展鸟体撞击铝合金薄板试验,并测得了撞击载荷、平板变形和应变数据;结合ISIGHT和PAM-CRASH软件,以薄板被撞击位置的法向变形为优化目标,提出一种鸟体本构参数识别方法。将识别得到的参数代入数值模型并开展数值分析,提取平板变形结果和撞击载荷数据,并与试验测试数据进行比较。结果表明,平板变形仿真结果与试验数据具有很好的一致性,撞击载荷分析结果与试验数据基本一致,证明了本文所建立的鸟体本构参数识别方法的有效性。     

基于筏化的镍基单晶合金高温蠕变模型

采用应力法以及界面能法对不同取向的镍基单晶合金在950℃下的筏化类型进行了预测.上述两种方法的筏化预测结果一致,[001]取向为N型筏化,[011]取向为P型筏化,[111]取向不筏化.进一步,根据筏化预测结果以及晶体滑移理论,结合Kachanov-Robotnov（K-R）损伤演化公式,建立了一个镍基单晶合金蠕变模型,采用该模型并结合商用有限元软件Abaqus的用户材料子程序（UMAT）二次开发接口,对[001],[011]和[111]取向下的CMSX-4镍基单晶合金,在950℃,180~450MPa应力条件下的蠕变变形行为进行了模拟.该模型能够准确预测镍基单晶合金的筏化类型以及滑移系开动规律,更加符合材料的蠕变变形物理机制,因此模型可以对镍基单晶合金的高温蠕变曲线的第2,3阶段进行很好的模拟,并得到了试验的验证.      

双脉冲发动机EPDM软隔层黏超弹本构模型

三元乙丙橡胶（EPDM）用于软隔层式双脉冲固体火箭发动机软质脉冲隔离装置（PSD）时,主要在Ⅰ脉冲工作时起到绝热抗烧蚀作用,同时保证Ⅱ脉冲能够可靠工作。根据橡胶类材料连续介质力学理论,建立了描述EPDM软隔层在有限变形下的黏超弹本构模型。模型由超弹部分和非线性黏弹性部分构成：超弹部分在Mooney-Rivlin模型的基础上进行了改进,使之能够描述大应变时的硬化现象;非线性黏弹性部分采用广义黏弹性模型,采用无量纲形式的KWW方程替代了传统的Prony级数,使得用2个参数就能预测较大应变率范围内的力学响应。利用万能材料实验机对EPDM软隔层进行了多步松弛实验和单轴等速率拉伸实验,然后根据实验结果,采用分步拟合的方法求出模型参数,利用所建立的本构模型对其余的实验结果进行预测并与实验结果进行比较,对比结果表明所建立的模型能较好地预测EPDM软隔层伸长比在800%以内的单轴等速拉伸响应。最后利用文献中的实验数据验证了所建立的模型能够较好地预测多种工况下的力学响应。     

形状记忆合金薄板低速冲击载荷下热力耦合行为分析

进行形状记忆合金(shape memory alloy,SMA)薄板不同冲击能量条件下低速冲击实验,以SMA 热力学本构模型为理论基础,通过数值手段对超弹性SMA 薄板的冲击响应进行仿真分析,探究冲击载荷作用下SMA 热力耦合行为特征。结果表明:数值仿真结果与实验数据吻合良好,有效表征了SMA 薄板冲击过程中的变形、相变、耗散、温变等热力耦合行为特征。     

基于非结构网格求解器的CHN-T1标模气动特性计算研究

高可信度计算流体力学(CFD)技术在大型运输机精细化气动设计中发挥着越来越重要的作用。虽然目前基于RANS方程的数值模拟技术已经应用于大型运输机复杂外形的气动性能计算与分析中,但其可信度仍需进一步确认。针对AeCW-1组委会发布的CHN-T1运输机标模,采用基于RANS方程的非结构网格求解器TNS进行了气动性能计算,并与组委会提供的风洞试验结果进行了对比。此外,还研究了网格量和湍流模型等因素对巡航气动特性计算结果的影响,以及支撑干扰和静气弹变形对跨声速抖振特性的影响。结果表明,TNS预测的巡航升阻比及跨声速抖振特性等气动性能与试验值吻合良好,计算结果完全位于试验值误差带内,验证了所采用的非结构网格求解器在大型运输机复杂外形气动力预测中的可靠性。