Loading [MathJax]/jax/output/SVG/jax.js

无机玻璃动态压缩破坏的离散元模拟

马棋棋 熊迅 郑宇轩 周风华

梁民族, 李翔宇, 覃金贵, 卢芳云, 王马法. PELE弹丸壳体的损伤断裂模型[J]. 高压物理学报, 2014, 28(6): 709-715. doi: 10.11858/gywlxb.2014.06.011
引用本文: 马棋棋, 熊迅, 郑宇轩, 周风华. 无机玻璃动态压缩破坏的离散元模拟[J]. 高压物理学报, 2019, 33(4): 044101. doi: 10.11858/gywlxb.20190719
LIANG Min-Zu, LI Xiang-Yu, QIN Jin-Gui, LU Fang-Yun, WANG Ma-Fa. Damage and Fracture Model for Shell of PELE[J]. Chinese Journal of High Pressure Physics, 2014, 28(6): 709-715. doi: 10.11858/gywlxb.2014.06.011
Citation: MA Qiqi, XIONG Xun, ZHENG Yuxuan, ZHOU Fenghua. Discrete Element Simulations of Dynamic Compression Failure of Inorganic Glass in SHPB Tests[J]. Chinese Journal of High Pressure Physics, 2019, 33(4): 044101. doi: 10.11858/gywlxb.20190719

无机玻璃动态压缩破坏的离散元模拟

doi: 10.11858/gywlxb.20190719
基金项目: 国家自然科学基金(11390361)
详细信息
    作者简介:

    马棋棋(1993-),男,硕士研究生,主要从事冲击动力学研究. E-mail:737488513@qq.com

    通讯作者:

    郑宇轩(1986-),男,博士,副教授,主要从事冲击动力学研究. E-mail:zhengyuxuan@nbu.edu.cn

  • 中图分类号: O347

