远场冲击波下螺旋桨毁伤与空化特征研究

王志凯; 郑景洲; 杨洋; 夏辉衡; 姚熊亮

doi:10.11883/bzycj-2023-0395

远场冲击波下螺旋桨毁伤与空化特征研究

doi: 10.11883/bzycj-2023-0395

1.
哈尔滨工程大学，船舶工程学院，黑龙江哈尔滨 150001
2.
中国船舶集团工程管理中心，北京 100081

基金项目: 国家重大专项（J2019I00170016）；国家自然科学基金（12472366）；黑龙江省优秀青年基金（YQ2021E009）

详细信息

作者简介:
王志凯（1989－　），男，博士，教授，wangzhikai@hrbeu.edu.cn

中图分类号: O383.1
计量
- 文章访问数: 190
- HTML全文浏览量: 16
- PDF下载量: 88
- 被引次数: 0
出版历程
- 收稿日期: 2023-10-31
- 修回日期: 2024-03-29
- 网络出版日期: 2024-04-01

Research on Damage and Cavitation Characteristics of Propellers under Far Field Shock Waves

1.
College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150001, China
2.
China Shipbuilding Group Engineering Management Center, Beijing 100081, China

摘要

摘要: 螺旋桨是舰船推进系统的核心部件，其运动稳定性和效率直接影响着舰船的性能。当前推进轴系抗冲击研究大多将螺旋桨等效成均质圆盘忽略其结构特征，不能准确得到水下爆炸瞬态冲击下螺旋桨的瞬态毁伤特征。故本文针对螺旋桨的结构特征，基于湿模态分析法得到实体建模优于壳体建模，开展了远场冲击波作用下螺旋桨物面空化冲击动响应及毁伤特征分析。并结合螺旋桨高速旋转状态下产生的水动力空化现象，进一步分析螺旋桨瞬态毁伤特征规律。得到了以下结论：在0°与90°攻角下，冲击波入射波作用于螺旋桨表面的物面载荷更高，但存在一个上限值，其与螺旋桨结构特征有关。在计及水动力空化状态下，桨叶的应力水平变化较为一致；桨叶主要塑性损伤区为叶根处，存在局部塑性和完全塑性两种模式。本文探讨了远场爆炸下螺旋桨毁伤与空化特征，研究结果可为推进轴系及螺旋桨抗冲击防护提供参考。
- 水下爆炸 /
- 螺旋桨 /
- 毁伤特性 /
- 空化效应
Abstract: The propeller is a critical component of a ship’s propulsion system that significantly influences the vessel’s performance through its stability and efficiency. Current research on the propulsion shaft system’s anti-shock properties often oversimplifies the propeller as a uniform circular disk, which disregards its structural intricacies and leads to inaccuracies in the transient damage characteristics during underwater explosions. This research focused on the propeller’s structural details and developed both an equivalent shell model and a more intricate solid model. Through structural wet modal numerical simulations, the study had determined that solid modeling outperforms shell modeling in accuracy. This finding is corroborated by coMParisons with empirical formulas, thereby validating the fluid-structure coupling analysis model.Building upon this foundation, the research examines the propeller’s transient Shock response and damage characteristics when subjected to far-field shockwaves. Utilizing the total wave algorithm in ABAQUS, the investigation extends to the cavitation and damage patterns of the propeller under such conditions, with confirmation provided by the one-dimensional Bleich-Sandler finite element model. To delve deeper into the phenomenon of hydrodynamic cavitation caused by the propeller’s high-speed rotation, the coupled Eulerian-Lagrangian (CEL) method was applied. Initially, a simplified propeller model was created to confirm the cavitation bubble layer’s fragmentation due to the flow field load resulting from explosive product expansion. Subsequent modifications to the propeller’s transient fluid-structure coupling calculation model allow for a more thorough analysis of its transient damage characteristics.The findings indicate that at attack angles of 0 and 90 degrees, the propeller surface experiences heightened shockwave loads, albeit with a threshold linked to the propeller’s structural properties. When hydrodynamic cavitation is factored in, the stress distribution on the propeller blade tends to be more uniform; the blade’s primary plastic damage is localized at the root, exhibiting both localized and complete plastic deformation patterns. This research elucidates the damage and cavitation effects on propellers due to far-field explosions, offering valuable insights for enhancing the anti-shock defenses of both the propulsion shaft system and the propeller itself.
- underwater explosion /
- propeller /
- damage characteristics /
- cavitation effect

HTML全文

图 1 水下兵器对舰船毁伤威力

