• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

散体介质SHPB被动围压试验体应力计算的理论修正方法

陈昊祥 李杰 邓树新 王德荣 王明洋

王维荣, 陈书杨, 王学军, 赵信文, 黄世璋, 李欣竹, 吴强. 一种适用于复杂加载状态下微喷射物质面密度测量的Asay膜方法[J]. 爆炸与冲击, 2024, 44(3): 034101. doi: 10.11883/bzycj-2023-0089
引用本文: 陈昊祥, 李杰, 邓树新, 王德荣, 王明洋. 散体介质SHPB被动围压试验体应力计算的理论修正方法[J]. 爆炸与冲击, 2022, 42(6): 064901. doi: 10.11883/bzycj-2021-0357
WANG Weirong, CHEN Shuyang, WANG Xuejun, ZHAO Xinwen, HUANG Shizhang, LI Xinzhu, WU Qiang. An improved Asay foil method for measuring areal density of ejecta under complex loading conditions[J]. Explosion And Shock Waves, 2024, 44(3): 034101. doi: 10.11883/bzycj-2023-0089
Citation: CHEN Haoxiang, LI Jie, DENG Shuxin, WANG Derong, WANG Mingyang. A theoretically-modified method for calculating the volumetric stresses in passive confined pressure SHPB tests of granular materials[J]. Explosion And Shock Waves, 2022, 42(6): 064901. doi: 10.11883/bzycj-2021-0357

散体介质SHPB被动围压试验体应力计算的理论修正方法

doi: 10.11883/bzycj-2021-0357
基金项目: 北京市自然科学基金(8222010);河南省特种防护材料重点实验室开放课题(SZKFKT202102);北京建筑大学内涵发展-青年教师科研能力提升计划(X2102080921019)
详细信息
    作者简介:

    陈昊祥(1992- ),男,博士,chx@stu.bucea.edu.cn

    通讯作者:

    王明洋(1966- ),男,博士,教授,博士生导师,wmyrf@163.com

  • 中图分类号: O347; O344

