Asay window for probing the microspall of materials
-
摘要: 微层裂是冲击波物理领域的重要基础问题,在工程上具有重要应用价值。近年来用于诊断样品多层层裂的传统Asay窗技术被用于诊断微层裂,但对其诊断能力和信号特征认识仍存在严重不足。为此,通过波系分析,揭示出在薄飞片击靶的微层裂实验中样品破碎存在1个“痂片”特征区、2个微层裂特征区以及1个“残体”特征区。实验表明,在样品窗口间隙合适的条件下,Asay窗不仅能够有效区分这些不同特征分区,而且能够灵敏探测样品表面发射的高速微喷粒子,从而实现对样品连续破碎过程的精密诊断。Abstract: Microspall is an essential problem in both theoretical investigation and engineering application in shock physics. The Asay window, originally developed to diagnose the multi-spall behavior of material, was recently employed to probe the microspall, but its ability for probing the problem calls for further demonstration, and the corresponding signal also needs further explanation. In this paper, wave propagation analysis indicates that the sample bearing microspall can be separated into several different characteristic regions, the experiments performed demonstrate that the Asay window can sensitively distinguish these regions as far as a reasonable experimental configuration is set up, and even the features of the micro jet particles can be detected. So the technique was proved of great value for dynamic fragmentation studies.
-
Key words:
- Asay window /
- wave propagation analysis /
- microspall /
- shock wave
-
脉冲爆轰发动机(pulse detonation engine, 简称PDE)是一种利用脉冲式爆轰波产生推力的新概念发动机,具有热效率高、结构简单等优点。多管PDE的构想最早出现在20世纪50年代,相对于单爆轰管的PDE,增加爆轰管数量具有增加推力和提高工作频率等优点。国内外对多爆轰管PDE研究颇有成效,2014年T.Morozumi等[1]采用乙烯为燃料,N2O为氧化剂,实现了四爆轰管PDE装置的成功起飞,工作过程中对称的两爆轰管同时起爆,在1500 ms的工作过程中测得平均推力为256 N,比冲为130.8 s。2011年,袁成等[2]以汽油为燃料,空气为氧化剂,成功实验了六个爆轰管的吸气式PDE的多循环起爆。然而作为一种动力装置,其在工作过程中产生较大的噪声会引起结构的声疲劳、对相关工作人员的听力造成损伤、对周围环境产生噪声污染等危害。因此,对多爆轰管PDE爆轰噪声的研究具有重要意义。
针对单爆轰管PDE爆轰噪声特性已有一系列的研究报道,L.Shaw等[3]对距离出口0~3.66 m范围内的PDE管外噪声进行了测量,发现距离管口越近,声压级越高。通过信号频率的处理,发现随着距管口的距离的增加,由于高频噪声在大气中的耗散作用,其声压的幅值迅速衰减。美国辛辛那提大学A.Glaser等[4]对距离出口方向不同角度、不同轴向距离的PDE爆轰噪声进行了实验,发现爆轰噪声辐射具有一定的指向性,并且得出了噪声幅值随轴向距离变化的衰减规律。国内对于PDE爆轰噪声的实验研究相对较少。郑龙席等[5]进行了气液两相PDE噪声实验研究,发现其噪声辐射呈现间歇性和周期性的特点,噪声源主要由单极子源和四极子源组成;噪声辐射的频谱为宽频带谱,其能量主要集中在低频部分,由气流脉动的基频和谐频组成;峰值声压级与爆轰频率无关,而脉冲声压级、声功率级和声功率随爆轰频率的增加而提高。许桂阳等[6]对加装不同喷管的PDE爆轰噪声进行了实验研究,发现加装喷管可以降低爆轰噪声和减小爆轰噪声的参考半径。
对于多管PDE爆轰噪声的研究,美国辛辛那提大学N.Caldwell等[7]对双爆轰管的PDE爆轰噪声进行了相关的实验研究,与单爆轰管PDE相比,双爆轰管PDE在距离管口2.89 m处产生的声压级要高出至少3 dB,并且爆轰噪声辐射的指向性更加明显。当前,国内多管PDE爆震噪声的研究尚未见到报道。
在前人研究的基础上,本文中着重于对正三角形这一典型多爆轰管组合方式的PDE爆轰噪声展开研究,通过实验探索多循环三爆轰管PDE近场噪声特性,分析其波形特征与衰减规律,得到爆轰噪声峰值衰减系数,并对噪声辐射指向性和持续时间进行研究。研究结果有助于探索多爆轰管PDE爆轰的噪声形成和传播机理以及多爆轰管PDE爆轰噪声的控制方法。
1. 实验系统
三爆轰管PDE实验系统如图 1所示,主要包括PDE、供气系统、供油系统、点火控制系统和测试系统等。
以压缩空气和压缩氧气作为氧化剂,燃料为93#汽油。3个爆轰管直径均为80 mm,摆放成正三角形样式,如图 1所示,中心轴间距为200 mm。PDE管头部以径向方式进气,通过精细雾化喷嘴和文氏管实现汽油雾化,并在混合室内实现氧化剂和燃料充分均匀混合。采用高能点火头,通过时序控制,实现三爆轰管同时起爆,点火频率为8 Hz,实验在常温常压下进行。在每个爆轰工作循环间隙内,会有新鲜燃料和氧化剂的混合物吹扫废气,同时达到完全填充PDE管,此时下一次点火开始,进行下一次工作循环。测试系统中的传感器为PCB公司高频动态传感器,传感器布置在垂直于爆轰管中心轴连线的水平面内0°、30°、45°和90°方向上,高度与图 1中正三角形中心点保持一致。在同一工况下,分别对距PDE出口400、600、800、1 000、1 200、1 600、2 400和3 000 mm处的噪声进行测试。测试得到的噪声电信号通过信号调理器处理后由同步数据采集系统采集,采样率为5×105 s-1。
2. PDE近场爆轰噪声传播特性分析
2.1 三管爆轰发动机爆轰噪声参考半径的变化
以r0为参考半径,将PDE噪声场划分为近场和远场,计算公式为[4]:
r0=(Ep0)13 (1) 式中:r0为参考半径, E为爆炸能量, p0为大气压力。定义r为距离三爆轰管名义“中点”的直线距离。当r<r0时,PDE爆轰噪声峰值的衰减与r的倒数满足3次方变化规律,而在参考半径r0之外时,PDE爆轰噪声峰值的衰减与r的倒数满足1次方变化规律:
p peak ={k(r0r)3r⩽r0k1r0rr>r0 (2) 式中:k、k1为衰减系数。
因此,参考半径r0对PDE爆轰噪声机理研究具有重要意义。图 2为不同方向上PDE爆轰噪声幅值随距离r的变化曲线。其中实线为负3次方拟合曲线,虚线为负1次方拟合曲线。
PDE近场爆轰噪声峰值ppeak与r的倒数满足一定的变化规律[4](公式(2))。
由曲线拟合情况可以看出:在距离管口较近的区域,三爆轰管PDE噪声场中的幅值衰减较快,衰减速度与1/r的3次方成正比。衰减速度逐渐放缓,在距离管口较远的地方,衰减速度与1/r成正比。这主要是因为在管口附近,管内爆轰波退化形成的高强度的激波和高温高压燃烧产物迅速喷出管口形成高压力区。伴随着高强度冲击波的迅速衰减和高温高压区的迅速扩散,压力幅值迅速下降。随着冲击波继续向前传播,冲击波强度逐渐减弱,其衰减速度也减慢。而中间距离的测试点,既不符合3次方,也不符合1次方的衰减规律。可将此区域定义为过渡区[8],即由强非线性声场向弱非线性声场衰减的过渡区域。该区域内,冲击波强度迅速衰减,并且衰减速度逐渐放缓,最终由衰减逐步蜕变为弱非线性声波,其声压的衰减速度较快且规律性较差。根据公式(2),以两条曲线的交汇点定义出参考半径r0,由0°、30°、60°和90°等4个方向的拟合曲线可以看出,参考半径r0≈0.7 m。本文中实验所用的爆轰管的直径d=80 mm,因此参考半径r0≈9倍管径。求得在0°、30°、60°和90°下,3次方衰减系数k=5.13×10-2、9.13×10-2、8.42×10-2和4.35×10-2; 1次方的衰减系数k1=5.16×10-2、8.26×10-2、6.41×10-2和4.53×10-2。
2.2 PDE近场爆轰噪声指向性分析
作为脉冲噪声,PDE近场爆轰噪声主要评价有脉冲声压级和峰值声压级。本文中着重研究峰值声压级,即爆轰噪声声压-时间曲线中声压峰值ppeak的分贝值,表达式为:
L peak =20lg(p peak /p0) (3) 式中:Lpeak为峰值声压级, dB;ppeak为爆轰噪声的峰值声压, Pa;p0=2×10-5 Pa, 为基准声压。
指向性是表示声源在不同的方向辐射声能量的差异性。由图 3可知,通过实验测得三管PDE近场爆轰噪声具有较强的指向性。三管PDE近场爆轰噪声峰值最大值均出现在30°方向。并且随着距离的增加,指向性更加明显。爆轰噪声,近似于高速喷流噪声的一种,主要由冲击波噪声和高速喷流噪声组成,冲击波噪声由爆轰波衰减而成的冲击波组成,高速喷流噪声由高温高压爆轰产物喷出PDE形成的高速喷流产生。文献[5]中对爆轰噪声的频谱计算分析认为,爆轰噪声的声源主要由单极子和四极子源组成。其中的四极子源,即由射流边界层中产生的湍流脉动引起的高频噪声,具有明显的指向性。文献[4]中,测量了单管PDE爆轰噪声指向性,得出了除1倍管径处,其余测量位置处的爆轰噪声辐射最大声压级均在30°方向上的结论,三爆轰管PDE虽不同于单爆轰管PDE,但其声源特性与单爆轰管PDE表现出了相似的规律性。
2.3 PDE噪声近场的时间特性
脉冲噪声的持续时间,也是脉冲噪声重要的评价标准。爆轰波出管口退化而成的冲击波在近场区域仍为较强的冲击波,具有很强的间断性。幅值上升快,下降也快,持续时间较短。本文中采用的判定标准有2种,即A持续时间和B持续时间[9]。按照国家标准[GJB 2A-1996]的定义,A持续时间为噪声在自由场的简单脉冲波,持续时间即指压力上升至主要正峰值,随后又迅速下降至环境压力所需要的时间。对于B持续时间,即在混合场或其他原因产生反射形成的一系列振荡波时,其持续时间为:在声压脉冲波峰值下降90%的声压振幅范围内的时间总和。图 4为0°方向,r=1 000 mm处的声压变化曲线,其中声压峰值为62 kPa, La和Lb分别为波动范围线的上下界限,大小分别为压力峰值ppeak的±1/10。图中O、A两点间的时间差为A持续时间,O、B两点间的时间差为B持续时间。
表 1为0°、30°、60°、90°等4个方向,测试距离分别为600、800、1 000、1 200 mm处的脉冲噪声的持续时间。
表 1 爆轰噪声的持续时间Table 1. Time of duration for detonation noiseφ/(°) r=600 mm r=800 mm r=1 000 mm r=1 200 mm A B A B A B A B 0 0.55 9.64 0.61 10.93 0.63 12.10 0.64 25.56 30 0.45 6.61 0.49 8.14 0.50 9.51 0.57 12.90 60 0.39 7.13 0.41 8.28 0.44 9.16 0.47 14.82 90 0.37 9.16 0.39 12.50 0.40 15.26 0.42 21.91 2.3.1 A持续时间
A持续时间随不同位置、不同方位角的变化情况如图 5(a)所示。在相同的位置r处,随着方位角φ的增加,A持续时间近似线性规律下降。在同一方位角情况下,随着距离r的增大,A持续时间也增大。分析认为,在相同的方位角上,随着轴向距离增大,一定量的脉冲声通过单位面积的声能量不断减小,因而峰压减小,压力-时间曲线拉宽,正脉冲时间也变长。在同一距离r处,由于三爆轰管PDE爆轰噪声的指向性,在30°方向能量辐射比较大,因而不仅压力值高,且A持续时间也长。而0°方向上能量分布较小,但是A持续时间较长,主要是因为三爆轰管虽然是同时点火,但由于三个爆轰管之间交汇流场的复杂性和波系之间的抵消作用,导致到达测试点的时间不一致。图 6为800 mm距离处测量所得压力随时间变化曲线,图中明显可以看出在压力上升区,存在3个峰值的叠加,导致其压力波形的正持续时间变长。
为验证对A持续时间规律分析的正确性,对单爆轰管PDE爆轰噪声近场噪声进行了实验测量。单爆轰管PDE爆轰噪声的最大值在30°方向,A持续时间随着方位角的变化曲线如图 5(b)所示:由图中可以看出,同一方位角下,随着距离的增加,A持续时间不断增加。同一轴向距离处,由于爆轰噪声辐射的指向性,在30°处,A持续时间较长,与上述分析结论一致。
2.3.2 B持续时间
B持续时间随着方位角φ的变化曲线如图 7所示。在600、800、1 000、1 200 mm处,30°和60°的B持续时间均小于同一距离下的0°和90°方向。图 8为800、1 000 mm处的爆轰噪声峰值指向性和B持续时间变化曲线。0°~30°,随着方位角的增加,峰值声压级增加,B持续时间相应减小,在60°~90°区间内,随着方位角的增加,峰值声压级下降,B持续时间增加。这两段区间内,峰值声压级的变化与B持续时间的关系成反比。在30°~60°处,B持续时间与方位角和峰值声压级没有明确的关系。在方位角相同的情况下,随着距离r的增大,B持续时间增加。
由上述分析可以看出,B持续时间与噪声峰值具有一定的关系,随着传播距离的增加,峰值逐渐减小,B持续时间增加。图 9为800和1 000 mm处0°方向上压力变化曲线,峰值为84、60 kPa。可以看出,峰值的衰减速度较快,但在小压力扰动区域,衰减速度较慢。3个流场之间的交汇作用较复杂,产生的波系也相对复杂。声场也并不是单管之间的简单叠加,在某些区域存在一定程度的抵消,某些区域存在一定程度的叠加[10],因此会出现峰值压力减小,B持续时间增大的现象。
3. 结论
通过对三爆轰管脉冲爆轰发动机同时起爆工况下爆轰噪声的研究,得出以下结论:
(1) 爆轰噪声随着传播距离的增加,迅速衰减,在3倍“名义管径”距离之前,按照距离r的负3次方的规律衰减,近场与远场的参考半径在3倍“名义管径”处,3倍“名义管径”之后一定距离处,按照距离r的负一次方的规律衰减。
(2) 爆轰噪声辐射具有明显的指向性,峰值均出现在30°方向。
(3) 在相同的距离r处,随着方位角φ的增大,A持续时间的下降规律近似线性的。在同一方位角情况下,不同距离r处,A持续时间随着r的增大而增加。
(4) B持续时间与声压峰值的大小成反比关系,并且随着距离r的增大,B持续时间增加。在不同的方位角上,峰值声压越大,B持续时间越小。
-
表 1 实验参数及结果
Table 1. Experimental parameters and results
实验编号 df/mm ds/mm d/mm v/(km·s-1) Shot 1 0.508 2.522 0.3 3.244 Shot 2 0.505 2.524 0.5 3.200 Shot 3 0.503 2.522 1.0 3.254 Shot 4 0.507 2.525 3.0 3.255 Shot 5 0.500 2.521 5.0 3.252 -
[1] Curran D R, Seaman L, Shockey D A. Dynamic failure of solids[J]. Physics Reports, 1987, 147(5):253-388. doi: 10.4028-www.scientific.net-KEM.261-263.239/ [2] 白以龙, 柯孚久, 夏蒙棼.固体中微裂纹系统统计演化的基本描述[J].力学学报, 1991, 23(3):290-298. http://www.cnki.com.cn/Article/CJFDTOTAL-LXXB199103004.htmBai Yilong, Ke Fujiu, Xia Mengfen. Formulation of statistical evolution of microcracks in solids[J]. Acta Mechanica Sinica, 1991, 23(3):290-298. http://www.cnki.com.cn/Article/CJFDTOTAL-LXXB199103004.htm [3] 黄筑平, 杨黎明, 潘客麟.材料的动态损伤和失效[J].力学进展, 1993, 23(4):433-467. http://d.old.wanfangdata.com.cn/Periodical/gthjjs201304023Huang Zhuping, Yang Liming, Pan Kelin. Dynamic damage and failure of materials[J]. Advances in Mechanics, 1993, 23(4):433-467. http://d.old.wanfangdata.com.cn/Periodical/gthjjs201304023 [4] Zhang L, Cai L C, Li Y L, et al. Simplified model for prediction of dynamic damage and fracture of ductile materials[J]. International Journal of Solids & Structures, 2004, 41(24):7063-7074. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=e7d9d8472fd17f65ef142c3ced83b26c [5] 张林.延性材料冲击响应: 动态损伤与断裂、结构相变的新模型[D].绵阳: 中国工程物理研究院, 2005. http://cdmd.cnki.com.cn/Article/CDMD-82818-2005136415.htm [6] 陈永涛, 唐小军, 李庆忠, 等.纯铁材料的冲击相变与"反常"层裂[J].爆炸与冲击, 2009, 29(6):638-641. doi: 10.11883/1001-1455(2009)06-0637-05Chen Yongtao, Tang Xiaojun, Li Qingzhong, et al. Phase transition and abnormal spallation in pure iron[J]. Explosion and Shock Waves, 2009, 29(6):638-641. doi: 10.11883/1001-1455(2009)06-0637-05 [7] 陈永涛, 唐小军, 李庆忠.Fe基α相合金的冲击相变及其对层裂行为的影响研究[J].物理学报, 2011, 60(4):486-494. http://d.old.wanfangdata.com.cn/Periodical/wlxb201104079Chen Yongtao, Tang Xiaojun, Li Qingzhong. Phase transition and influence of phase transition on spall in α phase Fe-based alloy[J]. Acta Physica Sinica, 2011, 60(4):486-494. http://d.old.wanfangdata.com.cn/Periodical/wlxb201104079 [8] 章杰, 苏少卿, 郑宇, 等.改进SPH方法在陶瓷材料层裂数值模拟中的应用[J].爆炸与冲击, 2013, 33(4):401-407. doi: 10.3969/j.issn.1001-1455.2013.04.011Zhang Jie, Su Shaoqing, Zheng Yu, et al. Application of modified SPH method to numerical simulation of ceramic spallation[J]. Explosion and Shock Waves, 2013, 33(4):401-407. doi: 10.3969/j.issn.1001-1455.2013.04.011 [9] 彭辉, 李平, 裴晓阳, 等.动态损伤演化的空间不连续性实验研究[J].物理学报, 2013, 62(22):226201. doi: 10.7498/aps.62.226201Peng Hui, Li Ping, Pei Xiaoyang, et al. Experimental study of the spatial discontinuity of dynamic damage evolution[J]. Acta Physica Sinica, 2013, 62(22):226201. doi: 10.7498/aps.62.226201 [10] 贺红亮.动态拉伸断裂的物理判据研究[J].高压物理学报, 2013, 27(2):153-161. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gywlxb201302001He Hongliang. Physical criterion of dynamic tensile fracture[J]. Chinese Journal of High Pressure Physics, 2013, 27(2):153-161. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gywlxb201302001 [11] 张林, 张祖根, 秦晓云, 等.D6A、921和45钢的动态破坏与低压冲击特性[J].高压物理学报, 2003, 17(4):305-310. doi: 10.3969/j.issn.1000-5773.2003.04.011Zhang Lin, Zhang Zugen, Qin Xiaoyun, et al. Dynamic fracture and mechanical property of D6A, 921 and 45 steels under low shock pressure[J]. Chinese Journal of High Pressure Physics, 2003, 17(4):305-310. doi: 10.3969/j.issn.1000-5773.2003.04.011 [12] 翟少栋, 李英华, 彭建祥, 等.平面碰撞与强激光加载下金属铝的层裂行为[J].爆炸与冲击, 2016, 36(6):767-773. doi: 10.11883/1001-1455(2016)06-0767-07Zhai Shaodong, Li Yinghua, Peng Jianxiang, et al. Spall behavior of pure aluminum under plate-impact and high energy laser shock loadings[J].Explosion and Shock Waves, 2016, 36(6):767-773. doi: 10.11883/1001-1455(2016)06-0767-07 [13] Barker L M, Hollenbach R E. Interferometer technique for measuring the dynamic mechanical properties of materials[J]. Review of Scientific Instruments, 1965, 36(11):1617-1620. doi: 10.1063/1.1719405 [14] 张林, 李英华, 陈大年, 等.金属材料层裂强度测量方法: GJB7368─2011[S].绵阳: 中国工程物理研究院, 2011. [15] 张万甲, 杨中正.93钨合金断裂特性研究[J].高压物理学报, 1995, 9(4):279-288. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199500150295Zhang Wanjia, Yang Zhongzheng. Studies on the fracture behaviour for 93 tungsten alloy[J]. Chinese Journal of High Pressure Physics, 1995, 9(4):279-288. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199500150295 [16] Signor L, De Rességuier T, Dragon A, et al. Investigation of fragments size resulting from dynamic fragmentation in melted state of laser shock-loaded tin[J]. International Journal of Impact Engineering, 2010, 37(8):887-900. doi: 10.1016/j.ijimpeng.2010.03.001 [17] De Rességuier T, Loison D, Lescoute E, et al. Dynamic fragmentation of laser shock-melted metals: Some experimental advances[J]. Journal of Theoretical & Applied Mechanics, 2010, 48(4):957-972. [18] Luo S N, An Q, Germann T C, et al. Shock-induced spall in solid and liquid Cu at extreme strain rates[J]. Journal of Applied Physics, 2009, 106(1):013502. doi: 10.1063/1.3158062 [19] Holtkamp D B, Clark D A, Ferm E N, et al. A survey of high explosive-induced damage and spall in selected metals using proton radiography[C]//Furnish M D, Gupta Y M, Forbes J M. Shock Compression of Confensed Matter-2003. New York: American Institute of Physics, 2004: 477-482. [20] Lescoute E, De Rességuier T, Chevalier J M, et al. Transverse shadowgraphy and new recovery techniques to investigate dynamic fragmentation of laser shock-loaded metals[C]//Elert M L, Buttler W T, Furnish M D, et al.Shock Compression of Confensed Matter-2009. New York: American Institute of Physics, 2009: 1043-1046. [21] McCluskey C W, Wilke M D, Anderson W W, et al. Asay window: A new spall diagnostic[J]. Review of Scientific Instruments, 2006, 77(11):477-482. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ022448623/ [22] Signor L, Roy G, Chanal P Y, et al. Debris cloud ejection from shock-loaded tin melted on release or on compression[C]// Elert M L, Buttler W T, Furnish M D, et al. Shock Compression of Confensed Matter-2009. New York: American Institute of Physics, 2009: 1065-1068. [23] 陈永涛, 任国武, 汤铁钢, 等.爆轰加载下金属样品的熔化破碎现象诊断[J].物理学报, 2013, 62(11):116202. doi: 10.7498/aps.62.116202Chen Yongtao, Ren Guowu, Tang Tiegang, et al. Experimental diagnostic of melting fragments under explosive loading[J]. Acta Physica Sinica, 2013, 62(11):116202. doi: 10.7498/aps.62.116202 [24] 张林, 李雪梅, 李英华.锡的微层裂破碎复合Asay窗实验研究[C]//第五届全国强动载效应及防护学术会议暨复杂介质/结构的动态力学行为创新研究群体学术研讨会论文集.北京: 北京理工大学, 2013: 58-63. [25] Zhang L, Li Y H, Yu Y Y, et al. General construction of mean-field potential and its application to the multiphase equations of state of tin[J]. Physica B, 2011, 406(22):4163-4169. doi: 10.1016/j.physb.2011.01.018 [26] Marsh S P. LASL shock Hugoniot data[M]. Berkeley:University of California Press, 1980. [27] De Rességuier T, Signor L, Dragon A, et al. On the dynamic fragmentation of laser shock-melted tin[J]. Applied Physics Letters, 2008, 92(13):013506. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=9666c07b6bdf4a3874acb4f34361ac0d 期刊类型引用(2)
1. 张凤国,刘军,何安民,赵福祺,王裴. 强冲击加载下延性金属卸载熔化损伤/破碎问题的物理建模及其应用. 物理学报. 2022(24): 284-292 . 百度学术
2. 贺年丰,任国武,陈永涛,郭昭亮. 爆轰加载下金属锡层裂破碎数值模拟. 爆炸与冲击. 2019(04): 91-98 . 本站查看
其他类型引用(1)
-