回转体高速倾斜入水的流场特性及结构响应

高英杰; 孙铁志; 张桂勇; 尤天庆; 殷志宏; 宗智

doi:10.11883/bzycj-2020-0014

回转体高速倾斜入水的流场特性及结构响应

doi: 10.11883/bzycj-2020-0014

高英杰^1,,
孙铁志^{1, 2, ,},
张桂勇^{1, 2, 3},
尤天庆⁴,
殷志宏¹,
宗智^{1, 2, 3}

1.
大连理工大学船舶工程学院辽宁省深海浮动结构工程实验室，辽宁大连 116024
2.
大连理工大学工业装备结构分析国家重点实验室，辽宁大连 116024
3.
高新船舶与深海开发装备协同创新中心，上海 200240
4.
北京宇航系统工程研究所，北京 100076

基金项目: 国家自然科学基金（51639003，51709042）；中央高校基本科研业务费专项（DUT2017TB05）；中国博士后科学基金（2018M631791，2019T120211）；工信部高技术船舶科研项目（2017-614）；辽宁省自然科学基金（20180550619）；辽宁省“兴辽英才计划”（XLYC1908027）；海洋工程国家重点实验室开放基金（1803）

详细信息

作者简介:
高英杰（1994－　），男，硕士研究生，yingjie_gao@foxmail.com

通讯作者:
孙铁志（1986－　），男，博士，副教授，suntiezhi@dlut.edu.cn

中图分类号: O352
计量
- 文章访问数: 1456
- HTML全文浏览量: 774
- PDF下载量: 93
- 被引次数: 15
出版历程
- 收稿日期: 2020-01-13
- 修回日期: 2020-05-27
- 刊出日期: 2020-12-05

Flow characteristics and structure response of high-speed oblique water-entry for a revolution body

GAO Yingjie^1
,,
SUN Tiezhi^{1, 2
, ,},
ZHANG Guiyong^{1, 2, 3},
YOU Tianqing⁴,
YIN Zhihong¹,
ZONG Zhi^{1, 2, 3}

1.
Liaoning Engineering Laboratory for Deep-Sea Floating Structures, School of Naval Architecture,Dalian University of Technology, Dalian 116024, Liaoning, China
2.
State Key Laboratory of Structural Analysis for Industrial Equipment,Dalian University of Technology, Dalian 116024, Liaoning, China
3.
Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration, Shanghai 200240, China
4.
Beijing Aerospace System Engineering Research Institute, Beijing 100076, China

摘要

摘要: 回转体高速入水过程涉及液体和固体的耦合作用，是一个复杂的非线性、非定常过程。为研究回转体高速入水的结构动响应及流场演变规律，本文中基于STAR-CCM+和ABAQUS平台，建立了回转体高速入水的双向流固耦合数值模型，开展了不同入水速度的回转体高速倾斜入水流固耦合数值计算。结果表明：数值计算的入水速度、位移曲线和空泡形态与实验结果良好吻合，验证了流固耦合方法的有效性；回转体倾斜高速入水的载荷先集中在触水部分边缘处，后集中于回转体底部中心处；流固耦合方法的入水冲击载荷峰值小于刚体的，弹性回转体的载荷曲线产生明显波动；撞水阶段，回转体空泡呈现不对称形态，随着入水加深，空泡不对称性变弱；入水速度60 m/s下，空泡发生表面闭合，回转体入水初速度越快，空泡表面闭合越晚；冲击载荷与入水速度有关，入水速度越大，峰值出现越早，震荡越明显，速度超过100 m/s时，回转体产生塑性形变。
- 回转体 /
- 高速入水 /
- 流固耦合 /
- 流场 /
- 载荷
Abstract: The high-speed water-entry process of a revolution body involves fluid and structure interaction, which is a complex nonlinear and unsteady process. The revolution body with high speed is subjected to extreme impact load at the moment of hitting water surface which would cause great deformation or even damage to the structure. In order to investigate the physical understanding of the structural strength of revolution bodies and the mechanism of the cavity dynamics during high-speed water-entry process, a fluid-structure interaction (FSI) model based on a co-simulation progress between STAR-CCM+ and ABAQUS was adopted. The model performed a two-way interaction analysis which can consider the influence of the deformation of structure to fluid into. In the form of CFD analysis, a three-dimensional simulation with a six-degree-of-freedom model was carried out, in which the Shear Stress Transfer (SST) turbulence model and the volume of fluid (VOF) technique were used for turbulence computation and air-liquid interface tracking, respectively. In the part of FEM research, the shell mesh form with a whole Johnson-Cook material model was implemented to give a full consideration of deformation process and the accuracy of structure computation, which can effectively reflect the plastic deformation of structure. Firstly, a comparison between FSI result and experimental result of the high-speed water-entry process was conducted. The results show that the velocity attenuation, displacement and the cavity features are in good agreement with the experimental result, which proves two-way FSI method can be effectively applied into the research of high-speed water-entry problem. Then a numerical simulation of revolution body oblique water-entry with different velocity was carried out. The results show that the stress initially focuses on the edge of the bottom side of the revolution body, then it transports to the central area, remaining steady in the end. Compared with the results of rigid body, the peak value of impact load of the flexible body appears smaller due to the repeated deformation for buffer, which also causes the fluctuation of the load curve. After the initial water impacting, the cavity presents an asymmetric shape. As water-entry time increases, the asymmetry of cavity becomes weaker. In the process of 60 m/s oblique water-entry, surface closure of the cavity occurs. With the increase of water-entry velocity, the time of cavity surface closure takes longer. The peak value of the impact load whose period is quite short appears immediately at very beginning of water-entry process. After entering the water surface, the impact load of the revolution body decreases dramatically and rapidly and fluctuates slightly. The peak value of the impact load is related to water-entry velocity. The higher the velocity is, the earlier peak value of the impact load occurs and more obviously it fluctuates. As water-entry velocity exceeds 100 m/s, plastic deformation appears in the central area of bottom of revolution body.
- revolution body /
- high-speed water entry /
- fluid-structure interaction /
- fluid field /
- load

