0.6m连续式跨声速风洞总压控制策略设计

黎壮声，杨鹏程，陈 旦，张永双，郭守春，陈天毅

0.6m连续式跨声速风洞总压控制策略设计

黎壮声，杨鹏程，陈 旦＊，张永双，郭守春，陈天毅

（中国空气动力研究与发展中心，四川绵阳 621000）

总压作为连续式风洞控制系统关键指标之一，其控制精度及快速性对提高风洞试验效率、降低能耗具有重要意义。0.6m风洞为国内首座具备负压试验能力的连续式跨声速风洞，其试验工况多，压力范围广，针对该风洞压力特性，设计了总压控制策略，根据不同的压力工况确定不同的阀门组合控制方式；同时针对模糊PID对连续式跨声速风洞宽压力范围、多调节工况下压力控制适应性较差的问题，提出分段变参数加模糊PID相结合的控制算法，即先根据目标总压确定不同分段区间下基本合理的基准P、I参数，再结合模糊控制算法对基准参数进行修正。风洞调试结果表明，总压控制精度优于0.1%，控制策略能够有效满足不同工况的控制要求。

连续式跨声速风洞；总压；控制策略；分段变参数；模糊PID

0 引 言

0.6 m×0.6m连续式跨声速风洞（以下简称0.6m风洞）为国内首座具备负压实验能力的连续式跨声速风洞，试验马赫数范围为0.2～1.6，压力范围为20～250kPa（绝对压力），具有3个试验段，该风洞通过阀门控制系统实现常压、增压和降压运行，其中增压、常压运行通过调节进排气系统实现，降压运行通过调节抽真空系统实现。

风洞的总压调节能力及精度是风洞的一个重要技术指标，对提高试验效率和试验数据精准度具有重要作用。目前，国内连续式跨声速风洞数量较少，NF－6风洞作为国内首座连续式跨声速风洞，其马赫数范围较窄，试验工况相对较少［1－2］。目前国内外风洞对马赫数控制方法描述较多，而对连续式风洞总压控制及马赫数和总压的耦合关系描述较少［3－11］。暂冲式风洞压力控制通过主调压阀和引射器调压阀等实现［9－11］，而0.6m连续式风洞压力主要通过进排气和抽真空阀门／泵组合来实现，且总压调节受马赫数调节手段影响，可借鉴的内容较少。各种改进的模糊PID控制方法在风洞中尤其是马赫数或风速控制的应用较多［12－13］，文献［13］对某暂冲式风洞主气流压力控制原理进行了简单分析，并引入模糊控制概念，总压控制精度达到0.5%，因此本文拟采用模糊PID来进行总压控制，但其控制参数大幅跳动时容易导致控制器输出跳动，进而导致压力波动，因而需要结合0.6m连续式风洞压力特性对控制算法进行优化。

0.6 m风洞阀门系统结构见图1。

图1 阀门系统结构示意图Fig.1 Structure sketch of valve system

不同于现有其他风洞，该风洞压力调节系统具有以下特性：

（1）总压调节手段多，包含9台阀门和3台泵，不同的压力工况采用不同的阀门组合进行控制。

（2）压力控制范围广，从绝对压力20kPa至250kPa，涉及常压、增压、负压3种不同的工况，且试验过程中可连续变总压、连续变马赫数和连续变迎角［］。

（3）压力控制精度要求高，设计指标为0.2%（常压和增压）［15］。

（4）作为大型连续式风洞的引导性风洞，风洞运行过程中迎角、二喉道、压缩机转速、驻室抽气流量等随时需要调节，这些因素都对总压有较大影响，从而导致压力场波动较大，需要总压控制策略具备很强的适应能力［14－15］。

本文主要研究内容为：分析0.6m连续式风洞的压力控制特性，摸索多试验工况、宽压力调节范围下具有较好适应性的压力控制策略，并结合控制对象特性通过对模糊PID算法分段变参数的方式进行改进，并实施该风洞的压力控制，为风洞控制策略优化及试验效率的提高奠定基础，同时为后续的大型风洞建设积累经验。

1 风洞总压和马赫数耦合特性

连续式跨声速风洞主要通过调节进气／排气／抽真空阀门的开度来实现总压的调节，由于总压和马赫数在调节时相互耦合，致使调节难度加大，且调节时间较长，经过分析，总压和马赫数调节手段的耦合关系如图2所示。

图2 总压、马赫数耦合关系Fig.2 The coupling relation between total pressure and Mach number

