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

具有恒定冲击载荷的梯度泡沫金属材料设计

常白雪 郑志军 赵凯 何思渊 虞吉林

引用本文:
Citation:

具有恒定冲击载荷的梯度泡沫金属材料设计

    作者简介: 常白雪(1992- ),女,博士研究生,bxchang@mail.ustc.edu.cn;
    通讯作者: 郑志军, zjzheng@ustc.edu.cn
  • 中图分类号: O347

Design of gradient foam metal materials with a constant impact load

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

  • 摘要: 多胞材料可通过大变形大量地吸收冲击能量,引入密度梯度可进一步提高其耐撞性。梯度多胞材料的宏观力学响应对材料密度分布极为敏感,不同类型的细观构型的影响也极为不同。已有的研究工作主要局限在对给定的密度梯度分析其动态响应,较少对耐撞性设计方法进行研究。本文针对梯度闭孔泡沫金属材料,基于非线性塑性冲击波模型发展了耐撞性反向设计方法,以维持冲击物受载恒定为目标,运用级数法获得了简化模型和渐近解。利用变胞元尺寸法构建了连续梯度变化的三维Voronoi细观有限元模型,并利用ABAQUS/Explicit有限元软件对理论设计进行数值验证。结果表明,反向设计理论简化模型的渐近解对于梯度闭孔泡沫金属材料的耐撞性设计是有效的,所提出的耐撞性设计方法在控制冲击吸能过程和冲击物受载方面具有指导意义。
  • 图 1(a)  不同相对密度下闭孔泡沫模型的准静态应力应变关系

    Figure 1(a).  Nominal stress–strain relations of closed-cell foam models with different relative densities

    图 1(b)  无量纲初始压溃应力和应变硬化参数与多胞模型相对密度的幂律关系[13]

    Figure 1(b).  Power-law fitting of material parameters of the R-PH idealization with relative density[13]

    图 2  质量块冲击梯度多胞杆

    Figure 2.  A diagram of 3D graded Voronoi models under mass impact

    图 3  相对密度分布的渐近解和冲击载荷历程的理论预测值

    Figure 3.  Theoretical predictions of asymptotic solutions of relative density distribution and the history curves of impact force

    图 4  冲击速度与冲击波波阵面位置的历史曲线

    Figure 4.  Evolution history of impact velocity and location of shock front

    图 5  密度梯度多胞杆细观有限元模型

    Figure 5.  Cell-based finite element models of density gradient cellular rods

    图 6  梯度多胞细观有限元模型中轴剖面变形图

    Figure 6.  Deformation patterns of cell-based finite element models in the middle section of the density gradient cellular rods perpendicular to x-axis

    图 7  冲击端和支撑端载荷理论预测与有限元计算结果的比较

    Figure 7.  Comparisons of impact force and support force curves history between theoretical predictions and finite element (FE) results

    图 8  有限元计算结果对比

    Figure 8.  Comparisons of finite element (FE) results

  • [1] 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.
    [2] 华云龙, 余同希. 多胞材料的力学行为 [J]. 力学进展, 1990, 21(4): 457–469. DOI: 10.6052/1000-0992-1991-4-J1991-052
    HUA Yunlong, YU Tongxi. Mechanical behavior of cellular solids [J]. Advances in Mechanics, 1990, 21(4): 457–469. DOI: 10.6052/1000-0992-1991-4-J1991-052
    [3] ELNASRI I, PATTOFATTO S, ZHAO H, et al. Shock enhancement of cellular structures under impact loading: part Ⅰ experiments [J]. Journal of the Mechanics and Physics of Solids, 2007, 55(12): 2652–2671. DOI: 10.1016/j.jmps.2007.04.005.
    [4] GIBSON L J, ASHBY M F. Cellular solids: structure and properties [M]. Cambridge: Cambridge University Press, 1999.
    [5] CUI L, KIERNAN S, GILCHRIST M D. Designing the energy absorption capacity of functionally graded foam materials [J]. Materials Science and Engineering: A, 2009, 507(1): 215–225. DOI: 10.1016/j.msea.2008.12.011.
    [6] WANG X K, ZHENG Z J, YU J L, et al. Impact resistance and energy absorption of functionally graded cellular structures [J]. Applied Mechanics and Materials, 2011, 69: 73−78. DOI: 10.4028/www.scientific.net/AMM.69.73.
    [7] WANG X K, ZHENG Z J, YU J L. Crashworthiness design of density-graded cellular metals [J]. Theoretical and Applied Mechanics Letters, 2013, 3(3): 031001. DOI: 10.1063/2.1303101.
    [8] SHEN C J, YU T X, LU G X. Double shock mode in graded cellular rod under impact [J]. International Journal of Solids and Structures, 2013, 50(1): 217–233. DOI: 10.1016/j.ijsolstr.2012.09.021.
    [9] ZHENG J, QIN Q, WANG T J. Impact plastic crushing and design of density-graded cellular materials [J]. Mechanics of Materials, 2016, 94: 66–78. DOI: 10.1016/j.mechmat.2015.11.014.
    [10] SHEN C J, LU G X, YU T X. Investigation into the behavior of a graded cellular rod under impact [J]. International Journal of Impact Engineering, 2014, 74: 92–106. DOI: 10.1016/j.ijimpeng.2014.02.015.
    [11] 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.
    [12] 蔡正宇, 丁圆圆, 王士龙, 等. 梯度多胞牺牲层的抗爆炸分析 [J]. 爆炸与冲击, 2017, 37(3): 396–404. DOI: 10.11883/1001-1455(2017)03-0396-09
    CAI Zhengyu, DING Yuanyuan, WANG Shilong, et al. Anti-blast analysis of graded cellular sacrificial cladding [J]. Explosion and Shock Waves, 2017, 37(3): 396–404. DOI: 10.11883/1001-1455(2017)03-0396-09
    [13] YANG J, WANG S L, DING Y Y, et al. Crashworthiness of graded cellular materials: a design strategy based on a nonlinear plastic shock model [J]. Materials Science and Engineering: A, 2017, 680: 411–420. DOI: 10.1016/j.msea.2016.11.010.
    [14] DING Y Y, WANG S L, ZHAO K, et al. Blast alleviation of cellular sacrificial cladding: a nonlinear plastic shock model [J]. International Journal of Applied Mechanics, 2016, 8(4): 1650057. DOI: 10.1142/S1758825116500575.
    [15] LI K, GAO X L, WANG J. Dynamic crushing behavior of honeycomb structures with irregular cell shapes and non-uniform cell wall thickness [J]. International Journal of Solids and Structures, 2007, 44(14−15): 5003–5026. DOI: 10.1016/j.ijsolstr.2006.12.017.
    [16] AJDARI A, NAYEB-HASHEMI H, VAZIRI A. Dynamic crushing and energy absorption of regular, irregular and functionally graded cellular structures [J]. International Journal of Solids and Structures, 2011, 48(3−4): 506–516. DOI: 10.1016/j.ijsolstr.2010.10.018.
    [17] 张新春, 刘颖. 密度梯度蜂窝材料动力学性能研究 [J]. 工程力学, 2012, 29(8): 372–377. DOI: 10.6052/j.issn.1000-4750.2010.12.0872
    ZHANG Xinchun, LIU Ying. Research on the dynamic crushing of honeycombs with density gradient [J]. Engineering Mechanics, 2012, 29(8): 372–377. DOI: 10.6052/j.issn.1000-4750.2010.12.0872
    [18] 吴鹤翔, 刘颖. 梯度变化对密度梯度蜂窝材料力学性能的影响 [J]. 爆炸与冲击, 2013, 33(2): 163–168. DOI: 10.3969/j.issn.1001-1455.2013.02.008
    WU Hexiang, LIU Ying. Influences of density gradient variation on mechanical performances of density-graded honeycomb materials [J]. Explosion and Shock Waves, 2013, 33(2): 163–168. DOI: 10.3969/j.issn.1001-1455.2013.02.008
    [19] FAN J H, ZHANG J J, WANG Z H, et al. Dynamic crushing behavior of random and functionally graded metal hollow sphere foams [J]. Materials Science and Engineering: A, 2013, 561: 352–361. DOI: 10.1016/j.msea.2012.10.026.
    [20] ZHANG J J, WANG Z H, ZHAO L M. Dynamic response of functionally graded cellular materials based on the Voronoi model [J]. Composites Part B: Engineering, 2016, 85: 176–187. DOI: 10.1016/j.compositesb.2015.09.045.
    [21] CHEN D, KITIPORNCHAI S, YANG J. Dynamic response and energy absorption of functionally graded porous structures [J]. Materials & Design, 2018, 140: 473–487. DOI: 10.1016/j.matdes.2017.12.019.
    [22] 常白雪, 郑志军, 赵凯, 等. 梯度多胞材料耐撞性设计的简化模型和渐近解 [J]. 中国科学: 物理学 力学 天文学, 2018, 48(9): 094615. DOI: 10.1360/SSPMA2018-00162
    CHANG Baixue, ZHENG Zhijun, ZHAO Kai, et al. A simplified model and its asymptotic solution for the crashworthiness design of graded cellular material [J]. Scientia Sinica: Physica, Mechanica & Astronomica, 2018, 48(9): 094615. DOI: 10.1360/SSPMA2018-00162
    [23] DING Y Y, WANG S L, ZHENG Z J, et al. Dynamic crushing of cellular materials: a unique dynamic stress–strain state curve [J]. Mechanics of Materials, 2016, 100: 219–31. DOI: 10.1016/j.mechmat.2016.07.001.
    [24] WANG P, ZHENG Z J, LIAO S F, et al. Strain-rate effect on initial crush stress of irregular honeycomb under dynamic loading and its deformation mechanism [J]. Acta Mechanica Sinica, 2018, 34(1): 117–129. DOI: 10.1007/s10409-017-0716-1.
    [25] HE S Y, ZHANG Y, DAI G, et al. Preparation of density-graded aluminum foam [J]. Materials Science and Engineering A, 2014, 618: 496–499. DOI: 10.1016/j.msea.2014.08.087.
  • [1] 王长峰郑志军虞吉林 . 泡沫杆撞击刚性壁的动态压溃模型. 爆炸与冲击, 2013, 33(6): 587-593. doi: 10.11883/1001-1455(2013)06-0587-07
    [2] 李松晏郑志军虞吉林 . 高速列车吸能结构设计和耐撞性分析. 爆炸与冲击, 2015, 35(2): 164-170. doi: 10.11883/1001-1455(2015)02-0164-07
    [3] 徐峰祥张锁武昆迎 . 厚度幂指数分布管状结构耐撞性设计准则与方法研究. 爆炸与冲击, 2019, 39(3): 035103-1-035103-9. doi: 10.11883/bzycj-2018-0013
    [4] 张健赵桂平卢天健 . 泡沫金属在冲击载荷下的动态压缩行为. 爆炸与冲击, 2014, 34(3): 278-284. doi: 10.11883/1001-1455(2014)03-0278-07
    [5] 宋延泽李志强赵隆茂 . 基于十四面体模型的闭孔泡沫材料动态力学性能的有限元分析. 爆炸与冲击, 2009, 29(1): 49-55. doi: 10.11883/1001-1455(2009)01-0049-07
    [6] 刘晶波姚玲 . 有限元畸变单元离散模型中的SH波动. 爆炸与冲击, 1993, 13(1): 9-17.
    [7] 陈二云赵改平戴韧马大为 . 可压缩多介质流动的间断有限元方法. 爆炸与冲击, 2010, 30(4): 419-423. doi: 10.11883/1001-1455(2010)04-0419-05
    [8] 宋顺成曹法和田时雨刘心德陈兆生 . 高速弹丸撞击地质材料的有限元计算. 爆炸与冲击, 1994, 14(2): 161-168.
    [9] 任波王乘 . 并行冲击/侵彻有限元数值模拟技术. 爆炸与冲击, 2005, 25(3): 260-264. doi: 10.11883/1001-1455(2005)03-0260-05
    [10] 谢宗蕻卞文杰昂海松乔新 . 峰窝夹芯结构雷达罩鸟撞有限元分析与模拟. 爆炸与冲击, 1999, 19(3): 235-242.
    [11] 周杰陶钢潘保青张洪伟 . 爆炸冲击波对人体胸部创伤机理的有限元方法研究. 爆炸与冲击, 2013, 33(3): 315-321. doi: 10.11883/1001-1455(2013)03-0315-06
    [12] 殷之平李玉龙黄其青 . 含诱导缺陷薄壁圆管耐撞性优化设计. 爆炸与冲击, 2011, 31(4): 418-422. doi: 10.11883/1001-1455(2011)04-0418-05
    [13] 钱七虎黄小平唐德高李志成 . 冲击荷载作用下有限元方法求解波动过程的精度研究. 爆炸与冲击, 1995, 15(1): 1-10.
    [14] 黄晨光段祝平 . 含热传导的冲击动力学有限元程序的研究与应用. 爆炸与冲击, 2001, 21(4): 253-258.
    [15] 本文用有限元方法分析了一定初缺陷杆受轴向冲击的弹塑性动态屈曲。由变形功相关的屈曲判据求出屈曲时间')">揭敏本文用有限元方法分析了一定初缺陷杆受轴向冲击的弹塑性动态屈曲。由变形功相关的屈曲判据求出屈曲时间计算了初缺陷及冲击载荷形状和大小对屈曲时间的影响。 . 一定初缺陷杆在轴向冲击下弹塑性动态屈曲有限元计算. 爆炸与冲击, 1991, 11(2): 153-160.
    [16] 黄延恩徐少文王乘 . 冲击载荷作用下结构动态响应分析的罚函数有限元法. 爆炸与冲击, 1998, 18(3): 226-230.
    [17] 朱秀云潘蓉林皋李亮 . 基于ANSYS/LS-DYNA的钢板混凝土墙冲击实验的有限元分析. 爆炸与冲击, 2015, 35(2): 222-228. doi: 10.11883/1001-1455(2015)02-0222-07
    [18] 刘海峰宁建国 . 冲击荷载作用下混凝土材料的细观本构模型. 爆炸与冲击, 2009, 29(3): 261-267. doi: 10.11883/1001-1455(2009)03-0261-07
    [19] 韩小平张元冲沈亚鹏张泰华赵壮华 . 冲击载荷下炸药装药动态响应的有限元分析及热点形成机理的数值模拟. 爆炸与冲击, 1997, 17(2): 143-152.
    [20] 王吉王肖钧卞梁 . 光滑粒子法与有限元的耦合算法及其在冲击动力学中的应用. 爆炸与冲击, 2007, 27(6): 522-528. doi: 10.11883/1001-1455(2007)06-0522-07
  • 加载中
