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

一种散体材料SHPB被动围压试验体应力修正方法

魏久淇 张春晓 曹少华 王世合 李磊

魏久淇, 张春晓, 曹少华, 王世合, 李磊. 一种散体材料SHPB被动围压试验体应力修正方法[J]. 爆炸与冲击, 2020, 40(12): 124201. doi: 10.11883/bzycj-2019-0411
引用本文: 魏久淇, 张春晓, 曹少华, 王世合, 李磊. 一种散体材料SHPB被动围压试验体应力修正方法[J]. 爆炸与冲击, 2020, 40(12): 124201. doi: 10.11883/bzycj-2019-0411
WEI Jiuqi, ZHANG Chunxiao, CAO Shaohua, WANG Shihe, LI Lei. A volume stress correction method for SHPB passiveconfined pressure of granular materials[J]. Explosion And Shock Waves, 2020, 40(12): 124201. doi: 10.11883/bzycj-2019-0411
Citation: WEI Jiuqi, ZHANG Chunxiao, CAO Shaohua, WANG Shihe, LI Lei. A volume stress correction method for SHPB passiveconfined pressure of granular materials[J]. Explosion And Shock Waves, 2020, 40(12): 124201. doi: 10.11883/bzycj-2019-0411

一种散体材料SHPB被动围压试验体应力修正方法

doi: 10.11883/bzycj-2019-0411
详细信息
    作者简介:

    魏久淇(1990- ),男,硕士,工程师,weijiuqi61489@163.com

    通讯作者:

    张春晓(1980- ),男,硕士,副研究员,cxz_007@163.com

  • 中图分类号: O341; O344