其中，G1（s）至G9（s）分别为每种控制量的传递函数，Gab（s）为总压对马赫数的影响函数，Gba（s）为马赫数对总压的影响函数，Ga（s）为总压控制传递函数，Gb（s）为马赫数控制传递函数。

为提高总压控制精度并减少调节稳定时间，理论上需要对总压和马赫数进行解耦控制，但目前总压和马赫数的解耦算法在工程应用中较少，更多是通过控制策略的优化并尽可能减少耦合的方式，因此，本系统设计时将多种马赫数调节手段调节时对总压的影响作为总压控制的扰动对待，同时风洞马赫数控制策略设计时减少总压和马赫数的耦合。0.6m风洞总压控制主要通过调节相应的调节阀开度实现，因此需首先建立调节阀开度和风洞总压的关系。对于可压缩流体，理想情况下，流量Q和压力之间关系为［16］：

其中，β为缩流处截面等效直径与管道截面等效直径之比；E为与流速有关的系数；A为缩流处截面积；ρ为密度；Q为体积流量；p1为入口静压；p2为出口静压，近似等于大气压；ε为流体膨胀系数。

同时，对于抛物线特性［17］的调节阀，在理想情况下其流量和阀门开度之间满足下式：

其中，Qmax为流经阀门的最大流量，R为阀门的可调比，为流经阀门的最大流量和最小流量的比值，国内调节阀R一般取为30或50，l／L为阀门开度的百分比。

由式（1）和（2）可得调节阀开度和静压的关系，而风洞总压由静压和动压组成，进而得到调节阀开度和风洞总压之间的关系：

其中，p0为总压，f1（n）为风速和压缩机转速的函数关系，考虑到稳定段气流速度较低，可忽略不计，f2（α，r，l，m）为马赫数调节对总压的影响量，a为迎角变化梯度，r为压缩机静叶角，l为二喉道位移，m为驻室抽气量。

由此可见，连续式风洞总压影响因素包括：马赫数调节方式、气流介质密度（由风洞内总压和总温确定）、调节阀开度、气源压力等，本文重点对气源压力、阀门开度、压缩机转速和目标总压进行考虑以设计总压控制策略。

2 总压控制策略

为保证总压控制精度并确保风洞运行安全，考虑到PID算法具有结构简单、稳定性好、工作可靠、调整方便的优点，但同时存在自适应能力差的缺点，因此拟采用模糊PID算法，其对压力波动有较强适应性，但由于该连续式风洞压力范围较广，且不同目标压力区间下压力特性不同，若单纯采用模糊PID进行控制，可能导致不同压力工况参数匹配度较差，从而带来震荡或超调，因此该连续式风洞总压控制的核心算法采用分段变参数加模糊PID控制。

结合该连续式跨声速风洞压力变化规律，以及连续式风洞运行时压缩机充气密封系统将向风洞洞体内漏气的特性，为节省气源耗气量，确定该风洞总压控制策略为：

（1）当工况为常压（控制压力目标在大气压至105kPa）时，由于压缩机充气密封往洞体漏气，因此仅利用排气调节阀进行总压控制；

（2）当工况为增压（控制压力目标大于105kPa）时，通过预置进气阀开度和调节排气阀进行总压控制；

（3）当工况为负压（控制压力目标低于94kPa）时，开启真空泵，通过预置进气调节阀开度、开启真空泵和调节抽真空调节阀进行总压控制；

（4）当压缩机转速大幅上升或下降时，压力会大幅波动，通过计算目标转速和当前转速误差以及转速误差变化率，对参与调节的阀门采用分段限幅跟随调节的方式，以减少积分饱和和超调，同时为减少总压和马赫数耦合，变马赫数时采用多阶梯升压缩机转速的控制方式。

（5）压力调节过程中，进气阀预置开度与气源压力和目标压力满足调试确定的函数关系，排气主调阀进行压力粗调（0.2%以上），待压力稳定后由排气辅调阀精调（0.2%以下至0.04%）。

3 分段变参数加模糊变参数PID算法

3.1模糊PID控制算法

模糊PID控制算法以常规PID控制为基础，系统首先计算被控量的误差e和误差变化率ec作为二维模糊控制器的输入变量，然后通过模糊化、模糊控制规则推理、解模糊得到KP，KI和KD的修正量，在控制过程中在线实时整定PID各参数，其原理已经非常成熟，这里不再累述。

3.2分段变参数加模糊变参数PI算法

