摘要
为分析倾转翼eVTOL在考虑电系统热管理对巡航性能的影响,以及不同工况下性能参数的选择,根据飞机二维动力学、阻力、能耗和电系统热管理计算模型,以飞机巡航速度作为控制变量,以巡航时间和电能消耗作为目标函数,建立倾转翼eVTOL巡航仿真模型。结合实际飞行条件设定不同权重,根据模型的具体特征,基于改进粒子群算法进行优化求解,以典型倾转翼eVTOL飞机Vahana为例,得出性能参数并对比分析。仿真结果表明:模型在各种权重组合下,对减少飞机巡航时间和降低电能消耗均有良好的优化效果;最快到达与最节能到达的飞行工况相比,可以节省约50%的时间,但同时会多消耗120%的电耗;与既考虑时间又考虑能耗的折中情况相比,节省了20%时间的同时能会多消耗28%的电能。
飞机推进系统的电气化已被确定为航空业降低碳足迹的潜在解决方案之一。然而,推进系统电气化仍然存在一些挑战,倾转翼eVTOL作为一种创新的低空交通工具,被视为未来城市空中交通(Urban air mobility,UAM)的重要组成部分。eVTOL必须重量轻,且电推进系统部件需提供足够的电力以持续飞行,电推进系统中的发热部件主要包括电池、电机和逆变器,各部件在飞行过程中处于相当大的热应力之下,为了保持它们正常运行,必须去除热量,否则可能会损坏部件,并导致灾难性故障。因此,在关注eVTOL运行参数优化的同时,考虑电系统热限制对其飞行性能的影响是十分重要的,这不仅可以帮助减少事故的发生,还可以降低运营商的运营成本。在未来的研究中,通过对电推进系统内部件参数的优化可能还会为eVTOL的设计改进以及UAM的政策实施和发展带来贡献。
国外学者对考虑电系统的飞行性能做了大量研究工作,Harrison
国内对于飞行性能优化也进行了大量研究,余樾
针对上述问题,首先以典型倾转翼eVTOL飞机为例构建了二维空气动力学和电系统热管理仿真计算模型,然后考虑速度、电系统热限制以及加速度等约束,以巡航时间和电能消耗为优化目标,利用改进粒子群算法对倾转翼eVTOL巡航段进行优化求解,并对不同工况下的巡航性能参数进行对比分析。
倾转翼eVTOL可以在不同飞行阶段通过改变拉力方向,实现起飞、巡航、悬停和降落。当倾转翼eVTOL巡航飞行时,可参照固定翼飞机的飞行状态进行性能分析并做出假设:整个巡航过程都是在标准大气温度下进行的;不考虑倾转旋翼的过渡阶段;不考虑电系统中各部件间电线的热效应;不考虑螺旋桨功率的损失;螺旋桨在产生推力的工作过程中不存在能量损耗。
在巡航状态飞行时,可以认为纵向和横向飞行动力学解耦,通过对姿态和位置控制,使飞机的横向和航向运动保持稳定。为简化问题,在整个巡航过程中不考虑偏航、滚转、侧滑角的变化,因此只分析二维平面内的受力情况,建立运动学模型,如

