大中型民用飞机起飞总重估算方法研究

推导了基于发动机推力估算起飞总重和基于旅客座位数估算起飞总重的公式。通过多种大中型喷气式民用客机相关参数的统计值,建立起飞总重对相关参数的多元回归分析模型,通过显著性检验,说明模型是有效的,并据此建立了新的起飞总重估算公式。通过算例结果分析表明,新的起飞总重估算公式具有很高精度。在新机的方案设计阶段若能给定此公式所需的相关参数,可以计算得到相对合理的设计起飞总重值。     

适用于气动学问题的DRP类有限差分方法

研究了一类低耗散低色散的高阶精度有限差分方法,目的是直接计算非定常欧拉方程用于气动声学问题.采用的数值方法是空间四阶、时间三阶精度的色散关系保持(DRP:dispersion-relation-preserving)类有限差分格式,通过波动方程算例验证了格式模拟波动问题的能力;采用加变幅值的高波数人工耗散项来抑制高频数值振荡,使得该格式可推广用于含激波、初始间断或非线性波动的问题.通过含激波、初始间断条件的算例验证了所用变幅值人工耗散的有效性,最后作为初步应用计算了二维亚、超音速均匀流中简单点声源辐射问题,得到了很好的结果,表明此类方法在气动声学问题计算上深有潜力.     

γ-Re_θ模式应用于高速边界层转捩的研究

应用γ-Reθ转捩模式对超声速、高超声速边界层转捩进行数值模拟。γ-Reθ模式通过求解关于当地雷诺数和间歇因子两个输运方程给出转捩起始位置和转捩区长度等信息。本文对马赫数3.5至7范围内的四种算例进行计算,研究了来流雷诺数、攻角变化和头部钝化半径等关键参数的变化对γ-Reθ模式预测流动转捩性能的影响。给出了壁面摩擦阻力系数、热流值与实验值的对比以及壁面附近间歇因子等值线分布等计算结果。γ-Reθ模式能够正确预测攻角和钝化半径变化时转捩位置和流动参数的变化趋势;在较大雷诺数时计算结果与实验值吻合很好。γ-Reθ模式对于超声速或高超声速边界层转捩的模拟仍需修正和改进。     

湍流模型及施密特数对点火位置的影响研究

准确预测燃料和空气湍流混合水平对精确模拟发动机性能至关重要.通过对Burrows&Kurkov氢气顺喷经典算例的数值模拟,研究了湍流模型和湍流关键参数(施密特数)对氢气和空气的混合程度以及点火位置的影响.数值结果表明:湍流模型对点火的位置具有一定的影响,其中Menter SST k-w湍流模型计算结果与实验值总体吻合最好;湍流施密特数严重影响着氢气和空气的混合程度以及燃烧流场的点火位置,其中Sct=0.5能满足大部分超燃冲压发动机燃烧算例的数值模拟.     

大坝渗流反问题初步研究

对大坝渗流反问题进行了初步研究,基于渗流场的有限元计算建立了反演坝体渗流参数的计算方法.该方法将参数反演问题转化为参数优化问题,利用大坝运行水位变化过程中坝体内某些测点上的水头测量数据,用牛顿-拉夫逊算法优化出坝体的渗流参数.通过算例计算分析表明,该反演方法具有推导简洁、反演精度高、适应性好的优点,有较好的工程实用价值.     

基于火焰面和射流模型的双模态冲压燃烧室一维性能计算

为了快速地得到双模态燃烧室中的参数分布,采用射流模型,获得燃料混合分数的时均值、脉动值以及当量标量耗散率的时均值在双模态冲压发动机燃烧室中的三维分布;采用火焰面模型得到燃烧组分和释热分布,将三维信息通过流道横截面面积平均加入到一维非稳态计算模型,求解流动控制方程得到双模态燃烧室一维参数分布。通过不同类型算例验证,表明该方法能较准确地预测燃烧内壁面压力的变化趋势,克服了燃料射流非一维因素对一维计算方法带来的不利影响。     

雷诺应力模型在三维湍流流场计算中的应用

