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

应力调整对石英玻璃珠低速冲击破碎行为的影响

宋一平 苗春贺 单俊芳 王鹏飞 徐松林

王维荣, 陈书杨, 王学军, 赵信文, 黄世璋, 李欣竹, 吴强. 一种适用于复杂加载状态下微喷射物质面密度测量的Asay膜方法[J]. 爆炸与冲击, 2024, 44(3): 034101. doi: 10.11883/bzycj-2023-0089
引用本文: 宋一平, 苗春贺, 单俊芳, 王鹏飞, 徐松林. 应力调整对石英玻璃珠低速冲击破碎行为的影响[J]. 爆炸与冲击, 2022, 42(7): 073103. doi: 10.11883/bzycj-2021-0244
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: SONG Yiping, MIAO Chunhe, SHAN Junfang, WANG Pengfei, XU Songlin. Effect of stress-state adjustment on fragmentation behavior of quartz glass beads subjected to low-velocity impact[J]. Explosion And Shock Waves, 2022, 42(7): 073103. doi: 10.11883/bzycj-2021-0244

应力调整对石英玻璃珠低速冲击破碎行为的影响

doi: 10.11883/bzycj-2021-0244
基金项目: 国家自然科学基金(11672286,11872361,11602267);中央高校基本科研业务费专项资金(WK2480000008);高压物理与地震科技联合实验室室开放基金(2019HPPES01)
详细信息
    作者简介:

    宋一平(1997- ),女,硕士,syiping@mail.ustc.edu.cn

    通讯作者:

    徐松林(1971- ),男,博士,研究员,博士生导师,slxu99@ustc.edu.cn

  • 中图分类号: O347

