倪章松, 刘森云, 王 桥, 王梓旭, 郭 龙

(中国空气动力研究与发展中心 空气动力学国家重点实验室, 四川 绵阳 621000)

0 引 言

结冰风洞[1-4]是开展飞行器结冰与防除冰研究的重要地面试验设备。历经70余年的建设发展，目前世界上已建成结冰风洞20余座[5]，主要分布在美国、意大利、中国、法国等。以美国IRT和意大利IWT大型结冰风洞为代表，世界上形成了欧美两大飞行器结冰试验研究中心[6-8]。IRT结冰风洞建成时间较早，试验技术发展十分完善，是目前公认最先进的结冰风洞之一。IWT结冰风洞是欧洲最主要的飞行器结冰试验基地，该风洞的试验技术也相对完善，国内多个型号在该风洞开展过结冰试验。为改变我国大型飞机防除冰设计依赖国外结冰风洞的现状，中国空气动力研究与发展中心(以下简称气动中心)于2013年建成了国内首座3 m×2 m大型多功能结冰风洞。该风洞由气动中心独立设计完成，拥有自主知识产权[9-12]。

结冰风洞包含制冷、云雾场、高度模拟等多个复杂的子系统。开展结冰与防除冰试验所涉及到的试验测试技术种类较多，例如风洞气流场、云雾场校测技术，常规结冰云雾环境与复杂结冰气象环境(SLD，冰晶等)模拟技术，防除冰试验技术，旋翼、发动机等复杂构件的结冰试验技术等[13-16]。结冰风洞试验技术是以多学科交叉理论为基础、以工程应用为目标的特种试验技术[17-20]。拥有先进可靠的试验技术是结冰风洞能力水平的重要体现，也是飞机气动性能评估和防除冰优化设计的必要保障。

为更好地发挥3 m×2 m结冰风洞在国家重大型号研制和国民经济建设中的作用，气动中心在结冰风洞试验技术领域开展了大量的研究工作。本文详细介绍了3 m×2 m结冰风洞云雾场参数校测技术，冰形特征捕获技术，热气防冰和电热除冰技术等关键技术，并对未来重点发展的进气模拟条件下的发动机结冰与防冰试验技术，直升机、螺旋桨结冰等特种试验技术进行了展望和探讨。本文的研究结果可为国内结冰风洞试验技术的发展规划提供参考。

1 3 m×2 m结冰风洞概况

1.1 总体情况

3 m×2 m结冰风洞是是一座闭口、高亚声速、回流式风洞，拥有主、次、高速3个可更换试验段，如图1所示。其中，主试验段宽3.0 m、高2.0 m、长6.5 m，次试验段宽4.8 m、高3.2 m、长9.0 m，高速试验段宽2.0 m、高1.5 m、长4.5 m。风洞配套有高度模拟系统、热气供气系统、进气模拟系统等各类试验系统及配套试验装置。该风洞是“十一五”国家重大科技基础设施之一，是现代飞机、高铁、风力机等型号研制必不可少的大型试验设施。

图1 3 m×2 m结冰风洞布局图

1.2 性能指标

表1给出了3 m×2 m结冰风洞的试验段尺寸和风速、湿度、温度、高度、云雾参数等主要性能指标参数[21]。需要注意的是，风洞的最大风速取决于所选用的试验段，其中在7000 m模拟高度条件下，主试验段的最大风速为210 m/s，次试验段的最大风速为78 m/s，高速试验段的最大风速为256 m/s。

表1 结冰风洞主要设计参数[21]Table 1 Main indexes of 3 m×2 m icing wind tunnel[21]

1.3 试验能力

3 m×2 m结冰风洞目前已形成的试验能力主要包括：(1) 飞机及其部件(包括机翼、尾翼、风挡、天线罩等)的结冰和冰脱落试验；(2) 飞行器结冰/防除冰试验与系统性能验证；(3) 发动机进气模拟条件下的结冰/防除冰试验；(4) 风力机、高速列车及地面设施的结冰特性试验；(5) 防/除冰新系统、新技术性能验证试验等。

