雷诺数三角翼绕流的影响

应用梁色液流动显示技术研究了雷诺数对６０°尖前缘三角翼前缘涡破裂位置、背风面流动结构等的影响，并详细分析了前风面流动随攻角的变化。     

三角翼前缘涡的某些破裂形式及特性研究

基于流动显示和PIV技术测量的实验结果，对三角翼前缘涡破裂的一些形式和破裂特性进行了分析和讨论。通过PIV测量所得到的涡量分布证实了在螺旋破裂的情况下，涡核的螺旋方向与前缘涡的旋转方向相反，及双螺旋破裂形式的存在等。进而对螺旋波的形成机理提出了与有关文献不同的看法。     

褶皱幅度对仿生翼型滑翔气动特性的影响

蜻蜓翅膀具有独特的褶皱状形貌．研究者们致力于利用仿生学原理,设计在低雷诺数条件下具有更优气动性能的褶皱翼型．本文采用计算流体力学方法,求解二维不可压Navier-Stokes方程组,探讨了四种翼型（平板翼型、流线翼型、小幅度褶皱翼型和大幅度褶皱翼型）的气动表现．在低雷诺数条件下得到以下结果：(1) 较小幅度的褶皱结构有利于增加升力和减小阻力．(2) 雷诺数变化时褶皱翼型的升力系数呈非线性变化;在特定雷诺数区间,幅度相近的褶皱翼型会发生相对气动优势的转变．(3) 褶皱结构内的回流区通过减小粘性阻力,使得翼型总阻力下降．(4) 翼型前缘的极小区域会产生脉冲高升力,对升力表现产生较大影响．这些结果表明,调整褶皱幅度是实现褶皱翼型气动优化的有效方案．     

尖顶弯折三角翼气动特性研究

对75°后掠三角翼在40%根弦处弯折引起的升阻特性及俯仰特性的变化进行了低速风洞实验研究,初步分析了尖顶弯折对三角翼气动特性的影响趋势. 得出尖顶弯折在低速大攻角下,可以使升力系数Cy增大、升阻比提高,并推迟失速攻角amax,从而说明顶部弯折对三角翼气动特性的改善作用是显著的.     

低雷诺数翼型蒙皮主动振动气动特性及流场结构数值研究

针对低雷诺数(Re)翼型气动性能差的特点,文章通过对翼型柔性蒙皮施加主动振动的方法,提高翼型低Re下的气动特性,改善其流场结构.采用带预处理技术的Roe方法求解非定常可压缩Navier-Stokes方程,对NACA4415翼型低Re流动展开数值模拟.通过时均化和非定常方法对比柔性蒙皮固定和振动两种状态下的升阻力气动特性和层流分离流动结构.初步研究工作表明在低Re下柔性蒙皮采用合适的振幅和频率,时均化升阻力特性显著提高,分离泡结构由后缘层流分离泡转变为近似的经典长层流分离泡,分离点后移,分离区缩小.在此基础上,文章更加细致研究了柔性蒙皮两种状态下单周期内的层流分离结构及壁面压力系数分布非定常特性和演化规律.蒙皮固定状态下分离区前部流场结构和压力分布基本保持稳定,表现为近似定常分离,仅在后缘位置出现类似于卡门涡街的非定常流动现象.柔性蒙皮振动时从分离点附近开始便产生分离涡,并不断向下游移动、脱落,表现为非定常分离并出现大范围的压力脉动.蒙皮振动使流体更加靠近壁面运动,大尺度的层流分离现象得到有效抑制.      

起始攻角附件三角翼涡破裂的演化特性研究

用数值方法求解Ｎ－Ｓ方程，并用“亚迭代”方法保证时间方向二阶精度，数值模拟了 ７５° 后掠三角翼非定常绕流，给出了详细的涡破裂起始攻角附近涡破裂的形成、发展和演化过程． 涡破裂从开始的螺旋破裂变为向泡状破裂演化的“过渡态”，然后发展到泡状破裂，再由泡状 破裂变为向螺旋破裂演化的“过渡态”，最终发展到典型的螺旋破裂形态．从计算结果中还发 现，当泡状破裂形成后破裂点前移速度明显减慢，同时当从泡状破裂向螺旋破裂演变时，破裂 点出现后移现象．可从最近相关的实验和数值模拟结果中发现这些现象的存在．     

