Modeling of axial compressor characteristic data based on symbolic regression
-
摘要:
为完整准确地描述不同压比轴流压气机的全转速范围的特性数据规律,在以前研究的基础上,利用符号回归机器学习方法对各种压比压气机峰值效率点参数之间的规律进行了进一步地探索,并提出了使用符号回归对多数据源多参数复杂模型的建模方法。使用PySR符号回归库,通过数据逐条拟合,多轮筛选验证,从47份压气机特性数据中,总结出了压气机峰值效率点参数之间的具有一定普适性的规律:总压比-转速、流量-转速之间具有双曲正切函数规律,而效率-转速具有分段二次及指数关系;利用构造的峰值效率点关系式及以前研究中提出的特性线近喘振区和近堵塞区的关联式构造了压气机特性非线性模型,可以对大部分压气机几乎全转速范围特性进行预测,且具有较高预测精度。
Abstract:To comprehensively and accurately characterize the performance characteristics of axial compressors with a variety of pressure ratios across the full rotational speed range, based on previous research, this study further investigated the parametric relationships at the peak efficiency points of compressors with all pressure-ratio using a symbolic regression which is a machine learning approach. A modeling framework leveraging symbolic regression for complex, multi-source, and multi-parameter systems is proposed. Using the PySR symbolic regression library, data points are fitted individually and subjected to multi-round screening and validation. From the characteristic data of 47 compressors, general correlations on the peak efficiency points are derived: both the total pressure ratio–rotational speed and mass flow rate–rotational speed relationships conform to a hyperbolic tangent function, whereas the efficiency–rotational speed relationship exhibits a piecewise quadratic and exponential behavior. By integrating these derived correlations on peak efficiency points with previously established equations for the near-surge and near-choke regions of the characteristic curves, a nonlinear compressor performance model was developed. The proposed model demonstrates the capability to predict the performance characteristics of a broad range of compressors over nearly the entire speed range with high fidelity.
-
表 1 流量-转速的表达式
Table 1. Expressions of mass flow with speed
序号 参数数量 表达式 εmse×104 1 3 $ y={c}_{2}+\tanh ({c}_{1}{x}^{{{c}_{3}}}) $ 2.45 2 3 $ y={c}_{2}+\tanh ({x}^{{{c}_{3}}} ({c}_{1}+{x}^{3}) ) $ 2.61 3 3 $ y={c}_{2}+\tanh ({x}^{{{c}_{3}}} ({c}_{1}+{x}^{2}) ) $ 2.86 4 4 $ y={c}_{1}\tanh ({c}_{2}{x}^{{{c}_{4}}}) +{c}_{3} $ 1.15 5 4 $ y={c}_{1}\tanh ({c}_{2}{x}^{{{c}_{4}}}+{c}_{3}) $ 1.16 6 4 $ y=\tanh (x ({c}_{1}{x}^{{{c}_{4}}}+{c}_{2}) ) +{c}_{3} $ 1.61 表 2 总压比-转速的表达式
Table 2. Expressions of total pressure ratio with speed
序号 参数数量 表达式 εmse×104 1 4 $ y={b}_{1}\tanh ({b}_{2}{x}^{{{b}_{4}}}) +{b}_{3} $ 1.21 2 4 $ \begin{gathered}y ={b}_{1}\exp (-{B}^{-1}) ,\\B ={b}_{2}\exp (\exp ({b}_{3}x) ) +{b}_{4}\end{gathered} $ 1.22 3 4 $ y={b}_{3}+ ({b}_{1}+1/x) \tanh ({b}_{2}{x}^{{{b}_{4}}}) $ 1.83 4 3 $ y=\tanh ({b}_{1}{x}^{{{b}_{3}}}) +{b}_{2} $ 3.52 5 3 $ y=\tanh ({b}_{2}{x}^{{{b}_{3}}}) /x+{b}_{1} $ 4.44 表 3 效率-转速的表达式
Table 3. Expressions of efficiency with speed
分段 表达式 εmse/105 左 $ y=\exp ({c}_{1}{x}^{2} ({c}_{2}+x+{ (x-{{c}_{3}}) }^{-1}) ) $ 4.0 左 $ y=-{c}_{1}{x}^{2}{ ({{c}_{2}}x+1) }^{2}+1 $ 5.3 左 $ y=-{c}_{1}{\mathrm{x}}^{2}\exp ({c}_{2}x) +1 $ 5.4 左 $ y=1- ({x}^{2}+x) ({x}^{2}+3x) ({c}_{1}x+{c}_{2}) $ 9.9 左 $ y=1+ ({x}^{3}+{x}^{2}) (-{c}_{1}x+{c}_{2}) $ 1.0 右 $ y=1-{c}_{1}{x}^{{{c}_{2}}} $ 7.0 右 $ y={c}_{1}{x}^{2} (-{c}_{2}x-{c}_{3}+\exp (\mathrm{x}) ) +1 $ 7.4 右 $ y=\exp (-{c}_{1}{x}^{{{c}_{2}}}) $ 7.6 右 $ y=1-x ({c}_{1}{x}^{{{c}_{3}}}-{c}_{2}x) $ 8.1 右 $ y=1-{a}_{1}{x}^{2} ({a}_{3}+x (2x-{a}_{2}) ) $ 8.1 -
[1] 吴施志, 石建成, 蔡建兵. 燃气涡轮轴和螺桨发动机总体设计[M]. 北京: 科学出版社, 2022. [2] 黄金泉, 周文祥, 潘慕绚, 等. 航空发动机建模、控制与健康管理[M]. 北京: 北京航空航天大学出版社, 2024. Huang Jinquan, Zhou Wenxiang, Pan Muxuan, et al. Aero-engine modelling, control and health management[M]. Beijing: Beihang University Press, 2024. (in ChineseHuang Jinquan, Zhou Wenxiang, Pan Muxuan, et al. Aero-engine modelling, control and health management[M]. Beijing: Beihang University Press, 2024. (in Chinese) [3] Agrawal R K, Yunis M. A generalized mathematical model to estimate gas turbine starting characteristics[J]. Journal of Engineering for Power, 1982, 104(1): 194-201. [4] 张娜, 林汝谋, 蔡睿贤. 压气机特性通用数学表达式[J]. 工程热物理学报, 1996(1): 21-24. Zhang Na, Lin Rumou, Cai Ruixian. General formulas for axial compressor performance estimation[J]. Journal of Engineering Thermophysics, 1996(1): 21-24. (in ChineseZhang Na, Lin Rumou, Cai Ruixian. General formulas for axial compressor performance estimation[J]. Journal of Engineering Thermophysics, 1996(1): 21-24. (in Chinese) [5] Kurzke J. Correlations hidden in compressor maps: GT2011-45519[R]. British Columbia, Canada: Proceedings of ASME Turbo Expo 2011, 2011. [6] Kurzke J. Generating compressor maps to simulate starting and windmilling: ISABE-2019-24299[R]. Canberra Australia: 24th ISABE Conference, 2019. [7] Kurzke J. How to get component maps for aircraft gas turbine performance calculations: 96-GT-164[R] . Birmingham UK: International Gas Turbine and Aeroengine Congress and Exhibition, 1996. [8] Kurzke J. Performance modeling methodology: efficiency definitions for cooled single and multistage turbines: GT-2002-30497[R]. Amsterdam Netherlands: ASME Turbo Expo: Power for Land, Sea, & Air, 2002. [9] Riegler C, Bauer M, Kurzke J. Some aspects of modeling compressor behavior in gas turbine performance calculations[J]. Journal of Turbomachinery, 2001, 123(2): 372-378. doi: 10.1115/1.1368123 [10] 王广, 李军, 牟牧. 航空发动机叶轮机部件低状态特性的指数预测方法[J]. 航空动力学报, 2005, 20(6): 1032-1036. Wang Guang, Li Jun, (mou/mu) Mu. Exponent extrapolation method of aero-engine turbomahinery low state characteristics[J]. Journal of Aerospace Power, 2005, 20(6): 1032-1036. (in Chinese doi: 10.3969/j.issn.1000-8055.2005.06.023Wang Guang, Li Jun, (mou/mu) Mu. Exponent extrapolation method of aero-engine turbomahinery low state characteristics[J]. Journal of Aerospace Power, 2005, 20(6): 1032-1036. (in Chinese) doi: 10.3969/j.issn.1000-8055.2005.06.023 [11] 刘建勋, 李应红, 张东方. 发动机起动过程部件特性的一种改进描述[J]. 燃气涡轮试验与研究, 2007, 20(1): 13-16, 31. Liu Jianxun, Li Yinghong, Zhang Dongfang. A modified description method for components characteristics of engine starting[J]. Gas Turbine Experiment and Research (Quarterly), 2007, 20(1): 13-16, 31. (in Chinese doi: 10.3969/j.issn.1672-2620.2007.01.003Liu Jianxun, Li Yinghong, Zhang Dongfang. A modified description method for components characteristics of engine starting[J]. Gas Turbine Experiment and Research (Quarterly), 2007, 20(1): 13-16, 31. (in Chinese) doi: 10.3969/j.issn.1672-2620.2007.01.003 [12] 王宇, 刘建勋, 李应红, 等. 基于抛物线外推的压气机低转速特性研究[J]. 航空动力学报, 2009, 24(5): 1136-1142. Wang Yu, Liu Jianxun, Li Yinghong, et al. Research on low speed characteristics of compressor based on parabola extrapolation[J]. Journal of Aerospace Power, 2009, 24(5): 1136-1142. (in ChineseWang Yu, Liu Jianxun, Li Yinghong, et al. Research on low speed characteristics of compressor based on parabola extrapolation[J]. Journal of Aerospace Power, 2009, 24(5): 1136-1142. (in Chinese) [13] Tsoutsanis E, Meskin N, Benammar M, et al. A component map tuning method for performance prediction and diagnostics of gas turbine compressors[J]. Applied Energy, 2014, 135: 572-585. doi: 10.1016/j.apenergy.2014.08.115 [14] Mikhailov A E, Mikhailova A B, Akhmedzyanov D A. New 1-D method for the prediction of axial-flow compressors off-design performance[J]. Procedia Engineering, 2016, 150: 155-160. doi: 10.1016/j.proeng.2016.06.740 [15] 施洋, 屠秋野, 严红明, 等. 一种航空发动机全状态性能模型[J]. 航空动力学报, 2017, 32(2): 373-381. Shi Yang, Tu Qiuye, Yan Hongming, et al. A full states performance model for aero engine[J]. Journal of Aerospace Power, 2017, 32(2): 373-381. (in Chinese doi: 10.13224/j.cnki.jasp.2017.02.015Shi Yang, Tu Qiuye, Yan Hongming, et al. A full states performance model for aero engine[J]. Journal of Aerospace Power, 2017, 32(2): 373-381. (in Chinese) doi: 10.13224/j.cnki.jasp.2017.02.015 [16] Li Xu, Yang Chuanlei, Wang Yinyan, et al. Compressor map regression modelling based on partial least squares[J]. Royal Society Open Science, 2018, 5(8): 172454. doi: 10.1098/rsos.172454 [17] 房友龙, 刘东风, 贺星, 等. 压气机特性曲线精确分步拟合方法研究[J]. 燃气涡轮试验与研究, 2019, 32(2): 21-27. Fang Youlong, Liu Dongfeng, He Xing, et al. Research on the precise step fitting method of compressor characteristic map[J]. Gas Turbine Experiment and Research, 2019, 32(2): 21-27. (in Chinese doi: 10.3969/j.issn.1672-2620.2019.02.005Fang Youlong, Liu Dongfeng, He Xing, et al. Research on the precise step fitting method of compressor characteristic map[J]. Gas Turbine Experiment and Research, 2019, 32(2): 21-27. (in Chinese) doi: 10.3969/j.issn.1672-2620.2019.02.005 [18] 尹大伟, 李本威, 王永华, 等. 基于Kriging方法的航空发动机压气机特性元建模[J]. 航空学报, 2011, 32(1): 99-106. Yin Dawei, Li Benwei, Wang Yonghua, et al. Aeroengine compressor characteristics metamodeling using Kriging method[J]. Acta Aeronautica et Astronautica Sinica, 2011, 32(1): 99-106. (in ChineseYin Dawei, Li Benwei, Wang Yonghua, et al. Aeroengine compressor characteristics metamodeling using Kriging method[J]. Acta Aeronautica et Astronautica Sinica, 2011, 32(1): 99-106. (in Chinese) [19] 涂环, 陈辉. 基于Kriging算法的压气机特性建模[J]. 内燃机学报, 2014, 32(4): 377-383. Tu Huan, Chen Hui. Modeling of compressor characteristics using Kriging method[J]. Transactions of CSICE, 2014, 32(4): 377-383. (in ChineseTu Huan, Chen Hui. Modeling of compressor characteristics using Kriging method[J]. Transactions of CSICE, 2014, 32(4): 377-383. (in Chinese) [20] Fei Jingzhou, Zhao Ningbo, Shi Yong, et al. Compressor performance prediction using a novel feed-forward neural network based on Gaussian kernel function[J]. Advances in Mechanical Engineering, 2016, 8(1): 1-14. doi: 10.1177/1687814016628396 [21] 饶高, 苏三买, 翟向博. 指数外推法和支持向量机相结合的压气机特性扩展方法[J]. 航空动力学报, 2017, 32(3): 749-755. Rao Gao, Su Sanmai, Zhai Xiangbo. Method of compressor characteristic extension combining exponent extrapolation method with support vector machine[J]. Journal of Aerospace Power, 2017, 32(3): 749-755. (in Chinese doi: 10.13224/j.cnki.jasp.2017.03.029Rao Gao, Su Sanmai, Zhai Xiangbo. Method of compressor characteristic extension combining exponent extrapolation method with support vector machine[J]. Journal of Aerospace Power, 2017, 32(3): 749-755. (in Chinese) doi: 10.13224/j.cnki.jasp.2017.03.029 [22] 刘锬韬, 高丽敏, 冯旭栋. 基于数据驱动的轴流压气机特性图规律及预测[J]. 航空动力学报, 2021, 36(5): 1072-1082. Liu Tantao, Gao Limin, Feng Xudong. Laws and prediction of axial compressor performance map based on data-driven[J]. Journal of Aerospace Power, 2021, 36(5): 1072-1082. (in Chinese doi: 10.13224/j.cnki.jasp.2021.05.018Liu Tantao, Gao Limin, Feng Xudong. Laws and prediction of axial compressor performance map based on data-driven[J]. Journal of Aerospace Power, 2021, 36(5): 1072-1082. (in Chinese) doi: 10.13224/j.cnki.jasp.2021.05.018 [23] 宋召运, 王宝潼, 范腾博, 等. 基于正反问题耦合的压气机特性快速预测方法[J]. 航空动力学报, 2022, 37(3): 619-628. Song Zhaoyun, Wang Baotong, Fan Tengbo, et al. Rapid prediction method of compressor characteristics with coupling direct and inverse problems[J]. Journal of Aerospace Power, 2022, 37(3): 619-628. (in Chinese doi: 10.13224/j.cnki.jasp.20210072Song Zhaoyun, Wang Baotong, Fan Tengbo, et al. Rapid prediction method of compressor characteristics with coupling direct and inverse problems[J]. Journal of Aerospace Power, 2022, 37(3): 619-628. (in Chinese) doi: 10.13224/j.cnki.jasp.20210072 [24] 王巍, 李哲, 刘祎阳, 等. 基于改进布谷鸟搜索算法的压气机特性曲线预测[J]. 推进技术, 2025, 46(1): 2402034. Wang Wei, Li Zhe, Liu Yiyang, et al. Compressor characteristic curve prediction based on improved cuckoo search algorithm[J]. Journal of Propulsion Technology, 2025, 46(1): 2402034. (in Chinese doi: 10.13675/j.cnki.tjjs.2402034Wang Wei, Li Zhe, Liu Yiyang, et al. Compressor characteristic curve prediction based on improved cuckoo search algorithm[J]. Journal of Propulsion Technology, 2025, 46(1): 2402034. (in Chinese) doi: 10.13675/j.cnki.tjjs.2402034 [25] 甄满, 董学智, 刘锡阳, 等. 基于物理增强的PSO-SVR压气机低转速特性预测方法[J]. 航空动力学报, 2025, 40(7): 20230303. Zhen Man, Dong Xuezhi, Liu Xiyang, et al. Physically enhanced prediction method of extrapolating low speed characteristics of compressor based on PSO-SVR[J]. Journal of Aerospace Power, 2025, 40(7): 20230303. (in Chinese doi: 10.13224/j.cnki.jasp.20230303Zhen Man, Dong Xuezhi, Liu Xiyang, et al. Physically enhanced prediction method of extrapolating low speed characteristics of compressor based on PSO-SVR[J]. Journal of Aerospace Power, 2025, 40(7): 20230303. (in Chinese) doi: 10.13224/j.cnki.jasp.20230303 [26] Cranmer M. Interpretable machine learning for science with PySR and SymbolicRegression. jl[PP/OL]. V3. (2023-05-05) [2026-05-15]. https://doi.org/10.48550/arXiv.2305.01582. [27] Cline S J, Fesler W, Liu H S, et al. High pressure compressor component performance report: NASA-CR-168245[R]. Cleveland, National Aeronautics and Space Administration, 1983. [28] Converse L, Giffin G. Extended parametric representation of compressor fans and turbines Volume I: CMGEN user's manual: NASA-CR-174645[R]. Cleveland, National Aeronautics and Space Administration, 1984. [29] Cornell W G. Experimental quiet engine program: NASA-CR-2519[R]. Cleveland, National Aeronautics and Space Administration, 1975. [30] Love C M, Box P. Design and development of an advanced F100 compressor: AGARD-CP-421[R]. Paris, NATO, 1987. [31] Converse L. Extended parametric representation of compressor fans and turbines volume Ⅲ - MODFAN user’s manual: NASA-CR-174647[R]. Cleveland: National Aeronautics and Space Administration, 1985. [32] Michael J, HALLE E. Energy efficient engine low-pressure compressor component test hardware detailed design report: NASA-CR-165354[R]. Cleveland, National Aeronautics and Space Administration, 1983. [33] Geiselhart K A, Caddy M J, Morris S J Jr. Computer program for estimating performance of air-breathing aircraft engines: NASA TM-4254[R]. Washington DC, National Aeronautics and Space Administration, 1991. [34] Doyle V L, Koch C C. Evaluation of range and distortion tolerance for high Mach number transonic stages: NASA-CR-72720[R]. Cleveland: National Aeronautics and Space Administration, 1970. [35] Flynn J T, Keenan M J, Sulam D H. Single-stage evaluation of highly-loaded high-Mach-number compressor stages. 2-Data and performance, multiple-circular-arc rotor: NASA-CR-72694[R]. Cleveland: National Aeronautics and Space Administration, 1970. [36] Seyler D R, Smith L H Jr. Single stage experimental evaluation of high Mach number compressor rotor blading. Part Ⅰ - Design of rotor blading: NASA CR-54581[R]. Cleveland: National Aeronautics and Space Administration, 1967. [37] Davis J V, Dellert E J. Advancement of small gas turbine component technology, advanced small axial compressor. volume II. addendum: USAAVLABS Technical Report 69-10B[R]. Detroit: Army Aviation Materiel Laboratories (USAAVLABS), 1970. [38] Michael C J, Halle J E. Energy efficient engine low-pressure compressor component test hardware detailed design report: NASA CR-165354[R]. Cleveland: National Aeronautics and Space Administration, 1981. [39] Sinnette J T Jr, Schey O W, King J A. Performance of NACA eight-stage axial-flow compressor designed on the basis of airfoil theory: NACA-TR-758[R]. Cleveland: National Advisory Committee for Aeronautics, 1943. [40] Lippett D, Woollatt G, Ivey P C, et al. The design, development and evaluation of 3D aerofoils for high speed axial compressors: part 1: test facility, instrumentation and probe traverse mechanism: ASME Paper GT2005-68792[R]. Reno: ASME, 2005: 291-302. [41] Tesch W, Pase R. Design and performance of a low aspect ratio, high tip speed multi-stage compressor: AIAA-1983-1161[R]. Seattle: AIAA, 1983. [42] Suzuki M, Kuno N. Research development of two-stage fan component in HYPR project: AlAA-95-2344[R]. San Diego, AIAA, 1995. [43] Dong Y, Mansour M, Hingorani S, et al. A new multistage axial compressor designed with the APNASA multistage CFD code: part 2-application to a new compressor design: ASME Paper 2001-GT-0350[R]. New Orleans : ASME, 2001. [44] Klapproth J, Miller M, Parker D. Aerodynamic development and performance of the CF6-6/LM2500 compressor: AIAA 1979-7030[R]. Orlando: AIAA, 1979. [45] Koff B L. F100-PW-229 higher thrust in same frame size[J]. Journal of Engineering for Gas Turbines and Power, 1989, 111(2): 187-192. doi: 10.1115/1.3240235 [46] Cumpsty N A. Jet Propulsion: A Simple Guide to the Aerodynamic and Thermodynamic Design and Performance of Jet Engines[M]. 2nd ed. Cambridge: Cambridge University Press, 2003. [47] Creason T, Baghdadi S. Design and test of a low aspect ratio fan stage: AIAA-1988-2816[R]. Boston: AIAA, 1988. [48] de la Rosa Blanco E, Hall C, Crichton D. Challenges in the silent aircraft engine design: AIAA 2007-454[R]. Reno: AIAA, 2007. [49] Neubert R J, Gendrich C P. HSCT forward swept fan performance: ASME Paper GT2003-38948[R]. Atlanta: ASME, 2003. [50] Otsuki Y, Nishihara Y, Hoshino E, et al. Kawasaki gas turbine engines for generator sets: ASME Paper 82-GT-6[R]. London: ASME, 1982. [51] Stevenson J D, Saravanamuttoo H I H. Simulating indirect thrust measurement methods for high bypass turbofans: ASME Paper 93-GT-335[R]. Cincinnati: ASME, 1993. [52] Yu Qing, Li Qiushi, Li Ling. The experimental researches on improving operating stability of a single stage transonic fan: ASME Paper GT2002-30640[R]. Amsterdam: ASME, 2002. [53] Veres J. Axial and centrifugal compressor mean line flow analysis method: AIAA-2009-1641[R]. Orlando: AIAA, 2009. [54] Wiedermann A, Frank D, Orth U, et al. Computational and experimental analysis of an industrial gas turbine compressor: ASME Paper GT2011-46336[R]. Vancouver: ASME, 2011. [55] Carchedi F, Wood G R. Design and development of a 12: 1 pressure ratio compressor for the ruston 6-MW gas turbine: ASME Paper 82-GT-20[R]. London: ASME, 1982. -

下载: