文章快速检索  
  高级检索
共轴刚性旋翼悬停状态地面效应气动特性
卢丛玲, 祁浩天, 徐国华, 史勇杰     
南京航空航天大学 直升机旋翼动力学国家级重点实验室, 南京 210016
摘要: 为了研究地面效应下共轴刚性旋翼的气动特性,建立了一套基于非定常雷诺平均Navier-Stokes方程的气动干扰数值方法,采用运动嵌套网格模拟双旋翼的反转运动。地面采用无滑移边界条件,并对旋翼和地面附近的网格进行加密,以更好地捕捉旋翼的流场细节和尾迹特征。计算结果与Lynx尾桨试验结果进行对比,验证了所建立方法的有效性。对地面效应下共轴刚性旋翼的气动性能和流场进行分析,结果发现:相对于单独的上下旋翼而言,共轴旋翼地面效应下的拉力增益更大,这是由于上下旋翼桨叶表面的压强干扰受地面高压的影响而减弱;地面的干扰主要影响双旋翼尾迹的径向位置,对其轴向位置影响不大,上下旋翼尾迹在地面附近相互融合、分裂,形成复杂的桨尖涡尾迹;双旋翼在地效下的尾迹径向扩张半径比单旋翼大,这是由于双旋翼的径向射流速度更大;随着旋翼距地面高度的增加,双旋翼间的气动干扰强度逐渐恢复,因此下旋翼拉力增益的下降速度比上旋翼更大;共轴旋翼桨尖涡相对卷起高度和扩张半径均随离地高度增加而减小。
关键词: 共轴刚性旋翼     地面效应     气动干扰     旋翼尾迹     计算流体力学     嵌套网格    
Aerodynamic characteristics of hovering coaxial rigid rotor in ground effect
LU Congling, QI Haotian, XU Guohua, SHI Yongjie     
National Key Laboratory of Science and Technology on Rotorcraft Aeromechanics, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
Abstract: In order to analyze the aerodynamic characteristics of hovering coaxial rigid rotors in ground effect, a numerical method based on unsteady Reynolds averaged Navier-Stokes equations is established. Moving overset grid is applied to simulate the reversal motion of the coaxial rotors. The non-slip boundary condition is adopted on the ground, and the grid around the rotors and the ground are refined to better capture the details of the flow field and the wake of the rotors. The comparison between the calculations and the experimental results of Lynx tail rotor verifies the effectiveness of the established method. The analyses of aerodynamic characteristics and the flow field of coaxial rigid rotors in ground effect show that the thrust increment of coaxial rotors in ground effect is larger than that of single upper and lower rotors. The increment is caused by the weakening of the interference between the surface of upper and lower blades due to high ground pressure. The ground interference mainly affects the radial position of the wake of the coaxial rotors, but it has little effect on its axial position. Upper and lower rotor wakes fuse and split with each other near the ground, forming complex tip vortex. The radius of the coaxial rotors wake in ground effect is larger than that of single rotor, because the radial jet velocity of the coaxial rotors is larger. With the increase of the height of the rotor from the ground, the aerodynamic interference intensity between the coaxial rotors gradually recovers, making the descending speed of the thrust increment of the lower rotor larger than that of the upper rotor. The relative roll height and the expansion radius of the coaxial rotors tip vortex decrease with the increase of off-ground height.
Keywords: coaxial rigid rotor     ground effect     aerodynamic interaction     rotor wake     CFD     overset grid    

直升机在垂直起降和贴地飞行时,必然会经历有地面效应的飞行状态。此时,地面对旋翼的性能有显著影响,在相同功率条件下,旋翼近地时产生的拉力与远离地面时相比明显增大,这种现象称之为“地面效应”(Ground Effect)。相对于传统的单旋翼直升机而言,共轴刚性旋翼直升机的地面效应更加复杂。这是由于共轴刚性旋翼由两副转向相反的旋翼构成,其下旋翼大部分区域处于上旋翼的下洗流和尾迹涡中,双旋翼间的气动干扰以及地面都会对旋翼的气动性能产生重要影响。旋翼靠近地面时,地面的存在也会使共轴旋翼间的气动干扰更加复杂。因此,针对共轴刚性旋翼地面效应状态进行气动性能和流场计算,分析其干扰特性和尾迹特征是十分必要的。