飞行器纵向俯仰气动特性计算研究

应用有限体积方法求解三维可压缩雷诺平均N-S方程,计算了巡航导弹外形飞行器作小振幅俯仰运动时的动态绕流流场和空气动力特性,开展了导弹绕不同转轴、以不同频率和在不同迎角范围内进行俯仰运动的非定常气动力迟滞特性研究。计算结果表明,当导弹作快速俯仰运动时,在上仰和下俯过程中的同一迎角瞬间,绕导弹流场流动明显不同,表现出明显的非定常迟滞特性。导弹的非定常气动力迟滞特性随俯仰运动频率的增大明显增强,且气动力迟滞曲线随着俯仰轴位置的变化而变化。在同一减缩频率下,导弹在不同迎角范围内作周期俯仰运动时,相同的运动相位角所对应的升力系数对迎角的导数是一致的,而不同减缩频率下升力系数对迎角的导数随运动相位角变化曲线明显不同。     

超高雷诺数下矩形截面气动特性研究

为研究工程常见的超高雷诺数下矩形柱流场特性和气动特性,以常见的起重机单箱梁截面(梁高h=1m,雷诺数Re=6.85×105)为例,采用雷诺平均RANS Realizable k-ε方法,对中等湍流强度来流下矩形柱绕流进行数值模拟。首先进行网格无关性验证,然后选取宽高比为0.6、0.8、1.0、1.5、2.0共5组不同截面的矩形柱结构进行流场分析与研究。通过对比不同宽高比柱体截面周围的气动特性参数,以及尾流区的尾涡、回流参数等特征量表明:阻力系数、平均背压系数、升力系数均方根、斯托罗哈数等均随着宽高比增大而降低,尾涡宽度和距离随着宽高比增大而减小,回流长度和最大回流速度随宽高比增大而增大。同时,通过与相关雷诺数(Re=104～105)的风洞实验和仿真结果进行对比,发现其在数值上和参考值相差很小,说明高雷诺数对绕流参数的影响很小,可为相关工程设计提供参考。     

偏转矢量喷流对双三角翼前缘涡破裂影响的实验研究

本文应用染色液流动显示技术对后缘偏转喷流情况下76°/40°双三角翼前缘涡破裂位置的变化进行了观测,实验结果表明偏转喷流主要推迟与喷流方向相同一侧前缘涡的破裂,而使另一侧前缘涡破裂略有提前.随着喷流偏转角度的增大,喷流使两前缘涡破裂位置差逐渐增大.另外,随着模型攻角的增大,前缘涡涡核与双三角翼翼面的夹角逐渐增大,导致偏转喷流的作用逐渐减弱.     

高速列车空气动力学研究技术综述

随着我国高速铁路的快速发展,高速列车运行速度越来越快,包括气动阻力、横风效应、会车效应、隧道效应和气动噪声等一系列空气动力学问题日益突出. 利用模型试验、实车测量和数值计算等不同的研究手段,开展全面的高速列车空气动力学研究显得十分重要. 本文比较全面和系统地介绍了国内外高速列车空气动力学研究在模型试验、实车测量、数值计算等方面的技术现状和进展情况,对今后的发展方向和内容进行了展望.     

结构变形对气动力影响的计算分析

本文介绍了笔者在静气动弹性效应的理论与实验研究中的若干进展。静气动弹性研究涉及作用于飞行器的气 力与结构变形的相互作用,是研究飞行器的风洞试验数据和飞行试验数据相关性的一个重要方面。文讨论了在进行静气动弹性分析时对一些总理２的处理并对风洞试验模型及某型导弹进行了计算,结果表明,对于所选算例,结构变形使升力线斜率降低,压心前移。     

基于绝对节点坐标的多柔体系统动力学高效计算方法

