Numerical simulation on the deflection behavior of large caliber conical nose projectile at oblique high-speed water entry
-
摘要: 结合某大口径锥头弹体高速倾斜入水试验,采用任意拉格朗日-欧拉(arbitrary Lagrange-Euler,ALE)流固耦合方法对弹体倾斜入水偏转行为进行数值模拟,研究了弹体以500 m/s高速倾斜入水过程中不同受力模式、载荷变化特征以及弹体发生偏转的力学机理,分析了入水角度对弹体偏转规律的影响。结果表明:在俯仰力矩作用下,弹体均发生抬头方向偏转,且偏转速度呈现先增大后减小的趋势,偏转程度在不同入水角度范围内呈现不同的变化趋势。当入水角度小于15°时,弹体会发生“跳弹”现象;当入水角度为30°~60°时,弹体偏转趋势基本一致,均由初始倾斜状态逐渐转动至水平状态、竖直状态并最终以弹头入水反方向的“出水”姿态向水下运动;当入水角度为75°时,弹体转动至水平状态后,并未继续偏转至竖直状态,弹头以朝斜上方的姿态向水下运动;弹体的入水侵深随入水角度的增大而增大,且增大趋势近似满足指数函数关系。Abstract: Integrated with a high-velocity oblique water entry test of a large caliber conical nose projectile, the deflection behavior of the projectile obliquely entering water was studied based on the arbitrary Lagrange-Euler (ALE) fluid-structure coupling method. Firstly, based on the experiment of a projectile impacting an inclined water tank at a high velocity, a finite element model was established to simulate the corresponding response characteristics, and the rationality of the numerical model and related method was verified. Secondly, the variation of the contact force mode and load characteristics of the projectile when the water entry velocity was 500 m/s was analyzed, and the corresponding mechanism was discussed. In addition, the influence of water entry angle on the deflection behavior of the projectile was investigated. The related analysis demonstrates that the projectile will deflect upward due to the effect of pitch moment, and the deflection velocity increases first and then decreases gradually during the entry process. The variation trends of deflection degrees are different within different entry angle ranges. When the entry angle is less than 15º, the projectile jumps usually out of the water. When the entry angle is in the range from 30º to 60º, the deflection trend of the projectile is almost the same, i.e., the projectile rotates from the initial tilted state to a horizontal state, then to a vertical state, and moves finally downwards with its nose in the opposite direction to the initial water entry direction. When the entry angle increases to 75º, the projectile cannot continue to rotate to a vertical state after it rotates to a horizontal state, but instead moves downwards in a tilted state with its nose facing upwards. Under different water entry angles, the axial force exerted on the projectile is negative, leading to a continuous decrease in the projectile velocity. Comparatively, the transverse force is positive, and the peak value decreases with increasing water entry angle. Moreover, the penetration depth of the projectile increases with the increase of entry angle, and it almost shows an exponential relationship.
-
地下爆炸广泛应用于工业、农业和国防等领域。因此,有关地下介质的爆破破碎机理、爆炸波数值模拟、地震波传播特征和地震效应等,已经进行了大量研究[1-10]。尤其是地下爆炸发生的地震效应, 虽然爆炸地震和天然地震的形成机理存在很大差别,但是他们都激发地震波,引起地震动,对地表和人类工程建筑物造成破坏,在某种程度上,破坏效果是类似的。而造成地震破坏效应的强弱主要取决于爆炸激发的地震波引起的地介质质点运动参数大小。因此,有关地下爆炸引起的地震波传播规律、介质质点运动参数计算方法,一直备受关注。
本文中,基于第四纪砂砾层中小当量(即当量在公斤级至百公斤级范围)的封闭爆炸的观测资料,针对波形中优势横波进行分析处理,分析爆炸地震波水平与垂向质点速度随当量和距离变化,并对几种模型在描述观测数据方面的残差进行分析对比,确定描述第四纪砂砾层介质的地震波质点运动速度模型。
1. 数学模型
介质的质点运动参数与爆炸当量、方式、源区介质和观测距离等众多因素有关。一般地, 爆炸当量越大, 激发的地震波就越强, 爆炸封闭越好, 爆炸能量耦合到介质中能量越高。爆炸能量耦合到介质中的比例(即能量耦合系数)主要和爆炸方式及爆炸所处的介质等因素有关。学者建立数学模型时,对地下封闭爆炸和抛洒爆炸、岩石与土层分别研究,从而大大简化了模型,使数学模型更有实用性[5]。归结起来,现有模型主要有以下几种。
(1) Sadauskas模型。Sadauskas模型[5]的地下爆炸质点速度为:
v=K(3√WR)α (1) 式中:质点速度v的单位为cm/s,当量W的单位为kg,距离R的单位为m。对于岩石,K=50,对于土壤,K=200,α=1~2。
(2) 显函模型。谢毓寿等[10]基于炸药量为0.05~2.0 kg TNT、距离在35~404 m范围内的数据,确定了振动质点速度为:
lgv=K0+αlgW+βlgR 上式可改写为:
v=KWαRβ (2) 式中:K=10K0,对于坚硬岩石,K0=2.2~2.6,α=0.60,β=1.80。
双极模型。用Sadauskas模型参数计算时,先计算比例距离,并且规定恒定的指数。而用显函模型计算α、β时,要基于较多样本参量,且需要各次爆炸空间位置不变等条件,否则将带来一定数据误差。为解决这个问题,双极模型[6]被提了出来。该模型在计算模型参数时不同于传统的上述两个模型计算方法,而是在一定当量指数变化区间内,查找最合适参数,使模型预测数据和实测数据具有更高的一致性。质点速度为:
v=K(WαR)β (3) 式中:对于黄土层,K为0.40,α为2.0,β为2.3。
2. 实验与检测
2.1 场地介质
野外爆炸实验场地相当开阔, 地形没有明显起伏,地面十分平坦。整个场地未见坚硬基岩,出露的地层为第四纪冲积层,为砂砾石混合体,砾石大小不一,多数砾石呈现长椭球和不规则形状,具有较低的磨圆度和筛选度。在地表以下至6 m之间,砂砾层中含砂比例较高, 6 m以下直至20 m深处, 砂砾层中含砂比例明显降低。实验场地构造相对简单,未发育大的断裂构造。砂砾石混合体结构十分松散,地表仅被一些植物稀疏覆盖,宏观上为一片植被不甚发育的荒漠区域。
2.2 检测仪器与布设
所用的检测地震波仪器是短周期型LE-3Dlite,机电换能灵敏度为400V/(m·s-1), 检波器噪声很低,能够检测1~80 Hz范围内的振动信号。该地震检波仪器为三分向组合体,可以在三维方向检测地震波质点振动幅值。
地表出露的为第四纪砂砾层,地势比较平坦。在爆炸点周围不同方向,距离在0.5~5 km范围内,布设了十几个地震波检测点,爆炸源与检测点分布如图 1所示。
系列爆炸源主要为20、100和300 kg的TNT炸药,从地表向下开挖约20 m深井,将炸药放在井底部,然后用挖出来的砂砾石进行回填,增强爆炸能量耦合,提高弹性能在爆炸能量中的比例。
3. 测量结果
3.1 测量概况
为确定本底噪声,在E1台站记录信号前截取了20s的数据。其本底波形比较稳定,振幅大约为0.5 μm/s,傅立叶振幅从低频到高频渐渐降低,最大傅立叶振幅为零频极限值,优势频带主要限制在3 Hz以下。针对这次实验布设十几个观测点,均记录到了良好的数据。地震波形持时很短,仅为数秒钟,两个水平方向的地震波形,在振动幅值大小、包络线形态和波形持时方面具有较好的一致性,而垂向和水平向波形在上述几方面明显不同,波形振幅远强于水平向振幅。随着距离增大,波形复杂性明显增加,如图 2所示。
3.2 质点运动速度
实验所用的LE-3Dlite地震计,能够测量东西、南北和垂向3个方向的地运动数据,通过灵敏度换算,可以直接得出地运动质点速度ve、vn、vv。表 1为爆炸当量为300 kg爆炸实验部分观测点质点速度。根据表 1,质点水平速度ve与南北速度vn互有大小,ven为水平均值,均大于垂直速度。统计上ve大于vn,垂向vv小于任意水平上幅值。数据表明,爆炸激发的地震波,水平上振幅强于垂向振幅。
表 1 300 kg实验部分观测点质点速度Table 1. Partial observation data of 300 kg explosion台站 R/km ve/(μm·s-1) vn/(μm·s-1) vv/(μm·s-1) ven/(μm·s-1) W2 1.038 0 288.00 619.29 179.55 453.65 W3 1.599 3 126.17 119.70 99.95 122.94 W4 2.144 7 170.35 92.06 51.20 131.20 W5 2.980 6 49.70 33.14 31.43 41.42 W6 4.967 7 13.10 8.73 4.77 10.92 E1 0.521 6 1 172.70 1 227.42 361.80 1 200.10 E3 1.151 8 221.67 122.48 96.97 172.08 E4 1.844 2 61.54 141.97 48.56 101.75 4. 模型确定
4.1 Sadauskas模型
根据式(1),基于一系列小当量地下爆炸不同观测点的距离和相应地表质点速度,先计算各观测点的数据 3√W/R ,然后,采用最小二乘法,进行回归计算。本文中对水平向质点速度vh和垂直向质点速度vv进行处理和回归,计算结果显示:水平向的线性相关因子0.960,残差0.245,模型中α为2.37,K为3.81;垂向的线性相关因子0.938,残差0.297,模型中α为2.26,K为2.28。拟合曲线如图 3所示,于是Sadauskas模型为:
{vh=3.81(3√WR)2.37vv=2.28(3√WR)2.26 (4) 根据上述回归参数,水平向质点速度相对于垂向质点速度,随3√W/R增长比垂直向快,二者的线性相关性更强。
4.2 显函模型
式(2)是距离R、W的二元函数,该式反映了地介质质点速度随距离增加而衰减、随当量增长而增长的关系,随距离的衰减指数可以通过一次实验不同距离的检测点速度计算确定,而随当量的指数可以通过同一观测点不同当量的实验数据得出。由于地下介质的复杂性,存在介质能量吸收和球面扩散作用,地震波在传播过程中,能量逐渐消耗,地震波幅值随着距离增加发生衰减,不同介质中衰减指数是不同的。本文中采用最小二乘法,对地震波的水平和垂直数据进行了回归计算,其计算结果如图 4所示。
图 4(a)为质点水平速度与距离的拟合曲线。图中3种不同当量爆炸的拟合曲线大致是平行的,距离指数为-2.07,相关系数为0.97。图 4(c)为质点垂直速度与距离拟合曲线, 距离指数为-1.57,相关系数为0.99。二者相比,质点垂直速度相比于水平速度,与距离的相关性更强。图 4(b)为质点水平速度与当量拟和曲线。图为20~300的4次爆炸,在0.51、1.03和4.99 km的距离测点上的数据拟和曲线,这几条曲线大致是平行,但存在一定差别,当量指数约为1.09,相关系数为0.99。图 4(d)为质点垂直速度与当量拟合曲线, 当量指数为0.77,相关系数为0.97。于是显函模型为:
{vh=1.79W1.09R−2.07vv=2.82W0.77R−1.57 (5) 4.3 双极模型
由于爆炸激发的质点振动参数和爆炸当量正相关,随着当量增大而增大,且呈显出指数增长关系。双极模型参数测定完全不同于传统的Sadauskas模型和显函模型模型测定方法,采用区间寻优方法,确定最佳的当量指数[4]。具体计算时,首先,给定一个当量指数,计算与距离的比值,利用最小二乘法进行回归计算,得出模型参数,采用模型反演质点速度数据,并和实测数据进行对比,确定实测和反演数据之间的残差值。然后, 用一定的步长,选择下一个当量指数,计算残差,步长越小,计算的结果越细致,但计算量也越大。将选定的区间内所有点计算完毕,从所有计算残差列表中选取最小残差所对应的参数,即为速度模型参数。本文中在1.0至2.5之间采用0.05步长进行计算,部分计算结果见表 2,如图 5所示。α、β、b、γ和σ分别表示当量开方次数、比例距离指数、回归截距、相关系数和残差。
表 2 双极模型部分计算数据Table 2. Partial calculating data from double extreme modelα βh bh γh σh βv bv γv σv 1.70 1.891 4 0.141 0 0.976 4 0.188 7 1.760 8 -0.289 8 0.930 0 0.314 2 1.80 1.957 6 0.078 7 0.978 7 0.179 6 1.827 1 -0.235 7 0.934 5 0.304 1 1.90 2.017 4 -0.016 3 0.980 0 0.174 1 1.887 5 -0.181 2 0.938 1 0.296 1 2.00 2.071 2 0.045 6 0.980 5 0.172 0 1.942 4 -0.126 6 0.940 7 0.289 9 2.10 2.119 4 0.106 5 0.980 2 0.173 0 1.992 1 -0.072 6 0.942 6 0.285 4 2.20 2.162 6 0.166 2 0.979 4 0.176 5 2.036 9 -0.019 3 0.943 8 0.282 5 2.30 2.201 1 0.224 3 0.978 1 0.182 1 2.077 3 0.033 0 0.944 3 0.281 1 2.40 2.235 3 0.280 8 0.976 3 0.189 2 2.113 5 0.084 0 0.944 4 0.280 9 2.50 2.265 6 0.335 6 0.974 2 0.197 4 2.146 1 0.133 6 0.944 1 0.281 8 2.60 2.292 4 0.388 4 0.971 7 0.206 4 2.175 2 0.181 8 0.943 3 0.283 6 根据表 2,水平向在当量指数1.7至2.0区间内,反演与实测数据之间的残差,从开始的逐渐减小,到2.0时达到极小值,为0.172,从2.0到2.6之间又逐渐回升增大。整个曲线为凹向上形态,如图 5(a)所示。而相关系数在在2.0达到极大值,相关性最强。垂直向在当量指数2.4时,同时达到残差最小和相关系数最大。残差变化如图 5(b)所示。图 5(c)~(d)为水平与垂直向质点速度拟合曲线。双极模型为:
{vh=1.11(W1/2.0R)2.07vv=1.21(W1/2.40R)2.11 (6) 5. 模型分析
Sadauskas模型中,质点水平向与垂向速度随3√W/R变化的指数分别为2.37和2.28,相关系数分别为0.960和0.938。数据表明,质点速度和3√W/R相关性是比较强的。
显函模型中,采用固定系列爆炸位置,确定质点速度随当量增大指数。事实上,爆炸位置尽管变化不大,但每次的地下爆炸不能完全重复,各爆炸位置之间依然有数十米的距离。严格来讲,位置不变,工程上是难以实现。因此当量指数测定存在一定误差。水平质点速度的当量指数具有一定分散性,均值为1.09、20~300 kg的质点速度随距离的衰减指数,一致性比较好,均值为-2.07。垂直向质点速度的当量指数为0.77、20~300 kg的质点速度随距离的衰减指数在-1.55至-1.58之间,均值为1.57。距离衰减指数相对于当量增长指数,精度要高得多。
双极模型中,水平与垂直向质点速度当量的指数分别是1/2.0和1/2.4,其残差分别0.172和0.280。而Sadauskas模型的水平与垂直向速度的残差分别为0.245和0.297,显函模型的水平与垂直向速度的残差分别为0.171和0.300。双极模型与Sadauskas模型和显函模型相比,在反演或预测质点速度参数方面,更加接近实际数据。
6. 结论
在第四纪砂砾介质中进行了公斤至百公斤级系列地下封闭爆炸,通过对0.5~5.0 km范围地震波测量和数据分析,得出如下结论。
地震波水平振幅强于垂向振幅, 质点速度随当量呈现指数增加,水平指数为1.09,垂直指数为0.77,质点速度随距离呈现指数衰减,水平指数为2.07,垂直指数为1.57。Sadauskas模型、显函模型和双极模型都能够描述小当量地下爆炸地震波质点速度。反演计算结果显示,双极模型的残差最小,显函模型次之,Sadauskas模型最大。换言之,采用双极模型反演的数据更加接近实际。
-
材料 ρ/(kg·m–3) E/GPa ν cp/(J·kg–1·K–1) Tr/K Tm/K ˙ε/s–1 A/MPa B/MPa n C G50 7 620 205 0.28 469.0 300 1765 1 1 445 1 326 0.356 0.005 7A04 2 850 69.35 0.33 921.0 293 878 1 602.5 732.1 0.753 0.014 材料 m D1 D2 D3 D4 D5 c0/(m·s–1) S1 γ0 a0 G50 1.120 0.100 0.760 1.57 0 0 4280 1.99 2.00 0.46 7A04 1.015 0.059 0.246 –2.41 –0.1 –0.1 5240 1.40 1.97 0.48 -
[1] VON KARMAN T H. The impact on seaplane floats during landing: NACA technical note No. 321 [R]. Washington: NACA, 1929. [2] LOGVINOVICH G V. Hydrodynamics of flows with free boundaries [M]. Kiev: Naukova Dumka, 1969. [3] MAY A, WOODHULL J C. Drag coefficients of steel spheres entering water vertically [J]. Journal of Applied Physics, 1948, 19(12): 1109–1121. DOI: 10.1063/1.1715027. [4] 陈先富. 弹丸入水空穴的试验研究 [J]. 爆炸与冲击, 1985, 5(4): 70–73.CHEN X F. Experimental studies on the cavitation phenomena as a pellet entering water [J]. Explosion and Shock Waves, 1985, 5(4): 70–73. [5] 张伟, 郭子涛, 肖新科, 等. 弹体高速入水特性实验研究 [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. [6] 郭子涛, 张伟, 郭钊, 等. 截卵形弹水平入水的速度衰减及空泡扩展特性 [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. [7] 刘思华, 王占莹, 李利剑, 等. 头型对射弹高速入水稳定性的影响 [J]. 航空学报, 2023, 44(21): 528437. DOI: 10.7527/S1000-6893.2023.28437.LIU S H, WANG Z Y, LI L J, et al. Influence of nose shapes on high-speed water entry stability of projectile [J]. Acta Aeronautica et Astronautica Sinica, 2023, 44(21): 528437. DOI: 10.7527/S1000-6893.2023.28437. [8] 王云, 袁绪龙, 吕策. 弹体高速入水弯曲弹道实验研究 [J]. 兵工学报, 2014, 35(12): 1998–2002. DOI: 10.3969/j.issn.1000-1093.2014.12.010.WANG Y, YUAN X L, LV C. Experimental research on curved trajectory of high-speed water-entry missile [J]. Acta Armamentarii, 2014, 35(12): 1998–2002. DOI: 10.3969/j.issn.1000-1093.2014.12.010. [9] SHI Y, HUA Y, PAN G. Experimental study on the trajectory of projectile water entry with asymmetric nose shape [J]. Physics of Fluids, 2020, 32(12): 122119. DOI: 10.1063/5.0033906. [10] 马庆鹏, 魏英杰, 王聪, 等. 不同头型运动体高速入水空泡数值模拟 [J]. 哈尔滨工业大学学报, 2014, 46(11): 24–29. DOI: 10.11918/j.issn.0367-6234.2014.11.004.MA Q P, WEI Y J, WANG C, et al. Numerical simulation of high-speed water entry cavity of cylinders [J]. Journal of Harbin Institute of Technology, 2014, 46(11): 24–29. DOI: 10.11918/j.issn.0367-6234.2014.11.004. [11] GAO J G, CHEN Z H, HUANG Z G, et al. Numerical investigations on the oblique water entry of high-speed projectiles [J]. Applied Mathematics and Computation, 2019, 362: 124547. DOI: 10.1016/j.amc.2019.06.061. [12] CHANG Y N, TONG A Y. A numerical study on water entry of cylindrical projectiles [J]. Physics of Fluids, 2021, 33(9): 093304. DOI: 10.1063/5.0059892. [13] 肖海燕, 罗松, 朱珠, 等. 高速射弹小角度入水弹道特性研究 [J]. 北京理工大学学报, 2019, 39(8): 784–791. DOI: 10.15918/j.tbit1001-0645.2019.08.003.XIAO H Y, LUO S, ZHU Z, et al. Trajectory and cavitation characteristics of high-speed projectiles at small angle of water entry [J]. Transactions of Beijing Institute of Technology, 2019, 39(8): 784–791. DOI: 10.15918/j.tbit1001-0645.2019.08.003. [14] 黄振贵, 王瑞琦, 陈志华, 等. 90°锥头弹丸不同速度下垂直入水冲击引起的空泡特性 [J]. 爆炸与冲击, 2018, 38(6): 1189–1199. DOI: 10.11883/bzycj-2018-0115.HUANG Z G, WANG R Q, CHEN Z H, et al. Experimental study of cavity characteristic induced by vertical water entry impact of a projectile with a 90° cone-shaped head at different velocities [J]. Explosion and Shock Waves, 2018, 38(6): 1189–1199. DOI: 10.11883/bzycj-2018-0115. [15] 胡明勇, 张志宏, 刘巨斌, 等. 低亚声速射弹垂直入水的流体与固体耦合数值计算研究 [J]. 兵工学报, 2018, 39(3): 560–568. DOI: 10.3969/j.issn.1000-1093.2018.03.018.HU M Y, ZHANG Z H, LIU J B, et al. Fluid-solid coupling numerical simulation on vertical water entry of projectile at low subsonic speed [J]. Acta Armamentarii, 2018, 39(3): 560–568. DOI: 10.3969/j.issn.1000-1093.2018.03.018. [16] GUO Z T, ZHANG W, XIAO X K, et al. An investigation into horizontal water entry behaviors of projectiles with different nose shapes [J]. International Journal of Impact Engineering, 2012, 49: 43–60. DOI: 10.1016/j.ijimpeng.2012.04.004. [17] SONG Z J, DUAN W Y, XU G D, et al. Experimental and numerical study of the water entry of projectiles at high oblique entry speed [J]. Ocean Engineering, 2020, 211: 107574. DOI: 10.1016/j.oceaneng.2020.107574. [18] 李佳川, 魏英杰, 王聪, 等. 不同扰动角速度高速射弹入水弹道特性 [J]. 哈尔滨工业大学学报, 2017, 49(4): 131–136. DOI: 10.11918/j.issn.0367-6234.201512058.LI J C, WEI Y J, WANG C, et al. Water entry trajectory characteristics of high-speed projectiles with various turbulent angular velocity [J]. Journal of Harbin Institute of Technology, 2017, 49(4): 131–136. DOI: 10.11918/j.issn.0367-6234.201512058. [19] 汪振, 吴茂林, 戴文留. 大口径弹体高速入水载荷特性研究 [J]. 弹道学报, 2020, 32(1): 15–22. DOI: 10.12115/j.issn.1004-499X(2020)01-003.WANG Z, WU M L, DAI W L. Study on load characteristics of high-speed water-entry of large caliber projectile [J]. Journal of Ballistics, 2020, 32(1): 15–22. DOI: 10.12115/j.issn.1004-499X(2020)01-003. [20] 孙玉松, 周穗华, 张晓兵, 等. 基于多介质ALE方法的大型弹体入水载荷特性研究 [J]. 海军工程大学学报, 2019, 31(6): 101–106. DOI: 10.7495/j.issn.1009-3486.2019.06.019.SUN Y S, ZHOU S H, ZHANG X B, et al. On water-impact load of heavy projectiles base on multi-material ALE method [J]. Journal of Naval University of Engineering, 2019, 31(6): 101–106. DOI: 10.7495/j.issn.1009-3486.2019.06.019. [21] 张斌, 李继承, 陈建良, 等. 构型弹体跌落冲击载荷及结构响应特性 [J]. 爆炸与冲击, 2023, 43(3): 033201. DOI: 10.11883/bzycj-2022-0098.ZHANG B, LI J C, CHEN J L, et al. Loading characteristics and structural response of a warhead during drop impact [J]. Explosion and Shock Waves, 2023, 43(3): 033201. DOI: 10.11883/bzycj-2022-0098. [22] 张伟, 肖新科, 魏刚. 7A04铝合金的本构关系和失效模型 [J]. 爆炸与冲击, 2011, 31(1): 81–87. DOI: 10.11883/1001-1455(2011)01-0081-07.ZHANG W, XIAO X K, WEI G. Constitutive relation and fracture model of 7A04 aluminum alloy [J]. Explosion and Shock Waves, 2011, 31(1): 81–87. DOI: 10.11883/1001-1455(2011)01-0081-07. 期刊类型引用(1)
1. 韩向东,王西泉,李超,陈俊彪,李军智,封慧勇,魏雪梅. 多因素下弹体跨介质入水试验系统构建. 兵工学报. 2024(S2): 199-207 . 百度学术
其他类型引用(3)
-