超声速进气道湍流流场的数值模拟

采用有限体积法数值模拟了超声速进气道的湍流流场,改进了离散纳维尔－斯托克斯方程中粘性项的方法,使之具有较好的守恒性,研究使用了变系数的隐式残差光顺等加速收敛技术,并讨论了超音来流、收敛型进气道复杂流场初场的给定方法,所得到的数值解正确地展示了流场的复杂激波结构。     

粘塑性材料结构的有限元分析方法

本文针对Ｗａｌｋｅｒ粘塑性本构方程,通过发展一自适应变步长积分算法,用户子程序的方法把该本构方程结合进了国际著名的非线性有限元软件ＡＢＡＱＵＳ,并进行了算例验证。它包括了四种单元类型,可考虑材料变形的率相关,加载的循环、非等温、变温等情况。适用于发动机热端部件非线性结构分析。     

基于混合动网格的二维非定常粘性流动数值模拟

利用基于Delaunay三角化的动网格方法,在混合网格上求解了带运动边界的二维非定常Navier-Stokes方程。在有限体积法格心格式的基础上,采用带人工耗散项的Jameson中心格式以及双时间步长推进方法对Navier-Stokes方程中的空间项和时间项分别进行了离散。湍流模型为Spalart-Allmaras方程模型。对在跨声速粘性绕流中进行俯仰振动以及后缘摆动的翼型进行了数值模拟,得到了令人满意的结果。     

基于边条优化外形的全机纵向响应与增升效果研究

用遗传基因算法优化设计了翼-身-尾全机模型机翼边条.通过求解全机纵向运动状态方程,分析了全机模型带与不带机翼边条时纵向动态响应问题.涉及的气动力、力距和纵向稳定导数等则由三维低阶板块法计算获得.通过优化设计,新的带有机翼边条的模型飞机,其升力线斜率可明显提高.当Ma=0.4-0.9时,斜率约提高13%-17%;当Ma=1.5时,约提高12%.数值计算结果表明:在装有机翼边条时,由于短周期频率和阻尼比的提高,飞机对在交战飞行时特别的偏转控制输入响应敏捷.最后,给出了部分基于板块法的气动力结果与Cy-20飞机飞行试验数据的比较.     

基于Hill方程的编队初始化误差分析

Hill方程具有形式简单、方便编队队形设计和队形保持控制律的设计等优点.然而,由于Hill方程自身模型误差的存在,导致基于Hill方程初始化所建立的队形具有一定的误差.针对这个问题,基于保留二阶非线性项的、改进形式的Hill方程和考虑偏心率的Lawden方程,分别在非线性和偏心率的影响下,分析了初始化队形的长期漂移和相对运动误差问题,分别给出了针对非线性和偏心率影响的修正初始化相对速度,从而消除队形的长期漂移,修正或减小相对运动误差.最后通过数值仿真进行了验证.     

测力仪校准方法的比较

介绍了中国ＪＪＧ１４４－９２、美国ＡＳＴＭ Ｅ７４－９５、国际标准化组织ＩＳＯ３７６－１９９９检定试验机所用测力仪的主要技术要求,并对此进行了比较。     

相关分析法辨识时域离散数学模型的频率特性

本文提出利用经典过程辨识法中的相关分析法,求取时域离散数学模型的频率特性。通过实例验证了此方法的有效性。此法误差小,能满足研究和设计上的要求。它适用于线性系统,也适用于非线性系统。     

基于厚板的弹性环式挤压油膜阻尼器建模及动力学特性系数识别

利用数值方法完善了弹性环式挤压油膜阻尼器（elastic ring squeeze film damper,ERSFD）的流固耦合计算模型,其中利用雷诺方程建立油膜的控制方程,利用厚板单元建立了弹性环的运动方程并采用分时迭代方法实现了弹性环-油膜的控制方程的耦合求解从而获得瞬时内外油膜的压力,并进一步识别了油膜以及ERSFD的动力学特性系数。结果表明ERSFD的阻尼系数受凸台高度影响最大,凸台高度从0.15mm增加到0.30mm, ERSFD的阻尼系数从5790(N·s)/m减小到718(N·s)/m; ERSFD刚度系数则主要取决于弹性环的厚度以及凸台数目,弹性环厚度从0.8mm增加到1.0mm,ERSFD的刚度从1.44×106N/m增加到2.51×106N/m。      

2024T4铝合金蠕变本构关系研究

设计了一种简易单向拉伸蠕变试验装置,对2024T4铝合金在175℃、多个应力水平下的单向拉伸蠕变进行了试验。根据试验数据的特点,将2024T4铝合金的蠕变本构关系分解为由一个双指数函数和一个一次函数组成的两分量表达式,并利用试验数据通过最小二乘拟合确定出本构关系中的待定参数,得到2024T4铝合金175℃时的本构关系。通过对拟合数据与试验数据的对比分析表明,所建立的蠕变本构关系能够有效表征2024T4铝合金175℃时在不同应力水平下的蠕变行为。     

A Comparison of Techniques for Solving the Poisson Equation in CFD

CFD is a ubiquitous technique central to much of computational simulation such as that required by aircraft design. Solving of the Poisson equation occurs frequently in CFD and there are a number of possible approaches one may leverage. The dynamical core of the MONC atmospheric model is one example of CFD which requires the solving of the Poisson equation to determine pressure terms. Traditionally this aspect of the model has been very time consuming and so it is important to consider how we might reduce the runtime cost. In this paper we survey the different approaches implemented in MONC to perform the pressure solve. Designed to take advantage of large scale, modern, HPC machines, we are concerned with the computation and communication behaviour of the available techniques and in this text we focus on direct FFT and indirect iterative methods. In addition to describing the implementation of these techniques we illustrate on up to 32768 processor cores of a Cray XC30 both the performance and scalability of our approaches. Raw runtime is not the only measure so we also make some comments around the stability and accuracy of solution. The result of this work are a number of techniques, optimised for large scale HPC systems, and an understanding of which is most appropriate in different situations.