3 m×2 m结冰风洞所有技术指标、结冰试验云雾场参数校测方法和流程均符合美国机械工程师学会颁布的相关标准(SAE ARP5905)[22]。2014和2018年，气动中心2次采用NACA 0012翼型开展结冰风洞标模试验，并与意大利CIRA-IWT结冰风洞进行了对比验证，同时将试验结果与SAE ARP5666参考文献[23]中其他同类型风洞的试验结果进行了对比，研究结果表明3 m×2 m结冰风洞结冰试验结果与国外主流工程性结冰风洞结果一致性良好，尤其是与IRT的结果较为接近(详见图2)。图中c表示弦长，图例CARDC为气动中心3 m×2 m结冰风洞试验结果，CIRA为意大利CIRA结冰风洞试验结果，Calc为数值计算结果，Cox为美国Cox结冰风洞试验结果，BRAIT为美国波音BRAIT结冰风洞试验结果，GKN为英国GKN结冰风洞试验结果，IRT为美国IRT结冰风洞试验结果，Body为NACA 0012翼型轮廓。

(a) 霜冰试验结果对比

(b) 明冰试验结果对比

2 试验技术现状

3 m×2 m结冰风洞已形成包括云雾场参数校测、冰形特征捕获、热气防冰和电热除冰等在内的常规结冰与防除冰试验技术，下面进行详细阐述。

2.1 云雾场参数校测技术

结冰风洞的云雾模拟技术以适航条例为指导，使风洞产生的云雾场与真实飞行条件下的云雾场相似。

针对云雾场均匀性，研制了结冰格栅。格栅整体尺寸为1955 mm(高)×2705 mm(宽)，格栅与上下洞壁相连，其面积覆盖了试验段中心88%的截面积。格栅单元格为150 mm×150 mm，深60 mm，迎风面宽度5 mm(见图3)。通过测量一定时间内格栅迎风面的结冰厚度，分析格栅表面结冰在风洞横截面的一致性，以此评估云雾场的均匀性。

针对云雾水滴直径参数，购置了马尔文粒径测量仪和相位多普勒干涉仪(PDI-FPDR)，如图4所示，并进行风洞内安装适应性配套，实现了云雾水滴粒径的测量。其中，马尔文粒径测量仪主要用于喷嘴的雾化特性和参数的地面测试，相位多普勒干涉仪主要用于风洞内云雾场水滴直径的直接测量。

图3 均匀度格栅结冰

图4 用于MVD校测的PDI-FPDR仪器

针对液态水含量参数，自主研制了气压驱动式冰刀装置和均匀度评估系统，提出了一种采用数值计算和结冰风洞试验相结合的手段测量结冰风洞中液态水含量的方法，实现了云雾场参数的高效、准确校测。

2.1.1 冰刀法

建立了基于气压驱动式冰刀的结冰风洞液态水含量校测装置和方法，气压驱动相比电驱动反应更迅速，且不易受低温环境的影响。冰刀法的基本原理是：将冰刀置于结冰风洞的来流中,一定时间后,工作面就会生长一定量的冰。试验结束后，通过测量冰的厚度和结冰时间等参数,采用式(1)计算液态水含量:

(1)

其中，LWC为液态水的含量值，C为单位转换常数，ρ为冰的密度，D为冰的生长厚度，Eb为冰刀收集效率(液态水收集系数)，v为试验风速，t为结冰时间。

冰刀法的关键是确保撞击到冰刀工作面的过冷液态水能够全部冻结,这就要求云雾场的温度要足够低(低于-18 ℃),此时冰刀上所结冰为霜冰。冰刀法测量液态水含量的方法为：将冰刀置于风洞中心线上,冰刀前装有防护罩(见图5)。防护罩准备状态处于关闭位置，将冰刀前缘保护起来。当结冰风洞流场稳定后,防护罩快速移开,暴露在云雾场中的冰刀前缘在液态水撞击下开始生长冰,试验结束后快速关闭防护罩保护生长的冰。防护罩作快速打开、关闭运动,开闭时间间隔即为结冰时间t。