HTML全文

镁合金是一种潜力巨大的轻质工程材料，有着广阔的应用前景，其中MB2是镁合金的典型代表。材料在复杂应力状态下的失效破坏行为是工程材料和结构设计的理论基础。目前，镁合金在动静态加载下的力学行为研究以单向应力状态为主^[1-3]，且能实现复杂应力动态加载的实验设备与技术尚不成熟。J.W.Hancock等^[4]研究发现金属的延性明显依赖于应力的三轴状态。近年来，O.S.Hopperstad等^[5-6]对结构钢进行了拉伸加载实验和数值分析，结果表明材料的等效破坏应变随应力三轴状态的变化趋势与J.W.Hancock等的研究结论一致；Y.B.Bao^[7]在单轴加载条件下对铝合金2024进行了不同应力状态的实验研究；D.Anderson等^[8]研究了应力三轴度和应变率对DP780钢破坏行为的影响。汤安民等^[9]、李智慧等^[10]研究了不同金属材料宏观断裂形式及断裂机理与应力状态的关系, 陈刚等^[11]建立了应力三轴度和应变率相关的45钢损伤失效模型, 朱浩等^[12]分析了应力三轴度和应变率对铝合金6063力学性能的影响, 张伟等^[13]以应力三轴度和应变率效应为基础给出了铝合金7A04的本构关系和失效模型。

所以，对于钢材、铝合金等常用金属材料，其失效破坏行为的研究已经很多，而对镁合金复杂应力状态下的失效破坏行为研究还很少。本文中将通过复杂应力下的动静态拉伸实验及数值模拟，分析应力状态和应变率对镁合金MB2失效破坏的影响，基于Johnson-Cook本构及失效模型建立适用于镁合金MB2的破坏准则，分析其宏观破坏模式及其微观失效机理，为镁合金MB2在实际工程中的应用提供强度设计依据和理论支撑。

1. 动静态拉伸实验方法

材料所受应力状态不同时，材料内产生的塑性变形与应力集中程度将不同，为了表征材料的复杂应力状态，引用应力三轴度η为应力状态参数，即平均应力与等效应力的比值σ_H/σ。在延性金属复杂应力状态的失效破坏研究中，依据Bridgman原理^[14]设计的缺口试件拉伸实验是最主要的研究方法，根据圆弧缺口颈部应力方程解，可以得到应力三轴度的计算公式：

$\eta =\frac{1}{3}+\text{ln}\left( 1+\frac{{{a}^{2}}-{{r}^{2}}}{2aR} \right)$

(1)

式中：a和R分别为最小横截面的半径及缺口半径，r为到横截面中心的距离，a₀为最小横截面初始半径。基于塑性不可压的假设，可以得到缺口处的等效应变定义公式如下：

