冲击载荷下仿贝壳珍珠层Voronoi结构的动态力学响应

武晓东 张海广 王瑜 孟祥生

引用本文:
Citation:

冲击载荷下仿贝壳珍珠层Voronoi结构的动态力学响应

    作者简介: 武晓东(1983-),男,博士,讲师,主要从事复合材料动力学研究. E-mail:wuxiaodong@tyut.edu.cn;
  • 中图分类号: O347.1

Dynamic Responses of Nare-Like Voronoi Structure under Impact Loading

  • CLC number: O347.1

  • 摘要: 贝壳珍珠层复合结构是一种有效的抗压结构系统,微观上具有Voronoi随机结构,具有良好的力学特性。为了研究仿贝壳珍珠层Voronoi结构在冲击载荷下的动态力学响应,建立了一种铝/乙烯基复合材料的三维Voronoi模型。首先,应用随机Voronoi技术建立仿贝壳珍珠层Voronoi随机模型,然后在随机多边形铝片之间引入黏结层,模拟黏结和分层过程,从最大变形、损伤分布和耗散能量等方面探讨Voronoi片板模型在弹丸冲击荷载作用下的力学性能,并与规则片板模型进行对比分析。结果显示:Voronoi结构更有利于冲击能量的扩散与吸收,减小应力集中,更好地发挥能量共享机制;而规则模型的冲击损伤主要集中在弹丸冲击点附近区域。最后讨论了黏结层厚度和分块尺寸对Voronoi模型力学性能的影响,结果表明:分块尺寸对Voronoi模型抗冲击性能的影响很小;黏结层对损伤耗散能和塑性能的影响很明显,黏结层越薄,模型的抗冲击性能越好。
  • 图 1  (a) Voronoi图的初始网格构型;(b)每个站点都在一个圆圈区域内随机移动;(c)新的Voronoi图是从新的站点系统中生成,通过矩形裁剪,形成有限区域的随机Voronoi图[25]

    Figure 1.  (a) Initial grid formation of Voronoi diagram; (b) each site moves randomly within a circle region; (c) a new Voronoi diagram is generated from the new site system, and by a rectangle cut, a finite Voronoi diagram is generated[25]

    图 2  规则8 × 8模型和随机Voronoi模型的示意图

    Figure 2.  Schematic of regular 8 × 8 model and random Voronoi models

    图 3  模型边界条件和加载条件

    Figure 3.  Boundary condition and loading of the model

    图 4  不同时刻规则片板模型的von Mises应力云图剖视图

    Figure 4.  Cutaway views of von Mises stress contours of regular plate model at different times

    图 5  不同时刻Voronoi片板模型的von Mises应力云图剖视图

    Figure 5.  Cutaway views of von Mises stress contours of Voronoi plate model at different times

    图 6  3.00 ms时规则片板模型的von Mises应力云图俯视图

    Figure 6.  Top view of von Mises stress contours of regular plate model at 3.00 ms

    图 7  3.00 ms时Voronoi片板模型的von Mises应力云图俯视图

    Figure 7.  Top view of von Mises stress contours of Voronoi plate model at 3.00 ms

    图 8  规则模型和不同分块尺寸Voronoi模型的损伤耗能

    Figure 8.  Damage energy of regular model and Voronoi models with different block sizes

    图 9  规则模型和不同分块尺寸Voronoi模型的塑性能

    Figure 9.  Plastic energy of regular model and Voronoi models with different block sizes

    图 10  具有不同黏性层厚度的Voronoi模型的损伤耗能

    Figure 10.  Damage energy of Voronoi models with different adhesive thicknesses

    图 11  具有不同黏性层厚度的Voronoi模型的塑性能

    Figure 11.  Plastic energy of Voronoi model with different adhesive thicknesses

    表 1  铝片的材料参数

    Table 1.  Parameters of aluminum plate

    Materialρ/(kg·m−3νE/GPa
    Aluminum27500.372
    下载: 导出CSV

    表 2  Cohesive模型的材料参数

    Table 2.  Parameters of Cohesive model

    $t_{\rm{n} }^{{0} }$/MPa$t_{\rm{s} }^{{0} }$/MPa$t_{\rm{t} }^{{0} }$/MPa$G_{\rm{n} }^{{0} }$/(kJ·m−2)$G_{\rm{s} }^{{0} }$/(kJ·m−2)$G_{\rm{t} }^{{0} }$/(kJ·m−2)ρ/(kg·m−3)Es/GPaEt/GPa
    8080801111 85041.5
    下载: 导出CSV

    表 3  Johnson-Cook模型参数[25]

    Table 3.  Parameters of Johnson-Cook model[25]

    MaterialA/MPaB/MPanCm
    Aluminum3916840.4360.009592
    下载: 导出CSV
  • [1] CURREY J D. Mechanical properties of mother of pearl in tension [J]. Proceedings of the Royal Society B: Biological Sciences, 1977, 196(1125): 443–463. doi: 10.1098/rspb.1977.0050
    [2] BARTHELAT F, TANG H, ZAVATTIERI P D, et al. On the mechanics of mother-of-pearl: a key feature in the material hierarchical structure [J]. Journal of the Mechanics and Physics of Solids, 2007, 55(2): 306–337. doi: 10.1016/j.jmps.2006.07.007
    [3] IMBALZANO G, TRAN P, NGO T D, et al. A numerical study of auxetic composite panels under blast loadings [J]. Composite Structures, 2016, 135: 339–352. doi: 10.1016/j.compstruct.2015.09.038
    [4] NGUYEN Q T, NGO T D, TRAN P, et al. Influences of clay and manufacturing on fire resistance of organoclay/thermoset nanocomposites [J]. Composites Part A: Applied Science and Manufacturing, 2015, 74: 26–37. doi: 10.1016/j.compositesa.2015.03.014
    [5] PRO J W, LIM R K, PETZOLD L R, et al. GPU-based simulations of fracture in idealized brick and mortar composites [J]. Journal of the Mechanics and Physics of Solids, 2015, 80: 68–85. doi: 10.1016/j.jmps.2015.03.011
    [6] ZHANG N, YANG S F, XIONG L M, et al. Nanoscale toughening mechanism of nacre tablet [J]. Journal of the Mechanical Behavior of Biomedical Materials, 2016, 53: 200–209. doi: 10.1016/j.jmbbm.2015.08.020
    [7] SARIKAYA M. An introduction to biomimetics: a structural viewpoint [J]. Microscopy Research and Technique, 1994, 27(5): 360–375. doi: 10.1002/jemt.1070270503
    [8] GAO H J, JI B H, JÄGER I L, et al. Materials become insensitive to flaws at nanoscale: lessons from nature [J]. Proceedings of the National Academy of Sciences of the United States of America, 2003, 100(10): 5597–5600. doi: 10.1073/pnas.0631609100
    [9] SONG F, SOH A K, BAI Y L. Structural and mechanical properties of the organic matrix layers of nacre [J]. Biomaterials, 2003, 24(20): 3623–3631. doi: 10.1016/S0142-9612(03)00215-1
    [10] BARTHELAT F. Nacre from mollusk shells: a model for high-performance structural materials [J]. Bioinspiration & Biomimetics, 2010, 5(3): 035001. doi: 10.1088/1748-3182/5/3/035001
    [11] TRAN P, NGO T D, MENDIS P. Bio-inspired composite structures subjected to underwater impulsive loading [J]. Computational Materials Science, 2014, 82: 134–139. doi: 10.1016/j.commatsci.2013.09.033
    [12] FLORES-JOHNSON E A, SHEN L M, GUIAMATSIA I, et al. A numerical study of bioinspired nacre-like composite plates under blast loading [J]. Composite Structures, 2015, 126: 329–336. doi: 10.1016/j.compstruct.2015.02.083
    [13] DUTTA A, TEKALUR S A. Crack tortuousity in the nacreous layer-topological dependence and biomimetic design guideline [J]. International Journal of Solids and Structures, 2014, 51(2): 325–335. doi: 10.1016/j.ijsolstr.2013.10.006
    [14] KOTHA S P, LI Y, GUZELSU N. Micromechanical model of nacre tested in tension [J]. Journal of Materials Science, 2001, 36(8): 2001–2007. doi: 10.1023/A:1017526830874
    [15] BARTHELAT F, RABIEI R. Toughness amplification in natural composites [J]. Journal of the Mechanics and Physics of Solids, 2011, 59(4): 829–840. doi: 10.1016/j.jmps.2011.01.001
    [16] SONG F, BAI Y L. Effects of nanostructures on the fracture strength of the interfaces in nacre [J]. Journal of Materials Research, 2003, 18(8): 1741–1744. doi: 10.1557/JMR.2003.0239
    [17] RABIEI R, BEKAH S, BARTHELAT F. Failure mode transition in nacre and bone-like materials [J]. Acta Biomaterialia, 2010, 6(10): 4081–4089. doi: 10.1016/j.actbio.2010.04.008
    [18] SHAO Y, ZHAO H P, FENG X Q, et al. Discontinuous crack-bridging model for fracture toughness analysis of nacre [J]. Journal of the Mechanics and Physics of Solids, 2012, 60(8): 1400–1419. doi: 10.1016/j.jmps.2012.04.011
    [19] GHAZLAN A, NGO T, VAN LE T, et al. Blast performance of a bio-mimetic panel based on the structure of nacre: a numerical study [J]. Composite Structures, 2020, 234: 111691. doi: 10.1016/j.compstruct.2019.111691
    [20] KO K, JIN S, LEE S E, et al. Impact resistance of nacre-like composites diversely patterned by 3D printing [J]. Composite Structures, 2020, 238: 111951. doi: 10.1016/j.compstruct.2020.111951
    [21] WU K J, ZHENG Z J, ZHANG S S, et al. Interfacial strength-controlled energy dissipation mechanism and optimization in impact-resistant nacreous structure [J]. Materials & Design, 2019, 163: 107532. doi: 10.1016/j.matdes.2018.12.004
    [22] 刘金义, 刘爽. Voronoi图应用综述 [J]. 工程图学学报, 2004, 25(2): 125–132. doi: 10.3969/j.issn.1003-0158.2004.02.023
    LIU J Y, LIU S. A survey on applications of Voronoi diagrams [J]. Journal of Engineering Graphics, 2004, 25(2): 125–132. doi: 10.3969/j.issn.1003-0158.2004.02.023
    [23] BARGMANN S, KLUSEMANN B, MARKMANN J, et al. Generation of 3D representative volume elements for heterogeneous materials: a review [J]. Progress in Materials Science, 2018, 96: 322–384. doi: 10.1016/j.pmatsci.2018.02.003
    [24] BAHMANI A, LI G, WILLETT T L, et al. Generating realistic representative microstructure of biomimetic composite materials for computational assessment of mechanical properties [J]. Materials Today: Proceedings, 2019, 7: 373–381. doi: 10.1016/j.matpr.2018.11.098
    [25] GHAZLAN A, NGO T D, TRAN P. Three-dimensional Voronoi model of a nacre-mimetic composite structure under impulsive loading [J]. Composite Structures, 2016, 153: 278–296. doi: 10.1016/j.compstruct.2016.06.020
  • [1] 王鹏李世强于国际吴桂英 . 轴向冲击载荷作用下梯度圆管的动力学行为. 高压物理学报, doi: 10.11858/gywlxb.2017.06.013
    [2] 黄涛陈鹏万张国新杨军 . 岩石Hopkinson层裂的流形元法模拟. 高压物理学报, doi: 10.11858/gywlxb.2006.04.003
    [3] 李汶蔚梅杰黄威 . 碳纤维增强复合材料层合板的抗冲击性能. 高压物理学报, doi: 10.11858/gywlxb.20190822
    [4] 刘静楠叶常青刘桂森沈耀 . 高温、高压、高应变速率动态过程晶体塑性有限元理论模型及其应用. 高压物理学报, doi: 10.11858/gywlxb.20190874
    [5] 王志华李志强任会兰赵隆茂 . 有机玻璃的几何形状对爆炸切割影响的研究. 高压物理学报, doi: 10.11858/gywlxb.2010.01.008
    [6] 胡昌明李雪梅彭建祥张祖根靳开诚王翔蔡灵仓刘仓理 . 冲击载荷下K9玻璃的光学特性. 高压物理学报, doi: 10.11858/gywlxb.2017.05.010
    [7] 梁瑞周文海余建平李珍宝杜超飞王敦繁 . 冲击载荷作用下岩体拉-压损伤破坏的边坡抛掷爆破模拟. 高压物理学报, doi: 10.11858/gywlxb.20180535
    [8] 傅华刘仓理王文强李涛 . 冲击动力学中离散元与有限元相结合的计算方法研究. 高压物理学报, doi: 10.11858/gywlxb.2006.04.007
    [9] 刘永贵沈玲燕 . 固定温度界面对相变波传播规律的影响. 高压物理学报, doi: 10.11858/gywlxb.20170559
    [10] 张志春强洪夫高巍然 . 光滑粒子流体动力学有限元法接触算法研究. 高压物理学报, doi: 10.11858/gywlxb.2011.02.001
    [11] 郑松林 . 晶体塑性有限元在材料动态响应研究中的应用进展. 高压物理学报, doi: 10.11858/gywlxb.20190725
    [12] 林华令 . 有限元法模拟混合物的冲击压缩特性. 高压物理学报, doi: 10.11858/gywlxb.1998.01.007
    [13] 彭克锋潘昊赵凯郑志军虞吉林 . 铜粉末动态压缩行为的多颗粒有限元分析. 高压物理学报, doi: 10.11858/gywlxb.20180665
    [14] 刘丽滨李海涛刁爱民王晓强 . 水下爆炸下有限尺度平板的载荷特性及结构响应试验研究. 高压物理学报, doi: 10.11858/gywlxb.20180516
    [15] 齐娟穆朝民 . 水射流对煤体冲击的有限元与光滑粒子耦合法数值模拟. 高压物理学报, doi: 10.11858/gywlxb.2014.03.016
    [16] 刘静楠叶常青陈开果俞宇颖沈耀 . <100> LiF高速冲击变形过程的晶体塑性有限元模拟. 高压物理学报, doi: 10.11858/gywlxb.20180551
    [17] 邓荣兵金先龙陈峻沈建奇陈向东 . 爆炸冲击波对玻璃幕墙破坏作用的多物质ALE有限元模拟. 高压物理学报, doi: 10.11858/gywlxb.2010.02.001
    [18] 王礼立朱珏赖华伟 . 冲击动力学研究中实测波信息的解读分析. 高压物理学报, doi: 10.11858/gywlxb.2010.04.007
    [19] 胡晓棉 . 可变元胞分子动力学方法对晶体结构稳定性的研究. 高压物理学报, doi: 10.11858/gywlxb.1989.02.005
    [20] 王永刚刘宏伟 . 强冲击载荷下含杂质的纯铝中微孔洞长大的动力学行为. 高压物理学报, doi: 10.11858/gywlxb.2010.04.002
  • 加载中
图(11)表(3)
计量
  • 文章访问数:  278
  • 阅读全文浏览量:  245
  • PDF下载量:  0
出版历程
  • 收稿日期:  2020-05-21
  • 录用日期:  2020-06-12
  • 网络出版日期:  2020-08-12

冲击载荷下仿贝壳珍珠层Voronoi结构的动态力学响应

    作者简介:武晓东(1983-),男,博士,讲师,主要从事复合材料动力学研究. E-mail:wuxiaodong@tyut.edu.cn
  • 1. 中国辐射防护研究院核应急与核安全研究所,山西 太原 030009
  • 2. 太原理工大学机械与运载工程学院,山西 太原 030024

摘要: 贝壳珍珠层复合结构是一种有效的抗压结构系统,微观上具有Voronoi随机结构,具有良好的力学特性。为了研究仿贝壳珍珠层Voronoi结构在冲击载荷下的动态力学响应,建立了一种铝/乙烯基复合材料的三维Voronoi模型。首先,应用随机Voronoi技术建立仿贝壳珍珠层Voronoi随机模型,然后在随机多边形铝片之间引入黏结层,模拟黏结和分层过程,从最大变形、损伤分布和耗散能量等方面探讨Voronoi片板模型在弹丸冲击荷载作用下的力学性能,并与规则片板模型进行对比分析。结果显示:Voronoi结构更有利于冲击能量的扩散与吸收,减小应力集中,更好地发挥能量共享机制;而规则模型的冲击损伤主要集中在弹丸冲击点附近区域。最后讨论了黏结层厚度和分块尺寸对Voronoi模型力学性能的影响,结果表明:分块尺寸对Voronoi模型抗冲击性能的影响很小;黏结层对损伤耗散能和塑性能的影响很明显,黏结层越薄,模型的抗冲击性能越好。

English Abstract

  • 仿贝壳珍珠层复合材料作为新型复合材料,有着优异的力学性能,不仅具有很高的强度,而且具有很好的韧性,近年来引起了学术界的不断关注[1-6]。很多发达国家非常重视贝壳结构材料和仿生材料研究,如美国等国家设置了专门的经费来研究贝壳生物材料的仿生设计和性能,用于装甲防弹衣和防爆装置。贝壳珍珠层复合材料的优异力学性能与其微观结构密切联系,为此研究人员对珍珠质的微观结构特征(体积分数、片剂长宽比、重叠长度等)进行了深入分析,试图将其与模型的力学性能联系起来[7-12]。Dutta等[13]研究了珍珠层中裂纹的萌生规律,评估了重叠长度对裂纹尖端驱动力的影响。Kotha等[14]的研究显示,低纵横比的文石片可以制造出具有高韧性的复合材料。Barthelat等[15]发现,珍珠层没有实现稳定状态的裂纹扩展,并将其归因于片层拔出增韧机制。其他学者也发现珍珠层内部和外部韧化机制阻止了裂纹的扩展[16-21]

    Barthelat等[15]通过观察发现,在每层贝壳珍珠层中,平板的排列与Voronoi图相似,从一个红色鲍鱼标本的光学图像中可以看到每个贝壳层压板都有矿物片的随机分布,并与其他珍珠层成键。基于这些光学图像,他们生成了一个由两层贝壳的平板结构组成的几何模型,用于有限元分析。自1907年Shamos和Hoey提出分治算法的最初定义和描述之后,Voronoi图便成为众多学科的中心主题之一。Voronoi图所具有的自然描述性和操纵能力,使其获得了广泛应用[22-24]。尽管Voronoi图对科学和工程中的各种应用具有重大的潜在影响,但是在很多领域包括仿生结构领域,Voronoi结构对材料力学性能的影响还未得到透彻的理解,为此本工作将探讨Voronoi结构的随机性对仿贝壳珍珠层结构力学性能的影响。

    为了研究仿贝壳珍珠层Voronoi随机模型结构的动态力学响应,首先建立一种铝/乙烯基复合材料结构的三维Voronoi模型,然后对模型在弹丸冲击载荷下的动态力学性能进行有限元模拟分析,最后讨论黏结层厚度和Voronoi模型分块尺寸对模型抗冲击力学性能的影响。

    • 利用文献[25]给出的随机Voronoi技术生成Voronoi随机模型。图1描述了仿贝壳珍珠层随机Voronoi结构的生成技术。图1(a)显示了由网格组成的Voronoi初始构型,每个网格内都包含1个站点;站点可以在圆内随机移动,如图1(b)所示,站点位置$\left( {x,y} \right)$由极坐标控制方程决定

      图  1  (a) Voronoi图的初始网格构型;(b)每个站点都在一个圆圈区域内随机移动;(c)新的Voronoi图是从新的站点系统中生成,通过矩形裁剪,形成有限区域的随机Voronoi图[25]

      Figure 1.  (a) Initial grid formation of Voronoi diagram; (b) each site moves randomly within a circle region; (c) a new Voronoi diagram is generated from the new site system, and by a rectangle cut, a finite Voronoi diagram is generated[25]

      $ x = {x_0} + r\cos \; \theta ,\;\;\;\;y = {y_0} + r\sin \;\theta $

      式中:$\left( {{x_0},{y_0}} \right)$为网格内站点的参考位置;$r$$\left( {0,R} \right)$之间的随机值,$R$为圆的半径;$\theta $$\left( {0,2{\text{π}} } \right)$之间的随机值。通过站点的随机移动,最终形成新的随机Voronoi图,如图1(c)所示。新生成的Voronoi图由一个无限大的区域组成,采用矩形切割使新生成的Voronoi图限定在有限区域内,丢弃有限区域外的站点和多边形网格,最终形成仿贝壳珍珠层随机Voronoi模型。

    • 图2给出了规则片板单元模型以及4种不同分块尺寸的不规则Voronoi片板模型,5种模型的总体几何尺寸相同,均为240 mm × 240 mm × 15 mm,模型总层数均为5层。图2(a)为规则片板模型,每个规则片板的几何尺寸为30 mm × 30 mm × 3 mm,每层由8 × 8共64个片板组成,图2(b)图2(e)分别给出了7 × 7、8 × 8、9 × 9和10 × 10分块的Voronoi不规则模型,每种模型均包括5层不同的随机单层结构,每层厚度为3 mm。

      图  2  规则8 × 8模型和随机Voronoi模型的示意图

      Figure 2.  Schematic of regular 8 × 8 model and random Voronoi models

      为了模拟仿贝壳珍珠层片层之间受冲击破坏时的脱黏现象,采用了内聚力Cohesive模型。通过合理的参数选择,内聚力Cohesive模型能够部分描述贝壳珍珠层内部层与层之间的变形和失效现象[2]。在片板之间以及板层之间插入Cohesive黏结层,考虑3种黏结层厚度0.1、0.2和0.3 mm,讨论黏结层厚度对模型冲击损伤的影响。

      内聚力Cohesive模型的牵引分离定律涉及黏性牵引应力矢量${{T}} = \left\{ {{t_{\rm{n}}},{t_{\rm{s}}},{t_{\rm{t}}}} \right\}$,其中下标n、s和t分别表示一个法向和两个切向分量。这些变量之间满足双线性二次黏聚律[25]

      $ {\left( {\dfrac{{{t_{\rm{n}}}}}{{t_{\rm{n}}^0}}} \right)^2} + {\left( {\dfrac{{{t_{\rm{s}}}}}{{t_{\rm{s}}^0}}} \right)^2} + {\left( {\dfrac{{{t_{\rm{t}}}}}{{t_{\rm{t}}^0}}} \right)^2} = 1 $

      式中:$t_{\rm{n}}^0$$t_{\rm{s}}^0$$t_{\rm{t}}^0$分别表示变形垂直于界面以及在第一、第二剪切方向上的最大应力。

      Cohesive黏结层刚度退化速率满足

      $ {\left( {\frac{{{G_{\rm{n}}}}}{{G_{\rm{n}}^0}}} \right)^2} + {\left( {\frac{{{G_{\rm{s}}}}}{{G_{\rm{s}}^0}}} \right)^2} + {\left( {\frac{{{G_{\rm{t}}}}}{{G_{\rm{t}}^0}}} \right)^2} = 1 $

      式中:$G_{\rm{n}}^{}$$G_{\rm{s}}^{}$$G_{\rm{t}}^{}$分别为正向和两个切向的断裂能,$G_{\rm{n}}^0$$G_{\rm{s}}^0$$G_{\rm{t}}^0$分别为正向和两个切向引起破坏所需的最大断裂能。

    • 仿贝壳珍珠层三维Voronoi结构模型包括两种材料模型,其中片层采用铝AA5083-H116,片层之间Cohesive黏性层使用乙烯树脂材料。表1表2列出了两种材料参数[19, 25],其中ρ为密度,ν为泊松比,E为弹性模量,EsEt为两个切向弹性模量。

      Materialρ/(kg·m−3νE/GPa
      Aluminum27500.372

      表 1  铝片的材料参数

      Table 1.  Parameters of aluminum plate

      $t_{\rm{n} }^{{0} }$/MPa$t_{\rm{s} }^{{0} }$/MPa$t_{\rm{t} }^{{0} }$/MPa$G_{\rm{n} }^{{0} }$/(kJ·m−2)$G_{\rm{s} }^{{0} }$/(kJ·m−2)$G_{\rm{t} }^{{0} }$/(kJ·m−2)ρ/(kg·m−3)Es/GPaEt/GPa
      8080801111 85041.5

      表 2  Cohesive模型的材料参数

      Table 2.  Parameters of Cohesive model

      在大应变情况下,铝合金的本构关系可采用Johnson-Cook模型描述

      $ \sigma = \left( {A + B{\varepsilon ^n}} \right)\left( {1 + C\ln \frac{{\dot \varepsilon }}{{{{\dot \varepsilon }_0}}}} \right)\left[ {1 - {{\left( {\frac{{T - {T_{\rm r}}}}{{T - {T_{\rm m}}}}} \right)}^{m}}} \right] $

      式中:$\sigma $为等效应力,$A$为材料在常温准静态下的屈服强度,$B$为应变强化系数,$n$为强化指数,$\varepsilon $为等效塑性应变,$\dot \varepsilon $为塑性应变率,${\dot \varepsilon _0}$为初始应变率,$C$为应变率敏感系数,$m$为温度软化系数,${T_{\rm{m}}}$为材料的熔点,$T{}_{\rm{r}}$为参考温度(取常温),$T$为当前温度。表3列出了Johnson-Cook模型参数[25]

      MaterialA/MPaB/MPanCm
      Aluminum3916840.4360.009592

      表 3  Johnson-Cook模型参数[25]

      Table 3.  Parameters of Johnson-Cook model[25]

    • 图3显示了模型的边界条件和加载条件。弹丸以18 m/s的初速度冲击Voronoi模型,弹丸速度属于中低速范围。弹丸模型上半部分是一个半径15 mm、长45 mm的圆柱体,下半部分是一个半径为15 mm的半球体,总长为60 mm,刚体属性。冲击载荷下仿贝壳珍珠层三维Voronoi模型的边界条件为4个侧边均完全固定,弹丸与复合结构模型的接触为通用接触,弹丸作用在复合结构模型中心。

      图  3  模型边界条件和加载条件

      Figure 3.  Boundary condition and loading of the model

      黏结层网格类型采用COH3D8,铝片网格类型采用C3D8R,黏结层网格大小为1 mm,铝片网格大小为2 mm。在该网格密度下,模型的网格单元总数达到167230。节点总数为567001时,最大应力值保持稳定,网格的收敛性较好。

    • 在珍珠层结构中,片板滑动机制被认为是激活内在和外在韧化机制的关键因素,可以阻止裂纹扩展。该机制分别引起内聚力和残余塑性应变,从而闭合裂纹。由于弹丸冲击载荷方向垂直于Voronoi板模型,冲击载荷破坏的主要形式是黏性层剥离,因此片板滑动引起的增韧机制在这种特定的冲击加载问题中不占主导地位。在冲击载荷下,损伤和变形耗散的能量比摩擦接触要多得多,Voronoi模型对珍珠层结构负载分配和能量吸收机制的影响是所要考虑的主要因素。本研究首先分析冲击载荷下模型的动态响应,在此基础上考察不规则Voronoi片板等几何因素对动态力学特性和能量分配的影响。

      图4图5分别为不同时刻规则片板模型与Voronoi片板模型受弹丸冲击时的应力云图剖视图。通过对比可以发现:在规则模型中,应力主要集中在弹丸冲击点及附近区域,远离冲击点区域的应力很小;在Voronoi模型中,应力分布区域更大,受力更加均匀。规则模型受冲击后很快就被冲破;而Voronoi模型的冲击模拟结果显示,其最大应力载荷小于规则模型,最终弹丸并未完全贯穿模型。

      图  4  不同时刻规则片板模型的von Mises应力云图剖视图

      Figure 4.  Cutaway views of von Mises stress contours of regular plate model at different times

      图  5  不同时刻Voronoi片板模型的von Mises应力云图剖视图

      Figure 5.  Cutaway views of von Mises stress contours of Voronoi plate model at different times

      图6图7分别为规则片板模型和Voronoi片板模型受弹丸冲击3.00 ms时应力云图的俯视图。从受弹丸冲击破坏情况来看:规则模型中脱黏现象主要集中在弹丸冲击点附近区域;而Voronoi模型的冲击影响区域更大,基本遍布整个模型。对于Voronoi模型,冲击区域发生变形时,其余片板受挤压后也发生了脱黏现象,吸收更多的冲击能量,从而有利于冲击能量的扩散与吸收,使模型的更多部分承担冲击负载,即增加承载区域,减小应力集中,更好地发挥能量共享机制。因此Voronoi片板模型抵抗冲击荷载的能力明显优于规则片板模型。

      图  6  3.00 ms时规则片板模型的von Mises应力云图俯视图

      Figure 6.  Top view of von Mises stress contours of regular plate model at 3.00 ms

      图  7  3.00 ms时Voronoi片板模型的von Mises应力云图俯视图

      Figure 7.  Top view of von Mises stress contours of Voronoi plate model at 3.00 ms

    • 图8图9分别给出了规则模型和不同分块Voronoi模型的损伤耗散能和塑性耗散能对比。在模型总尺寸相同的情况下,Voronoi模型的损伤耗散能远远高于规则片板模型,而塑性耗散能则小于规则片板模型,说明在冲击载荷作用下Voronoi模型抵抗冲击的能力优于规则片板模型。不同分块Voronoi模型的损伤耗散能和塑性耗散能差别不大,分块尺寸对Voronoi模型抗冲击性能的影响很小。

      图  8  规则模型和不同分块尺寸Voronoi模型的损伤耗能

      Figure 8.  Damage energy of regular model and Voronoi models with different block sizes

      图  9  规则模型和不同分块尺寸Voronoi模型的塑性能

      Figure 9.  Plastic energy of regular model and Voronoi models with different block sizes

      图10图11分别给出了不同黏结层厚度(h)的Voronoi模型的损伤耗散能和塑性耗散能。可以看出,黏结层对损伤耗散能和塑性耗散能的影响很明显。黏结层越薄,模型整体吸能越大,越薄的黏结层使模型具有更高的抗弯刚度,抗冲击性能越强。由此可见,Voronoi模型的不规则性是仿贝壳珍珠层复合结构模型抗冲击性能的影响因素,对贝壳结构韧性的提升发挥着重要作用。

      图  10  具有不同黏性层厚度的Voronoi模型的损伤耗能

      Figure 10.  Damage energy of Voronoi models with different adhesive thicknesses

      图  11  具有不同黏性层厚度的Voronoi模型的塑性能

      Figure 11.  Plastic energy of Voronoi model with different adhesive thicknesses

    • 通过有限元数值模拟研究了仿贝壳珍珠层Voronoi模型在弹丸冲击载荷下的动态力学响应,得到如下主要结论。

      (1)从冲击破坏受损情况来看,不规则Voronoi模型的冲击影响区域比规则模型更大,基本遍布整个模型。对于Voronoi模型,当冲击区域发生变形时,其余片板受挤压后发生脱黏现象,从而吸收更多的冲击能量,有利于冲击能量的扩散与吸收,让模型的更多部分承担冲击负载,即增加承载区域,并且减小应力集中,更好地发挥共享机制。规则模型的脱黏现象主要集中在弹丸冲击点及其附近区域。

      (2)在模型总尺寸相同的情况下,Voronoi片板模型的损伤耗散能远远高于规则片板模型,而塑性耗散能则小于规则模型。在冲击载荷作用下,不规则Voronoi片板模型抵抗冲击的能力优于规则片板模型。

      (3)分块尺寸对Voronoi模型抗冲击性能的影响很小,而黏结层对损伤耗散能和塑性耗散能的影响很明显,黏结层越薄,模型整体吸能越大,抗冲击性能越强。

      鉴于目前制备具有Voronoi结构的金属/高分子材料复合实验模型具有一定困难,因此未对模拟结果进行实验验证。随着3D打印技术的进一步发展,可采用金属/高分子材料混合3D打印技术制备仿贝壳珍珠层Voronoi实验模型,届时即可对铝/乙烯基复合三维Voronoi模型的数值模拟工作进行验证,进而开展更深入的研究。

参考文献 (25)

目录

    /

    返回文章
    返回