由于风洞总压控制范围较宽，在不同的压力区间，风洞的总压特性不同，尤其当总压较高时（p0≥200kPa），风洞排气调节阀开度增大后总压下降明显，同时目前调节阀存在响应时间及执行速度较慢的问题，当压力快速变化时阀门的调节速度跟不上“节奏”，因此若不根据压力控制范围对控制参数进行调整，而仍使用常压时的基准调节参数，必然导致PID算法适应性差，产生较大的超调。因此，以模糊自适应PID为基础，增加目标总压作为输入变量，首先根据目标总压进行分段，基于不同的总压控制目标结合简单PI调试规律初步确定PI参数基准值，尔后由模糊算法进行变参数控制，以提高不同试验状态的自适应能力，原理框图见图3。

图3 分段变参数加模糊PID控制框图Fig.3 Block diagram of valve control system using segmented parameter and fuzzy PID

3.3模糊控制参数及规则表

阀门控制系统中，将计算出的目标压力、实际压力和给定压力的误差e和误差变化率ec作为模糊控制器的输入变量，PI控制器的2个参数KP、KI的修正量ΔKP、ΔKI作为输出变量。误差e基本论域取［－0.75，0.75］，其偏差变化率ec的基本论域取［－0.5，0.5］，输出变量ΔKP、ΔKI基本论域分别取为［－0.5，0.5］、［－1.5，1.5］，E、EC、ΔKP、ΔKI分别是输入变量e、ec、输出变量ΔKP、ΔKI所对应的模糊语言变量。E、EC、ΔKP、ΔKI都分别取7个模糊子集，其中NB、NM、NS、ZO、PS、PM、PB分别代表负大、负中、负小、零（又分成正零ZP、负零ZN）、正小、正中、正大。模糊语言变量选择三角隶属函数。

调试过程先摸索出不同压力区间段的大体P、I参数，同时为防止控制输出突跳，将控制参数进行拟合处理后再作为模块控制的输入，见式7。

式中KP、KI为分段后基准比例系数和积分系数，p为目标压力。根据不同目标压力分段，能有效缩小适合当前目标压力的基准P、I值范围，基准P、I波动范围缩小，也就缩小了最终的P、I参数波动的绝对值，进而避免了控制参数大幅跳动的问题。

表1给出了结合风洞压力特性的分段变参数模糊PID控制规则表。

PID的实时控制参数是在基准P、I参数的基础加上模糊算法得到的P、I参数增量。

表1 阀门控制系统ΔKP、ΔKI参数Table 1 ΔKP、ΔKIparameter of valve control system

3.4进气阀门开度预置

由于连续式风洞有大量用于气密封的空气往洞内泄露，因此常压时采用气密封漏气加排气阀调节的方式（进气阀关闭）；而在增压时尤其高目标压力段，当压力欠压时，若仅通过气密封漏气的方式进行补气则速度过慢，需通过进气阀门对风洞进行补气，调试时为尽量减少气源压力波动对总压控制的影响，结合阀门性能曲线，减少气源消耗，进气阀以一个相对固定的进气量往风洞灌气，进气阀预置开度与气源压力和目标压力满足下列函数关系：

4 软件流程设计

程序设计时首先基于模块化编程方法，将阀门和泵进行组合以实现不同压力工况控制，风洞运行时，根据工况命令调用不同的模块，以满足多种工况下阀门开关组合的控制逻辑需求，并通过压力控制模块实现闭环控制，压力闭环由阀门系统PLC 315实现，风洞主控系统和阀门控制系统通过Profinet工业以太网进行通信，阀门控制系统软件流程见图4。

5 试验验证

为测试控制算法性能，对变总压、变马赫数、变迎角等不同工况均进行了调试，迎角范围为－4°～14°，图5和6列举了不同工况下的典型测试结果，其中图5为风洞在常压和增压工况下变马赫数变迎角调试结果，图6为负压20kPa变马赫数变迎角试验结果。

图5 变马赫数变迎角调试结果Fig.5 The debug result with varying Mach numbers and attack angles

由图5和6可见，趋于稳态后，常压100kPa时，总压最大值100.071，最小值99.923，平均值100.002，［18］计算控制精度为0.079%，增压200kPa时，总压最大值200.152，最小值199.83，平均值200.01，控制精度为0.09%；负压20kPa时，总压最大值20.028，最小值19.972，平均值19.998，控制精度为0.15%；具体测试结果见文献［19］。因此，总压控制精度在常压增压工况时优于0.1%，在负压工况时优于0.2%，达到国军标先进指标［20］，图中马赫数上升过程的超调主要由于调节阀的响应速度较慢导致。