关于地面干扰对旋翼气动性能的影响,在试验和理论分析方面已有诸多研究。在试验方面,Sheridan和Weisner[1]在风洞中测量了不同飞行状态模型直升机贴地飞行时气动力和力矩的变化。Curtiss[2]、Light[3]和Ganesh[4]等分别以烟流法、宽场流动显示法和粒子图像测速(PIV)技术对地面效应下单旋翼尾迹涡线的几何形状进行了测量。在理论研究方面,何承健和高正[5]使用涡流理论研究了地面效应状态下单旋翼地面涡的变化。Curtiss等[2, 6]基于自由尾迹方法采用镜像法求解地面效应下的旋翼前飞流场。计算流体力学(CFD)方法近年来逐步被应用于旋翼地面效应下的气动特性研究[7-16]。康宁和孙茂[16]提出了一个简化的单旋翼贴地流场数值模拟方法,采用动量源模型模拟旋翼对流场的作用。Bhattacharyya和Conlisk[7]采用结构的嵌套网格方法进行了旋翼/地面干扰流场的数值计算。Kutz等[12]通过求解雷诺平均Navier-Stokes(RANS)方程来对地面效应下的单旋翼流场进行求解并针对旋翼离地高度对尾迹空间位置的影响进行了研究;朱明勇等[14]建立了一套基于嵌套网格的流场求解方法,运用该方法计算了单旋翼直升机小前进比贴地飞行时的流场。

目前,国内外针对双旋翼地面效应的研究较为稀少。康宁和孙茂[17]采用动量源方法,通过求解RANS方程,得到前飞地面效应状态下横列式、纵列式和共轴双旋翼的流场形态。覃燕华等[18]用自由尾迹方法计算了地面效应下的共轴双旋翼流场及诱导速度。陆陶冶等[19]研究了共轴双旋翼的地面效应及参数影响。谭剑锋等[20]采用涡粒子方法,计算了纵列式旋翼地面效应下的气动性能和流场特征,并与试验结果进行了比较。目前,关于共轴双旋翼地面效应的研究多采用动量源方法和尾迹方法,这些方法无法模拟双旋翼地面效应干扰状态下的流场细节。针对共轴双旋翼地面效应下的气动干扰,仍有很多方面需要研究。

鉴于此,本文在国内外研究工作的基础上,建立了一套适用于共轴双旋翼地面效应的流场分析方法,通过与地面效应下的试验结果进行对比,验证了该方法的有效性。在此基础上研究了共轴刚性旋翼地面效应下的气动性能和流场特性;通过与单旋翼地面效应的对比,分析了地面对双旋翼间气动干扰的影响;研究了双旋翼拉力增益特性及尾迹的几何形状特征;分析了旋翼离地高度对共轴旋翼地面效应下气动性能和流场特征的影响。

1 计算方法 1.1 共轴旋翼网格系统

本文采用运动嵌套网格方法对共轴旋翼地面效应的流场进行数值模拟。图 1为共轴旋翼地面效应网格系统示意图。计算中对旋翼和地面附近的背景网格进行加密处理,以便更好地捕捉旋翼尾迹。图 2为双旋翼对齐时嵌套网格系统展向截面示意图。计算中采用的笛卡尔背景网格节点数为264×320×264,分别对应XYZ方向。同时,对桨叶贴体网格奖尖区域进行加密,每片桨叶网格节点数为215×94×124,分别对应桨叶周向、法向和展向。采用的嵌套网格总网格量约为3.2千万。

图 1 共轴刚性旋翼嵌套网格系统示意图 Fig. 1 Sketch map of overset grid system for coaxial rigid rotor
图 2 共轴刚性旋翼网格展向截面 Fig. 2 Spanwise cross-section of coaxial rigid rotor grid
1.2 流场求解方法 1.2.1 控制方程及求解方法

采用惯性坐标系下的雷诺平均Navier-Stokes方程进行流场求解,其积分形式的表达式为

$\frac{\partial}{\partial t} \iiint\limits_{\partial V} \boldsymbol{W} \mathrm{d} V+{\mathop{{\int\!\!\!\!\!\int}\mkern-21mu \bigcirc}\limits_{\partial S} {\left( {{\mathit{\boldsymbol{F}}_{\rm{c}}} - {\mathit{\boldsymbol{F}}_{\rm{v}}}} \right) \cdot \mathit{\boldsymbol{n}}{\rm{d}}\mathit{S}} = {\bf{0}}} $ (1)
 

式中:SV分别为控制体表面积和体积;W为守恒变量;FcFv分别为对流通量和黏性通量,其表达式分别为

