飞鸟、飞行昆虫、蝙蝠等自然界的飞行生物均采用扑翼方式飞行，尺度缩小后，扑翼飞行是一种可有效提高气动性能的飞行方式。微型扑翼飞行MicroAirVehicle，简称FMAV）是一种仿生飞行器，但是目前人造的微型扑翼飞行器和飞行生物的飞行能力相比，还相差甚远。

　　近年来，FMAV受到广泛关注，包括美国、欧盟、日本以及中国在内的世界各国竞相对其展开研究，同时，也得到了诸多研究成果，涉及气动计算、风洞，常见的有“°”字形和“8”字形。

　　不同飞行生物的翼稍运动轨迹本文针对包含前后扫掠的扑动运动的气动特性展开研究，考虑到现有技术条件下仿鸟尺度飞行器具有更强的可实现性，选择仿鸟尺度的微型扑翼为研究对象，通过对其三个自由度的运动设计不同的参数，研究三维扑翼采用仿生“°”字形和“8”字形（为了研究方便，假设理想的“°”字形和“8”字形规律）运动时的气动特性。

　　1扑动方式定义按照转轴的不同，扑翼的运动可分为：上下扑动、弦向扭转和前后掠动（以下简称扑动、扭转和掠动），如所示，其中，扑动运动为绕轴的转动，扭转运动为绕轴的转动，掠动运动为绕轴的转动。扑动与扭转被认为是基本运动方式，当加入不同的掠动频率后，扑翼翼稍轨迹将体现出不同的形态，本节将分别定义。

　　1.1基本运动基本运动包括扑动和扭转。扑动运动规律为扭转运动规律为幅度为扑动频率；为在Z时间步的瞬时迎角，也称之为扭转角（包含扭转变形导致的迎角变化）；a.为初始迎角；am为大扭转角；z为展向站位；6为半展长。

　　扑翼按照基本运动规律的运动如所示，扭转与扑动的相位差为9.°，该相位差下推进效率（a）扑动和扭转的角度函数包含扑动与扭转规律的基本运动律为动规律下时刻在展向站位z的瞬时掠动角度；/，mS0字形或“8”字形运动规律下，在展向站位z的大掠动角度；/，m，lip为“0”字形或“8”字形运动规律下，翼稍大掠动角度。

　　1.2 o字形与8字形扑动在基本运动规律上，增加定的掠动规律，即可形成“0”字形或者“8”字形运动规律，这两种规律的区别在于掠动的频率不同。若掠动频率等于扑动频率，则翼尖轨迹体现为“0”字形；若掠动频率为扑动频率的两倍，则翼尖轨迹体现为“8”字形。当掠动频率更高时，翼尖轨迹则更为复杂，该扑动方式也比较少见，不具有代表性，本文不讨论其影响与规律。

　　律为（b）“0”字形掠动的运动规律“0”字形扑动时的掠动规律与扑翼运动示意字形掠动规律与翼运动示意图如（a）“8”字形掠动的角度函数（b）“8”字形掠动的运动规律“8”字形扑动时的掠动规律与扑翼运动示意数值方法采用数值模拟方法研究微型扑翼仿生复杂扑动情况下的气动特性。微型扑翼飞行器速度低、尺度小，其绕流具有低马赫数、低雷诺数的特点，受空气粘性的影响很大，故采用雷诺平均N-S（Reynolds-Averaged为基本控制方程。

　　在连续介质假设下，忽略彻体力和热源的三维非定常可压缩RANS方程可写为如下守恒形式二（u9v-，v）、p、E分别为流体的密度、速度在直角坐标系下的三个分量、压力、单位质量总能。

　　单位质量总能（E）和单位质量总焓（H）分别为完全气体的压力和温度状态方程为采用中心格式有限体积法求解积分形式的N-S方程，认为守恒流动变量在网格单元中心，时间推进方法采用显式Runge-Kutta混合多步法和隐式LU-SGS方法，非定常计算采用双时间法。采用预处理方法求解低马赫数下的收敛问题。湍流模型采用SST二方程模型。

