• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊
高级检索 E-mail Alert

多胞材料的动态应力应变状态及其一致近似关系

王鹏 朱长锋 郑志军 虞吉林

引用本文:
Citation:

多胞材料的动态应力应变状态及其一致近似关系

    作者简介: 王鹏(1988-), 男, 博士, 助理研究员;
    通讯作者: 郑志军, zjzheng@ustc.edu.cn
  • 基金项目: 国家自然科学基金 11772330
    国家自然科学基金 11372308

  • 中图分类号: O347.3

Dynamic stress-strain states of cellular materials and a uniformly approximated relation

    Corresponding author: ZHENG Zhijun, zjzheng@ustc.edu.cn ;
  • CLC number: O347.3

  • 摘要: 多胞材料在高速冲击下呈现出逐层压溃的变形模式,塑性冲击波模型可以用来表征这种集中变形带的传播行为。本文中采用截面应力计算方法得到了随机蜂窝在恒速冲击下的一维应力分布,进而对冲击波的传播规律进行了分析。比较了高速冲击下由不同方法得到的冲击波速度与冲击速度的关系,结果表明R-PP-L(率无关,刚性-理想塑性-锁定)模型高估了冲击波速度,但R-PH(率无关,刚性-塑性硬化)模型以及一维冲击波理论得到的冲击波速度与有限元结果比较接近。冲击波速度与冲击速度在高速情形下趋于线性关系,但随着冲击速度的减小,冲击波速度不断减少并趋于常数。根据这一特征和塑性冲击波模型,发展了可以表征冲击波速度与冲击速度的关系、动态应力应变关系的一致近似模型。
  • 图 1  Voronoi蜂窝试件细观有限元模型

    Figure 1.  A cell-based finite element model of a Voronoi honeycomb specimen

    图 2  两种冲击速度下的应力分布图

    Figure 2.  One-dimensional stress distributions at two impact velocities

    图 3  一维的应力分布及其应力梯度分布

    Figure 3.  One-dimensional stress distributions and the corresponding stress gradients

    图 4  不同冲击速度下冲击波位置随时间的关系

    Figure 4.  Variation of shock front position with impact time at different impact velocities

    图 5  一维应力和应变分布及其理想化

    Figure 5.  One-dimensional stress and strain distributions and the idealizations

    图 6  不同方法得到的冲击波速度的比较

    Figure 6.  Comparison of shock wave speeds obtained using different methods

    图 7  冲击波速度与冲击速度的关系

    Figure 7.  Variations of the shock wave speed with the impact velocity

    图 8  波后应变和波后应力随冲击速度的变化

    Figure 8.  Variations of strain and stress behind shock front with impact velocity

    图 9  不同冲击速度下的速度分布

    Figure 9.  Velocity distributions at different impact velocities

    图 10  准静态及动态应力应变关系

    Figure 10.  Quasi-static and dynamic stress-strain relations

  • [1] REID S R, BELL W W, BARR R A. Structural plastic shock model for one-dimensional ring systems[J]. International Journal of Impact Engineering, 1983, 1(2):175-191.DOI:10.1016/0734-743X(83)90005-2.
    [2] REID S R, PENG C. Dynamic uniaxial crushing of wood[J]. International Journal of Impact Engineering, 1997, 19(5/6):531-570. DOI:10.1016/S0734-743X(97)00016-X.
    [3] ZHENG Z J, YU J L, LI J R. Dynamic crushing of 2D cellular structures:A finite element study[J]. International Journal of Impact Engineering, 2005, 32(1):650-664. DOI:10.1016/j.ijimpeng.2005.05.007.
    [4] LIU Y D, YU J L, ZHENG Z J, et al. A numerical study on the rate sensitivity of cellular metals[J]. International Journal of Solids and Structures, 2009, 46(22):3988-3998. DOI:10.1016/j.ijsolstr.2009.07.024.
    [5] 张新春, 刘颖, 章梓茂.组合蜂窝材料面内冲击性能的研究[J].工程力学, 2009, 26(6):220-225.
    ZHANG Xinchun, LIU Ying, ZHANG Zimao. Research on dynamic properties of suppercell honeycomb structures under in-plane impact loading[J]. Engineering Mechanics, 2009, 26(6):220-225.
    [6] 胡玲玲, 蒋玲.胞孔构型对金属蜂窝动态力学性能的影响机理[J].爆炸与冲击, 2014, 34(1):41-46. DOI:10.3969/j.issn.1001-1455.2014.01.008.
    HU Lingling, JIANG Ling. Mechanism of cell configuration affecting dynamic mechanical properties of metal honeycombs[J]. Explosion and Shock Waves, 2014, 34(1):41-46. DOI:10.3969/j.issn.1001-1455.2014.01.008.
    [7] ZOU Z, REID S R, TAN P J, et al. Dynamic crushing of honeycombs and features of shock fronts[J]. International Journal of Impact Engineering, 2009, 36(1):165-176. DOI:10.1016/j.ijimpeng.2007.11.008.
    [8] LIAO S F, ZHENG Z J, YU J L. Dynamic crushing of 2D cellular structures:Local strain field and shock wave velocity[J]. International Journal of Impact Engineering, 2013, 57(1):7-16. DOI:10.1016/j.ijimpeng.2013.01.008.
    [9] TAN P J, REID S R, HARRIGAN J J, et al. Dynamic compressive strength properties of aluminium foams. Part Ⅱ:"Shock"theory and comparison with experimental data and numerical models[J]. Journal of the Mechanics and Physics of Solids, 2005, 53(10):2206-2230. DOI:10.1016/j.jmps.2005.05.003.
    [10] TAN P J, REID S R, HARRIGAN J J, et al. Dynamic compressive strength properties of aluminium foams. Part Ⅰ:Experimental data and observations[J]. Journal of the Mechanics and Physics of Solids, 2005, 53(10):2174-2205. DOI:10.1016/j.jmps.2005.05.007.
    [11] ZHENG Z J, YU J L, WANG C F, et al. Dynamic crushing of cellular materials:A unified framework of plastic shock wave models[J]. International Journal of Impact Engineering, 2013, 53(1):29-43. DOI:10.1016/j.ijimpeng.2012.06.012.
    [12] 王长峰, 郑志军, 虞吉林.泡沫杆撞击刚性壁的动态压溃模型[J].爆炸与冲击, 2013, 33(6):587-593. DOI:10.3969/j.issn.1001-1455.2013.06.005.
    WANG Changfeng, ZHENG Zhijun, YU Jilin. Dynamic crushing models for a foam rod striking a rigid wall[J]. Explosion and Shock Waves, 2013, 33(6):587-593. DOI:10.3969/j.issn.1001-1455.2013.06.005.
    [13] ZHENG Z J, WANG C F, YU J L, et al. Dynamic stress-strain states for metal foams using a 3D cellular model[J]. Journal of the Mechanics and Physics of Solids, 2014, 72:93-114. DOI:10.1016/j.jmps.2014.07.013.
    [14] BARNES A T, RAVI-CHANDA K, KYRIAKIDES S, et al. Dynamic crushing of aluminum foams:Part Ⅰ:Experiments[J]. International Journal of Solids and Structures, 2014, 51(9):1631-1645. DOI:10.1016/j.ijsolstr.2013.11.019.
    [15] SUN Y L, LI Q M, MCDONALD S A, et al. Determination of the constitutive relation and critical condition for the shock compression of cellular solids[J]. Mechanics of Materials, 2016, 99:26-36. DOI:10.1016/j.mechmat.2016.04.004.
    [16] WANG P, WANG X K, ZHENG Z J, et al. Stress distribution in graded cellular materials under dynamic compression[J]. Latin American Journal of Solids and Structures, 2017, 14(7):1251-1272. DOI:10.1016/j.ijimpeng.2017.10.006.
    [17] WANG L L. Foundations of stress waves[M]. Amsterdam:Elsevier Science Ltd., 2007. DOI:10.1002/prep.200790014.
  • [1] 廖华林李根生 . 孔隙流体耦合效应对射流冲击应力分布的影响分析. 爆炸与冲击, 2006, 26(1): 84-90. doi: 10.11883/1001-1455(2006)01-0084-07
    [2] 王荣波田建华何莉华李泽仁赵剑衡吴廷烈 . 石英光纤探针在非金属材料冲击实验中的应用. 爆炸与冲击, 2006, 26(3): 284-288. doi: 10.11883/1001-1455(2006)03-0284-04
    [3] 郭昭亮范诚刘明涛汤铁钢刘仓理 . 膨胀环受力历史对应力应变关系的影响. 爆炸与冲击, 2016, 36(6): 819-824. doi: 10.11883/1001-1455(2016)06-0819-06
    [4] 蔡正宇丁圆圆王士龙郑志军虞吉林 . 梯度多胞牺牲层的抗爆炸分析. 爆炸与冲击, 2017, 37(3): 396-404. doi: 10.11883/1001-1455(2017)03-0396-09
    [5] 韩会龙张新春王鹏 . 负泊松比蜂窝材料的动力学响应及能量吸收特性. 爆炸与冲击, 2019, 39(1): 013103-1-013103-11. doi: 10.11883/bzycj-2017-0281
    [6] 姜楠徐全军龙源廖瑜林卫 . 大孔径静态破碎膨胀压力特性及布孔参数分析. 爆炸与冲击, 2015, 35(4): 467-472. doi: 10.11883/1001-1455(2015)04-0467-06
    [7] 朱兆祥胡时胜陆菁骆利民王礼立 . 冲击载荷下钢-铅层合体的弹粘塑性应力应变关系研究. 爆炸与冲击, 1989, 9(1): 1-10.
    [8] 经福谦 . 冲击波速度—粒子速度关系式的一个简单推导及其直线表达式适用范围的讨论. 爆炸与冲击, 1982, 2(3): 17-24.
    [9] 江厚满张若棋张寿齐 . 用冲击波应力历史曲线确定材料物态方程参数. 爆炸与冲击, 1998, 18(2): 177-181.
    [10] 郭靖华王祥林孙晓明 . 动态光测法研究冲击作用下的应力分布问题. 爆炸与冲击, 1987, 7(2): 146-153.
    [11] 孙悦慈佳祥孙东 . 一种冲击波速度传感器性能的检测装置. 爆炸与冲击, 2007, 27(5): 473-476. doi: 10.11883/1001-1455(2007)05-0473-04
    [12] 余尚江李科杰 . 混凝土结构内冲击波应力传感器设计及其行为. 爆炸与冲击, 2005, 25(4): 350-354. doi: 10.11883/1001-1455(2005)04-0350-05
    [13] 寇绍全 . 波动问题的一个逆问题从测得的应变波形建立材料的应力应变应变率关系. 爆炸与冲击, 1983, 3(1): 44-51.
    [14] 黄霞汤文辉蒋邦海 . 平面应变各向异性本构关系及在应力波传播模拟中的应用. 爆炸与冲击, 2010, 30(4): 383-389. doi: 10.11883/1001-1455(2010)04-0383-07
    [15] 金解放李夕兵常军然陶伟邱灿 . 循环冲击作用下岩石应力应变曲线及应力波特性. 爆炸与冲击, 2013, 33(6): 613-619. doi: 10.11883/1001-1455(2013)06-0613-07
    [16] 范春雷胡金伟陈大年王焕然谢书港 . 无氧铜平面冲击波实验的横向应力测试及屈服强度确定. 爆炸与冲击, 2008, 28(2): 110-115. doi: 10.11883/1001-1455(2008)02-0110-06
    [17] 张宝超潘景龙 . FRP约束混凝土快速荷载下应力应变关系初探. 爆炸与冲击, 2003, 23(5): 466-471.
    [18] 乔登江 . 空中爆炸冲击波(Ⅰ)基本理论. 爆炸与冲击, 1985, 5(4): 78-85.
    [19] 周佩杰王坚陶钢周杰 . 泡沫材料对冲击波的衰减特性. 爆炸与冲击, 2015, 35(5): 675-681. doi: 10.11883/1001-1455(2015)05-0675-07
    [20] 王礼立胡时胜 . 钛合金TB-2在高应变率下的动态应力应变关系. 爆炸与冲击, 1985, 5(1): 9-15.
  • 加载中