图5 风洞中的冰刀装置[24]

冰刀装置的主要技术指标为：(1) 液态水收集系数Eb大于0.95；(2) 结冰时间控制精度要求防护罩作0°～90°角行程时间小于0.2 s；(3) 最大试验风速下关键部件的变形量：冰刀小于0.05 mm，防护罩小于1 mm；(4)冰刀装置阻塞度小于5%；(5)冰刀表面粗糙度为0.8 μm[25]。

冰刀法的优点是简便快捷，缺点是未考虑结冰过程中液态水蒸发的影响，同时为保证收集系数大于特定值，要求风洞云雾场的温度足够低。为克服冰刀法的缺点，提出了采用数值计算和试验相结合的手段测量结冰风洞中液态水含量的方法[26]，下文称其为混合法。

2.1.2 混合法

该方法无需专门的LWC测量装置，只要有常规的质量测量工具即可。其基本步骤为：(1) 选取典型外形(如圆柱)进行结冰风洞试验，测量一定速度和结冰时间下冰的质量；(2) 用数值模拟的手段计算不同液态水含量条件所对应的结冰质量，建立LWC与结冰质量的关系曲线LWC=f(Mi)；(3) 将试验测量的结冰质量在曲线上进行插值，获得试验LWC的值。

以直径为L的圆柱为例，位于极限轨迹之间的过冷水滴会与物面碰撞并结冰，而在极限轨迹范围之外的水滴将绕过物体。极限轨迹定义为分别与圆柱的上下表面相切的水滴运动轨迹。设流动速度为v，则单位时间、单位展长物面上冰的质量Mi为:

Mi=Em·LWC·v·L-Me

(2)

其中，Em=Δy/L为水滴总收集系数，Δy为2条极限轨迹之间的距离。蒸发质量Me的计算采用经验公式[26]。为了计算水滴的极限轨迹，需要完成空气的流场计算，水滴的运动方程求解，以及迭代寻迹等过程[26]。图6给出了风速为25和35 m/s条件下液态水含量与结冰质量的关系。

利用混合法测量LWC，操作方便，成本低廉，且测量结果的精度较高，克服了常规LWC测量方法的不足。为进一步提高云雾场参数测量的精准度，3 m×2 m结冰风洞未来主要开展云雾发生技术、云雾液滴运动过程研究、非接触式结冰云雾场诊断技术研究等。

图6 不同风速下的结冰质量[26]

2.2 冰形特征捕获技术

在冰形测量方面，获取冰形的传统方法是利用卡板及“热刀”装置(见图7)，手工测量冰形的二维轨迹线[27-28]，通过数据拟合的方式获取冰形的三维信息。该方法操作简便，容易实现，但操作过程中对冰形信息有所破坏，影响了测量的准确度，且所得结果为有限的二维冰形，不能完全反映冰形的真实三维外形。

图7 传统热刀法冰形测量技术

考虑到热刀法在复杂结冰表面冰形提取的局限性，发展了基于立体成像和激光扫描[29]的冰形非接触测量技术，通过基于金属粉末的显影方法，解决了冰形三维测量中的冰透光性问题；通过选用具有防水等级和耐低温的测量产品，解决了低温高湿环境的影响问题，实现了冰形的三维激光扫描，同时建立了三维数字冰形提取、重构与分析技术[30-32]，见图8。

2.3 热气防冰和电热除冰试验技术

热气供气系统是热气防除冰试验的关键[33-36]，3 m×2 m风洞目前已完成热气供气系统建设，见图9。结合基础设备平台，配套了热气温度、压力、流量测控系统，形成了热气防冰试验能力。热气供气系统的主要技术指标如下：