$\boldsymbol{W}=\left[\begin{array}{lllll}{\rho} & {\rho u} & {\rho v} & {\rho w} & {\rho E}\end{array}\right]^{\mathrm{T}} $
$\boldsymbol{F}_{\mathrm{c}}=\left[\begin{array}{c}{\rho \boldsymbol{V}_{\mathrm{r}}} \\ {\rho u \boldsymbol{V}_{\mathrm{r}}+\boldsymbol{n}_{x} p} \\ {\rho v \boldsymbol{V}_{\mathrm{r}}+\boldsymbol{n}_{y} p} \\ {\rho w \boldsymbol{V}_{\mathrm{r}}+\boldsymbol{n}_{z} p} \\ {\rho H \boldsymbol{V}_{\mathrm{r}}}\end{array}\right] $
$\boldsymbol{F}_{\mathrm{v}}=\left[\begin{array}{c}{0} \\ {\boldsymbol{n}_{x} \tau_{x x}+\boldsymbol{n}_{y} \tau_{x y}+\boldsymbol{n}_{z} \tau_{x z}} \\ {\boldsymbol{n}_{x} \tau_{y x}+\boldsymbol{n}_{y} \tau_{y y}+\boldsymbol{n}_{z} \tau_{y z}} \\ {\boldsymbol{n}_{x} \tau_{z x}+\boldsymbol{n}_{y} \tau_{z y}+\boldsymbol{n}_{z} \tau_{z z}} \\ {\boldsymbol{n}_{x} \Phi_{x}+\boldsymbol{n}_{y} \Phi_{y}+\boldsymbol{n}_{z} \Phi_{z}}\end{array}\right] $

其中:uvw表示3个方向的速度分量;ρp分别表示流体密度和压强;EH为单位质量气体的总能和总焓;Vr为相对速度;n=[nx  ny  nz]T为单元表面法矢量;τijΦi(i=x, y, z; j=x, y, z)为黏性相关项。

为减小非物理耗散的影响,提高流场计算精度,采用二阶MUSCL(Monotone Upstream-centered Schemes for Conservation Laws)格式对单元内流场进行重构,并采用逆风格式(Roe格式)[21]计算网格交接面上的对流通量。由于地面效应下双旋翼干扰流场是非定常的,需要采用双时间法进行时间推进,在伪时间方向上使用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式[22]进行时间推进以提高流场的计算效率。湍流模型采用Spalart-Allmaras一方程模型[23]

1.2.2 边界条件

本文在使用CFD方法计算共轴旋翼地面效应下的流场时,地面采用无滑移、不可穿透的物面边界条件,远场边界条件采用一维Riemann不变量进行设置。在无地效状态下,旋翼流场底面边界与旋翼的距离增大,且与远场边界设置相同。

2 算例验证

图 3图 4分别为Lynx尾桨有地面效应(In Ground Effect, IGE)下桨尖涡等涡量图和桨尖涡空间位置分布。Lynx尾桨由4片无扭转的矩形桨叶构成,其半径R为1.105 m、弦长c为0.18 m、翼型为NPL9615。尾桨在悬停状态下,旋翼距离地面的高度h/R=0.52、总距为15°、桨尖马赫数为0.56。为了研究网格尺寸对桨尖涡计算结果的影响,本文生成了3套不同尺寸的嵌套网格来计算地效状态的尾桨桨尖涡。由于旋翼桨尖区域和尾迹区域的网格尺寸对其尾迹计算影响较大,因此在这两个区域进行加密。网格1的桨尖部分网格尺寸为0.08c,尾迹区域的背景网格尺寸为0.15c;网格2对应的尺寸分别为0.05c和0.08c;网格3为0.025c和0.05c。由图 3可见,网格尺寸对桨尖涡空间位置分布影响较小,主要影响桨尖涡的涡强及耗散情况,更密的网格能够减小桨尖涡的耗散,捕捉到更长的尾迹。根据本文对桨尖涡的计算需求,选择网格2的尺寸来对后续的地面效应下共轴旋翼流场进行计算。如图 4所示,其中r为桨尖涡半径,桨尖涡位置的计算值在100°涡龄角前与试验吻合良好,随着涡龄角增大,桨尖涡空间位置分散程度增加,这表明尾迹在地面附近的非定常性增加,这可能是由于地面阻塞改变流场中气流的流动方向,进而影响桨尖涡的空间位置;且地面的黏性作用会增大桨尖涡的耗散,使得桨尖在地面附近拉伸扩展[20]。因此,地面附近桨尖涡的空间位置分布更加分散。

图 3 Lynx尾桨地面效应下桨尖涡等涡量图 Fig. 3 Blade tip vortex of Lynx tail rotor IGE
图 4 Lynx尾桨地面效应下桨尖涡轨迹 Fig. 4 Trajectories of tip vortex of Lynx tail rotor IGE

图 5给出了Lynx尾桨[3]拉力增益随离地高度的变化曲线。对比试验值可见,计算值与其基本吻合,说明本文所建立的旋翼地面效应计算方法是可靠的。同时由图可见,当离地高度h/R>2时,地面对其拉力基本没有影响。

图 5 Lynx尾桨拉力增益随地面高度的变化 Fig. 5 Variation of Lynx tail rotor thrust increments with heights
3 共轴旋翼地面效应计算与分析 3.1 拉力增益与流场

本文以Harrington-2共轴旋翼为计算模型,计算了其在有/无地面效应状态下的悬停干扰流场和气动性能。该旋翼参数如表 1所示,其采用对称翼型,厚度沿展向呈线性分布。计算中采用本文课题组[24]所建立的共轴双旋翼配平方法对无地面效应下的双旋翼进行扭矩配平,得到上旋翼总距为9.71°,下旋翼为10.41°。地面效应状态下,上旋翼离地面高度为0.5R,在满足上下旋翼扭矩平衡、双旋翼总功率与无地效状态相同两个条件下,对共轴旋翼进行配平得到上旋翼总距为9.63°,下旋翼为10.12°。

表 1 Harrington-2共轴双旋翼参数 Table 1 Parameters of Harrington-2 coaxial rotor
参数 数值
旋翼半径/m 3.81
旋翼间距/m 0.609 6
桨叶弦长/m 0.457 2
桨叶根切 0.2R
桨叶片数 2+2

图 6为共轴旋翼在有地效(In Ground Effect, IGE)和无地效(Out Ground Effect, OGE)两种状态下旋翼瞬时拉力系数CT对比。可见,在有/无地效两种状态,上下旋翼瞬时拉力均呈现明显的周期性波动。这是因为双旋翼旋转一周,上下桨叶相遇4次,拉力也随之产生4次剧烈波动。在无地效状态,由于下旋翼一直处于上旋翼的尾流之中,有效迎角减小,所以虽然下旋翼总距比上旋翼大,其拉力却较小。在地效状态,上下旋翼拉力除了周期性波动外,还有地面干扰引起的小幅度非定常波动。而且下旋翼的平均拉力略高于上旋翼,这是因为下旋翼更接近地面,拉力增益更为明显。

图 6 共轴旋翼有/无地效悬停状态旋翼瞬时拉力 Fig. 6 Temporal variations of coaxial rotor thrust at IGE/OGE hovering status

为了捕捉共轴旋翼桨盘内的瞬时轴向诱导速度Vy,在桨叶展向0.4R~0.9R的位置设置监测点,监测点的位置如图 7所示。其中,采用桨尖速度Vtip对其进行无量纲化。在0°和180°方位角附近,由于桨叶穿过监测点,因此产生诱导速度断点。如图 8所示,诱导速度随着展向位置的增大而增加;且上旋翼的诱导速度分布更为均匀,而受到上旋翼尾流干扰,下旋翼的诱导速度随展向位置变化更为剧烈。在无地面干扰时,下旋翼的桨盘诱导速度受到上旋翼下洗流的作用,其桨叶展向外侧(0.6R~0.9R)诱导速度明显大于上旋翼。地面的干扰使得上下旋翼诱导速度整体降低,桨叶的有效迎角因诱导速度的减小而增加,引起旋翼拉力增加。且地面的干扰主要降低桨叶0.4R~0.8R处的诱导速度,而桨尖部分(0.9R)却有所增加,这主要是由地面卷起的尾迹涡对桨尖的干扰所引起的。

图 7 诱导速度监测点位置示意图 Fig. 7 Sketch map of induced velocity monitoring points in different positions
图 8 共轴旋翼有/无地效悬停状态桨盘诱导速度 Fig. 8 Disc induced velocity of coaxial rotor at IGE/OGE hovering status

旋翼拉力增益现象是直升机旋翼贴地飞行时出现的重要特征。在功率相同的条件下,旋翼有地效状态下的拉力与其在无地效状态下的拉力比值即为拉力增益。表 2给出了悬停时有/无地效状态下,单旋翼与共轴旋翼的拉力系数对比。其中,单独的上(下)旋翼是指去除共轴旋翼的下(上)旋翼后得到的旋翼。单独的上下旋翼地效下的拉力增益分别为1.309和1.318。共轴旋翼的上旋翼拉力增益为1.390,下旋翼拉力增益为1.494。与单旋翼相比,共轴旋翼地效下的拉力增益更大,且下旋翼的拉力增益增幅更加明显。无地效状态下,由于上下旋翼间存在气动干扰,因此双旋翼的拉力均低于单旋翼,且受上旋翼下洗流的影响,双旋翼间的气动干扰对下旋翼拉力影响更大。地效状态下,地面的干扰使得双旋翼间的气动干扰减弱,因此双旋翼拉力增益比单旋翼更大,且下旋翼拉力增益比上旋翼更明显。

表 2 有/无地效状态旋翼拉力系数 Table 2 CT of rotor at IGE/OGE status
旋翼 拉力系数 拉力增益
无地效 有地效
单独上旋翼 0.009 4 0.012 3 1.309
单独下旋翼 0.010 7 0.014 1 1.318
共轴上旋翼 0.008 2 0.011 4 1.390
共轴下旋翼 0.007 7 0.011 5 1.494