　　扑翼计算网格基于嵌套网格方法，该方法不要求各子域共享边界，可有效减轻各区自身网格生成的难度，避免了计算网格扭曲过大、网格疏密分布不合理的状况，非常适合复杂外形或者有相对运动的情况。各个重叠子域的流场解在内边界要求必须光滑，通常由插值来实现，从而得到整个流场的解。洞点识别和插值单元寻址是重叠网格的两个关键技术。本文采用“距离减缩法”高效确定不同网格之间的对应单元，进而快速进行洞点识别，该方法效率高、易于实现，应用于插值单元寻址时，能够快速锁定插值贡献单元。本文所用的数值方法，已经过验证，其过程详见。

　　3结果分析作为基本规律研究，采用半展弦比为4的直机翼，翼型为NACA0008.考虑到仿鸟的研究范畴，取马赫数（Ma）为0. 03，雷诺数（D为7X104，迎角（吻）为5°，扑动幅度为20°，减缩频率（々）为0.2，其余参数详见各小节。

　　3.1升力和阻力特性计算无掠动情况下的气动特性，扭转角变化范围为0°25°，间隔3°。扭转角变化对升力系数和阻力系数的影响如所示。

　　扭转角对升力和阻力的影响规律Fig.7Theliftanddragcoefficientsvs.从可以看出：当扭转角在18范围内时，随着扭转角的增大，升力系数和推力系数均增加（即阻力系数减小，因为推力系数和阻力系数为相反数）当扭转角大于18°时，推力减小。18°是个较大的数值，实际中的扭转角通常不会超过18°因此，增加扑动时的扭转角可有效增加升力和推力，通过增加扭转角可改变飞行状态，例如改变飞行速度或高度等。

　　计算“0”字形扑动运动方式下的气动特性。掠动角度变化范围为0°20°间隔2°。“0”字形扑动方式下掠动角度对升力和阻力特性的影响规律如所示，可以看出：当掠动角度很小时（小于4），增加掠动角会减小升力，但却增加推力；当掠动角度在中等范围时（46），增加掠动角会使升力和推力同时增大；当掠动角度继续增加时（大于6），增加掠动角度会增加升力，但却减小推力。

　　计算“8”字形扑动运动方式下的气动特性。掠动角度变化范围为0°20°间隔2°。“8”字形扑动方式下掠动角度对升力和阻力特性的影响规律如所示。

　　从可以看出：8字形扑动规律下，随着掠动角的增加，升力增加但推力减小；当掠动角度在较小范围变化时（小于7），升力的增加比例大于推力的减小比例；当掠动角度在较大范围变化时（大于7），升力的增加比例小于推力的减小比例。

　　力为代价。

　　3.2流场分析展向50截面处的流场压力云图，如表1所示。

　　表1不同扑动方式流场压力变化云图扑动过程纯扑动扑动加扭转“”字形扑动“8”字形扑动下扑过程上扑过程从表1可以看出：纯扑动“产生的压力波动范围大，表明‘’纯扑动”产生的气动力波动范围也大；在“纯扑动”过程中，有相当部分气动力在上扑和下扑的过程中被抵消，造成运动的浪费；扑动加扭转“的压力波动情况优于”纯扑动“，其压力变化范围减小，且后缘压力增大，有利于增加推力：”“字形扑动方式在下扑过程中，由于有向前的掠动，增强了升力，但同时也增加了阻力：”8“字形扑动方式由于掠动时速度更快，对升力有较强的增加作用，但同时也显著增加了阻力。

　　4结论要增加升力，“8”字形扑动方式比“0”字形扑动方式在增强升力方面效果更好。

　　从运动机理上讲，掠动运动其本质是在定程度上改变了扑翼的合速度方向，从而影响了有效攻角，导致不同的气动特性。

　　增加扑动时的扭转角可有效增加升力和推力，从而改变飞行状态；通过利用掠动角也可达到增升的目的。

　　本文的研究结果可为复杂运动规律的微型扑翼飞行器设计提供了。