异物穿越螺旋桨概率分析方法研究及应用

朱嘉伟; 杨雪鹤; 陈垦伦; 雷柏茂

doi:10.13224/j.cnki.jasp.20230476

异物穿越螺旋桨概率分析方法研究及应用

doi: 10.13224/j.cnki.jasp.20230476

朱嘉伟^{1, 2, 3,},
杨雪鹤^{1, 2},
陈垦伦^{1, 3},
雷柏茂^{1, 2, 3,*, ,}

1.
工业和信息化部电子第五研究所，广州 511370
2.
工业和信息化部电子第五研究所广东省电子信息产品可靠性技术重点实验室，广州 511370
3.
工业和信息化部电子第五研究所民用飞机及航空发动机质量与可靠性工程技术工业和信息化部重点实验室，广州 511370

基金项目: 民用飞机专项科研技术研究项目（MJ-2020-F-10）

详细信息

作者简介:
朱嘉伟（1985-），男，高级工程师，硕士，研究领域为航空发动机排异特性、可靠性与适航性。E-mail：zhujiawei@ceprei.com

通讯作者:
雷柏茂（1987-），男，高级工程师，博士，研究领域为航空发动机结构强度与适航性。E-mail：leibaimao@126.com

中图分类号: V235.12
计量
- 文章访问数: 457
- HTML浏览量: 222
- PDF量: 34
- 被引次数: 0
出版历程
- 收稿日期: 2023-07-24
- 网络出版日期: 2024-08-22

Research and application of a novel simulation method for evaluating the traversal chance of flying object through the propeller

ZHU Jiawei^{1, 2, 3
,},
YANG Xuehe^{1, 2},
CHEN Kenlun^{1, 3},
LEI Baimao^{1, 2, 3,*
, ,}

1.
The Fifth Electronics Research Institute of Ministry of Industry and Information Technology，Guangzhou 511370，China
2.
Guangdong Provincial Key Laboratory of Electronic Information Products Reliability Technology，The Fifth Electronics Research Institute of Ministry of Industry and Information Technology，Guangzhou 511370，China
3.
The Ministry of Industry and Information Technology Key Laboratory of Quality and Reliability Engineering Technology of Civil Aircraft and Aero-Engine，The Fifth Electronics Research Institute of Ministry of Industry and Information Technology，Guangzhou 511370，China

摘要

摘要:
螺旋桨形状和外来异物飞行姿态复杂多样，导致异物穿越螺旋桨概率定量分析存在着困难，目前行业缺乏有效的技术手段。提出了一种可针对圆柱体异物，在螺旋桨真实形貌和异物飞行姿态可调等多种工程实际应用需求情况下开展异物穿越螺旋桨概率分析的技术方法，通过对比异物正撞击穿越长方体桨叶概率的理论计算值和仿真计算值的差异，验证了该方法的准确性，最后选取了某型螺旋桨开展异物穿越概率仿真分析。研究结果表明：该技术方法分析精度高，理想模型的仿真结果与理论穿越概率的最大绝对误差和最大相对误差为分别为0.92%和1.10%；异物穿越某型螺旋桨的概率所呈现的变化趋势与异物正撞击穿越长方体桨叶的概率相一致；在相同异物撞击条件下，穿越概率会随着异物长径比的增加而降低，随着姿态角的增加而增加。
- 外来异物 /
- 螺旋桨 /
- 碰撞检测 /
- 穿越概率 /
- 蒙特卡洛算法
Abstract:
Chance quantitative analysis of foreign object crossing through the propeller is difficult due to the complexity of propeller shape and foreign object flight attitude, and there is lack of effective technical means for it. A technical method to analyze the chance of cylindrical foreign object crossing the propeller under multiple engineering applications was presented, such as realistic three-dimensional propeller shape and adjustable foreign object flight attitude, etc. The accuracy of the technical method was verified by comparing the difference between the theoretical and simulated chance of foreign object crossing the rectangular blade propeller. Finally, a certain type propeller was selected to carry out the traversal chance simulation. The results showed that, the technical method presented herein had high analytical accuracy, the maximum absolute error and relative error between traversal chance results of simulation and theoretical calculation under ideal model were 0.92% and 1.10%. Respectively, the traversal chance variation trend of a certain type propeller was consistent with the rectangular blade propeller. Under the same impact conditions of foreign object, traversal chance decreased with the increase of aspect ratio and increased with the increase of impact attitude angle.
- foreign object /
- propeller /
- collision detection /
- traversal chance /
- Monte Carlo algorithm

