某型客机风挡雨刷系统故障分析

针对某型客机风挡雨刷系统在飞行过程中被气流吹离正常停放位置的故障,利用CFD软件对该型客机风挡雨刷系统进行了气动力分析。同时利用气动力分析结果,找出了导致风挡雨刷系统发生故障的主要原因并设计了加装导流板的排故方案。最终,通过试飞对故障原因和排故方案进行了验证。     

非结构网格生成方法的改进及气膜冷却三维静子流场的求解

讨论五面体网格下,涡轮静子叶片含多排冷却孔时三维粘性流场雷诺平均Navier-Stokes方程组的数值求解问题。考虑到叶片含多排冷却孔的工程背景,本文发展了生成五面体非结构网格的快速有效方法即仅在一个拟S₁面上生成非结构网格,而沿叶高方向仍采用结构网格的生成办法。文中还采用了新的网格优化措施以改善生成网格的质量。流场的计算采用了非结构网格下的高分辨率迎风格式。几个典型算例表明:本文发展的网格生成方法快速有效,而且网格生成质量较高;发展的计算格式能够获得较为满意的三维流场数值解,并具有较高的激波分辨率。      

二维间断有限元混合网格混合算法研究

针对二维流场进行结构／非结构混合网格划分,在靠近物面处生成结构化网格,在结构化网格以外和远场之间以非结构化网格填充。使用间断有限元方法（DGM）离散Euler方程。在混合网格混合算法中,结构化网格区域使用结构DGM,非结构化网格区域使用非结构DGM,保证结构／非结构对接面处的通量守恒,并加入moment限制器消除激波前后伪震荡。对NACA0012和RAE2822翼型跨音速无粘流动的数值模拟结果说明了混合算法的有效性以及在捕捉间断方面的能力。     

应用混合网格分析复杂积冰翼型气动性能

通过发展一种点-点对接的结构/非结构混合网格生成方法,避免了复杂积冰外形难以划分结构化网格的问题,并对复杂积冰翼型的气动性能进行了分析计算.在靠近积冰边界的内层采用了非结构网格以拟合复杂的边界,在非结构网格外采用结构化网格以节省存储空间和计算时间.为了检验网格对选用湍流模型的影响,整个流动区域分别采用SST和SA湍流模型,求解了雷诺平均N-S方程.数值计算结果表明SST模型更适于模拟复杂分离流动,积冰对翼型的气动性能造成了严重的影响.     

飞机风挡雨刷系统地面试验方法研究

提出了一种进行飞机风挡雨刷系统地面试验的试验方法。提出等效降雨强度的概念,并按照雨量相等的原则,把飞机在飞行过程中主风挡玻璃遭遇的降雨强度转换成风挡雨刷系统地面试验所需的等效降雨强度,为风挡雨刷系统的地面试验提供一个较为真实的降雨环境,解决了进行风挡雨刷系统地面试验的关键问题。通过工程实践,证明本方法是可行且有效的。     

类天宫航天器迎风面积建模及变化特性分析

针对类天宫航天器复杂结构外形迎风面积的精确计算问题,提出了一种在轨飞行期间迎风面积的计算模型。该模型考虑航天器姿态变化和帆板转动因素,在航天器本体系下建立经纬度网格,将三维的投影方向矢量转化为二维的经度、纬度;基于图形学几何投影方法,对投影经度、投影纬度、帆板转角模型参数按照一定颗粒度离线穷举,计算各种参数状态的迎风面积,离线建立航天器迎风面积模型;轨道计算过程中,根据航天器在轨飞行模式确定该时刻的姿态和帆板转动角度,通过迎风面积模型插值获得迎风面积。应用该方法计算天宫一号长期在轨飞行阶段的迎风面积,对不同时间尺度的变化规律进行了分析,验证了方法的有效性。     

生成非结构网络的改进阵面推进法