图(9)
计量
  • 文章访问数:  192
  • HTML全文浏览量:  608
  • PDF下载量:  62
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-10-31
  • 录用日期:  2018-12-07
  • 网络出版日期:  2019-04-25
  • 刊出日期:  2019-04-01

具有恒定冲击载荷的梯度泡沫金属材料设计

    作者简介:常白雪(1992- ),女,博士研究生,bxchang@mail.ustc.edu.cn
    通讯作者: 郑志军, zjzheng@ustc.edu.cn
  • 1. 中国科学技术大学近代力学系中国科学院材料力学行为和设计重点实验室,安徽 合肥 230026
  • 2. 西安交通大学航天学院机械结构强度与振动国家重点实验室,陕西 西安 710049
  • 3. 东南大学生物科学与医学工程学院生物电子学国家重点实验室,江苏 南京 210096

摘要: 多胞材料可通过大变形大量地吸收冲击能量,引入密度梯度可进一步提高其耐撞性。梯度多胞材料的宏观力学响应对材料密度分布极为敏感,不同类型的细观构型的影响也极为不同。已有的研究工作主要局限在对给定的密度梯度分析其动态响应,较少对耐撞性设计方法进行研究。本文针对梯度闭孔泡沫金属材料,基于非线性塑性冲击波模型发展了耐撞性反向设计方法,以维持冲击物受载恒定为目标,运用级数法获得了简化模型和渐近解。利用变胞元尺寸法构建了连续梯度变化的三维Voronoi细观有限元模型,并利用ABAQUS/Explicit有限元软件对理论设计进行数值验证。结果表明,反向设计理论简化模型的渐近解对于梯度闭孔泡沫金属材料的耐撞性设计是有效的,所提出的耐撞性设计方法在控制冲击吸能过程和冲击物受载方面具有指导意义。