$\bar{\varepsilon }=2\text{ln(}{{a}_{0}}/a)$

(2)

当a取试件断裂时的横截面直径a_f时，式(2)计算得到的即为试件的等效破坏应变ε_f。静态拉伸实验试件采用直径为5 mm的光滑圆柱以及最小横截面直径为6 mm，缺口曲率半径分别为9、6、3 mm的缺口圆柱试件。动态拉伸实验试件采用直径为3 mm的光滑圆柱以及最小横截面直径为2 mm，缺口曲率半径分别为1.0、1.5、2.0、3.0 mm的缺口圆柱试件。静态拉伸加载由材料试验机MTS 810完成，而动态拉伸加载则通过分离式Hopkinson拉杆(SHTB)实验装置完成，如图 1所示。

图 1 SHTB实验系统

Figure 1. SHTB testing system

下载: 全尺寸图片幻灯片

2. 镁合金MB2的动静态本构关系

2.1 应变率及绝热效应

对于光滑圆柱试件的动态拉伸实验，按照一维应力波理论进行数据处理分析可知，试件在不同应力幅值的入射波作用下表现出了明显的应变率效应，不同应变率下的真实应力应变曲线如图 2所示。

图 2 镁合金MB2真实应力应变曲线

Figure 2. True stress-strain curves of magnesium alloy MB2

下载: 全尺寸图片幻灯片

在高速冲击载荷作用下，加载过程为绝热过程，材料在该过程中的温度变化ΔT可由塑性变形进行求解：

$\Delta T=\frac{\beta }{\rho {{c}_{p}}}\int_{0}^{{{\varepsilon }_{\text{p}}}}{\sigma \text{d}{{\varepsilon }_{\text{p}}}}$

(3)

式中：β为摩擦能量转换系数，ρ为材料密度，c_p为比定压热容，ε_p为等效塑性应变。对于镁合金MB2，摩擦能量转换系数通常取为0.925，密度为1 800 kg/m³，比定压热容为1 040~1 148 J/(kg·K)。

2.2 本构模型参数拟合

鉴于镁合金MB2的应变率效应，这里引用经典的J-C本构模型对实验结果进行拟合。在室温、准静态(0.001 s^-1)拉伸条件下，J-C本构模型的应变率项和温度项均为1，此时材料的本构模型退化为：

$\sigma =A+B\varepsilon _{\text{p}}^{n}$

(4)

采用最小二乘法对参数A、B、n进行拟合可以得到：A=192 MPa、B=218.3 MPa、n=0.370 56。由于在300 s^-1下试件未被拉断，所以这里以800和1 800 s^-1应变率下的抗拉强度及其对应的等效塑性应变作为比较点与参考应变率0.001 s^-1进行比较分析，获得参数C的算术平均值为0.015。对于反映温度效应的参数m，由于本研究未涉及高温下的动态拉伸实验，所以此处引用文献^[15]中的研究结论，给出参数m的值为0.95。

3. 失效破坏行为

3.1 宏观失效破坏

对于光滑圆柱试件和缺口圆柱试件，由于其加载过程中应力状态的差异，导致其动静态加载条件下的破坏形式也有所不同，如图 3所示。

图 3 不同试件破坏模式

Figure 3. Macro fracture patterns of different specimens

下载: 全尺寸图片幻灯片

光滑圆柱试件在整体上发生剪切破坏，且在破坏前有一定的颈缩，试件断口处有一定的杯锥特征。缺口试件在准静态条件下整体上发生正拉断，断口边缘有杯锥特征；在动态加载下整体表现为杯锥形断裂。通过测量断裂后缺口处的最小横截面直径，得到了动态拉伸试件不同加载条件下的等效破坏应变，如表 1所示。

表 1 不同类型试件的等效破坏应变

Table 1. Equivalent fracture strain of different specimens

入射波幅值/MPa	等效破坏应变
入射波幅值/MPa	光滑	R=3.0 mm	R=2.0 mm	R=1.5 mm	R=1.0 mm
50	未断裂	0.333 8	0.250 2	0.254 1	0.205 3
120	0.311 7	0.307 1	0.225 7	0.239 7	0.197 7
350	0.289 3	0.249 2	0.180 0	0.175 8	0.177 6

下载: 导出CSV

| 显示表格

3.2 J-C失效模型