为了进一步研究地面效应对双旋翼气动干扰的影响,图 9给出了上下旋翼桨叶相遇时刻的展向剖面压力系数Cp分布。无地效状态下,由于下旋翼桨叶上方的低压与上旋翼桨叶下方的高压相互影响,使得两者的桨叶上下表面压强差均减小,因此拉力减小。地效状态下,由于地面的阻塞及黏性作用,地面附近气流速度减小,压强增大,桨叶下表面的压强受地面产生的高压影响而增加,因此桨叶上下表面压强差增大,拉力增加。由于下旋翼距地面更近,因此其拉力增益比上旋翼更大。如图 9所示,受地面高压影响,双旋翼桨叶间的压强干扰也随之减弱,可见地面的存在使上下旋翼间的干扰减弱,因此地面效应下的旋翼拉力增益相对于单旋翼更明显。

图 9 有/无地效状态单/共轴旋翼展向剖面压力系数分布 Fig. 9 Pressure coefficient distribution of spanwise cross-sections of single rotor/coaxial rotor at IGE/OGE status

图 10为共轴旋翼有/无地效两种状态桨叶表面压力系数Cp分布。其中,双旋翼位置与图 8相同。对比可见,地效状态下共轴旋翼的上旋翼桨叶下表面压强较无地效有所增大,而下旋翼的上表面压强则比无地效状态降低。这也表明了上旋翼下表面的高压与下旋翼上表面的低压间的干扰减弱,上下桨叶表面压强均向无干扰状态靠近,这与图 9所述结论一致,表明地面效应减弱了共轴旋翼间的气动干扰强度。

图 10 有/无地效状态共轴旋翼桨叶表面压力系数分布 Fig. 10 Surface pressure coefficient distribution of coaxial rotor at IGE/OGE status
3.2 地面效应尾迹特征

图 11对比了有/无地效下的双旋翼桨尖涡径向位置和轴向位置。在无地效状态,上下旋翼尾迹随着涡龄角增加均逐渐收缩。上旋翼桨尖涡半径收缩速度大于下旋翼,在180°涡龄角位置与下旋翼桨叶发生干扰,桨尖涡半径发生剧烈变化。地效状态,上下旋翼桨尖涡半径均是先收缩后扩张,上旋翼桨尖涡在140°涡龄角位置收缩至最小半径,140°~170°涡龄角间半径逐渐扩张,由于其在180°位置处会与下旋翼桨叶进行干扰,因此在170°~200°涡龄角间桨尖涡半径产生小幅度抖动,200°涡龄角后,其半径继续收缩。由于下旋翼离地面更近,其半径收缩速度小于上旋翼,在130°涡龄角位置收缩至最小值。在上旋翼220°、下旋翼140°涡龄角附近,由于上下旋翼桨尖涡靠近、相互诱导,上下旋翼尾迹半径迅速扩张、变化剧烈。

图 11 共轴刚性旋翼桨尖涡轨迹 Fig. 11 Trajectories of blade tip vortex of coaxial rigid rotor

对比两种状态下的尾迹轴向位置分布可见,在无地效状态下,上旋翼尾迹下降速度较为稳定,在经过下旋翼桨盘时(180°涡龄角)被加速,下降速度增加。下旋翼在110°~150°涡龄角位置,受到附近的上旋翼桨尖涡诱导,其桨尖涡轴向位置略有上升。在地效状态下,上旋翼桨尖涡与无地效状态基本一致,由于下旋翼桨盘的诱导速度减小,因此其加速作用减弱,桨尖涡下降速度相对无地效状态减小。由于在上旋翼220°、下旋翼140°涡龄角附近,上下旋翼桨尖涡发生融合,因此二者的桨尖涡在该位置的轴向高度均基本不变。可见,地面的干扰作用对于双旋翼桨尖涡的半径影响显著,而对其轴向位置影响较小。地面效应下,双旋翼尾迹的收缩速度明显减小,且上、下旋翼分别在140°和130°涡龄角位置半径收缩到最小。

图 12为共轴旋翼地面效应下桨尖涡示意图,由图可见,上下旋翼桨尖涡在地面附近相互融合、分裂,形成复杂的桨尖涡尾迹。随着旋翼的旋转,桨尖涡从桨尖脱落螺旋向下发展。由图 13可见,无地效状态下桨尖涡在下移过程中沿径向逐渐收缩。在有地效状态下,桨尖涡在向下发展的过程中受到地面阻碍,其径向位置先收缩后扩张,桨尖涡靠近地面后向上卷起,其径向位置继续扩张。单旋翼的尾迹涡从地面卷起后的扩张半径约为1.65R,最大卷起高度约为0.28R;双旋翼扩张半径约为2.2R,最大卷起高度约为0.24R。可见,单旋翼与共轴旋翼的桨尖涡卷起高度差别较小,但共轴旋翼桨尖涡的扩张半径远大于单旋翼。旋翼下洗流由于受到地面的阻挡,形成向外的径向射流,空气经过共轴旋翼的两次加速,其下洗流速度大于单旋翼,因此径向射流速度也随之较大。图 14对比了不同径向位置处,单旋翼和共轴刚性旋翼在有/无地效悬停状态的径向速度Vr。可以明显看出,地效下的旋翼径向速度均大于无地效状态的径向速度。在0.8R径向位置处,无地效下的旋翼径向速度为负,表明桨尖涡处于收缩状态。在1.2R处,共轴旋翼地面附近的径向速度约为单旋翼的两倍,其桨尖涡扩张半径也较大。