Discrete Element Simulations of Dynamic Compression Failure of Inorganic Glass in SHPB Tests

  • 摘要: 利用离散元软件PFC2D(Particle Flow Code)建立了分离式霍普金森压杆(SHPB)系统,模拟了无机玻璃圆柱和圆盘试件在冲击压缩下的动态力学行为和失效破坏模式。结果表明:无机玻璃作为典型的脆性材料,其抗压强度具有明显的应变率效应,而杨氏模量则对应变率不敏感;无机玻璃圆柱的破坏过程受纵向压力、端面摩擦力以及横向惯性力的影响,初期微裂纹呈三角状分布,随着纵向应力水平的提高,出现明显的泊松效应,产生横向张应力,致使微裂纹沿纵向扩展,最终试件发生沿轴向的劈裂断裂;摩擦系数和泊松比对试件破坏模式及强度有一定影响。将建立的SHPB数值实验平台用于模拟无机玻璃巴西圆盘试验,揭示了圆盘发生中心开裂的拉伸特征及拉伸强度的应变率相关性。

     

  • 自1996年German-French Research Institute Saint Louis (ISL)实验室第一次提出横向效应增强型侵彻体(Penetrator with Enhanced Lateral Effect, PELE)概念以来, 国内外学者对PELE弹丸的作用过程进行了广泛研究。在PELE弹丸壳体断裂的机理方面, Kesberg等人[1]认为壳体侵彻目标的内核因密度较低而被压缩并沿径向膨胀, 导致壳体内壁压力剧增至吉帕级, 弹丸穿透目标后能量瞬间释放, 导致壳体破裂成大量破片。Paulus等人[2]结合Recht-Ipson模型[3]和Mott模型[4]计算出弹体的剩余速度、装填物的压力以及壳体破片的径向飞散速度, 并且用弱激波理论描述了PELE弹丸壳体的断裂过程。对于PELE弹丸壳体断裂的研究, 除了可从力学和能量的角度定性地解释PELE弹丸壳体断裂过程外, 还可以从损伤断裂的角度, 具体分析其在侵彻过程中的断裂机理。

    在研究材料动态拉伸断裂方面有许多模型, Seaman、Curran和Shockey[5]从实验出发, 基于微损伤的统计性分析结果, 提出了成核与长大模型(Nucleation and Growth, NAG)。但是NAG模型计算需要大量的实验参数, 实用性比较差。封加波等人[6-7]在NAG模型的基础上利用单孔洞近似, 从能量平衡的角度出发, 提出了简单实用的损伤度函数模型。本研究基于损伤度模型, 结合PELE弹丸侵彻过程, 将PELE弹丸壳体的膨胀断裂过程分为加速膨胀和匀速膨胀两个阶段, 推导出这两个阶段PELE弹丸壳体的损伤断裂方程, 并根据损伤断裂方程分析影响PELE弹丸横向效应的内在机理。

    PELE弹丸壳体一般为高密度金属或合金。金属或合金材料在高应变率下的变形过程比较复杂, 理论上通常把动态变形过程划分为几个阶段, 如弹性阶段、塑性阶段和失效阶段, 对不同变形阶段需要用不同的本构关系表达式进行描述。动态塑性本构方程可以写为

    σ=φ(ep,˙ep) (1)

    式中:σ为材料的动态屈服应力, ep为径向塑性自然应变, 为径向塑性自然应变率。

    根据文献[7], 金属壳体的动态本构关系可以表示为

    σ=σ0+Kep+η˙ep (2)

    式中:σ0为静态屈服强度, K为材料强化系数, η为材料黏度系数。

    对于柱壳结构有

    ep=ln(R(t)R0) (3)
    ˙ep=u(t)/R(t) (4)
    εp=R(t)R0R0 (5)

    式中:R0为壳体初始半径, R(t)为壳体实时半径, R0R(t)均是壳体内径和外径的平均值, εp为径向塑性应变, u(t)为壳体径向膨胀速度。结合(3)式和(5)式可以得到epεp的关系

    ep=ln(1+εp) (6)

    将(6)式代入(2)式, 可以将动态本构关系改写成以下形式

    σ=σ0+Kln(1+εP)+η˙ep (7)

    本研究所采用的断裂准则建立在封加波提出的损伤度模型的基础上。金属材料的损伤度包括在应力波作用下材料初始损伤缺陷及其发展, 以及新产生的损伤缺陷及其发展, 而且材料的损伤度具有累积效应。在PELE弹丸壳体动态断裂的过程中, 为了考虑材料的损伤累积效应, 引入损伤变量D(t), D(t)定义为材料某一截面上微孔洞所占的面积比。参照文献[6], 有如下形式的方程

    dDdt=[(1D)(σσ0)]2ηK (8)

    求解方程(8)式需要知道σ(或εp随时间变化的规律。从(4)式和(5)式看出, 需要知道壳体膨胀的实时半径R(t)和壳体径向膨胀速度u(t)的函数形式。在处理实际问题时, 很难求出解析表达式, 一般需要结合实验结果求解。

    PELE弹丸是一种基于物理效应、不需要装填炸药和引信的新概念弹药, 主要用于防空、反导、城市作战等领域。PELE弹丸由壳体和内核组成, 壳体一般为密度较大的金属材料, 内核是密度及强度比较低的金属或聚合物。高密度金属壳体对靶板具有良好的侵彻性能, 而内核材料在侵彻过程中被封闭压缩, 巨大的压力在内核中迅速积聚, 弹丸穿透靶板后能量瞬间释放导致外壳破裂解体。

    PELE的横向增强效应作用过程分为3个主要时段和一个后效作用时段[1]。如图 1所示, 第1时段, 当弹丸击中目标后, 外层高密度壳体侵彻目标, 内核因材料强度低而无法侵彻, 将会在靶前停滞而被压缩; 第2时段, 封闭压缩过程中产生的高压将使内核材料沿径向膨胀, 导致壳体内壁压力急剧增加, 可达到数吉帕; 第3时段, 当应力应变达到壳体的破坏极限时, 出靶后的卸载应力波致使外层弹体破裂成大小不同的破片。后效作用时段则是壳体破片的径向飞散以及内核和塞块间继续相互作用的过程。

    图  1  PELE的3个阶段
    Figure  1.  Three stages of penetration process of the PELE

    图 2给出了PELE弹丸计算模型, 其中壳体材料为35CrMnSi, 长度为105 mm, 外径为30 mm, 内径为18 mm, 弹底部厚为15 mm; 内核材料为聚乙烯, 长度为90 mm, 直径为18 mm; 靶板材料为装甲钢, 厚度为8 mm。网格平均尺寸为壳体1 mm, 内核1.2 mm, 靶板采用变步长网格(0.5~10 mm)以提高计算效率。为了获得侵彻过程中的时间历程参数, 在弹丸模型的壳体中心线沿轴向均匀分布4个观测点(Gauge点), 每个观测点间隔20 mm, 如图 2所示。

    图  2  PELE弹丸计算模型
    Figure  2.  Calculation model of the PELE

    表 1表 2分别为材料模型和材料参数。数值计算采用拉格朗日方法, 材料均采用侵蚀算法, 靶板采用材料应变失效模型, 壳体使用有效应力失效模型。采用Stochastic模型[8]模拟壳体随机断裂, 其中材料弱化点的破坏概率在[0, 1]范围为

    表  1  材料模型
    Table  1.  Material model
    Component Material species Equation of state Strength model Invalidation model
    Shell 35CrMnSi Linear Johnson-Cook Principal stress
    Core Polyethylene Shock von-Mises -
    Target plate RHA Shock von-Mises Material strain
    下载: 导出CSV 
    | 显示表格
    表  2  材料参数
    Table  2.  Material parameters
    Material ρ/(g/cm3) E/(GPa) σ/(MPa) G/(GPa) ν
    35CrMnSi 7.8 210 1 270 81 0.33
    Polyethylene 0.94 1.5 0.26 0.55 0.45
    RHA 7.8 171 1400 64.1 0.33
    下载: 导出CSV 
    | 显示表格
    P=1exp(ceγε/γ) (9)

    式中:P为应变为ε时的网格破坏概率, cγ取决于材料性质, 根据文献[8], 取c=0.467, γ=10。

    PELE弹丸以700 m/s的着靶速度垂直侵彻8 mm厚的装甲钢靶板。壳体中Gauge点沿壳体径向速度变化曲线如图 3所示。从图 3中看出, PELE侵彻过程第2时段, 壳体径向先经历短暂加速膨胀阶段, 膨胀速度达到一个较大值; PELE侵彻过程第3时段, 壳体膨胀速度小幅震荡并趋于稳定。

    图  3  PELE弹丸壳体Gauge点的径向速度
    Figure  3.  The radical velocity of gauge points in the PELE

    PELE实验弹丸参数:外径为30 mm, 内径为18 mm, 长度为105 mm; 壳体材料采用35CrMnSi合金钢, 内核为聚乙烯, 长度为90 mm, 靶板为装甲钢, 厚度为8 mm; 装配全弹质量为350 g。

    实验弹采用轻气炮发射。弹丸着靶和出靶分别设置了测速点探针, 以测量着靶速度和弹丸穿透靶板的时间。为了获得弹丸头部出靶时的膨胀速度, 设计了一套环形电探针测量装置。实验装置示意图如图 4所示, 实验装置各部分如图 5所示。

    图  4  实验装置示意图
    Figure  4.  Experimental setup
    图  5  实验装置
    Figure  5.  Experimental device

    表 3分别统计了3发有效实验的弹丸着靶速度、弹丸头部从着靶到出靶的穿靶时间, 以及出靶时弹丸头部的径向速度。应该指出的是, 由于弹丸穿靶时间极短, 弹丸出靶时头部的径向速度比较大, 探针测试存在一定的误差。

    表  3  着靶速度与穿靶时间统计
    Table  3.  Statistics of the impact velocity and the penetrating time
    Experimental No. Velocity/(m/s) Time/(μs) Radical velocity/(m/s)
    1 643 8.4 127
    2 694 7.2 142
    3 762 6.5 153
    下载: 导出CSV 
    | 显示表格

    壳体膨胀断裂主要处于PELE侵彻过程的第2时段和第3时段。壳体动态断裂过程可以描述为:应力波从内核传入壳体内壁, 壳壁开始加速膨胀, 在应力波多次反射效应下膨胀速度振荡衰减, 直至最终趋于平稳。其衰减规律基本上按余弦函数变化[9-10], 所以径向速度公式可以写为

    u(t)={Ctt0t<t1Aeλtcos(ψt)+u0t1t<tc (10)

    式中:AψCλ为与总体模型有关的参数, u0为稳态时的径向速度值。t0为应力波传到壳体外壁的时刻, t1为壳体膨胀加速结束的时刻, tc为壳体发生贯穿性断裂的时刻。

    根据侵彻过程壳体的膨胀规律, 可以将PELE壳体断裂过程分为壳体加速膨胀阶段和壳体匀速膨胀阶段。调整数值模拟中弹丸的着靶速度, 使其工况与3发实验相吻合, 对比数值模拟结果与实验结果发现, PELE壳体加速膨胀的时间与壳体从着靶到出靶的时间相吻合。可以看出当壳体处于靶板内时, 壳体处于加速状态; 当壳体出靶后, 壳体匀速膨胀。

    根据PELE侵彻过程壳体的膨胀规律, 将PELE壳体断裂过程分为两个阶段:壳体加速膨胀阶段和壳体匀速膨胀阶段。

    第1阶段:壳体加速膨胀阶段。假定自然应变率可以近似看作常数

    ˙ep=u(t)/R(t)=c (11)

    另外

    dt=dep/˙ep (12)

    对(12)式积分得

    t1t0=ep1/˙ep (13)

    将(7)代入(8), 在[t0, t1]内积分得

    D1=13ω1(ω1+ep1)3+3ω1ω31 (14)

    式中:D1t1时刻的损伤度。(14)式中

    ω1=(K/η)˙ep (15)

    通过实验结果可以得到代入(13)式可以求出代入(14)式求出D1

    第2阶段:壳体匀速膨胀阶段。第1阶段结束时刻即为第2阶段初始时刻。依据PELE壳体膨胀规律, 给出描述壳体膨胀速度的变化规律为

    u(t)=Aeλtcos(ψt)+u0 (16)

    (16) 式反映应力波多次反射引起的衰减振荡, 代表了壳体膨胀速度的真实情况, 但直接代入损伤度方程难以求解。为了方便积分, 这里按照能量等效原理取均值, 假设柱壳膨胀速度u(t)等于常数

    u(t)=u0 (17)

    由于

    ˙ep=u(t)/R(t) (18)
    ˙ep2=u(t)/R2 (19)
    u(t)=dR(t)/dt (20)

    将(5)式、(18)式、(19)式代入(20)式, 得

    dt=dR(t)/u(t)=dεp/˙ep2 (21)

    式中: 为第2阶段的初始径向塑性自然应变率, R2为第2阶段壳体的初始半径。

    将(7)式、(17)式、(21)式代入(8)式, 在[t1, tc]积分得

    DcD1dD(1D)2=εpcεp11ω2[ln(1+εp)+ω21+εp]dεp (22)

    式中: 为壳体发生贯穿性破坏的损伤度。经过整理, 最终得

    εpcεp1(1+εpc)(1+εp1)ω22+[(epc)2(ep1)2DcD1(1Dc)(1D1)]ω2+(1+εpc)[(epc)22epc+2](1+εp1)[(ep1)22ep1+2]=0 (23)

    式中: 为第一阶段壳体材料径向塑形应变。由(6)式知满足以下关系式

    ep1=ln(1+εp1) (24)

    (23) 式为PELE弹丸壳体的损伤断裂方程。第1阶段的自然应变率和第2阶段的初始半径R2、径向速度u0可以通过实验得到, 通过R2u0可以求出而壳体材料参数σ0KηDc均可以通过查找手册和文献获得, 最终可以得到一个确定的临界应变作为损伤断裂判据。

    将(17)式、(19)式代入(21)式, 在[t1, tc]内积分可得壳体发生贯穿断裂的时刻tc

    tct1=R2(εpcεp1)/u0 (25)

    对于35CrMnSi, 参照文献[10], 壳体特征参数取η=-12 kPa·s, K=1.66 GPa, ωλ=3.5 kJ/m2, D0=10-5, Dc=0.2, σ0=0.72 GPa。当初始半径R0=15 mm的PELE弹丸以1 km/s速度着靶时, 结合实验数据, 得到第1阶段的自然应变率和第2阶段的初始半径R2及径向膨胀速度u0, 见表 4

    表  4  3个时刻膨胀断裂过程的参数
    Table  4.  Parameters of expansion and fracture at three time points
    t/(μs) ˙ep/(104 s1) R2/(mm) u0/(m/s)
    8.4 2.62 17.7 127
    7.2 1.91 18.6 142
    6.5 1.71 19.9 153
    下载: 导出CSV 
    | 显示表格

    将第2阶段的初始半径R2及径向膨胀速度u0代入(19)式, 得到第2阶段的初始应变率将第一阶段代入(13)式可以求出代入(14)式求出D1; 同时将代入(24)式求出代入(23)式, 得到3个典型时刻临界应变图 6所示。从图 6中可以看出, 在一定速度范围内, 随着弹丸着靶速度的增加, 壳体断裂的临界应变明显降低。

    图  6  PELE壳体断裂的临界应变
    Figure  6.  Critical strain of the PELE

    (1) PELE弹丸在侵彻过程中壳体膨胀断裂可以分为加速膨胀和匀速膨胀两个阶段, 在封加波提出的损伤度模型的基础上, 得到了PELE弹丸壳体的损伤断裂方程和断裂判据。

    (2) 根据PELE弹丸壳体损伤断裂方程, 计算出PELE弹丸壳体膨胀断裂过程的参数, 分析了弹丸着靶速度与壳体断裂之间的关系, 可为PELE弹丸设计提供理论指导。

    致谢: 感谢胡玉涛博士和曹雷博士在公式推导和实验方面给予的帮助。
  • 图  SHPB离散元模型结构

    Figure  1.  Structure of discrete element model of SHPB

    图  子弹的形状和尺寸

    Figure  2.  Shape and dimensions of the projectile

    图  锥形子弹撞击速度为18 m/s时入射和透射波形

    Figure  3.  Incident and transmitted waves created by the impact of the conical projectile at 18 m/s

    图  Hopkinson杆上记录的应力波波形

    Figure  4.  Stress waves recorded on the Hopkinson bars

    图  玻璃试样中测量圆的分布

    Figure  5.  Distribution of probing circles in the glass specimen

    图  应力和应变率随应变的变化

    Figure  6.  Dependence of stress and strain rate on strain

    图  试件两侧的应力时程曲线以及内部应力不均匀系数

    Figure  7.  Temporal profile of stress on the surfaces facing the incident bar and the transmitted bar and the stress inhomogeneity coefficient

    图  试件内部应力时程曲线(μ=0.1˙ε=700s1

    Figure  8.  Stress history curve of specimen (μ=0.1˙ε=700s1

    图  不同时刻试件内部裂纹分布(a)和破碎形态(b)(μ=0.1˙ε=700s1

    Figure  9.  Crack distributions (a) and fragmentation morphologies (b) of the specimen at different time (μ=0.1˙ε=700s1

    图  10  不同应变率下无机玻璃的应力-应变曲线

    Figure  10.  Stress-strain curves of inorganic glass at different strain rates

    图  11  无摩擦时试件内部的应力时程曲线

    Figure  11.  Temporal profile of stress inside a specimen without the boundary friction

    图  12  无摩擦时试件内部的破坏模式

    Figure  12.  Evolution of internal damage process inside a specimen without the boundary friction

    图  13  ν0 时材料的失效破坏演化

    Figure  13.  Evolution of material failure when ν0

    图  14  不同摩擦力和泊松比条件下试件的宏观应力时程曲线(˙ε≈700 s–1

    Figure  14.  Macroscopic stress histories of the specimen with various boundary friction and Poisson’s ratio (˙ε≈700 s–1)

    图  15  SHPB离散元模型结构

    Figure  15.  Structure of discrete element model of SHPB

    图  16  v0=9 m/s时无机玻璃动态巴西圆盘劈裂破坏过程:(a)压力时程曲线,(b)试样在各个时刻的破坏形貌

    Figure  16.  Typical Brazilian disk splitting process under dynamic loading (v0=9 m/s): (a) the pressure history; (b) failure patterns of the disk at different time

    图  17  不同冲击速度下巴西圆盘试件的加载压力时程曲线

    Figure  17.  Pressure histories of the specimens under different impact velocities

    图  18  压缩与拉伸强度的动态增强因子与应变率的关系

    Figure  18.  DIF of compressive and tensile strengths as a function of strain rate

    表  1  SHPB数值实验的离散元模型的主要微观参数

    Table  1.   Main microscopic parameters of discrete element model in numerical experiments of SHPB

    MaterialEffective modulus
    of linear contact/
    GPa
    Normal-to-shear
    stiffness ratio of
    linear contact
    Minimum radius of
    particles/mm
    Size ratio of maximum
    and minimum
    particles
    Porosity
    Steel bar1904.00.1001.50.15
    Inorganic glass 632.10.0261.50.10
    MaterialEffective modulus of
    flat-joint contact/
    GPa
    Normal-to-shear
    stiffness ratio of
    flat-joint contact
    Tensile strength of
    flat-joint contact/
    GPa
    Shear strength of
    flat-joint
    contact/GPa
    Density of
    particles/
    (kg·m–3)
    Steel bar1904.0100010008800
    Inorganic glass 632.10.0730.352444
    下载: 导出CSV

    表  2  石英玻璃宏观参数的数值模拟结果与文献数据的对比

    Table  2.   Macroscopic parameters of quartz glass: comparison of simulation results with ones published in the literatures

    MethodE/GPaρe/(kg·m–3)σc/MPaσt/MPaσb/MPaKIc/(MPa·m1/2)ν
    DEM simulation72.522006104767.50.960.17
    Refs.[1617]72.52200500–11004960–700.810.17
    下载: 导出CSV
  • [1] ZHANG X, ZOU Y, HAO H, et al. Laboratory test on dynamic material properties of annealed float glass [J]. International Journal of Protective Structures, 2012, 3(4): 407–430. doi: 10.1260/2041-4196.3.4.407
    [2] PERONI M, SOLOMOS G, PIZZINATO V, et al. Experimental investigation of high strain-rate behaviour of glass [C]// Applied Mechanics and Materials. Trans Tech Publications, 2011, 82: 63–68.
    [3] 王振, 张超, 王银茂, 等. 飞机风挡无机玻璃在不同应变率下的力学行为 [J]. 爆炸与冲击, 2018, 38(2): 295–301.

    WANG Z, ZHANG C, WANG Y M, et al. Mechanical behaviours of aeronautical inorganic glass at different strain rates [J]. Explosion and Shock Waves, 2018, 38(2): 295–301.
    [4] NIE X, CHEN W W, WERESZCZAK A A, et al. Effect of loading rate and surface conditions on the flexural strength of borosilicate glass [J]. Journal of the American Ceramic Society, 2009, 92(6): 1287–1295. doi: 10.1111/jace.2009.92.issue-6
    [5] NIE X, CHEN W W. Rate and surface treatment effect on the strength of boro-glass [C]//11th International Congress and Exhibition on Experimental and Applied Mechanics. Orlando, Florida: Society for Experimental Mechanics, 2008: 122–123.
    [6] NIE X, CHEN W W, SUN X, et al. Dynamic failure of borosilicate glass under compression/shear loading experiments [J]. Journal of the American Ceramic Society, 2007, 90(8): 2556–2562. doi: 10.1111/jace.2007.90.issue-8
    [7] ZHANG X, HAO H, MA G. Dynamic material model of annealed soda-lime glass [J]. International Journal of Impact Engineering, 2015, 77: 108–119. doi: 10.1016/j.ijimpeng.2014.11.016
    [8] SUN X, LIU W, CHEN W, et al. Modeling and characterization of dynamic failure of borosilicate glass under compression/shear loading [J]. International Journal of Impact Engineering, 2009, 36(2): 226–234. doi: 10.1016/j.ijimpeng.2008.01.014
    [9] 臧孟炎, 李军, 雷周. 基于DEM的两层结构夹层玻璃冲击破坏特性研究 [J]. 科学技术与工程, 2009, 9(3): 549–553. doi: 10.3969/j.issn.1671-1815.2009.03.005

    ZANG M Y, LI J, LEI Z. Study on impact fracture behavior of Bi-layer laminated glass based on DEM [J]. Science Technology and Engineering, 2009, 9(3): 549–553. doi: 10.3969/j.issn.1671-1815.2009.03.005
    [10] 臧孟炎, 雷周, 尾田十八. 汽车玻璃的静力学特性和冲击破坏现象 [J]. 机械工程学报, 2009, 45(2): 268–272.

    ZANG M Y, LEI Z, ODA J. Static characteristic and impact fracture behavior of automobile glass [J]. Journal of Mechanical Engineering, 2009, 45(2): 268–272.
    [11] CUNDALL P A. A computer model for simulating progressive, large-scale movement in blocky rock system [C]//Proceedings of the International Symposium on Rock Mechanics, 1971: 129–136.
    [12] POTYONDY D O, CUNDALL P A. A bonded-particle model for rock [J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(8): 1329–1364. doi: 10.1016/j.ijrmms.2004.09.011
    [13] CHEN W W, SONG B. Split Hopkinson (Kolsky) bar: design, testing and applications [M]. New York: Springer Science & Business Media, 2011: 37–49.
    [14] LI X, ZOU Y, ZHOU Z. Numerical simulation of the rock SHPB test with a special shape striker based on the discrete element method [J]. Rock Mechanics and Rock Engineering, 2014, 47(5): 1693–1709. doi: 10.1007/s00603-013-0484-6
    [15] 熊迅, 李天密, 马棋棋, 等. 石英玻璃圆环高速膨胀碎裂过程的离散元模拟 [J]. 力学学报, 2018, 50(3): 622–632.

    XIONG X, LI T M, MA Q Q, et al. Discrete element simulations of the high velocity expansion and fragmentation of quartz glass rings [J]. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(3): 622–632.
    [16] 王玉芬, 刘连城. 石英玻璃 [M]. 北京: 化学工业出版社, 2007.

    WANG Y F, LIU L C. Quartz glass [M]. Beijing: Chemical Industry Press, 2007.
    [17] 王承遇, 卢琪, 陶瑛. 玻璃的脆性(一) [J]. 玻璃与搪瓷, 2011, 39(6): 37–43. doi: 10.3969/j.issn.1000-2871.2011.06.009

    WANG C Y, LU Q, TAO Y. Brittleness of glass [J]. Glass & Enamel, 2011, 39(6): 37–43. doi: 10.3969/j.issn.1000-2871.2011.06.009
    [18] 宋力, 胡时胜. SHPB测试中的均匀性问题及恒应变率 [J]. 爆炸与冲击, 2005, 5(3): 207–216. doi: 10.3321/j.issn:1001-1455.2005.03.003

    SONG L, HU S S. Stress uniformity and constant strain rate in SHPB test [J]. Explosion and Shock Waves, 2005, 5(3): 207–216. doi: 10.3321/j.issn:1001-1455.2005.03.003
    [19] 邓志方, 黄西成, 谢若泽. SHPB实验入射波形分析[C]//中国计算力学大会, 2010.

    DENG Z F, HUANG X C, XIE R Z. Analysis of incident wave in SHPB experiments [C]//Chinese Conference on Computational Mechanics, 2010.
    [20] 陶俊林, 田常津, 陈裕泽, 等. SHPB系统试件恒应变率加载实验方法研究 [J]. 爆炸与冲击, 2004, 24(5): 413–418. doi: 10.3321/j.issn:1001-1455.2004.05.006

    TAO J L, TIAN C J, CHEN Y Z, et al. Investigation of experimental method to obtain constant strain rate of specimen in SHPB [J]. Explosion and Shock Waves, 2004, 24(5): 413–418. doi: 10.3321/j.issn:1001-1455.2004.05.006
    [21] SHEIKH M Z, WANG Z, DU B, et al. Static and dynamic Brazilian disk tests for mechanical characterization of annealed and chemically strengthened glass [J]. Ceramics International, 2019, 45(6): 7931–7944. doi: 10.1016/j.ceramint.2019.01.106
  • 加载中
图(18) / 表(2)
计量
  • 文章访问数:  6003
  • HTML全文浏览量:  3610
  • PDF下载量:  75
出版历程
  • 收稿日期:  2019-01-22
  • 修回日期:  2019-03-07
  • 刊出日期:  2019-04-25

目录

/

返回文章
返回