J-C失效模型类似于J-C本构模型，采用多项乘积的形式对应力状态参数、应变率参数以及温度效应参数进行解耦。由于在本研究中未涉及等效破坏应变的温度效应，所以此处给出简化后的表达形式如下：

${{{\bar{\varepsilon }}}_{\text{f}}}=[{{D}_{1}}+{{D}_{2}}\text{exp(}{{D}_{3}}\eta )][1+{{D}_{4}}\text{ln(}\dot{\varepsilon }/{{{\dot{\varepsilon }}}_{0}})]$

(5)

式中：D₁、D₂、D₃、D₄为失效模型系数。研究发现，试件最小横截面处的应力三轴度是随等效应变的增加而不断变化的，整个加载过程的应力三轴度计算需要考虑到应变的累积效应^[16]，所以D.M.Goto等^[17]、Y.Bao等^[18]在研究不同延性金属材料时，对整个加载过程的平均应力三轴度η_av进行了定义：

${{\eta }_{\text{av}}}=\frac{1}{{{{\bar{\varepsilon }}}_{\text{f}}}}\int_{0}^{{{{\bar{\varepsilon }}}_{\text{f}}}}{\eta (\bar{\varepsilon })\text{d}\bar{\varepsilon }}$

(6)

采用ABAQUS/Explicit分析软件对镁合金MB2不同类型试件的动态拉伸实验过程进行模拟计算，鉴于SHTB和试件的对称性，以轴对称单元进行建模，以拉伸本构关系和实验中采集的入射波为输入条件，模拟得到了试件在加载过程中应力三轴度的变化情况，如图 4所示。

图 4 不同类型试件应力三轴度变化

Figure 4. Stress triaxiality curves of different specimens

下载: 全尺寸图片幻灯片

同时，由于缺口试件动态加载过程中，其塑性变形集中于缺口区域，并不满足一维应力波原理，无法通过理论公式对应变率进行求解。所以这里根据缺口试件等效应变的定义，借助数值模拟得到最小横截面直径的变化和应力波作用时程，通过时间积分得到试件的平均应变率：

${{{\dot{\varepsilon }}}_{\text{av}}}=\frac{1}{\Delta t}\int_{{{t}_{0}}}^{{{t}_{\text{f}}}}{\dot{\varepsilon }(t)\text{d}t=\frac{2}{\Delta t}}\text{ln}\left( \frac{{{a}_{0}}}{{{a}_{\text{f}}}} \right)$

(7)

式中：Δt为试件发生拉伸变形的时间，t₀为载荷作用开始的时刻，t_f为载荷作用结束的时刻。通过数值计算得到不同类型试件入射波在试件上的作用时间，进而计算得到其平均应变率，如表 2所示。

表 2 缺口试件的平均应变率

Table 2. Average strain rates of different specimens

入射波幅值/MPa	平均应变率/s^-1
入射波幅值/MPa	R=3 mm	R=2 mm	R=1.5 mm	R=1 mm
50	828	968	1 047	1 239
120	2 315	2 573	2 569	2 694
350	3 905	4 129	4 914	4 981

下载: 导出CSV

| 显示表格

在准静态( $\dot{\varepsilon }={{\dot{\varepsilon }}_{0}}$ )拉伸条件下，J-C失效模型的应变率项为1，且镁合金MB2的等效破坏应变存在明显的转折点^[19]，通过拟合得到，应力三轴度小于0.66时，参数D₁、D₂、D₃的值分别为1.04×10^-5、0.073 7以及0.323；当应力三轴度大于0.66时，其值分别为105.7、-0.103以及0.236。对于镁合金MB2，J-C失效模型中的(ε_f/ε₀-1)与ln( $\dot{\varepsilon }/{{\dot{\varepsilon }}_{0}}$ )并不满足线性关系，所以这里对模型中的应变率项进行了修正，修正后的失效模型为：

${{{\bar{\varepsilon }}}_{\text{f}}}=[{{D}_{1}}+{{D}_{2}}\text{exp(}{{D}_{3}}\eta )]\{1+{{D}_{5}}\text{ }\!\![\!\!\text{ 1-exp(}{{D}_{6}}\text{ln(}\dot{\varepsilon }/{{{\dot{\varepsilon }}}_{0}}))]\}$

(8)

式中：D₅、D₆为应变率项参数。通过拟合分析发现，修正模型中(ε_f/ε₀-1)与ln( $\dot{\varepsilon }/{{\dot{\varepsilon }}_{0}}$ )的关系与实验数据吻合，D₅、D₆的值分别为8.87×10^-20、2.801，如图 5所示。