图 12 共轴旋翼地面效应状态桨尖涡示意图 Fig. 12 Blade tip vortex schematic of coaxial rotor at IGE status
图 13 有/无地效状态桨尖涡展向截面位置 Fig. 13 Spanwise cross-sections station of blade tip vortex at IGE/OGE status
图 14 有/无地效状态旋翼径向速度 Fig. 14 Radial velocity of rotor at IGE/OGE status
3.3 旋翼离地高度影响

图 15为共轴旋翼平均拉力系数随高度的变化曲线,本文所述共轴旋翼离地高度为其上旋翼离地高度。如图所示,随着旋翼距地面高度增加,双旋翼拉力均减小,且下旋翼拉力下降速度更大。这主要是由于随着地面效应的减弱,地面对双旋翼间的气动干扰的影响减小,双旋翼间的气动干扰强度逐渐恢复。因此,下旋翼在地面效应减弱的同时,受到上旋翼的干扰增大,拉力较上旋翼减小更多。

图 15 共轴旋翼拉力随高度的变化 Fig. 15 Variation of coaxial rotor thrust with heights

图 16给出了共轴旋翼在不同离地高度下瞬时拉力变化。共轴旋翼拉力在每个相遇方位角附近(93°)有脉冲式拉力波动,这是由于桨叶的“厚度效应”使得相遇时刻桨叶间的气流速度增大、压强减小,因此上旋翼拉力突降而下旋翼拉力激增[25]。由图可见,随着旋翼离地高度的减小,“厚度效应”引起的上下旋翼的脉冲拉力波动也随之减弱,这与图 9所述的共轴旋翼间气动干扰减弱的结论一致。

图 16 不同离地高度下共轴旋翼有/无地效状态的瞬时拉力 Fig. 16 Temporal variations of coaxial rotor thrust with different heights at IGE/OGE status

图 17图 18分别为不同离地高度下,共轴旋翼展向截面的尾迹和诱导速度图。图 17中:Vor表示涡强度(Vorticity Magnitude),由图可见,随着旋翼离地高度h/R从1.5降至0.5,旋翼地面效应逐渐增强,桨根涡强度增加明显,桨尖涡的扩张半径也随之从1.54R增加至2.2R。旋翼尾迹在耗散前最终的卷起高度也随旋翼离地高度的减小而降低。如图 18所示,随着旋翼离地高度的减小,流场中的轴向诱导速度逐渐减小,而径向诱导速度逐渐增大。地面效应下桨尖涡的扩张半径随径向速度的增加而增加。随着离地高度的增加,共轴旋翼桨尖涡的最终卷起高度也随之增加。但以离地高度为基准,桨尖涡的相对卷起高度却在逐渐降低,旋翼离地高度从0.5R增至1.5R,桨尖涡相对卷起高度从0.48减少至0.35。这与流场中轴向诱导速度随离地高度增加而增加的分布特征一致。

图 17 不同离地高度共轴旋翼有地面效应下展向截面尾迹 Fig. 17 Blade tip vortex of spanwise cross-sections of coaxial rotor with different heights at IGE status
图 18 不同离地高度共轴旋翼有地面效应下展向截面诱导速度 Fig. 18 Spanwise cross-sections induced velocity of coaxial rotor with different heights at IGE status
4 结论

本文基于CFD方法和运动嵌套网格技术,研究了共轴刚性旋翼地面效应状态下的气动性能和流场特征。分析了地面对共轴旋翼间的气动干扰以及共轴旋翼尾迹的影响,在此基础上研究了旋翼离地高度对共轴旋翼地面效应的影响,得到以下结论:

1) 由于地面的阻塞及黏性作用,地面附近气流速度减小、压强增大;地面干扰引起的高压对双旋翼的表面压强产生影响,使双旋翼间的气动干扰减弱。因此在地面效应状态下,共轴旋翼的上下旋翼的拉力增益均大于单独的上旋翼或下旋翼,且其下旋翼拉力增益增加幅值比上旋翼更大。