绝对节点坐标法已经被广泛应用于柔性多体系统的动力学研究之中, 但是其计算效率问题尚未得到很好的解决. 基于绝对节点坐标方法计算弹性力及其对广义坐标的偏导数矩阵(Jacobi矩阵), 通常是基于第二类Piola-Kirchhoff应力张量来完成, 计算效率不高.根据虚功原理并采用第一类Piola-Kirchhoff应力张量的方法直接推导得到了弹性力及其Jacobi矩阵的解析表达式. 基于不同方法所得的数值算例结果对比研究表明, 该方法可使计算效率大大提高.      

AERODYNAMIC OPTIMIZATION OF HIGH-SPEED TRAIN BASED ON RBF MESH DEFORMATION

基于局部型函数三维参数化方法、改进的蚁群算法和改进的克里金（Kriging）代理模型,开展了列车头型的三维气动减阻优化设计研究。为了避免复杂几何外形大变形情况下千万量级网格的重复生成,提高高速列车头型优化设计的效率,引入了缩减控制点的径向基函数网格变形技术。优化结果表明：径向基函数网格变形技术在不降低网格质量的情况下可以有效缩短网格变形的时间消耗,能够用于复杂几何外形的气动优化设计;在给定的设计空间内,控制鼻锥外形的6个关键设计参数对列车气动阻力的影响呈单调递增关系;优化后,在满足约束条件的情况下,简化外形列车的整车气动阻力减小5.41%,头尾车减阻效果明显,中间车气动阻力基本不变。     

基于非递推公式的机械臂正动力学的并行计算策略

利用Ｊｏｕｒｄａｉｎ原理来消除约束反力，并通过引入Ｌａｇｒａｎｇｅ乘子释放约束，得到了机械臂正动力学非递推形式的计算模型．由于在这种模型中引进了一些冗余计算，因此减弱了方程中各计算量之间的依赖性，从而提高了模型的内在并行性．为了尽量减少并行计算所需处理器的数目，本文采用了求解规模缩减技术，并基于这个模型提出了一种面向Ｏ（ｎ）个处理器的机械臂正动力学的并行计算策略．最后以ＰＵＭＡ５６０机器人的前三个臂为例进行了计算效率分析     

MULTI-OBJECTIVE OPTIMIZATION AND AERODYNAMIC PERFORMANCE ANALYSIS OF THE UPPER SURFACE FOR HYPERSONIC VEHICLES

为分析小攻角巡航条件下吸气式高超声速飞行器上壁面的变化对其气动性能和容积的影响, 以参数化后的飞行器上壁面对称面型线为设计变量, 在飞行马赫数6.5, 飞行高度27 km, 飞行攻角为4°的条件下, 采用计算流体力学为性能分析工具, Pareto多目标遗传算法为优化设计方法, 开展了二维条件下的升阻比/容积双目标优化设计. 在此基础上, 选择典型的二维优化结果, 重构生成对应的三维构型并进行数值分析, 获得了飞行器气动性能和容积间的相互关系. 结果表明在巡航条件下, 尽管二维/三维条件下飞行器的气动参数数值有较大差别, 但在这2种条件下, 飞行器的升阻比和容积间的关系均近似呈线性反比例关系. 同时, 对于三维构型而言, 在给定容积不变的条件下, 通过改变上壁面对称面型线的形状仅能使升阻比获得较小的增量(约0.36%). 相比之下, 当给定升阻比基本不变的条件下, 飞行器容积可调空间相对较大, 约为1.93%. 此外, 计算结果还表明, 在飞行器的容积基本不变情况下, 通过调节上壁面对称面型线, 可使飞行器的俯仰力矩获得5%左右的调节空间, 且其升阻比基本不变.     

蜻蜓滑翔时柔性褶皱前翅气动特性分析