HTML

图 1 异物穿越螺旋桨概率分析技术过程

Figure 1. Technical process of traversal chance analysis for foreign object crossing through the propeller

下载: 全尺寸图片幻灯片

图 2 圆柱体内、外接球示意图

Figure 2. Inner and outer ball model of cylinder

下载: 全尺寸图片幻灯片

图 3 螺旋桨空间表面数字离散化构建

Figure 3. Discretized construction of propeller’s surface

下载: 全尺寸图片幻灯片

图 4 异物的运动轨迹方程构建过程

Figure 4. Process of coordinate system transformation

下载: 全尺寸图片幻灯片

图 5 异物正几何曲面坐标旋转变换过程

Figure 5. Process of coordinate system transformation before edge collision detection

下载: 全尺寸图片幻灯片

图 6 异物穿越螺旋桨概率仿真计算分析软件界面

Figure 6. Operating interface of simulation software for traversal chance of flying object through the propeller

下载: 全尺寸图片幻灯片

图 7 长方体桨叶几何尺寸示意图

Figure 7. Geometric relationship between foreign object and blade when they are tangent

下载: 全尺寸图片幻灯片

图 8 异物穿越桨叶过程

Figure 8. Process of foreign object fly across the blade

下载: 全尺寸图片幻灯片

图 9 相邻桨叶的旋转过程

Figure 9. Process of blade rotation

下载: 全尺寸图片幻灯片

图 10 长方体桨叶螺旋桨数字离散化构建

Figure 10. Discretized construction of cuboid propeller

下载: 全尺寸图片幻灯片

图 11 螺旋桨三维结构图

Figure 11. 3D model of a certain type of propeller

下载: 全尺寸图片幻灯片

图 12 螺旋桨空间表面数字离散化构建图

Figure 12. Discretized construction for spatial surface of a certain type of propeller

下载: 全尺寸图片幻灯片

图 13 不同异物撞击速度的穿越概率对比

Figure 13. Comparison of traversal chance for different foreign object impact velocities

下载: 全尺寸图片幻灯片

图 14 不同长径比的穿越概率对比

Figure 14. Comparison of traversal chance for different aspect ratio of h/d

下载: 全尺寸图片幻灯片

图 15 异物撞击姿态角示意图

Figure 15. Illustration of foreign object’s impact attitude angle

下载: 全尺寸图片幻灯片

表 1 异物形状及其内、外接球构建方式

Table 1. Construction way of cylindrical foreign object and its inner and outer ball

参数	数值及详情
异物形状	圆柱体
异物几何特征参数	h，d
内接球半径r	min（h/2，d/2）
外接球半径R	$ {{\sqrt {{h^2} + {d^2}} } / 2} $

下载: 导出CSV

表 2 螺旋桨参数

Table 2. Parameters of propeller

参数	数值及说明
桨叶数量n/片	4，6，8
转速ω/（r/min）	920
桨叶形状	长方体
桨叶宽度a/mm	400
桨叶厚度b/mm	200

下载: 导出CSV

表 3 圆柱形异物参数

Table 3. Parameters of cylindrical foreign object

参数	数值
异物撞击速度v/（m/s）	80
高h/mm	165.7
直径d/mm	82.85

下载: 导出CSV

表 4 临界速度下的穿越概率仿真情况（n=4）

Table 4. Simulation results of traversal chance under critical velocity （n=4）

径向位置R_d/mm	临界速度v₁/（m/s）	仿真穿越概率P_s/%
1000	32.53	0.67
1500	28.24	0.83
2000	26.51	0.92
2500	25.58	0.72
3000	24.99	0.79
3500	24.59	0.84

下载: 导出CSV

表 5 临界速度下的穿越概率仿真情况（n=6）

Table 5. Simulation results of traversal chance under critical velocity （n=6）