2) 共轴旋翼的上下旋翼桨尖涡在地面附近会相互融合、分裂,形成复杂的桨尖涡尾迹。地面的干扰作用对于双旋翼桨尖涡的径向位置影响较为显著,而对其轴向位置影响较小。地面效应状态下,双旋翼尾迹的收缩速度明显减小,在本文计算状态上、下旋翼尾迹分别在140°和130°涡龄角位置半径收缩到最小值。

3) 双旋翼的桨尖涡在地面附近的扩张半径远大于单旋翼,桨尖涡卷起高度与单旋翼差别不大。这是由于流场中向下的诱导速度受到地面的阻挡转变成沿径向向外的射流速度。与单旋翼相比,双旋翼的诱导速度较大,因此地面附近的径向速度随之增大,导致其桨尖涡扩张半径也较大。

4) 随着离地高度的增加,地面对双旋翼间的气动干扰的影响减小,双旋翼间的气动干扰强度逐渐恢复。双旋翼拉力增益随之减小,且下旋翼拉力增益减小更快。旋翼距地面高度从0.5R增至1.5R,桨尖涡扩张半径从2.2R减小至1.54R,相对卷起高度从0.48减少至0.35。

参考文献
[1] SHERIDAN P F, WEISNER W. Aerodynamics of helicopter flight near ground[C]//Proceedings of 33rd Annual Forum of American Helicopter Society, 1977.
[2] CURTISS H C, SUN M, PUTNAM W F, et al. Rotor aerodynamics in ground effect at low advance ratios[J]. Journal of the American Helicopter Society, 1984, 29(1): 48-55.
Click to display the text
[3] LIGHT J S. Tip vortex geometry of a hovering helicopter rotor in ground effect[J]. Journal of the American Helicopter Society, 1993, 38(2): 34-42.
Click to display the text
[4] GANESH B, KOMERATH N. Study of ground vortex structure of rotorcraft in ground effect at low advance ratios[C]//24th AIAA Applied Aerodynamics Conference. Reston, VA: AIAA, 2006.
[5] 何承健, 高正. 贴地飞行的旋翼尾迹研究[J]. 航空学报, 1986, 7(4): 325-331.
HE C J, GAO Z. A study of the rotor wake in nap-of-the-earth[J]. Acta Aeronautica et Astronautica Sinica, 1986, 7(4): 325-331. (in Chinese)
Cited By in Cnki (6) | Click to display the text
[6] 孙茂, CURTISS H C. 近地低速飞行时旋翼尾涡系的畸变及其诱导速度[J]. 空气动力学学报, 1989, 7(1): 35-42.
SUN M, CURTISS H C. Rotor vortex wake distortion and its induced velocity in ground effect at low speed[J]. Acta Aerodynamica Sinica, 1989, 7(1): 35-42. (in Chinese)
Cited By in Cnki (3) | Click to display the text
[7] BHATTACHARYYA S, CONLISK A T. The structure of the rotor wake in ground effect[C]//41th Aerospace Sciences Meeting and Exhibit. Reston, VA: AIAA, 2003.
[8] HWANG J Y, CHOI J H, KWON O J. Assessment of S-76 rotor hover performance in ground effect using an unstructured mixed mesh method[C]//54th AIAA Aerospace Sciences Meeting. Reston, VA: AIAA, 2016.
[9] LAKSHMINARAYAN V K, KALRA T S, BAEDER J D. Detailed computational investigation of a hovering microscale rotor in ground effect[J]. AIAA Journal, 2013, 51(4): 893-909.
Click to display the text
[10] CROZON C, STEIJL R, BARAKOS G N. Numerical study of helicopter rotors in a ship airwake[J]. Journal of Aircraft, 2014, 51(6): 1813-1832.
Click to display the text
[11] SUGIURA M, TANABE Y, SUGAWARA H, et al. Numerical simulations and measurements of the helicopter wake in ground effect[J]. Journal of Aircraft, 2017, 54(1): 209-219.
Click to display the text
[12] KUTZ B M, BENSING F, KEßLER M, et al. CFD calculation of a helicopter rotor hovering in ground effect[C]//DILLMANN A, HELLER G, KREPLIN H P, et al. New Results in Numerical and Experimental Fluid Mechanics Ⅷ: Contributions to the 17th STAB/DGLR Symposium Berlin, Germany 2010. Berlin Heidelberg: Springer, 2013: 297-304.
[13] 叶靓, 招启军, 徐国华. 基于非结构嵌套网格方法的旋翼地面效应数值模拟[J]. 航空学报, 2009, 30(5): 780-786.
YE L, ZHAO Q J, XU G H. Numerical simulation on flowfield of rotor in ground effect based on unstructured embedded grid method[J]. Acta Aeronautica et Astronautica Sinica, 2009, 30(5): 780-786. (in Chinese)
Cited By in Cnki (10) | Click to display the text
[14] 朱明勇, 招启军, 王博. 基于CFD和混合配平算法的直升机旋翼地面效应模拟[J]. 航空学报, 2016, 37(8): 2539-2551.
ZHU M Y, ZHAO Q J, WANG B. Simulation of helicopter rotor in ground effect based on CFD method and hybrid trim algorithm[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(8): 2539-2551. (in Chinese)
Cited By in Cnki (2) | Click to display the text
[15] 康宁, 孙茂. 旋翼近地飞行时诱导速度的N-S方程计算[J]. 空气动力学学报, 1998, 16(2): 221-225.
KANG N, SUN M. Navier-Stokes calculations of induced velocity of a rotor in forward flight with ground effects[J]. Acta Aerodynamica Sinica, 1998, 16(2): 221-225. (in Chinese)
Cited By in Cnki (14) | Click to display the text
[16] 康宁, 孙茂. 旋翼尾流与地面干扰时地面涡现象的研究[J]. 力学学报, 1998, 30(5): 615-620.
KANG N, SUN M. Investigation of the ground vortex phenomenon due to the interaction between rotor's wake and the ground[J]. Chinese Journal of Theoretical and Applied Mechanics, 1998, 30(5): 615-620. (in Chinese)
Cited By in Cnki (1) | Click to display the text
[17] KANG N, SUN M. Simulated flowfields in near-ground operation of single-and twin-rotor configurations[J]. Journal of Aircraft, 2000, 37(2): 214-220.
Click to display the text
[18] 覃燕华, 朱清华, 邵松. 共轴双旋翼悬停地面效应气动特性分析[J]. 南京航空航天大学学报, 2015, 47(2): 266-274.
QIN Y H, ZHU Q H, SHAO S. Aerodynamic characteristics analysis for hovering coaxial rotors in ground effect[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2015, 47(2): 266-274. (in Chinese)
Cited By in Cnki | Click to display the text
[19] 陆陶冶, 陈仁良, 吉洪蕾, 等. 共轴双旋翼悬停地面效应分析[J]. 哈尔滨工业大学学报, 2017, 49(10): 45-52.
LU T Y, CHEN R L, JI H L, et al. Performance analysis of coaxial-rotorfor hovering in ground effect[J]. Journal of Harbin Institute of Technology, 2017, 49(10): 45-52. (in Chinese)
Cited By in Cnki | Click to display the text
[20] TAN J F, SUN Y M, BARAKOS G N. Vortex approach for downwash and outwash of tandem rotors in ground effect[J]. Journal of Aircraft, 2018, 55(6): 2491-2509.
Click to display the text
[21] ROE P L. Approximate Riemann solvers, parameter vectors, and difference schemes[J]. Journal of Computational Physics, 1981, 43: 357-372.
Click to display the text
[22] LUO H, BAUM J D, LÃHNER R. A fast, matrix-free implicit method for computing low Mach number flows on unstructured grids[J]. International Journal of Computational Fluid Dynamics, 2000, 14(2): 133-157.
Click to display the text
[23] SPALART P, ALLMARAS S. A one-equation turbulence model for aerodynamic flows[C]//30th Aerospace Sciences Meeting and Exhibit. Reston, VA: AIAA, 1992.
[24] 祁浩天, 史勇杰, 徐国华, 等. 共轴刚性旋翼气动干扰及操纵特性分析[J]. 航空动力学报, 2017, 32(12): 3004-3012.
QI H T, SHI Y J, XU G H, et al. Analysis of aerodynamic interaction andtrim characteristics of rigid coaxial rotor[J]. Journal of Aerospace Power, 2017, 32(12): 3004-3012. (in Chinese)
Cited By in Cnki | Click to display the text
[25] LAKSHMINARAYAN V K, BAEDER J D. High-resolution computational investigation of trimmed coaxial rotor aerodynamics in hover[C]//63rd American Helicopter Society Annual Forum, 2007.
http://dx.doi.org/10.7527/S1000-6893.2019.23055
中国航空学会和北京航空航天大学主办。
0

文章信息

卢丛玲, 祁浩天, 徐国华, 史勇杰
LU Congling, QI Haotian, XU Guohua, SHI Yongjie
共轴刚性旋翼悬停状态地面效应气动特性
Aerodynamic characteristics of hovering coaxial rigid rotor in ground effect
航空学报, 2019, 40(12): 123055.
Acta Aeronautica et Astronautica Sinica, 2019, 40(12): 123055.
http://dx.doi.org/10.7527/S1000-6893.2019.23055

文章历史

收稿日期: 2019-04-03
退修日期: 2019-04-12
录用日期: 2019-04-29
网络出版时间: 2019-05-09 10:37

相关文章

工作空间