一种缩比模型自由飞试验失速判断方法_粟建波.pdf

收稿日期:20210301修回日期:20210307第 40 卷第 2 期计算机仿真2023 年 2 月文章编号:10069348(2023)02003904一种缩比模型自由飞试验失速判断方法粟建波1，张甲奇2，邱宇豪2，王锦锦2(1 桂林电子科技大学，广西 桂林 541704;2 中国飞行试验研究院中航工业飞行仿真航空科技重点实验室，陕西 西安 710089)摘要:针对某通航类飞机缩比模型自由飞失速特性试验中，难以准确判断飞机模型进入失速的问题，根据中国民用航空规章第 23 部中失速相关条款的要求，提出了依据俯仰角速率变化判断飞机模型出现不可控下俯运动趋势，进而判断飞机模型失速的方法，通过建立该飞机六自由度非线性仿真模型进行仿真验证，同时利用缩比模型自由飞飞行试验对此方法进行了验证。试验结果表明:上述方法能够准确的判断飞机进入失速，进而提高了缩比模型自由飞试验的精准度。关键词:模型自由飞;失速特性;失速判断中图分类号:V249.1文献标识码:BA Judgment Method for Stall in FreeFlight Testof educed atio ModelSU Jianbo1，ZHANG Jiaqi2，QIU Yuhao2，WANG Jinjin2(1 Guilin University of Electronic Technology，Guilin Guangxi 541704，China;2 AVIC Aeronautical Science and Technology Key Laboratory of Flight Simulation，CFTE，Xi an Shanxi 710089，China)ABSTACT:Aiming at the problem that it is difficult to accurately judge whether the airplane model entering thestall during the freeflying stall characteristics test of a general aviation aircraft scaled model，according to the re-quirements of the stall related clauses in Part 23 of the Civil Aviation egulations of China，a method was proposed forjudging the uncontrollable pitching motion trend of the airplane model based on the change of the pitch angle rateThen the method of judging the stall of the aircraft model was verified by establishing a sixdegreeoffreedom non-linear simulation model of the aircraft At the same time，the method was verified by the free flight test of the scaleddown model The test results show that this method can accurately determine that the aircraft has entered a stall，thereby improving the accuracy of the scaled model free flight testKEYWODS:Model freeflight;Stall characteristics;Stall judgment1引言缩比模型自由飞试验是遵循几何、质量特性、动力相似准则关系，使用全尺寸飞机缩比模型，在真实大气环境中开展的飞行试验，主要用于飞机失速速度、失速特性以及尾旋特性研究，其试验结果与全尺寸飞机具有很好的相关性，在国内外航空研究院所获得广泛的认可1。美国 NASA 曾利用缩比模型自由飞试验完成了 F15、F16、F117、F18 等多型飞机的大迎角特性研究23，国内航空工业试飞中心也先后完成了运10、AJ、HO300、C919、AG600 等飞机的缩比模型自由飞试验研究，试验结果降低了新机研制风险和成本，为全尺寸飞机的试飞和适航提供了重要的技术支持。CCA23 中，23.201(b)条款规定了机翼水平失速要求:飞机的机翼水平失速特性必须按下述要求在飞行中进行演示:在至少高于失速速度 10 节开始，必须先拉升降舵操纵器件使减速率不超过每秒一节，直到失速发生，可用下列任一表明:飞机出现不可控制的下俯运动;防失速装置(如:推杆器)激发了飞机的下俯运动;操纵器件达到止动点45。条款中规定“飞机出现不可控的下俯运动”即可判断飞机失速，在全尺寸飞机的失速试验中，驾驶员可以实时感知飞机的过载、姿态、声音、机体抖震等变化，分析飞机出现不可控的下俯运动，进而判断飞机失速。而在模型自由飞93试验中，地面操纵手主要依靠目视和下传的飞行参数对飞机模型状态进行判断，存在较大的视觉误差和时间延迟，不能及时准确地判断出飞机模型出现了不可控的下俯运动。在多个飞机型号的模型自由飞试验中发现，依靠操纵手目视和下传的飞行参数进行飞机失速的判断，再进行改出操作，此过程存在 3 5s 的时间延迟，飞机可能进入过失速或深失速5。因此，国内外各研究院所和机构均在探索利用飞行控制系统自主实现模型失速判断的方法。本文依据 CCA23 中，23.201(b)条款要求，创新性的设计了一种基于俯仰角速率判断飞机模型失速的方法，并通过飞行仿真和缩比模型飞行试验对此方法进行了验证，获得了可信的试验结果。2失速判断方法俯仰角速率判断模型失速的方法是基于中国民用航空规章第 23 部 23.201(b)条款中“飞机出现不可控制的下俯运动”而设计的。当持续缓慢拉杆使模型进入失速的过程中，俯仰角速率为正值，俯仰角持续增大，迎角持续增大。在该过程中，如模型进入失速，则模型出现“低头”现象，此时俯仰角速率由正值变为负值，即持续拉杆的过程中出现不可控的下俯运动。同时，考虑到模型大迎角状态时可能出现震动现象以及传感器精度，故不能以单个俯仰角速率值作为失速判断的依据，需采集 t 秒时间内 n 个俯仰角速率的平均值?q作为判断依据，如式(1)?q=q0+q1+qnt(1)当 t 秒内俯仰角速率平均值?q小于设置的失速俯仰角速率均值 qs时，认为模型在持续拉杆的情况下出现了低头现象，即飞机出现不可控的下俯运动，从而判断模型已失速。为了更好的模拟飞行员失速操纵，判断模型失速后继续拉杆1 秒再推杆改出失速，模拟飞行员 1 秒反应时间。当飞机失速特性不佳时，失速后模型可能会先出现剧烈滚转的现象，而不会先出现不可控的下俯运动。CCA23中，23.201(d)条款规定了机翼水平失速:在进入和改出机动时，必须有可能使用正常的操纵手段就能防止大于 15的滚转和偏航4。自由飞模型的姿态角与真实飞机对应关系为1:1，因此可以沿用该标准作为模型失速的一项判断依据。迎角是飞机失速试飞的一个重要参数，通过风洞试验可以得到飞机的失速迎角，而失速一般提前或略晚于达到该失速迎角时发生。飞机失速后，可能会出现异常的机头上仰。此时，通过俯仰角速率和滚转角均无法准确判断失速。因此，在模型自由飞试验中，为保证飞行安全，可以将迎角作为失速判定的一项基本条件，即在迎角达到某个预设值(一般大于失速迎角 3左右)，认为飞机模型已经进入失速，可开始执行失速改出动作，迎角预设超限值可根据前期飞行试验结果进行调整。某通航类飞机模型自由飞失速判断如图 1 所示。图 1模型失速判断流程图3控制律设计3.1纵向控制律以往模型自由飞失速试飞方法是由地面操纵手缓慢拉杆，控制模型缓慢抬头进入失速。依靠操纵手控制进入失速过程中，受操纵人员能力限制和遥控器驾驶杆的操纵特点限制，无法充分地实现精准而缓慢的控制，因此实际拉杆过程一般比理论拉杆过程更加剧烈，导致减速率和过载偏大。为此，本文采用自主拉杆失速控制律6，首先采集模型失速试验前配平时的升降舵值作为初始值 de0，输入为升降舵随时间变化的舵面增量，反馈为俯仰角速率，输出为升降舵偏转指令，模拟地面操纵手缓慢拉杆的过程，实现自主拉杆进入失速。控制律构型如图 2 所示。图 2自主拉杆控制律结构框图3.2横航向控制律在机翼水平失速的过程中，通过副翼控制滚转角符合 23部条款中不超过 15的要求。该控制律构型以滚转角指令04作为输入信号，引入滚转角 和滚转角速率 作为反馈信号，解算得到副翼输出指令，实现了滚转角的保持功能。于此同时，在方向舵回路引入了偏航角速率，用以增大飞机模型的荷兰滚阻尼比以及降低滚转时产生的侧滑角。控制律构型见图 3。图 3滚转角保持回路控制律结构框图4飞行仿真及试飞验证4.1相似参数缩比模型自由飞试验为了能够真实反映全尺寸飞机本体的失速特性，必须满足相似准则关系。缩比模型与全尺寸飞机不仅要满足几何外形相似，同时还要满足弗劳德数 相似和质量分布相似准则，关键参数的相似关系见表 1。表 1缩比模型与全尺寸飞机参数相似关系质量惯量距速度拉力K3KK5KK5KK3K迎角侧滑角姿态角角速率舵面偏角11K1表中:K 为缩比系数，K为模型试验高度和飞机飞试验高度的密度比。4.2试验对象根据某通航类飞机三维数模，按照全尺寸飞机缩比进行模型的制造，并与全尺寸飞机具有几何相似的桨叶、襟翼、升降舵、副翼和方向舵，模型螺旋桨按照动力相似关系进行模拟。缩比模型的基本参数为机身长 L=2.5565m，质量 M=50kg，翼展 b=2.7705m，机翼面积 S=0.9287m2，平均气动弦长 MAC=0.3455m。4.3飞行仿真依据某型飞机缩比模型气动数据，由力方程组、力矩方程组、运动方程组和导航方程组十二个状态方程，建立了高精度六自由度非线性仿真模型。1)力方程组:V=(T cos cos D+Gxa)/m?=(T sin L+Gza)/mVcos+(p cossin+q cos r sin sin)/cos=(T cos sin+Y+Gya)/mV+psin r cos2)力矩方程组?p=(c1r+c2p)q+c3?L+c4N?q=c5pr c6(p2 r2)+c7M?r=(c8r c2r)q+c4?L+c9N3)运动方程组=p+(r cos+q sin)tan=q cos r sin=(r cos+q sin)/cos 4)导航方程组?xg=u cos cos+v(sin sin cos cos sin)+w(sin sin+cos sin cos)?yg=u cos sin+v(sin sin sin+cos cos)+w(sin cos+cos sin sin)h=u sin v sin cos w cos cos(5)使用 Matlab 2015a 建立仿真模型，电脑处理器为 Intel()Core(TM)i56500CPU。基于建立的六自由度非线性仿真模型进行纵向和横航向控制律设计，并通过反复迭代确定控制律参数。同时，通过飞行仿真充分验证了俯仰角速率判断缩比模型失速方法的可行性。4.4试飞流程某通航类飞机模型自由飞试验流程如下7 1)飞行前，预设自动开伞高度、滚转角超限值、迎角超限值、失速改出俯仰角速率值;2)通过载机将缩比模型运载升空，到达指定的拉距空域后，缩比模型上电，地面操纵手操纵完成空中拉距试验，确认遥控遥测系统和各个舵面出舵正常;3)载机到达指定的投放空域后，地面操纵手纵杆推杆指令俯仰角为20，同时侧杆中立指令滚转角为 0，在模型与载机安全分离后，纵杆回中指令俯仰角为 0。4)待模型配平后，地面操纵手切换到失速模态，此时飞行控制系统纵向自主拉杆进入失