图1 eVTOL自由体示意图
Fig.1 Free body diagram of eVTOL
因此,倾转翼eVTOL在气动参考系中的运动方程为
(1) |
(2) |
(3) |
(4) |
式中:为巡航时间,为巡航距离;为巡航高度;是由第i个螺旋桨产生的拉力;为螺旋桨个数;为倾转翼eVTOL质量。
倾转翼eVTOL飞机的巡航性能基于传统的二次阻力曲线,使用的阻力模型基于文献[
(5) |
式中:为基于总阻力模型的寄生阻力系数(机翼与机身寄生阻力系数之和);为升力系数;为倾转翼eVTOL等效展弦比;为机翼的总参考面积;为跨度效率;为巡航高度处的空气密度。其中和可表示为
(6) |
(7) |
式中:为翼展;为飞机的最大升力系数;为失速速度。升力可表示为
(8) |
倾转翼eVTOL飞行器在巡航阶段的DEP功耗模型是由标准升力、阻力和拉力模型导出
(9) |
式中为螺旋桨效率。因此,倾转翼eVTOL在巡航阶段消耗的总能量可表示为
(10) |
电动飞机的电推进系统一般由蓄电池、电机、逆变器和螺旋桨构

图2 电推进系统构成示意图
Fig.2 Composition diagram of electric propulsion system
为了计算电池的热效应,采用文献[
(11) |
式中:为戴维南电压瞬时变化率;、和分别为戴维南电压、电阻和比热容;为电池容量;为充电状态;为开路电压;、和分别为导线电流、电压和电阻;、和分别为单个电池的效率、质量和比热容。
电机部件基于螺旋桨所需转速和扭矩来确定电机的效率,电机效率用于确定电机部件上的功率损耗,该功率损耗转换为热量。电机允许空气通过机舱的散热器被动冷却,空气通过散热器后,被输送到电机后面的逆变器,就形成了电机和逆变器之间的循环关系:逆变器为电机提供功率,电机将热量排入气流,气流反过来冷却逆变器。电机的产热率和散热率及电机排出的冷却气流温度可表示为
(12) |
(13) |
式中:为电机输入功率;为电机内部导热系数;为电机对流换热系数;为大气温度;为电机平均温度;为冷却空气质量流率;为气流比热容。因此,电机温度变化率可表示为
(14) |
式中为电机热容量。
针对不同工况的飞行,考虑巡航时电系统中电池和电机发热限制,使巡航时间最短或能量消耗最小,以及二者在不同权重下的性能指标。可以通过控制倾转翼eVTOL巡航速度让飞机在巡航时加速或减速,并利用启发式算法进行优化求解。
根据建立的倾转翼eVTOL巡航飞行模型,在适当的初始状态和必要的约束条件下,建立性能优化模型。鉴于未来低空经济环境,eVTOL运营商通常关注飞行任务的成本、服务水平以及环境影响等方面,即主要考虑运营时间和电能消耗的影响,并根据实际情况调整它们之间的权衡,用来提高飞行任务的经济性和满足乘客需求。因此,将巡航时间和巡航电能消耗作为优化目标函数,对二者赋予不同的权重值相加形成一个新的目标函数,由于两个目标的单位不同,需进行归一化以确保它们可以在相同的尺度上计算。将巡航速度作为控制变量,可得到目标函数
(20) |
式中、为归一化后的巡航时间和电能消耗的函数,表示方法为
(21) |
(22) |
式中:、、、分别表示为倾转翼eVTOL以最小和最大巡航速度飞行时需要的时间和电能消耗。、分别表示巡航时间和电能消耗的权重,二者之和为1。当为1时,即为0,则表示求解巡航时间的最小值;同理,当为0,为1时,求得的结果为巡航时电能消耗的最小值。
粒子群优化算法是一种基于群体智能的启发式优化算法,其在计算效率上的优势使其成为诸多学者关注的重点,并已广泛应用于轨迹和参数优化问题
对于有n个待优化参数的问题,粒子的搜索空间为n维空间,其位置代表该问题1个可能参数取值,记某粒子在第w次迭代后的位置和速度分别为xi和vi,则该粒子在w+1次迭代后的位置为
(23) |
设该粒子历史上达到最佳适应性时对应的位置为xP,该粒子所在的粒子群历史上达到最佳适应性时对应的位置为xG,则粒子的速度为
(24) |
式中:kG和kP为粒子个体评价和群体评价对速度贡献的权重;r1和r2为区间(0,1)上服从均匀分布的随机数。
为提高算法的收敛性,平衡算法在全局最优和局部最优之间搜索的需求,向
(25) |
进一步地,为了使算法具备更好的收敛性,同时提高算法的全局搜索能力,kM、kP和kG均采用基于适应性的动态调节,c11、c12、c21、c22、c31和c32为常数,为第w+1次迭代的适应度函数,、为个体评价和群体评价的适应度函数,即
(26) |
对于搜索空间存在边界的情况,则通过对粒子适应性的调整进行限制。若粒子超出搜索空间所限制的边界,则其适应性将被记为1个不利的惩罚值,同时重新设置kM、kP、kG为常数,从而通过粒子迭代更新机制实现对搜索空间边界的回避。
通过上述的改进粒子群优化算法确定初始状态、约束条件和目标函数再利用MATLAB软件对问题进行求解。具体流程按以下步骤进行优化。
如

图3 改进粒子群算法优化巡航性能流程图
Fig.3 Flowchart of improved particle swarm optimization algorithm for cruise performance
以典型倾转翼eVTOL⁃Vahana为研究对象,飞机采用前后两个可转动的机翼和鸭式布局,可以改变由电力推进驱动的8个螺旋桨的推力方向,其主要参数见
基本参数 | 数值 |
---|---|
起飞质量m/kg | 725 |
机翼总平面面积S/ | 9 |
螺旋桨半径R/m | 0.75 |
平均叶片弦长cb/m | 0.1 |
电池质量Mbatt/kg | 272 |
电机质量Mmotor/kg | 49 |
将倾转翼eVTOL实际可飞行的数据代入电推进系统热管理模型的电机和逆变器模块中,假设巡航高度为1 000 ft(305 m),巡航速度54 m/s,电机初始温度为90 ℃,逆变器初始温度为45 ℃,电机效率为95%,可以得到电机和逆变器模块的热特性如

图4 电机和逆变器模块热特性
Fig.4 Thermal characteristics of motor and inverter modules
倾转翼eVTOL在经历起飞爬升后,螺旋桨和机翼由直升机模式完全过渡到固定翼模式,进入巡航平飞阶段。假设倾转翼eVTOL电池初始温度为50 ℃,电机初始温度为90 ℃。其他参数设置为:失速速度为35 m/s,最大巡航速度为80 m/s,最大加速度为0.5 m/
根据参考倾转翼eVTOL数据和巡航参数初始值,代入2.3节的基于改进粒子群算法的巡航性能模型进行优化求解。分别列举在不同权重、下归一化的巡航时间和巡航电能消耗函数,并进行对比分析。不同工况下的粒子群算法进化图如

图5 改进粒子群算法进化图
Fig.5 Evolutionary graph of improved PSO algorithm
倾转翼eVTOL巡航性能参数和电推进系统温度变化曲线如


图6 参考倾转翼eVTOL巡航性能参数对比
Fig.6 Comparison of cruise performance parameters for reference tilt⁃wing eVTOL
为了确保优化结果的正确性,将
基于改进粒子群算法的 eVTOL巡航性能优化 | 考虑RTA约束的倾转翼 eVTOL能量最优速度剖 | ||
---|---|---|---|
巡航速度Vcurise/(m· | 巡航时间tcurise/s | 对地速度V/(m· | 飞行时间t/s |
45.5 | 1 081 | 45.5 | 约1 100 |
— | — | 50 | 约1 000 |
65 | 760 | 60 | 约850 |
— | — | 70 | 约710 |
80 | 630 | 80 | 约600 |
由
不同工况下的具体的目标函数值、巡航时间和电能消耗值如
权重、 | 目标函数值 | 巡航时间/s | 电能消耗 /(kW·h) |
---|---|---|---|
0f1 1f2 | 0.608 | 1 081 | 14.9 |
0.2f1 0.8f2 | 0.677 | 1 066 | 15.0 |
0.5f1 0.5f2 | 0.626 | 760 | 23.7 |
0.8f1 0.2f2 | 0.347 | 630 | 33.1 |
1f1 0f2 | 0.010 | 630 | 33.1 |
由
针对考虑电系统热限制的倾转翼eVTOL巡航性能优化问题,建立了二维巡航飞行和电系统热管理的仿真模型,并提出了一种改进粒子群优化模型,针对不同目的飞行工况,对巡航过程中的时间和电能消耗的权衡进行优化计算分析,研究表明:
(1)倾转翼eVTOL巡航在考虑电系统热管理的条件下,采用改进粒子群算法的性能参数优化是可行的。为兼顾不同飞行工况目标的飞机性能,可选用不同的权重进行优化研究。
(2)电系统的温度主要与飞机功率有关,倾转翼eVTOL在确定的巡航高度和距离飞行时,速度越大会导致越多的电能消耗,速度越小能耗越小,但需要更长时间的飞行。
(3)在权重组合为1和0的条件下,飞机将加速到最大巡航速度80 m/s飞行到目标距离;在权重组合为0.5和0.5的条件下,飞机将先以65 m/s左右的巡航速度保持600 s左右后加速飞行到目标距离;在权重组合为0和1的条件下,飞机将以最小巡航速度45 m/s保持840 s左右后加速飞行到目标距离。
因此,通过基于改进粒子群算法的倾转翼eVTOL巡航性能优化,可以为未来eVTOL运营商对飞机的飞行性能参数选择提供一定的帮助,从而提高飞行任务的经济性、效率和可持续性,从而创造更大的价值和竞争优势。
参考文献
HARRISOn J, CHARLE D, ZENKER J, et al. Using multi-physics system simulation to predict battery pack thermal performance and risk of thermal runaway during eVTOL aircraft operations[C]//Proceedings of 2019 AIAA/IEEE Electric Aircraft Technologies Symposium (EATS). [S.l.]: IEEE, 2019: 1-13. [百度学术]
FALCK R D, Chin J C, SCHNULO S L, et al. Trajectory optimization of electric aircraft subject to subsystem thermal constraints[C]//Proceedings of the 18th AIAA/ISSMO Multidisciplinary Analysis and Optimization Conference. [S.l.]: AIAA, 2017. [百度学术]
QIAO X, CHEN G, LIN W, et al. The impact of battery performance on urban air mobility operations[J]. Aerospace, 2023, 10(7): 631. [百度学术]
PRADEEP P, WEI P. Energy optimal speed profile for arrival of tandem tilt-wing eVTOL aircraft with RTA constraint[C]//Proceedings of 2018 IEEE CSAA Guidance, Navigation and Control Conference (GNCC). [S.l.]: IEEE, 2018. [百度学术]
SáNCHEZ C N, SANCHEZ J C, RUIZ M Á V, et al. BADA eVTOL performance model for UTM traffic simulation and analysis[C]//Proceedings of the 11th SESAR Innovation Days. Brussels, Belgium: SESAR, 2021: 91. [百度学术]
WARREN M, GARBO A, HERNICZEK M K, et al. Effects of range requirements and battery technology on electric VTOL sizing and operational performance[C]//Proceedings of AIAA SciTech 2019 forum. [S.l.]: AIAA, 2019: 527. [百度学术]
余樾,牛欣怡,索晓杰.基于遗传算法的飞机垂直剖面优化研究[J].电脑编程技巧与维护,2023(5):132-134,173. [百度学术]
YU Yue, NIU Xinyi, SUO Xiaojie. Research on optimization of aircraft vertical profile based on genetic algorithm[J]. Computer Programming Skills & Maintenance, 2023(5): 132-134,173. [百度学术]
魏志强,张文秀,韩博.考虑飞机排放因素的飞机巡航性能参数优化方法[J].航空学报,2016,37(11):3485-3493. [百度学术]
WEI Zhiqiang, ZHANG Wenxiu, HAN Bo. Optimization method of aircraft cruise performance parameters considering pollution emissions[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(11): 3485-3493. [百度学术]
回忆,魏志强,魏路欢.基于飞机性能影响的路径优化算法[J].南京航空航天大学学报,2022,54(6): 1145-1151. [百度学术]
HUI Yi, WEI Zhiqiang, WEI Luhuan. A trajectory optimization algorithm based on impact of aircraft performance[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2022, 54(6): 1145-1151. [百度学术]
LOVERING Z. Vahana configuration trade study part I[EB/OL]. (2016-12-14)[2023-10-01]. https://acubed.airbus.com/blog/vahana/vahana-configuration-trade-study-part-i/. [百度学术]
KLESH A T, KABAMBA P T. Solar-powered aircraft: Energy-optimal path planning and perpetual endurance[J]. Journal of Guidance, Control, and Dynamics, 2009, 32(4): 1320-1329. [百度学术]
王书礼,孙金博,康桂文,等.一种电动飞机电推进系统的能效优化方法[J].航空学报,2021,42(3): 53-61. [百度学术]
WANG Shuli, SUN Jinbo, KANG Guiwen, et al. Energy efficiency optimization method for electric aircraft propulsion system[J]. Acta Aeronautica et Astronautica Sinica, 2021, 42(3): 53-61. [百度学术]
CHIN J, SCHNULO S L, MILLER T, et al. Battery performance modeling on sceptor X-57 subject to thermal and transient considerations[C]//Proceedings of AIAA Scitech 2019 Forum. [S.l.]: AIAA, 2019: 784. [百度学术]
罗鸣,左林玄.基于改进粒子群算法的高超声速飞机爬升轨迹优化[J].飞机设计,2023,43(4): 7-11. [百度学术]
LUO Ming, ZUO Linxuan. Climb trajectory optimization of hypersonic aircraft based on improved particle optimization[J]. Aircraft Design, 2023, 43(4): 7-11. [百度学术]
LEE B S, YUN J Y, HWANG H Y. Flight range and time analysis for classification of eVTOL PAV[J]. Journal of Advanced Navigation Technology, 2020, 24(2): 73-84. [百度学术]