CK-IBNN surrogate model for aerodynamic performance prediction of multi-stage axial flow compressor
-
摘要:
CoKriging模型在百千维变量问题上的建模代价往往难以承受,以致多级压气机叶片形状优化的百千维问题无法求解。为打破CoKriging模型的维数诅咒问题,建立了一种基于无限宽度贝叶斯神经网络(IBNN)的关联函数及基于IBNN关联函数的CoKriging模型。理论分析显示IBNN关联函数不利用欧式空间距离评价任两样本点之间的关联程度,且关联函数的超参数对任意维度问题均为3个,与变量个数无关。因此,所提方法的建模代价可显著降低。为验证所提方法的有效性和高效性,求解了5维、31维、144维和
1512 维的压气机气动建模问题,并与多保真度深度神经网络、分层Kriging等方法相比,结果显示所提方法可快速建立高精度代理模型,仅需0.1 s和7 s即可完成144维和1512 维问题的压气机气动性能预测模型的建立。-
关键词:
- 代理模型 /
- CoKriging模型 /
- 维数诅咒 /
- 多级压气机 /
- 气动建模
Abstract:The prohibitive computational cost of the model construction of CoKriging for hundreds or thousands of variables prevents the practical application of CoKriging. Therefore, the surrogate-based optimization design of multi-stage axial flow compressor with hundreds or thousands of variables cannot be conducted. To break the curse of dimensionality of the CoKriging model, a correlation function based on infinite-width Bayesian neural network (IBNN) along with CoKriging modeling method utilizing the IBNN correlation function was established. The theoretical analysis demonstrated that the IBNN correlation function was not dependent on the spatial distance to measure the similarity of any two points. Furthermore, the number of hyperparameters of the IBNN function was fixed as 3, being independent of the number of modeling variables like the existing correlation functions. Therefore, the modeling efficiency can be decreased significantly. To verify the effectiveness and efficiency of the proposed method, aerodynamic problems of axial flow compressor with 5, 31, 144 and
1512 variables were solved. Performance of the proposed method was compared with multi-fidelity deep neural network, and Hierarchical Kriging, etc. The proposed method can build the surrogate models efficiently with adequate accuracy. The surrogate models for the aero-dynamic performance prediction with 144 and1512 variables can be tuned within 0.1 and 7 seconds. -
表 1 常见关联函数
Table 1. Representative correlation functions
函数名称 表达式 EXP $ \exp \left(\displaystyle\sum_{k=1}^{d}-\theta_{k} \left| x_{k}-x_{k}^{\prime}\right|\right) $ EXPG $ \exp \left(\displaystyle\sum_{k=1}^{d}-\theta_{k} \left| x_{k}-x_{k}^{\prime}\right|^{\theta_{d+1}}\right)\quad\quad 0<\theta_{d+1}<2 $ Gaussian $ \exp \left(\displaystyle\sum_{k=1}^{d}-\theta_{k}\left| x_{k}-x_{k}^{\prime}\right|^2\right) $ Linear $\displaystyle\prod_{k=1}^{d}\max\{0,1-\theta_k|x_k-x'_k|\} $ Cubic $ \begin{array}{l}\displaystyle\prod_{k=1}^{d} 1-3 \xi_{k}^{2}+2 \xi_{k}^{3} \\\xi_{k}=\min \left\{0, \theta_{k} \mid x_{k}-x_{k}^{\prime}\mid\right\}\end{array} $ Spline $\displaystyle\prod_{k=1}^{d}\zeta (\xi_k) ,\xi_k=\theta_k|x_k-x'_k| $ Matérn 5/2 $ \begin{array}{l}\left(1+\sqrt{5}a+\dfrac{5a^2}{3}\right)\exp\left(-\sqrt{5}a\right) \\a=\sqrt{\displaystyle\sum_{k=1}^d\theta_k\left|x_k-x'_k\right|^2}\end{array} $ 表 2 测试算例信息汇总
Table 2. Information of test problems
$d$ ${N_{\text{l}}}$ $N$ ${N_{\text{t}}}$ 问题名称 5 50 10 10 10StgVGV-$\pi $ 10StgVGV-$\eta $ 10StgVGV-$ \dot m $ 31 296 164 285 Rotor 37-$\eta $ Rotor 37-$\pi $ 144 1168 520 384 3StgDef-$\eta $ 3StgDef-$\pi $ 3StgDef-$ \dot m $ 1512 15120 7560 3024 10StgDef-$\pi $ 10StgDef-$\eta $ 10StgDef-$ \dot m $ 表 3 MFDNN训练参数
Table 3. Training parameters of MFDNN
层数 宽度 最大步数 问题名称 3 40 30000 10StgVGV-$\pi $ 10StgVGV-$\eta $ 10StgVGV-$ \dot m $ 10 100 40000 Rotor 37-$\eta $ Rotor 37-$\pi $ 12 600 40000 3StgDef-$\eta $ 3StgDef-$\pi $ 3StgDef-$ \dot m $ 15 800 40000 10StgDef-$\pi $ 10StgDef-$\eta $ 10StgDef-$ \dot m $ 表 4 R2指标对比
Table 4. Comparison of R2 metric
问题名称 MFDNN K-DIC HK-DIC K-IBNN CK-IBNN 10StgVGV-$\pi $ 0.9647 0.0053 0.9855 0.8684 0.9841 10StgVGV-$\eta $ 0.3624 − 0.1303 0.7684 0.7411 0.8669 10StgVGV-$ \dot m $ 0.4581 − 0.0093 0.9797 0.8594 0.9512 Rotor 37-$\eta $ 0.6179 0.8694 0.9686 0.8457 0.9554 Rotor 37-$\pi $ 0.9131 0.9834 0.9953 0.9710 0.9935 3StgDef-$\eta $ 0.9508 0.9610 0.9736 0.9054 0.9589 3StgDef-$\pi $ 0.8574 0.8715 0.9048 0.8244 0.8789 3StgDef-$ \dot m $ 0.9252 0.9386 0.9575 0.8630 0.9371 10StgDef-$\pi $ 0.8726 0.9404 0.9661 10StgDef-$\eta $ − 9.8144 0.6649 0.6808 10StgDef-$ \dot m $ 0.9118 0.9193 0.9406 表 5 ERMSE指标对比
Table 5. Comparison of ERMSE metric
问题名称 MFDNN K-DIC HK-DIC K-IBNN CK-IBNN 10StgVGV-$\pi $ 0.1784 0.9462 0.1141 0.3441 0.1197 10StgVGV-$\eta $ 0.7575 1.0086 0.4566 0.4827 0.3461 10StgVGV-$ \dot m $ 0.6983 0.9531 0.1351 0.3558 0.2097 Rotor 37-$\eta $ 0.6170 0.3607 0.1769 0.3921 0.2108 Rotor 37-$\pi $ 0.2942 0.1285 0.0681 0.1699 0.0803 3StgDef-$\eta $ 0.2216 0.1971 0.1622 0.3071 0.2024 3StgDef-$\pi $ 0.3771 0.3579 0.3082 0.4185 0.3475 3StgDef-$ \dot m $ 0.2731 0.2475 0.2060 0.3697 0.2505 10StgDef-$\pi $ 0.3569 0.2440 0.1842 10StgDef-$\eta $ 3.2880 0.5788 0.5649 10StgDef-$ \dot m $ 0.2969 0.2840 0.2437 表 6 EMAE指标对比
Table 6. Comparison of EMAE metric
问题名称 MFDNN K-DIC HK-DIC K-IBNN CK-IBNN 10StgVGV-$\pi $ 0.4505 1.8576 0.2100 0.5880 0.2091 10StgVGV-$\eta $ 1.6343 1.4028 1.1866 0.9220 0.7485 10StgVGV-$ \dot m $ 1.8306 1.9491 0.2299 0.6632 0.3987 Rotor 37-$\eta $ 2.6992 1.9322 1.7223 1.9705 2.0445 Rotor 37-$\pi $ 0.9830 0.7303 0.4902 0.7794 0.6312 3StgDef-$\eta $ 0.6965 0.7204 0.5286 1.1033 0.7095 3StgDef-$\pi $ 2.0361 1.9271 1.9778 1.9177 2.0978 3StgDef-$ \dot m $ 1.1077 0.9200 0.8396 1.3683 0.8606 10StgDef-$\pi $ 1.5141 1.1538 1.2180 10StgDef-$\eta $ 14.3835 4.2887 3.8824 10StgDef-$ \dot m $ 1.6787 1.7340 1.9766 表 7 建模时间对比
Table 7. Comparison of modeling time
问题名称 MFDNN K-DIC HK-DIC K-IBNN CK-IBNN 10StgVGV-$\pi $ 234.6345 0.0028 0.0087 0.0003 0.0004 10StgVGV-$\eta $ 259.8951 0.0019 0.0093 0.0004 0.0003 10StgVGV-$ \dot m $ 259.7941 0.0021 0.0097 0.0001 0.0004 Rotor 37-$\eta $ 712.4516 0.0511 0.1825 0.0051 0.0068 Rotor 37-$\pi $ 718.8171 0.0621 0.2174 0.0013 0.0073 3StgDef-$\eta $ 797.9772 2.5133 13.5913 0.0156 0.0956 3StgDef-$\pi $ 793.5154 2.2602 12.9230 0.0120 0.0932 3StgDef-$ \dot m $ 801.3981 2.1465 12.9916 0.0122 0.0953 10StgDef-$\pi $ 6723.5427 3.1202 6.5156 10StgDef-$\eta $ 6466.9472 2.9386 6.5598 10StgDef-$ \dot m $ 6801.4526 3.0376 6.6351 表 8 不同M值的CK-IBNN模型R2值
Table 8. R2 metric value of CK-IBNN model with various M
问题名称 M=1 M=2 M=3 M=4 M=5 10StgVGV-$\pi $ 0.9823 0.9829 0.9841 0.9844 0.9844 10StgVGV-$\eta $ 0.8624 0.8561 0.8669 0.8730 0.8759 10StgVGV-$ \dot m $ 0.9354 0.9470 0.9512 0.9536 0.9554 Rotor 37-$\eta $ 0.9463 0.9531 0.9554 0.9565 0.9571 Rotor 37-$\pi $ 0.9923 0.9933 0.9935 0.9934 0.9932 3StgDef-$\eta $ 0.9560 0.9580 0.9589 0.9592 0.9593 3StgDef-$\pi $ 0.8726 0.8767 0.8789 0.8800 0.8807 3StgDef-$ \dot m $ 0.9328 0.9354 0.9371 0.9381 0.9388 10StgDef-$\pi $ 0.9691 0.9679 0.9661 0.9638 0.9614 10StgDef-$\eta $ 0.6790 0.6813 0.6808 0.6786 0.6756 10StgDef-$ \dot m $ 0.9443 0.9429 0.9406 0.9378 0.9348 表 9 不同M值的CK-IBNN模型EMAE值
Table 9. EMAE metric value of CK-IBNN model with various M
问题名称 M=1 M=2 M=3 M=4 M=5 10StgVGV-$\pi $ 0.2690 0.1986 0.2091 0.2358 0.2570 10StgVGV-$\eta $ 0.7283 0.7683 0.7485 0.7192 0.6884 10StgVGV-$ \dot m $ 0.4798 0.4392 0.3987 0.3630 0.3658 Rotor 37-$\eta $ 1.9332 2.0129 2.0445 2.0585 2.0645 Rotor 37-$\pi $ 0.5648 0.6011 0.6312 0.6562 0.6777 3StgDef-$\eta $ 0.7090 0.7089 0.7095 0.7103 0.7115 3StgDef-$\pi $ 2.1331 2.1117 2.0978 2.0883 2.0817 3StgDef-$ \dot m $ 0.8754 0.8690 0.8606 0.8522 0.8445 10StgDef-$\pi $ 1.3230 1.2661 1.2180 1.1755 1.1375 10StgDef-$\eta $ 3.8006 3.8451 3.8824 3.9146 3.9429 10StgDef-$ \dot m $ 1.9297 1.9557 1.9766 1.9941 2.0093 表 10 不同${\boldsymbol{\sigma}} _{\bf{w}}^{\bf{2}}$值的CK-IBNN模型精度指标值
Table 10. Metric value of CK-IBNN model with various ${\boldsymbol{\sigma}} _{\bf{w}}^{\bf{2}}$
问题名称 R2 MAE $\sigma _{\text{w}}^2$=5 $\sigma _{\text{w}}^2$=10 $\sigma _{\text{w}}^2$=15 $\sigma _{\text{w}}^2$=5 $\sigma _{\text{w}}^2$=10 $\sigma _{\text{w}}^2$=15 10StgVGV-$\pi $ 0.9851 0.9841 0.9834 0.1935 0.2091 0.2146 10StgVGV-$\eta $ 0.8717 0.8669 0.8635 0.7312 0.7485 0.7600 10StgVGV-$ \dot m $ 0.9538 0.9512 0.9494 0.3966 0.3987 0.4022 Rotor 37-$\eta $ 0.9548 0.9554 0.9555 2.0493 2.0445 2.0425 Rotor 37-$\pi $ 0.9936 0.9935 0.9934 0.6268 0.6312 0.6325 3StgDef-$\eta $ 0.9588 0.9589 0.9589 0.7181 0.7095 0.7057 3StgDef-$\pi $ 0.8777 0.8789 0.8792 2.1007 2.0978 2.0969 3StgDef-$ \dot m $ 0.9361 0.9371 0.9374 0.8663 0.8606 0.8586 10StgDef-$\pi $ 0.9671 0.9661 0.9656 1.2385 1.2180 1.2104 10StgDef-$\eta $ 0.6813 0.6808 0.6804 3.8591 3.8824 3.8913 10StgDef-$ \dot m $ 0.9419 0.9406 0.9400 1.9707 1.9766 1.9784 表 11 不同${\boldsymbol{\sigma}} _{\bf{b}}^{\bf{2}}$值的CK-IBNN模型精度指标值
Table 11. Metric value of CK-IBNN model with various ${\boldsymbol{\sigma}} _{\bf{b}}^{\bf{2}}$
问题名称 R2 MAE $\sigma _{\text{b}}^2$=1.2 $\sigma _{\text{b}}^2$=1.6 $\sigma _{\text{b}}^2$=2.0 $\sigma _{\text{b}}^2$=1.2 $\sigma _{\text{b}}^2$=1.6 $\sigma _{\text{b}}^2$=2.0 10StgVGV-$\pi $ 0.9837 0.9841 0.9843 0.2127 0.2091 0.2058 10StgVGV-$\eta $ 0.8650 0.8669 0.8681 0.7550 0.7485 0.7437 10StgVGV-$ \dot m $ 0.9501 0.9512 0.9519 0.4003 0.3987 0.3978 Rotor 37-$\eta $ 0.9555 0.9554 0.9552 2.0429 2.0445 2.0459 Rotor 37-$\pi $ 0.9934 0.9935 0.9935 0.6321 0.6312 0.6303 3StgDef-$\eta $ 0.9589 0.9589 0.9589 0.7069 0.7095 0.7116 3StgDef-$\pi $ 0.8791 0.8789 0.8786 2.0972 2.0978 2.0984 3StgDef-$ \dot m $ 0.9373 0.9371 0.9368 0.8592 0.8606 0.8619 10StgDef-$\pi $ 0.9657 0.9661 0.9663 1.2128 1.2180 1.2227 10StgDef-$\eta $ 0.6805 0.6808 0.6809 3.8885 3.8824 3.8769 10StgDef-$ \dot m $ 0.9402 0.9406 0.9409 1.9778 1.9766 1.9753 -
[1] 黄磊, 张军, 何旭东, 等. 核心驱动风扇与压气机匹配特性研究[J]. 航空动力学报, 2024, 39(10): 20220798. HUANG Lei, ZHANG Jun, HE Xudong, et al. The matching performance research of core driven fan stage and compressor[J]. Journal of Aerospace Power, 2024, 39(10): 20220798. (in ChineseHUANG Lei, ZHANG Jun, HE Xudong, et al. The matching performance research of core driven fan stage and compressor[J]. Journal of Aerospace Power, 2024, 39(10): 20220798. (in Chinese) [2] 范忠岗, 巴顿, 邱佳慧, 等. 轴流压气机叶片与机匣处理一体化优化设计[J]. 航空动力学报, 2024, 39(7): 20220069. FAN Zhonggang, BA Dun, QIU Jiahui, et al. Integrated design optimization of blade and casing treatment in axial compressor[J]. Journal of Aerospace Power, 2024, 39(7): 20220069. (in ChineseFAN Zhonggang, BA Dun, QIU Jiahui, et al. Integrated design optimization of blade and casing treatment in axial compressor[J]. Journal of Aerospace Power, 2024, 39(7): 20220069. (in Chinese) [3] 李斌, 严红明, 李方刚, 等. 基于骨架特性的压气机可调叶片模型特性修正[J]. 航空动力学报, 2024, 39(7): 20220459. LI Bin, YAN Hongming, LI Fanggang, et al. Compressor variable stator vane model performance correction based on backbone map[J]. Journal of Aerospace Power, 2024, 39(7): 20220459. (in ChineseLI Bin, YAN Hongming, LI Fanggang, et al. Compressor variable stator vane model performance correction based on backbone map[J]. Journal of Aerospace Power, 2024, 39(7): 20220459. (in Chinese) [4] 韩忠华, 许晨舟, 乔建领, 等. 基于代理模型的高效全局气动优化设计综述方法研究进展[J]. 航空学报, 2020, 41(5): 623344. HAN Zhonghua, XU Chenzhou, QIAN Jianlin, et al. Recent progress of efficient global aerodynamic shape optimization using surrogate-based approach[J]. Acta Aeronautica et Astronautica Sinica, 2020, 41(5): 623344. (in ChineseHAN Zhonghua, XU Chenzhou, QIAN Jianlin, et al. Recent progress of efficient global aerodynamic shape optimization using surrogate-based approach[J]. Acta Aeronautica et Astronautica Sinica, 2020, 41(5): 623344. (in Chinese) [5] 吕利叶, 鲁玉军, 王硕, 等. 代理模型技术及其应用: 现状与展望[J]. 机械工程学报, 2024, 60(3): 254-281. LÜ Liye, LU Yujun, WANG Shuo, et al. Survey and prospect of surrogate model technique and application[J]. Journal of Mechanical Engineering, 2024, 60(3): 254-281. (in ChineseLÜ Liye, LU Yujun, WANG Shuo, et al. Survey and prospect of surrogate model technique and application[J]. Journal of Mechanical Engineering, 2024, 60(3): 254-281. (in Chinese) [6] 赵欢, 高正红, 夏露. 基于新型多可信度代理模型的多目标优化方法[J]. 航空学报, 2023, 44(6): 126962. ZHAO Huan, GAO Zhenghong, XIA Lu. Novel multi-fidelity surrogate model assisted many-objective optimization method[J]. Acta Aeronautica et Astronautica Sinica, 2023, 44(6): 126962. (in ChineseZHAO Huan, GAO Zhenghong, XIA Lu. Novel multi-fidelity surrogate model assisted many-objective optimization method[J]. Acta Aeronautica et Astronautica Sinica, 2023, 44(6): 126962. (in Chinese) [7] WIEGAND M, PROTS A, MEYER M, et al. Robust design optimization of a compressor rotor using recursive cokriging based multi-fidelity uncertainty quantification and multi-fidelity optimization[J]. Journal of Turbomachinery, 2025, 147(6): 061009. doi: 10.1115/1.4067076 [8] YANG Jinguang, ZHANG Min, PENG Cheng, et al. Stator re-stagger optimization in multistage axial compressor[J]. Propulsion and Power Research, 2021, 10(2): 107-117. doi: 10.1016/j.jppr.2021.03.002 [9] WANG Qineng, SONG Liming, GUO Zhendong, et al. A novel multi-fidelity surrogate for efficient turbine design optimization[J]. Journal of Turbomachinery, 2024, 146(4): 041011. doi: 10.1115/1.4064228 [10] DU Qiuwan, LI Yunzhu, YANG Like, et al. Performance prediction and design optimization of turbine blade profile with deep learning method[J]. Energy, 2022, 254: 124351. doi: 10.1016/j.energy.2022.124351 [11] HE Youwei, LUO Jinliang. Efficient hierarchical Kriging modeling method for high-dimension multi-fidelity problems[J]. Chinese Journal of Mechanical Engineering, 2024, 37(1): 151. doi: 10.1186/s10033-024-01136-z [12] MAGRIS M, IOSIFIDIS A. Bayesian learning for neural networks: an algorithmic survey[J]. Artificial Intelligence Review, 2023, 56(10): 11773-11823. doi: 10.1007/s10462-023-10443-1 [13] KENNEDY M. Predicting the output from a complex computer code when fast approximations are available[J]. Biometrika, 2000, 87(1): 1-13. doi: 10.1093/biomet/87.1.1 [14] 周奇, 杨扬, 宋学官, 等. 变可信度近似模型及其在复杂装备优化设计中的应用研究进展[J]. 机械工程学报, 2020, 56(24): 219-245. ZHOU Qi, YANG Yang, SONG Xueguan, et al. Survey of multi-fidelity surrogate models and their applications in the design and optimization of engineering equipment[J]. Journal of Mechanical Engineering, 2020, 56(24): 219-245. (in Chinese doi: 10.3901/JME.2020.24.219ZHOU Qi, YANG Yang, SONG Xueguan, et al. Survey of multi-fidelity surrogate models and their applications in the design and optimization of engineering equipment[J]. Journal of Mechanical Engineering, 2020, 56(24): 219-245. (in Chinese) doi: 10.3901/JME.2020.24.219 [15] PALAR P S, PARUSSINI L, BREGANT L, et al. On kernel functions for bi-fidelity Gaussian process regressions[J]. Structural and Multidisciplinary Optimization, 2023, 66(2): 37. doi: 10.1007/s00158-023-03487-y [16] LI Y, RUDNER T G J, WILSON A. A study of Bayesian neural network surrogates for Bayesian optimization[C]//Proceedings of International Conference on Learning Representations. San Diego, US: OpenReview, 2023: 1-39. [17] COUCKUYT I, DHAENE T, DEMEESTER P. ooDACE toolbox: a flexible object-oriented Kriging implementation[J]. Journal of Machine Learning Research, 2014, 15(1): 3183-3186. [18] WILLIAMS C K I. Computing with infinite networks[C]//Advances in Neural Information Processing Systems. Red Hook, US: Curran Associates Incorporation, 1997: 295-301. [19] LEE J, BAHRI Y, NOVAK R, et al. Deep neural networks as Gaussian processes[C]//Proceedings of the 6th International Conference on Learning Representations. San Diego, US: OpenReview, 2018: 1-17. [20] NOVAK R, XIAO L, HRON J, et al. Neural tangents: fast and easy infinite neural networks in Python[C]//Proceedings of the 8th International Conference on Learning Representations. San Diego, US: OpenReview, 2020: 1-19. [21] DANIELY A, FROSTIG R, SINGER Y. Toward deeper understanding of neural networks: The power of initialization and a dual view on expressivity[C]//Advances in Neural Information Processing Systems. Red Hook, US: Curran Associates Incorporation, 2016: 2261-2269. [22] DENTON J D. Multall: an open source, computational fluid dynamics based, turbomachinery design system[J]. Journal of Turbomachinery, 2017, 139(12): 121001. doi: 10.1115/1.4037819 [23] WANG Xuesong, SUN Jinju, BENINI E, et al. An overall blockage attenuation-based aerodynamic performance and stability design optimization method for transonic axial flow compressors[J]. International Journal of Numerical Methods for Heat and Fluid Flow, 2023, 33(5): 1853-1885. doi: 10.1108/HFF-07-2022-0437 [24] DASADHIKARI K, KUWAMURA Y. Rapid algorithmic blade design applying machine learning from shape optimization to satisfy multidisciplinary constraints[J]. Journal of Turbomachinery, 2024, 146(4): 041010. doi: 10.1115/1.4064227 [25] HE Youwei, GUI Qingwen, LUO Jinliang. An efficient parallel multi-fidelity multi-objective Bayesian optimization method and application to 3-stage axial compressor with 144 variables[J]. Aerospace Science and Technology, 2024, 150: 109235. doi: 10.1016/j.ast.2024.109235 [26] MENG Xuhui, KARNIADAKIS G E. A composite neural network that learns from multi-fidelity data: Application to function approximation and inverse PDE problems[J]. Journal of Computational Physics, 2020, 401: 109020. doi: 10.1016/j.jcp.2019.109020 [27] WU Xiaojing, ZUO Zijun, MA Long, et al. Multi-fidelity neural network-based aerodynamic optimization framework for propeller design in electric aircraft[J]. Aerospace Science and Technology, 2024, 146: 108963. doi: 10.1016/j.ast.2024.108963 [28] FU Chongbo, WANG Peng, ZHAO Liang, et al. A distance correlation-based Kriging modeling method for high-dimensional problems[J]. Knowledge-Based Systems, 2020, 206: 106356. doi: 10.1016/j.knosys.2020.106356 -

下载: