夹层板精密成形的数值模拟

为了实现基于“点阵钉模、真空负压、蜂窝夹层”的双曲度面板精密成形过程的有限元数值模拟,分析了成形 2个阶段的变形过程,建立了简化的力学模形,以一种广义协调薄板弯曲四边形位移单元为基础构造了具有公共横向位移自由度的几何非线性四边形四节点分别考虑和不考虑横向剪切的多层板单元,通过牛顿 -拉斐逊法结合载荷、位移约束量的增量法求解几何非线性方程,较好地再现了真实成形过程     

飞行器水载荷结构完整性数值模拟现状与展望-Part I:水上迫降和水上漂浮

现代飞行器面临水上迫降、水上漂浮、贮箱晃动和投汲水等复杂水载荷的结构完整性和乘员安全性分析问题日趋重要,随着科学技术的发展,数值模拟已经成为飞行器设计、分析和适航取证的重要手段。以固定翼飞机、水陆两栖飞机、直升机、火箭和卫星等现代航空航天飞行器为对象,围绕适用于飞行器水载荷分析的数值模拟方法进行综述,根据外流（水上迫降和水上漂浮）和内流（贮箱晃动和投汲水）的不同将综述内容分为Part I和Part II两部分。Part I的主要工作为：首先,归纳水上迫降和水上漂浮的事故和试验,总结水气两相流和流固耦合算法的发展现状和优缺点;随后,结合工程实际,介绍飞行器水上迫降和水上漂浮的范畴、水载荷分析要点、适用的数值模拟方法和软件的国内外发展情况,其中,水上迫降的总结包括飞行参数、波浪水面和弹性体对迫降性能的影响研究,水上漂浮的总结涵盖了飞行器构型参数、破舱和波浪对漂浮性能的影响研究;最后,指出复杂风浪情况下水上迫降和漂浮的水气固三相耦合工程应用难点和解决途径,并探讨飞行器水载荷数值分析的技术挑战和未来的发展方向。     

防冰热载荷计算的一种新方法

提出了一种计算飞机防冰热载荷的新方法,该方法利用CFD软件FLUENT及其自带的用户自定义函数UDF进行数值模拟。首先用FLUENT的两相流欧拉模型模拟空气-水滴两相流动,进而使用FLUENT内部的数据结构,通过UDF二次开发计算水滴撞击量和防冰热载荷,并且利用FLUENT的后处理功能对计算结果进行可视化图形显示。以某3D迎角传感器为例,进行了三维空气-水滴两相流和防冰热载荷计算,防冰所需热功率计算结果和实验数据误差为7.5%,在工程计算允许范围内。     

基于数值风洞技术的民用飞机系留气动载荷计算研究

基于数值风洞技术,针对在低速风洞进行的民用运输飞机的地面系留载荷风洞试验进行了模拟研究。使用CFD软件建立了风洞试验段和飞机风洞模型的数值模型。根据风洞试验条件设置数值计算条件,计算了侧滑角在0毅~ -180毅范围内的飞机模型气动力。由于数值模型包括了洞壁、地板,计算结果和试验结果吻合良好。结果分析表明:当-120毅臆茁臆-90毅时,前起落架系留装置将承受较大的载荷,该载荷主要源于大偏航角时平尾部件产生的抬头力矩引起的纵向载荷,以及偏航力矩引起的横向载荷;Cl 比Cm 小两个数量级,对系留载荷的影响不占主导地位;将平尾、升降舵或方向舵预偏可降低系留载荷。为民用飞机低速风洞试验数值建模提供了参考。     

控制体有限元方法中的一种高阶插值格式

针对控制体有限元方法对流动方向不平行于坐标轴方向的回流问题数值模拟采用各种迎风格式都容易产生假扩散的问题,本文给出了一种由插值单元及其迎风相邻单元所确定的高阶插值格式。通过对Re=1000时三种不同形状的二维四边形空穴流模拟,表明该格式与局部斜迎风格式的适当组合可有效地减小控制体有限元方法对回流问题数值模拟的假扩散。      

铝锂合金加筋壁板剪切屈曲性能

为研究剪切载荷下2A97铝锂合金加筋壁板的屈曲与后屈曲行为,设计了加筋壁板和夹具,完成了壁板的剪切试验;得到了加筋壁板的失稳载荷、破坏载荷以及破坏模式;采用受剪板屈曲与张力场理论计算了加筋壁板的剪切屈曲失稳载荷;建立有限元数值计算模型对加筋壁板屈曲行为进行计算分析,并将数值结果与试验结果对比。结果表明：加筋壁板的屈曲模式为筋条间蒙皮的局部屈曲;加筋壁板的破坏模式为沿加载对角线方向蒙皮的凸起,破坏原因为蒙皮的塑性变形、撕裂以及筋条的扭转变形;利用张力场理论可以得到较准确的屈曲失稳载荷,与试验误差为6.56%;数值模拟得到的屈曲与破坏模式与试验吻合,失稳载荷和极限载荷与试验结果误差分别为1.22%和11.52%。     

基于激励盘理论的螺旋桨滑流数值模拟研究

针对螺旋桨滑流流动复杂、非定常数值模拟难度大的问题,对激励盘模型进行深入研究,提高其对滑流模拟的精度。首先根据非定常数值模拟单独螺旋桨流动结果确定激励盘载荷。分析了考虑旋向载荷分布对滑流数值模拟的重要性。采用激励盘模型得到的定常滑流速度分布与螺旋桨非定常滑流的时间平均结果非常吻合,且采用激励盘理论能够明显地缩短模拟时间,在工程应用中具有很大的优势。     

转、静子轴向间距对转子叶片振动应力的影响

采用数值模拟方法,分析了转、静子轴向间距对转子叶片振动应力的影响.利用CFD程序对进口静叶尾迹作用下的压气机转子叶片排非定常流动进行了数值模拟,并采用所编写的程序将气动载荷转化为有限元模型上的动态压力载荷,进行了结构瞬态响应分析.结果表明:在所考查的轴向间距下,随轴向间距的增大,转子叶片振动应力呈减小趋势,这主要是由于造成压力载荷波动的尾迹强度沿流向是减弱的.      

亚声速谐振升力面的点偶极子法

本文介绍计算亚声速谐振薄翼非定常气动载荷的点偶极子法(DPM)。网格划分方式与偶极子网格法(DLM)相类似。此法将网格上的载荷分布集中置于网格中央剖面的1/4弦点上,控制点取在3/4弦点上,无需作数值积分计算。很容易用于计算形状复杂机翼的非定常气动载荷。 文中给出了矩形翼和后掠翼定常和非定常气动载荷的数值结果,其精度数值与DLM的相当,计算效率提高一倍。     

固体火箭发动机内气粒两相流动的SPH-FVM耦合方法数值模拟

为解决现有数值模型在求解固体火箭发动机内颗粒的燃烧流动等问题时所遇到的难题,基于拟流体模型,采用光滑粒子流体动力学方法 (SPH)求解离散颗粒相,追踪颗粒运动轨迹,采用有限体积方法 (FVM)对气体连续相进行描述,捕捉流场特性,两相间通过曳力、压力梯度、热传导、体积分数等参量进行耦合,建立了两种不同坐标系下方法间的耦合框架。耦合方法模拟了喷气推进实验室(JPL)喷管中气粒两相流动过程,得到了气相和颗粒相的参数分布规律,与相关实验结果及离散颗粒模型(DPM)数值模拟结果对比,吻合较好。针对Φ315mm实验发动机工作状态和结构特点,进行了发动机燃烧室内过载条件下的气粒两相流动数值模拟,分析了纵、横加速度载荷对发动机燃烧室内粒子场和聚集带的影响,与已有的数值模拟结果基本吻合。结果表明,新方法求解发动机气粒两相流问题准确可靠,适用于发动机内更为复杂的两相流问题数值模拟。     

超声速前机身及翼－身组合体的全速势方程数值计算

通过数值求解全速势方程，计算了超声速来流的前机身及翼－身组合体。当流场存在亚声速区时，在此区域内采用中心差分格式，迭代求解，并引入多重网格技术，加快收敛；当流场中某一区域沿某一方向是超声速时，在此区域内采用沿该方向的推进解法。计算结果表明，本文的方法可靠，结果准确，可以向工程应用方面推广。     

亚超音速机翼最佳弯扭综合设计计算方法

本文介绍一种亚超音速机翼最佳弯扭综合设计的计算方法,它应用了有限基本解方法。分别在亚超音速各选取一个设计点(M数和C_L),进行机翼弯扭设计,其目的是减小与升力相关的阻力。在此基础上,顾及亚音速和超音速这两个设计点的气动力特性,还要兼顾到飞机其它性能和结构上实现的可能性,进行机翼的综合设计。本文分别给出了亚音速最佳弯扭设计,超音速最佳弯扭设计和综合设计的计算结果。经过分析表明,计算结果是合理的。     

二维带动力吸气式高超声速飞行器绕流的PNS-NS混合求解

为了提高吸气式高超声速飞行器绕流的求解效率,采用空间推进求解抛物化PNS(parabolizedNavier-Stokes)方程和时间迭代求解Navier-Stokes(N-S)方程的混合计算流体动力学(CFD)方法来求解高超声速飞行器整机绕流.在超声速占主导的流动区域采用空间推进求解抛物化N-S方程的方法,在亚声速和分离区采用时间迭代求解N-S方程的方法.对于求解二维带化学反应的吸气式高超声速飞行器绕流,混合CFD方法和完全时间迭代方法相比,可得到同等准确的数值模拟结果,并且求解效率提高了数倍.      

叶轮机平面叶栅的气动设计及统流计算

通过全速势方程的数值积分，利用ＡＦ２格式，给出了一种用于叶轮机中叶栅型面的气动设计和绕流计算的便捷方法，对于稳态亚音速、跨音速及超音速流动均适用。与目前普及很广的、在超音区应用人工压缩性的方案相比，该算法有明显的优势。     

跨音速叶栅流场的无粘-边界层迭代计算

二维跨音速叶栅流场的计算是在无粘时间相关法计算程序的基础上采用无粘-边界层迭代而得。马赫数从亚音速到超音速范围的几种数值试验均与测试结果吻合一致。在边界层无强烈脱流时可获得很好的结果。     

旋成体轴对称跨声速全位势流的高效差分算法

本文用守恒型全位势方程,贴体坐标网格,对旋成体轴对称跨声速绕流的差分数值计算方法进行了研究;根据最佳收敛准则,提出了轴对称情形的AF2迭代算法,并将此算法应用于半球头柱体、弹体等各种外形的旋成体。与一般方法仅适用于亚声速自由流不同,本计算可从亚声速、跨声速自由流一直到低超声速自由流。计算结果表明,本文方法收敛快,与实验及其它方法的结果符合较好。     

再人体复杂外形无粘流场数值计算研究

本文在文[1,2,3]的基础上研究了非定常Euler方程的推进迭代方法,并将无波动、无自由参数的耗散差分格式(下称NND格式)发展为隐式、迎风格式,用以计算有差曲控制翼的再入复杂飞行器的育攻角、有侧滑角的超声速无粘绕流流场。其结果是准确的。本文的方法具有以下特点: 1.可以在同一计算程序中实现亚跨超流场的计算:在超声速区域,采用推进技术;在亚跨声速区,则自动选择推进迭代求解技术。 2.NND格式不仅具有公式简单,无自由参数的特点,而且过激波时满足熵条件,不会产生非物理解,也不发生振荡。 3.隐式差分格式的无条件稳定和Gauss-Seidel迭代的快速收敛,计算省时。 本文的方法曾用于烧蚀头部凹陷外形和航天飞机简化外形的计算,能自动选择迭代和推进区域。可以预言,如计算机条件许可,本方法可以用于大攻角的无粘流场计算。     

飞机亚跨超绕流的数值模拟研究

本文以Euler方程为数学模型,采用一种高精度的TVD(Total Variation Dimishing)离散格式及一种含近似因式分解的推进迭代方法,求解亚跨超绕流’流场。通过若干算例的试算,证明方法是可行的,可以用来模拟飞机的复杂流场。     

亚跨超音速喷管内流场的整体矢通量分裂算法

从N－S方程出发，采用先进的矢通量分裂算法，对亚跨超音速喷管内流场进行了整体模拟计算。文中对几种不同计算格式的计算效率进行了比较，给出了超音速占优混合流动的直接分割求解方法。计算采用代数湍流模型，跨音速段的计算结果同其它实验和计算结果进行了比较。计算中发现对出口边界条件进行亚音速、超音速和有回流分别处理很有必要；边界层内的M数沿壁面法向变化非常快；不同壁温条件和不同的流动Re数均会改变边界层的厚度。     

复杂构形的亚、跨、超音速定常及非定常气动力的Green函数算法综述

对亚、跨、超音速定常及非定常气动力的 Green函数算法作了简要回顾 ,包括该方法的基本原理、近年来在国内的发展及其在气动弹性和气动载荷计算中的应用等。算例表明该方法适用范围广、计算效率高且比较准确