海鸥翼折转运动的数值模拟及分析

扑翼的非定常流动及气动特性分析对微型飞行器设计具有重要意义。本文采用动态混合网格技术以及非定常数值计算方法,对海鸥翼的拍动过程进行了数值研究。采用基于径向基函数的插值技术实现网格的变形,为了提高插值效率,发展了基于最大物面误差的参考点选择算法。基于文献观测数据建立了海鸥翼的三维模型,并设计了相应的拍动以及变形规律,对拍动角、折转角的影响进行了分析。分析结果表明,折转角可以减小上拍过程的不利影响,对提高整个扑翼周期的时均升力、减少时均阻力和能耗是有益的。进一步通过调整折转的相位来增加翼的折转时间,并减少其展开的时间,对提高扑翼过程的平均气动力特性是有益的。     

多段翼型局部主动变形流动控制的非定常数值模拟

对30P30N三段翼型失速攻角附近的分离流动进行了数值模拟研究。为了抑制大攻角时背风区的流动分离,在主翼段上表面引入了行波壁变形模型和抛物型局部主动振动模型,利用作者以往发展的动态混合网格技术和相应的非定常计算方法,对变形过程中的非定常分离流动进行了数值模拟,分析了各种变形参数对流动分离的影响。计算结果表明,在适当的条件下,局部主动变形能够抑制翼型背风区的分离,由此可以起到增升减阻的作用,改善翼型的气动性能。     

粘弹性蜻蜓型翼挥拍的流固耦合求解途径初探

昆虫翼挥拍时的动态变形与受到的惯性力和气动力是相互影响的,这是一个复杂的流固耦合过程,本文利用简化的非定常气动力估算方法、简化的昆虫翼结构模型和简化的拍动方式,对该问题的耦合求解途径做了初次探索.首先,在比较了昆虫翼拍动中受到的惯性力和气动力大小的基础上,阐明了耦合求解蜻蜓翼变形过程的必要性.其次,提出一个迭代求解该耦合问题的简化方案,并利用有限元数值模拟了粘弹性蜻蜓模型翼绕定轴往复挥拍的变形过程,比较了耦合气动力前后的模型翼变形,并初步分析了能耗变化.     

微型飞行器低雷诺数矩形扑翼非定常气动特性的数值模拟

采用双时间步方法求解三维可压缩非定常N-S方程,数值模拟了微型飞行器低雷诺数矩形扑翼的非定常绕流,首先将得到的结果与文献进行了对比,数据间具有较好的一致性.然后针对不同的展弦比、减缩频率及初始攻角,计算了矩形扑翼的非定常气动特性及表面流态和动态压力分布,并分析了翼尖涡对扑翼非定常气动特性的影响.     

可变形翼战术导弹飞行动力学联合仿真研究

基于集中质量假设发展的一套适于可变形翼导弹的多体运动模型,以附加惯性项及快变气动参数的形式来反映弹翼变形对导弹动态特性的影响,涉及交叉求解不同弹翼外形对应的气动特性及运动方程.由此提出的一种基于Missile DATCOM与MATLAB/SIMULlNK的联合仿真方法,应用于典型轴对称导弹的变翼动态特性研究中,快速有效...     

Z型翼变体飞机的纵向多体动力学特性

机翼变形时,变体飞机的翼面积、惯性特性、全机焦点和重心位置等均会发生较大的变化,从而引起飞机的动态特性也随之改变。为此对机翼变形过程中的Z型翼变体飞机进行了纵向多体动力学建模仿真;推导了变形过程中变体飞机的六自由度非线性动力学方程,并通过简化得到了解耦后的纵向动力学方程。机翼折叠动态过程的气动特性数值模拟结果表明,不同折叠角速度下飞机的气动力相差不大。在机翼折叠角速度较小且忽略非定常气动效应的情况下,采用气动力准定常假设对变形过程中不同机翼折叠角速度下变体飞机的纵向响应进行了数值仿真,并研究了重心位置移动和气动特性变化对飞机变形过程动态特性的影响规律。结果表明,折叠过程中气动特性的变化是影响飞机动态特性的主要因素,机翼折叠后飞机的速度和迎角增加,且飞行高度下降较大。     