从雷诺应力模型出发,通过求解雷诺平均N-S方程组获得三维湍流流场的数值解。计算中比较了多种湍流模式,并进行了相应的流场计算。本文完成了两个典型算例。从算例1的计算结果与实验值比较中发现:采用雷诺应力模型(RSM)计算的三维流场比采用k-ε模型更贴近实验值;算例2采用了RSM模型及三维非结构网格,对一典型内流问题进行了三维流场计算。算例的数值实践表明:采用雷诺应力模型可以有效的计算各向异性的湍流流场;另外,发展非结构网格有助于模拟壁面附近的流动,并节省计算机内存。      

基于SMP机群的航空CFD并行技术应用研究

针对SMP机群的双层并行体系结构,研制出一种双层并行技术--在计算节点间采用MPI消息传递模式并行,而在计算节点内部采用OpenMP编译器指令模式并行.通过算例测试证明,采用该技术的航空CFD程序,并结合贪婪负载平衡分配技术进行并行计算,其效率可达94%,计算结果与实验值吻合.     

小子样场合下估算母体百分位值置信下限和可靠度置信下限的Bootstrap方法

将Bootstrap方法引入到小子样场合下母体百分位值置信下限的计算中,并与传统的单侧容限系数法和新单侧容限系数法进行了大量的对比计算,算例表明Bootstrap方法明显优于其他两种方法,在较高可靠度和置信度要求下母体百分位值置信下限不会出现负值的情况,而且计算得到的结果更接近真值。发展了一种半参数Bootstrap方法用于计算可靠度的置信下限,模拟计算表明半参数Bootstrap方法很好地克服了置信度较高时,新旧单侧容限系数计算出的可靠度置信下限过低的局限性,而且半参数Bootstrap方法的计算结果均有较高的精度。     

固体火箭发动机稳态燃速二维模型参数最优化辨识

1引言发动机参数辨识常用方法是最小二乘法[1],即将参数估计问题转化为多元函数求极值问题,通过迭代选择待估参数值,使燃烧室压强计算值和实验值的残差平方和最小。采用最小二乘法进行参数辨识存在局限性[2],需要通过台劳级数展开式近似模型替代真实数学模型计算残差平方和对待     

基于特征正交分解法的翼型结冰冰形快速预测

为加快翼型结冰冰形的计算速度,提出了一套基于特征正交分解(POD)法的结冰冰形快速预测算法.通过计算流体动力学(CFD)数值模拟计算得到的冰形结果作为样本,以结冰温度单参数变化时为例,详细介绍了POD法预测结冰冰形的实现步骤.考虑结冰温度、结冰时间以及液态水含量的影响,完成了单参数、两参数与三参数变化时的POD结冰冰形快速预测.通过算例结果发现:完成POD预测计算只需几秒钟且POD法与CFD法得到的冰形吻合较好,仅在冰角附近有较小差别,表明POD方法能够快速、准确地得到翼型的结冰冰形.      

弧齿锥齿轮误差敏感性优化设计

从齿面结构上提出了通过优化差曲面全曲率来改善弧齿锥齿轮的安装误差敏感性问题的方法.推导了啮合点处沿齿线方向的两啮合齿面全曲率作为敏感性系数,分析了局部综合参数和参考点位置参数对参考点处的敏感系数的影响.提出了通过优化传动比一阶导数、接触迹线方向、二阶变性系数和三阶变性系数,获得对安装误差敏感性低的小轮加工参数.算例表明:优化后的齿轮副在啮合过程中的敏感性系数控制在较小范围之内,传动误差的幅值和对称性均满足设计要求,改善了齿轮副的啮合质量.      

螺旋桨滑流-机翼干扰效应的数值计算方法

本文用有限基本解数值计算方法来计算螺旋桨滑流对机翼气动特性的影响,并编制了FORTRAN语言程序。对于带有一个、两个或四个螺旋桨的升力面算例,计算结果与实验值很接近,比国外资料中的解析理论结果大有改善。     

一种设计参数型的有限元模型修正方法

