基于多目标优化的跨介质航行体水面滑跳初始运动参数设计

doi:10.12382/bgxb.2022.0704

摘要/Abstract

摘要：

基于耦合欧拉-拉格朗日方法,开展了初始运动参数对带滑板跨介质航行体近水面滑跳过程的影响规律研究。采用参数回归方法分别建立最大浸深、最大过载和俯仰角变化3个变量随航行体初始速度、入水角、滑板倾角变化的关系式,得到了跨介质航行体稳定滑跳的参数包络,在入水角2°~5°、滑板倾角10°~25°、速度60~80m/s范围之间,且在范围内各参数之间会相互制约。在此基础上,基于多目标优化算法获得了最适宜航行体稳定滑跳的初始参数为入水速度60.1m/s、滑板倾角18.8°、入水角2°,在此参数范围内,跨介质航行体的水面滑跳过程可以在不超过航行体过载要求的条件下获得相对最小的入水浸深和俯仰角变化。

关键词: 跨介质, 水面滑跳, 多目标优化, 耦合欧拉-拉格朗日方法

Abstract:

The influence of the initial kinematic parameters on the skipping process of a trans-media vehicle over water surface is numerically studied based on coupled Eulerian-Lagrangian method. The relationships among three parameters, including the maximum immersion depth, the maximum overload and the pitching angle variation, and the initial motion parameters, including the initial velocity and incident angle of trans-media vehicle and the inclination angle of skipping plate, are analyzed by means of the parameter regression method. The parameter envelope of stable skipping process of the trans-media vehicle is obtained at the incident angles ranging from 2 ° to 5 °, the inclination angles of skipping plate ranging from 10 ° to 25 °, and the initial velocities ranging from 60 to 80m/s. Based on the multi-objective optimization algorithm, the most suitable target parameters for the stable skipping process of the vehicle are obtained as follows: the initial velocity of 60.1m/s, the skipping plate’s inclination angle of 18.8 °, and the incident angle of 2 °. Within this parameter range, the relative minimum values of immersion depth and pitching angle variation can be obtained during the skipping process of the trans-media vehicle without exceeding the overload requirements.

Key words: trans-medium, skipping motion, multi-objective optimization, coupled Euler-Lagrangian method

中图分类号:

O368

孙士明, 郁伟, 王晓辉, 李振旺, 刘彩连. 基于多目标优化的跨介质航行体水面滑跳初始运动参数设计[J]. 兵工学报, 2024, 45(2): 541-551.

SUN Shiming, YU Wei, WANG Xiaohui, LI Zhenwang, LIU Cailian. Initial Kinematic Parameters Design of Trans-media Vehicle Skipping over Water Surface Based on Multi-objective Optimization[J]. Acta Armamentarii, 2024, 45(2): 541-551.

图/表 16

参考文献 23