径向位置R_d/mm	临界速度v₁/（m/s）	仿真穿越概率P_s/%
1000	62.97	0.59
1500	48.67	0.88
2000	43.76	0.00
2500	41.27	0.79
3000	39.76	0
3500	38.75	0

下载: 导出CSV

表 6 临界速度下的穿越概率仿真情况（n=8）

Table 6. Simulation results of traversal chance under critical velocity （n=8）

径向位置R_d/mm	临界速度v₁/（m/s）	仿真穿越概率P_s/%
1000	118.34	0.41
1500	76.25	0.56
2000	64.84	0.68
2500	59.52	0.65
3000	56.44	0.73
3500	54.43	0.75

下载: 导出CSV

表 7 不同径向撞击位置穿越概率对比情况（n=4）

Table 7. Simulation results of traversal chance under different impact radius （n=4）

径向位置R_d/mm	理论穿越概率P_x/%	仿真穿越概率P_s/%	相对误差E_r/%
1500	51.38	51.96	0.58
1750	54.34	54.07	0.27
2000	56.56	56.27	0.29
2250	58.27	57.92	0.35
2500	59.65	59.12	0.53
2750	60.77	60.43	0.34
3000	61.71	61.85	0.14
3250	62.50	63.39	0.89
3500	63.17	64.27	1.10
3750	63.38	63.76	0.38

下载: 导出CSV

表 8 不同径向撞击位置穿越概率对比情况（n=6）

Table 8. Simulation results of traversal chance under different impact radius （n=6）

径向位置R_d/mm	理论穿越概率P_x/%	仿真穿越概率P_s/%	相对误差E_r/%
1500	27.07	26.68	0.39
1750	31.51	31.65	0.14
2000	34.83	35.40	0.57
2250	37.41	37.13	0.28
2500	39.47	39.56	0.09
2750	41.16	41.69	0.53
3000	42.56	43.44	0.88
3250	43.74	43.57	0.17
3500	44.76	44.31	0.45
3750	45.64	46.27	0.63

下载: 导出CSV

表 9 不同径向撞击位置穿越概率对比情况（n=8）

Table 9. Simulation results of traversal chance under different impact radius （n=8）

径向位置R_d/mm	理论穿越概率P_x/%	仿真穿越概率P_s/%	相对误差E_r/%
1500	2.76	2.69	0.07
1750	8.68	8.92	0.24
2000	13.11	13.33	0.22
2250	16.55	17.03	0.48
2500	19.30	19.21	0.09
2750	21.54	21.59	0.05
3000	23.41	23.75	0.34
3250	24.99	25.29	0.30
3500	26.35	26.25	0.10
3750	27.52	28.15	0.63

下载: 导出CSV

表 10 穿越概率仿真条件

Table 10. Traversal chance simulation conditions for a certain type of propeller

参数	数值及说明
螺旋桨转速/（r/min）	920
异物撞击姿态	正撞击
异物撞击速度v/（m/s）	50，70，80，120
撞击径向位置R_d/mm	600～1700 （每100 mm一个撞击点，共12个撞击点）

下载: 导出CSV

表 11 不同径向撞击位置穿越概率仿真结果（v=50 m/s）

Table 11. Simulation results of traversal chance under different impact radius （v=50 m/s）

撞击径向位置R_d/mm	仿真穿越概率P_s/%	撞击径向位置R_d/mm	仿真穿越概率P_s/%
600	2.92	1200	24.14
700	12.48	1300	25.10
800	16.78	1400	23.80
900	20.70	1500	24.02
1000	23.56	1600	23.94
1100	24.48	1700	23.28

下载: 导出CSV

表 12 不同径向撞击位置穿越概率仿真结果（v=70 m/s）

Table 12. Simulation results of traversal chance under different impact radius （v=70 m/s）

撞击径向位置R_d/mm	仿真穿越概率P_s/%	撞击径向位置R_d/mm	仿真穿越概率P_s/%
600	11.88	1200	34.56
700	19.38	1300	34.64
800	25.90	1400	34.30
900	28.66	1500	35.26
1000	32.16	1600	35.60
1100	34.46	1700	36.32

下载: 导出CSV