图 5 应变率项参数拟合曲线

Figure 5. Fitted curve for the item of strain rates

下载: 全尺寸图片幻灯片

结合参数拟合的结果，可以得到镁合金MB2最终的破坏准则为以应力三轴度0.66为转折点的分段函数形式，如图 6所示。图中实心圆点为实验结果，曲面为破坏准则的拟合结果。

图 6 镁合金MB2的破坏准则曲面

Figure 6. Fracture criterion surface of magnesium alloy MB2

下载: 全尺寸图片幻灯片

通过失效破坏曲面可以发现，镁合金MB2的断裂延性(等效破坏应变)随着应力三轴度的增加先增大后减小，随着应变率的增大而不断减小。从破坏模式上来看，随着应力三轴度的增大，材料先发生剪切破坏，随后又发生正拉断破坏；随着应变率的增大，光滑圆柱及缺口圆柱试件在动静态载荷作用下的破坏形式基本一致。所以，应力状态对断裂延性和破坏模式的影响存在明显的转折点，而应变率的影响则不存在转折点。

为了对该失效模型进行验证，这里通过ABAQUS/Explicit分析软件和单元失效法对不同类型试件的动态失效破坏形式进行了数值模拟分析，其模拟结果如图 7所示。

图 7 不同试件失效破坏模式的数值模拟结果

Figure 7. Fracture patterns of different specimens from numerical simulation

下载: 全尺寸图片幻灯片

从数值模拟结果可以看出，无论是光滑圆柱试件还是缺口圆柱试件，其最终破坏形式均呈现出杯锥形特征，这与实验现象是基本吻合的；但对于光滑圆柱试件，该模型不能体现试件的整体剪切断裂特征，这与J-C失效模型的物理涵义以及镁合金MB2在不同应力状态下的微观损伤机理有关。

3.3 微观失效机理

宏观破坏行为与材料的微观损伤变形机制密切相关，为进一步认识应力状态对镁合金MB2破坏行为的影响，对准静态条件下光滑圆柱及缺口试件的断口形貌进行微观扫描观察，如图 8所示。

图 8 试件微观断口形貌

Figure 8. Micro fracture morphology of specimens

下载: 全尺寸图片幻灯片

对于光滑拉伸试件，微观断口上能看到河流花样特征即解理断裂，同时也呈现出了一定的微孔洞特征，即有韧性断裂的特点，如图 8(a)所示，属韧脆性混合断裂；对于缺口试件，不同缺口程度均表现出较明显的微孔洞断裂的特征，在图 8(b)中能清晰地看到韧窝聚合、连接的特点。这些特征与试件在宏观破坏形式上的表现是一致的，即光滑圆柱试件破坏是由混合损伤机制引起的，而缺口试件破坏则是由微孔洞损伤演化机制决定的。J-C失效模型在微观上反映的是材料的微孔洞损伤演化机理，在宏观上体现为该微观损伤机理导致的材料断裂延性变化，不能体现其他变形机制的作用，这就是模拟结果可以描述光滑圆柱试件的杯锥形破坏特点而无法描述其整体剪切特征的原因。

4. 结论

(1) 通过缺口圆柱试件设计能够实现镁合金MB2不同程度的多向应力状态，结合SHTB装置可以实现其在复杂应力状态下的动态拉伸加载研究。

(2) 在研究的实验条件范围内，应力状态对镁合金MB2破坏行为的影响存在明显的转折点，而应变率(5 000 s^-1以内)则不存在转折点。通过模型修正和参数拟合得到了镁合金MB2的等效破坏应变准则，且其数值模拟结果可以较好地体现拉伸试件的杯锥形破坏特征。

(3) 镁合金MB2不同应力状态下的破坏行为与其微观损伤机理密切相关，随着应力三轴度的增加，其微观损伤机制由混合断裂转变为韧窝断裂，这是导致其宏观延性和破坏模式发生变化的重要原因。

图 1 双向耦合求解过程

Figure 1. Computing process of fluid-structure interaction method

下载: 全尺寸图片幻灯片

图 2 计算域

Figure 2. Computational domain

下载: 全尺寸图片幻灯片

图 3 流体网格模型

Figure 3. Mesh of fluid domain

下载: 全尺寸图片幻灯片

图 4 速度衰减曲线