(1) 防冰热气最大流量：1.5 kg/s，精度±1%；

(2) 防冰热气最高温度：400 ℃，精度±2 ℃；

(3) 机翼模型防冰热气最大压力：0.4 MPa，精度±1.0%；

(4) 发动机模型防冰热气最大压力：1.4 MPa，精度±1.0%。

(a) 实际冰形

(b) 扫描冰形

Fig.8Iceshapenon-contactmeasurementtechnologyusinglaserscanning

图9 热气供气系统

系统通过控制热气防冰模型入口热气的温度、压力和流量参数，实现不同防冰热气条件的模拟。该系统已应用于C919飞机热气防冰试验。

3 m×2 m结冰风洞已完成电热防除冰试验系统的研制。由于需要承接各种飞机型号试验，该电加热防除冰系统配套了直流28 V、270 V、交流115 V/400 Hz中频电源，供电主回路采用PLC控制固态继电器方案，每套电源可以实现20个加热区温度控制，每个加热区加热功率不低于2 kW，控制方式可以选择通断控制间断加热或调节电压连续加热，满足型号单位不同试验需求[37-38]。该装置已应用于某直升机电热除冰技术研究。

3 试验技术展望

在结冰风洞试验技术研究方面，下一步将重点关注发动机结冰与防除冰、直升机结冰与防除冰，过冷大水滴和冰晶雨雾模拟等结冰风洞试验新技术，并形成相应的试验能力。

3.1 发动机结冰与防除冰试验技术

发动机结冰与防除冰试验技术主要包含发动机进气道模拟技术、复杂构型结冰测量技术、结冰特性捕获技术、结冰危害评估技术等。

图10给出了结冰试验在发动机防除冰系统设计中的地位和作用。由于发动机结冰数值模拟研究具有高度的复杂性和结果预测的不确定性，因此结冰风洞试验在发动机设计环节中是必不可少的。

图10 发动机结冰与防除冰研究关系图

Fig.10Researchmapoficingandde-icingtechniquesonaircraftengine

3 m×2 m结冰风洞目前已完成发动机进气模拟系统主体建设，建设了主、辅2条引气支路分别用于模拟大小流量，抽吸气流通过进气道与发动机舱，实现发动机进气道入口的流量模拟[39]，从而模拟发动机的不同功率状态。结合基础研究平台配套了动力系统和测控系统，形成了发动机进气道模拟试验能力(见图11)，已顺利完成某型发动机进气道结冰试验。

图11 发动机进气模拟系统

发动机结冰与防除冰试验技术研究未来的重点工作包括：(1) 加强对结冰试验参数的控制，特别是LWC和MVD的控制；(2) 积极研制物理尺寸小、测量精度高、测量范围大、对流场干扰小的测量装置，特别是利用光学或声学原理的测量手段[40]。

3.2 直升机结冰与防除冰试验技术

在直升机结冰与防除冰试验技术研究方面，需要建设旋翼结冰试验台，并构建高速旋转模拟与控制技术。2019年，3 m×2 m结冰风洞将建成2 m直径直升机旋翼结冰试验台并投入使用。试验台各系统的技术指标如下：

(1) 桨毂中心至试验段底端的高度范围：1.60～2.14 m；

(2) 旋翼直径：2 m；

(3) 最大驱动功率：90 kW；

(4) 额定转速：2100 r/min；

(5) 最大拉力：3500 N；

(6) 主轴倾角范围：-25°～15°。

通过控制直升机旋翼的不同转速、主轴倾角等模拟旋翼的不同工作状态。受风洞尺寸限制，采用缩比旋翼进行结冰试验研究还面临诸多问题，因此，直升机旋翼的结冰与防除冰试验技术未来重点的发展内容包括[41-42]：

(1) 旋翼结冰试验相似准则研究；

(2) 旋翼结冰风洞试验模型设计技术；