图(10)
计量
  • 文章访问数:  94
  • HTML全文浏览量:  204
  • PDF下载量:  14
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-08-09
  • 录用日期:  2018-03-06
  • 刊出日期:  2019-01-05

多胞材料的动态应力应变状态及其一致近似关系

    作者简介:王鹏(1988-), 男, 博士, 助理研究员
    通讯作者: 郑志军, zjzheng@ustc.edu.cn
  • 1. 中国科学技术大学近代力学系中国科学院材料力学行为和设计重点实验室, 安徽 合肥 230026
  • 2. 中国工程物理研究院总体工程研究所, 四川 绵阳 621999
基金项目:  国家自然科学基金 11772330国家自然科学基金 11372308

摘要: 多胞材料在高速冲击下呈现出逐层压溃的变形模式,塑性冲击波模型可以用来表征这种集中变形带的传播行为。本文中采用截面应力计算方法得到了随机蜂窝在恒速冲击下的一维应力分布,进而对冲击波的传播规律进行了分析。比较了高速冲击下由不同方法得到的冲击波速度与冲击速度的关系,结果表明R-PP-L(率无关,刚性-理想塑性-锁定)模型高估了冲击波速度,但R-PH(率无关,刚性-塑性硬化)模型以及一维冲击波理论得到的冲击波速度与有限元结果比较接近。冲击波速度与冲击速度在高速情形下趋于线性关系,但随着冲击速度的减小,冲击波速度不断减少并趋于常数。根据这一特征和塑性冲击波模型,发展了可以表征冲击波速度与冲击速度的关系、动态应力应变关系的一致近似模型。

English Abstract

  • 多胞材料在高速冲击下具有变形局部化的典型特征。Reid等[1-2]研究了由圆环组成的一维系统的动态冲击行为,提出了“结构冲击波”的概念,并在对木材的动态压溃实验中证实了压溃带的冲击波传播特征。Zheng等[3]、Liu等[4]在对随机蜂窝的动态冲击中观察到3种变形模式,即均匀模式、过渡模式和冲击模式。张新春等[5]、胡玲玲等[6]探讨了胞元构型对蜂窝结构动态冲击性能的影响。Zou等[7]认为蜂窝中传播的冲击波波阵面在高速冲击下具有胞元尺寸的宽度。Liao等[8]采用了局部应变场计算方法表征了蜂窝铝中一维冲击波的传播并得到了冲击波速度的大小。一系列的一维冲击波模型被发展来表征多胞材料中塑性压溃行为,如R-PP-L模型(率无关、刚性-理想塑性-锁定假设)[2, 9-10]和R-PH模型(率无关,刚性-塑性硬化假设)[11-13]。然而,大部分的冲击波模型都是基于准静态名义应力-应变曲线得到的,直接得到动态的应力应变关系在实验上是困难的。

    多胞材料在动态冲击下存在变形局部化和应力增强等现象,名义应力应变曲线失去了物理意义[4]。Zheng等[13]通过对泡沫铝的三维细观有限元模型的分析得到了波后应力应变状态点,并提出了D-R-PH(动态,刚性-塑性硬化)模型来表征动态应力-应变行为。Barnes等[14]、Sun等[15]分别对开孔泡沫和闭孔泡沫进行实验和数值研究,也得到了多胞材料的动态应力应变状态。但这些研究都是针对冲击模式下的应力应变状态,对过渡模式下的应力应变状态的认识仍然不清楚。此外,文献[13-15]中波后应力都是由撞击端的名义应力得到的,缺少对材料内部局部应力信息的直接表征。

    认识多胞材料在中速过渡模式下的应力应变行为,对于完整地理解多胞材料的动态本构关系有重要意义。本文中将采用截面应力计算方法研究基于细观有限元模型的随机蜂窝结构在恒速冲击下的应力分布,分析冲击波在随机蜂窝中的传播特性,表征随机蜂窝的动态应力应变关系。

    • 采用二维Voronoi技术构造随机蜂窝结构,详细建模过程参见文献[3]。如图 1,蜂窝试件的长和宽均为50 mm,试件的不规则度为0.3,由500个胞元组成。试件相对密度ρ0/ρs= 0.1,其中ρ0为蜂窝试样变形前的初始密度,ρs为基体材料的密度。蜂窝试件的平均胞元尺寸约为2.5 mm,其值定义为与平均胞元面积相等的圆的直径。

      图  1  Voronoi蜂窝试件细观有限元模型

      Figure 1.  A cell-based finite element model of a Voronoi honeycomb specimen

      采用ABAQUS/Explicit有限元软件对蜂窝材料的恒速压溃过程进行数值模拟。胞壁材料由S4R壳单元(4节点,减缩积分、沙漏控制、有限膜应变)进行模拟。通过网格收敛性分析将模型划分了10 480个壳单元,平均单元长度约为0.2 mm。基体材料属性设置为弹性-理想塑性,材料参数设置为泊松比ν = 0.3,弹性模量E= 66 GPa,屈服应力σy = 175 MPa,基体密度ρs = 2 700 kg/m3。所有接触面均定义为通用接触,摩擦因数取为0.02。为了模拟面内变形,模型中所有节点的面外变形均被限制。

      蜂窝试件的左右两端设置两个刚性面,左侧刚性面以恒定不变的速度v沿x方向向右压缩试件。在压缩过程中右侧刚性板为固支端,如图 1所示。

    • 对于蜂窝试件在承受匀速压缩的情况下,一维近似是适用的,因此可以通过截面工程应力来表征局部应力信息,详见文献[16]。不同的拉格朗日截面位置的工程应力由截面内力除以截面初始面积得到。截面的内力由两部分组成,一部分是通过基体材料的单元节点所传递的力,称为节点力,记为Fx(n);另一部分是由变形中胞壁接触所引起的接触力,记为Fx(c)。这两部分力由下式可以计算得到:

      $ \left\{ \begin{array}{l} F_x^{\left( {\rm{n}} \right)} = \sum\limits_{i = 1}^M {F_{xi}^{\left( {\rm{n}} \right)}} \\ F_x^{\left( {\rm{c}} \right)} = \sum\limits_{i = 1}^K {\left( {F_{xi}^{\left( {{\rm{cn}}} \right)} + F_{xi}^{\left( {{\rm{ct}}} \right)}} \right)} \end{array} \right. $

      式中:M为所考察的截面所贯穿的单元在参考构型中左侧的节点数,K是从撞击端到所考察截面的总的节点数目。对于恒速压缩情形,K包含了刚性端的节点。接触力Fx(c)的简化计算公式的算法得益于牛顿第三定律,完整的接触对的接触力大小相等方向相反,互相抵消之后剩余的部分为被截面所隔开的接触对的接触力。由于刚性面与试件撞击端位置重合,因此撞击端面上的接触力由刚性面输出接触力得到。Fxi(n)Fxi(cn)Fxi(ct)分别为任意节点编号i的节点力的x分量以及节点的法向和切向接触力的x分量,均可以直接从ABAQUS的输出文件提取得到。截面的工程应力可以由下式计算为:

      $ \sigma = \sigma _x^{\left( {\rm{n}} \right)} + \sigma _x^{\left( {\rm{c}} \right)} $

      式中:σx(n)=Fx(n)/A0,表示由节点所传递的应力;σx(c)=Fx(c)/A0,表示由胞壁接触所引发的应力。其中A0为截面的初始面积。

    • 冲击速度分别为200和50 m/s时的应力分布如图 2所示。在任意名义应变时刻下,应力会在某个位置发生陡然减小,并且这个位置随着名义应变的增大往试件右侧移动。这说明在高速和中速撞击下试件内均存在冲击波的传播,应力发生急剧变化的位置即为冲击波波阵面的位置,中速撞击下的波阵面区域相较于高速撞击更加平缓。从图 2可以明显地观察到应力增强现象,跨过冲击波阵面,应力由波前应力迅速增加到波后应力并几乎保持不变。

      图  2  两种冲击速度下的应力分布图

      Figure 2.  One-dimensional stress distributions at two impact velocities

    • 可以由应力分布梯度来确定冲击波波阵面Φ的位置,如图 3所示。应力分布梯度反映了一维应力在加载方向上的变化,在应力梯度绝对值达到最大的拉格朗日坐标位置,应力变化最快,这个位置即为当前时刻冲击波波阵面的位置。中速和高速压溃下波阵面位置均可以采用此方法得到。不同冲击速度下得到的波阵面位置随冲击时间变化关系如图 4所示。在恒速压缩情形下,冲击波的波阵面位置与冲击时间关系为线性关系。可通过对线性关系拟合得到的斜率来预估冲击波速度。从图 4中可以看出,冲击速度越大,斜率越大,则相应的冲击波速度也越大。

      图  3  一维的应力分布及其应力梯度分布

      Figure 3.  One-dimensional stress distributions and the corresponding stress gradients

      图  4  不同冲击速度下冲击波位置随时间的关系

      Figure 4.  Variation of shock front position with impact time at different impact velocities

    • 除了由有限元结果直接得到的应力信息来得到冲击波速度之外,一维冲击波理论也可以对蜂窝中的冲击波传播行为进行描述。由应力波理论[17],跨波阵面的质量和动量守恒关系分别为:

      $ {v_{\rm{B}}}\left( t \right) - {v_{\rm{A}}}\left( t \right) = {v_{\rm{s}}}\left( t \right)\left( {{\varepsilon _{\rm{B}}}\left( t \right) - {\varepsilon _{\rm{A}}}\left( t \right)} \right) $

      $ {\sigma _{\rm{B}}}\left( t \right) - {\sigma _{\rm{A}}}\left( t \right) = {\rho _0}{v_{\rm{s}}}\left( t \right)\left( {{v_{\rm{B}}}\left( t \right) - {v_{\rm{A}}}\left( t \right)} \right) $

      式中:vs(t)是冲击波速度, σBεBvB为波后方的应力、应变和速度,σAεAvA为波前方的应力、应变和速度。联立这两个守恒关系可以得到冲击波速度与波前波后应力、应变的关系为:

      $ {v_{\rm{s}}}\left( t \right) = \sqrt {\frac{{{\sigma _{\rm{B}}}\left( t \right) - {\sigma _{\rm{A}}}\left( t \right)}}{{{\rho _0}\left( {{\varepsilon _{\rm{B}}}\left( t \right) - {\varepsilon _{\rm{A}}}\left( t \right)} \right)}}} $

      如果能直接从有限元结果中得到波前和波后应力/应变信息,就可以由式(5)计算得到冲击波速度的大小。为了简便,这里将这种结合冲击波理论和数值模拟结果的方法叫做一维冲击波理论方法。

      通过对一维应力和应变分布求平均值,可以得到波前和波后的应力和应变,如图 5所示。图中一维应变分布采用文献[8]提出的局部应变计算方法得到。在求平均的过程中,为了消除冲击波波阵面平均效应的影响,将紧邻波阵面附近1.5倍胞元半径Rc内的数据排除掉。波后应力和应变信息由冲击端到Φ-1.5Rc区域内平均得到,波前应力应变信息由Φ+1.5Rc到支撑端区域平均得到。再由式(5)就可以得到不同冲击速度下的冲击波速度。

      图  5  一维应力和应变分布及其理想化

      Figure 5.  One-dimensional stress and strain distributions and the idealizations

    • 通过截面应力和局部应变场等有限元方法以及一维冲击波理论方法都可以得到冲击波速度。此外,也可由已有的冲击波模型得到冲击波速度的大小。由R-PP-L模型[2]得到的冲击波速度为:

      $ {v_{\rm{s}}} = v/{\varepsilon _{\rm{L}}} $

      式中:εL为多胞材料的锁定应变。由R-PH模型[13]得到的冲击波速度为:

      $ {v_{\rm{s}}} = v + \sqrt {C/{\rho _0}} $

      式中:C为多胞材料的塑性硬化参数。在本文中,εL=0.64,C=0.226 MPa。

      由不同方法得到的冲击波速度如图 6所示。由一维冲击波理论方法以及R-PH冲击波模型得到的冲击波速度比较接近有限元方法得到的冲击波速度,并且与冲击速度的差值均为一个常值,即材料冲击参数。然而由R-PP-L模型预估得到的冲击波速度明显高于有限元结果。这是由于R-PP-L模型采用的是锁定不变的压实应变,实际上其值小于动态压溃的压实应变。由质量守恒关系可以发现,压实应变越小,冲击波速度越大。R-PH模型中的压实应变是随冲击速度变化的动态应变,因此更加适合描述多胞材料在动态压溃下的冲击波传播。

      图  6  不同方法得到的冲击波速度的比较

      Figure 6.  Comparison of shock wave speeds obtained using different methods

    • 前文中通过应力分布得到了高速冲击下的冲击波速度与冲击速度的关系。在中等冲击速度下,仍然可以观察到冲击波传播的现象[4, 8],因此仍然采用应力梯度分布方法来确定冲击波速度的大小,结果见图 7。当冲击速度很高的时候,冲击波速度与冲击速度呈线性关系并且两者之间的差值几乎为一个常数[13]。但随着冲击速度的减小,冲击波速度趋近于常数。因此,冲击波速度对冲击速度的一阶导数在低速和高速下分别趋于0和1,从低速到高速的完整曲线具有明显的S型特征。利用S函数(sigmoid function),冲击波速度对冲击速度的一阶导数可以写作:

      $ \frac{{{\rm{d}}{v_{\rm{s}}}}}{{{\rm{d}}v}} = \frac{1}{{1 + {{\rm{e}}^{\frac{{b - v}}{a}}}}} $

      图  7  冲击波速度与冲击速度的关系

      Figure 7.  Variations of the shock wave speed with the impact velocity

      式中:ab为待定参数。上式的积分形式可以写作:

      $ {v_{\rm{s}}} = a\ln \left( {1 + {{\rm{e}}^{\frac{{v - b}}{a}}}} \right) + c $

      式中:c为待定参数。可用式(9)对有限元结果进行拟合,从图 7可见拟合效果较好,其中拟合参数为a=3.603 m/s, b=40.68 m/s和c=54.91 m/s。若式(8)右边的值较小,对应于低速冲击,此时并不能观察到冲击波传播的现象,因此不妨取该值为0.05,可以估算从均匀模式向过渡模式转变的临界速度vc1=baln19=30.07 m/s。若式(8)右边的值较大,对应于高速冲击,不妨取该值为0.95,可以估算从过渡模式向冲击模式转变的临界速度vc2=b+aln19=51.29 m/s。可见参数ab的线性组合反映了模式转变速度的大小。当0<vb时,冲击波速度vs近似等于c。当vb时,$1 + {{\rm{e}}^{\frac{{v - b}}{a}}} \approx {{\rm{e}}^{\frac{{v - b}}{a}}}$,此时冲击波速度vs=v-b+c,和R-PH模型中得到的冲击波速度与冲击速度变化关系一致,即c1=c-b。因此,式(9)表征了冲击波速度与冲击速度的一致近似关系。

    • 在本文中考察的恒速冲击情形中,若取εA=0, vA=0, vB=vσA=σ0,由跨波阵面的守恒关系式(3)和(4)有:

      $ \left\{ \begin{array}{l} {\varepsilon _{\rm{B}}} = v/{v_{\rm{s}}}\\ {\sigma _{\rm{B}}} = {\sigma _0} + {\rho _0}v{v_{\rm{s}}} \end{array} \right. $

      将式(9)代入上式,可得:

      $ \left\{ \begin{array}{l} {\varepsilon _{\rm{B}}} = v/\left[ {a\ln \left( {1 + {{\rm{e}}^{\frac{{v - b}}{a}}}} \right) + c} \right]\\ {\sigma _{\rm{B}}} = {\sigma _0} + {\rho _0}v\left[ {a\ln \left( {1 + {{\rm{e}}^{\frac{{v - b}}{a}}}} \right) + c} \right] \end{array} \right. $

      与有限元结果的比较如图 8所示,图中虚线对应于式(11)。可见式(11)与有限元结果存在一定的误差。事实上,在恒速冲击下,冲击波波阵面前方的应变和粒子速度不能直接忽略。

      图 9给出了不同冲击速度下试件中的速度分布情况,局部粒子速度在冲击波波阵面位置附近约两个胞元区间内急剧变小。当冲击速度为中等大小时,波阵面位置处的粒子速度变化更为平缓。整个波前区域处于低速状态,近似呈线性分布,与冲击速度的大小关系不大。冲击波前方的粒子速度和变形不能忽略,可以通过图 9估算vA≈10 m/s,初始压溃应变εA取准静态实验中常用的ε0 =0.02。因此,由跨波阵面的守恒关系式(3)和(4)有:

      $ \left\{ \begin{array}{l} {\varepsilon _{\rm{B}}} = {\varepsilon _{\rm{0}}} + \left( {v - {v_{\rm{A}}}} \right)/{v_{\rm{s}}}\\ {\sigma _{\rm{B}}} = {\sigma _0} + {\rho _0}\left( {v - {v_{\rm{A}}}} \right)/{v_{\rm{s}}} \end{array} \right. $

      图  8  波后应变和波后应力随冲击速度的变化

      Figure 8.  Variations of strain and stress behind shock front with impact velocity

      图  9  不同冲击速度下的速度分布

      Figure 9.  Velocity distributions at different impact velocities

      将式(9)代入上式,可得:

      $ \left\{ \begin{array}{l} {\varepsilon _{\rm{B}}} = {\varepsilon _{\rm{0}}} + \left( {v - {v_{\rm{A}}}} \right)/\left[ {a\ln \left( {1 + {{\rm{e}}^{\frac{{v - b}}{a}}}} \right) + c} \right]\\ {\sigma _{\rm{B}}} = {\sigma _0} + {\rho _0}\left( {v - {v_{\rm{A}}}} \right)/\left[ {a\ln \left( {1 + {{\rm{e}}^{\frac{{v - b}}{a}}}} \right) + c} \right] \end{array} \right. $

      图 8所示,这个结果与有限元结果吻合较好。引入无量纲参数$\xi = {{\rm{e}}^{\frac{{v - b}}{a}}}$,式(13)可以进一步写作:

      $ \left\{ \begin{array}{l} {\varepsilon _{\rm{B}}} = {\varepsilon _{\rm{0}}} + \frac{{\ln \xi + \beta - \alpha }}{{\ln \left( {1 + \xi } \right) + \gamma }}\\ {\sigma _{\rm{B}}} = {\sigma _0} + B\left( {\ln \xi + \beta - \alpha } \right)\left[ {\ln \left( {1 + \xi } \right) + \gamma } \right] \end{array} \right. $

      式中:α=vA/a, β =b/a, γ=c/aB=ρ0a2。因此,式(14)给出恒速冲击下动态应力应变状态的一致近似关系。准静态和动态应力应变关系的比较如图 10所示。在图 10中,准静态应力应变曲线及动态应力应变状态点由有限元计算得到,动态应力应变曲线由式(14)得到。在动态应力应变关系的曲线上的每个点都对应着特定的一个冲击速度。在高应变情形下,动态应力应变关系位于准静态曲线右侧,这与已有文献中的结论一致。随着应变的减小,动态的应力应变关系逼近准静态应力应变曲线。

      图  10  准静态及动态应力应变关系

      Figure 10.  Quasi-static and dynamic stress-strain relations

    • 采用细观有限元模型研究了随机蜂窝结构在恒速冲击下的动态性能,采用截面应力计算方法得到了随机蜂窝试件内的应力分布和冲击波速度,比较了高速冲击下由不同方法得到的冲击波速度与冲击速度的关系。结果表明,R-PP-L模型高估了冲击波速度,但R-PH模型以及一维冲击波理论得到的冲击波速度与有限元模拟得到的结果接近。在冲击速度很高的时候,冲击波速度与冲击速度的关系趋于线性,但随着冲击速度的减小,冲击波速度不断减少并趋于常数。最后,发展了可以表征冲击波速度与冲击速度的关系的一致近似模型,并基于一维冲击波理论发展了动态应力应变关系的一致近似模型。这个模型可以对多胞材料在过渡模式及冲击模式下的动态压溃行为进行表征。

参考文献 (17)

目录

    /

    返回文章
    返回