![点击分享此内容可以赚币 分享](/master/images/share_but.png)
一种组合动力飞行器模态转换过程轨迹优化与控制方案.pdf
《一种组合动力飞行器模态转换过程轨迹优化与控制方案.pdf》由会员分享,可在线阅读,更多相关《一种组合动力飞行器模态转换过程轨迹优化与控制方案.pdf(9页珍藏版)》请在咨信网上搜索。
1、第 45 卷第 3 期 2024 年 3 月宇 航 学 报Journal of AstronauticsNo.32024MarchVol.45一种组合动力飞行器模态转换过程轨迹优化与控制方案刘凯1,张永亮2,聂聆聪2(1.大连理工大学力学与航空航天学院,大连 116024;2.北京动力机械研究所,北京 100074)摘要:为应对组合动力飞行器涡轮/冲压发动机模态转换过程中容易出现的推力陷阱问题,开展了组合动力飞行器模态转换阶段飞行/推进一体化保护控制研究。该研究利用组合动力爬升段飞行轨迹优化方法,解决模态转换过程中组合发动机总推力无法满足平飞加速需求的推力陷阱问题。结合轨迹线性化控制方法,完成
2、了宽速域飞行器/发动机一体化轨迹跟踪控制设计,利用迎角、油门等变量的协同调节实现轨迹跟踪控制。仿真结果表明,轨迹优化策略能够在一定程度上克服模态转换过程中的推力不足问题,一体化轨迹线性化控制方法可以有效实现爬升轨迹跟踪,避免模态转换过程对组合动力飞行器飞行任务造成的不利影响。关键词:宽速域飞行器;涡轮/冲压发动机;飞行/推进一体化控制;轨迹优化;轨迹跟踪控制中图分类号:V249.122+.2 文献标识码:A 文章编号:1000-1328(2024)03-0443-09 DOI:10.3873/j.issn.1000-1328.2024.03.011A Trajectory Optimizati
3、on and Control Scheme for Mode Conversion Process of Turbine-based Combined Cycle VehicleLIU Kai1,ZHANG Yongliang2,NIE Lingcong2(1.College of Mechanics and Aerospace,Dalian University of Technology,Dalian 116024,China;2.Beijing Power Machinery Institute,Beijing 100074,China)Abstract:To deal with the
4、 problem of thrust trap during the process of turbine/ramjet mode conversion of the turbine based combined cycle(TBCC),the integrated control of turbine/scramjet mode conversion has been studied.The combined power climbing flight trajectory optimization method is used to solve the thrust trap proble
5、m that the total thrust of the combined engine cannot meet the acceleration demand of horizontal flight during the modal conversion process.The control variables and state variables are discretized to complete the trajectory optimization.The integrated control strategy is designed to complete the jo
6、int adjustment of the angle of attack and the throttle during the mode conversion.The simulation results show that the trajectory optimization method can overcome the thrust trap through a special trajectory design.The integrated protection control method improves the control quality and avoids the
7、adverse effects of the modal conversion process.Key words:Wide-speed-range vehicle;Turbine engine/Ramjet;Integrated flight/propulsion control;Trajectory optimization;Trajectory tracking control0引言为了适应低成本天地往返运载需求,可重复使用的自主水平起降宽速域飞行器成为各国研究热点。因其宽速域特点,传统涡扇发动机不能满足高马赫数下的推力需求,冲压发动机无法在低马赫数下稳定工作。涡轮基组合循环(Turbi
8、ne-based combine cycle,TBCC)发动机方案通过组合涡扇与收稿日期:2023-09-08;修回日期:2024-01-26宇航学报第 45 卷冲压动力系统,兼顾飞行器宽速域内推进性能需求,成为设计速域在Ma 07内的宽速域飞行器理想动力方案1-3。宽速域飞行器起飞爬升加速过程中,亚声速与低马赫状态涡扇单独工作,当飞行器临近涡扇发动机工作包线边界时,涡扇/冲压共同工作(以下称模态转换),最后达到冲压单独工作状态。其中模态转换阶段需要调整进气道分流板与尾喷管,将一部分进气道压缩流从低速涡扇流道分流到高速冲压流道并产生推力4-5。此时冲压流道内部燃烧不稳定,对来流状态与燃油当量比
9、的要求和限制十分严格,极易发生进气道不起动5;此时宽速域飞行器的飞行速度接近涡扇系统的设计极限,又尚未达到冲压系统的有效工作区间,导致推进系统提供的总推力不足,形成“推力陷阱”现象,给飞行控制带来新的挑战。因此,对涡轮基组合循环发动机方案模态转换过程的相关研究有重要的工程实践意义。近年来,一些学者对TBCC模态转换过程与推力陷阱问题进行了研究。文献 7 介绍了包括变循环方案、预冷方案和火箭助力方案在内的3种提升涡扇发动机工作上限的发动机设计方法。文献 8-9分别提出了多催化剂接力催化的自由基接力燃烧(FRRC)策略来大幅降低煤油点火温度和优化进气道回流区附面层抽吸设计,从而使冲压发动机实现低马
10、赫数起动。文献 10 采用动网技术分析攻角对外二元并联TBCC模态转换期间进气道出口气流状态的影响。文献 11 以进气道流量分配比例变化为条件,分析了模态转换段TBCC性能参数的变化规律,阐述了发动机性能衰减的主要原因。目前国内外学者关于模态转换与推力陷阱的研究集中于产生机理分析,提升涡扇发动机工作上界,降低冲压发动机工作下界以及飞行状态对发动机性能影响分析等方面,鲜有从轨迹优化与飞行/推进一体化控制的角度给出模态转换阶段的解决方案。宽速域飞行器采用飞推一体化构型设计,单独对推进系统或飞行状态进行解耦分析设计无法系统性地克服推力陷阱。因此,本文针对一型TBCC宽速域飞行器进行飞行动力学建模、定
11、量分析模态转换过程的推进系统能力与推力需求,通过轨迹设计克服模态转换阶段的推力陷阱,采用轨迹线性化制导方法实现飞行/推进一体化控制,为解决宽速域飞行器推力陷阱问题提供一条新的技术途径。1宽速域飞行器动力学建模1.1宽速域飞行器质点动力学模型在飞行器航迹坐标系下,依据经典牛顿定律,飞行器的质心运动可以用如式(1)的矢量方程形式表述:mdVdt=m(Vt+V)=F+(1)式中,m是飞行器质量,V为飞行速度,是转动角速度,F是飞行器合外力矢量,T是推力矢量。忽略地球自转影响,根据飞行器航迹坐标系与地面坐标系之间的转换关系,整理可得飞行器三变量质心运动的动力学方程如式(2)所示:dVdt=Tcos-D
12、m-gsinddt=Tsincos+LcosmV-gcosVddt=-Tsinsin+LsinmVcosdhdt=VsinL=fL()h,z,e,VD=fD()h,z,e,r,VT=fT()h,V(2)式中:,h分别表示飞行器航迹倾角、飞行器航迹偏角、飞行高度;D,L为阻力、升力;,为攻角、倾侧角;z,e,r分别为升降舵偏角、副翼偏角、方向舵偏角。1.2模态转换阶段轨迹策略分析本文基于激波理论构建二维进气道模型,得到进气道出口的流场参数。采用 T-mats工具包建立涡扇发动机模型,依据质量守恒、动量守恒和能量守恒,建立准一维冲压发动机模型。采用面元法分析宽速域飞行器气动特性12-13。对于经过
13、优化设计后的TBCC模态转换过程,转换速域能够降低至Ma 2.12.6。此时涡扇发动机效能不能满足推力需求。另一方面,飞行速度低导致冲压发动机进气道内无法形成足够强的激波系从而对自由流进行有效压缩,抗反压能力很弱。燃油当量比受限于进气道不起动保护约束和贫油熄火约束之间,导致其值和调节范围都很小。冲压发动机也不能提供足够的推力满足飞行器推力需求。综上所述,模态转换阶段飞行器的核心需求是444第 3 期刘凯等:一种组合动力飞行器模态转换过程轨迹优化与控制方案加速,使飞行马赫数达到冲压发动机的正常工作区间,完成模态转换。首先对模态转换区间内的TBCC性能与宽速域飞行器推力需求进行分析,对比飞行器平飞
14、加速最小推力需求与TBCC最大推力能力。宽速域飞行器在设定飞行状态的配平攻角如图1所示。在得到宽速域飞行器配平攻角的基础上,可以求解得到飞行器在设定飞行状态下的阻力。结合涡扇与冲压发动机性能,可以对比阻力与总推进能力。图1为考虑宽速域飞行器飞推耦合效应的气动特性情况下,配平攻角与飞行状态之间的关系,图2为配平攻角下推力需求与推进系统能够提供最大推力的对比。可以看出,宽速域飞行器组合动力发动机在Ma 2.252.55速域内,即冲压发动机尚未达到有效工作区间的阶段,出现最大推力不满足平飞加速所需最小推力的推力陷阱现象。在此阶段,需要采用特殊轨迹设计,使宽速域飞行器跨越推力陷阱,恢复机动能力,如图3
15、所示。具体轨迹策略是把帮助宽速域飞行器完成模态转换的轨迹设计分为2个阶段:第1阶段是爬升预备阶段,此阶段需要宽速域飞行器利用涡扇发动机达到足够的高度和速度,给模态转换和可能的俯冲加速做动能和势能储备;第2阶段为模态转换阶段,此阶段涡扇发动机和冲压发动机共同工作,利用优化后的飞行轨迹进行进一步加速,达到冲压发动机的正常工作速域。在轨迹规划上,第1阶段设定飞行目标是使飞行器达到涡扇发动机的极限高度和冲压启动速度;第2阶段的主要目标是使飞行器达到冲压发动机正常工作的速度要求,而相应地放松对飞行器的高度限制。2宽速域飞行器飞推一体化控制器设计2.1基于自适应高斯伪谱法的模态转换轨迹优化Gauss伪谱法
16、14-15由于在计算效率上表现出显著的优势广泛应用于求解最优控制问题,同时在航空航天领域常用于轨迹优化。该方法通过对控制量和状态量进行离散,直接进行寻优解算,从而完成非线性规划任务16。考虑一般形式的非线性系统动力学方程为:x(t)=f x(t),u(t),t(3)式中:状态变量x(t)Rn,控制变量u(t)Rm,时间t t0,tf。若最优控制的时间范围是t0,tf,而Gauss伪谱法时间范围是-1,1,则对时间t进行变换可以表示为:=2ttf-t0-tf+t0tf-t0(4)Gauss 伪 谱 法 经 过 离 散 处 理 的 点 为 K 阶1.52.02.53.03.54.0Ma012345
17、67配平攻角与飞行状态关系攻角/()14 km推进性能裕度16 km推进性能裕度18 km推进性能裕度20 km推进性能裕度22 km推进性能裕度图1配平攻角Fig.1The trim angle of attack1.52.02.53.03.54.0Ma21012345678推力裕度与飞行状态关系推力裕度/N10414 km推进性能裕度16 km推进性能裕度18 km推进性能裕度20 km推进性能裕度22 km推进性能裕度0裕度参考线图2推力裕度与飞行状态关系Fig.2Relationship between thrust margin and flight stateMaH/m1.01.5
18、0.52.02.53.03.02.52.01041.51.00.53.54.0ACB(Ma0.9爬升至10 km,平飞跨声速)E1(俯冲完成模态转换)E(平飞完成模态转换)D(发动机模态转换)F常规轨迹优化轨迹图3模态转换阶段的轨迹优化策略Fig.3Trajectory optimization strategy in mode transition445宇航学报第 45 卷Legendre-Gauss 点,等价于 K 阶 Legendre 项式的根k=1,K,。增加额外点 0=-1,作为第K+1个插值点,然后以K+1个 Lagrange 插值多项式Li()(i=0,1,K)为基函数,进行状态
19、变量的近似:x()X()=i=1KLi()x(i)(5)最后一个点 f=1,与之相对应的状态变量x(f)同样符合动力学方程关系。离散末端状态的动力学方程并通过Gauss积分来近似,如式(6):X(f)=X(0)+tf-t02wkfX(k),U(k),;t0,tf(6)式 中:wk=-11Li()d为 Gauss 权 重 系 数,k为Legendre-Gauss点。同样地,用 Lagrange 插值多项式为基函数来近似控制变量,其中U(k),(k=1,2,K)为离散点上的控制变量:i(i=1,2,K)为Legendre-Gauss点。对式(6)求导得到状态变量的导数,可以将动力学方程约束变为代数
20、约束,再反代入动力学方程,则得到状态变量在配点处应受的代数方程约束:i=0KDki(k)X(i)-tf-t02fX(k),U(k),k;tf,t0=0(7)边界条件约束满足:(X0,t0,Xf,tf)=0(8)过程约束满足:C(Xk,Uk,k;t0,tf)0,(k=1,2,K)(9)用Gauss积分近似最优控制问题中性能指标函数的积分项,可以得到性能指标函数:J=F(X0,t0,Xf,tf)+tf-t02k=1Kwkg(Xk,Uk,k;t0,tf)(10)此时的最优控制问题已经被转换为:使得经过离散化后的状态变量Xi、控制变量Uk以及终端时刻tf,在满足动力学方程约束的前提下,求解J函数最小,
21、即式(10)最小。2.2模态转换阶段轨迹策略对比仿真第1阶段宽速域飞行器以爬升为目标,设定第1阶段的状态量和控制量约束如表1所示。采用高斯伪谱法优化且经过积分验证的优化爬升轨迹如图4所示。表1爬升段状态量及控制量约束Table 1State quantity and control quantity constraints in climbing stage物理量速度 V/(ms-1)弹道倾角/()弹道偏角/()高度 H/m质量 m/kg攻角/()涡扇油门开度初始条件4900011 00015 000终端条件5700,12021 0001,1220%,120%01020304050607080
22、90 100 1.01.21.41.61.82.02.2 最少油耗最短时间最少油耗最短时间最少油耗最短时间最少油耗最短时间最少油耗最短时间4805005205405605806001001020304050 1.4801.4851.4901.4951.500 051015高度/m 104104时间/s(a)涡扇爬升段高度曲线0102030405060708090 100 时间/s(b)涡扇爬升段速度曲线0102030405060708090 100 时间/s(c)涡扇爬升段弹道角曲线0102030405060708090 100 时间/s(d)涡扇爬升段质量变化曲线01020304050607
23、08090 100 时间/s(e)涡扇爬升段攻角变化曲线 速度/(ms1)弹道角/()质量/kg 攻角/()446第 3 期刘凯等:一种组合动力飞行器模态转换过程轨迹优化与控制方案由图4优化爬升轨迹可知:在最终高度差小于200 m,最终速度差小于 4 m/s的情况下对比,最短时间路径对比最少油耗路径少用时间39.4 s,多消耗油料66.4 kg。将该阶段的末端飞行状态作为第2阶段(模态转换段)的初始状态,进行第2阶段的轨迹优化。由表2可知模态转换阶段马赫数较小,冲压发动机燃烧不稳定,过大的燃油当量比容易引起进气道不起动和点火困难17-18。因此将燃油当量比限制为常值0.35。采用高斯伪谱法优化
24、且经过积分验证的模态转换阶段轨迹如图5所示。由图5可以看到,传统轨迹在达到710 m/s时,遭遇推力陷阱,推力和推力需求此时达到平衡,无法加速。优化轨迹均可以克服推力陷阱,达到目标速度。最大高度轨迹耗时61.5 s,下降高度1.37 km;最少时间轨迹耗时36.8 s,下降高度5.3 km,比最大高度轨迹少消耗燃油60 kg。最小油耗轨迹比最少时间轨迹多耗时7.2 s;少消耗燃油5.4 kg;高度少下降615.1 m。3宽速域飞行器轨迹线性化控制策略本节基于前文所建的质点动力学模型,对宽速10401020304050607080901.41.61.82.02.2最大高度最少时间最少油耗传统轨迹
25、最大高度最少时间最少油耗传统轨迹最大高度最少时间最少油耗传统轨迹最大高度最少时间最少油耗传统轨迹最大高度最少时间最少油耗传统轨迹最大高度最少时间最少油耗传统轨迹最大高度最少时间最少油耗传统轨迹600700800900 2010010 1.4551.4601.4651.4701.475 202468708090100110120 涡扇油门开度/%0.250.300.35 冲压燃油当量比 高度/m 时间/s(a)共同工作段高度曲线0102030405060708090时间/s(b)共同工作段速度曲线0102030405060708090时间/s(c)共同工作段弹道角曲线01020304050607
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 一种 组合 动力 飞行器 转换 过程 轨迹 优化 控制 方案
![提示](https://www.zixin.com.cn/images/bang_tan.gif)
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。