表 13 不同径向撞击位置穿越概率仿真结果（v=80 m/s）

Table 13. Simulation results of traversal chance under different impact radius （v=80 m/s）

撞击径向位置R_d/mm	仿真穿越概率P_s/%	撞击径向位置R_d/mm	仿真穿越概率P_s/%
600	14.44	1200	36.78
700	23.22	1300	38.34
800	27.36	1400	39.68
900	30.76	1500	39.04
1000	33.92	1600	40.80
1100	35.58	1700	39.54

下载: 导出CSV

表 14 不同径向撞击位置穿越概率仿真结果（v=120 m/s）

Table 14. Simulation results of traversal chance under different impact radius （v=120 m/s）

撞击径向位置R_d/mm	仿真穿越概率P_s/%	撞击径向位置R_d/mm	仿真穿越概率P_s/%
600	19.20	1200	43.88
700	25.96	1300	45.08
800	33.30	1400	46.36
900	37.22	1500	47.08
1000	39.36	1600	47.60
1100	41.50	1700	48.62

下载: 导出CSV

表 15 不同撞击姿态角穿越概率仿真结果

Table 15. Simulation results of traversal chance under different impact attitude angles

撞击径向位置 R_d/mm	撞击速度 v/（m/s）	撞击姿态角/ （°）	仿真穿越概率P_s/%
1400	120	0	46.36
		45	49.16
		90	57.44

下载: 导出CSV

参考文献(22)