Figure 4. Velocity attenuation curves

下载: 全尺寸图片幻灯片

图 5 侵蚀位移曲线

Figure 5. Erosion displacement curves

下载: 全尺寸图片幻灯片

图 6 不同时刻的空泡形态

Figure 6. Cavity features at different times

下载: 全尺寸图片幻灯片

图 7 回转体轴向、底部径向的单元分布

Figure 7. Element distribution of revolution body

下载: 全尺寸图片幻灯片

图 8 回转体的应力沿轴向分布曲线

Figure 8. Equivalent stress curves of elements along axis

下载: 全尺寸图片幻灯片

图 9 不同时刻回转体底部应力分布

Figure 9. Stress distributions of revolution body bottom at different times

下载: 全尺寸图片幻灯片

图 10 回转体底部的应力沿径向分布曲线

Figure 10. Equivalent stress curves of bottom elements along radius

下载: 全尺寸图片幻灯片

图 11 回转体底部的应变沿径向分布曲线

Figure 11. Strain curves of bottom elements along radius

下载: 全尺寸图片幻灯片

图 12 不同时刻的空泡演变

Figure 12. Cavity features at different times

下载: 全尺寸图片幻灯片

图 13 不同时刻的应力分布

Figure 13. Stress distributions at different times

下载: 全尺寸图片幻灯片

图 14 速度60 m/s倾斜入水的空泡演变

Figure 14. Cavity features for oblique water-entry at v₀=60 m/s

下载: 全尺寸图片幻灯片

图 15 入水速度80 m/s时回转体底部中心点压力曲线

Figure 15. Pressure curve of bottom central point at v₀=80 m/s

下载: 全尺寸图片幻灯片

图 16 速度80 m/s入水时速度曲线

Figure 16. Velocity curves at v₀=80 m/s

下载: 全尺寸图片幻灯片

图 17 归一化的速度衰减曲线

Figure 17. Normalized velocity attenuation curves

下载: 全尺寸图片幻灯片

图 18 回转体中心单元的有效应力曲线

Figure 18. Equivalent stress curves of central elements

下载: 全尺寸图片幻灯片

图 19 空泡形态对比

Figure 19. Comparisons of cavity features

下载: 全尺寸图片幻灯片

表 1 不同入水速度的压力峰值

Table 1. Peak pressures at different initial water entry velocities

入水速度/(m·s⁻¹) 压力峰值/MPa
CFD FSI

60 9.93 9.42
80 16.03 15.59
100 27.92 25.12

下载: 导出CSV

参考文献(24)