通过对非结构化网格生成的阵面推进方法进行分析,对该算法进行了优化设计,采用笛卡尔直角坐标网格来代替背景非结构化网格,其上网格尺度参数通过控制体有限元法求解Possion方程的Dirchlet问题得到,提高了网格生成的效率和质量。同时利用Laplacian滤波算子和对角交换方法来改善三角形网格的质量,尤其是对角交换技术对消除畸形三角形单元效果明显。对固壁边界采用三次样条拟合后重新离散的方式形成初始阵面,改进了网格的贴体性,提高了网格的光滑性。在编制通用程序时利用链表结构和派生类型的数据变量,节省了计算内存,提高了网格生成效率。     

高超声速热化学非平衡钝体绕流数值模拟

对高超声速热化学非平衡钝体绕流进行研究,发展了高温热化学非平衡条件下的高超声速气动热数值模拟方法。采用二阶精度迎风TVD格式的N-S方程有限体积法多块结构网格求解器,考虑了空气5组分及化学非平衡效应的影响,对组分生成源项及振动能源项采用了点隐式的处理,提高了计算效率。对球锥模型和RAM-CII飞行实验等算例进行了数值气动热模拟,分析了不同非平衡模型在热化学非平衡条件下流场的影响。计算结果表明建立的数值模拟方法具有较高的精度。     

基于混合网格的声爆/气动一体化设计方法研究

探索了基于混合网格的近场预测以及基于波形参数法的远场预测相结合的声爆预测方法,综合利用了结构化网格计算效率高以及非结构网格对复杂外形适应性好的优点,在保证声爆预测精度的基础上提高了计算效率。在此基础上,建立了基于声爆预测方法、响应面模型和Pareto遗传算法的声爆/气动一体化设计方法。分析表明,发动机短舱对超声速升阻特性和声爆特征都有较大影响,利用声爆/气动一体化设计方法对发动机短舱位置进行了一体化设计研究。结果表明,建立的响应面模型的精度完全满足要求,根据不同的设计偏向选择了三种优化方案,优化后方案的性能较原始方案都有不同程度的提高。     

三维分块非结构化网格上的可压缩流计算方法

提出了分块非结构化网格生成方法,通过子块划分、子块网格划分、整体网格拓扑连接生成可局部加密的分块非结构化网格,既具备了非结构化网格对复杂区域的处理能力,又避免了分块结构化网格块间信息传递的复杂计算过程.通过将密度修正值引入压力修正方程,发展了非结构化网格上可压缩流动的SIMPLE算法.对孤立翼型及透平叶栅的流动进行了数值计算,压力系数、速度、折转角的数值结果与已有的实验值吻合较好,验证了数值方法的精度.      

过渡状态的倾转旋翼气动特性计算分析(英文)

针对特殊的旋翼倾转运动,建立了一个过渡状态的旋翼非定常气动力数值计算方法.为正确模拟旋翼倾转运动使桨叶受到的附加惯性力及哥氏力作用,重新推导了旋翼倾转时的桨叶气动力模型和挥舞运动方程;为了适合于旋翼倾转时的入流和气动力计算,入流模型中考虑了倾转运动引起的旋翼尾迹弯曲影响.应用建立的方法,首先进行了旋翼配平计算,以验证计算模型,并给出了倾转旋翼的操纵量.然后,着重计算了旋翼在倾转过渡时的拉力.俯仰和滚转力矩随倾转角的变化,分析了倾转飞行时的前飞速度.倾转时间等对旋翼气动力的影响,得出了一些新的结论.     

NNW-HyFLOW高超声速流动模拟软件框架设计

