涡桨飞机襟翼载荷计算研究

魏浩格*，温 庆，王艳艳

（中航通飞华南飞机工业有限公司研发中心，广东 珠海）

引言

襟翼是机翼上用来改善气流状态和增加升力的装置。在飞机起飞、着陆、爬升或低速机动飞行阶段，偏转襟翼增加机翼剖面弯曲度及有效迎角[1]，增加机翼最大升力系数，实现增大升力的作用。涡轮螺旋桨发动机（简称涡桨发动机）以其功率大、运转稳定性好、寿命长、费用低[1]的优势被广泛应用于现代军民用飞机中。

然而螺旋桨桨叶旋转时会引起气流扰动，对飞机位于螺旋桨之后并处于滑流之中的部件产生滑流效应[2]，国内外学者对滑流影响开展了大量研究工作，左岁寒[3]、张小莉[4]等人采用准定常N-S 方法对螺旋桨产生的滑流影响进行了分析。

本文以某4 发螺旋桨飞机为例，通过工程算法对涡桨飞机考虑滑流影响的后缘襟翼载荷进行计算，将考虑滑流喷流影响的CFD 仿真计算结果与工程算法计算结果进行对比，得出：采用工程算法可较好反映涡桨飞机襟翼载荷情况，满足工程使用要求。

1 滑流效应工程算法

为了实现高效工程评估的目的，简化螺旋桨对气流扰动的影响效果，采用等效盘模型的无厚度圆盘代替真实螺旋桨[5]，获得近似真实螺旋桨的滑流效果，能够快速评估襟翼等部件滑流影响效果及规律。

螺旋桨通过增加流经螺旋桨桨盘区域流管的动量产生推力，可以比拟为翼面产生升力[2]。基于动量定理，采用一个等效盘代表螺旋桨建立螺旋桨滑流简化模型。模型考虑滑流收缩，不计滑流旋转，假定在滑流各截面轴向速度相等。

被滑流覆盖的襟翼部分产生更高的动压，从而比邻近襟翼部分产生更大的升力。如图1 所示。

图1 滑流覆盖的襟翼部分动压增加

桨盘前方远处动压：

桨盘后方被滑流覆盖的襟翼部分动压：

动量定理描述如下：

式（3）中：T 为螺旋桨拉力，N；

则式（3）可转换为：

桨盘平面处相对远前方来流的速度增量系数为：

被滑流覆盖的襟翼部分相对远前方来流速压比为：

滑流收缩截面面积与螺旋桨桨盘面积比为：

收缩截面直径与桨盘直径比为：

本文以4 发涡轮螺旋桨飞机为例，滑流对襟翼载荷的影响可假定为使其收缩截面覆盖的襟翼部分的速压从原始的Q 增至Q'，则单侧（两发）襟翼载荷的滑流修正系数可得：

代入式（6），式（10）可转化为：

将式(11)得到的滑流修正系数KS乘以襟翼无动力情况下载荷即可得到考虑滑流影响的襟翼载荷，将滑流覆盖襟翼部分（收缩截面直径D'）内的襟翼压力分布乘以速压变化比即可得到考虑滑流的襟翼压力分布。

2 计算输入

2.1 计算对象

本文研究飞机构型为4 发涡桨飞机，两侧各两个涡桨发动机，图2 为襟翼偏转后飞机外形示意。

图2 飞机外形示意

X 轴沿飞机水平基准线，逆航向为正；Y 轴垂直于机身对称面，逆航向向左为正；Z 轴在飞机对称面内垂直于横轴指向上方。

襟翼形式为定轴固定子翼双缝襟翼，在展向某截面处分为内、外襟翼，共4 片襟翼，图3 为襟翼形式切面示意。

图3 襟翼形式示意

2.2 计算工况

通过工程算法计算考虑滑流襟翼载荷，比较筛选出的襟翼载荷严重工况如表1 所示。

表1 严重工况

采用CFD 对严重工况进行考虑滑流、喷流计算，用以对比工程算法结果。CFD 验证计算使用的网格如图4 所示。

图4 襟翼局部网格

3 计算结果及比较分析

工程算法计算襟翼载荷以风洞试验的压力分布数据为基础，计算无动力情况的气动力分布，对于襟翼而言，气动分布载荷主要是指：

（1） 总压力中心(压心)位置。

（2） 展向力分布。

（3） 剖面压力中心。

（4） 剖面弦向载荷分布。

根据风洞得到的飞机表面压力分布（Cp）数据库，通过积分，计算襟翼的气动分布载荷，再按前文工程计算方法计算考虑滑流的滑流修正系数。另外，工程算法计算时认为左右襟翼载荷对称。

采用数值模拟方法对全机考虑滑流及喷流进行仿真计算，提取两侧襟翼CFD 计算气动压力分布数据，通过积分计算载荷情况。

3.1 展向力分布对比

将工程算法结果与基于CFD 方法考虑滑流、喷流影响仿真后计算所得的展向载荷分布结果进行对比，如图5- 图8 所示。

图5 主内襟翼气动力比较

图6 主外襟翼气动力比较

图7 子内襟翼气动力比较

图8 子外襟翼气动力比较

从图气动力分布可以看出，在螺旋桨后方，受到滑流影响的襟翼站位区域，气动力相对较大，在边缘区域的气动力较小。

3.2 工程算法与仿真对比

将滑流工程算法与CFD 计算进行对比。如图9 所示为工程算法与仿真计算气动力对比，分别为计算状态序号1 和2 的工程算法与CFD 计算右、左襟翼的主内、主外、子内、子外襟翼及合力结果。

图9 工程算法与CFD 气动力对比

如表2 所示，为对应图9 的CFD 与工程算法气动力的比值。

表2 CFD 与工程算法比值

从表2 可以看出，工程算法计算的气动力合力均大于采用CFD 计算，说明工程算法所算气动力可覆盖CFD 结果。

如图10 所示为工程算法与仿真计算铰链力矩对比，分别为计算状态序号1 和2 的工程算法与CFD 计算右、左襟翼的主内、主外、子内、子外襟翼及合力矩结果。

图10 工程算法与CFD 力矩对比

如表3 所示，为对应图10 的CFD 与工程算法铰链力矩的比值。

表3 CFD 与工程算法比值

从表3 可以看出，工程算法计算的铰链力矩小于CFD 计算结果，但考虑到CFD 中补充了喷流的影响，工程算法已经较好的反映出襟翼载荷情况，基本满足工程需求。

4 结论

本文使用风洞试验压力数据通过考虑螺旋桨滑流影响的工程算法计算襟翼载荷，筛选出严重工况。基于严重工况采用CFD 方法考虑滑流、喷流影响因素仿真，得到的飞机襟翼表面压力分布（Cp），通过积分，计算得到两侧襟翼的气动载荷。将仿真结果与考虑滑流工程算法的襟翼载荷严重情况进行对比，合力最大偏差为仿真结果较工程算法小4.8%，合力矩最大偏差为仿真结果较工程算法大5.4%。

采用工程算法更加便捷高效，可基本满足工程需求，较好的反映襟翼载荷情况，但在后续设计计算中，仍需通过CFD 仿真计算、试验对飞机进行校验。