(3) 旋翼升力模式控制技术(定升力、定总矩)；

(4) 旋翼结冰过程冰形在线捕获技术；

(5) 旋翼参数、气动性能测量技术(形变、升力、扭矩)；

(6) 旋翼冰脱落测量与分析技术；

(7) 旋翼防除冰试验技术。

3.3 过冷大水滴云雾模拟试验技术

过冷大水滴具有平均直径大、分布范围广、液态水含量高、结冰环境温度较高的特点。由于过冷大水滴造成的结冰会导致飞机的气动性能和操纵性能急剧下降，对飞行安全造成极大的威胁，2014年11月4日，FAA颁布了有关SLD结冰的补充条例。在FAR 25部增加了附录O，明确了冻细雨(Freezing Drizzle)、冻雨(Freezing Rain)等SLD云雾条件，该修正案于2015年1月5日生效。同年3月12日，EASA在16号修正案中增加了相关内容，SLD与冰晶的分布高度及危害如图12所示[43]。

SLD试验技术的主要内容是冰晶条件生成、测试仪器校准、大粒径喷雾装置研制、云雾模拟与校测、云雾水滴分布符合性评估以及结冰/防除冰试验验证。未来这些技术的建立与以下基础性物理问题的解决密不可分，主要包括：

图12 SLD和冰晶大气特征示意[43]

(1) 喷嘴大粒径水滴喷雾机理与实现；

(2) 水滴沿程热力学特性(直径、温度)变化规律；

(3) 水滴沿程动力学特性(变形、破碎、轨迹)变化规律；

(4) 水滴特征参数(粒径、温度)非接触式光学测量；

(5) 喷嘴分布、控制策略与云雾粒径、温度分布影响规律。

SLD试验技术的核心是需要模拟适航条例FAR25部附录O中过冷大水滴的粒径分布，并在有限条件下建立满足符合适航符合性验证要求的验证方法[44-45]。

3.4 冰晶云雾模拟试验技术

自2002年，全球已有超过14起因冰晶造成的飞机发动机失效的严重事故[46]。针对冰晶问题，美国联邦适航条例33部附录D规定了典型的冰晶结冰云雾条件。在此基础上，美国NASA的推进试验室(PSL)于2013年形成了发动机冰晶试验能力,并在同年开展了第一次试验。由此看来，我国自主研发的发动机在未来适航取证的过程中,同样将会面临冰晶结冰的适航合格审定问题。

冰晶结冰试验技术研究需要回答3个问题[46]：

(1) 冰晶是如何在发动机内部结上的？

(2) 冰晶结冰是如何影响发动机正常工作的？

(3) 如何才能防止或者减轻发动机发生冰晶结冰？

针对这3个问题，主要开展以下研究内容，如图13所示。其中，对于冰晶试验技术研究，需要重点突破的关键技术包括[47-48]：

(1) 冰晶喷雾系统研制；

(2) 冰晶的风洞探测技术；

(3) 冰晶结冰条件模拟技术(冰晶热力学、相变、动力学、侵蚀等物理特性研究)；

(4) 冰晶结冰试验相似方法；

(5) 冰晶结冰测量仪器评估技术；

(6) 冰晶结冰非接触式测量技术。

在技术指标方面，未来可以参考NASA PSL风洞[49]，如表2所示。

图13 发动机冰晶结冰研究内容

Fig.13Researchcontentofcrystal-icingtesttechnologiesforaircraftengine

表2 冰晶模拟技术指标Table 2 Crystal-icing simulation technical indexes

4 结 论

(1) 介绍了云雾参数校测、冰形测量、热气/电热防除冰试验等技术，并在型号试验中得到了应用；

(2) 提出了3 m×2 m结冰风洞试验技术的近期的发展方向，主要包括：发动机结冰与防除冰，直升机旋翼结冰与防除冰，过冷大水滴、冰晶云雾模拟试验技术等。