NNW-HyFLOW在国家数值风洞（NNW）工程支持下,由风雷开源软件提供基础框架支撑,将建设成为基于结构/非结构混合网格,面向高超声速应用领域的国产自主工业CFD软件,具备高温气体热化学非平衡效应模拟及其相关气动力、气动热和气动物理特性计算分析等主要特色功能。本文从框架设计、数据结构、耦合方法、并行计算方法以及接口设计等方面对软件设计思想和框架特点进行了介绍,给出了求解器采用的理论模型、核心数值算法及其实现方法,结合HEG风洞试验、RAM-C飞行试验、Electre飞行试验以及航天飞机OV102飞行试验等典型算例开展了数值模拟。研究表明：NNW-HyFLOW具有底层代码复用、功能兼容性好、扩展能力强和接口灵活等优点,其当前测试版本已经具备了较好的高超声速非平衡流动数值模拟能力,在热化学非平衡效应及其影响的气动力特性、气动热环境和等离子体分布特性预测与评估方面,具有较高的数值计算精准度,初步满足了高超声速复杂飞行器高温非平衡流动数值模拟的需求。     

弧翼飞行器柔性飞行动力学研究

以一种弧翼布局飞行器为背景,介绍了飞行器的特性;建立了线性分布式气动力模型,添加到全机的三维有限元柔性模型上,进行飞行动力学的仿真分析。研究了分布式气动力模型对弧翼布局飞行器的柔性飞行动力学特性造成的影响。通过对不同响应情况下的飞行状态的仿真计算与分析,总结出了采用分布式气动力建模的大展弦比弧翼飞行器的柔性飞行动力学的变化规律。     

基于非线性试验气动力的飞机静气动弹性响应分析

分别使用线性和非线性静气动弹性分析方法,对某飞机纵向静气动弹性响应特性随飞行动压、攻角、平尾偏度和纵向过载变化的趋势进行了分析,并将分析结果和飞行试验进行了比较。线性方法的气动力计算使用亚音速偶极子格网法。非线性方法由风洞试验提供刚体气动力,并使用线性气动力影响系数矩阵对其进行弹性化处理。通过对两种方法的计算结果进行对比分析可以看出:①使用刚体试验气动力的非线性分析方法能够获得和飞行试验比较一致的结果;②在线性方法所提供的结果中,翼面的弯曲变形及剪力、弯矩和扭矩等部分参数能为初步设计提供大致参考,但气动力系数的弹性增量、翼面的扭转变形、翼面的剪力、弯矩和扭矩的弹性部分可能不能为型号设计提供正确的指导;③不论刚体气动力为线性还是非线性,气动力系数随动压变化的趋势均呈线性。     

一种适合于旋翼前飞非定常流场计算的新型运动嵌套网格方法

针对旋翼流场计算中关键的运动嵌套网格技术,提出了一个新的网格单元关系快速判断方法。该方法通过计算特殊网格点在背景网格上的位置来建立桨叶网格在背景网格中所占据范围的关系,并用于嵌套网格洞单元的识别,又进一步将该方法推广用于Inverse Map的建立,在此基础上,建立了一套用于前飞流场计算的运动嵌套方法。通过多种前飞旋翼测试算例,分别评估了背景网格、桨叶网格尺寸改变对计算时间的影响,同时验证了该方法的可靠性,结果表明构建网格嵌套关系所用时间随桨叶网格和背景网格尺寸的变化幅度较小,且具有良好的鲁棒性,能够满足旋翼前飞状态计算的需要。最后,在基于前飞非定常N-S/Euler流场求解器中,采用该方法模拟了UH-60A和7A旋翼悬停/前飞状态的气动特性,通过计算时间的对比体现出该方法的高效性,并且流场计算结果与试验数据对比表明该方法可以有效地用于旋翼前飞非定常流场的计算。     

转静干涉对叶片非定常表面压力的影响

为了研究转静干涉对叶片表面压力的影响,通过数值模拟的方法,计算了一级低速轴流压气机流场,得到了充分收敛的周期性非定常流场解。在此基础上分析了叶片表面压力及气动力的变化规律,并通过改变背压,研究了流量对叶片表面压力和气动力的影响。最后在低速轴流压气机实验器上进行了实验验证,结果表明计算和实验的结果吻合得较好。     

Mechanism analysis of hysteretic aerodynamic characteristics on variable-sweep wings