[1]	JOHNSON W, REID S R. Ricochet of spheres off water[J]. Journal of Mechanical Engineering Science, 1975, 17(2): 71-81. doi: 10.1243/JMES_JOUR_1975_017_013_02 URL
[2]	HUTCHINGS I M. The ricochet of spheres and cylinders from the surface of water[J]. International Journal of Mechanical Sciences, 1976, 18(5): 243-247. doi: 10.1016/0020-7403(76)90006-0 URL
[3]	MILOH T, SHUKRON Y. Ricochet off water of spherical projectiles[J]. Journal of Ship Research, 1991, 35(2): 91-100. doi: 10.5957/jsr.1991.35.2.91 URL
[4]	BELDEN J, HURD R, JANDRON M. et al. Elastic spheres can walk on water[J]. Nature Communications, 2016, 7: 10551. doi: 10.1038/ncomms10551 pmid: 26842860
[5]	BOCQUET L. The physics of stone skipping[J]. American Journal of Physics, 2003, 71(2): 150-155. doi: 10.1119/1.1519232 URL
[6]	CLANET C, HERSEN F, BOCQUET L. Secrets of successful stone-skipping[J]. Nature, 2004, 427(6969): 29. doi: 10.1038/427029a
[7]	NAGAHIRO S, HAYAKAWA Y. Theoretical and numerical approach to “magic angle” of stone skipping[J]. Physical review letters, 2005, 94(17): 174501. doi: 10.1103/PhysRevLett.94.174501 URL
[8]	SUN S M, YAN K, CHEN W Q. Numerical study on kinetic characteristics of disk skipping over water surface[J]. Journal of Ship Mechanics, 2021, 25(6): 716-725.
[9]	ROSELLINI L, HERSEN F, CLANET C, et al. Skipping stones[J]. Journal of Fluid Mechanics, 2005, 543: 137-146. doi: 10.1017/S0022112005006373 URL
[10]	HEWITT I J, BALMFORTH N J, MCELWAINE J N. Continual skipping on water[J]. Journal of Fluid Mechanics, 2011, 669: 328-353. doi: 10.1017/S0022112010005057 URL
[11]	MURALI V, NAIK S D. Skipping stone to projectile ricochet[J]. International Review on Modelling and Simulations, 2015, 8(1): 104-110.
[12]	YAN R, MONAGHAN J J. SPH simulation of skipping stones[J]. European Journal of Mechanics B/Fluids, 2017, 61: 61-71. doi: 10.1016/j.euromechflu.2016.10.001 URL
[13]	PARK M, JUNG Y, PARK W. Numerical study of the impact force and ricochet behaviour of high speed water entry bodies[J]. Comput Fluids, 2003, 32: 939-951. doi: 10.1016/S0045-7930(02)00087-7 URL
[14]	杨超, 姜宇, 吴志刚. 基于边界元法的三维结构体滑跳运动数值仿真[J]. 北京航空航天大学学报, 2022, 48(9):1678-1691.
	YANG C, JIANG Y, WU Z G. Numerical simulation of skipping motion of three-dimensional structure based on boundary element method[J]. Journal of Beijing University of Aeronautics and Astronautics, 2022, 48(9):1678-1691. (in Chinese)
[15]	裴譞, 张宇文, 李闻白, 等. 跨介质飞行器气/水两相弹道仿真研究[J]. 工程力学, 2010, 27(8): 223-228.
	PEI X, ZHANG Y W, LI W B, et al. Simulation and analysis on the gas/water two-phase ballistics of trans-media aircraft[J]. Engineering Mechanics, 2010, 27(8): 223-228. (in Chinese)
[16]	王伟, 张宇文, 朱灼. 跨介质飞行器弹道仿真分析[J]. 计算机仿真, 2011, 28(12): 1-4.
	WANG W, ZHANG Y W, ZHU Z. Simulation and analysis on ballistic trajectory of trans-media aircraft[J]. Computer Simulation, 2011, 28(12): 1-4. (in Chinese)
[17]	侯东伯. 运动体触水滑跳过程运动特性研究[D]. 哈尔滨: 哈尔滨工业大学. 2020.
	HOU D B. Research on kinematic characteristics of moving object during water-skipping[D]. Harbin:Harbin Institute of Technology, 2020. (in Chinese)
[18]	田北晨, 李达钦, 刘涛涛, 等. 跨介质飞行器触水滑跳运动特性数值模拟[J]. 兵工学报, 2022, 43(3): 586-598. doi: 10.12382/bgxb.2021.0164
	TIAN B C, LI D Q, LIU T T, et al. Numerical simulation on kinematic characteristics of trans-media aircraft during water-skipping[J]. Acta Armamentarii, 2022, 43(3): 586-598. (in Chinese) doi: 10.12382/bgxb.2021.0164
[19]	周可, 黄振贵, 陈志华, 等. 跨介质航行器高速斜入水跳弹现象研究[J]. 装备环境工程, 2022, 19(5):39-48.
	ZHOU K., HUANG Z G, CHEN Z H, et al. Ricochet phenomenon of tans-media vehicle at high-speed oblique water-entry[J]. Equipment Environmental Engineering, 2022, 19(5): 39-48. (in Chinese)
[20]	刘祥, 李天雄, 李林, 等. 入水参数对“水漂式”航行体跨介质运动影响数值仿真分析[J]. 宇航总体技术, 2020, 4(3): 62-70.
	LIU X, LI T X, LI L, et al. Anlysis of influence of water entry parameters on cross-medium movement of “water drift” aerial-auqatic vehicle[J]. Astronautical System Engineering Technology, 2020, 4(3): 62-70. (in Chinese)
[21]	王晓辉, 褚学森, 冯光. 基于Abaqus显式CEL方法的球体入水数值研究[J]. 船舶力学, 2018, 22(7): 838-844.
	WANG X H, CHU X S, FENG G. Numerical simulation of sphere’s water-entry base on coupled Euler-Lagrangian method[J]. Ship Mechanics, 2018, 22(7): 838-844. (in Chinese)
[22]	周连第. 多元回归分析方法在船舶科研设计中的应用[M]. 北京: 国防工业出版社, 1979.
	ZHOU L D. Application of multiple regression analysis method in ship scientific research and design[M]. Beijing: National Defense Industry Press, 1979. (in Chinese)
[23]	张火明, 陆萍蓝, 张晓菲. 基于INSGA-Ⅱ的等效水深截断系泊系统优化设计[J]. 船舶力学, 2014, 18(11): 1292-1298.
	ZHANG M H, LU P L, ZHANG X F. Investigation on optimization design of equivalent water depth truncated mooring system based on INSGA-Ⅱ[J]. Ship Mechanics, 2014, 18(11): 1292-1298. (in Chinese)

编号	v/(m·s^-1)	β/(°)	α/(°)	d_max/m	N_max/g	Δθ/(°)	编号	v/(m·s^-1)	β/(°)	α/(°)	d_max/m	N_max/g	Δθ/(°)
1	70	5	2	0.036	13.9	-5.6	31	80	15	6	0.124	46.8	-18.5
2	70	5	4	0.098	17.7	-10.6	32	80	15	8	0.175	55.1	-20.6
3	70	5	6	0.185	23.0	-13.5	33	80	20	2	0.046	25.4	0.8
4	70	5	8	0.340	28.1	-13.2	34	80	20	4	0.076	42.4	-4.3
5	70	10	2	0.046	14.8	-4.5	35	80	20	6	0.112	52.9	-7.7
6	70	10	4	0.086	22.0	-12.3	36	80	20	8	0.157	63.6	-19.8
7	70	10	6	0.141	28.4	-15.7	37	80	25	2	0.047	23.8	1.3
8	70	10	8	0.238	33.6	-19.9	38	80	25	4	0.080	25.1	-2.5
9	70	15	2	0.045	19.6	-0.5	39	80	25	6	0.114	40.9	-7.8
10	70	15	4	0.081	27.2	-9.2	40	80	25	8	0.149	55.5	-15.1
11	70	15	6	0.124	36.7	-15.6	41	90	5	2	0.035	23.8	-7.8
12	70	15	8	0.179	45.8	-18.6	42	90	5	4	0.097	29.8	-11.2
13	70	20	2	0.046	19.5	0.7	43	90	5	6	0.248	36.6	-10.2
14	70	20	4	0.077	32.1	-4.2	44	90	5	8	0.310	40.2	-24.1
15	70	20	6	0.112	40.4	-5.9	45	90	10	2	0.045	23.5	-6.1
16	70	20	8	0.154	48.9	-19.8	46	90	10	4	0.087	35.6	-14.9
17	70	25	2	0.048	21.9	1.1	47	90	10	6	0.160	43.5	-19.9
18	70	25	4	0.081	30.9	-0.9	48	90	10	8	0.248	61.8	-21.7
19	70	25	6	0.115	42.2	-6.8	49	90	15	2	0.044	27.8	0.1
20	70	25	8	0.148	52.4	-14.7	50	90	15	4	0.080	44.9	-11.6
21	80	5	2	0.035	18.6	-6.5	51	90	15	6	0.128	60.1	-21.2
22	80	5	4	0.100	23.3	-11.5	52	90	15	8	0.180	70.3	-24.8
23	80	5	6	0.208	28.6	-14.3	53	90	20	2	0.045	32.0	1.0
24	80	5	8	0.386	37.2	-14.5	54	90	20	4	0.076	51.6	-4.7
25	80	10	2	0.045	18.0	-5.2	55	90	20	6	0.110	65.4	-10.6
26	80	10	4	0.087	27.4	-13.5	56	90	20	8	0.157	75.2	-19.7
27	80	10	6	0.160	32.0	-17.9	57	90	25	2	0.048	34.8	1.3
28	80	10	8	0.254	36.1	-20.8	58	90	25	4	0.080	54.0	-3.1
29	80	15	2	0.045	20.9	-0.5	59	90	25	6	0.114	67.9	-5.3
30	80	15	4	0.079	34.8	-10.4	60	90	25	8	0.148	83.2	-20.5

编号	v/(m·s^-1)	β/(°)	α/(°)	d_max/m	N_max/g	Δθ/(°)	编号	v/(m·s^-1)	β/(°)	α/(°)	d_max/m	N_max/g	Δθ/(°)
1	70	5	2	0.036	13.9	-5.6	31	80	15	6	0.124	46.8	-18.5
2	70	5	4	0.098	17.7	-10.6	32	80	15	8	0.175	55.1	-20.6
3	70	5	6	0.185	23.0	-13.5	33	80	20	2	0.046	25.4	0.8
4	70	5	8	0.340	28.1	-13.2	34	80	20	4	0.076	42.4	-4.3
5	70	10	2	0.046	14.8	-4.5	35	80	20	6	0.112	52.9	-7.7
6	70	10	4	0.086	22.0	-12.3	36	80	20	8	0.157	63.6	-19.8
7	70	10	6	0.141	28.4	-15.7	37	80	25	2	0.047	23.8	1.3
8	70	10	8	0.238	33.6	-19.9	38	80	25	4	0.080	25.1	-2.5
9	70	15	2	0.045	19.6	-0.5	39	80	25	6	0.114	40.9	-7.8
10	70	15	4	0.081	27.2	-9.2	40	80	25	8	0.149	55.5	-15.1
11	70	15	6	0.124	36.7	-15.6	41	90	5	2	0.035	23.8	-7.8
12	70	15	8	0.179	45.8	-18.6	42	90	5	4	0.097	29.8	-11.2
13	70	20	2	0.046	19.5	0.7	43	90	5	6	0.248	36.6	-10.2
14	70	20	4	0.077	32.1	-4.2	44	90	5	8	0.310	40.2	-24.1
15	70	20	6	0.112	40.4	-5.9	45	90	10	2	0.045	23.5	-6.1
16	70	20	8	0.154	48.9	-19.8	46	90	10	4	0.087	35.6	-14.9
17	70	25	2	0.048	21.9	1.1	47	90	10	6	0.160	43.5	-19.9
18	70	25	4	0.081	30.9	-0.9	48	90	10	8	0.248	61.8	-21.7
19	70	25	6	0.115	42.2	-6.8	49	90	15	2	0.044	27.8	0.1
20	70	25	8	0.148	52.4	-14.7	50	90	15	4	0.080	44.9	-11.6
21	80	5	2	0.035	18.6	-6.5	51	90	15	6	0.128	60.1	-21.2
22	80	5	4	0.100	23.3	-11.5	52	90	15	8	0.180	70.3	-24.8
23	80	5	6	0.208	28.6	-14.3	53	90	20	2	0.045	32.0	1.0
24	80	5	8	0.386	37.2	-14.5	54	90	20	4	0.076	51.6	-4.7
25	80	10	2	0.045	18.0	-5.2	55	90	20	6	0.110	65.4	-10.6
26	80	10	4	0.087	27.4	-13.5	56	90	20	8	0.157	75.2	-19.7
27	80	10	6	0.160	32.0	-17.9	57	90	25	2	0.048	34.8	1.3
28	80	10	8	0.254	36.1	-20.8	58	90	25	4	0.080	54.0	-3.1
29	80	15	2	0.045	20.9	-0.5	59	90	25	6	0.114	67.9	-5.3
30	80	15	4	0.079	34.8	-10.4	60	90	25	8	0.148	83.2	-20.5

初始变量			目标函数
v/ (m·s^-1)	β/(°)	α/(°)	d_max/m	N_max/g	Δθ/(°)
60.1	18.8	2.0	0.0245	25.8	0.0015

初始变量			目标函数
v/ (m·s^-1)	β/(°)	α/(°)	d_max/m	N_max/g	Δθ/(°)
60.1	18.8	2.0	0.0245	25.8	0.0015