A volume stress correction method for SHPB passiveconfined pressure of granular materials

  • 摘要: 本文利用有限元仿真给出了一种修正方法,并用数值仿真和试验验证了该方法的可靠性。研究表明:散体材料SHPB被动围压试验中,试样厚度远小于厚壁圆筒长度时,端部效应会导致厚壁圆筒不均匀凸出变形,计算材料的体应力-应变关系不能将厚壁圆筒应力状态简化为平面应力问题;厚壁圆筒处于弹性状态下,通过厚壁圆筒理论计算出的径向力与真实径向力存在一定比例关系,在一定范围内,折算系数与试样实时厚度呈二次函数关系。
  • 散体介质是一种非平衡态复杂的能量耗散体系,具有优良的削波耗能特性[1],在地震工程、采矿工程和防护工程领域受到广泛应用,如砂土常作为防护结构分配层的填充材料。该类工程在地震、爆炸冲击等强荷载作用下的响应与材料在高应变率下的动态力学特性密切相关。分离式霍普金森压杆(SHPB)是研究材料中高应变率力学特性的常用设备,散体材料往往需要厚壁圆筒和垫块控制材料的初始条件,从而实现该材料的SHPB被动围压试验,国内外诸多学者[2-16]利用该方法对砂土等散体材料的动态力学特性开展了大量研究。SHPB被动围压试验中,试样及厚壁圆筒内部的应力场和应变场很难通过现有的技术直接测量,目前得到材料体应力-应变关系的主要手段是将厚壁圆筒的受力状态简化为平面应力问题[17],测量厚壁圆筒外表面(如图1所示)的环向应变εc,由下式计算得到,

    图  1  散体材料受约束的几何结构
    Figure  1.  Geometric structure constrained by SHPB test for bulk material
    σrr=0.5(a21)Ecεc
    (1)
    εrr=0.5εc[(1νc)+(1+νc)a2]
    (2)
    p=[σzz(t)+2σrr]/3
    (3)
    εV=εzz+2εrr
    (4)

    式中:σrr为极坐标下径向应力;εrr为极坐标下的径向应变;εzz为轴向应变;Ecνc分别为厚壁圆筒的弹性模量和泊松比,εc为应变片测到的厚壁圆筒环向应变,a为厚壁圆筒的外径与内径的比;p为体应力,εV为体应变。

    对于固体材料而言,因不存在装样困难,在满足SHPB均匀性假定的条件下,试样的初始厚度可以做到与厚壁圆筒长度等长,在小应变范围内,材料的体应力应变利用式(1)~(4)计算并无不妥;当试样厚度与厚壁圆筒不等长时,根据Forquin等[18]分析厚壁圆筒会发生不均匀凸出变形,图2样厚度较长时,试验时厚壁圆筒的几何变形图。由圣维南原理易知,厚壁圆筒L段受力状态仍可以看成平面应力问题,式(1)~(4)仍有效。然而由于散体材料波阻抗远小于压杆的波阻抗,试样厚度很小才能满足SHPB均匀性假定条件,加之受装样限制,试样厚度远小于厚壁圆筒的长度,此时端部效应可能会覆盖整个不均匀凸出区,公式(1)~(4)的有效性需要予以探讨。以往研究石英砂的动态力学特性时,大多只分析砂子的轴向应力应变,Luo等[11]没有考虑端部效应的影响直接引用式(1)~(4)计算砂子的体应力-应变关系,这有失严谨。魏久淇等[5]和文祝等[6]认为由式(1)计算的材料径向应力与真实应力存在折算系数关系,但并未给出合理的验证。Bragov等[15]注意到了该问题,并设计了一个阶梯型的圆筒(如图3所示),装样时将试样置于厚壁对应位置,该方法可以极大限度弱化端部效应,计算材料体应力应变时可直接引用式(1)~(4)。然而该方法存在两个弊端:一是圆筒薄壁越薄,端部效应对试验结果的影响越小,但圆筒加工难度越大;二是试验时无法改变试样的初始厚度。

    图  2  试验时厚壁圆筒的几何变形图
    Figure  2.  Geometric deformation diagram of thick-walled cylinder in test
    图  3  阶梯型套筒[15]
    Figure  3.  Ladder sleeve[15]

    本文中,不改变厚壁圆筒结构,使用传统厚壁圆筒固定装样的方法,利用Autodyn软件分析了试样实时厚度对材料真实体应力应变的影响,基于厚壁圆筒理论提出一种修正方法,并对非饱和钙质砂进行准一维应变冲击试验和数值仿真,验证修正方法的可靠性,以期对开展散体材料SHPB被动围压进行指导。

    Autodyn可用于分析气体、液体和固体及它们之间耦合作用的高速、瞬态动力学的问题,本文用该软件辅助分析特殊工况下(试样厚度远小于厚壁圆筒长度),式(1)~(4)的有效性,并给出科学计算材料体应力应变的思路。

    假定试样为准一维应变状态,图4为试样应力均匀时厚壁圆筒的受力示意图,此模型忽略了试样与厚壁圆筒内壁摩擦力。由图4可知,当试样应力均匀时,厚壁圆筒内壁承受沿材样实时厚度Ls均布径向内推力σss,根据连续性条件可知σss=σrr,由式(1)可知εc=2σrr/[Ec(α21)]。利用Autodyn建立厚壁圆筒模型,对厚壁圆筒内壁试样实时厚度Ls区施加给定σss(t),输出厚壁圆筒外中点环向应变εs,对比εsεc,若εs=εc,则说明式(1)~(4)有效,若εsεc,则反之。图5为建好的厚壁圆筒模型,选用朗格朗日算法,厚壁圆筒为自由边界条件。

    图  4  厚壁圆筒受力分析
    Figure  4.  Force analysis of thick-walled cylinder
    图  5  厚壁圆筒模型
    Figure  5.  Thick-walled cylinder model

    图6为给定的σss时程曲线,由图6可知σss为动态力,可以真实模拟试验时厚壁圆筒受到的动态径向内推力。

    图  6  径向力σss时程曲线
    Figure  6.  Time-history curve of radial force σss

    对于不同的LcLs图7给出了公式值εc和数值计算值εs(Lc,Ls)对比图,其中LcLs的单位为mm。公式值εcσss带入式(1)计算出的厚壁圆筒外表面的环向应变。由图7可知,当试样实时厚度大于30 mm,式(1)有效,当试样厚度小于20 mm时,式(1)无效。研究发现:当砂样实时厚度小于20 mm,公式值εc与数值计算值εs(Lc,Ls)存在折算关系:k=εc/εs(Lc,Ls)k为折算系数,试样实时厚度仅影响折算系数的大小。

    图  7  应变公式值(εc)与数值计算值(εs(Lc, Ls))的对比
    Figure  7.  Comparison diagram between strain by fomula (εc) and that by simulation (εs(Lc, Ls))

    本文试验中,试样初始厚度为10 mm,且最大应变约为0.25,因此需着重分析试样实时厚度10 mm至7 mm的折算系数。图8为折算系数k与试样实时厚度Ls的关系图,根据图8中的数据拟合得

    图  8  折算系数与试样实时厚度的关系图
    Figure  8.  relation diagram between conversion coefficient and real-time thickness of samples
    k=0.023L2s0.593Ls+5.207(7mmLs10mm)
    (5)

    因此,当7mmLs10mm时,计算钙质砂的体应力应变时应采用以下公式

    σrr=0.5k(α21)Ecεc
    (6)
    εrr=0.5k[(1νc)+(1+νc)α2]εc
    (7)
    p=1/3(σzz+2kσrr)
    (8)
    εV=εzz+2kεrr
    (9)

    为了验证上述仿真方法的可靠性,本文开展了数值模拟,图9为试样实时厚度与厚壁圆筒等长(Ls=Lc)时,本文模拟计算值εs(Lc, Ls)与εc的(式(1)理论值)对比图。由图9可知,当试样实时厚度与厚壁圆筒相等时,仿真计算值与理论值吻合较好,因此本文Autodyn有限元分析方法具有一定的可信度。

    图  9  验证性模拟
    Figure  9.  Verification simulation

    为进一步验证本文方法的可靠性,本文开展了含水率为25%和30%非饱和钙质砂SHPB被动围压高速冲击试验,每种含水率重复3次试验。含水率为25%和30%非饱和钙质砂受高速冲击压缩可转变为饱和钙质砂。水在饱和钙质砂中对应力应变关系起主导作用,由于水的等向性,理论上饱和钙质砂的轴向应力等于径向应力。非饱和钙质砂受高速冲击压缩转变为饱和钙质砂时,若式(6)计算的径向应力和轴向应力相等,则能进一步验证本文方法的可靠性。

    本文所用材料为钙质砂,图10为钙质砂试样,呈乳白色,主要成分为CaCO3,掺杂有少量红色颗粒(碳酸镁)。原始砂样经过24 h充分烘干后,经标准筛筛分称重后的级配曲线见图11,考虑到试样直径及试样厚度,大颗粒的存在无法满足试样均匀性,极小颗粒量少不易称量,故将砂样中粒径小于0.15 mm和大于1.18 mm的粒径剔除,剔除后剩余粒径的砂样组分不变,钙质砂的比重为2.82,最大最小干密度分别为1.378和1.183 g/cm3,实际工况如表1所示。

    表  1  试验工况表
    Table  1.  Summary of SHPB tests
    编号试验材料含水率/%砂质量/g装样厚度/mm干密度/(g∙cm−3)相对密实度/%气压/MPa
    G60-0.2-01钙质砂2513.9010.021.2958.730.2
    G60-0.2-022510.001.2960
    G60-0.2-0325 9.981.3061.31
    G60-0.2-043010.021.2958.73
    G60-0.2-053010.001.2960
    G60-0.2-063010.021.2958.73
    下载: 导出CSV 
    | 显示表格
    图  10  试验砂样
    Figure  10.  Sand specimens tested
    图  11  砂样颗粒级配曲线
    Figure  11.  Grain size distribution

    本文所用SHPB试验装置为直径37 mm的铝杆,如图12(a)所示。子弹长度、入射杆以及透射杆分别为400 mm、2 000 mm和2 000 mm。由于钙质砂为松散介质,为精确控制试样的初始条件,试验时用厚壁圆筒和垫块束缚成试件。垫块和厚壁圆筒如图12(b)图12(c)所示。厚壁圆筒材质为4340钢,弹性模量为210 GPa, 泊松比为0.3,内径37.1 mm,外径为47 mm,长为94 mm;垫块材质跟压杆材质完全相同,直径为37 mm,厚度为30 mm。

    图  12  试验设备
    Figure  12.  Test equipment

    试样装样步骤参考Song等[2]研究成果,整形器为厚1 mm直径10 mm橡胶片。图13是试验时的原始典型波形。砂样动态应力平衡是评判试验数据有效的重要方法,图14为钙质砂试样的动态应力平衡图,由图可知试样基本满足应力平衡条件。

    图  13  试验原始波形
    Figure  13.  Test the original waveform
    图  14  砂样动态应力平衡
    Figure  14.  Dynamic stress equilibrium in the sand specimens

    图15为含水率25%、30%钙质砂的应力-应变曲线。厚壁圆筒的屈服强度为835 MPa,由拉梅公式[19]易知厚壁圆筒内应力(等于砂子的径向应力)大于325 MPa时,厚壁圆筒才会屈服。由图14可知钙质砂的最大轴向应力为102 MPa,易知厚壁处于弹性状态。本文砂样的轴向应变为10−1量级,由式(5)~(7)易知径向应变为10−4量级,故钙质砂可以看成准一维应变状态。因此可得,

    图  15  含水率25%、30%钙质砂的应力应变
    Figure  15.  Stress and strain of calcareous sand with 25%, 30% water content
    ρ=ρ0/(1εz)
    (10)

    式中:ρ0为初始干密度,εz为轴向应变,ρ为砂样应变εz的实时密度。由土力学[20]易知,

    e=wsatGs
    (11)
    e=Gsρ1
    (12)

    式中:e为孔隙比,wsat为饱和含水率,Gs为骨架密度。由式(10)、(11)和(12)易得含水率25%和30%非饱和钙质砂受高速冲击压缩转变为饱和时的应变拐点,分别记为Sr1Sr2(见图15)。由图14可知:钙质砂含水率越高,受冲击压缩更易饱和,达到饱和点后,由于水不可压缩,其应力-应变曲线迅速趋陡峭。图14中钙质砂的径向力σrr是由式(5)~(6)计算得到的,由图可知,饱和状态时,钙质砂的轴向应力和修正过的径向应力基本相等,这主要是因为水在饱和钙质砂中起主导作用,水对四周力相等。综上可知,本文方法具有可靠性。

    本文利用有限元分析给出了一种解决办法,并用数值仿真和试验两种方法验证本文方法的可靠性,主要得出以下结论。

    (1)散体材料SHPB被动围压试验时,试样厚度远小于厚壁圆筒长度时,由于端部效应导致厚壁圆筒不均匀凸出变形,求解材料的体应力-应变关系时不能将厚壁圆筒应力状态简化为平面应力问题。

    (2)厚壁圆筒处于弹性状态下,通过厚壁圆筒理论计算出的材料径向力与材料真实径向力存在一定比例关系,在一定范围内,折算系数与试样实时厚度存在二次函数关系。

    (3)试样厚度、厚壁圆筒长度、厚壁圆筒内外径之比以及厚壁圆筒材料特性都会影响折算系数,各因素之间还存在多种耦合关系,解析解也尚未给出,仍需进一步探讨。开展类似试验时,针对具体情况,可利用本文提供的方法分析折算系数。

  • 图  1  散体材料受约束的几何结构

    Figure  1.  Geometric structure constrained by SHPB test for bulk material

    图  2  试验时厚壁圆筒的几何变形图

    Figure  2.  Geometric deformation diagram of thick-walled cylinder in test

    图  3  阶梯型套筒[15]

    Figure  3.  Ladder sleeve[15]

    图  4  厚壁圆筒受力分析

    Figure  4.  Force analysis of thick-walled cylinder

    图  5  厚壁圆筒模型

    Figure  5.  Thick-walled cylinder model

    图  6  径向力σss时程曲线

    Figure  6.  Time-history curve of radial force σss

    图  7  应变公式值(εc)与数值计算值(εs(Lc, Ls))的对比

    Figure  7.  Comparison diagram between strain by fomula (εc) and that by simulation (εs(Lc, Ls))

    图  8  折算系数与试样实时厚度的关系图

    Figure  8.  relation diagram between conversion coefficient and real-time thickness of samples

    图  9  验证性模拟

    Figure  9.  Verification simulation

    图  10  试验砂样

    Figure  10.  Sand specimens tested

    图  11  砂样颗粒级配曲线

    Figure  11.  Grain size distribution

    图  12  试验设备

    Figure  12.  Test equipment

    图  13  试验原始波形

    Figure  13.  Test the original waveform

    图  14  砂样动态应力平衡

    Figure  14.  Dynamic stress equilibrium in the sand specimens

    图  15  含水率25%、30%钙质砂的应力应变

    Figure  15.  Stress and strain of calcareous sand with 25%, 30% water content

    表  1  试验工况表

    Table  1.   Summary of SHPB tests

    编号试验材料含水率/%砂质量/g装样厚度/mm干密度/(g∙cm−3)相对密实度/%气压/MPa
    G60-0.2-01钙质砂2513.9010.021.2958.730.2
    G60-0.2-022510.001.2960
    G60-0.2-0325 9.981.3061.31
    G60-0.2-043010.021.2958.73
    G60-0.2-053010.001.2960
    G60-0.2-063010.021.2958.73
    下载: 导出CSV
  • [1] 钱七虎, 王明洋. 三相介质饱和土自由场中爆炸波的传播规律 [J]. 爆炸与冲击, 1994, 14(2): 97–104.

    QIAN Q H, WANG M Y. Propagation of explosive wave in the free-field of three-phase saturated soil [J]. Explosion and Shock Waves, 1994, 14(2): 97–104.
    [2] SONG B, CHEN W, 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.
    [3] 赵章泳, 邱艳宇, 紫民, 等. 含水率对非饱和钙质砂动力特性影响的试验研究 [J]. 爆炸与冲击, 2019, 40(2): 023102. DOI: 10.11883/bzycj-2019-0066.

    ZHAO Z Y, QIU Y Y, ZI Min, et al. Experimental study on dynamic compression of unsaturated calcareous sand [J]. Explosion and Shock Waves, 2019, 40(2): 023102. DOI: 10.11883/bzycj-2019-0066.
    [4] 魏久淇, 王明洋, 邱艳宇, 等. 钙质砂动态力学特性试验研究 [J]. 振动与冲击, 2018, 37(24): 7–12. DOI: 10.13465/j.cnki.jvs.2018.24.002.

    WEI J Q, WANG M Y, QIU Y Y, et al. Impact compressive response of calcareous sand [J]. Journal of Vibration and Shock, 2018, 37(24): 7–12. DOI: 10.13465/j.cnki.jvs.2018.24.002.
    [5] 魏久淇, 吕亚茹, 刘国权, 等. 钙质砂一维冲击响应及吸能特性试验 [J]. 岩土力学, 2019, 40(1): 191–198, 206. DOI: 10.16285/j.rsm.2017.1235.

    WEI J Q, LÜ Y R, LIU G Q. et al One-dimensional impact responses and energy absorption of calcareous sand [J]. Rock and Soil Mechanics, 2019, 40(1): 191–198, 206. DOI: 10.16285/j.rsm.2017.1235.
    [6] 文祝, 邱艳宇, 紫民, 等. 钙质砂的准一维应变压缩试验研究 [J]. 爆炸与冲击, 2019, 39(3): 033103. 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): 033103. DOI: 10.11883/bzycj-2018-0015.
    [7] 于潇, 陈力, 方秦. 珊瑚砂中应力波衰减规律的实验研究 [J]. 岩石力学与工程学报, 2018, 37(6): 1520–1529. DOI: 10.13722/j.cnki.jrme.2018.0147.

    YU X, CHEN L, FANG Q. Experimental study on the attenuation of stress wave in coral sand [J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(6): 1520–1529. DOI: 10.13722/j.cnki.jrme.2018.0147.
    [8] BARR A D, CLARKE S D, TYAS A, et al. Effect of moisture content on high strain rate compressibility and particle breakage in loose sand [J]. Experimental Mechanics, 2018, 58(8): 1331–1334. DOI: 10.1007/s11340-018-0405-4.
    [9] ROSS C A, THOMPSON P Y, CHARLIE W A, et al. Transmission of pressure waves in partially saturated soils [J]. Experimental Mechanics, 1989, 29(1): 80–83. DOI: 10.1007/BF02327786.
    [10] MARTIN B E, KABIR Md E, CHEN W. Undrained high-pressure and high strain-rate response of dry sand under triaxial loading [J]. International Journal of Impact Engineering, 2013, 54: 51–63. DOI: 10.1016/j.ijimpeng.2012.10.008.
    [11] 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.
    [12] OMIDVAR M, ISKANDER M, BLESS S. Stress-strain behavior of sand at high strain rates [J]. International Journal of Impact Engineering, 2012, 49: 192–213. DOI: 10.1016/j.ijimpeng.2012.03.004.
    [13] 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.
    [14] MARTIN B E, CHEN W N, SONG B, et al. Moisture effects on the high strain-rate behavior of sand [J]. Mechanics of Materials, 2009, 41(6): 786–798. DOI: 10.1016/j.mechmat.2009.01.014.
    [15] 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.
    [16] BRAGOV A M, GRUSHEVSKY G M, LOMUNOV A K. Use of the Kolsky method for confined tests of soft soils [J]. Experimental Mechanics, 1996, 36(3): 237–242. DOI: 10.1007/BF02318013.
    [17] RAVI-CHANDAR K, Ma Z. Inelastic deformation in polymers under multiaxial compression [J]. Mechanics of Time-Dependent Materials, 2000, 4(4): 333–357. DOI: 10.1023/a:1026570826226.
    [18] FORQUIN P, GARY G, GATUINGT F. A testing technique for concrete under confinement at high rates of strain [J]. International Journal of Impact Engineering, 2008, 35(6): 425–446. DOI: 10.1016/j.ijimpeng.2007.04.007.
    [19] 徐秉业. 应用弹塑性力学[M]. 北京: 清华大学出版社, 1995.
    [20] 国家质量技术监督局, 中华人民共和国建设部. GBT/50123 土工试验方法标准[S]. 北京: 中国计划出版社, 1999.
  • 期刊类型引用(1)

    1. 陈昊祥,李杰,邓树新,王德荣,王明洋. 散体介质SHPB被动围压试验体应力计算的理论修正方法. 爆炸与冲击. 2022(06): 115-121 . 本站查看

    其他类型引用(1)

  • 加载中
图(15) / 表(1)
计量
  • 文章访问数:  1558
  • HTML全文浏览量:  798
  • PDF下载量:  78
  • 被引次数: 2
出版历程
  • 收稿日期:  2019-10-24
  • 修回日期:  2020-01-19
  • 刊出日期:  2020-12-05

目录

/

返回文章
返回