A theoretically-modified method for calculating the volumetric stresses in passive confined pressure SHPB tests of granular materials

  • 摘要: SHPB被动围压试验为探究散体介质在爆炸和冲击荷载作用下的力学行为提供了一个行之有效的方法。针对相关试验设计和计算中存在的弊端和不足,借助经典板壳理论将SHPB被动围压试验中用于约束散体介质的刚性套筒简化为受均匀带状内压作用的圆柱形壳体。理论计算了套筒径向位移、环向应变与均匀带状内压及套筒几何、力学参数的关系,得到了套筒径向位移、环向应变沿其轴向的分布规律;分析了套筒长度、厚度、内外径以及均匀带状内压宽度之间等无量纲几何参数对计算结果的影响;将理论计算结果与试验和数值模拟结果进行对比,验证了理论计算结果的正确性。本文中提出的理论修正方法可为指导散体介质SHPB被动围压试验提供参考。
  • 当材料表面机械加工等微缺陷经历冲击波作用后,会产生速度显著高于自由面的微米级颗粒,这就是所谓的微喷射现象[1]。微喷射物质的质量、速度、空间分布会受到样品表面粗糙度、冲击波强度和加卸载历程、充气条件等一系列因素的影响[2-6],定量给出微喷射物质分布与影响因素之间的关系是研究微喷射形成机理的重要途径。基于动量守恒原理的Asay膜方法[7-8]在研究微喷射物质方面取得了巨大成功,被广泛应用,并被不断改进[9-10],相继出现了测量更高密度的Asay窗[8, 11-12]和尺寸更小的Asay膜[10]

    传统Asay膜方法测量微喷射物质的原理如图1所示,高速飞片撞击样品后在样品自由面产生微喷射,放置与样品自由面有确定间距的薄膜(Asay膜)用于收集样品自由面的微喷射物质,微喷射物质与薄膜碰撞后发生动量传递,并一起运动。采用光子多普勒速度计(photonic Doppler velocimetry, PDV)精确测量薄膜的运动过程,并按照如下4个基本假设,可获得微喷射物质的速度和质量分布。(1)瞬时喷射假设:当冲击波在自由面卸载时,全部微喷射物质是在足够短时间内产生的,即可近似认为微喷射是瞬时的,由此通过薄膜与样品自由面的距离及微喷射物质穿过该距离的时间,可计算碰撞到薄膜上的微喷射物质的速度。(2)微喷射物质分布均匀假设:即平板运动可以采用一维运动处理。(3)完全非弹性碰撞假设:微喷物质与薄膜碰撞后无反向溅射。(4)薄膜刚体运动假设:即不考虑冲击波在薄膜中的传播。

    图  1  Asay膜法微喷测量原理
    Figure  1.  The principle of the traditional Asay foil ejecta diagnostic

    基于上述假设,可计算不同时刻到达薄膜的微喷射粒子群的速度ve[9]

    ve=(d0+t0vf dt)/(d0+t0vf dt)tt
    (1)

    式中:d0为薄膜到样品的初始距离,t的时间零点为冲击波从自由面卸载的时刻,vf为薄膜的运动速度。

    再利用动量守恒,可计算任意时刻到达薄膜的微喷射物质的质量me

    me=mfvf0t dvfd0vf0t dvf
    (2)

    式中:mf为薄膜的质量。对各个时刻的me积分,则得到作用于薄膜上的微喷射物质总量分布。前述各质量都为单位面积质量,即面密度。

    随着对微喷射理论认识的深入,复杂加卸载条件下的微喷射现象被广泛关注,特别是样品经历连续2次冲击加载后的微喷射物质是目前较活跃的研究方向[4, 13-14]。在样品经历2次加载的情况下,自由面及微喷物质的速度随时间变化的典型状态如图2所示,首次加载后微喷射物质的分布和Asay膜常用场景是一致的。但是经历第2次加载后,首次加载后微喷粒子群速度ve1和二次加载后微喷粒子群速度ve2的值存在跳变,不再满足瞬时喷射假设和匀速运动假设,利用空腔距离和作用时间不能准确给出二次加载后某时刻微喷颗粒速度,即式(1)不成立,传统方法失效。因此,针对复杂加载微喷射物质测量,需要对Asay膜方法进行改进。

    图  2  典型二次加载微喷射速度示意图
    Figure  2.  Schematic of ejecta velocity under double shock loading

    为了解决前述复杂加载条件下传统Asay膜方法应用失效的问题,本文中,提出基于Asay膜探头相同位置PDV测速信号[15]给出微喷物速度分布信息的改进Asay膜方法(以下简称新方法),针对PDV探头对微喷物质穿透深度有限、微喷物速度上下限提取存在偏差等实际应用中的问题,采用数值实验对新方法的测量偏差及原因进行分析,通过轻气炮实验对新方法与传统Asay膜方法的测量效果进行对比评估。

    图3所示,t0时刻微喷射物质开始作用在Asay膜上,此时自由面与Asay膜的距离用l表示,h为膜的厚度,ˉve为某时刻微喷物质的平均速度;t为微喷物质作用于膜上的任意时刻,其中se为微喷射物质稠密主体区的运动距离,l为未作用于膜上的微喷物质的厚度,d为作用于膜上的微喷物质压缩后的厚度,sf为膜的运动距离,vf为膜的运动速度。

    图  3  微喷物质与Asay膜作用过程示意图
    Figure  3.  Schematic of ejecta impacting the Asay foil

    t时刻,假设有n个颗粒作用于Asay膜上,且质点间完全非弹性碰撞,根据动量守恒,微喷物质的速度和质量满足:

    {(me+mf)vf=ni=1(meivei)me=ni=1mei
    (3)

    测试过程中不可能追踪每个颗粒的速度和质量,假设t时刻微喷射粒子群的速度ve和质量me满足:

    ni=1(meivei)=meve
    (4)

    在单颗粒速度一致的情况下,式(4)严格成立。结合式(3)~(4)以及图3所示的运动关系:

    Δl=ll=se+dsf
    (5)

    并假设作用在Asay膜上的微喷射物质被压实,可以给出微喷物质的总量及分布满足关系:

    {me=mfvfvevfΔl=ve(tt0)+1ρ0mfvfvevft0vfdt
    (6)

    在实际测量过程中,mfvfve都可以通过相应测试方法给出,Δl无需初始空腔距离l即可求得,因此通过式(6)便可以获得微喷物质的总量及分布。

    为评估上述方法的数据处理不确定度,基于完全非弹性碰撞及瞬时喷射假定推导了Asay膜响应模型,如图4所示。为简化推导过程,假定单颗粒质量mp相同(可证明颗粒质量有分布,但不影响结果),忽略碰撞产生的动量矩。令喷射产生后t0时刻第1个颗粒碰膜,可获得如下的颗粒速度分布和膜速度的控制方程组分别为:

    图  4  微喷颗粒冲击Asay膜响应示意图
    Figure  4.  Schematics of ejecta particles impacting the Asay foil
    mpvp,i=(mp+mf)vf,i,tf,i=0,sf,i=0i=1
    (7)
    {vf,i={[(i1)mp+mf]vf,i1+mpvp,i}/(imp+mf)dti=(t0+ti1)(vp,i1vp,i)/(vp,ivf,i1)tf,i=tf,i1+dtisf,i=vf,i1dti+sf,i1i=2,3,,N
    (8)

    式中:vp,i为第i个颗粒的速度,vf,i为第i个颗粒碰膜后的膜速,tf,i为第i个颗粒碰膜的时间,sf,i为第i个颗粒碰膜后膜的运动距离。

    通过生成有随机速度分布的颗粒群,基于每个颗粒都符合匀速运动的假设,则可以获得t0时刻颗粒空间速度位置分布。将其代入式(4),则可以获得vf-t曲线,ve-t,即Asay膜数值实验数据。

    为验证上述新Asay膜方法在实际应用中的效果,设计了对比实验,如图5所示,气体炮发射搭载飞片的弹丸,撞击样品产生微喷射。相较于传统的Asay膜方法,新方法增加了直接测量微喷物速度的PDV探头。由此可以获得的实验数据包括Asay膜测速结果和PDV测微喷射速度结果,结合样品与Asay膜探头的初始距离及Asay膜的面密度,即可分别采用传统Asay膜方法和本文中提出的新Asay膜方法对微喷物的面密度进行分析,对比评估新方法的测量效果。

    图  5  基于PDV测速的Asay膜实验示意图
    Figure  5.  Schematic of the improved Asay foil experiment based on PDV velocity measurement

    采用Matlab随机函数,生成了3种具有典型特征的颗粒群速度分布,颗粒群运动(速度方向为水平向右)t0时间后,其分布如图6所示,其中指数分布(exponential distribution, ED)和反指数分布(anti-exponential distribution, AED)都是利用Matlab的指数分布随机函数生成的,线性分布(linear distribution,LD)则是利用Matlab的均匀分布随机函数生成的。ED分布中颗粒主要集中在低速区,使用的生成函数为:vp,i=2350+random(‘exp’, 2400, 1, N)/50;而AED分布中颗粒集中在高速区,使用的生成函数为:vp,i=2650+random(‘exp’, 2400, 1, N)/50。LD分布中每种速度的颗粒数量基本一致,使用的生成函数为:vpi=2350+300×rand (N, 1)。以上生成函数中‘exp’为生成方法,N为颗粒总数,在生成过程中N=10 000,总面密度为100 mg/cm2,单颗粒等效面密度为0.01 mg/cm2。本文中讨论的碰膜时刻喷射粒子速度分布满足高速粒子在前、低速粒子在后这一基本特征,忽略了在复杂加载后较短时间内粒子速度分布不满足该特征的特殊情形,而常见应用场景中都是在距离复杂加载较长时间后进行测量的,因此本文中暂不对这一特殊情形进行讨论。

    图  6  数值实验获得的微喷射分布
    Figure  6.  Ejecta particles distributions obtained from numerical experiments

    当设定膜的面密度为160 mg/cm2时,利用控制方程(8)获得的膜速曲线vf-t图7(a)所示。图7(b)~(c)则分别展示了理论ve-t、Me-t曲线,Me代表各个时刻作用于膜上的微喷物累积面密度。

    图  7  数值实验数据
    Figure  7.  Numerical experiment data

    精确反演面密度分布时,准确的vf-t曲线和ve-t曲线缺一不可,而在实际复杂加载实验条件下,几乎不可能获得精确的ve-t曲线。以下将讨论使用近似ve-t曲线的处理方法。

    (1)定值法:在实验测量中,直接测量微喷物速度的PDV探头至少可以给出微喷物速度分布中的某个速度值ˉve,使用ˉve-t进行数据处理。

    (2)线性插值法:假定已知速度上限(t1, ve1)和下限(t2, ve2),利用两点线性插值给出ve-t曲线。

    分别采用ve取定值(2350、2450、2550、2650 m/s)及取速度上限(t1, 2650 m/s)和下限(t2, 2350 m/s),结合3种分布下的t1、t2值,线性插值给出ve-t曲线。利用式(6)对数值实验数据处理分析结果的Me-t曲线及数据处理获得的面密度值Me与理论值Met的相对偏差如图810所示,其中Met通过累积数值实验中每一时刻碰膜的粒子质量精确给出。

    图  8  指数分布条件下Asay膜数值实验数据处理结果及相对偏差分析
    Figure  8.  Numerical experiment data processing results and deviation analysis using the new Asay foil method for exponential distribution of ejecta
    图  9  反指数分布条件下Asay膜数值实验数据处理结果及相对偏差分析
    Figure  9.  Numerical experiment data processing results and deviation analysis using the new Asay foil method for anti-exponential distribution of ejecta
    图  10  线性分布条件下Asay膜数值实验数据处理结果及偏差分析
    Figure  10.  Numerical experiment data processing results and deviation analysisusing the new Asay foil method for linear distribution of ejecta

    Asay膜方法常用于较大面密度(10 mg/cm2以上)微喷物质定量测量,以下将针对实际使用测量范围对新方法的测量误差进行评估分析。

    当微喷颗粒呈指数分布时,除ve取值为2650 m/s时最大偏差达到15%外,其余3种取值情况下在5~100 mg/cm2面密度范围内的偏差都小于10%;当微喷颗粒呈反指数分布时,当ve取值为2350和2450 m/s时,偏差大于10.0%,分别为18.5%和11%,其余2种取值时在5~100 mg/cm2面密度范围内的偏差都小于10%;当微喷颗粒呈线性分布时,除ve取值为2350 m/s时最大偏差达到14%外,其余3种取值情况下在5~100 mg/cm2面密度范围内的偏差都小于10%。

    使用线性插值方法对3种微喷物质分布的数值实验结果进行数据分析都能取得较好的效果,即在5~70 mg/cm2面密度范围内的偏差不超过5%,在5~100 mg/cm2面密度范围内的偏差不超过10%。

    由上述分析可见,在可以给出微喷物质速度分布的上限和下限的情况下,无论微喷物质分布如何,采用PDV测速结合Asay膜速度进行数据分析,可以获得较高精度的微喷物质分布和累积面密度结果。然而,实验条件下,PDV测速探头可以给出微喷物质速度上限(t1, ve1),但由于PDV测速探头使用的激光对微喷物穿透深度存在极限,通常无法给出真实的速度下限。

    针对PDV信号无法给出微喷物质速度下限的情况,可以先根据定值法给出初步的Me-t分布曲线,然后根据PDV探头穿透深度在Me-t曲线上给出极限穿深面密度Me, l时的t2,并取此时PDV测速信号对应的速度下限为ve2。根据PDV测速探针使用的1550 nm激光在典型微喷射物质中穿透深度约20 mg/cm2,取Me,l为10、20和40 mg/cm2,利用数值实验结果对考虑实验条件下的PDV探头穿透深度进行测量偏差分析,分析使用的(t1, ve1)、(t2, ve2)如表1所示。

    表  1  不同数值实验条件下的数据处理参数
    Table  1.  Data processing parameters under different numerical experiment conditions
    面密度
    分布情况
    t1/μst2/μsve1/(m·s−1)ve2/(m·s−1)
    Me,l=10 mg/cm2Me,l=20 mg/cm2Me,l=40 mg/cm2Me,l=10 mg/cm2Me,l=20 mg/cm2Me,l=40 mg/cm2
    指数02.9503.3553.7852772.702460.832426.112393.23
    反指数00.0460.1010.2402649.992644.862639.212625.85
    线性00.2680.5871.2512649.992620.402588.022528.43
    下载: 导出CSV 
    | 显示表格

    图11给出了假定PDV探头穿透深度分别为10、20、40 mg/cm2时数据分析结果与理论值的偏差,3种微喷射物质分布情况下,随着穿透深度的增大,数据分析结果与理论值的偏差都呈现出减小趋势,且数据分析结果较理论值都偏大;经分析,如图7(b)所示,ve-t曲线是下凹曲线,采用曲线靠前部分的两点进行线性插值后,用于计算碰膜微喷射质量的ve将被低估,且低估值越来越大,由式(6)可知,ve偏小将导致me偏大。

    图  11  模拟实际测试条件Asay膜数值实验结果与理论偏差分析
    Figure  11.  Analysis of numerical experimental results and theoretical deviation of Asay foil under simulated actual test conditions

    当微喷物质分布集中靠近Asay膜方向(ED)或样品方向(AED)时,采用线性插值法数据处理的面密度结果在5~100 mg/cm2范围内与理论值偏差小于10%;而当微喷射物质为均匀分布(LD)时,采用线性插值方法数据处理的面密度结果在5~100 mg/cm2范围内的与理论偏差小于20%。针对LD分布时偏差较其余2种分布大的现象进一步分析,通过图7(b)可知,处于插值点附近区域的ve与理论值差较小,而越远离插值区域偏差越大,ED和AED分布中主要微喷射物质与插值区域更靠近,而LD分布中各个区域微喷物都一样,因此,新的Asay膜方法对均匀分布的微喷射累积面密度测量误差较大。

    2.1节中通过数值实验分析了新Asay膜方法进行微喷物面密度测量的偏差及其原因,为检验新方法在实验中的应用效果,利用轻气炮开展了验证实验。采用飞片冲击表面预置微米粉末样品开展实验,如图5所示,飞片为尺寸13.2 mm×2 mm的304不锈钢,飞片速度为1923.1 m/s,样品为20 mm×2 mm不锈钢片表面均匀预置粒径分布为3~6 μm的球形锡颗粒,粉末初始实测面密度约44 mg/cm2。PDV探头中心和Asay膜探头中心位于同一圆周上。测得的Asay膜速度及相应的PDV测速数据如图12所示。

    图  12  实验Asay膜测速频谱信号和PDV测速频谱信号
    Figure  12.  Asay foil velocity spectrum signal and PDV velocity spectrum signal in experiment

    结合测试结果及已知的初始条件,分别利用传统的空腔法(传统Asay膜法)和本文中提出的新方法(新Asay膜法)对实验数据进行了处理,相应的处理参数如表2所示。其中新方法中速度上限(t1, ve1)和下限(t2, ve2)的选取分别为:膜有速度响应时刻为t1,膜速陡升时刻为t2t1时刻对应的频谱速度上限为ve1t2时刻对应的速度频谱下限为ve2

    表  2  两种Asay膜数据处理方法使用的参数
    Table  2.  Parameters used in two different Asay foil data processing methods
    方法 l/mm t0/μs mf/(mg·cm−2) t1/μs t2/μs ve1/(m·s−1) ve2/(m·s−1)
    传统Asay膜法 19.5 28.33 156
    新Asay膜法1 156 34.05 39.33 4180 2026
    新Asay膜法2 156 34.05 39.33 4380 1826
    新Asay膜法3 156 34.05 39.33 3980 2226
    下载: 导出CSV 
    | 显示表格

    表2中新方法选取了3组数据,其中新Asay膜法1是按图12(b)中PDV频谱对应的t1t2时刻分别提取速度上限和下限。考虑到在频谱上提取速度时引入的误差,新Asay膜法2和新Asay膜法3中考虑了2种极端情况,即频谱数据提取速度上下限存在±200 m/s的误差。图13(a)对比了分别采用传统Asay膜方法和本文中提出的新Asay膜方法的数据处理结果,在考虑PDV频谱速度上下限误差的情况下,2种方法给出的微喷物质分布及累积面密度有较好的一致性。需特别说明的是,Asay膜方法测得的面密度比预置粉末的面密度小的主要原因是粉末未完全从304钢样品表面脱落。图13(b)分析了2种方法的测量结果的相对偏差,采用新方法测得的累积面密度在初始阶段较传统方法测得的小,但随着累积面密度的增大,新Asay膜方法与传统Asay膜方法的差异逐渐减小,并且新Asay膜法1和新Asay膜法2得到的累积面密度较传统方法的大。综合来看,在累积面密度大于5 mg/cm2的情况下,在考虑PDV穿透深度、微喷物速度上下限提取偏差等因素后,使用新Asay膜方法得到的面密度空间分布和累积面密度相对传统方法的偏差都小于20%。

    图  13  两种Asay膜数据处理方法实验结果及相对偏差
    Figure  13.  Comparison of experimental results and relative deviation of two Asay foil data processing methods

    针对复杂加载条件下微喷射物质测量这一具体问题,分析了传统Asay膜方法存在的不足,特别是在某些多次加载条件下无法给出传统Asay膜数据处理方法所必需的初始空腔距离及自由面起动时间等参数,提出了基于PDV测速结合Asay膜测速结果进行微喷射物质质量分布分析的新方法。采用传统方法和新方法对数值实验及气炮加载实验数据进行对比分析,验证了新方法的有效性,并分析了测量偏差,获得的主要认识如下。

    (1)基于PDV测速结合Asay膜测速处理微喷射物质质量分布在无法获得自由面起跳时间及膜与自由面初始空腔距离时有效,扩充了Asay膜方法的适用范围;新Asay膜方法的测量值在累积面密度大于5 mg/cm2时,相对传统Asay膜方法测量偏差小于20%,即与传统Asay膜方法测试精度基本相当。

    (2)在考虑PDV探头对微喷射物质穿透深度有限的条件下,新Asay膜方法对微喷射物质空间分布具有一定的敏感性,当微喷射物质集中于高速区或低速区时,微喷物累积面密度测试结果与理论值有较小的偏差;当微喷射物质均匀分布时,结果与参考值偏差较大;对于已知被测对象大致分布的情况,新Asay膜方法通过改变对ve的拟合方式,理论上测试精度仍可提升。

    (3)微喷射速度测试PDV探头要求与Asay膜探头的被测对象具有相同物理状态或处于相同位置,因此,当微喷射状态空间分布极不均匀时,如果无法保证2种探头的被测对象基本一致,该方法的测试误差可能被放大,应当慎用。

    本文中,仅对针对激光干涉测速信号和Asay膜信号的本身特点开展了分析,分析中不包括Asay膜探头自身的测量偏差、实际物理过程偏离基本假设等。针对此类因素,后续将开展进一步研究。本文中讨论的喷射粒子碰膜时都具有高速粒子先碰膜、低速粒子后碰膜的特点,针对第2次加载强度较低或间隔时间较长时,理论上存在高速粒子先碰膜、低速粒子后碰膜这一特殊情况,后续将进一步开展数值实验分析,拓宽本文方法的应用场景。

  • 图  1  钢套筒端部效应的示意图

    Figure  1.  The configuration for the edge effect in the steel sleeve

    图  2  套筒受力

    Figure  2.  Stress on the steel sleeve

    图  3  套筒环向应变沿轴线的分布

    Figure  3.  Distribution of hoop strain on the outer surface along with the steel sleeve

    图  4  理论计算得到的修正系数与数值模拟结果的比较

    Figure  4.  Comparison of theoretical correction factors with numerical results

    图  5  模型网格划分

    Figure  5.  The element meshes for the numerical models

    图  6  理论计算得到的修正系数与数值模拟结果的比较(算例1)

    Figure  6.  Comparison of theoretical correction factors with numerical results (case 1)

    图  7  理论计算得到的修正系数与数值模拟结果的比较(算例2)

    Figure  7.  Comparison of theoretical correction factors with numerical results (case 2)

    图  8  理论轴向应力-应变理论曲线与试验结果[4]的比较

    Figure  8.  Comparison of theoretical axial stress-stain curves with experimental results[4]

  • [1] YAMAMURO J A, ABRANTES A E, LADE P V. Effect of strain rate on the stress-strain behavior of sand [J]. Journal of Geotechnical and Geoenvironmental Engineering, 2011, 137(12): 1169–1178. DOI: 10.1061/(ASCE)GT.1943-5606.0000542.
    [2] LUO H, LU H, COOPER W L, et al. Effect of mass density on the compressive behavior of dry sand under confinement at high strain rates [J]. Experimental Mechanics, 2011, 51(9): 1499–1510. DOI: 10.1007/s11340-011-9475-2.
    [3] LUO H Y, COOPER W L, LU H B. Effects of particle size and moisture on the compressive behavior of dense Eglin sand under confinement at high strain rates [J]. International Journal of Impact Engineering, 2014, 65: 40–55. DOI: 10.1016/j.ijimpeng.2013.11.001.
    [4] 魏久淇, 张春晓, 曹少华, 等. 一种散体材料SHPB被动围压试验体应力修正方法 [J]. 爆炸与冲击, 2020, 40(12): 124201. DOI: 10.11883/bzycj-2019-0411.

    WEI J Q, ZHANG C X, CAO S H, et al. A volume stress correction method for SHPB passive confined pressure of granular materials [J]. Explosion and Shock Waves, 2020, 40(12): 124201. DOI: 10.11883/bzycj-2019-0411.
    [5] 赵章泳, 邱艳宇, 紫民, 等. 含水率对非饱和钙质砂动力特性影响的试验研究 [J]. 爆炸与冲击, 2020, 40(2): 023102. DOI: 10.11883/bzycj-2019-0066.

    ZHAO Z Y, QIU Y Y, ZI M, et al. Experimental study on dynamic compression of unsaturated calcareous sand [J]. Explosion and Shock Waves, 2020, 40(2): 023102. DOI: 10.11883/bzycj-2019-0066.
    [6] 文祝, 邱艳宇, 紫民, 等. 钙质砂的准一维应变压缩试验研究 [J]. 爆炸与冲击, 2019, 39(3): 033101. DOI: 10.11883/bzycj-2018-0015.

    WEN Z, QIU Y Y, ZI M, et al. Experimental study on quasi-one-dimensional strain compression of calcareous sand [J]. Explosion and Shock Waves, 2019, 39(3): 033101. DOI: 10.11883/bzycj-2018-0015.
    [7] HOPKINSON B. A method of measuring the pressure produced in the detonation of high explosives or by the impact of bullets [J]. Proceedings of the Royal Society of London: Series A, 1914, 89(612): 411–413. DOI: 10.1098/rspa.1914.0008.
    [8] KOLSKY H. An investigation of the mechanical properties of materials at very high rates of loading [J]. Proceedings of the Physical Society: Section B, 1949, 62(11): 676–700. DOI: 10.1088/0370-1301/62/11/302.
    [9] BRAGOV A M, LOMUNOV A K, SERGEICHEV I V, et al. Determination of physicomechanical properties of soft soils from medium to high strain rates [J]. International Journal of Impact Engineering, 2008, 35(9): 967–976. DOI: 10.1016/j.ijimpeng.2007.07.004.
    [10] MA Z, RAVI-CHANDAR K. Confined compression: a stable homogeneous deformation for constitutive characterization [J]. Experimental Mechanics, 2000, 40(1): 38–45. DOI: 10.1007/BF02327546.
    [11] SONG B, CHEN W N, LUK V. Impact compressive response of dry sand [J]. Mechanics of Materials, 2009, 41(6): 777–785. DOI: 10.1016/j.mechmat.2009.01.003.
    [12] 曹志远. 板壳振动理论 [M]. 北京: 中国铁道出版社, 1989: 304–305.
    [13] 王敏中, 王炜, 武际可. 弹性力学教程 [M]. 北京: 北京大学出版社, 2011: 264–266.
  • 加载中
图(8)
计量
  • 文章访问数:  386
  • HTML全文浏览量:  224
  • PDF下载量:  55
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-08-23
  • 修回日期:  2021-12-31
  • 网络出版日期:  2022-04-15
  • 刊出日期:  2022-06-24

目录

/

返回文章
返回