运输机短舱挂架纵向位置优选数值研究

针对运输机的翼/身/架/舱干扰问题,本文基于Euler方程,采用笛卡尔直角网格,应用有限体积法对全机绕流流场进行数值计算,分析和比较了短舱挂架不同前伸量、下沉量的影响,提出了基于最佳前伸量和下沉量的优选配合位置,并与初始位置的全机气动特性进行对比。结果表明,获得的短舱挂架优选配合位置使全机气动特性得到了明显改善和提高。     

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

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

基于统计学检验指标的燃气轮机部件特性方程拟合次数的选择

针对燃气轮机在数学建模过程中难以确定其部件特性方程最优拟合次数的问题,提出了采用统计学检验指标确定特性方程最优拟合次数的方法。最优拟合次数的选择,使得拟合方程既防止由于拟合次数过高,将测量数据中的噪声也纳入拟合模型,又避免了拟合次数太低所引起的结果粗糙。结果表明：综合运用拟合检验指标,避免了由经验选择拟合次数的局限性,并为特性方程最优拟合次数的选取提供了理论依据。     

跨声速颤振计算方法研究

采用改进的弹簧近似方法,以ALE(Arbitrary Lagrangian-Eulerian)描述的Euler方程为控制方程,耦合结构运动方程,计算了Isogai机翼模型的颤振边界。气动力求解采用双时间推进,外时间为物理时间,与结构运动方程同步;结构运动方程求解采用Rayleigh-Ritz方法,时间推进也采用双时间法,节省了大量计算时间。研究表明,计算结果与参考文献提供的结果吻合良好,证明了所提出的跨声速颤振计算方法的准确性。     

八节点Hamiltonian等参元列式

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

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

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

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

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

机构磨损可靠性

首先简述了磨损的常规失效判据,然后提出了磨损可靠性的安全裕量方程,给出了可靠性计算公式。用算例阐明基本概念及方法思路。最后对防止粘着磨损咬卡提出了合理的设计准则。对极限磨损量的代表性作了说明。     

用全位势方程计算机翼的亚声速,跨声速和超声速绕流

对大后掠小展弦比细长机翼,本文对机翼纵轴垂直的每一横流截面生成O型网格,形成对机翼流场的H-O型网格,用守恒型全位势方程、差分和隐式近似因式分解迭代算法计算绕机翼的可压缩位流。自由流可从亚声速直到低超声速的全部跨声速范围。本算法要求机翼前缘有大后掠角,后缘可稍许后掠或前掠。本文算例表明,所研制的计算程序已可提供工程实用。     

二维可压湍流边界层改进的积分方法

采用含高阶项的边界层动量积分方程,同时对常规边界层卷吸方程及延迟方程进行相应的高阶影响修正,得到改进的边界层积分方程组。应用此方程组对多种高亚音速及跨音速翼型边界层流动作了计算,并与一阶Ｇｒｅｅｎ方法计算结果及实验结果作了比较。结果表明,在边界层积分方程中保留法向压力梯度项及雷诺法向应力项,明显改进了翼面边界层接近分离区域处参数的计算精度。