八节点Hamiltonian等参元列式

为了使平面八节点等参元的优越性在弹性力学Hamilton正则方程的半解析法得到应用。结合弹性材料修正后的Hellinger—Reissner（H—R）变分原理和二次插值函数表达平面外应力和位移函数,建立了Hamilton正则方程的八节点等参元列式。首先简要地介绍了弹性材料修正后的H—R变分原理．然后用二次插值函数表示平面外应力和位移变量,并详细地推导了Hamilton正则方程的八节点等参元列式。数值实例结果证明了本文等参元列式的正确性。     

跨流域三维复杂绕流问题的气体运动论并行计算

通过研究求解描述跨流域三维绕流问题的Boltzmann模型方程气体运动论耦合迭代数值格式,分析气体运动论数值计算方法的内在并行性;从变量依赖关系、数据通信与并行可扩展性三方面开展基于离散速度空间区域分解计算的研究,发展求解稀薄流到连续流跨流域三维复杂绕流问题的并行算法。通过对不同Knudsen数、不同马赫数、不同攻角三维球体及返回舱绕流的并行计算,计算结果与实验数据和理论分析吻合较好。研究表明,该并行算法负载平衡和并行可扩展性较好,对不同并行计算机系统具有很好适应性,显示直接求解分子速度分布函数的气体运动论计算方法有良好的并行计算特性。     

径向基函数在三维Euler方程数值计算中的应用

基于三维Euler方程,对非结构四面体网格给出了一种基于紧支径向基函数重构的ENO型有限体积格式,方法的主要思想是先对每一个四面体单元构造插值径向基函数,而在计算交界面的流通量采用高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法。最后用该格式对一些典型算例进行了数值模拟,结果表明该方法计算速度快,对间断有很好的分辨能力。     

含弱粘结界面复合材料层合板的半解析法

采用线性弹簧模型表征界面弱粘结特性．并运用了Hamilton正则方程的半解析方法处理界面弱粘结对复合材料层合板弯曲问题的影响。首先简要地介绍了弹性材料修正后的H—R变分原理和Hamilton元,然后着重推导了弱粘结界面的线弹簧模型的平面元素列式及正则方程的半解析法的控制方程。最后根据界面传递关系建立了弱粘结层合板的控制方程。数值实例结果证明了本文方法的正确性。     

同轴双剪切气-气喷嘴试验研究

为研究应用于全流量补燃循环发动机的新型结构气-气喷注器,对以气氢和气氧为推进剂的同轴双剪切喷嘴进行试验研究,分析喷嘴不同中心氢流量比例和流量变化对推进剂燃烧效率的影响。结果表明:中心氢流量与外圈氢流量的比例是同轴双剪切喷嘴重要的设计参数。当中心氢流量比例为0.3至0.4之间时,同轴双剪切喷嘴使推进剂高效燃烧;在喷嘴流量相当于航天飞机主发动机单喷嘴8倍的工况下,同轴双剪切喷嘴能使推进剂达到高的燃烧效率。     

水平非均匀大气条件下的抛物方程模型研究

在分析传统分步傅立叶(SSFT)和中心差分离散混合傅立叶变换(DMFT)算法的基础上,针对算法复杂度及数值稳定性问题,提出了前向-后向混合差分DMFT算法及其具体实现步骤,提高计算效率的同时解决了传统算法在复杂条件下的盲点问题.给出了算法中水平非均匀大气折射剖面及复杂条件下最大传播仰角θmax等关键参数的计算方法,并与传统模型及经典文献实验数据进行了仿真比较,验证了该算法及模型在复杂大气条件下的正确性和有效性.     

空间拦截防御的射击诸元快速解算算法

针对空间平台发射非制导毁伤单元对近程高速非合作目标进行拦截防御的射击诸元快速解算问题,建立了命中方程,并基于命中方程的一阶解,给出了2种解析算法确定发射时刻、发射速度矢量与提前点坐标之间的函数关系.提出了计算命中时刻、发射速度矢量与发射时刻间函数关系的二分数值方法.最后通过仿真分析了算法的方法误差及适用范围.     

两种物态方程在超高速碰撞数值模拟中的比较

材料在超高速碰撞下将发生熔化甚至气化相变。为反映超高速碰撞中材料相变带来的影响,文章采用GRAY三相物态方程与Tillotson物态方程,对超高速碰撞进行数值模拟对比研究。研究表明:当碰撞速度在3 km/s以下时,除了靶板穿孔直径外,两种物态方程所给出的碎片云的结果基本一致;但当碰撞速度在3 km/s以上时,两种物态方程给出的数值模拟结果有较大差异,这说明在超高速碰撞中相变的产生对碎片云形状参数有较大影响。     

四边简支压电热弹性体层合开口壳的精确解

根据压电热弹性体的本构关系和热传导关系的相似性,通过增维方式将热传导关系写进压电热弹性体的本构关系中,直接得到了一个新的增维的压电热弹性体本构关系。依据此增维的本构关系,结合压电热弹性体平衡方程和热平衡方程,成功地导出了压电热弹性体层合开口壳机、电、热耦合问题的齐次状态方程。这类方法大大简化了人们在分析压电热弹性体耦合问题时通常要求解的非齐次微分方程和关于温度的平衡方程和导热方程的二阶微分方程的繁琐方法,同时,齐次状态方程的提出大大减少了数值计算的工作量。     

桥联重积分方程一种新的数值方法

纤维增强复合材料断裂力学中的桥联问题常常可归结为一个二重积分方程。常见的求解方法有矩阵迭代法和积分迭代法，但这两种方法都产复杂。本文给出一种新的分段插值方法，对此积分方程进行求解。数值计算结果表明，该方法效果令人满意。