Figure 1. The damage power of underwater weapons to ships

下载: 全尺寸图片幻灯片

图 2 螺旋桨水动力空泡

Figure 2. Propeller hydrodynamic cavitation

下载: 全尺寸图片幻灯片

图 3 螺旋桨模型

Figure 3. Propeller model

下载: 全尺寸图片幻灯片

图 4 螺旋桨一阶湿模态

Figure 4. The first-order wet mode of propeller

下载: 全尺寸图片幻灯片

图 5 Bleich—Sandler平板模型

Figure 5. Bleich-Sandler plate model

下载: 全尺寸图片幻灯片

图 6 冲击波压力曲线

Figure 6. Shock wave pressure curve

下载: 全尺寸图片幻灯片

图 7 一维Bleich—Sandler有限元模型

Figure 7. One-Dimensional Bleich-Sandler Finite Element Model

下载: 全尺寸图片幻灯片

图 8 平板速度随时间的变化曲线的理论值与计算值对比

Figure 8. Comparison of the plate velocity-time curve between theoretical and calculation values

下载: 全尺寸图片幻灯片

图 9 爆源位置示意图

Figure 9. Detonation source location diagram

下载: 全尺寸图片幻灯片

图 10 不同攻角物面空化区域

Figure 10. Cavitation area at different angles of attack

下载: 全尺寸图片幻灯片

图 11 空化区域大小随攻角变化

Figure 11. The size of cavitation area changes with the angle of attack.

下载: 全尺寸图片幻灯片

图 12 不同冲击因子物面空化区域

Figure 12. Cavitation zones at the material surface under different shock factors

下载: 全尺寸图片幻灯片

图 13 空化区域大小随冲击因子变化曲线

Figure 13. Curve of cavitation zone size variation with shock factor

下载: 全尺寸图片幻灯片

图 14 物面空化溃灭二次加载

Figure 14. Cavitation collapse secondary loading on object surface

下载: 全尺寸图片幻灯片

图 15 不同攻角二次加载峰值

Figure 15. Secondary loading peak values at different angles of attack

下载: 全尺寸图片幻灯片

图 16 90°攻角下二次加载峰值随网格尺寸变化

Figure 16. At 90 Degrees angle of attack, the variation of the secondary loading peak with mesh size

下载: 全尺寸图片幻灯片

图 18 冲击因子为1.0时不同攻角下不同测点位移响应

Figure 18. Deflection response of different measuring points at different angles of attack when the shock factor is 1.0

下载: 全尺寸图片幻灯片

图 17 测点位置

Figure 17. Position of measuring point

下载: 全尺寸图片幻灯片

图 19 冲击因子为1.3时不同攻角下不同测点塑性变形

Figure 19. Plastic deformation of different measuring points at different angles of attack when the shock factor is 1.3,

下载: 全尺寸图片幻灯片

图 20 爆源-空化层-物面冲击波传递数值模拟示意图

Figure 20. Simulation of explosion source-cavitation layer-surface shock wave transmission

下载: 全尺寸图片幻灯片

图 21 流场压力图

Figure 21. Flow field pressure contour

下载: 全尺寸图片幻灯片

图 22 空化泡层体积变化

Figure 22. Variation of cavitation bubble layer volume

下载: 全尺寸图片幻灯片

图 23 板面应力波传递云图

Figure 23. Stress wave propagation contour on the plate surface

下载: 全尺寸图片幻灯片

图 24 空化泡层示意图

Figure 24. Cavitation bubble layer diagram

下载: 全尺寸图片幻灯片

图 25 冲击载荷作用下桨叶轴向位移

Figure 25. Axial displacement of blade under shock load

下载: 全尺寸图片幻灯片

图 26 相同攻角不同冲击因子下典型测点位移响应

Figure 26. Deflection response of typical measuring points under different shock factors at the same angle of attack

下载: 全尺寸图片幻灯片

图 27 冲击因子为1.3时不同攻角下不同测点塑性变形

Figure 27. Part of the plastic deformation cloud diagram when the shock factor is 1.3,

下载: 全尺寸图片幻灯片

表 1 螺旋桨材料参数

Table 1. Propeller material parameters

密度/ (kg∙m⁻³)	杨氏模量/ MPa	泊松比	强度极限/ MPa
7600	1.17×10⁵	0.3	630

下载: 导出CSV

表 2 螺旋桨部分设计参数

Table 2. Design parameters of propeller part

直径D/ mm	叶数Z	切面形状	螺距H/ mm	螺距比 H/D	后倾角α	毂径比 d/D
6000	5	弓型	4200	0.7	15°	0.183

下载: 导出CSV

表 3 螺旋桨有限值与计算值误差

Table 3. Error table of propeller experimental value and calculated value

	实体一阶湿模态固有频率/ Hz		壳体一阶湿模态固有频率/Hz
	一阶挥舞	一阶扭转	一阶挥舞	一阶扭转
数值实验	14.87	35.81	12.01	31.91
经验公式	14.09	33.38	14.09	30.54
误差/%	+5.6%	+7.3%	−14.7%	−8.5%

下载: 导出CSV

表 4 计算工况

Table 4. Working condition

工况	攻角/(°)	冲击因子/（$ {\mathrm{k}\mathrm{g}}^{\frac{1}{2}}\cdot \mathrm{m}^{-1} $）	工况	攻角/(°)	冲击因子/（$ {\mathrm{k}\mathrm{g}}^{\frac{1}{2}}\cdot \mathrm{m}^{-1} $）
1	0	0.4	9	60	0.4
2		0.7	10		0.7
3		1.0	11		1.0
4		1.3	12		1.3
5	30	0.4	13	90	0.4
6		0.7	14		0.7
7		1.0	15		1.0
8		1.3	16		1.3

下载: 导出CSV

表 5 不同尺寸网格数量表

Table 5. Table of mesh quantities for different sizes

轴体尺寸/ mm	桨叶尺寸/ mm	网格数量
200	100	23 343
160	80	57 599
100	50	153 536
90	45	207 462
80	40	302 294
70	35	576 104
60	30	654 490
50	25	1 105 909

下载: 导出CSV

表 6 不同攻角塑性损伤情况

Table 6. Plastic damage at different angles of attack

攻角/(°)	根部塑性损伤叶片数量	3/4半径叶缘处塑性损伤峰值范围	3/4半径叶缘处塑性损伤平均值范围	塑性损伤范围总值
0	5	3.4%	3.0%	13.1%
30	4	1.9%	1.3%	7.5%
60	4	1.0%	0.5%	3.5%
90	5	3.2%	2.7%	14.4%

下载: 导出CSV

表 7 不同攻角塑性损伤情况

Table 7. Plastic damage at different angles of attack

攻角/（°）	叶根塑性损伤模式	塑性损伤范围
0	叶根局部塑性	4.4%
30	叶根局部塑性	4.2%
60	叶根整体塑性	5.6%
90	叶根整体塑性	7.9%

下载: 导出CSV

参考文献(17)

[1]	Materials Processing Technology, 2008, 206(1/2/3): 275-281. DOI: 10.1016/j.jmatprotec.2007.12.022.
[2]	许斐, 周力, 宗智. 铝合金舰艇在水下冲击波作用下动态响应的数值研究 [J]. 舰船科学技术, 2011, 33(1): 31–40. DOI: 10.3404/j.issn.1672-7649.2011.01.006. XU F, ZHOU L, ZONG Z. A numerical analysis on dynamic responses of aluminum vessels subjected to underwater shock wave [J]. Ship Science and Technology, 2011, 33(1): 31–40. DOI: 10.3404/j.issn.1672-7649.2011.006. DOI: 10.3404/j.issn.1672-7649.2011.01.006.
[3]	宗智, 陈岗, 叶帆, 等. 水下爆炸冲击波作用下的空化效应数值仿真研究 [J]. 船舶力学, 2014, 18(4): 385–394. DOI: 10.3969/j.issn.1007-7294.2014.04.005. ZONG Z, CHEN G, YE F, et al. Numerical simulation of cavitation caused by underwater explosion [J]. Journal of Ship Mechanics, 2014, 18(4): 385–394. DOI: 10.3969/j.issn.1007-7294.2014.04.005.
[4]	PAN J Y, FARAG N, LIN T, et al. Propeller induced structural vibration through the thrust bearing [C]// Innovation in Acoustics and Vibration Annual Conference of the Australian Acoustical Society. Adelaide, 2002.
[5]	KOBAYASHI S. Effects of shaft vibration on occurrence of asymmetric cavitation in inducer [J]. JSME International Journal Series B Fluids and Thermal Engineering, 2006, 49(4): 1220–1225. DOI: 10.1299/jsmeb.49.1220.
[6]	孙文林. 冰区航行船舶螺旋桨强度研究 [D]. 哈尔滨: 哈尔滨工程大学, 2016. DOI: 10.7666/d.D01105449. SUN W L. Research on the propeller strength of ice-going ships [D]. Harbin: Harbin Engineering University, 2016. DOI: 10.7666/d.D01105449.
[7]	徐学伟, 刘寅桐, 张耀林. 爆炸声源对舰船螺旋桨效能损伤评估研究 [J]. 舰船电子工程, 2018, 38(10): 153–155,164. DOI: 10.3969/j.issn.1672-9730.2018.10.038. XU X W, LIU Y T, ZHANG Y L. Research on evaluation of damage to ship propeller effectiveness for explosion sound source [J]. Ship Electronic Engineering, 2018, 38(10): 153–155,164. DOI: 10.3969/j.issn.1672-9730.2018.10.038.
[8]	UCAR H. Dynamic response of a catamaran-hull ship subjected to underwater explosions [D]. Monterey: Naval Postgraduate School, 2006.
[9]	孙帅, 常欣, 叶礼裕, 等. 四桨推进船舶螺旋桨负荷数值计算分析 [J]. 哈尔滨工程大学学报, 2017, 38(9): 1351–1358. DOI: 10.11990/jheu.201607004. SUN S, CHANG X, YE L Y, et al. Numerical analysis of the propeller load of a four-propeller propulsion ship [J]. Journal of Harbin Engineering University, 2017, 38(9): 1351–1358. DOI: 10.11990/jheu.201607004.
[10]	刘强, 王永生, 董蕾, 等. 船用螺旋桨振动特性数值计算与分析 [J]. 海军工程大学学报, 2016, 28(3): 17–23. DOI: 10.7495/j.issn.1009-3486. 2016.03.005. DOI: 10.7495/j.issn.1009-3486.2016.03.005. LIU Q, WANG Y S, DONG L, et al. Numerical analysis on vibration of marine propeller [J]. Journal of Naval University of Engineering, 2016, 28(3): 17–23. DOI: 10.7495/j.issn.1009-3486.2016.03.005.
[11]	娄本强, 嵇春艳. 船用螺旋桨流固耦合振动特性分析 [J]. 大连理工大学学报, 2019, 59(2): 154–161. DOI: 10.7511/dllgxb2019.02.007. LOU B Q, JI C Y. FSI vibrations characteristics analysis of marine propeller [J]. Journal of Dalian University of Technology, 2019, 59(2): 154–161. DOI: 10.7511/dllgxb2019.02.007.
[12]	李程辉. 基于声固耦合法的水下爆炸冲击波在复杂边界附近空化特性研究[D]. 哈尔滨: 哈尔滨工程大学, 2018: 26–30. LI C H. Investigation on cavitation of underwater explosion shockwave near complex boundaries with coupled CEL method [D]. Harbin: Harbin Engineering University, 2018: 26–30.
[13]	陈高杰, 贾则, 宗智, 等. 考虑空化效应的舰艇水下爆炸仿真研究 [J]. 爆破, 2012, 29(2): 1–3,25. DOI: 10.3963/j.issn.1001-487X.2012.29.0213. CHEN G J, JIA Z, ZONG Z, et al. Investigation on simulation of underwater explosion cavitation surface ship [J]. Blasting, 2012, 29(2): 1–3,25. DOI: 10.3963/j.issn.1001-487X.2012.29.0213.
[14]	BLEICH H H, SANDLER I S. Interaction between structures and bilinear fluids [J]. International Journal of Solids and Structures, 1970, 6(5): 617–639. DOI: 10.1016/0020-7683(70)90034-X.
[15]	马艳, 陆芳, 黄红波, 等. 集装箱船模型螺旋桨空泡与脉动压力比对试验研究[C]//第二十五届全国水动力学研讨会暨第十二届全国水动力学学术会议文集(上册). 北京: 海洋出版社, 2013.
[16]	SALVATORE F. Description of the INSEAN E779A propeller experimental dataset[R]. Technical Report INSEAN 2006-085 INSEAN-Italian Ship Model Basin. 2006.
[17]	寇晓枫, 王高辉, 卢文波, 等. 空气隔层对水下爆炸冲击波的缓冲效应 [J]. 振动与冲击, 2017, 36(3): 7–13. DOI: 10.13465/j.cnki.jvs.2017.03.002. KOU X F, WANG G H, LU W B, et al. Mitigation effects of air interlayer on underwater explosion shock wave [J]. Journal of Vibration and Shock, 2017, 36(3): 7–13. DOI: 10.13465/j.cnki.jvs.2017.03.002.