On ductile fracture of 30CrMnSiNi2A steel considering effects of stress triaxiality
-
摘要: 30CrMnSiNi2A钢是一种在军工领域应用广泛的低合金高强度钢。针对结构完整性的评估问题,采用试验和数值计算结合的方法研究了30CrMnSiNi2A钢的韧性断裂特性。对光滑圆棒试件在不同温度下进行准静态和动态拉伸试验,并通过有限元迭代方法标定了材料的Johnson-Cook动态本构模型参数,分析了温度和应变率对30CrMnSiNi2A钢断裂行为的影响。开展了缺口圆棒拉伸、缺口平板剪切和圆柱压缩试验,计算了各试件对应的平均应力三轴度和断裂应变,给出了应力三轴度在−1/3~1.5区间内的断裂应变变化曲线,分别确定了Johnson-Cook和Bao-Wierzbicki失效模型参数。研究表明,30CrMnSiNi2A钢的断裂应变与应力状态密切相关,且在不同的应力三轴度区间内曲线单调性差异较大,Bao-Wierzbicki失效模型较好地描述了这种钢在不同应力状态下的断裂特性。
-
关键词:
- 韧性断裂 /
- 动态力学性能 /
- 应力三轴度 /
- 30CrMnSiNi2A钢
Abstract: 30CrMnSiNi2A steel is a high-strength low-alloy steel and it is widely used in the defense industry. Its ductile fracture properties were investigated by using a combined experimental-numerical approach to address the need for assessment of structural integrity. The parameters of the Johnson-Cook constitutive model were obtained through quasi-static and dynamic tensile tests on round bar specimens at different temperatures and relevant finite element method iterations. The effects of strain rate and temperature on fracture were studied. The tensile tests on notched round bar specimens were performed to calibrate the fracture strain in the range of high positive stress triaxiality. The tensile tests on butterfly plates and compression tests on short cylindrical specimens covered the fracture properties in the range of low and negative stress triaxiality. The finite element models were computed and the fracture loci in the space of the effective plastic strain to fracture and the stress triaxiality in a wide range from −1/3 to 1.5 were constructed. The parameters of the Johnson-Cook fracture model and Bao-Wierzbicki fracture model were calibrated. It is shown that stress triaxiality has a significant effect on the fracture of 30CrMnSiNi2A steel and the monotonicity of the fracture loci varies in different stress triaxiality range. The Bao-Wierzbicki model is capable of predicting the fracture patterns of the 30CrMnSiNi2A steel in different stress states.-
Key words:
- ductile fracture /
- dynamic mechanical behavior /
- stress triaxiality /
- 30CrMnSiNi2A steel
-
随着混凝土含水率的增加, 混凝土的力学性能发生改变, 研究表明, 饱和状态下的混凝土结构在复杂受力状态下的力学性能与干燥状态下的差别较大[1-2]。探讨饱水混凝土的力学性能对于水下爆破开挖、水中建筑物拆除、大型水下建筑物抗震安全评价等有着重要的意义。
在水环境下混凝土类材料的力学性能研究方面, 目前已有一些学者开展了相关的研究, 但主要局限于准静态情况, 其动态力学性能的研究方面还存在不足:P.Rossi等[3-4]讨论了动态拉伸应力下混凝土中的自由水对混凝土的作用, 从机制上对湿态混凝土力学性能的变化进行解释, 但对湿态混凝土的抗压强度没有研究; P.K.Mehta等[5]和X.H.Vu等[6]研究了多种含水率和干燥混凝土的强度特性, 但仅描述了这一实验现象, 没有深入探讨水的作用机理; 王海龙等[7]对比了不同加载速率条件下的干燥和饱和混凝土抗压强度变化, 但其加载速率均小于10-4 s-1, 仍属于准静态加载范畴。
材料的力学性能在应变率比较高时通常会发生较大的变化, 此时建立材料的本构关系就需要考虑惯性力和应变率对材料性能的影响。目前应用于应变率在10-6~104 s-1的材料试验装置及其测量技术已趋成熟, 大口径分离式Hopkinson装置是实现对混凝土材料中等应变率加载的有效途径之一, 而饱和混凝土在中等应变率100~102 s-1条件下压缩加载的动态力学特性研究极少。本文中, 开展干燥和饱和混凝土的冲击压缩动态力学特性实验, 对2种状态混凝土在中等应变率条件下力学性能进行研究, 并进行理论上的分析探讨。
1. 实验
1.1 实验设备
静载单轴压缩实验在MTS810型材料试验机上完成。动态实验是在∅75 mm的大口径分离式Hopkinson装置上进行的, 如图 1所示。该装置的入射杆长5 461 mm, 透射杆长3 485 mm, 子弹长50 mm。输入杆和传送杆上贴有应变片用以记录波形。
SHPB方法建立在2个基本假定基础上:(1)杆中一维应力波假定, 即不考虑非轴向应力; (2)试件均匀假定, 即应力应变沿试件长度均匀分布。根据假设, 有:
εi(t)+εr(t)=εt(t) (1) 式中:εi(t)、εr(t)和εt(t)分别为入射波、反射波和透射波产生的应变。基于式(1)可得到如下关系:
˙ε=2c0Ls[εi(t)−εt(t)] (2) ε=2c0Ls∫t0[εi(t)−εt(t)]dt (3) σ=(A0/As)E0εt(t) (4) 式中:ε和σ为试件的应变和应力,
为应变率, c0是杆中弹性波速度, E0是杆的弹性模量, A0是杆的横截面积, As和Ls分别是试件的原始截面积和原始长度。
由实验得到输入杆和输出杆中的应变波形, 通过式(3)~(4)可求得材料在某一应变率下的应力应变曲线, 相应的应变率可由式(2)得到。
1.2 实验材料和方法
MTS实验所用试样为∅50 mm×100 mm的圆柱体, 两端表面平行度在0.05 mm以内, 表面平整度在0.02 mm以内, 试样如图 2(a)所示。SHPB实验所用试样为∅75 mm×37 mm的圆柱体, 其表面平行度和平整度与静态实验试样相同, 试样如图 2(b)所示。2个实验中的试样均分为2种状态:(1)干燥状态, 试样制备好自然风干, 然后开展实验; (2)饱和状态, 试样在水中自由吸水至饱和后开展实验。共开展干燥和饱和状态的准静态单轴压缩实验各5个, 以及不同速率的SHPB实验各12个。
准静态实验中加载速率为4×10-4 s-1。SHPB实验中, 将干燥和饱和混凝土分别以4个为一组进行分批, 通过调整气压控制加载速率, 依次进行加载速率从低到高的动态冲击压缩实验。
2. 实验结果及分析
2.1 静态力学特性
在准静态加载条件下, 对实验结果取平均值, 得到干燥和饱和混凝土的强度和峰值应变。实验结果见表 1, 其中:σ0为平均强度, ε0为平均峰值应变。表 1的结果表明:饱和混凝土的静态强度比干燥状态的静态强度低约37.8%, 而其峰值应变(峰值应力处的应变)小23.6%。
表 1 准静态条件下的实验结果Table 1. Results under quasi-static conditions状态 σ0/MPa ε0 /s-1 干燥 49.84 0.006 19 0.000 4 饱和 36.17 0.005 01 0.000 4 2.2 动态变形特性
选取不同应变率条件下干燥和饱和混凝土的实验数据, 绘制应力应变曲线, 如图 3~4所示, 其中:
应的曲线为准静态应力应变曲线。从实验结果可以看出, 干燥和饱和混凝土的动态应力应变曲线与准静态下的应力应变曲线有显著差别, 上升段都较陡, 表明动态加载条件下的弹性模量大于静态加载条件下的弹性模量, 且在中等应变率条件下的峰值应变均小于准静态条件下的峰值应变。
为进一步对比分析相近应变率下水对混凝土动态力学特性的影响, 结合实验的加载应变率范围, 选取应变率相近的试样数据进行分析, 得到的应力应变曲线如图 5所示。从图 5可以看出, 干燥和饱和混凝土的动态应力应变曲线的形状是基本相同的, 但在峰值应力之后, 饱和混凝土的曲线下降较缓, 这说明饱和混凝土具有一定的延性。
对有效SHPB实验试样的峰值应变εp进行统计, 得到它们的关系趋势线, 如图 6所示。在中等应变率条件下, 干燥和饱水混凝土的峰值应变均随着应变率的增大而增大, 并且饱和混凝土的峰值应变比干燥混凝土的峰值应变大, 而随着应变率的增加两者之间的差值逐渐减小。
2.3 动态强度特性
实验得到的干燥与饱和混凝土的强度σb与应变率
的关系如图 7所示。结合己有混凝土动态特性的应变率依赖特性研究成果[8-11], 以单轴压缩平均抗压强度σ0及其对应的应变率
0为归一化因子将应力和应变率进行量纲一化, 得到干燥与饱和混凝土的动态强度-应变率对数曲线, 如图 8所示。混凝土的干燥和饱和状态下的动态强度-应变率关系拟合结果如下:
σdσd0=0.139ln˙εd˙εd0−0.191 (5) σsσs0=0.242ln˙εs˙εs0−0.826 (6) 式中:σd和σs分别为干燥和饱和混凝土的动态强度, σd0和σs0分别为干燥和饱和混凝土在准静态条件下的单轴压缩平均抗压强度,
d和
s分别为干燥和饱和混凝土在冲击荷载作用下的应变率,
d0和
s0分别为干燥和饱和混凝土在准静态条件下的应变率。
根据图 7~8可知, 随着应变率的增大, 干燥和饱和混凝土强度均有提高的趋势, 其中, 干燥混凝土动态强度提高了32.2%, 饱和混凝土动态强度提高了72.4%, 增加幅度为干燥混凝土的近2倍。这表明:不论干燥或是饱水情况, 混凝土均具有明显的应变率效应, 并且饱和混凝土比干燥混凝土具有更强的应变率敏感性。
从图 7可以看出, 应变率存在一临界值
c(本实验中为(17.74±2) s-1), 当应变率
c时, 饱和混凝土的动态强度大于干燥混凝土的的动态强度, 即在低应变率条件下, 干燥混凝土强度比饱水混凝土大, 而在高应变率条件下, 则相反。李赞成等[10]在研究C10混凝土的动态强度时, 也发现在应变率100~500 s-1的条件下, 饱和混凝土的动态强度比干燥混凝土高25%~100%。
3. 水对混凝土强度影响机理探讨
饱和混凝土在微破裂稳定发展阶段会出现塑性变形和体胀, 这为孔隙中的自由水对混凝土的破坏作用提供有利条件。在准静态或应变率较低(小于
c)的情况下, 由于裂纹拓展的速率较小, 自由水有足够的时间沿裂隙扩散, 并在一定程度上润滑裂隙接触面, 从而促进裂纹拓展, 最终使混凝土破坏。根据目前的实验和理论研究可知, 混凝土是在裂纹尖端的拉应力作用下破坏的, 当应变率较大(大于
c)时, 裂纹动态扩展的速率很大, 混凝土中自由水无法在瞬间扩散到新形成的裂隙中, 水的作用反而会阻碍裂纹的扩展, 导致饱水混凝土在较高强度下破裂, 从而使饱和混凝土的应变率效应更明显。P.Rossi等[11-12]和D.Zheng等[13]认为在动态加载条件下湿混凝土强度提高的原因可由物理学中的Stefan效应来解释, 即假设在两半径为r的平行圆形平板中间有黏性液体存在, 当两平板以相对速度v分离时, 液体会产生一个黏聚力L来阻止平板的分离, 由Stefan效应产生的阻力L表示为:
L=3πμr42h3v (7) 式中:μ为液体的黏度; h为两圆形平板的间距。动态荷载作用下的应变率就是该模型分离速度的具体表征量。该黏聚力在饱和混凝土内部裂纹表面处产生的应力为:
σc=3μr22h3v (8) D.Zheng等[13]认为在低加载速率下考虑自由水的作用时, 动态增强因子会大于1.0(不考虑水的作用时动态增强因子为1.0)。
图 7的实验结果和其他实验研究表明[6-7-14]:对混凝土类脆性材料而言, 在准静态或低加载速率下, 材料在含水状态下的强度低于干燥状态下的强度。这与文献[11-13]中的观点矛盾, 主要原因是文献[11-13]中在考虑Stefan效应时没有考虑材料的应变率临界值, 只有当应变率足够大时, Stefan效应才能表现出来。
实际上, 水对混凝土的影响可以分为2部分(如图 9所示)。
(1) 对混凝土基体的影响。
经充分浸泡后, 基体由于水的作用被弱化, 饱和基体的弱化效应主要是受与基体结合水的影响, 若将弱化因子定义为D, D与混凝土的含水量有关, 物理意义上与抗压强度存在一定关系:
D=σs0σd0 (9) (2) 混凝土空隙或裂隙中的自由水对混凝土基体的力学作用。
饱和混凝土中的Stefan效应、楔入效应及自由水毛细作用均是由裂隙中的自由水作用产生的, 各效应产生的作用力如图 10所示。楔入效应主要是由在孔隙裂隙压密时所形成的孔隙水压力p产生的, 并且与裂纹的扩展速度有关[14-15]。在准静态或应变率较低(见图 10(a))时, 自由水能沿新形成的微裂隙继续扩散, 孔隙水压力p不断减小, 直至消失, 其影响作用包含在弱化因子中。在较高应变率条件下(见图 10(b)), 楔入效应的过程极其短暂, 由于裂纹的迅速扩展, 因表面张力和毛细作用pc对裂隙表面的吸附, 自由水不容易达到缝尖, 反而以相对速率v向裂隙内运动, 产生显著的Stefan效应。自由水因黏性产生的Stefan效应阻碍裂纹的扩展。
水的毛细作用pc的计算公式[15]为:
pc=4σ∗ρcosθ (10) 式中:σ*为界面张力; ρ为裂纹宽度(毛细管的直径); θ为接触角。
由图 10(b)可知, 毛细作用pc作用方向平行于劈裂面, 对裂纹的扩展基本没有影响, 其作用可以忽略。因此不同饱和度下混凝土的强度与应变率的关系为:
σb=f(˙ε)+[g(˙ε)−Df(˙ε)]φ(w) (11) 式中:
)为干燥状态下的应变率函数, 可由实验拟合式(5)求得;
)为Stefan效应L关于应变率的函数, 在冲击荷载产生的中等应变率条件下, 根据式(8)可求得
/(2h3), φ(w)是关于混凝土饱和度w的函数:
φ(w)={0w=0aw0<w<11w=1 (12) 式中:a为实验参数。在中等应变率条件下, 干燥和饱和混凝土的动态强度相等时:
f(˙ε)=1Dg(˙ε) (13) 此时得到的应变率即为临界应变率
c。当
c时, σs < σd。仅当
c时, 饱和混凝土的强度σs才大于干燥混凝土的强度σd, 混凝土中因自由水的黏性产生的Stefan效应才能良好地展现出来, 因此考虑材料中的Stefan效应时, 必须考虑应变率这个条件, 这也是目前各学者实验研究结果不一致的原因。根据式(11), 只需知道干燥混凝土的强度与应变率的关系就可预测出不同饱和度下的混凝土的强度与应变率的关系。若将该关系进一步量化, 还需要更深入的研究和探索。
4. 结论
(1) 在中等应变率条件下, 干燥和饱和混凝土的应力应变曲线的形状是基本相同的, 动态加载条件下的弹性模量大于静态加载条件下的弹性模量; 饱和混凝土的峰值应变比干燥混凝土的峰值应变大, 随着应变率的增加两者之间的差值逐渐减小。
(2) 干燥和饱和混凝土均具有明显的应变率效应, 随着应变率的提高, 饱和混凝土动态强度提高的幅度为干燥混凝土的近两倍, 前者具有更强的应变率敏感性。
(3) 存在应变率临界值
c:当
c时, 饱和混凝土的强度小于干燥混凝土的的强度; 而当
c时, 则相反。因此考虑材料中的Stefan效应时, 必须考虑应变率。
-
表 1 异形试件
Table 1. Abnormal specimens
序号 试件 载荷 数量 应力三轴度范围 a R4缺口圆棒 拉伸 3 [1/3,+∞] b R2缺口圆棒 拉伸 3 [1/3,+∞] c 45°缺口平板 拉伸 4 [0,1/3] d 0°缺口平板 拉伸 4 [0,1/3] e 短粗圆柱 压缩 3 [−1/3,0] f 细长圆柱 压缩 3 [−1/3,0] 表 2 30CrMnSiNi2A钢的化学成分(质量分数,%)
Table 2. Composition of 30CrMnSiNi2A steel (mass fraction, %)
C Cr Mn Si Ni Cu P S 0.300 1.030 1.160 1.040 1.680 0.080 0.014 0.020 表 3 光滑圆棒试件的屈服应力和断裂应变
Table 3. Yield strength and fracture strain of round bars
序号 温度/K 应变率/s−1 屈服应力/MPa 断裂应变 1 298 0.001 1 163 0.753 2 298 1 000 1 591 0.775 3 298 2 000 1 716 0.881 4 298 4 000 1 883 0.962 5 473 0.001 1 115 0.827 6 773 0.001 698 1.513 表 4 平均应力三轴度和断裂应变对应表
Table 4. A summary of average stress triaxiality and equivalent plastic fracture strain
序号 应力三轴度理论值 平均应力三轴度 断裂应变 1 −1/3 −0.246 3 1.203 2 −1/3 −0.097 2 0.437 3 0 0.013 6 0.281 4 1/6 0.086 6 0.455 5 1/3 0.609 4 0.753 6 0.652 1.069 0 0.382 7 0.893 1.353 0 0.337 -
[1] 武海军, 姚伟, 黄风雷, 等. 超高强度钢30CrMnSiNi2A动态力学性能实验研究 [J]. 北京理工大学学报, 2010, 30(3): 258–262. DOI: 10.15918/j.tbit1001-0645.2010.03.002.WU H J, YAO W, HUANG F L, et al. Experimental study on dynamic mechanical properties of ultrahigh strength 30CrMnSiNi2A steel [J]. Transactions of Beijing Institute of Technology, 2010, 30(3): 258–262. DOI: 10.15918/j.tbit1001-0645.2010.03.002. [2] 周义清, 张治民. 30CrMnSiNi2A钢在不同应变率下的力学性能研究 [J]. 兵器材料科学与工程, 2010, 33(4): 46–50. DOI: 10.14024/j.cnki.1004-244x.2010.04.021.ZHOU Y Q, ZHANG Z M. Mechanical property of 30CrMnSiNi2A steel under various strain rates [J]. Ordnance Material Science and Engineering, 2010, 33(4): 46–50. DOI: 10.14024/j.cnki.1004-244x.2010.04.021. [3] 李磊, 张先锋, 吴雪, 等. 不同硬度30CrMnSiNi2A钢的动态本构与损伤参数 [J]. 高压物理学报, 2017, 31(3): 239–248. DOI: 10.11858/gywlxb.2017.03.005.LI L, ZHANG X F, WU X, et al. Dynamic constitutive and damage parameters of 30CrMnSiNi2A steel with different hardnesses [J]. Chinese Journal of High Pressure Physics, 2017, 31(3): 239–248. DOI: 10.11858/gywlxb.2017.03.005. [4] HANCOCK J W, MACKENZIE A C. On the mechanisms of ductile failure in high-strength steels subjected to multi-axial stress-states [J]. Journal of the Mechanics and Physics of Solids, 1976, 24(2−3): 147–169. DOI: 10.1016/0022-5096(76)90024-7. [5] MIRZA M S, BARTON D C, CHURCH P. The effect of stress triaxiality and strain-rate on the fracture characteristics of ductile metals [J]. Journal of Materials Science, 1996, 31(2): 453–461. DOI: 10.1007/BF01139164. [6] CLAUSEN A H, BORVIK T, HOPPERSTAD O S, et al. Flow and fracture characteristics of aluminium alloy AA5083-H116 as function of strain rate, temperature and triaxiality [J]. Materials Science and Engineering: A, 2004, 364(1−2): 260–272. DOI: 10.1016/j.msea.2003.08.027. [7] BAO Y B, WIERZBICKI T. On fracture locus in the equivalent strain and stress triaxiality space [J]. International Journal of Mechanical Sciences, 2004, 46(1): 81–98. DOI: 10.1016/j.ijmecsci.2004.02.006. [8] 李营, 吴卫国, 朱海清, 等. 基于大范围应力三轴度的船用低碳钢断裂特性研究 [J]. 中国造船, 2016, 57(3): 54–64. DOI: 10.3969/j.issn.1000-4882.2016.03.007.LI Y, WU W G, ZHU H Q, et al. Fracture properties of ship-build mild steel in wide range of stress triaxiality [J]. Shipbuilding of China, 2016, 57(3): 54–64. DOI: 10.3969/j.issn.1000-4882.2016.03.007. [9] 张朴, 赵明艳, 李营, 等. 复杂应力状态下船用低碳钢断裂应变的“断崖现象” [J]. 船舶力学, 2020, 24(2): 191–196. DOI: 10.3969/j.issn.1007-7294.2020.02.008.ZHANG P, ZHAO M Y, LI Y, et al. Failure strain sharp-down of marine mild steel under complex stress states [J]. Journal of Ship Mechanics, 2020, 24(2): 191–196. DOI: 10.3969/j.issn.1007-7294.2020.02.008. [10] 桂良进, 高付海, 范子杰. 先进高强度钢的断裂失效准则研究 [J]. 固体力学学报, 2012, 33(4): 395–403. DOI: 10.19636/j.cnki.cjsm42-1250/o3.2012.04.007.GUI L J, GAO F H, FAN Z J. Study of failure criterion of advanced high strength steel [J]. Chinese Journal of Solid Mechanics, 2012, 33(4): 395–403. DOI: 10.19636/j.cnki.cjsm42-1250/o3.2012.04.007. [11] MAE H, TENG X, BAI Y, et al. Calibration of ductile fracture properties of a cast aluminum alloy [J]. Materials Science and Engineering: A, 2007, 459(1/2): 156–166. DOI: 10.1016/j.msea.2007.01.047. [12] 朱浩, 朱亮, 陈剑虹. 应力三轴度和应变率对6063铝合金力学性能的影响及材料表征 [J]. 材料科学与工程学报, 2007, 25(3): 358–362. DOI: 10.3969/j.issn.1673-2812.2007.03.009.ZHU H, ZHU L, CHEN J H. Influence of stress triaxiality and strain rate on the mechanics behavior of 6063 aluminum alloy and material characterization [J]. Journal of Materials Science & Engineering, 2007, 25(3): 358–362. DOI: 10.3969/j.issn.1673-2812.2007.03.009. [13] BAO Y B, WIERZBICKI T. A comparative study on various ductile crack formation criteria [J]. Journal of Engineering Materials and Technology, 2004, 126(3): 314–324. DOI: 10.1115/1.1755244. 期刊类型引用(10)
1. 张阿漫,明付仁,刘云龙,李帅,王诗平. 水下爆炸载荷特性及其作用下的舰船毁伤与防护研究综述. 中国舰船研究. 2023(03): 139-154 . 百度学术
2. 孙裴俊,刘宜平. 水库大坝在爆炸作用下的工程防护措施研究进展. 防护工程. 2023(06): 50-57 . 百度学术
3. Jian Liu,Feng-jiang An,Cheng Wu,Sha-sha Liao,Ming-xue Zhou,Dong-yu Xue,Huan Guo. Experimental investigations on small-and full-scale ship models with polyurea coatings subjected to underwater explosion. Defence Technology. 2022(07): 1257-1268 . 必应学术
4. 马上,陈叶青,王振清,汪剑辉,胡晶. 水下爆炸离心机试验中爆炸荷载的缩比关系研究. 土木工程学报. 2021(10): 20-29 . 百度学术
5. 宋歌,龙源,钟明寿,王敏,吴建宇. 离心机和密闭加压容器中水下爆炸相似关系. 爆炸与冲击. 2019(02): 95-103 . 本站查看
6. 张阿漫,王诗平,彭玉祥,明付仁,刘云龙. 水下爆炸与舰船毁伤研究进展. 中国舰船研究. 2019(03): 1-13 . 百度学术
7. 陈材,石全,尤志锋,郭驰名,戈洪宇. 圆柱形弹药空气中爆炸相似性规律. 爆炸与冲击. 2019(09): 37-43 . 本站查看
8. HU Jing,CHEN ZuYu,ZHANG XueDong,WEI YingQi,LIANG XiangQian,LIANG JianHui,MA GuoWei,WANG QiuSheng,LONG Yuan. Underwater explosion in centrifuge part Ⅰ:Validation and calibration of scaling laws. Science China(Technological Sciences). 2017(11): 1638-1657 . 必应学术
9. 胡晶,陈祖煜,魏迎奇,张雪东,梁建辉,黄志杰. 挡水结构水下爆炸作用的离心模型试验. 水利学报. 2017(09): 1073-1081 . 百度学术
10. 金键,朱锡,侯海量,陈鹏宇,吴林杰. 水下爆炸载荷下舰船响应与毁伤研究综述. 水下无人系统学报. 2017(06): 396-409 . 百度学术
其他类型引用(12)
-