图6 20kPa变马赫数变迎角调试结果Fig.6 The debug result with varying Mach numbers and attack angles at the pressure of 20kPa

6 结 论

经过试验验证，可以得出以下结论：

（1）利用分段变参数加模糊PID控制方法实现了0.6m风洞总压精确控制，总压控制精度优于设计指标要求，该控制方法对较宽压力范围或调节阀性能比较敏感的风洞压力控制是有效的；

（2）如果进一步提高各调节阀的控制精度及响应速度，总压控制精度和稳定速度将进一步提升；

（3）本文在0.6m风洞探讨的风洞总压控制策略可为后续大型风洞及其他连续式风洞的压力控制方法提供参考。

下一步还将重点对总压和马赫数之间，不同马赫数调节手段与总压控制之间的耦合关系等方面展开更深入的研究。

［1］张永双，陈旦，陈娇.NF－6连续式跨声速风洞马赫数控制方式比较与研究［J］.实验流体力学，2013，27（2）：95－99.

Zhang Y S，Chen D，Chen J.Comparison and research on the Mach number control methods for the NF－6continuous transonic wind tunnel［J］.Journal of Experiments in Fluid Mechanics，2013，27（2）：95－99.

［2］郝礼书，乔志德，张永双，等.NF－6风洞马赫数闭环控制系统设计研究［J］.实验流体力学，2010，24（4）：85－88.

Hao L S，Qiao Z D，Zhang Y S，et al.Design research on the Mach number closed－loop control system in the NF－6wind tunnel［J］.Journal of Experiments in Fluid Mechanics，2010，24（4）：85－88.

［3］Jonathan R Potts1，Ian Lunnon，William J Crowther.Development of a transonic wind tunnel test bed for MEMS flow control actuators and sensors［R］.AIAA－2009－319，2009.

［4］Nhan Nguyen，Mark Ardema.Adjoint method and predictive control for 1－D flow in NASA ames 11－foot transonic wind tunnel［R］.AIAA－2006－1433，2006.

［5］Edward J Mickle，Capt Carmen Evans.Restoration of AEDC′s 16－ft transonic propulsion test capability［R］.AIAA－2005－7614，2005.

［6］Tirres C，Baxter L E.AEDC test facility planning“toward new horizons”to today［R］.AIAA－2004－2731，2004.

［7］廖达雄，陈吉明，彭强，等.连续式跨声速风洞设计关键技术［J］.实验流体力学，2011，25（4）：74－78.

Liao D X，Chen J M，Peng Q，et al.Key design techniques of the low noise continuous transonic wind tunnel［J］.Journal of Experiments in Fluid Mechanics，2011，25（4）：74－78.

［8］连晓飞.基于遗传算法优化BP网络的风洞马赫数控制研究［D］.沈阳：东北大学，2011.

Lian X F.Wind tunnel mach number control based on genetic algorithm optimized bp－neural network［D］.Shenyang：Northeastern University，2011.

［9］高川，周波，蒋婧妍，等.基于labview的大型超声速风洞总压测控系统设计与应用［J］.测控技术，2014，33（8）：84－87.

Gao C，Zhou B，Jiang J Y，et al.Design and implementation of total pressure measurement and control system of large－scare supersonic wind tunnel based on LabVIEW［J］.Measurement＆Control Technology，2014，33（8）：84－87.

［10］芮伟，易凡，杜宁，等.2.4m跨声速风洞颤振试验流场控制技术研究［J］.实验流体力学，2012，26（6）：83－86.

Rui W，Yi F，Du N，et al.Study on flow field control technique of flutter test in 2.4mtransonic wind tunnel［J］.Journal of Experiments in Fluid Mechanics，2012，26（6）：83－86.

［11］褚卫华，汤更生，王帆.2m×2m超声速风洞流场控制策略研究与实现［J］.实验流体力学，2012，26（5）：98－102.

Chu W H，Tang G S，Wang F.Research and realization on the control strategies of the 2m×2msupersonic wind tunnel［J］.Journal of Experiments in Fluid Mechanics，2012，26（5）：98－102.

［12］李可，刘旺开，王浚.专家＿模糊PID在低速风洞风速控制系统中的应用［J］.北京航空航天大学学报，2007，33（12）：1387－1390.

