Study on the scaling law of geometrically-distorted thin-walled cylindrical shells subjected to axial impact
-
摘要: 对于受轴向冲击载荷作用的薄壁圆管动态响应的相似律问题,由于圆管的薄壁特性导致厚度无法与高度和半径按相同的比例进行结构缩放,从而产生模型的几何畸变,此时传统的相似律已无法描述原型与畸变模型之间的动态响应规律。基于薄壁圆管轴向冲击问题的控制方程,通过能量守恒和量纲分析,推导了考虑几何畸变条件下轴向冲击载荷作用的理想弹塑性薄壁圆管动态响应的相似律。通过在给定应变与应变率区间上建立比例模型预测的流动屈服应力与原型流动屈服应力的最佳逼近关系,将几何畸变相似律进一步推广至包含应变率和应变硬化的材料。通过数值方法验证了提出的几何畸变模型相似律的适用性。分析结果表明,提出的考虑厚度畸变的受轴向冲击薄壁圆管的相似律可用于预测原型结构的冲击动态响应,并显著降低比例模型与原型结构平均载荷和能量的偏差。Abstract: In scaling the dynamic responses of thin-walled cylindrical shells subjected to axial impact loading, the thickness cannot be adjusted according to the same scale as the radius and height due to the thin wall characteristics. Hence, geometrically-distorted models would be used, and the traditional scaling law cannot describe the relationship between the dynamic responses of the prototype and the geometrically-distorted model. In this paper, the scaling law for this case was derived for elastic-ideal plastic thin-walled cylindrical shells under axial impact loading. For strain hardening and strain-rate hardening material, based on the average load, deformation energy, and displacement of the shell in the axisymmetric deformation mode, the dimensionless numbers of three key design parameters, namely the stress, mass, and displacement, were obtained through the law of energy conservation. Then, the optimal approximation of the flow stress predicted by the distorted scaled model to the flow stress of the prototype was established on a given strain and strain rate interval. In this way, the derived scaling law can be applied to the case considering the coupling effects of geometric distortion, strain-rate sensitivity, and strain hardening. Finally, several finite element models of thin-walled cylindrical shell models subject to axial mass impact were established. These models use the elastic-ideal plastic material model and the general material model with strain-rate hardening and strain hardening effects. The modified impact velocity and impact mass were obtained by the present method using the geometrically-distorted model, which verified the effectiveness and correctness of the proposed scaling law. The results show that the geometrically distorted model corrected by the method proposed in this article can quite accurately predict the dynamic responses of the prototype, and significantly reduce the errors in the dynamic responses of the thin-walled cylindrical shell subjected to axial impact loading, especially the average load and deformation energy.
-
随着计算机技术的发展,数值模拟已经成为复杂系统研究及设计的重要手段。以往,研究人员常常将材料参数、几何尺寸等作为确定性的输入参数,代入程序进行数值模拟,从而得到确定性的结果。然而实际上,数值模拟过程中存在各种不确定性来源。近十几年,不确定度量化(uncertainty quantification, UQ)方法发展迅速,逐渐成为数值模拟研究的热点[1-2]。不确定度量化主要关注如何表征、量化、减小数值模拟过程中的各种不确定性来源,是实现高置信度数值模拟的有效手段。在进行不确定度量化分析时,为了研究输入参数引入的不确定性,每个输入参数常常需要人为给定概率分布。输入参数概率分布的选择很可能会对分析结果有直接影响,为了减少分布选择的主观性,输入参数的概率分布最好能够由试验数据提供。
目前,圆筒试验是评估炸药的做功能力、确定爆轰产物状态方程参数的基准试验,已经得到了广泛应用[3-8]。标准的圆筒试验方法主要采用电离探针测量炸药爆速,用高速扫描相机记录定常滑移爆轰驱动下圆筒壁的径向膨胀距离随时间的变化关系,并利用经验公式计算圆筒壁的膨胀速度和比动能等表征炸药做功能力的特征量,最后利用这些试验结果标定JWL(Jones-Wilkins-Lee)状态方程的参数[9-10]。本文中,利用贝叶斯方法对圆筒试验JWL状态方程参数进行标定,给出标定参数的估计值,得到标定参数的后验分布,然后将该后验分布作为不确定度量化分析时输入参数的初始分布,以减小输入参数分布选择引入的不确定性。
1. 圆筒试验
圆筒试验的试验装置示意图如图 1所示,高速相机VISAR从狭缝扫描记录的是圆筒壁外表面的径向膨胀过程,通过对照片底片的处理获得圆筒壁的径向膨胀位移和时间的关系。得到的观测数据如图 2所示。炸药爆轰产物JWL状态方程形式和等熵条件如下:
{p=A[1−ω/(R1V)]e−R1V+B[1−ω/(R2V)]e−R2V+ωe/VpS=Ae−R1V+Be−R2V+C/Vω+1 (1) 式中:p为压力,下标“S”表示熵,V为相对比容,e为初始体积能量,A、B、C、R1、R2和ω为待标定参数。利用JWL状态方程在C-J点的特性,由圆筒试验和相关试验的观测数据得到爆压、爆速、密度、内能和参数ω,再构造代数方程组,给定R1和R2,联立求解A、B和C[11]。
2. 贝叶斯标定方法
M.Kennedy和A.O‘Hagan[12]提出用贝叶斯方法进行参数标定和预测。由此发展出来一系列方法,统称为KOH方法[13-15]。KOH方法用高斯过程(Gaussian process, GP)对计算结果和观测结果进行建模,同时对模型偏差进行建模。KOH方法的主要思想是用高斯过程建立代理模型,然后用代理模型计算后验分布和标定参数。下面简单介绍KOH方法的框架,更多技术细节参见文献[12]。
2.1 定义和统计模型
令X=(X(1), …, X(qX))为qX维控制变量,θ=(θ(1), …, θ(qθ))为qθ维标定参数,Y为数值模拟计算的输出结果,Z为试验观测结果。KOH模型如下:
yi=fC(x∗i,si)i=1,⋯,M (2) zj=ρfC(xj,θ)+b(xj)+εjj=1,⋯,N (3) 式中:fC(·)为计算模型,b(·)为模型偏差,εj为观测误差,ρ为回归系数,M为数值模拟样本量,N为试验样本量,xi*为控制变量在第i次数值模拟中的取值,xj为控制变量在第j次实验中的取值(按照数理统计中约定,采用大写字母表示随机变量,相应的小写字母表示该随机变量的取值,本文中,“X”、“Y”、“Z”、“T”、“R1”、“R2”均采用该约定)。为了得到数值模拟的结果,需要给标定参数赋值,令s表示数值模拟时输入的标定参数,θ表示标定参数的真实值,以便于区分。
假设εj服从均值为零、方差为λ的正态分布;fC(·)用高斯过程建模,均值函数为m1(x, s),协方差函数为c1(·, ·);模型偏差b(x)也用高斯过程建模,均值函数为m2(x),协方差函数为c2(·, ·)。假设均值函数和协方差函数的形式如下:
mi(⋅)=hi(⋅)Tβii=1,2 (4) ci((xk,sk),(xl,sl))=σ2iexp{−qX∑j=1ωXij(x(j)k−x(j)l)2−qθ∑r=1ωθir(s(r)k−s(r)l)2}i=1,2 (5) 式中:hi(·)为统计建模时任意给定的连接函数,σi2、ωXi、ωθi为协方差函数建模中需要估计的未知参数。令ψ=(ψ1, ψ2),ψi是协方差函数ci(·, ·)中的超参数(σi2, ωXi, ωθi),i=1, 2。KOH模型的超参数包括β=(β1, β2)和φ=(ρ, λ, ψ), 则全部要估计的参数为(θ, β, φ)。
2.2 数据结构
数据包括观测数据z和数值模拟结果y, 即dT=(zT, yT)。数值模拟过程中输入变量的取值记为D1={(x1*, s1), …, (xM*, sM)},试验数据控制变量的取值记为D2={x1, …, xN},与D1相对应,可以将D2扩充为D2(θ)={(x1, θ), …, (xN, θ)}。
2.3 后验分布
令先验分布P(β)∝1。根据贝叶斯理论,在给定先验分布的条件下,后验分布为:
P(θ,β,φ|d)∝P(θ)P(φ)f(d;md(θ),Vd(θ)) (6) 式中:f(d; md(θ), Vd(θ))是多元正态分布N(md(θ), Vd(θ))的密度函数,且
md(θ)=E(d|θ,β,φ)=H(θ)β=(H1(D1)0ρH1(D2(θ))H2(D2))β (7) Vd(θ)=V(d|θ,β,φ)=(V1(D1)ρC1(D1,D2(θ))TρC1(D1,D2(θ))λIN+ρ2V1(D2(θ))+V2(D2)) (8) 式中:H1、H2为由连接函数hi构成的向量,V1(D1)为由集合D1中数据构成的方差矩阵,V2(D2)为由集合D2中数据构成的方差矩阵,C1为D1和D2中数据构成的协方差矩阵,IN为N维单位矩阵。
2.4 计算步骤
直接计算式(6)中的后验分布较困难,因此将计算过程分为4个步骤:步骤1和步骤2是建立GP模型,步骤3是用建立的GP模型计算后验分布,步骤4是利用步骤3得到的后验分布标定参数。
步骤1:建立GP代理模型,估计β1和c1(·, ·)的超参数ψ1。用数值模拟所得数据y估计β1和ψ1,得到估计值ˆβ1和ˆψ1。步骤2:建立观测数据的GP代理模型,估计ρ、λ、β2和c2(·, ·)的超参数ψ2。在ˆβ1和ˆψ1给定的条件下,用观测数据z估计ρ、λ、β2和c2(·)的超参数ψ2,从而得到估计值ˆβ=(ˆβ1,ˆβ2)和ˆφ=(ˆρ,ˆλ,ˆψ)。步骤3:计算标定参数θ的后验分布。由于P(θ|φ=ˆφ,d)∝P(θ,ˆφ|d),因此通过计算P(θ,ˆφ|d)即可得到θ的后验分布。步骤4:利用θ的后验分布标定参数。标定参数的估计可以定义为后验分布的均值ˆθ(mean)、中位数ˆθ(mid)或最大概率值ˆθ(max),其中ˆθ(max)也称为最大后验估计(maximum posterior estimates, MPE)。
3. JWL状态方程参数的标定
采用KOH贝叶斯方法标定圆筒试验JWL状态方程参数。将标定参数记为(R1, R2);将时间视为设计变量,记为T;同时令M=30, N=10。数值计算过程中,R1和R2的抽样区间分别为[4.3, 5.5]和[0.8, 1.5],用拉丁超立方抽样(latin hypercube sample, LHS)方法得到M组参数并代入模拟程序进行计算,得到M条径向半径随时间变化的曲线。对设计变量T进行LHS抽样得到M个样本。将T和(R1, R2)的样本随机配对得到模拟计算的输入数据D1={(t1*, R11, R21), …, (tM*, R1M, R2M)}。D1对应的计算输出结果为y={y1, …, yM}。需要注意的是,每次数值计算的结果都是径向半径随时间变化的曲线,这里只是从曲线上随机取一个抽样点作为输出结果。对时间变量T随机抽样得到N个样本,即D2={t1, …, tN}。将D2扩充为D2(r1, r2)={(t1, r1, r2), …, (tN, r1, r2)},其中(r1, r2)是标定参数的真实值。D2(r1, r2)对应的径向半径观测数据为z={z1, …, zN}。假设(R1, R2)的先验分布为均匀分布,即对(R1, R2)没有任何知识;超参数ψ的先验分布为标准正态分布。数值模拟使用的计算程序为由北京应用物理与计算数学研究所自主研制的相容拉氏流体动力学工具包CHAP (compatible hydrodynamics analysis program)[16]。
下面分别采用两个模型进行建模,其主要区别在于模型中是否包含偏差。模型1是KOH模型的通用形式,如式(2)~(3)所示;模型2是KOH模型的特殊形式,相当于在模型1的基础上令ρ=1,b(xj)≡0,此时式(3)可以简化为:
zj=fC(xj,θ)+εjj=1,⋯,N (9) 按照第2.4节中的4个步骤依次进行计算。第1步是建立数值模拟的GP代理模型。GP代理模型中参数的估计值分别为ˆβ1=3.090,ˆσ21=2.296,ˆωX1=0.603,ˆωR11=0.064,ˆωR21=0.059。为检验GP代理模型预测的能力,重新随机抽样4组(R1, R2)样本代入数值模拟程序进行计算,将计算得到的结果和采用GP模型预测的结果进行比较,如图 3所示。从图 3可以看出:代理模型预测曲线和数值计算的曲线基本重合,说明代理模型的预测足够精确。第2步是建立观测数据的GP代理模型。模型1中模型参数的估计值分别为ˆρ=1.057,ˆλ=1.810,ˆβ2=−0.607,ˆσ22=1.585,ˆωX2=2.220。由于模型2忽略了模型偏差并且ρ=1,所以该步骤中只需要估计λ,估计值为ˆλ=2.608×10−5。第3步是计算标定参数的后验分布。图 4给出了标定参数(R1, R2)的后验概率分布等高线,其中后验概率共分为9个概率区间,分别由9种颜色进行区分,从小到大依次用Pk (k=1, …, 9)表示,P9表示最大概率区间的值域。图 4(a)为采用模型1得到的后验概率分布,可以看出,概率分布从右下角向左上角逐级递增,概率最大区域P9位于左上角;图 4(b)为采用模型2得到的后验概率分布,可以看出,概率分布从右下角和左边向中间区域逐级递增,P9近似位于中间位置。第4步是标定参数。由图 4可知:模型1的MPE为ˆR(max)1=4.68258,ˆR(max)2=1.31066;模型2的MPE为ˆR(max)1=4.74985,ˆR(max)2=1.10314。两个模型的MPE均近似位于P9区域的中心位置。
图 5为模型1和模型2的预测效果,其中标定结果(红线)是将MPE代入数值模拟程序计算得到,预测结果(绿线)采用观测数据的GP代理模型计算得到,试验结果用黑线表示。图 5显示:两个模型的预测结果均与试验结果符合很好,说明建立的代理模型都能够比较精确地描述试验数据。模型2中试验结果和标定结果符合得比较好,表明标定的参数值是最佳拟合的估计;但是,当t>39μs时,模型1的标定结果和试验结果差异逐渐增大。由于数值模拟和基于观测数据代理模型的预测(即ˆfC(xj,ˆθ)和ˆfC(xj,ˆθ)+ˆb(xj))都足够精确,所以两者的差异即为模型偏差的预测ˆb(x),说明模型1标定的参数值是“调和”了模型偏差ˆb(x)和标定参数^θ 的最佳估计。
4. 结论
(1) 用KOH贝叶斯方法标定圆筒试验JWL状态方程的参数R1和R2,得到了标定参数的估计值和后验分布,并将后验分布作为子系统级(或系统级)复杂系统数值模拟不确定度量化分析的输入参数的初始分布,以尽可能地减少分布选择引入的认知不确定性。(2)分析了两种KOH统计模型假设对标定参数估计值和后验分布的影响。从参数标定的角度分析,模型2得到的估计值更好;从不确定度量化的角度分析,由于模型1还分析了模型偏差引入的不确定性,因此模型1更合适。
-
表 2 受轴向冲击的理想弹塑性薄壁圆管比例因子
Table 2. Scaling factors of the elastic-ideal plastic thin-walled cylindrical shell under axial impact loading
变量 比例因子 变量 比例因子 长度L β=Lm/Lp 位移δ βδ=β 密度ρ βρ=ρm/ρp 应力σd βσd=βρβ2v√β/βh 速度v βv=vm/vp 应变ε βε=√βh/β 质量m βm=βρβ2βh 应变率˙ε β˙ε=(βv/β)√βL/βh 时间t βt=β/βv 载荷P βP=βρββhβ2v 加速度 a βa=β2v/β 动能Ek βEk=βρβhβ2β2v ρ/(g∙cm−3) E/GPa μ A/MPa B/MPa C n ˙ε0/s−1 7.89 207 0.3 350 275 0.022 0.36 1 表 4 理想弹塑性模型几何畸变比例因子
Table 4. Scaling factors of geometrically-distorted models of elastic-ideal plastic material
β η βh βM βv βt βP βδ βε 0.1 1.2 0.12 0.12 1.0466 0.0955 0.0131 0.1 1.0954 0.1 1.5 0.15 0.15 1.1067 0.0904 0.0184 0.1 1.2247 0.1 1.8 0.18 0.18 1.1583 0.0863 0.0241 0.1 1.3416 0.1 2.0 0.20 0.20 1.1892 0.0841 0.0283 0.1 1.4142 表 5 考虑应变率效应和应变硬化效应几何畸变模型比例因子
Table 5. Scaling factors of geometrically-distorted models considering strain-rate sensitivity and strain hardening
β η βh βM βv βt βP βδ βε 0.1 1.2 0.12 0.12 1.0756 0.0930 0.0139 0.1 1.0954 0.1 1.5 0.15 0.15 1.1442 0.0874 0.0196 0.1 1.2247 0.1 1.8 0.18 0.18 1.2038 0.0831 0.0261 0.1 1.3416 0.1 2.0 0.20 0.20 1.2396 0.0807 0.0306 0.1 1.4142 表 6 比例模型与原型的峰值位移和平均载荷相对误差
Table 6. Relative errors in the peak displacement and average force between the scale models and prototype
模型 (δ/βδ)/mm 相对误差/% (P/βP)/kN 相对误差/% 原型 124.258 − 71.95 − β=1/10, η=1.2 121.669 2.084 73.43 2.057 β=1/10, η=1.5 120.651 2.903 69.48 3.433 β=1/10, η=1.8 118.642 4.520 75.85 5.420 β=1/10, η=2.0 117.678 5.295 72.25 0.417 表 7 考虑应变率与应变硬化效应的比例模型位移与平均载荷相对误差
Table 7. Relative errors in the peak displacement and average force of the scaled models considering strain-rate sensitivity and strain hardening
模型 (δ/βδ) /mm 相对误差/% (P/βP)/kN 相对误差/% 原型 90.740 − 97.79 − β=1/10, η=1.2 89.619 1.235 98.42 0.644 β=1/10, η=1.5 89.222 1.673 99.47 1.718 β=1/10, η=1.8 86.201 5.002 101.06 3.344 β=1/10, η=2.0 86.232 4.968 105.14 7.516 -
[1] JONES N. Structural impact [M]. 2nd ed. New York: Cambridge University Press, 2012. [2] 徐海斌, 张德志, 谭书舜, 等. 轴向压缩的金属薄壁圆管相似律的实验研究 [C] // 第20届全国结构工程学术会议论文集. 浙江宁波: 中国力学学会工程力学编辑部, 2011: 554–559. [3] ALEXANDER J M. An approximate analysis of the collapse of thin cylindrical shells under axial loading [J]. The Quarterly Journal of Mechanics and Applied Mathematics, 1960, 13(1): 10–15. DOI: 10.1093/qjmam/13.1.10. [4] ABRAMOWICZ W, JONES N. Dynamic axial crushing of circular tubes [J]. International Journal of Impact Engineering, 1984, 2(3): 263–281. DOI: 10.1016/0734-743X(84)90010-1. [5] KARAGIOZOVA D, JONES N. Influence of stress waves on the dynamic progressive and dynamic plastic buckling of cylindrical shells [J]. International Journal of Solids and Structures, 2001, 38(38/39): 6723–6749. DOI: 10.1016/S0020-7683(01)00111-1. [6] KARAGIOZOVA D, NURICK G N, YUEN S C K. Energy absorption of aluminium alloy circular and square tubes under an axial explosive load [J]. Thin-Walled Structures, 2005, 43(6): 956–982. DOI: 10.1016/j.tws.2004.11.002. [7] LU G, YU J L, ZHANG J J, et al. Alexander revisited: upper- and lower-bound approaches for axial crushing of a circular tube [J]. International Journal of Mechanical Sciences, 2021, 206: 106610. DOI: 10.1016/j.ijmecsci.2021.106610. [8] CASABURO A, PETRONE G, FRANCO F, et al. A review of similitude methods for structural engineering [J]. Applied Mechanics Reviews, 2019, 71(3): 030802. DOI: 10.1115/1.4043787. [9] COUTINHO C P, BAPTISTA A J, RODRIGUES J D. Reduced scale models based on similitude theory: a review up to 2015 [J]. Engineering Structures, 2016, 119: 81–94. DOI: 10.1016/j.engstruct.2016.04.016. [10] OSHIRO R E, ALVES M. Scaling impacted structures [J]. Archive of Applied Mechanics, 2004, 74(1/2): 130–145. DOI: 10.1007/BF02637214. [11] OSHIRO R E, ALVES M. Scaling of cylindrical shells under axial impact [J]. International Journal of Impact Engineering, 2007, 34(1): 89–103. DOI: 10.1016/j.ijimpeng.2006.02.003. [12] 王帅, 徐绯, 代震, 等. 结构冲击畸变问题的直接相似方法研究 [J]. 力学学报, 2020, 52(3): 774–786. DOI: 10.6052/0459-1879-19-327.WANG S, XU F, DAI Z, et al. A direct scaling method for the distortion problems of structural impact [J]. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(3): 774–786. DOI: 10.6052/0459-1879-19-327. [13] WANG S, XU F, ZHANG X Y, et al. Material similarity of scaled models [J]. International Journal of Impact Engineering, 2021, 156: 103951. DOI: 10.1016/j.ijimpeng.2021.103951. [14] 李肖成, 徐绯, 杨磊峰, 等. 薄板在冲击载荷下线弹性理想塑性响应的相似性研究 [J]. 爆炸与冲击, 2021, 41(11): 113103. DOI: 10.11883/bzycj-2020-0374.LI X C, XU F, YANG L F, et al. Study on the similarity of elasticity and ideal plasticity response of thin plate under impact loading [J]. Explosion and Shock Waves, 2021, 41(11): 113103. DOI: 10.11883/bzycj-2020-0374. [15] 秦健, 张振华. 原型和模型不同材料时加筋板冲击动态响应的相似预报方法 [J]. 爆炸与冲击, 2010, 30(5): 511–516. DOI: 10.11883/1001-1455(2010)05-0511-06.QIN J, ZHANG Z H. A scaling method for predicting dynamic responses of stiffened plates made of materials different from experimental models [J]. Explosion and Shock Waves, 2010, 30(5): 511–516. DOI: 10.11883/1001-1455(2010)05-0511-06. [16] ALVES M, OSHIRO R E, CALLE M A G, et al. Scaling and structural impact [J]. Procedia Engineering, 2017, 173: 391–396. DOI: 10.1016/j.proeng.2016.12.036. [17] MAZZARIOL L M, ALVES M. Similarity laws of structures under impact load: geometric and material distortion [J]. International Journal of Mechanical Sciences, 2019, 157/158: 633–647. DOI: 10.1016/j.ijmecsci.2019.05.011. [18] WANG S, XU F, DAI Z. Suggestion of the DLV dimensionless number system to represent the scaled behavior of structures under impact loads [J]. Archive of Applied Mechanics, 2020, 90(4): 707–719. DOI: 10.1007/s00419-019-01635-9. [19] WANG S, XU F, ZHANG X Y, et al. A directional framework of similarity laws for geometrically distorted structures subjected to impact loads [J]. International Journal of Impact Engineering, 2022, 161: 104092. DOI: 10.1016/j.ijimpeng.2021.104092. [20] 李志斌, 虞吉林, 郑志军, 等. 薄壁管及其泡沫金属填充结构耐撞性的实验研究 [J]. 实验力学, 2012, 27(1): 77–86.LI Z B, YU J L, ZHENG Z J, et al. An experimental study on the crashworthiness of thin-walled tubes and their metallic foam-filled structures [J]. Journal of Experimental Mechanics, 2012, 27(1): 77–86. [21] 朱文波, 杨黎明, 余同希. 薄壁圆管轴向冲击下的动态特性研究 [J]. 宁波大学学报(理工版), 2014, 27(2): 92–96.ZHU W B, YANG L M, YU T X. Study on dynamic properties of thin-walled circular tubes under axial compression [J]. Journal of Ningbo University (Natural Science & Engineering Edition), 2014, 27(2): 92–96. [22] 余同希, 卢国兴, 张雄. 能量吸收: 结构与材料的力学行为和塑性分析 [M]. 北京: 科学出版社, 2019. [23] 白以龙, 黄筑平, 虞吉林, 等. 材料和结构的动态响应 [M]. 合肥: 中国科学技术大学出版社, 2005. [24] JOHNSON G R, COOK W H. A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures [C] // Proceedings of the 7th International Symposium on Ballistics. Hague, Netherlands, 1983: 541–547. 期刊类型引用(4)
1. Songmao Zhao,Haotian Jian,Ke Wang,Zheng Ning,Peng Zhu,Ruiqi Shen. Investigation on the exploding foil initiators ignition enhanced by explosion-electricity coupling. Defence Technology. 2024(12): 1-11 . 必应学术
2. 薛创,宁成,彭先觉. 水中金属丝电爆炸动力学过程的零维模型. 爆炸与冲击. 2023(05): 150-157 . 本站查看
3. 王万军,孙秀娟,张雷,雷凡,郭菲,杨爽,付秋菠. 亚微秒内金丝电爆炸的激光探针测试及三维动力学仿真. 含能材料. 2019(06): 473-480 . 百度学术
4. 张永民,姚伟博,邱爱慈,汤俊萍,王宇,呼义翔. 金属丝电爆炸现象研究综述. 高电压技术. 2019(08): 2668-2680 . 百度学术
其他类型引用(1)
-