Loading [MathJax]/jax/output/HTML-CSS/jax.js
  • ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

受轴向冲击薄壁圆管的几何畸变相似律研究

杨磊峰 常新哲 徐绯 王帅 刘小川 惠旭龙 李肖成

陈华, 周海兵, 刘国昭, 孙占峰, 张树道. 圆筒试验JWL状态方程参数的贝叶斯标定[J]. 爆炸与冲击, 2017, 37(4): 585-590. doi: 10.11883/1001-1455(2017)04-0585-06
引用本文: 杨磊峰, 常新哲, 徐绯, 王帅, 刘小川, 惠旭龙, 李肖成. 受轴向冲击薄壁圆管的几何畸变相似律研究[J]. 爆炸与冲击, 2022, 42(5): 053205. doi: 10.11883/bzycj-2021-0452
Chen Hua, Zhou Haibing, Liu Guozhao, Sun Zhanfeng, Zhang Shudao. Bayesian calibration for parameters of JWL equation of state in cylinder test[J]. Explosion And Shock Waves, 2017, 37(4): 585-590. doi: 10.11883/1001-1455(2017)04-0585-06
Citation: YANG Leifeng, CHANG Xinzhe, XU Fei, WANG Shuai, LIU Xiaochuan, XI Xulong, LI Xiaocheng. Study on the scaling law of geometrically-distorted thin-walled cylindrical shells subjected to axial impact[J]. Explosion And Shock Waves, 2022, 42(5): 053205. doi: 10.11883/bzycj-2021-0452

受轴向冲击薄壁圆管的几何畸变相似律研究

doi: 10.11883/bzycj-2021-0452
基金项目: 国家自然科学基金(11972309);中央高校基本科研业务费专项资金(310201901A012);高等学校学科创新引智计划(111计划)(BP0719007)
详细信息
    作者简介:

    杨磊峰(1997- ),男,硕士研究生,409332963@mail.nwpu.edu.cn

    通讯作者:

    徐 绯(1970- ),女,博士,教授,xufei@nwpu.edu.cn

  • 中图分类号: O347.3;V214.4