Li K，Liu W K，Wang J.Parameters self－tuning fuzzy－PID combined with expert control on wind velocity control system of wind tunnels at home［J］.Journal of Beijing University of Aeronautics and Astronautics，2007，33（12）：1387－1390.

［13］杨海滨，张伟，罗承友，等.模糊控制在风洞主气流压力自动调节系统中的应用［J］.兵工自动化，2015，34（4）：39－43.

Yang H B，Zhang W，Luo C Y，et al.Application of fuzzy control in wind tunnel main airflow pressure auto－adjust system［J］.Ordnance Industry Automation，2015，34（4）：39－43.

［14］张永双.某型连续式跨声速风洞测控系统初步设计方案［R］.中国空气动力研究与发展中心，2008.

Zhang Y S.The design scheme of measurement and control system for one type of continuous transonic wind tunnel［R］.China Aerodynamics Research and Development Center，2008.

［15］陈旦.某型连续式跨声速风洞阀门控制系统技术设计［R］.中国空气动力研究与发展中心，2012.

Chen D.The technological design of measurement and control system for one type of continuous transonic wind tunnel［R］.China Aerodynamics Research and Development Center，2012.

［16］安延涛.大型压力调节阀的动态分析及故障检测研究［D］.山东：山东大学，2012.

An Y T.Dynamic analysis and fault detection of large pressure control valve［D］.Shandong：Shan dong University，2012.

［17］明赐东.调节阀计算选型使用［M］.成都：成都科技大学出版社，1999.

Ming C D.Adjusting valve calculation lectotype use［M］.Chengdu：Chengdu University of Science ＆Technology publishing house，1999.

［18］国防科工委.GJB 1179－91.高速风洞和低速风洞流场品质规范［S］.北京，1992.

State Commission of Science and Technology for National Defense Industry.GJB 1179－91.Specification for flow quality of high and low speed wind tunnels［S］.Beijing，1992.

［19］陈吉明，柳新民，张永双.某型连续式跨声速风洞设计技术总结［R］.中国空气动力研究与发展中心，2015.

Chen J M，Liu X M，Zhang Y S.The summarization of design for one type of continuous transonic wind tunnel［R］.China Aerodynamics Research and Development Center，2015.

［20］总装备部.GJB 5221－2004.风洞控制系统设计及检定准则［S］.北京，2004.

General Armament Department.GJB 5221－2004.Design and calibration criteria of control systems for wind tunnels［S］.Beijing，2004.

The design of total pressure control strategy for 0.6mcontinuous transonic wind tunnel

Li Zhuangsheng，Yang Pengcheng，Chen Dan＊，Zhang Yongshuang，Guo Shouchun，Chen Tianyi
（China Aerodynamics Research and Development Center，Mianyang Sichuan 621000，China）

The total pressure is one of the key indicators of the continuous tunnel control system，the control precision and speed of which are very important to improve the efficiency of wind tunnel test and reduce energy consumption.The 0.6mwind tunnel is the first domestic continuous transonic wind tunnel with the negative pressure test ability，which can provide lots of test conditions covering wide total pressures.Based on its pressure condition，the total pressure control strategy is designed，and the different control ways of valve combination are ascertained depending on the different pressure conditions.Meanwhile，the control algorithm of segmented parameters combining with the fuzzy PID is proposed，to improve the poor adaptability of fuzzy PID to the wide pressure range and multi－adjustment condition in continuous transonic wind tunnels.In our algorithm，the reasonable benchmark parameters P and I under the different segmentation interval due to the target total pressure should been determined firstly.Then it would been modified by fuzzy control algorithm.The result shows that the total pressure control accuracy is better than 0.1%，and the control strategy can adapt to the control requirement of different working conditions.

continuous transonic wind tunnel；total pressure；control strategy；segmented variable parameter；fuzzy PID

V211.74

：A

（编辑：李金勇）

1672－9897（2016）04－0087－06

10.11729／syltlx20150129

2015－11－04；

2016－03－25

＊通信作者E－mail：cd0606＠163.com

Li Z S，Yang P C，Chen D，et al.The design of total pressure control strategy for 0.6mcontinuous transonic wind tunnel.Journal of Experiments in Fluid Mechanics，2016，30（4）：87－92.黎壮声，杨鹏程，陈旦，等.0.6m连续式跨声速风洞总压控制策略设计.实验流体力学，2016，30（4）：87－92.

黎壮声（1975－），男，广西玉林人，高级工程师。研究方向：风洞控制。通信地址：四川省绵阳市二环路南段6号403信箱（621000）。E－mail：livdug＠163.com。