English Abstract

    • 多胞材料具有高孔隙率、可变形量大和轻量化特征以及超强的能量吸收能力,广泛地应用于汽车、航空航天和铁路等领域。在动态压缩情形下,多胞材料具有明显的变形局部化特征和应力增强现象[1-3],可以很好地提高结构的耐撞性。如何优化设计多胞材料以满足结构耐撞性要求,是近几年结构冲击领域的研究热点。

      多胞材料的准静态名义应力应变关系呈现明显的三段式特征:弹性段、长而稳定的平台段和压实段。基于此特征,Reid等[1]提出了刚性–塑性–锁定(R-PP-L)模型,并发展了一维冲击波模型以表征多胞材料的动态力学行为。Gibson等[4]定量化给出了R-PP-L模型中多胞材料的平台应力和锁定应变与相对密度和基体材料屈服应力的关系。该理论已被应用于研究不同相对密度和强度分布下梯度多胞材料的吸能能力和峰值应力[5-10]。然而,R-PP-L模型只满足对多胞材料力学性能的一阶近似,Zheng等[11]提出了更为精确的刚性–塑性硬化(R-PH)模型,其应力应变关系写作

      $ \sigma ={\sigma _0} + \frac{{C\varepsilon }}{{{{(1 - \varepsilon )}^2}}} $

      式中: σ0C分别为初始压溃应力和应变硬化参数,这两个材料参数与多胞材料的相对密度相关[12-13]。Ding等[14]发现基于R-PH模型可以更好地指导爆炸牺牲层的设计。

      随着更为精确的多胞材料本构模型的提出,梯度多胞材料的动态变形机理的研究也趋于完善。但多数研究停留在讨论不同梯度多胞材料的动态力学性能上,分析具有不同密度和强度分布的多胞模型的动态力学性能[15-21]。只通过观察比较其力学响应,选择满足耐撞性的要求的梯度模型,无法直接达到最优化的目的。Yang等[13]首次提出了梯度多胞材料的耐撞性反向设计理论,直接设定物体冲击力要求,反向设计多胞材料的密度分布,并利用三维细观有限元模拟验证其有效性。不过,该耐撞性反向设计理论无显式解,Yang等[13]运用四阶龙格库塔法获得了数值解。常白雪等[22]利用该理论对二维Voronoi模型进行耐撞性反向设计,并提出了该理论的简单而方便的渐近解。由于二维随机蜂窝与三维闭孔胞元的结构差异,针对于二维模型提出的渐近解,无法适用于三维多胞构型的耐撞性设计需求,而对于闭孔多胞构型的耐撞性设计理论的近似解尚未提出。

      本文中进一步研究了质量块以特定的初速度撞击梯度闭孔泡沫金属杆的理论模型和有限元模拟:通过讨论物体冲击过程中受载恒定的耐撞性要求,对梯度多胞材料耐撞性反向设计理论进行简化。针对闭孔泡沫金属,拟合R-PH模型中初始压溃应力和应变硬化参数与相对密度的关系;运用级数法,求解耐撞性反向设计理论,获得梯度多胞材料的密度分布;利用变胞元尺寸法和三维随机Voronoi技术,生成特定密度分布的梯度闭孔多胞模型;利用有限元软件ABAQUS/Explicit模块计算,获取冲击力结果验证设计的有效性,并与均匀密度分布的多胞材料结果进行对比,讨论合理的密度分布设计的必要性。

    • 本文中采用的均匀密度的闭孔泡沫模型与文献[11]中的一致。在10 m/s的恒速压缩下,利用ABAQUS/Explicit模块进行有限元模拟,获取了不同相对密度的Voronoi模型准静态压缩下的力学性能,其名义应力应变曲线如图1(a)所示。进而,利用Zheng等 [11]提出的R-PH模型对准静态名义应力应变曲线进行最小二乘法拟合,如图1(a)所示。所得的初始压溃应力和应变硬化参数与相对密度呈幂律关系[13],如图1(b)所示。为方便理论分析,本文采用幂次为3/2的幂律关系:

      图  1(a)  不同相对密度下闭孔泡沫模型的准静态应力应变关系

      Figure 1(a).  Nominal stress–strain relations of closed-cell foam models with different relative densities

      图  1(b)  无量纲初始压溃应力和应变硬化参数与多胞模型相对密度的幂律关系[13]

      Figure 1(b).  Power-law fitting of material parameters of the R-PH idealization with relative density[13]

      $ \left\{ \begin{gathered} {\sigma _0}(\rho )={k_1}{\sigma _{{\rm{ys}}}}{\rho ^{3/2}} \hfill \\ C(\rho )={k_2}{\sigma _{{\rm{ys}}}}{\rho ^{3/2}} \hfill \\ \end{gathered} \right. $

      式中:ρ为相对密度,σys为基体材料的屈服应力,k1k2为拟合参数,分别为1.166和0.115。

    • 在质量块初速度冲击梯度泡沫金属材料的过程中,为保护物体安全,且平稳减速,须控制物体在冲击过程中受冲击载荷恒定且低于物体耐受载荷。因而,梯度泡沫金属杆需根据耐撞性要求进行特定的相对密度分布设计。针对二维随机蜂窝模型的耐撞性反向设计中,常白雪等[22]提出了控制质量块冲击梯度多胞金属杆过程中受载恒定的相对密度分布渐近解,通过二维细观有限元模拟验证了设计的有效性,并指出特定的相对密度分布只适用于特定的冲击情形。对于闭孔梯度泡沫材料,Yang等[13]提出了耐撞性反向设计理论,通过四阶龙格库塔法,针对于控制冲击过程中质量块受载恒定、线性递增和线性递减的耐撞性要求,获得了相对密度分布的数值解,但未给出可方便实际工程应用的公式解。本文中针对物体受载恒定的耐撞性要求,对耐撞性反向设计理论进行简化,运用级数法,求得相对密度分布的渐近解。

      考虑一质量为M的物块以冲击速度v0撞击横截面面积为A0的闭孔泡沫金属杆,被保护质量块在整个冲击过程中受载恒定,额定受载力为F0。冲击过程中质量块的速度为v=v0a0t,其中a0 =F0/M为冲击物体的加速度,t为冲击时间。冲击的总时间T0 =Mv0/F0,特征长度L0=Mv02/F0。根据Yang等[13]的耐撞性反向设计理论,当冲击载荷为恒定值时,该理论可以简化为

      $ \left\{ \begin{gathered} \frac{{\rm{d} \rho }}{{\rm{d}} t}=\frac{{{a_0}{\rho _{\rm{s}}}\rho \left( {3v + 2c(\rho )} \right)}}{{{{\sigma_0 '}}(\rho ) + {\rho _{\rm{s}}}{v^2} + {\rho _{\rm{s}}}v\left( {c(\rho ) + \rho c'(\rho )} \right)}} \hfill \\ \frac{{\rm{d} \varPhi }}{{\rm{d}} t}=v + c(\rho ) \hfill \end{gathered} \right. $

      式中: ρ为泡沫金属材料相对密度,ρs为基体材料密度,Φ为冲击波阵面位置,撇号表示对ρ的导数,$c(\rho )=\sqrt {{{C(\rho )} \mathord{\left/ {\vphantom {{C(\rho )} {{\rho _{\rm{s}}}\rho }}} \right.} {{\rho _{\rm{s}}}\rho }}} ={c_0}{\rho ^{1/4}}$${c_0}{\rm{=}}\sqrt {{k_2}{\sigma _{{\rm{ys}}}}/{\rho _{\rm{s}}}} $。考虑到相对密度、长度、冲击速度和时间的特征量,定义无量纲量:

      $ \bar \rho =\frac{\rho }{{{\rho _0}}},\quad \bar \varPhi =\frac{\varPhi }{{{L_0}}},\quad \bar v=\frac{v}{{{v_0}}},\quad \bar t=\frac{t}{{{T_0}}} $

      式中: ρ0为冲击端的相对密度。此时,式(3)可以进一步写作:

      $ \left\{ \begin{gathered} \frac{{\rm{d} \bar \rho }}{{\rm{d} \bar \varPhi }}=\frac{{4\bar \rho \left( {3\bar v + 2{{\bar c}_0}{{\bar \rho }^{1/4}}} \right)}}{{\left( {6{{\bar s}_0}{{\bar \rho }^{1/2}} + 4{{\bar v}^2} + 5\bar v{{\bar c}_0}{{\bar \rho }^{1/4}}} \right)\left( {\bar v + {{\bar c}_0}{{\bar \rho }^{1/4}}} \right)}} \hfill \\ \frac{{{\rm{d}} \bar v}}{{\rm{d} \bar \varPhi }}= - {\left( {\bar v + {{\bar c}_0}{{\bar \rho }^{1/4}}} \right)^{ - 1}} \hfill \\ \end{gathered} \right. $

      式中: $\bar v=1 - \bar t$${\bar s_0}={s_0}{\rho _0}^{1/2}/({\rho _{\rm{s}}}v_0^2)$${\bar c_0}={c_0}{\rho _0}^{1/4}/{v_{\rm{0}}}$s0=k1σys。利用级数法,可以求得渐近解:

      $ \left\{ \begin{gathered} \bar \rho (\bar \varPhi )=1 + {a_1}\bar \varPhi + {a_2}{{\bar \varPhi }^2} + O({{\bar \varPhi }^3}) \\ \bar v(\bar \varPhi )=1 + {b_1}\bar \varPhi + {b_2}{{\bar \varPhi }^2} + O({{\bar \varPhi }^3}) \\ \end{gathered} \right. $

      式中:一阶系数和二阶系数分别为:

      $ \left\{ \begin{gathered} {a_1}=\frac{{12 + 8{{\bar c}_0}}}{{\left( {1 + {{\bar c}_0}} \right)\left( {4 + 5{{\bar c}_0} + 6{{\bar s}_0}} \right)}} \hfill \\ {b_1}= - \frac{1}{{1 + {{\bar c}_0}}} \hfill \\ \end{gathered} \right. $

      $ \left\{ \begin{aligned} &{a_2}=\frac{{2\left( {1 + {{\bar c}_0}} \right)\left( {240 + 543{{\bar c}_0} + 404\bar c_0^2 + 110\bar c_0^3} \right) + 12\times\left( {42 + 86{{\bar c}_0} + 61\bar c_0^2 + 18\bar c_0^3 - 6{{\bar c}_0}{{\bar s}_0}} \right){{\bar s}_0}}}{{{{\left( {1 + {{\bar c}_0}} \right)}^3}{{\left( {4 + 5{{\bar c}_0} + 6{{\bar s}_0}} \right)}^3}}} \\& {b_2}=\frac{{(1 + {{\bar c}_0})({{\bar c}_0} - 2) - 3{{\bar s}_0}}}{{{{\left( {1 + {{\bar c}_0}} \right)}^3}\left( {4 + 5{{\bar c}_0} + 6{{\bar s}_0}} \right)}} \\ \end{aligned} \right. $

      梯度泡沫金属杆冲击端的初始相对密度ρ0可以由一维塑性冲击波模型给出:

      $ {\sigma _{\rm{B}}}(0)={s_0}\rho _0^{3/2} + {\rho _{\rm{s}}}{\rho _0}{v_0}\left( {{v_0} + {c_0}\rho _0^{1/4}} \right)={F_0}/{A_0} $

      在质量块冲击过程中,相对密度单调递增变化的梯度泡沫金属杆中塑性压溃波的传播满足一维冲击波假设,且质量块随波后压实区同步运动,即速度一致[13, 22],波后压实区应力记为σB。根据应力波理论和牛顿运动定律,如文献[22],可以获得:

      $ \left\{ \begin{gathered} \frac{{{\rm{d}} v(t)}}{{{\rm{d}} t}}= - \frac{{F(t)}}{M}= - \frac{{{\sigma _0}(\rho ) + {\rho _{\rm{s}}}\rho v(t)\left( {v(t) + c(\rho )} \right)}}{{{M \mathord{\left/ {\vphantom {M {{A_0}}}} \right. } {{A_0}}} + {\rho _{\rm{s}}}\int_0^{\varPhi (t)} {\rho (x)} {\rm{d}}x}} \hfill \\ \frac{{{\rm{d}} \varPhi }}{{{\rm{d}} t}}=v + c(\rho ) \hfill \\ \end{gathered} \right. $

      根据渐近求解的相对密度分布和对应的冲击情形代入上式,可以获得冲击力等参量随时间的变化。

    • 变胞元尺寸法可用于生成壁厚均一且密度分布连续变化的二维或三维的随机Voronoi细观有限元模型[6-7, 13]。与基于六角蜂窝构型的二维多胞模型的撒点规则不同,三维多胞结构的撒点规则是基于卡尔文十四面体提出的。Yang等 [13]基于均匀三维Voronoi结构引入了新的撒点规则,在一定区域内随机撒点,通过已知的相对密度分布和壁厚,控制任意两个相邻成核点ij之间的间距不小于当前位置的最小许可间距,记为

      $ {\delta _{{{ij}}}} \geqslant (1 - k) \times 2.9h/\rho ({x_{{{ij}}}},{y_{{{ij}}}},{z_{{{ij}}}}) $

      式中: h为胞元胞壁厚度,k为胞元不规则度,ρxijyijzij)为ij两核点连线中点位置对应的相对密度。本文中只考虑x方向上的相对密度分布,xij=(xi + xj)/2,其他两个方向相对密度均匀分布。

      梯度泡沫金属的基体材料模型和参数与文献[13]一致。基体材料假设为弹性–理想塑性,其密度ρs=2 700 kg/m3,杨氏模量Es=69 GPa,泊松比ν=0.3,屈服应力σys=165 MPa。采用壳单元S3R划分胞壁网格,通过网格收敛性分析[11],其壳单元的特征尺寸大小为约0.3 mm。梯度泡沫材料细观有限元模型的横截面为正方形,其面积A0=30×30 mm2。模型中所有可能的接触设置为摩擦因数为0.02的通用接触,冲击加载示例见图2

      图  2  质量块冲击梯度多胞杆

      Figure 2.  A diagram of 3D graded Voronoi models under mass impact

    • 上文已经给出了恒定冲击力耐撞性反向设计理论的一阶和二阶近似解,本文中采用与Yang等[13]冲击情形一致:冲击物体质量M=50 g,初始冲击速度v0=100 m/s,额定受力为F0=7.5 kN。由此可得,特征时间T0 =66.7 ms,特征长度L0 =66.7 mm。冲击端相对密度为0.086 9。参数${\bar c_0}=0.453$${\bar s_0}=2.10$,则有渐近解中的系数为a1=0.570和a2=0.191。

      通过已获得的上述参数表征相对密度分布,并与Yang等[13]的结果进行比较验证。根据初始压溃应力与相对密度的关系,获得相对密度梯度分布的渐近解,将Yang等[13]的数值结果作为精确解进行对比,如图3(a)所示。渐近解与精确解在冲击端的相对密度有小量差距。渐近解的结果0.086 9与精确解的结果0.085 0相比,相对误差约为2%,可以忽略。比较渐近解与精确解给出的相对密度分布表明,一阶近似解只是线性近似,冲击端附近差异较小,越靠近远端与精确解差异越大;二阶近似解的整体差异相比较小,冲击端附近稍大于精确解,远端附近稍小于精确解。

      图  3  相对密度分布的渐近解和冲击载荷历程的理论预测值

      Figure 3.  Theoretical predictions of asymptotic solutions of relative density distribution and the history curves of impact force

      结合应力波理论和牛顿第二运动定律,理论预测特定密度分布的梯度泡沫金属材料的力学行为。将相对密度分布和特定冲击情形代入公式进行理论预测,运用四阶龙格库塔法求解,可以获得冲击力、速度和波阵面位置随时间的演化过程。通过观察冲击力随时间的演化过程,二阶近似和精确解模型的结果与耐撞性额定冲击力值非常吻合,而一阶近似模型结果初期吻合较好,后期差别较大,如图3(b)所示。通过冲击速度降为0时确定波阵面传播停止的位置,获得梯度泡沫金属杆的长度分别为:69.4 mm(一阶近似模型)和65.9 mm(二阶近似模型)。精确解给出的梯度泡沫金属杆的长度为66.3 mm,如图4(b)所示。通过理论预测,二阶近似解具有足够的精度,可用于设计具有恒定冲击载荷的梯度泡沫金属材料细观有限元模型的相对密度分布。

      图  4  冲击速度与冲击波波阵面位置的历史曲线

      Figure 4.  Evolution history of impact velocity and location of shock front

    • 基于已获得的相对密度梯度分布,利用变胞元尺寸法,生成梯度泡沫金属杆的细观有限元模型,并生成均匀密度分布的泡沫金属细观有限元模型作对比。均匀密度模型与二阶近似密度梯度模型的整体尺寸一致,平均相对密度一致,相对密度为0.117,三维Voronoi模型如图5所示。取中轴面x=15mm处剖面图观察其变形过程,可以发现一阶近似和二阶近似模型在刚性板的冲击过程中,以逐层压溃变形模式为主,在密度较小端胞元开始发生压溃,并形成一道较为稳定的塑性冲击波向远端传播,如图6(a)(b)。而对于均匀密度分布的细观有限元模型,其撞击初始变形以冲击端胞元压溃为主,形成一道向远端传播的冲击波。随着刚性板的继续推进,由于端部效应,杆内应力波在固支端反射加强,其应力超过了支撑端部分胞元的初始压溃应力,随即发生压溃变形。由于冲击情形为质量块初速度撞击,冲击能量有限,端部反射波只能使得少量胞元发生压溃,并不能使得该压溃波传播很远,最终压溃带宽度约为3~4个胞元直径,如图6(c)。比较密度均匀和梯度分布多胞模型变形图,发现合理地设计相对密度分布,可以有效地控制压溃波的形成和传播。

      图  5  密度梯度多胞杆细观有限元模型

      Figure 5.  Cell-based finite element models of density gradient cellular rods

      图  6  梯度多胞细观有限元模型中轴剖面变形图

      Figure 6.  Deformation patterns of cell-based finite element models in the middle section of the density gradient cellular rods perpendicular to x-axis

      梯度泡沫金属杆冲击端发生局部化压溃,产生塑性冲击波向远端传播。波前未压实区应力和速度分别为为σ0和0,其中初始压溃应力为σ0(ρ(Φ(t)))=s0ρ3/2,即当前t时刻波阵面位置处对应的初始压溃应力,则支撑端力Fsup=A0σ0。有限元计算的支撑端和冲击端力的结果与相对密度分布精确解模型的理论预测结果比较见图7。有限元计算的结果出现震荡,其原因为:刚性板初始接触阶段,试件和刚性板间的罚函数接触算法会导致震荡;后期频率较高的震荡为胞壁内部的弹性波来回传播及交互作用导致的;后期低频且幅值较大的震荡是由于胞元的失稳压溃产生的。一阶近似下的有限元模型的计算结果与理论预测值比较:初始阶段有限元计算的冲击力结果偏高于额定冲击力,后逐渐趋于额定值并重合;有限元计算的支撑端力结果与理论预测结果初期呈递增趋势,初始阶段有限元计算的结果高于理论预测值,上升斜率较小,后逐渐低于理论预测值,整体相比差异明显。二阶近似下的有限元模型的计算结果与理论预测值比较表明,冲击力结果整体维持较为稳定的状态,在某一恒定值附近震荡,其值稍大于额定设计值;支撑端力结果与理论预测值初期都呈现递增趋势,初始阶段有限元计算结果稍高于理论预测值,后逐渐趋于理论预测值,整体偏高。二阶近似模型相比于一阶近似模型可以很好地控制冲击力在某一恒定值附近震荡,其值偏高于额定值,与Yang等[13]的结果基本一致,如图8(a)所示。这表明二阶近似解可以很好地表征精确解,而一阶近似解的精度不够。初始压溃应力与相对密度3/2次幂的关系可以应用于指导梯度泡沫材料的耐撞性设计。

      图  7  冲击端和支撑端载荷理论预测与有限元计算结果的比较

      Figure 7.  Comparisons of impact force and support force curves history between theoretical predictions and finite element (FE) results

      图  8  有限元计算结果对比

      Figure 8.  Comparisons of finite element (FE) results

      与二阶近似模型结果对比,均匀密度分布的泡沫金属模型的冲击力结果具有明显的峰值且高于额定载荷,后随时间衰减。由于冲击初期,质量块冲击速度最高,冲击端相对密度(0.117)大于二阶近似模型(0.086 9),其初始压溃应力和应变硬化参数较高,使得初期冲击力明显高于二阶近似模型,如图8(a)所示。随冲击速度的降低,其对应的冲击力也相应减小,出现明显的冲击力峰值后衰减,而二阶近似模型相对密度单调递增,补偿了因速度衰减对冲击力的影响,很好地控制了冲击力的稳定演化。观察物体冲击速度随时间的衰减过程发现,有限元模拟结果要比理论预测衰减偏快,用时稍短。以冲击速度降为0时表征冲击结束,一阶和二阶近似模型实际冲击时间分别为0.64 ms和0.62 ms,而均匀密度分布模型速度衰减得最快且不稳定,冲击时间为0.56 ms,耗时最短,如图8(b)所示。合理地优化梯度泡沫金属材料相对密度的分布,对于控制质量块冲击力是必要的。

      图7图8中有限元模拟的结果偏高于理论预测值,可能原因是本文中采用的冲击波模型未考虑多胞材料的率敏感性影响,初始压溃应力和应变硬化参数均使用准静态应力应变关系拟合得到,导致理论值低估有限元模拟的结果。已有研究表明,即使在未考虑基体材料应变率效应的影响下,多胞材料的动态应力应变关系与准静态的应力应变关系有明显的不同。Zheng等[11]发现了泡沫金属一种特有的率敏感性——“冲击速率敏感性”,即多胞材料动态压实区的应力应变状态对应于特定的冲击速度,不同冲击速度下的状态点可以用一条曲线描述,该曲线可用于表征多胞材料动态下的应力应变关系。通过比较动态应力应变关系和准静态的应力应变关系发现,动态初始压溃应力高于准静态初始压溃应力,而动态应变硬化参数小于准静态应变硬化参数。Ding等[23]和Wang等[24]进一步研究发现了泡沫金属的初始压溃应力具有应变率敏感性,在应变率高到一定程度后二者有幂率关系。因而,后期研究中可以考虑泡沫金属材料率敏感性的影响,以提高设计的准确度。在实验研究方面,He等[25]提出了通过控制水冷前发泡时间可以制备连续密度梯度变化的闭孔泡沫金属,将有助于后期耐撞性反向设计理论的实验验证。

    • 本文中基于一维率无关刚性–塑性硬化(R-PH)模型,获得了泡沫金属的初始压溃应力和应变硬化参数与相对密度呈3/2次幂的关系。针对于恒定冲击力的耐撞性要求,对梯度泡沫金属的耐撞性反向设计理论进行简化,运用级数法,获得了闭孔梯度泡沫金属材料相对密度分布的渐近解,与Yang等[13]的精确解进行了比较。通过理论结果的比较发现:两者初始相对密度具有微小差异,但可以忽略;一阶近似解是对精确解的线性近似,误差较大;二阶近似与精确解吻合较好。通过冲击力理论预测判断,二阶近似密度分布可以达到恒定冲击力要求,而一阶近似结果在后期无法确保冲击力维持恒值。根据相对密度分布的渐近解,设计三维梯度细观有限元模型并进行计算验证。通过与数值模拟结果的比较发现:二阶近似模型冲击力一直维持较为稳定的恒定冲击力,很好地满足了恒定冲击载荷的耐撞性要求;一阶近似模型冲击力无法维持较为稳定的载荷;均匀密度分布的多胞模型冲击力具有明显的峰值,后随时间衰减。相比均匀密度多胞材料,单调递增的二阶近似密度梯度设计,合理地分配密度的强弱,控制载荷稳定演化,使得速度平稳降低,提高了撞击过程物体的稳定性。综上所述,耐撞性反向设计理论简化模型的渐近解对于梯度闭孔泡沫材料的耐撞性设计是有效的,所提出的耐撞性设计方法在控制冲击吸能过程和冲击物受载方面具有指导意义。

参考文献 (25)

目录

    /

    返回文章
    返回