23-悬停状态直升机桨叶的气动优化设计—刘国强-6

190

第二十六届（2010）全国直升机年会论文

悬停状态下旋翼桨叶的气动外形优化设计

刘国强1 高卓飞2 唐正飞1

（1.南京航空航天大学直升机旋翼动力学重点实验室 南京 210016）

（2.中国人民解放军61267部队，北京通州 101114）

摘 要：本文在改进的动量-叶素理论基础上建立了遗传算法优化方法，并且验证了方法的有效性。针对旋翼桨叶气动外形设计的有关问题，介绍了桨叶气动外形优化设计方法，以旋翼的悬停效率为目标函数，设计变量主要包括桨叶的弦长、扭转角、翼型选择和配置，在满足约束条件下悬停性能有所提高，表明了优化方法的可行性。

关键词：桨叶气动优化；动量-叶素理论；悬停性能

1 引言

旋翼是直升机的主要气动面，直升机的升力、方向操纵都与旋翼有关，旋翼气动特性的好坏与旋翼桨叶的气动外形有着密切的关系，先进的桨叶外形对提高直升机的气动性能，延缓失速、降低功率损失、降低振动和噪声都有明显的作用，所以桨叶的气动外形设计以成为直升机研制的关键技术之一。

由于桨叶气动环境的复杂性，涉及设计变量之多，以及设计变量之间的相互影响关系，固定一个变量而调整其他变量的设计方法已经不再适用，1981年，Bingham, Gene J 等人对桨叶的根梢比、扭转角、翼型分布等变化对悬停和前飞性能的影响做了详细分析，其方法就是运用上述改变单一变量法，使得计算量比较大。1991年，Joanne L. Walsh 采用CAMRAD 计算程序，计算加入了尾流模型，通过CONMIN 优化程序在满足约束的情况下计算目标函数的最值，建立了一种桨叶优化设计方法。

本文采用改进的动量叶素理论求解旋翼的气动特性，优化方法采用遗传算法的优化，以悬停效率为目标函数，桨叶沿展向不同位置的弦长和扭转角以及不同种翼型位置为设计变量，加入约束，在满足条件的情况下使得桨叶性能有所提高，达到了预期设计目的。

2 改进的动量叶素理论气动特性计算

首先采用“微观”的动量理论求解旋翼的拉力及其扭矩[5]

，对于微观动量理论，即把滑流分成很多同心基元滑流圆环处理，此时的旋翼拉力公式转化为：

101122(2)()d T d m v r d r V v v ρπ==⨯⋅+ （1） 加入尾流旋转对诱导速度的影响，定义b 为弦向入流因子即尾流旋转速度与弦向速度的比值。由角动量守恒得到桨叶微段扭矩为：

3

1()4t dM dm v r r v b dr πρ==Ω⋅ （2） 然而在叶素理论中，垂直飞行时，旋翼以角速度一方面以角速度Ω旋转，一方面以爬升速度0V 向上运动，以距离桨榖中心r 处的叶素为例，k 片桨叶在半径r 处产生的拉力和扭矩之和为：

21

(cos sin )2y x dT W kc C C dr ρφφ=

- （3） 2

1(sin cos )2

y x dM W kc C C rdr ρφφ=+ （4）

191

如果加入桨尖修正，考虑普朗特叶尖修正因子： 2

arccos f F e π=

，其中，2sin k R r

f r φ

-=-⋅。 （5） 令式（1）和式（3）相等即可求出旋翼在半径r 处的诱导速度，由于等式两端都与诱导速度有关，所以本次采用迭代的方法求解诱导速度，用σ/T C 配平悬停时所需的总距0θ，直至满足收敛条件，进而得到悬停效率3/2/2

FM C m =。如下（图1）所示，具体的流程为：

图1 旋翼桨叶气动力计算流程

通过采用上述方法，按照文献[1]给出的样例直升机数据，对桨叶的气动特性进行了计算，验证气动模型的正确性。模型的主要参数如下表：

通过编程计算，计算结果和给定的样例直升机对比如下表二：

3 遗传算法在桨叶优化中的应用

遗传算法是一种借鉴生物界自然选择和自然遗传机制的随机搜索算法。它与传统的算法不同，大多数古典的优化算法是基于一个单一的度量函数（评估函数）的梯度或者较高次统计，以产生一个确定性的试验解序列；遗传算法不依赖于梯度信息，而是通过模拟自然进化过程来搜索最优解,它利用某种编码技术，作用于称为染色体的数字串，模拟由这些串组成的群体的进化过程。由于遗传算法是模拟自然界的生物进化过程，与传统的优化算法相比较，遗传算法具有自组织、自适应和自学习性、搜索点是并行的，无需求导，这样可以解决很多复杂性、多极点的问题，具有极强的搜索能力和鲁棒性。

3.1旋翼桨叶优化问题描述

扭转角无论是在悬停或者是前飞都起着很大作用，好的扭转角可以提高悬停效率约为3%—5%，一般为了考虑制造上的难题将桨叶做成线性扭转的，原则上真正的最佳扭转是非线性的，因此我们将

192

桨叶的扭转设成非线性的，为了优化的灵活性，桨叶的预扭被定义为如下[2]

，沿展向变化如下：

()1(1)r r r δ

θθτ⎡⎤=⨯+⋅-⎣⎦ （6）

其中r θ是桨根处的扭角，δ为扭角外形参数，τ为桨根与桨尖扭角比例系数。弦长沿径向分布如下： []

1

()1(1)1p

r c r c r r ε

λ⎡⎤=⨯+-⋅-⎢⎥⎣

⎦

（7）

其中r c 为桨根处弦长，λ为根梢比，,p ε皆为弦长分布参数。

3.2 数学模型建立

对于一般的优化算法，其表述如下：

目标函数：max ()i F x , 1,2,...,i M =

约束条件：()0

()0j i j i i i i

g x h x a x b ≥⎧⎪=⎨⎪

≤≤⎩

1,2,...,j N =

1， 优化目标:

以悬停效率为最大目标函数 max ：FM （X ）。 2， 设计变量：

12(,,,,,,,,)r r pos pos X c p r r λεθδτ=

1）弦长变量：r c 、λ、ε、p 2）扭角变量：r θ、δ、τ 3）翼型配置位置：1pos r ，2pos r

1pos r 为第二种翼型的起始位置，2pos r 是第三种翼型的起始位置，这里优化选用了沿桨叶径向用三

种翼型，在现在的多数先进旋翼设计中一般三种翼型布置是常见的。

3，约束条件：

旋翼实度的约束，min σσ≥，为不使扭转角过大可以根据自己的设计需要对扭转角的上下限加以限制iL i iU θθθ≤≤，弦长约束iL i iU c c c ≤≤其中1,2,3...i =。在每次优化过程中，每一步迭代均用拉力系数T C 配平悬停时的总距0θ，保证每步迭代都会有正确解。

3.3 优化过程

优化过程流程图如下图2所示。本文采用二进制编码，取种群交叉概率为0.8，变异概率为0.01。种群规模为80，遗传代数为200代。

4 优化算例及结果

为验证优化的可行性，首先以计算的样例直升机优化模型，在只改变弦长和扭转角情况下，悬停效率由0.7413提高到0.7796，优化后的桨叶弦长和扭转角如下图3、4所示：

图3样例直升机桨叶弦长分布 图4样例直升机桨叶扭转角分布

弦长（m ）

半径站位 （r/R）

预扭角（度）

半径站位（r/R）