本文提出了一种设计参数型的有限元模型修正方法。依据正交条件,采用系统矩阵对设计参数的台劳展开式,建立参数修改量的线性方程式,用奇异值分解方法求解。在本方法中充分考虑了试验模态参数误差的影响,及试验振型自由度扩充等问题。计算机仿真及实际算例的结果表明,本方法计算量小、精确度较好。     

燃烧室中辐射热流分布的蒙特卡罗计算

本文用蒙特卡罗法计算了燃烧室壁面辐射热流密度分布。计算中考虑了燃气参数如:压力、温度、燃气组分及燃气的吸收和发射特性沿燃烧室轴向、径向和周向的变化。详细地说明了计算过程并进行了算例计算。计算所得热流密度值与分布规律是合理的。      

螺旋锥齿轮振动研齿系统的拍击响应

建立了位移激励下的螺旋锥齿轮振动研齿系统的扭转振动模型,推导出考虑摩擦的轮齿斜碰撞前后状态间的局部映射关系式,得到了螺旋锥齿轮振动研齿系统的Lyapunov指数计算方法.通过算例,运用相图、Poincaré映射图、Lyapunov指数和系统分叉图,分析了研齿系统随激励参数变化的动态响应,结果表明,随激励幅值的增大,研齿系统由单边碰撞状态经倍周期分叉进入混沌,并在混沌状态中有周期窗口的出现.继续增大激励幅值,系统进入了双边碰撞状态.      

基于UKF的实时弹道解算参数稳健设计

在基于UKF(Unscented Kalman Filter,无迹卡尔曼滤波)的实时弹道解算应用中,外部环境不稳定的影响以及观测过程中野值的存在,使得基于经验参数设置的滤波易于发散.因此,需要针对指定型号调节和设置滤波器参数的取值,使其克服上述现象.针对这一问题,提出了基于理论弹道先验的参数稳健设计方法,对滤波器参数最优值进行预测.从滤波器的架构出发,梳理出影响滤波器性能的相关参数,结合稳健设计思想,外设计部分基于理论弹道先验信息,通过模拟弹道特征点构造弹道观测模型;而内设计部份,构造模型的测元数据进行正交试验,判断参数的显著性,获取相应的最优参数取值,从而为实时解算弹道的参数设置提供指导.实际数据试验结果表明,该方法可显著提高滤波器的实算稳健性和收敛性能.     

基于FFD的气动外形参数化方法

主要介绍了FFD气动外形参数化方法的相关理论,给出二维算例。针对FFD方法的不足,又提出了隐式直接操作FFD方法和显式直接操作FFD方法,并计算了三维算例。通过对比,得出各种参数化方法的优缺点。对气动外形优化设计过程中的变形问题有一定的参考价值。     

基于均值迭代法的稳态数据处理方法

参数不确定度的分析及合成是动力装置性能试飞中最为关键的环节之一,而参数名义值的获取以及剔除异常数据是不确定度分析的基础。本文介绍了一种基于均值迭代的稳态数据处理算法,详细分析了该方法的工作原理、名义值求取步骤。算例计算表明,本方法对数据样本的概率分布无要求,适用性强,避免了一些传统算法的复杂计算,而且可以实时地评估稳定段数据的抖动程度,适合于参数样本分布未知的航空发动机关键截面参数的数据预处理。     

高升阻比翼型和机翼的多点优化设计

研究并给出了一种适用于CFD计算的尾迹面气动力计算方法,并讨论了实现尾迹面法的气动力积分技术.数值算例的结果与传统的表面积分法的结果比较表明,尾迹面法的结果更接近于实验数据.结果还表明此计算方法利用CFD可给出工程允许的精度范围的阻力值,因此也可用于复杂外形的阻力优化计算.讨论了多目标优化及综合目标函数构造的方法;给出了基于N-S方程流场求解的高气动性能翼型的优化设计方法和算例;讨论了基于欧拉和N-S方程流场求解的三维机翼优化设计方法和算例.翼型和机翼的优化设计算例结果表明,优化设计计算方法可有效地提高翼型的升力系数和降低阻力系数,具有更高的综合气动性能;三维机翼厚度非线性分布和负扭角会改善机翼流场的流动状态、提高机翼的升阻比.算例结果也表明所提出的优化设计方法是正确和有效的.