[1]	王永虎, 石秀华. 入水冲击问题研究的现状与进展 [J]. 爆炸与冲击, 2008, 28(3): 276–282. DOI: 10.11883/1001-1455(2008)03-0276-07. WANG Y H, SHI X H. Review on research and development of water-entry impact problem [J]. Explosion and Shock Waves, 2008, 28(3): 276–282. DOI: 10.11883/1001-1455(2008)03-0276-07.
[2]	Von KARMAN T H. The impact on seaplane floats during landing [R]. Washington DC: National Advisory Committee on Aeronautics, 1929: 309−313. DOI: 10.1115/1.4023571.
[3]	WAGNER H. Über stoß-und gleitvorgänge an der oberfläche von flüssigkeiten [J]. Zeitschrift für Angewandte Mathematik und Mechanik, 1932, 12(4): 193–215. DOI: 10.1002/zamm.19320120402.
[4]	GAVRILENKO V N, KUBENKO V D. Plane problem of rigid body penetration into a compressible fluid [J]. Soviet Applied Mechanics, 1985, 21(4): 345–352. DOI: 10.1007/BF00886581.
[5]	WORTHINGTON A M, COLE R S. Impact with a liquid surface studied by the aid of instantaneous photography. Paper Ⅱ [J]. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 1900, 194: 175–199. DOI: 10.1098/rsta.1900.0016.
[6]	MAY A. Effect of surface condition of a sphere on its water-entry cavity [J]. Journal of Applied Physics, 1951, 22(10): 1219–1222. DOI: 10.1063/1.1699831.
[7]	MAY A, WOODHULL J C. The virtual mass of a sphere entering water vertically [J]. Journal of Applied Physics, 1950, 21(12): 1285–1289. DOI: 10.1063/1.1699592.
[8]	ABELSON H I. Pressure measurements in the water-entry cavity [J]. Journal of Fluid Mechanics, 1970, 44(1): 129–144. DOI: 10.1017/S0022112070001738.
[9]	PANAHI R. Simulation of water-entry and water-exit problems using a moving mesh algorithm [J]. Journal of Theoretical and Applied Mechanics, 2012, 42(2): 79–92.
[10]	ERFANIAN M R, ANBARSOOZ M, RAHIMI N, et al. Numerical and experimental investigation of a three dimensional spherical-nose projectile water entry problem [J]. Ocean Engineering, 2015, 104: 397–404. DOI: 10.1016/j.oceaneng.2015.05.024.
[11]	NEAVES M D, EDWARDS J R. All-speed time-accurate underwater projectile calculations using a preconditioning algorithm [J]. Journal of Fluids Engineering, 2006, 128(2): 284–296. DOI: 10.1115/1.2169816.
[12]	孙辉, 卢炽华, 何友声. 二维楔形体冲击入水时的流固耦合响应的实验研究 [J]. 水动力学研究与进展, 2003, 18A(1): 104–109. DOI: 10.3969/j.issn.1000-4874.2003.01.018. SUN H, LU C H, HE Y S. Experimental research on the fluid-structure interaction in water entry of 2D elastic wedge [J]. Journal of Hydrodynamics, 2003, 18A(1): 104–109. DOI: 10.3969/j.issn.1000-4874.2003.01.018.
[13]	张伟, 郭子涛, 肖新科, 等. 弹体高速入水特性实验研究 [J]. 爆炸与冲击, 2011, 31(6): 579–584. DOI: 10.11883/1001-1455(2011)06-0579-06. ZHANG W, GUO Z T, XIAO X K, et al. Experimental investigations on behaviors of projectile high-speed water entry [J]. Explosion and Shock Waves, 2011, 31(6): 579–584. DOI: 10.11883/1001-1455(2011)06-0579-06.
[14]	郭子涛, 张伟, 郭钊, 等. 截卵形弹水平入水的速度衰减及空泡扩展特性 [J]. 爆炸与冲击, 2017, 37(4): 727–733. DOI: 10.11883/1001-1455(2017)04-0727-07. GUO Z T, ZHANG W, GUO Z, et al. Characteristics of velocity attenuation and cavity expansion induced by horizontal water-entry of truncated-ogive nosed projectiles [J]. Explosion and Shock Waves, 2017, 37(4): 727–733. DOI: 10.11883/1001-1455(2017)04-0727-07.
[15]	陈诚, 袁绪龙, 党建军, 等. 超空泡航行器20°角倾斜入水冲击载荷特性试验研究 [J]. 兵工学报, 2018, 39(6): 1159–1164. DOI: 10.3969/j.issn.1000-1093.2018.06.016. CHEN C, YUAN X L, DANG J J, et al. Experimental investigation into impact load during oblique water-entry of a supercavitating vehicle at 20° [J]. Acta Armamentarii, 2018, 39(6): 1159–1164. DOI: 10.3969/j.issn.1000-1093.2018.06.016.
[16]	YAN G X, PAN G, SHI Y, et al. Experimental and numerical investigation of water impact on air-launched AUVs [J]. Ocean Engineering, 2018, 167: 156–168. DOI: 10.1016/j.oceaneng.2018.08.044.
[17]	钱铖铖, 余春华, 穆青, 等. 发射速度和发射角度对射弹高速入水流动的影响 [J]. 兵器装备工程学报, 2019, 40(7): 35–39; 50. DOI: 10.11809/bqzbgcxb2019.07.008. QIAN C C, YU C H, MU Q, et al. Numerical research of effects of launch speed and launch angle on water entry of high-speed projectile [J]. Journal of Ordnance Equipment Engineering, 2019, 40(7): 35–39; 50. DOI: 10.11809/bqzbgcxb2019.07.008.
[18]	张佳悦, 李达钦, 吴钦, 等. 航行体回收垂直入水空泡流场及水动力特性研究 [J]. 力学学报, 2019, 51(3): 803–812. DOI: 10.6052/0459-1879-18-364. ZHANG J Y, LI D Q, WU Q, et al. Numerical investigation on cavity structures and hyrodynamics of the vehicle during vertical water-entry [J]. Chinese Journal of Theoretical and Applied Mechanics, 2019, 51(3): 803–812. DOI: 10.6052/0459-1879-18-364.
[19]	何春涛, 王聪, 闵景新, 等. 回转体匀速垂直入水早期空泡数值模拟研究 [J]. 工程力学, 2012, 29(4): 237–243. HE C T, WANG C, MIN J X, et al. Numerical simulation of early air-cavity of cylinder cone with vertical water-entry [J]. Engineering Mechanics, 2012, 29(4): 237–243.
[20]	侯昭, 孙铁志, 张桂勇, 等. 回转体倾斜入水空泡试验及六自由度数值计算研究 [J]. 宇航总体技术, 2017, 1(4): 38–45. HOU Z, SUN T Z, ZHANG G Y, et al. Experimental investigation and 6-DOF simulation of oblique water-entry cavity of revolution body [J]. Astronautical Systems Engineering Technology, 2017, 1(4): 38–45.
[21]	黄志刚, 孙铁志, 杨碧野, 等. 平头锥型回转体高速入水结构强度数值分析 [J]. 爆炸与冲击, 2019, 39(4): 043201. DOI: 10.11883/bzycj-2017-0330. HUANG Z G, SUN T Z, YANG B Y, et al. Numerical analysis on structural strength of a cone-shaped flatted revolution body during high-speed water-entry [J]. Explosion and Shock Waves, 2019, 39(4): 043201. DOI: 10.11883/bzycj-2017-0330.
[22]	陈宇, 曹红松, 王智军, 等. 一种新型防空MEFP的设计与仿真 [J]. 现代防御技术, 2013, 41(2): 202–206; 211. DOI: 10.3969/j.issn.1009-086x.2013.02.037. CHEN Y, CAO H S, WANG Z J, et al. New type of air defense MEFP [J]. Modern Defense Technology, 2013, 41(2): 202–206; 211. DOI: 10.3969/j.issn.1009-086x.2013.02.037.
[23]	AGMELL M, AHADI A, STÅHL J E. A numerical and experimental investigation of the deformation zones and the corresponding cutting forces in orthogonal cutting [J]. Advanced Materials Research, 2011, 223: 152–161. DOI: 10.4028/www.scientific.net/AMR.223.152.
[24]	CHEN T, HUANG W, ZHANG W, et al. Experimental investigation on trajectory stability of high-speed water entry projectiles [J]. Ocean Engineering, 2019, 175: 16–24. DOI: 10.1016/j.oceaneng.2019.02.021.

