09-基于傅立叶法的旋翼尾迹悬停状态稳定性分析-吕维梁-7.docVIP

第二十六届（2010）全国直升机年会论文 基于傅立叶方法的旋翼尾迹悬停状态稳定性分析 吕维梁徐国华 （） 摘 要： 关键词：傅立叶方法；尾迹；悬停；稳定性；特征值分析 直升机的尾迹涡线与固定翼飞行器不同，呈螺旋状下降，并以其畸变的尾迹形状影响整个诱导速度场，从而对旋翼的气动载荷、噪声和性能产生重要的影响。尾迹涡线的运动过程受很多因素的影响，包括旋翼、前方来流、自诱导、相互诱导等，研究尾迹的稳定性对于认识旋翼流场的物理本质和改进尾迹数值计算具有重要的意义。 旋翼尾迹涡线流场的不稳定现象包括局部和全局的不稳定[1]。局部不稳定现象限于流场一小范围空间区域内，主要与桨尖涡有关，其不稳定性在一定条件下会扩散到整个流场，出现全局不稳定现象[2]。在1972年，Landgrebe[3]就已观察到尾迹不稳定或非周期尾迹的现象。同一时期，Widnall[4]理论分析了单根无限长的理想螺旋状涡线某一点受到正弦位移扰动后的稳定性，表明螺旋涡线存在不稳定模态。随后，在1974年，Gupta[5]将Widnall的方法拓展用于多条无限长螺旋状涡线的稳定性分析，该研究考虑了涡线间的干扰，但没有考虑到尾迹的畸变及桨叶的影响。1999年，Caradonna[6]通过实验发现了尾迹涡线的“涡对”现象。在2000年，Jain[7]在对两片桨叶旋翼在轴向飞行状态下桨尖涡的干扰进行研究时指出，桨尖涡干扰所形成的“涡对”现象与两无粘涡环所产生的蛙跳运动类似，并非一种内在不稳定现象，引起了对尾迹涡线稳定性的争论。同在2000年，Bhagwat在文[1]中使用自由尾迹方法，同样是假设在涡线上存在正弦扰动，对悬停状态的畸变尾迹涡线进行了特征值分析，得出了与无限长理想螺旋尾迹类似的稳定性结论，再次印证了旋翼尾迹涡线内在不稳定的结论。在国内，2007年，王海、徐国华[8]在文[1]的基础上，进一步分析了前飞状态下的尾迹稳定性，认为前飞尾迹也存在不稳定性。 前述论文在尾迹稳定性的研究上做出了卓有成效的工作，但并未指明正弦扰动的理论来源及意义，为进一步完整稳定性分析方法的理论依据，本文从傅立叶方法出发，将正弦扰动深化为任意扰动下的傅立叶谐波，分析了在任意扰动下的尾迹涡线上任一点的稳定性，假定涡线上存在任意的扰动分布，并将其展开成为一组傅立叶级数，通过对尾迹运动的线化扰动方程左右端取变分从而得到任意扰动下的尾迹涡线扰动方程，再对其进行特征值分析从而得到瞬时状态下尾迹上任一点的稳定性结果。另外，前述论文停留在孤立分析规则螺旋或旋翼尾迹上某一点的瞬时稳定性，本文在深化为任意扰动的基础上，以H-34桨叶模型为例，详细计算分析了悬停状态下，尾迹条数（桨叶片数）、涡核分布以及尾迹畸变对尾迹随寿命角变化的尾迹整体瞬时稳定性的影响，从而得出一些重要结论。 2 方法 2.1 傅立叶稳定性分析方法 假设涡线上存在任意大小及分布的一系列扰动，其中，寿命角为零的起始端由于固连桨叶，扰动为零。则这组扰动可以分解为一组傅立叶级数，如下所示： (1) 其中，为尾迹寿命角，对应各阶傅立叶级数，为各阶谐波放大系数，考虑到，式同时表达了正弦和余弦，但其中只有实数部分代表扰动量。为当前谐波的波数，即整个尾迹所能包含的整波数量。由于涡线末端扰动量不为零，为了补齐扰动函数周期性，所以波数可以为小数。 随时间 (2) 其中，为各阶谐波的权重，的增长率系数，可见，当时，该阶谐波发散不稳定，当时，该阶谐波收敛稳定，当时，该阶谐波处于中性状态。 根据涡的运动定理，旋翼尾迹涡线的运动方程为： (3) 式中，为为该点的当地速度，包括该点的总诱导速度和来流速度，通过该方程可以迭代求解出尾迹的收敛几何形状。 ，则扰动后的尾迹涡线的运动方程为： (4)式减去(3)式，得到扰动量的运动方程为： (5) 将上式右端用泰勒级数展开并略去高阶小量，得到线化后的扰动方程： (6) 2.3 诱导速度扰动项 方程(6)右端的速度项包括该点的来流速度和总诱导速度，目前只考虑尾迹涡线的扰动所带来的影响，因此，小扰动只对诱导速度项产生作用。对尾迹涡线进行离散处理，推导涡线上任意一点由于尾迹扰动量所产生的诱导速度项的解析式。 由诱导速度公式，即毕奥-萨瓦定律，一小段环量强度为的对任意一点 (7) 式中，为诱导速度的方向向量，到涡元的距离，，。其中、分别为向量、与向量的夹角，，，，。 点在靠近涡线处会得出的不正确的无穷大诱导速度，在(7)式中引入涡核修正因子，其中，为涡核半径，得到： 对上式左右两端取变分，得到： (9) 式中，，展开各个变分项，得到： (10) 式中，，，，其中，，，。 式中，