柔性扑翼全机气动结构耦合特性及稳定性研究

由于微型扑翼飞行器具有体积小、重量轻等特点,在飞行过程中会产生明显的柔性变形,因此在对扑翼气动特性进行数值模拟时,有必要考虑柔性变形的影响。基于结构动力学理论,发展了一种适用于扑翼全机气动结构耦合特性的数值模拟方法与一种适用于扑翼的结构动力学求解方法。方法利用哈密顿原理,对动能和应变能进行变分,进而得到扑翼动力学方程,采用结构有限元方法对运动方程进行离散并求解;采用基于嵌套网格的CFD求解器对微型扑翼全机非定常绕流进行数值模拟。采用CFD/CSD耦合求解器对微型扑翼飞行器全机气动结构耦合特性进行数值模拟,分析了结构变形对柔性扑翼气动特性的影响,并分别对扑动过程中刚性扑翼和柔性扑翼的压心变化范围进行研究,分析了结构变形对柔性扑翼稳定性的影响。     

蝙蝠飞行的空气动力学研究进展

蝙蝠飞行的空气动力学研究是近二十年来关注的热点问题,其中蕴含着翼主动大变形的流动控制机理,对其进行研究可以揭示不同于鸟类和昆虫的非定常空气动力产生的机制。首先,根据研究工具的进步程度,把对蝙蝠飞行的研究分为了三个历史阶段,即在第一和第二阶段通过获得部分"静态"的数据来研究,在第三阶段对蝙蝠翅膀飞行中的动态变形进行测量和流场结构测量。然后陈述了蝙蝠翼的特征,飞行中挥拍运动及变形规律特征,并简述了特征的四类变形模型(扭转变形、弦向弓形变形、展向弯曲变形和翼面积改变)。有了初步的变形模型后,对主动变形产生非定常空气动力产生机制和气动性能调控进行了探讨。最后提出了当前对蝙蝠飞行的空气动力学研究面临的困难和挑战。     

基于动态混合网格的鱼体巡游数值模拟

对鱼类等柔性变形体的数值模拟需要解决如下几个问题(1) 柔性体动态网格生成;(2) 基于动态网格的不可压流非定常算法.本文基于Delaunay背景网格插值方法,发展了一套快速生成柔性体变形的动态混合网格生成方法;对于不可压流问题,采用虚拟压缩方法进行求解,结合双时间步方法来处理非定常问题;最后数值模拟了鱼体巡游的非定常流动,分析了鱼体的摆动频率对力、功耗等流体动力学特性的影响.     

飞机着陆过程中非定常地面效应数值模拟方法

针对飞机着陆过程的大位移运动,提出适于多块结构动网格生成的改进谢别德-无限插值混合方法,改进谢别德方法准确计算各网格块的边界网格变形,无限插值方法快速计算内部网格点变形,结合局部重构保证动态网格的合理分布.混合方法有效解决了单一无限插值法存在的动网格交错问题.进而将二维非定常地面效应N-S方程计算方法发展至三维情形,应用于带增升装置的翼身组合体着陆过程中非定常地面效应计算分析.结果表明:模型的升力系数随离地高度的减小而减小,非定常和准定常计算方法的结果存在一定差异,而考虑了非定常效应的算法更为准确.     

跨声速机翼非定常气动力的全位势粘位迭代计算

采用Ｃ－Ｈ型网格,守恒型非定常全位势方程的时间精确近似因式分解差分地计算二维,三维的跨声速非定常位势流,用准定常,准二维方法计算边界层位移厚度,通过粘位迭代得到的跨声速翼型,机翼的非定常气动力,所得结果与实验数据吻合很好。     