Effect of stress-state adjustment on fragmentation behavior of quartz glass beads subjected to low-velocity impact

  • 摘要: 结合高速摄影技术,应用SHPB加载装置,分别使用钢制、铝制和有机玻璃制3种透射杆,对直径约7.90、11.80、15.61 mm 3种尺寸的石英玻璃珠进行了低速冲击实验。根据不同透射杆条件下的玻璃珠破碎过程中的载荷-位移曲线,结合有限元软件计算玻璃珠在冲击作用下载荷的变化情况以及实验过程中玻璃珠的应变,探讨了应力调整对玻璃珠破碎过程的影响。结果表明:相同冲击条件作用下,改变透射杆的材料,会改变玻璃珠破碎过程中的载荷分布,即透射端边界波阻抗的改变会导致反射波发生改变,从而导致玻璃珠内部载荷发生变化;透射杆为铝材和有机玻璃材质时,玻璃珠在破碎过程中的载荷明显下降,在加载过程中伴随着垫块的变形,玻璃珠内部的应力调整时间变长;透射杆为钢杆时,玻璃珠的应变主要表现为两端最大,越靠近中间应变越小,对于透射杆为铝杆和有机玻璃杆的玻璃珠,透射端局部出现了卸载行为。采用有机玻璃透射杆之后,局部应力和变形降低的结果使得玻璃珠在经受较大的变形之后发生破碎,表明玻璃珠的破碎行为由局部变形和局部变形梯度共同控制。
  • 当材料表面机械加工等微缺陷经历冲击波作用后,会产生速度显著高于自由面的微米级颗粒,这就是所谓的微喷射现象[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.  Force-displacement curves and breakage patterns of quartz glass spheres

    图  2  改进的SHPB实验装置

    Figure  2.  Schematic diagram of a modified split Hopkinson pressure bar device

    图  3  直径7.90 mm玻璃珠的透射载荷-位移关系曲线

    Figure  3.  Transmitted load-displacement curves of glass sphere with diameter 7.90 mm

    图  4  直径7.90 mm玻璃珠的压缩过程

    Figure  4.  Compression processes of glass spheres with diameter 7.90 mm

    图  5  直径15.61 mm玻璃珠的透射载荷-位移关系曲线

    Figure  5.  Transmitted load-displacement curves of glass sphere with diameter 15.61 mm

    图  6  直径15.61 mm玻璃珠的压缩过程

    Figure  6.  Compression processes of glass spheres with diameter 15.61 mm

    图  7  玻璃珠的计算模型和有限元模型

    Figure  7.  Calculation model and finite element method model of glass sphere

    图  8  直径12.00 mm玻璃珠内部不同截面载荷-时间曲线

    Figure  8.  Load-time curves of different cross sections of glass beads with diameter 12.00 mm

    图  9  直径12.00 mm玻璃珠内部不同截面载荷-位置关系

    Figure  9.  Load-position relations of different cross sections of glass beads with diameter 12.00 mm

    图  10  直径12.00 mm玻璃珠内部剪应力分布

    Figure  10.  Shear stress distributions of glass beads with diameter 12.00 mm

    图  11  结合高速摄影计算的直径7.90 mm玻璃珠中应变演化

    Figure  11.  Strain evolutions in glass sphere with diameter 7.90 mm based on the high-speed photography

    图  12  基于高速摄影计算的直径15.61 mm玻璃珠中应变演化

    Figure  12.  Strain evolutions in glass spheres with diameter 15.61 mm based on the high-speed photography

    图  13  破碎阵面

    Figure  13.  Failure wave fronts

    图  14  基于剪切扩散理论计算的直径7.90 mm玻璃珠中应变演化

    Figure  14.  Strain evolutions in glass spheres with diameter 7.90 mm based on the shear activation diffusion theory

    图  15  基于剪切扩散理论计算的直径15.61 mm玻璃珠中应变演化

    Figure  15.  Strain evolution in glass sphere with diameter 15.61 mm based on shear activation diffusion theory

    表  1  玻璃珠两端载荷差统计

    Table  1.   Statistics of load differences between two ends of glass sphere

    直径/mm透射杆材料平均应变率/s−1两端载荷差/%
    7.907003.7~4.1
    5006.1~8.2
    有机玻璃45012.8~25.5
    11.806007.9~8.3
    450 8.0~15.8
    有机玻璃30062.9~77.8
    15.614505.1~9.1
    40048.3~54.1
    有机玻璃30074.3~82.7
    下载: 导出CSV

    表  2  材料参数

    Table  2.   Parameters of materials

    材料密度/(kg·m−3)声速/(m·s−1)弹性模量/GPa泊松比
    7 8005 100203.00.30
    2 7005 090 70.00.25
    有机玻璃1 8001 270 2.90.32
    石英玻璃珠2 5306 027 89.00.24
    下载: 导出CSV
  • [1] 徐松林, 单俊芳, 王鹏飞. 脆性材料高应变率压缩失效机制综述与研究进展 [J]. 现代应用物理, 2020, 11(3): 030101. DOI: 10.12061/j.issn.2095-6223.2020.030101.

    XU S L, SHAN J F, WANG P F. Review and research progress of dynamic failure mechanism for brittle materials under high strain rate [J]. Modern Applied Physics, 2020, 11(3): 030101. DOI: 10.12061/j.issn.2095-6223.2020.030101.
    [2] HUANG J Y, XU S L, HU S S. Influence of particle breakage on the dynamic compression responses of brittle granular materials [J]. Mechanics of Materials, 2014, 68: 15–28. DOI: 10.1016/j.mechmat.2013.08.002.
    [3] LIU C H, NAGEL S R, SCHECTER D A, et al. Force fluctuations in bead packs [J]. Science, 1995, 269(5223): 513–515. DOI: 10.1126/science.269.5223.513.
    [4] MAJMUDAR T S, BEHRINGER R P. Contact force measurements and stress-induced anisotropy in granular materials [J]. Nature, 2005, 435(7045): 1079–1082. DOI: 10.1038/nature03805.
    [5] GOLDENBERG C, GOLDHIRSCH I. Force chains, microelasticity, and macroelasticity [J]. Physical Review Letters, 2002, 89(8): 084302. DOI: 10.1103/physrevlett.89.084302.
    [6] HARTLEY R R, BEHRINGER R P. Logarithmic rate dependence of force networks in sheared granular materials [J]. Nature, 2003, 421(6926): 928–931. DOI: 10.1038/nature01394.
    [7] HUANG J Y, LU L, FAN D, et al. Heterogeneity in deformation of granular ceramics under dynamic loading [J]. Scripta Materialia, 2016, 111: 114–118. DOI: 10.1016/j.scriptamat.2015.08.028.
    [8] HUANG J Y, XU S L, YI H S, et al. Size effect on the compression breakage strengths of glass particles [J]. Powder Technology, 2014, 268: 86–94. DOI: 10.1016/j.powtec.2014.08.037.
    [9] 易洪昇, 徐松林, 单俊芳, 等. 不同加载速度下脆性颗粒的破坏特性 [J]. 爆炸与冲击, 2017, 37(5): 913–922. DOI: 10.11883/1001-1455(2017)05-0913-10.

    YI H S, XU S L, SHAN J F, et al. Fracture characteristics of brittle particles at different loading velocities [J]. Explosion and Shock Waves, 2017, 37(5): 913–922. DOI: 10.11883/1001-1455(2017)05-0913-10.
    [10] SHAN J F, XU S L, LIU Y G, et al. Dynamic breakage of glass sphere subjected to impact loading [J]. Powder Technology, 2018, 330: 317–329. DOI: 10.1016/j.powtec.2018.02.009.
    [11] MCDOWELL G R, AMON A. The application of Weibull statistics to the fracture of soil particles [J]. Soils and Foundations, 2000, 40(5): 133–141. DOI: 10.3208/sandf.40.5_133.
    [12] CHEONG Y S, SALMAN A D, HOUNSLOW M J. Effect of impact angle and velocity on the fragment size distribution of glass spheres [J]. Powder Technology, 2003, 138(2/3): 189–200. DOI: 10.1016/j.powtec.2003.09.010.
    [13] SALMAN A D, REYNOLDS G K, FU J S, et al. Descriptive classification of the impact failure modes of spherical particles [J]. Powder Technology, 2004, 143/144: 19–30. DOI: 10.1016/j.powtec.2004.04.005.
    [14] 方继松, 王珠, 熊迅, 等. 石英玻璃球撞击刚性壁的破碎过程 [J]. 高压物理学报, 2020, 34(1): 014101. DOI: 10.11858/gywlxb.20190764.

    FANG J S, WANG Z, XIONG X, et al. Fragmentation process of quartz glass spheres impacting rigid wall [J]. Chinese Journal of High Pressure Physics, 2020, 34(1): 014101. DOI: 10.11858/gywlxb.20190764.
    [15] POTAPOV A V, CAMPBELL C S. The two mechanisms of particle impact breakage and the velocity effect [J]. Powder Technology, 1997, 93(1): 13–21. DOI: 10.1016/S0032-5910(97)03242-7.
    [16] 简世豪, 苗春贺, 张磊, 等. 双石英玻璃珠的低速冲击破碎行为 [J]. 高压物理学报, 2021, 35(2): 024202. DOI: 10.11858/gywlxb.20200629.

    JIAN S H, MIAO C H, ZHANG L, et al. Fragmentation of double quartz glass spheres subjected to lower-velocity impact [J]. Chinese Journal of High Pressure Physics, 2021, 35(2): 024202. DOI: 10.11858/gywlxb.20200629.
    [17] SHIPWAY P H, HUTCHINGS I M. Fracture of brittle spheres under compression and impact loading. Ⅰ. Elastic stress distributions [J]. Philosophical Magazine A, 1993, 67(6): 1389–1404. DOI: 10.1080/01418619308225362.
    [18] CHAU K T, WEI X X, WONG R H C, et al. Fragmentation of brittle spheres under static and dynamic compressions: experiments and analyses [J]. Mechanics of Materials, 2000, 32(9): 543–554. DOI: 10.1016/S0167-6636(00)00026-0.
    [19] 黄俊宇. 冲击载荷下脆性颗粒材料多尺度变形破碎特性研究 [D]. 合肥: 中国科学技术大学, 2016: 82–84.

    HUANG J Y. Dynamic multiscale deformation behavior and particle-breakage properties of granular materials subjected to impact loading [D]. Hefei: University of Science and Technology of China, 2016: 82–84.
    [20] 苗春贺, 陈丽娜, 单俊芳, 等. 水泥砂浆抗弹性能研究 [J]. 高压物理学报, 2021, 35(2): 024205. DOI: 10.11858/gywlxb.20200609.

    MIAO C H, CHEN L N, SHAN J F, et al. Research on the ballistic performance of cement mortar [J]. Chinese Journal of High Pressure Physics, 2021, 35(2): 024205. DOI: 10.11858/gywlxb.20200609.
    [21] FENG R. Formation and propagation of failure in shocked glasses [J]. Journal of Applied Physics, 2000, 87(4): 1693–1700. DOI: 10.1063/1.372079.
    [22] JIANG H B, XU S L, SHAN J F, et al. Dynamic breakage of porous hexagonal boron nitride ceramics subjected to impact loading [J]. Powder Technology, 2019, 353: 359–371. DOI: 10.1016/j.powtec.2019.05.028.
  • 加载中
图(15) / 表(2)
计量
  • 文章访问数:  279
  • HTML全文浏览量:  128
  • PDF下载量:  44
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-06-22
  • 录用日期:  2022-06-01
  • 修回日期:  2021-08-20
  • 网络出版日期:  2022-06-07
  • 刊出日期:  2022-07-25

目录

/

返回文章
返回