Study on the scaling law of geometrically-distorted thin-walled cylindrical shells subjected to axial impact

  • 摘要: 对于受轴向冲击载荷作用的薄壁圆管动态响应的相似律问题,由于圆管的薄壁特性导致厚度无法与高度和半径按相同的比例进行结构缩放,从而产生模型的几何畸变,此时传统的相似律已无法描述原型与畸变模型之间的动态响应规律。基于薄壁圆管轴向冲击问题的控制方程,通过能量守恒和量纲分析,推导了考虑几何畸变条件下轴向冲击载荷作用的理想弹塑性薄壁圆管动态响应的相似律。通过在给定应变与应变率区间上建立比例模型预测的流动屈服应力与原型流动屈服应力的最佳逼近关系,将几何畸变相似律进一步推广至包含应变率和应变硬化的材料。通过数值方法验证了提出的几何畸变模型相似律的适用性。分析结果表明,提出的考虑厚度畸变的受轴向冲击薄壁圆管的相似律可用于预测原型结构的冲击动态响应,并显著降低比例模型与原型结构平均载荷和能量的偏差。
  • 随着计算机技术的发展,数值模拟已经成为复杂系统研究及设计的重要手段。以往,研究人员常常将材料参数、几何尺寸等作为确定性的输入参数,代入程序进行数值模拟,从而得到确定性的结果。然而实际上,数值模拟过程中存在各种不确定性来源。近十几年,不确定度量化(uncertainty quantification, UQ)方法发展迅速,逐渐成为数值模拟研究的热点[1-2]。不确定度量化主要关注如何表征、量化、减小数值模拟过程中的各种不确定性来源,是实现高置信度数值模拟的有效手段。在进行不确定度量化分析时,为了研究输入参数引入的不确定性,每个输入参数常常需要人为给定概率分布。输入参数概率分布的选择很可能会对分析结果有直接影响,为了减少分布选择的主观性,输入参数的概率分布最好能够由试验数据提供。

    目前,圆筒试验是评估炸药的做功能力、确定爆轰产物状态方程参数的基准试验,已经得到了广泛应用[3-8]。标准的圆筒试验方法主要采用电离探针测量炸药爆速,用高速扫描相机记录定常滑移爆轰驱动下圆筒壁的径向膨胀距离随时间的变化关系,并利用经验公式计算圆筒壁的膨胀速度和比动能等表征炸药做功能力的特征量,最后利用这些试验结果标定JWL(Jones-Wilkins-Lee)状态方程的参数[9-10]。本文中,利用贝叶斯方法对圆筒试验JWL状态方程参数进行标定,给出标定参数的估计值,得到标定参数的后验分布,然后将该后验分布作为不确定度量化分析时输入参数的初始分布,以减小输入参数分布选择引入的不确定性。

    圆筒试验的试验装置示意图如图 1所示,高速相机VISAR从狭缝扫描记录的是圆筒壁外表面的径向膨胀过程,通过对照片底片的处理获得圆筒壁的径向膨胀位移和时间的关系。得到的观测数据如图 2所示。炸药爆轰产物JWL状态方程形式和等熵条件如下:

    图  1  圆筒试验的试验装置示意图
    Figure  1.  Device of cylinder test
    图  2  圆筒试验径向半径和时间的关系
    Figure  2.  Radius vs. time in cylinder test
    {p=A[1ω/(R1V)]eR1V+B[1ω/(R2V)]eR2V+ωe/VpS=AeR1V+BeR2V+C/Vω+1 (1)

    式中:p为压力,下标“S”表示熵,V为相对比容,e为初始体积能量,ABCR1R2ω为待标定参数。利用JWL状态方程在C-J点的特性,由圆筒试验和相关试验的观测数据得到爆压、爆速、密度、内能和参数ω,再构造代数方程组,给定R1R2,联立求解ABC[11]

    M.Kennedy和A.O‘Hagan[12]提出用贝叶斯方法进行参数标定和预测。由此发展出来一系列方法,统称为KOH方法[13-15]。KOH方法用高斯过程(Gaussian process, GP)对计算结果和观测结果进行建模,同时对模型偏差进行建模。KOH方法的主要思想是用高斯过程建立代理模型,然后用代理模型计算后验分布和标定参数。下面简单介绍KOH方法的框架,更多技术细节参见文献[12]。

    X=(X(1), …, X(qX))为qX维控制变量,θ=(θ(1), …, θ(qθ))为qθ维标定参数,Y为数值模拟计算的输出结果,Z为试验观测结果。KOH模型如下:

    yi=fC(xi,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{qXj=1ωXij(x(j)kx(j)l)2qθr=1ωθir(s(r)ks(r)l)2}i=1,2 (5)

    式中:hi(·)为统计建模时任意给定的连接函数,σi2ωXiωθi为协方差函数建模中需要估计的未知参数。令ψ=(ψ1, ψ2),ψi是协方差函数ci(·, ·)中的超参数(σi2, ωXi, ωθi),i=1, 2。KOH模型的超参数包括β=(β1, β2)和φ=(ρ, λ, ψ), 则全部要估计的参数为(θ, β, φ)。

    数据包括观测数据z和数值模拟结果y, 即dT=(zT, yT)。数值模拟过程中输入变量的取值记为D1={(x1*, s1), …, (xM*, sM)},试验数据控制变量的取值记为D2={x1, …, xN},与D1相对应,可以将D2扩充为D2(θ)={(x1, θ), …, (xN, θ)}。

    令先验分布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)

    式中:H1H2为由连接函数hi构成的向量,V1(D1)为由集合D1中数据构成的方差矩阵,V2(D2)为由集合D2中数据构成的方差矩阵,C1D1D2中数据构成的协方差矩阵,INN维单位矩阵。

    直接计算式(6)中的后验分布较困难,因此将计算过程分为4个步骤:步骤1和步骤2是建立GP模型,步骤3是用建立的GP模型计算后验分布,步骤4是利用步骤3得到的后验分布标定参数。

    步骤1:建立GP代理模型,估计β1c1(·, ·)的超参数ψ1。用数值模拟所得数据y估计β1ψ1,得到估计值ˆβ1ˆψ1。步骤2:建立观测数据的GP代理模型,估计ρλβ2c2(·, ·)的超参数ψ2。在ˆβ1ˆψ1给定的条件下,用观测数据z估计ρλβ2c2(·)的超参数ψ2,从而得到估计值ˆβ=(ˆβ1,ˆβ2)ˆφ=(ˆρ,ˆλ,ˆψ)。步骤3:计算标定参数θ的后验分布。由于P(θ|φ=ˆφ,d)P(θ,ˆφ|d),因此通过计算P(θ,ˆφ|d)即可得到θ的后验分布。步骤4:利用θ的后验分布标定参数。标定参数的估计可以定义为后验分布的均值ˆθ(mean)、中位数ˆθ(mid)或最大概率值ˆθ(max),其中ˆθ(max)也称为最大后验估计(maximum posterior estimates, MPE)。

    采用KOH贝叶斯方法标定圆筒试验JWL状态方程参数。将标定参数记为(R1, R2);将时间视为设计变量,记为T;同时令M=30, N=10。数值计算过程中,R1R2的抽样区间分别为[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×105。第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区域的中心位置。

    图  3  GP代理模型预测的结果和数值计算结果的比较
    Figure  3.  Comparison of predictions of the GP surrogate model with numerical simulations
    图  4  两模型的后验概率等高线和最大后验估计
    Figure  4.  Posterior distribution and MPE with different models

    图 5为模型1和模型2的预测效果,其中标定结果(红线)是将MPE代入数值模拟程序计算得到,预测结果(绿线)采用观测数据的GP代理模型计算得到,试验结果用黑线表示。图 5显示:两个模型的预测结果均与试验结果符合很好,说明建立的代理模型都能够比较精确地描述试验数据。模型2中试验结果和标定结果符合得比较好,表明标定的参数值是最佳拟合的估计;但是,当t>39μs时,模型1的标定结果和试验结果差异逐渐增大。由于数值模拟和基于观测数据代理模型的预测(即ˆfC(xj,ˆθ)ˆfC(xj,ˆθ)+ˆb(xj))都足够精确,所以两者的差异即为模型偏差的预测ˆb(x),说明模型1标定的参数值是“调和”了模型偏差ˆb(x)和标定参数^θ 的最佳估计。

    图  5  径向半径的预测、试验和模拟结果的比较
    Figure  5.  Comparison of radius histories obtained by prediction, experiment and simulation

    (1) 用KOH贝叶斯方法标定圆筒试验JWL状态方程的参数R1R2,得到了标定参数的估计值和后验分布,并将后验分布作为子系统级(或系统级)复杂系统数值模拟不确定度量化分析的输入参数的初始分布,以尽可能地减少分布选择引入的认知不确定性。(2)分析了两种KOH统计模型假设对标定参数估计值和后验分布的影响。从参数标定的角度分析,模型2得到的估计值更好;从不确定度量化的角度分析,由于模型1还分析了模型偏差引入的不确定性,因此模型1更合适。

  • 图  1  轴向受压薄壁圆管的轴对称压溃模式

    Figure  1.  Axisymmetric crushing mode of the thin-walled cylindrical shell under axial compression

    图  2  金属材料应力-应变曲线

    Figure  2.  Stress-strain curves of metal materials

    图  3  受轴向冲击的薄壁圆管示意图

    Figure  3.  Schematic diagram of a thin-walled cylindrical shell under axial impact

    图  4  原型与比例模型的动态响应-时间曲线

    Figure  4.  Dynamic response-time curves of the scale models and the prototype

    图  5  修正后的比例模型与原型的动态响应-时间曲线

    Figure  5.  Dynamic response-time curves of the modified scaled models and the prototype

    表  1  纯几何相似律比例因子[10-11]

    Table  1.   Scaling factors by pure geometric similarity[10-11]

    变量比例因子变量比例因子
    长度L[10] β=Lm/Lp 位移δ[10] βδ=β
    密度ρ[10] βρ=ρm/ρp 应力σd[10] βσd=βρβ2v
    速度v[10] βv=vm/vp 应变ε[10] βε=1
    质量m[10] βm=βρβ3 应变率˙ε[10] β˙ε=βv/β
    时间t[10] βt=β/βv 载荷P[11] βP=βρβ2β2v
    加速度a[10] βa=β2v/β 动能Ek[11] βEk=βρβ3β2v
    下载: 导出CSV

    表  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
    下载: 导出CSV

    表  3  1006 钢的材料参数[24]

    Table  3.   Material parameters of 1006 steel[24]

    ρ/(g∙cm−3)E/GPaμA/MPaB/MPaCn˙ε0/s1
    7.892070.33502750.0220.361
    下载: 导出CSV

    表  4  理想弹塑性模型几何畸变比例因子

    Table  4.   Scaling factors of geometrically-distorted models of elastic-ideal plastic material

     βη βh βM βv βt βP βδ βε
    0.11.20.120.121.04660.09550.01310.11.0954
    0.11.50.150.151.10670.09040.01840.11.2247
    0.11.80.180.181.15830.08630.02410.11.3416
    0.12.00.200.201.18920.08410.02830.11.4142
    下载: 导出CSV

    表  5  考虑应变率效应和应变硬化效应几何畸变模型比例因子

    Table  5.   Scaling factors of geometrically-distorted models considering strain-rate sensitivity and strain hardening

     βη βh βM βv βt βP βδ βε
    0.11.20.120.121.07560.09300.01390.11.0954
    0.11.50.150.151.14420.08740.01960.11.2247
    0.11.80.180.181.20380.08310.02610.11.3416
    0.12.00.200.201.23960.08070.03060.11.4142
    下载: 导出CSV

    表  6  比例模型与原型的峰值位移和平均载荷相对误差

    Table  6.   Relative errors in the peak displacement and average force between the scale models and prototype

    模型δ/βδ)/mm相对误差/%P/βP/kN相对误差/%
    原型124.25871.95
    β=1/10, η=1.2121.6692.08473.432.057
    β=1/10, η=1.5120.6512.90369.483.433
    β=1/10, η=1.8118.6424.52075.855.420
    β=1/10, η=2.0117.6785.29572.250.417
    下载: 导出CSV

    表  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.74097.79
    β=1/10, η=1.289.6191.23598.420.644
    β=1/10, η=1.589.2221.67399.471.718
    β=1/10, η=1.886.2015.002101.063.344
    β=1/10, η=2.086.2324.968105.147.516
    下载: 导出CSV
  • [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)

  • 加载中
图(5) / 表(7)
计量
  • 文章访问数:  480
  • HTML全文浏览量:  227
  • PDF下载量:  102
  • 被引次数: 5
出版历程
  • 收稿日期:  2021-11-02
  • 修回日期:  2022-01-04
  • 网络出版日期:  2022-04-24
  • 刊出日期:  2022-05-27

目录

/

返回文章
返回