二维非定常跨音速时间线化积分方程的推导

1.引言 求解非定常跨音速流动的主要困难是非线性问题。对于微幅翼型振动问题(同时引起激波的微幅振动)可做时间线化简化处理,得到时间线化微分方程定解问题。时间线化积分方程是由时间线化微分方程推导出来的。     

跨声速翼型绕流的Euler／边界层方程干扰数值解

本文利用Euler方程和可压缩湍流边界层积分方程研究绕跨声速翼型的有粘与无粘强干扰流动。应用有限差分法在贴体的网格上求解时间相关的Euler方程,以剪功积分方法求解翼面贴附和分离湍流边界层流动,并引入一个松弛方程描述剪应力对上游湍流历程的延迟响应。有粘/无粘干扰采用表面源模型。计算结果表明,对翼面存在强干扰流动情况,获得了与实验值基本吻合的结果。     

计算机翼跨声速绕流的全速位积分方程数值解

本文利用格林定理将全速位方程转换为全速位积分方程，应用高斯定理对沿激波面的显式积分进行简化。将全速位积分方程的解表示成为面元项和场元项之和，然后离散进行数值求解。本文将激波捕捉方法和激波装配技术结合起来，对机翼跨声速绕流进行数值实验，计算结果令人满意。     

计算绕薄翼型跨音速非定常流的积分方程法

 本文从小扰动方程出发推导出绕薄翼型的跨音速非定常流动的积分方程,为使此方程适用于具有激波的流场引入人工粘性,并对其数值求解。通过算例讨论了人工粘性和计算域大小对计算结果的影响,与实验及其它数值计算结果比较,表明本方法的准确度令人满意。计算量小,收敛性好是此法的特点。     

翼-身组合体跨音速绕流全位势积分方程数值计算

 用积分方程方法求解 Prandtl- Glauert算子表示的全位势方程 ,并计算了翼 -身组合体跨音速绕流。用 Murman- Cole差分格式计算空间场源强度 ,以捕捉激波。计算结果与相应的实验结果符合良好     

Euler方程与边界层积分方程耦合求解跨音速翼型绕流

 应用Euler方程求解跨音速翼型特性时考虑了粘性影响,粘性影响是通过边界层动量和能量积分方程求解的,即粘流/无粘流迭代方法。其中Euler方程采用LU-ADI方法求解;边界层方程均由正解法过渡到反解法,以解决强激波干扰区出现小分离泡的计算问题。计算中使用了贴体C网格,通过一定变换使其保持基本正交。计算结果表明,压力分布、摩阻系数分布与实验结果符合较好。     

亚音速升力面理论的一种新解法

本文提供了亚音速升力面理论的一种新解法,可求解亚音速流中有限翼展机翼的载荷分布。计算中采用传统的弦向级数和展向分条。证明了Mangler二重奇异积分可交换积分次序,从而使积分可解析地求出。对弦向积分的被积函数作多项式分段逼近,可使核函数的表达式成为封闭形式的有限和,导致性态很好的系数矩阵。计算结果与实验及其它方法的结果作了比较,并给出了计算效率。     

一种求解核函数方程的新方法—辅助核法

本文提出了求解核函数方程的辅助核法。此法克服了偶极子栅格法对网格分布的限制,具有更大的适用范围。数值计算表明,此法对求解绕翼型的跨声速激波流动是很有利的。     

时间推进积分方程法求解非定常跨声速流动

本文提出了一种求解非定常跨声速流动的新方法——时间推进积分方程法,此法克服了时间线化积分方程法的限制,能较好地模拟激波的运动。本文首先用一维波(?)问题——模型问题阐明此法的基本思想,然后将它应用于二维非定常跨声速流动中。本文还首次引入拟速度位的概念,使时间推进积分方程式得到简化,尾涡条件和Kutta条件更易处理。数值计算表明时间推进积分方程法是合理可靠的。