[1]	顾晨轩,苏艳,王辉. 国内外运输类飞机鸟撞适航条款及其修订背景研究分析[J]. 科技创新导报,2017,14(7): 247-249. GU Chenxuan,SU Yan,WANG Hui. Research and analysis on airworthiness clause of bird impact of transport aircraft at home and abroad and its revision background[J]. Science and Technology Innovation Herald,2017,14(7): 247-249. (in Chinese GU Chenxuan, SU Yan, WANG Hui. Research and analysis on airworthiness clause of bird impact of transport aircraft at home and abroad and its revision background[J]. Science and Technology Innovation Herald, 2017, 14(7): 247-249. (in Chinese)
[2]	罗刚,张海洋,吴春波,等. 航空发动机吞鸟与鸟撞飞机适航通用分析方法[J]. 航空发动机,2019,45(6): 90-96. LUO Gang,ZHANG Haiyang,WU Chunbo,et al. General airworthiness analysis method of aeroengine bird ingestion and aircraft bird impact[J]. Aeroengine,2019,45(6): 90-96. (in Chinese LUO Gang, ZHANG Haiyang, WU Chunbo, et al. General airworthiness analysis method of aeroengine bird ingestion and aircraft bird impact[J]. Aeroengine, 2019, 45(6): 90-96. (in Chinese)
[3]	何光宇,李应红,柴艳,等. 航空发动机压气机叶片砂尘冲蚀防护涂层关键问题综述[J]. 航空学报,2015,36(6): 1733-1743. HE Guangyu,LI Yinghong,CHAI Yan,et al. Review of key issues on coating against sand erosion of aero-engine compressor blade[J]. Acta Aeronautica et Astronautica Sinica,2015,36(6): 1733-1743. (in Chinese HE Guangyu, LI Yinghong, CHAI Yan, et al. Review of key issues on coating against sand erosion of aero-engine compressor blade[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(6): 1733-1743. (in Chinese)
[4]	汪洋,李玉龙. 冰雹冲击复合材料层合板仿真研究[J]. 振动与冲击,2015,34(2): 187-190,203. WANG Yang,LI Yulong. Modeling of high velocity hailstone impact onto composite material panel[J]. Journal of Vibration and Shock,2015,34(2): 187-190,203. (in Chinese WANG Yang, LI Yulong. Modeling of high velocity hailstone impact onto composite material panel[J]. Journal of Vibration and Shock, 2015, 34(2): 187-190, 203. (in Chinese)
[5]	乔艳江,宋显宏,胡锦翔,等. 道面外来物撞击航空发动机叶片现场鉴定方法探究[C]//全国失效分析学术会议论文集. 北京: 中国机械工程学会失效分析分会,2013: 49: 274-277. QIAO Yanjiang,SONG Xianhong,HU Jinxiang,et al. Research on field appraise method of aero-engine blade damaged by foreign object from runway[C]//Proceedings ofNational Failure Analysis Academic Conference. Beijing: Chinese Failure Analysis Institution of Chinese Mechanical Engineering Society (CMES),2013: 49: 274-277. (in Chinese QIAO Yanjiang, SONG Xianhong, HU Jinxiang, et al. Research on field appraise method of aero-engine blade damaged by foreign object from runway[C]//Proceedings ofNational Failure Analysis Academic Conference. Beijing: Chinese Failure Analysis Institution of Chinese Mechanical Engineering Society (CMES), 2013: 49: 274-277. (in Chinese)
[6]	张柱国,郝一鸣. 运输类飞机鸟撞及符合性验证综述[J]. 航空科学技术,2013,24(6): 1-4. ZHANG Zhuguo,HAO Yiming. Review about bird strike and substantiation for transport category airplane[J]. Aeronautical Science and Technology,2013,24(6): 1-4. (in Chinese ZHANG Zhuguo, HAO Yiming. Review about bird strike and substantiation for transport category airplane[J]. Aeronautical Science and Technology, 2013, 24(6): 1-4. (in Chinese)
[7]	金镭,刘友丹. 运输类飞机鸟撞适航性技术研究[J]. 航空标准化与质量,2012(6): 11-15,27. JIN Lei,LIU Youdan. Research on airworthiness technology of bird impact of transport aircraft[J]. Aeronautic Standardization and Quality,2012(6): 11-15,27. (in Chinese JIN Lei, LIU Youdan. Research on airworthiness technology of bird impact of transport aircraft[J]. Aeronautic Standardization and Quality, 2012(6): 11-15, 27. (in Chinese)
[8]	张清,崔博,霍枫. 涡扇发动机吞雹适航符合性验证试验方法[J]. 沈阳航空航天大学学报,2021,38(5): 23-29. ZHANG Qing,CUI Bo,HUO Feng. Test method of hail ingestion airworthiness conformance for turbofan engine[J]. Journal of Shenyang Aerospace University,2021,38(5): 23-29. (in Chinese ZHANG Qing, CUI Bo, HUO Feng. Test method of hail ingestion airworthiness conformance for turbofan engine[J]. Journal of Shenyang Aerospace University, 2021, 38(5): 23-29. (in Chinese)
[9]	魏雪霞,李广利. 沙粒对飞机发动机叶片的冲击碰撞[J]. 北京理工大学学报,2005,25(6): 471-474. WEI Xuexia,LI Guangli. The impact of sand on the blades of airplane engines[J]. Transactions of Beijing Institute of Technology,2005,25(6): 471-474. (in Chinese WEI Xuexia, LI Guangli. The impact of sand on the blades of airplane engines[J]. Transactions of Beijing Institute of Technology, 2005, 25(6): 471-474. (in Chinese)
[10]	李玉龙,石霄鹏. 民用飞机鸟撞研究现状[J]. 航空学报,2012,33(2): 189-198. LI Yulong,SHI Xiaopeng. Investigation of the present status of research on bird impacting on commercial airplanes[J]. Acta Aeronautica et Astronautica Sinica,2012,33(2): 189-198. (in Chinese LI Yulong, SHI Xiaopeng. Investigation of the present status of research on bird impacting on commercial airplanes[J]. Acta Aeronautica et Astronautica Sinica, 2012, 33(2): 189-198. (in Chinese)
[11]	李旦,赵廷渝,王永虎. 鸟撞飞机风挡数值模拟研究[J]. 四川理工学院学报(自然科学版),2017,30(4): 36-40. LI Dan,ZHAO Tingyu,WANG Yonghu. Numerical simulation of bird impact on aircraft windshield[J]. Journal of Sichuan University of Science & Engineering (Natural Science Edition),2017,30(4): 36-40. (in Chinese LI Dan, ZHAO Tingyu, WANG Yonghu. Numerical simulation of bird impact on aircraft windshield[J]. Journal of Sichuan University of Science & Engineering (Natural Science Edition), 2017, 30(4): 36-40. (in Chinese)
[12]	张宇,王彬文,白春玉,等. 典型航空器结构抗冰雹、仿真鸟弹冲击特性及其对比研究[J]. 振动与冲击,2021,40(18): 307-314. ZHANG Yu,WANG Binwen,BAI Chunyu,et al. Anti-impact characteristics and contrastive study of typical aircraft structures against hail and simulated bird projectile[J]. Journal of Vibration and Shock,2021,40(18): 307-314. (in Chinese ZHANG Yu, WANG Binwen, BAI Chunyu, et al. Anti-impact characteristics and contrastive study of typical aircraft structures against hail and simulated bird projectile[J]. Journal of Vibration and Shock, 2021, 40(18): 307-314. (in Chinese)
[13]	王利敏,张彦军,米百刚,等. 涡桨飞机发动机进气道排除异物特性数值研究[J]. 航空工程进展,2020,11(2): 264-271,278. WANG Limin,ZHANG Yanjun,MI Baigang,et al. Numerical simulation on excluding foreign objects from engine inlet of turboprop aircraft[J]. Advances in Aeronautical Science and Engineering,2020,11(2): 264-271,278. (in Chinese WANG Limin, ZHANG Yanjun, MI Baigang, et al. Numerical simulation on excluding foreign objects from engine inlet of turboprop aircraft[J]. Advances in Aeronautical Science and Engineering, 2020, 11(2): 264-271, 278. (in Chinese)
[14]	TUCKER V A. A mathematical model of bird collisions with wind turbine rotors[J]. Journal of Solar Energy Engineering,1996,118(4): 253-262. doi: 10.1115/1.2871788
[15]	DE LUCAS M,JANSS G,FERRER M. Birds and wind farms: risk assessment and mitigation[M]. 1st ed. Madrid,Spain: Quercus,2007.
[16]	MASDEN E A,COOK A S C P. Avian collision risk models for wind energy impact assessments[J]. Environmental Impact Assessment Review,2016,56: 43-49. doi: 10.1016/j.eiar.2015.09.001
[17]	PODOLSKY R. Method of and article of manufacture for determining probability of avian collision: US7315799[P]. 2008-01-01.
[18]	HOLMSTROM L A,HAMER T E,COLCLAZIER E M,et al. Assessing avian-wind turbine collision risk: an approach angle dependent model[J]. Wind Engineering,2011,35(3): 289-312. doi: 10.1260/0309-524X.35.3.289
[19]	BOLKER E D,HATCH J J,ZARA C. Modeling how windfarm geometry affects bird mortality[EB/OL]. (2014-08-06)[2023-10-28]. https://arxiv.org/abs/1408.1580.pdf
[20]	CALAVERA L P R,ARROYO S C,MONTANES L B. Spinner shed ice engine ingestion assessment methodology [C]//Proceedings of the 29th Congress of the International Council of the Aeronautical Sciences (ICAS). St. Petersburg,Russia: ICAS,2014: 1-10.
[21]	马力,姜甲玉,薛庆增. 航空发动机第1级风扇叶片鸟撞研究[J]. 航空发动机,2014,40(2): 65-69. MA Li,JIANG Jiayu,XUE Qingzeng. Research on bird impact of aeroengine first stage fan[J]. Aeroengine,2014,40(2): 65-69. (in Chinese MA Li, JIANG Jiayu, XUE Qingzeng. Research on bird impact of aeroengine first stage fan[J]. Aeroengine, 2014, 40(2): 65-69. (in Chinese)
[22]	蒋盼盼,韩涛锋,黄琪,等. 一种飞鸟越桨撞击风挡概率的评估方法: CN108090257B[P]. 2021-05-04. JIANG Panpan,HAN Taofeng,HUANG Qi,et al. An Evaluation Method for the Probability of a Bird Striking the Windshield by Crossing the Propeller: China Patent CN108090257B[P]. 2021-05-04. (in Chinese JIANG Panpan, HAN Taofeng, HUANG Qi, et al. An Evaluation Method for the Probability of a Bird Striking the Windshield by Crossing the Propeller: China Patent CN108090257B[P]. 2021-05-04. (in Chinese)