施引文献

期刊类型引用(4)

1.	祁武超，刘炳宏，田素梅. 镁合金带解锁装置热冲击断裂特性及拓扑优化. 兵器装备工程学报. 2023(11): 47-56 . 百度学术
2.	周伦，苏兴亚，敬霖，邓贵德，赵隆茂. 6061-T6铝合金动态拉伸本构关系及失效行为. 爆炸与冲击. 2022(09): 113-124 . 本站查看
3.	衣海娇，甄莹，曹宇光，张士华，史永晋. 6061-T6铝合金断裂应变与应力三轴度关系研究. 机械强度. 2020(03): 551-558 . 百度学术
4.	Dewang Zhao，Daxin Ren，Kunmin Zhao，Pan Sun，Xinglin Guo，Liming Liu. Ultrasonic Welding of Magnesium–Titanium Dissimilar Metals:A Study on Thermo-mechanical Analyses of Welding Process by Experimentation and Finite Element Method. Chinese Journal of Mechanical Engineering. 2019(06): 191-201 . 必应学术

其他类型引用(11)

资源附件(0)

访问统计

图(19) / 表(1)

计量

文章访问数: 1456
HTML全文浏览量: 774
PDF下载量: 93
被引次数: 15

1. 动静态拉伸实验方法
2. 镁合金MB2的动静态本构关系
2.1 应变率及绝热效应
2.2 本构模型参数拟合
3. 失效破坏行为
3.1 宏观失效破坏
3.2 J-C失效模型
3.3 微观失效机理
4. 结论

入水速度/(m·s⁻¹)	压力峰值/MPa
入水速度/(m·s⁻¹)	CFD	FSI
60	9.93	9.42
80	16.03	15.59
100	27.92	25.12

回转体高速倾斜入水的流场特性及结构响应

doi: 10.11883/bzycj-2020-0014

作者简介: 高英杰（1994－ ），男，硕士研究生，yingjie_gao@foxmail.com

通讯作者: 孙铁志（1986－ ），男，博士，副教授，suntiezhi@dlut.edu.cn

计量

出版历程