钝体俯仰阻尼导数数值计算

分别采用：（1）Euler方程的谐振摄动法；（2）薄层近似的非定常N－S方程模拟有粘扰动流场两种方法计算了钝体外形的高超声速俯仰阻尼导数。通过计算表明：对于复杂外形，方法（1）适用于零攻角附近的流动计算，而方法（2）考虑了流场的粘性作用，能够应用于一般流态（3＜M＜20，0＜a＜20°）的计算。此外，本文发展了文献［6］的动导数转换公式，使之适用于二维转轴情况，并在方法（2）中通过引入“亚迭代”过程提高了非定常流场的计算精度。     

CONSTITUTIVE MODEL AND DYNAMIC MECHANICAL PROPERTIES OF GLASSCLOTH/EPOXY LAMINATES AT DIFFERENT TEMPERATURES

ＣＯＮＳＴＩＴＵＴＩＶＥＭＯＤＥＬＡＮＤＤＹＮＡＭＩＣＭＥＣＨＡＮＩＣＡＬＰＲＯＰＥＲＴＩＥＳＯＦＧＬＡＳＳ┐ＣＬＯＴＨ／ＥＰＯＸＹＬＡＭＩＮＡＴＥＳＡＴＤＩＦＦＥＲＥＮＴＴＥＭＰＥＲＡＴＵＲＥＳＨＡＮＸｉａｏｐｉｎｇ（韩小平）１,ＨＡＮＳｈｅｎｌｉ...     

矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术,其中非结构网格仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外形;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。混合网格系统由修正的四分树（２Ｄ）／八分树（３Ｄ）方法生成。本文将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲面边界的台阶效应,同时保证了网格间的通量守恒。该算法被成功地运用于ＮＡＣＡ００１２翼型的跨、低超声速绕流,二维弯管内运动激波的非定常流,以及三维捆绑火箭的超声速绕流的数值模拟。计算结果表明本方法在模拟复杂无粘流场方面的灵活性和高效性。     

基于运动嵌套网格的前飞旋翼绕流N-S方程数值计算

 通过求解 Navier-Stokes方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间上采用中心平均的有限体积法进行离散,时间方向应用含隐式子迭代的双时间法推进求解。为了验证程序的正确性,数值计算了一有升力悬停流场,旋翼桨叶表面压力分布的计算值与实验值吻合很好。     

本构关系对连续介质损伤力学基本方程适定性的影响——一维弹性损伤情况

许多工程材料当变形超过某一极限时由于损伤都会出现应变软化行为。首先从非局部理论出发,推导了应变梯度损伤本构方程;然后利用一阶拟线性偏微分方程组的特征理论,在一维弹性损伤情况下分析了两种不同的本构模型,即 Kachanov损伤本构方程与应变梯度损伤本构方程,对连续介质损伤力学基本方程适定性的影响。结果表明,当损伤发展时,与 Kachanov损伤本构模型相关的连续介质损伤力学的基本方程是不适定的;而与应变梯度损伤本构模型相关的,则无论损伤是否发展,其基本方程始终是适定的。这说明在连续介质损伤力学的本构方程中必须考虑材料内部微结构的尺度效应。     

HTPB复合固体推进剂非线性本构研究

采用Ｓｃｈａｐｅｒｙ蠕变型非线性粘弹性本工用Ｍ－Ｃ方法,分析蠕变实验数据,研究了ＨＴＰＢ复合固体推进剂非线性本构关系。结果表明,该方法可以减少加卸载时数据波动的影响,降低参数拟合对实验条件的要求。Ｓｃｈａｐｅｒｙ非线性模型较好地反映了ＨＴＰＢ复合固体推进剂的力学性能,可以用于固体推进剂的力学分析。     

边界层转换阶段直接数值模拟的影响矩阵法

提出一种求解时间相关、三维不可压完全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的影响矩阵法。考虑流动区域为上半空间的平极边界层流动,在水平方向,利用周期性条件,采用Ｆｏｕｒｉｅｒ离散,在半无穷方向,定义了一种特殊映像并利用Ｃｈｅｂｙｓｈｅｖ多项式展开。由此将Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程转化为以水平波数（ｋ１,ｋ２）为参数的一组法向Ｈｅｌｍｈｏｌｔｚ方程组,同时可更有效地求解压力。本文讨论了连续问题的影响矩阵法与原问题的等价性。影响矩阵法的离散形式,以及影响矩阵法的算法实现等。最后给出了一维Ｈｅｌｍｈｏｌｔｙ问题的数值实验结果。     

双时间法求解大迎角翼型绕流非定常N-S方程

本文采用双时间推进法、顶点格式有限体积法求解非定常ＮＳ方程、应用多层网格技术计算了ＮＡＣＡ００１２翼型大迎角振荡绕流,并对伪时间推进收敛精度,物理时间步长,多层网格迭代方式等对计算结果与计算ＣＰＵ时间的影响进行了讨论,另外还讨论了采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ湍流模型时用固定转捩或全湍流模式对计算结果的影响,并与实验结果进行了比较。     

螺旋桨流场Euler/NS数值解对比分析

采用多块点对点网格生成技术,生成螺旋桨流场空间计算网格。针对不同前进比,分别求解Euler和NS方程对螺旋桨流场进行数值模拟,并与实验结果进行了对比分析。研究发现无粘和有粘的计算结果均接近于实验结果,因此基于Euler方程的流场数值模拟更适用于工程设计。     

稀疏非结构网格上的亚声速流高精度数值模拟

采用高阶间断有限元法在非结构网格上数值求解二维亚声速Euler方程。数值结果表明,尽管采用的非结构网格非常稀疏,但通过采用真实曲线物面边界和高阶的基函数仍可得到高精度的数值解。另外,对于超低速情况,方程无需经过任何特殊处理就可以得到收敛的数值解。由于采用牛顿一般最小余量法(Newton-GMRES)时通常需要较好的初始值,本文设计了一种阶谱循环过程来提高数值求解的鲁棒性。