蜻蜓是自然界优秀的飞行家,滑翔是其常见且有效的飞行模式.蜻蜓优异的飞行能力来源于其翅膀的巧妙结构,褶皱是蜻蜓翅膀上最为显著的结构之一,不仅提高了翅膀的刚度,还改变了其气动特性,而飞行过程中柔性翅膀会产生变形是蜻蜓翅膀的另一特性.为揭示蜻蜓在滑翔时,柔性褶皱前翅的变形,探究褶皱和柔性的共同作用对其气动特性的影响,基于逆向工程,依据前人的测量数据和研究成果,通过三维建模软件建立了蜻蜓三维褶皱前翅的计算流体力学(computational fluiddynamics,CFD)模型和计算结构力学(computational structuralmechanics,CSD)模型,并通过模态分析验证了此模型有足够的精度.基于CFD方法和CFD/CSD双向流固耦合计算方法分别对蜻蜓滑翔飞行时刚性和柔性褶皱前翅的气动特性进行了数值模拟,结果表明,柔性褶皱前翅受气动载荷后,翅脉和翅膜产生形变,柔性前翅上下表面压力差相较于刚性前翅减小了,从而其升力和阻力也减小了,而在大攻角时,变形后的前缘脉诱导出比刚性前翅更强的前缘涡.因此在攻角小于10$^\circ$时刚性前翅的气动特性优于柔性前翅,继续增大攻角,柔性前翅的气动特性则优于刚性前翅.前翅受载后气动响应时间短,翅尖的变形最大,仅仅产生了垂直于翅膀所在平面方向上的变形,而没有发生扭转,翼根处受到应力最大,褶皱上凸部分承受蜻蜓滑翔时前翅的主要载荷.      

基于响应面法的新型排翼式飞艇的气动优化设计

采用均匀试验设计的方法, 以后掠排翼式飞艇为初始气动外形, 后翼的轴、法向坐标为设计变量, 升力、阻力系数和升阻比为目标函数, 通过CFD结合响应面近似理论, 建立了气动参数的二阶响应面优化模型, 并针对排翼式飞艇的升阻特性进行了优化设计. 研究结果表明: 响应面方法是解决气动优化问题的一种行之有效的方法;优化后排翼式飞艇的气动性能得到明显改善, NACA0020, NACA0030, NACA0040排翼式飞艇的升阻比分别提高了9.61%, 6.08%, 13.08%; 飞翼的轴向和法向位置对各气动参数的影响不尽相同, 升力系数C_rm l和升阻比k对轴向位置更为敏感, 而轴向和法向位置对阻力系数C_d的影响则基本相当.      

CFD模拟方法的发展成就与展望

本文综述了计算流体力学(computational fluid dynamics, CFD),尤其是计算空气动力学的发展概况.从计算方法、网格技术、湍流模型、大涡模拟等方面分别总结了CFD所取得的成就,分析了当前存在的问题、困惑,展望了其发展趋势.在CFD计算方法中,主要介绍了中心格式、迎风格式、TVD格式、WENO格式、紧致格式以及间断Galerkin有限元方法,对不同方法的原理和特性进行了系统阐述.网格技术包括结构网格、非结构网格、混合网格以及重叠网格,重点讨论了重叠网格的若干关键技术.在湍流模型中,对目前的模型进行分类介绍,包括线性涡黏性模型、二阶矩模型、非线性模型等,还介绍了转捩模型、DES方法以及SAS方法等.在大涡模拟方法中,就其中若干相关的研究方向进行了探讨,包括滤波方法、亚格子模型、收敛标准、数值格式等. 文中还包含了作者在相关领域的若干研究成果.     

INVESTIGATIONS ON THE MECHANICAL CHARACTERISTICS OF NANOBEAMS BASED ON THE NONLOCAL STRAIN GRADIENT THEORY1)

基于非局部应变梯度理论,建立了一种具有尺度效应的高阶剪切变形纳米梁的力学模型. 其中,考虑了应变场和一阶应变梯度场下的非局部效应. 采用哈密顿原理推导了纳米梁的控制方程和边界条件,并给出了简支边界条件下静弯曲、自由振动和线性屈曲问题的纳维级数解. 数值结果表明,非局部效应对梁的刚度产生软化作用,应变梯度效应对纳米梁的刚度产生硬化作用,梁的刚度整体呈现软化还是硬化效应依赖于非局部参数与材料特征尺度的比值. 梁的厚度与材料特征尺度越接近,非局部应变梯度理论与经典弹性理论所预测结果之间的差异越显著.