用于棘轮变形预测的棘轮演化统一模型研究

 在大量单轴棘轮实验基础上,研究了均值、幅值、峰值和谷值应力对 304不锈钢的饱和棘轮应变的影响规律,首次提出了棘轮应力σ_r和棘轮门槛值σ_rth 的概念,建立了基于单参数控制的、用于饱和棘轮应变预测的 SRM饱和棘轮本构模型和用于独立循环应力工况下棘轮应变演化预测的 REM棘轮演化模型,并由此发展了全面描述任意循环应力工况下棘轮应变演化规律的 URM棘轮演化统一模型。抛物律模型 SRM和幂律模型 REM对实验数据拟合的相关系数均超过 0.98。URM建模容易,只需 4～ 6个试样的单轴棘轮实验数据。此外还讨论了获得 SRM本构模型的单试样实验法。     

Ni基单晶叶片粘塑性分析

从位错滑移机理出发，建立了一种适用于Ｎｉ基高温单晶结构分析用的晶体滑移粘塑性本构模型。该模型将八面体、十二面体及六面体滑移系作为潜在的开动滑移系，考虑阻应力和背应力两种内应力状态参量。利用某Ｎｉ基单晶叶片材料在７００℃和９５０℃的试验结果对模型进行了考核并标定了模型参数。进而将所建模型编入有限元结构分析软件中，对某型发动机单晶叶片进行了计算分析。     

动静圆柱坐标系下N-S方程的分区算法及子域干涉处理研究

耦合了相对和绝对圆柱坐标系下三维Navier-Stokes方程,统一了动子域和静子域的控制方程;在子域计算中,采用高分辨率的三阶ENN格式及LU-SGS隐式解法;在子域干涉处理中,研究并对比了动静子域干涉面的重合处理与重叠处理;最终发展了一种新的既高效又准确的分区算法。通过对某带进口导叶的三级跨声速轴流压气机叶栅内三维紊流流场的数值模拟,验证了该分区算法的有效性,计算结果与实验数据吻合较好。     

马赫数对侧压式高超音速进气道及等直隔离段三维内流场的影响的数值分析

用 N-S方程和 RNG k-ε紊流模型计算了 RBCC用侧压式高超音速进气道三维内流场 ,重点分析了马赫数对流场的影响。从出口截面上气流参数的均匀程度来看 ,Ma3 时比 Ma6时的更均匀 ;从压缩气流的主要部位来看 ,Ma6情况下在支板前肩点之后到喉部截面之前的这个等宽度收敛通道内 ,而 Ma3 时主要是在前肩点之前的收敛通道内 ;从对气流的压缩程度来看 ,Ma6比 Ma3 更大 ;从唇口激波的强度来看 ,Ma6比 Ma3 更强 ;从进气道的流量捕获率来看 ,Ma6比 Ma3 更高 ;高超音速进气道内靠近侧壁以及侧壁和顶板的角区附近是最有可能出现亚音速的区域。     

超燃冲压发动机燃烧室的准一维计算与分析

在考虑有限速率化学反应的准一维Euler方程基础上，通过增加截面面积变化、壁面摩擦和添质的源项，发展了适用于超燃燃烧室性能分析的准一维计算方法。依次以中国空气动力研究与发展中心（CARDC）和日本国家航空与航天实验室（NAL）的氢燃料燃烧室模型作为验证算例，分别采用传统的一步反应模型和发展的有限速率反应模型，模拟了燃烧室流场，并基于NAL燃烧室，计算分析了不同当量比和进口压强对燃烧室流动特性的影响。结果表明：两种方法都能得到与实验数据吻合良好的结果；和一步反应模型相比，有限速率反应模型不仅可以更细致地捕捉流场细节，而且能够初步分析化学非平衡效应的影响；对于NAL燃烧室，当量比≥0.6时，压强随当量比的升高而增大，当达到1.0时，反压已推进隔离段，且推进速度随当量比增大而增加；进口压强不大于110.444kPa时，反压随进口压强增大而升高，且当反压不小于82.833kPa时，反压被隔离在等直段燃烧室入口处；过小的当量比和过大的进口压强均会导致燃烧室出口马赫数严重下降，甚至出现亚声速出流状态。     

跨声速开槽壁超临界翼型洞壁干扰的N—S方程模拟

本文采用二维可压缩非定常Ｎ－Ｓ方程模拟了带有上下开槽壁的跨声速二维管风洞中超临界翼型绕流的流场。网格生成中内层采用解双曲型偏微分方程，外层采用代数生成法。开槽壁采用只有压力梯度没有顶板运动的Ｃｏｕｅｔｔｅ流动的解析解近似。推导出槽间流动与开闭比的三次方关系。上下壁边界层可由无滑移条件直接模拟，也可用平板边界层近似。本文对驻室压力作了几种近似处理，结果表明，所给Ｎ－Ｓ方程的边界条件与风洞实验多接近一     

解决欧拉方程奇异性的方法探讨

提出了一种解决欧拉方程奇异性的工程算法,从飞行器运动方程入手,推导出具体计算公式,并进行了数值计算验证。该方法无论是理论分析还是工程计算都证明是可行的,没有发现原理缺陷和方法误差,不会产生奇异性,在飞行器仿真、模拟和姿态控制中具有普遍的应用和参考价值。     

虚构解方法程序验证

对求解偏微分方程问题的程序的严格验证,一直以来由于其精确解有限,较难进行,针对这种情况,本文给出了一种新的程序验证方法,虚构解方法,该方法旨在解决大型科学计算程序如何较严格地进行程序验证这一问题;该方法通过构造虚构解,修改原控制方程,然后通过对计算结果的分析处理达到对程序的精确验证;文中给出了一组通用的虚构解的选取办法,解决了虚构解方法验证过程中确定虚构解的问题;同时还给出了运用该方法对非结构二维Euler计算程序进行验证的具体例子,网格收敛分析结果表明该方法是一种有效的程序验证方法,并且具有较好的通用性。     

扑动翼型的低雷诺数气动特性分析

通过求解引入拟压缩项的不可压Navier-Stokes方程,数值模拟了绕扑动翼型的低雷诺数非定常流动。针对厚度在4%-12%之间的NACA对称翼型,分析了翼型厚度等参数对扑动翼型气动特性的影响。在低雷诺数条件下,对于纯俯仰运动,随着翼型厚度的减小,平均阻力系数也变小。而对于纯沉浮运动,发现翼型厚度对气动特性的影响和俯仰运动有很大的差别,平均阻力系数随着翼型厚度的减小而变大。通过对沉浮运动一个周期流线图的分析,认为这是翼型前缘涡的影响造成的。由于前缘涡的影响,翼型厚度增加,平均压差阻力系数变小,甚至会出现负值。雷诺数的影响研究表明,随着雷诺数的增加,扑动翼型的阻力系数减小的趋势越缓慢。     

再入弹抛壳气动特性研究

本文用有限差分方法求解再入母弹头与其所抛壳体之间的无粘干扰流场,获得了体—壳的干扰气动力,并在准定常假设的基础上近似模拟了弹体与壳体的分离过程。数值研究的结果表明,体—壳间存在的干扰会导致壳的压心位置急剧变化,危及母弹安全;抛壳速度是体—壳安全分离的关键参数。