Variable-sweep wings have large shape-changing capabilities and wide flight envelops, which are considered as one of the most promising directions for intelligent morphing UAVs. Aerodynamic investigations always focus on several static states in the varying sweep process, which ignore the unsteady aerodynamic characteristics. However, deviations to static aerodynamic forces are inevitably caused by dynamic sweep motion. In this work, first, unsteady aerodynamic characteristics on a typical variable-sweep UAV with large aspect ratio were analyzed. Then, deep mechanism of unsteady aerodynamic characteristics in the varying sweep process was studied. Finally, numerical simulation method integrated with structured moving overset grids was applied to solve the unsteady fluid of varying sweep process. The simulation results of a sweep forward-backward circle show a distinct dynamic hysteresis loop surrounding the static data for the aerodynamic forces. Compared with the static lift coefficients , at the same sweep angles, dynamic  lift coefficient in sweep forward process are all smaller, while dynamic sweep backward  lift coefficient are all larger. In addition, dynamic deviations to static  lift coefficient are positively related with the varying sweep speeds. Mechanism study on the unsteady aerodynamic characteristics indicates that three key factors lead to the dynamic hysteresis loop in varying sweep process. They are the effects of additional velocity caused by varying sweep motion, the effects of flow hysteresis and viscosity. The additional velocity induced by sweep motion affects the transversal flow direction along the wing and the effective angle of attack at the airfoil profile. The physical properties of flow, the hysteresis and viscosity affect the unsteady aerodynamic characteristics by flow separation and induced vortexes.     

部分三代战机“迎角悬挂”现象的机理分析与仿真计算

参考有关风洞试验结果,探讨了部分三代战机"迎角悬挂"现象的机理,对比分析了不同机型"迎角悬挂"的特点。并以某型三代战机为例,根据其纵向气动和操稳特性曲线及飞控系统特性,利用时域动态响应的方法,仿真计算了该型机大迎角下"迎角悬挂"的有关现象,计算结果与其技术资料的叙述相符。分析与计算表明,大迎角下"迎角悬挂"现象,是由该类型三代战机气动布局决定的一个共性问题;同时,"迎角悬挂"现象的特点又与具体机型的气动布局特点以及重心位置等因素有关。     

近距编队飞行气动影响建模与分析

针对近距编队飞行面临的相互气动影响问题,采用等效气动效应法建立了尾流气动影响的等效模型.该模型可引入编队飞行六自由度方程.通过仿真验证表明,由等效模型计算出的诱导力与诱导力矩和风洞数据具有较好的符合度,等效模型可直接用于编队飞行建模及控制器的设计与仿真.     

基于目标压力分布的旋翼先进气动外形反设计分析方法

基于雷诺平均Navier-Stokes（RANS）方程、运动嵌套网格、目标压力分布及余量修正方法,构建了一套直升机旋翼桨叶先进气动外形反设计方法。为避免由桨叶气动外形变化导致的网格畸变,发展了一套基于Poisson方程求解的旋翼桨叶结构化贴体正交网格的快速、自动化生成方法,为提高运动嵌套网格的生成质量和通用性,采用剖面间网格插值与桨尖翻折相结合的方法,同时建立了基于“Top Map”和“Inverse Map”相结合的洞边界划定和贡献单元搜寻的新方法。基于Navier-Stokes方程和双时间法建立了旋翼非定常流场模拟方法,通量求解采用Roe-MUSCL格式,并使用低速预处理法来克服前飞旋翼流场收敛中遇到的刚性问题。在计算流体力学（CFD）方法基础上,基于旋翼翼型压力系数余量联立各方位角处的反设计MGM（Modified Garabedia-McFadden）超定方程组,并依据激波分离、失速等约束设置了各方位角处的反设计权重系数,创建了基于MGM超定方程组最小二乘解的旋翼气动外形（翼型）设计方法。应用所建立的方法,分别针对多目标、多状态和前飞时的旋翼（翼型）气动外形进行反设计分析,验证了本文方法的有效性。最后,将该方法拓展应用到旋翼桨尖气动外形设计中,设计得到与UH-60A直升机旋翼气动特性相似的矩形桨叶外形。