碳纤维-泡沫铝夹芯板低速冲击响应

刘姗姗 刘亚军 张英杰 李志强

引用本文:
Citation:

碳纤维-泡沫铝夹芯板低速冲击响应

    作者简介: 刘姗姗(1993-),女,硕士研究生,主要从事冲击动力学研究. E-mail:liushanshan0104@163.com;
    通讯作者: 李志强, lizhiqiang@tyut.edu.cn
  • 中图分类号: O344.1

Low-Velocity Impact Response of Carbon Fiber-Aluminum Foam Sandwich Plate

    Corresponding author: LI Zhiqiang, lizhiqiang@tyut.edu.cn
  • CLC number: O344.1

  • 摘要: 为研究夹芯结构的低速冲击响应,以碳纤维(T700)/环氧树脂复合材料层合板为上下面板,以闭孔泡沫铝为芯层,模拟夹芯板落锤冲击时的损伤演化过程。复合材料层合板采用三维实体单元建模,基于有限元软件ABAQUS中的用户子程序VUMAT,引入三维Hashin失效准则模拟复合材料的损伤破坏;采用二次应力准则,Cohesive单元模拟黏结层的层间失效;闭孔泡沫铝芯层采用3D Voronoi细观模型建模。分析复合材料夹芯结构在落锤冲击下的损伤起始、损伤扩展和最终破坏模式,通过锤头的接触力、位移、夹芯板的内能、后面板的最大位移研究夹层结构的能量吸收情况及抗冲击特性,得出了在质量保持不变的情况下,5种芯层相对密度和厚度的耦合关系中的最优设计是芯层相对密度15.0%,厚度为10 mm,为满足实际工程中的需求提供了设计依据。
  • 图 1  夹芯板结构示意图

    Figure 1.  Diagram of sandwich plate structure

    图 2  三维有限元模型

    Figure 2.  Three-dimensional finite element model

    图 3  能量-时间曲线

    Figure 3.  Energy-time curves

    图 4  能量-位移曲线

    Figure 4.  Energy-displacement curve

    图 5  不同冲击能量下的冲击载荷-时间曲线

    Figure 5.  Force-time curves under different impact energy

    图 6  33.0 J冲击能量下夹芯板的破坏模式

    Figure 6.  Failure mode of sandwich panel under 33.0 J impact energy

    图 7  不同冲击能量下锤头位移-时间曲线

    Figure 7.  Displacement-time curves of impactor at different impact energy

    图 8  不同冲击能量下锤头动能-时间曲线

    Figure 8.  Kinetic energy-time curves of impactor at different impact energy

    图 9  不同冲击能量下夹芯板内能-时间曲线

    Figure 9.  Internal energy-time curves of the sandwich board at different impact energy

    图 10  不同冲击能量下后面板的最大位移

    Figure 10.  Maximum displacement of rear panel at different impact energy

    图 11  不同夹芯结构的冲击载荷-时间曲线

    Figure 11.  The impact force-time curves for different sandwich structures

    图 12  不同夹芯结构的锤头位移-时间曲线

    Figure 12.  The displacement of impactor-time curves for different sandwich structures

    图 13  不同夹芯结构的锤头动能-时间曲线

    Figure 13.  Kinetic energy of impactor-time curves of impactor for different sandwich structures

    图 14  不同夹芯结构的芯层塑性耗散能-时间曲线

    Figure 14.  Plastic dissipation energy of core layer-time curves for different sandwich structures

    图 15  不同结构后面板最大位移与芯层厚度的比值曲线

    Figure 15.  The ratio curves of the maximum displacement of back layer to core thickness for different sandwich structures

    图 16  不同夹芯结构后面板各层的最大应力变化曲线

    Figure 16.  Maximum stress of each back layer for different sandwich structures

    图 17  冲击载荷-时间曲线

    Figure 17.  Force-time curves

    图 18  锤头位移-时间曲线

    Figure 18.  Displacement-time curves of impactor

    图 19  后面板应力云图

    Figure 19.  Stress plot of rear panel

    表 1  复合材料力学性能参数[11]

    Table 1.  Mechanical parameters of the composite materials[11]

    E1/GPaE2/GPaνG12/GPaG13/GPaG23/GPa
    180100.282.63.93.9
    Xt/MPaXc/MPaYt/MPaYc/MPaS12/MPaρ/(g·cm–3)
    2 5001 25060186851.95
    下载: 导出CSV

    表 2  Cohesive单元材料参数[14]

    Table 2.  Material parameters of cohesive elements[14]

    Knn/(GPa·mm–1)Kss(= Ktt)/(GPa·mm–1)N/MPaS(= T)/MPaG1/(J·m–2)G2(= G3)/(J·m–2)
    120433080520970
    下载: 导出CSV

    表 3  Al6061-T6材料参数

    Table 3.  Material parameters of Al6061-T6

    ρ/(g·cm–3)E/GPaνA/MPaB/MPaNm
    2.7700.282654260.341
    下载: 导出CSV

    表 4  不同能量下后面板的撕裂程度

    Table 4.  Tear degree of rear panel under different energy

    Impact energy/JTear layers
    33.00
    58.73
    91.75
    下载: 导出CSV

    表 5  5种不同的夹芯结构

    Table 5.  Five different sandwich structures

    Structure typePlane size of specimen/
    (mm × mm)
    Stacking sequenceUpper (lower) panel thickness /mmCore relative density/%The thickness of the core layer/mmDiameter of impactor/mmImpact energy/J
    1#100 × 100[45°/0°/−45°/90°]s110.015.012.533.0
    2#12.512.0
    3#15.010.0
    4#17.5 8.6
    5#20.0 7.5
    下载: 导出CSV

    表 6  5种不同结构后面板的撕裂程度

    Table 6.  Tear degree of rear panel for five sandwich structures

    Structure No.Tear layers
    1#0
    2#3
    3#5
    4#7
    5#8
    下载: 导出CSV
  • [1] 王巍, 安子军, 彭春彦, 等. 泡沫铝填充钢/铝复合管轴向抗冲击吸能特性 [J]. 哈尔滨工程大学学报, 2017, 38(7): 1093–1099.
    WANG W, AN Z J, PENG C Y, et al. Simulative research on the energy absorption characteristics of aluminum foam-filled steel/Al clad tube under axial impact loading [J]. Journal of Harbin Engineering University, 2017, 38(7): 1093–1099.
    [2] 骆伟, 谢伟, 刘敬喜. 芯层几何构形对复合材料波纹夹层结构冲击特性的影响 [J]. 江苏科技大学学报(自然科学版), 2018, 32(1): 21–26.
    LUO W, XIE W, LIU J X. Research on dynamic characteristics of a sandwich structures with various core shapes under impact loads [J]. Journal of Jiangsu University of Science and Technology (Natural Science Edition), 2018, 32(1): 21–26.
    [3] TITA V, CARVALHO J D, VANDEPITTE D. Failure analysis of low velocity impact on thin composite laminates: experimental and numerical approaches [J]. Composite Structures, 2008, 83(4): 413–428. doi: 10.1016/j.compstruct.2007.06.003
    [4] 韩守红, 吕振华. 铝泡沫夹层结构抗爆炸性能仿真分析及优化 [J]. 兵工学报, 2010, 31(11): 1468–1474.
    HAN S H, LÜ Z H. Numerical simulation of blast-resistant performance of aluminum foam sandwich structures and optimization [J]. Acta Armamentarii, 2010, 31(11): 1468–1474.
    [5] 李志斌, 卢芳云. 泡沫铝夹芯板压入和侵彻性能的实验研究 [J]. 振动与冲击, 2015(4): 1–5.
    LI Z B, LU F Y. Tests for indentation and perforation of sandwich panels with aluminium foam core [J]. Journal of Vibration and Shock, 2015(4): 1–5.
    [6] 赵金华, 曹海琳, 晏义伍, 等. 泡沫铝夹层结构复合材料低速冲击性能 [J]. 材料工程, 2018, 46(1): 92–98. doi: 10.11868/j.issn.1001-4381.2015.001295
    ZHAO J H, CAO H L, YAN Y W, et al. Low velocity impact properties of aluminum foam sandwich structural composite [J]. Journal of Materials Engineering, 2018, 46(1): 92–98. doi: 10.11868/j.issn.1001-4381.2015.001295
    [7] 荣誉. 梯度泡沫金属力学性能的Lagrangian分析 [D]. 太原: 太原理工大学, 2018.
    [8] HASHIN Z. Failure criteria for unidirectional fiber composites [J]. Journal of Applied Mechanics, 1980, 47: 329–334. doi: 10.1115/1.3153664
    [9] 谭开忍, 肖熙. 含有腐蚀缺陷海底管道极限载荷分析 [J]. 海洋工程, 2006, 24(3): 63–67. doi: 10.3969/j.issn.1005-9865.2006.03.010
    TAN K R, XIAO X. Analysis on limit load of corroded submarine pipelines [J]. The Ocean Engineering, 2006, 24(3): 63–67. doi: 10.3969/j.issn.1005-9865.2006.03.010
    [10] 沈鋆. 极限载荷分析法在压力容器分析设计中的应用 [J]. 石油化工设备, 2011, 40(4): 35–38. doi: 10.3969/j.issn.1000-7466.2011.04.010
    SHEN J. Limit load analysis application in pressure vessel analytical design [J]. Petro-Chemical Equipment, 2011, 40(4): 35–38. doi: 10.3969/j.issn.1000-7466.2011.04.010
    [11] 肖先林, 王长金, 赵桂平. 碳纤维复合材料-泡沫铝夹芯板的冲击响应 [J]. 振动与冲击, 2018, 37(15): 110–117.
    XIAO X L, WANG C J, ZHAO G P. Dynamic responses of carbon fiber composite sandwich panels with aluminum foam core subjected to impact loading [J]. Journal of Vibration and Shock, 2018, 37(15): 110–117.
    [12] JOHNSON G R, COOK W H. Fracture characteristics of three metals subjected to various strains, strains rates, temperatures and pressures [J]. Engineering Fracture Mechanics, 1985, 21(1): 31–48. doi: 10.1016/0013-7944(85)90052-9
    [13] 熊明洋, 向忠, 胡旭东, 等. 基于ABAQUS的CCF300碳纤维层合板低速冲击破坏数值模拟 [J]. 轻工机械, 2017, 35(4): 27–32. doi: 10.3969/j.issn.1005-2895.2017.04.006
    XIONG M Y, XIANG Z, HU X D, et al. Numerical simulation of low velocity impact failure of CCF300 carbon fiber laminate based on ABAQUS [J]. Light Industry Machinery, 2017, 35(4): 27–32. doi: 10.3969/j.issn.1005-2895.2017.04.006
    [14] 陈县辉. 基于内聚力单元的层合板低速冲击响应模拟研究 [D]. 太原: 中北大学, 2014.
    [15] STUBSS C. Compilation strategies: alternate approaches to achieve low power consumption [J]. Electronic Component News, 2008, 52(4): 11–113.
    [16] FOO C C, SEAH L K, CHAI G B. Low-velocity impact failure of aluminium honeycomb sandwich panels [J]. Composite Structures, 2008, 85(1): 20–28. doi: 10.1016/j.compstruct.2007.10.016
    [17] SAHU S, MONDAL D P, CHO J U, et al. Low-velocity impact characteristics of closed cell AA2014-SiCp composite foam [J]. Composites Part B: Engineering, 2019, 160: 394–401. doi: 10.1016/j.compositesb.2018.12.054
  • [1] 丁铁张晓晴姚小虎 . 筋条形状对复合材料加筋壁板低速冲击动态响应的影响. 高压物理学报, 2017, 31(6): 769-777. doi: 10.11858/gywlxb.2017.06.012
    [2] 廖斌斌周建武林渊贾利勇王栋亮花争立郑津洋顾超华 . CFRP层合板低速冲击响应及损伤特性研究. 高压物理学报, 2019, 33(4): 044202-1-044202-9. doi: 10.11858/gywlxb.20180699
    [3] 吴和成肖毅华 . 陶瓷和仿珍珠母陶瓷/聚脲复合结构的冲击损伤对比. 高压物理学报, 2020, 34(2): 024201-1-024201-7. doi: 10.11858/gywlxb.20190808
    [4] 王光勇裴晨浩林加剑 . 在顶爆作用下含裂隙锚固洞室的损伤演化及动态响应规律. 高压物理学报, 2018, 32(6): 064103-1-064103-9. doi: 10.11858/gywlxb.20180561
    [5] 郭扬李永池曹结东李建光 . 铝锂合金材料损伤演化方程和层裂准则的研究. 高压物理学报, 2005, 19(4): 312-318 . doi: 10.11858/gywlxb.2005.04.005
    [6] 朱锡李海涛牟金磊黄晓明 . 水下近距爆炸作用下船体梁的动态响应特性. 高压物理学报, 2010, 24(5): 343-350 . doi: 10.11858/gywlxb.2010.05.005
    [7] 郑松林 . 晶体塑性有限元在材料动态响应研究中的应用进展. 高压物理学报, 2019, 33(3): 030108-1-030108-21. doi: 10.11858/gywlxb.20190725
    [8] 贾宪振胡毅亭董明荣许学忠刘家骢 . 深水爆炸冲击波作用下圆柱壳动态响应影响因素的数值模拟研究. 高压物理学报, 2008, 22(2): 208-214 . doi: 10.11858/gywlxb.2008.02.016
    [9] 王永刚贺红亮王礼立经福谦 . 延性材料动态拉伸断裂早期连通过程的逾渗描述. 高压物理学报, 2006, 20(2): 127-132 . doi: 10.11858/gywlxb.2006.02.003
    [10] 李汶蔚梅杰黄威 . 碳纤维增强复合材料层合板的抗冲击性能. 高压物理学报, 2020, 34(2): 024101-1-024101-8. doi: 10.11858/gywlxb.20190822
    [11] 李平李大红宁建国经福谦 . 冲击载荷下Al2O3陶瓷的动态响应. 高压物理学报, 2002, 16(1): 22-28 . doi: 10.11858/gywlxb.2002.01.004
    [12] 李焰张向荣谭红梅刘晓玲裴明敬 . 国产PVDF压电薄膜的冲击加载及卸载响应研究. 高压物理学报, 2004, 18(3): 261-266 . doi: 10.11858/gywlxb.2004.03.011
    [13] 段卓平关智勇 . 冲击波中横向、纵向应力计响应初步研究. 高压物理学报, 2002, 16(4): 265-270 . doi: 10.11858/gywlxb.2002.04.005
    [14] 展婷变宁建国王志华 . 冲击载荷下钢筋混凝土力学行为的研究. 高压物理学报, 2016, 30(2): 109-115. doi: 10.11858/gywlxb.2016.02.004
    [15] 罗景润李大红张寿齐赵方芳赵玉华 . 简单拉伸下高聚物粘结炸药的损伤测量及损伤演化研究. 高压物理学报, 2000, 14(3): 203-208 . doi: 10.11858/gywlxb.2000.03.008
    [16] 贺红亮 . 动态拉伸断裂的物理判据研究. 高压物理学报, 2013, 27(2): 153-161. doi: 10.11858/gywlxb.2013.02.001
    [17] 李汶蔚黄威 . 金属蜂窝夹芯结构抗水下冲击性能. 高压物理学报, 2020, 34(3): 035102-1-035102-9. doi: 10.11858/gywlxb.20190790
    [18] 宿华祥易伟建 . 钢筋混凝土墙抗冲击性能的数值模拟分析. 高压物理学报, 2020, 34(1): 014201-1-014201-14. doi: 10.11858/gywlxb.20190772
    [19] 徐明利张若棋张光莹 . 用元胞法研究层合材料中的应力波. 高压物理学报, 2000, 14(3): 182-188 . doi: 10.11858/gywlxb.2000.03.005
    [20] 王光勇曹安生余锐林加剑 . 顶爆和拱腰侧爆同时作用下锚固洞室的动态响应. 高压物理学报, 2020, 34(2): 025202-1-025202-8. doi: 10.11858/gywlxb.20190812
  • 加载中
图(19)表(6)
计量
  • 文章访问数:  1939
  • 阅读全文浏览量:  979
  • PDF下载量:  17
出版历程
  • 收稿日期:  2019-12-23
  • 录用日期:  2020-02-04
  • 网络出版日期:  2020-05-15
  • 刊出日期:  2020-06-01

碳纤维-泡沫铝夹芯板低速冲击响应

    作者简介:刘姗姗(1993-),女,硕士研究生,主要从事冲击动力学研究. E-mail:liushanshan0104@163.com
    通讯作者: 李志强, lizhiqiang@tyut.edu.cn
  • 1. 太原理工大学机械与运载工程学院应用力学研究所,山西 太原 030024
  • 2. 太原理工大学材料强度与结构冲击山西省重点实验室,山西 太原 030024
  • 3. 力学国家级实验教学示范中心(太原理工大学),山西 太原 030024

摘要: 为研究夹芯结构的低速冲击响应,以碳纤维(T700)/环氧树脂复合材料层合板为上下面板,以闭孔泡沫铝为芯层,模拟夹芯板落锤冲击时的损伤演化过程。复合材料层合板采用三维实体单元建模,基于有限元软件ABAQUS中的用户子程序VUMAT,引入三维Hashin失效准则模拟复合材料的损伤破坏;采用二次应力准则,Cohesive单元模拟黏结层的层间失效;闭孔泡沫铝芯层采用3D Voronoi细观模型建模。分析复合材料夹芯结构在落锤冲击下的损伤起始、损伤扩展和最终破坏模式,通过锤头的接触力、位移、夹芯板的内能、后面板的最大位移研究夹层结构的能量吸收情况及抗冲击特性,得出了在质量保持不变的情况下,5种芯层相对密度和厚度的耦合关系中的最优设计是芯层相对密度15.0%,厚度为10 mm,为满足实际工程中的需求提供了设计依据。

English Abstract

  • 纤维增强复合材料以其高比强度、高比刚度、耐疲劳和耐腐蚀等优异性能被广泛应用于航空航天、航海、汽车等重要领域。泡沫铝以其密度小、抗冲击能力强、耐高温、抗腐蚀、易加工、成型精度高等优异的物理和化学性能,在工厂机房、机器设备、建筑、汽车防护等领域被广泛应用。而以纤维增强复合材料为面板、泡沫金属为芯层的夹芯结构,具有轻质、耐高温、隔热阻燃以及可以较好地吸收冲击能量等优良的特性,已大量应用于汽车防冲档、机械装置保护壳、飞机外壳夹层和太空飞行器的防护层等领域[1-2]

    近年来,研究人员对纤维增强复合材料夹芯结构的力学性能进行了大量研究。Tita等[3]通过进行低速冲击试验、压痕试验与准静态有限元模拟,将压痕试验结果与动力试验结果进行对比,验证了复合材料薄板在低速冲击时的惯性效应。韩守红等[4]对几种铝泡沫夹层结构的抗爆炸性能进行了有限元分析,从变形模式、运动响应和吸能特性等方面对比研究了夹层结构的动态响应特性,并找到了其中抗爆炸性能最优的夹层复合结构,对抗爆炸结构的工程设计具有一定指导意义。李志斌等[5]利用材料试验机(MTS)和落锤试验机研究了夹芯板结构在压入和侵彻时的变形和失效行为,得到不同压头/锤头形状和边界条件对泡沫铝夹芯板的压入和侵彻响应以及能量吸收性能的影响。赵金华等[6]使用泡沫铝芯层、玄武岩纤维(BF)和超高分子量聚乙烯纤维(UHMWPE)面板的复合夹芯结构,通过Instron 9350落锤式冲击试验机对夹层结构复合材料进行了冲击试验,并与蜂窝铝夹层结构的性能进行了对比,得到纤维类型、铺层结构设计、芯材对夹层结构复合材料冲击性能和损伤模式的影响规律。目前对复合材料渐进损伤分析大多采用试验和有限元方法,在研究吸能特性时,多采用比能量来衡量能量的吸收,忽略了芯层相对密度和高度的耦合作用对夹芯板结构在低速冲击下的力学性能影响。

    本研究考虑了芯层相对密度和高度耦合作用的影响,在保证夹芯板质量一定的情况下,通过分析复合材料夹芯结构在低速冲击下锤头的接触力、位移、动能,以及芯层的塑性耗散能、后面板的最大位移来研究夹层结构的能量吸收能力及抗冲击特性,从而得到5种不同芯层的相对密度与厚度的耦合关系中最优的设计,为满足实际工程需求提供设计依据。

    • 冲击试件由碳纤维-闭孔泡沫铝夹芯板组成,面板为T700碳纤维板,基体材料为环氧树脂,上下面板均为8层,单层厚度为0.125 mm,具体结构如图1所示。

      图  1  夹芯板结构示意图

      Figure 1.  Diagram of sandwich plate structure

    • 图2给出了有限元模拟的三维模型。碳纤维面板的网格划分采用渐进过渡的疏密网格,网格密度从板中间到两边依次降低,网格单元类型为C3D8R,层间界面采用Cohesive单元,单层厚度为0.01 mm,网格单元类型为COH3D8。芯层闭孔泡沫铝使用3D Voronoi随机模型建模,首先根据芯层的体积、相对密度等进行相关计算得出形核点的个数,给定空间内任意一个胞元核点坐标,采用Matlab软件中的随机函数并经过多次循环得到所有的核点坐标,通过MPT_Voronoi工具箱得到随机模型,随机模型中所有胞元的节点和面的坐标信息以ANSYS命令流的格式生成txt文本,将其导入ANSYS中,得到闭孔泡沫铝的几何模型,最后再以iges格式输出,在Hypermesh中划分网格,最终将其导入ABAQUS软件与碳纤维面板及锤头装配后进行相应计算[7]。闭孔泡沫铝芯层基体材料为Al6061-T6,有限元网格单元类型为S4R,面板与芯层之间采用通用接触。

      图  2  三维有限元模型

      Figure 2.  Three-dimensional finite element model

      夹芯板四周固支,冲击部位位于板中心,整个模型采用通用接触。面板为100 mm × 100 mm的正方形,芯层厚度以及冲击能量视结构类型而定,夹芯板的上下面板均为[45°/0°/−45°/90°]s铺层的碳纤维板。

    • 采用三维Hashin损伤失效准则[8-10]模拟落锤冲击下T700碳纤维面板层内的失效模式,具体的失效准则如下。

      纤维拉伸失效(σ11 ≥ 0)

      ${F_{{\rm {ft}}}} = {\left( {\frac{{{\sigma _{11}}}}{{{X_{\rm t}}}}} \right)^2} + {\left( {\frac{{{\tau _{12}}}}{{{S_{12}}}}} \right)^2} + {\left( {\frac{{{\tau _{13}}}}{{{S_{13}}}}} \right)^2} \geqslant 1$

      纤维压缩失效(σ11 < 0)

      ${F_{{\rm {fc}}}} = {\left( {\frac{{{\sigma _{11}}}}{{{X_{\rm c}}}}} \right)^2} \geqslant 1$

      基体拉伸失效(σ22 + σ33 ≥ 0)

      ${F_{{\rm{mt}}}} =\left(\frac{{{\sigma _{22}} + {\sigma _{33}}}}{{{Y_{\rm t}}}}\right)^2+ \frac{{{\sigma _{23}^2 - {\sigma _{22}}{\sigma _{33}}} }}{{S_{23}^2}} + \left( \frac{{{\tau _{12}}}}{{{S_{12}}}} \right)^2 + \left(\frac{{{\tau _{13}}}}{{{S_{13}}}}\right)^2 \geqslant 1$

      基体压缩失效(σ22 + σ33 < 0)

      ${F_{\rm {mc}}} = \frac{1}{{{Y_{\rm c}}}}\left[ {{{\left( {\frac{{{Y_{\rm c}}}}{{2{S_{12}}}}} \right)}^2} - 1} \right]\left( {{\sigma _{22}} + {\sigma _{33}}} \right) + {\left( {\frac{{{\sigma _{22}} + {\sigma _{33}}}}{{2{S_{12}}}}} \right)^2} + \frac{{ {\tau _{23}^2 - {\sigma _{22}}{\sigma _{33}}} }}{{S_{23}^2}} + {\left( {\frac{{{\tau _{12}}}}{{{S_{12}}}}} \right)^2} + {\left( {\frac{{{\tau _{13}}}}{{{S_{13}}}}} \right)^2} \geqslant 1$

      式中:Fft表示纤维拉伸断裂,Ffc表示纤维挤压断裂,Fmt表示基体拉伸断裂,Fmc表示基体挤压断裂,XtXc分别表示1方向的拉伸强度和压缩强度,YtYc表示2方向的拉伸强度和压缩强度,S12表示12层间剪切强度,S23表示23层间剪切强度,σ11σ22σ33为3个方向的正应力,τ12τ23τ13为3个方向的剪应力。

      以上4种损伤模式的相关材料参数如表1所示,其中E1E2分别表示1、2方向的弹性模量,G12G13G23为剪切模量分量,ν为泊松比,ρ为密度。

      E1/GPaE2/GPaνG12/GPaG13/GPaG23/GPa
      180100.282.63.93.9
      Xt/MPaXc/MPaYt/MPaYc/MPaS12/MPaρ/(g·cm–3)
      2 5001 25060186851.95

      表 1  复合材料力学性能参数[11]

      Table 1.  Mechanical parameters of the composite materials[11]

      层间的损伤采用内聚力(Cohesive)单元[12]模拟。内聚力单元可以有效模拟分层损伤的起始、扩展以及最终分层的发生[13]。具体材料参数见表2[14],表中KnnKssKtt为内聚力模型刚度,NST分别表示3个方向的拉伸强度,G1G2G3分别表示3个方向的断裂能。

      Knn/(GPa·mm–1)Kss(= Ktt)/(GPa·mm–1)N/MPaS(= T)/MPaG1/(J·m–2)G2(= G3)/(J·m–2)
      120433080520970

      表 2  Cohesive单元材料参数[14]

      Table 2.  Material parameters of cohesive elements[14]

      闭孔泡沫铝芯层基体材料为Al6061-T6,其应力-应变关系采用Johnson-Cook本构方程[15]描述。Johnson-Cook本构方程表达式为

      $\sigma = \left( {A + B{\varepsilon ^n}} \right)\left( {1 + C{\rm ln}\;{{\dot \varepsilon }^*}} \right)\left( {1 - {{ {{T^*}} }^m}} \right)$

      基于Johnson-Cook本构模型的失效表达式为

      $\varepsilon _{\rm u}^{\rm f} = \left[ {{D_1} + {D_2}\exp \left( {{D_3}\eta } \right)} \right]\left( {1 + {D_4}\ln {{{\dot \varepsilon }^*}} } \right)$

      式中:${\varepsilon } $为应变,${{\dot \varepsilon }^*} $为等效应变率,$\varepsilon _{\rm u}^{\rm f} $为失效应变,A为静态屈服强度,B为应变硬化模量,C为应变率敏感系数,T*为无量纲温度,n为应变硬化因子,m为温度软化因子,D1D4为材料参数。材料参数见表3

      ρ/(g·cm–3)E/GPaνA/MPaB/MPaNm
      2.7700.282654260.341

      表 3  Al6061-T6材料参数

      Table 3.  Material parameters of Al6061-T6

    • 以冲击能量为91.7 J、上下碳纤维面板铺层为[45°/0°/-45°/90°]s、闭孔泡沫铝芯层相对密度为10%、边界条件为四周固支、冲击部位为板中心为例来验证模型的有效性。在进行网格划分时,经过反复计算尝试,碳纤维面板的网格划分采用渐进过渡的疏密网格来提高有限元的计算效率。

      图3给出了夹芯板的能量-时间曲线、系统动能-时间曲线和沙漏能-时间曲线。夹芯板的内能和系统动能峰值分别为92.3 J和91.2 J,都与冲击总能量相差不超过1%,说明整个系统的能量是守恒的,夹芯板的吸能曲线以及能量-位移曲线(见图4)都与文献[16-17]中的曲线趋势一致,验证了模型及模拟过程的正确性。伪应变能的峰值为4.48 J,占冲击总能量的比值小于5%,说明沙漏能是可以控制的。

      图  3  能量-时间曲线

      Figure 3.  Energy-time curves

      图  4  能量-位移曲线

      Figure 4.  Energy-displacement curve

    • 为研究不同冲击能量下碳纤维-泡沫铝夹芯板的变形模式,选取3种冲击能量33.0、58.7、91.7 J,分别对相同的夹芯板进行冲击。通过控制冲击速度为6、8、10 m/s实现3种不同能量的冲击,泡沫铝芯层的相对密度为10.0%。

      图5给出了不同冲击能量下夹芯结构的冲击载荷-时间曲线,图6给出了碳纤维面板和芯层的破坏模式。结合图5图6可以看出,在夹芯板的冲击过程中,锤头作用区域的面板和芯层均发生局部破坏。在冲击起始阶段,碳纤维面板和芯层共同承受冲击力,冲击力迅速上升;随着冲击的继续,上层面板发生破坏直至被撕裂(见图6(a)图6(c)),芯层的泡沫铝被压缩并且出现明显的塑性变形,塑性变形的区域如图6(b)所示;冲击载荷达到峰值并出现一个较长的“平台区”,然后开始下降。冲击能量的大小决定了“平台区”冲击力的大小,冲击能量越大,冲击载荷的峰值就越高。表4给出了冲击过程中后面板的撕裂层数,冲击能量越大,后面板撕裂的层数越多。

      Impact energy/JTear layers
      33.00
      58.73
      91.75

      表 4  不同能量下后面板的撕裂程度

      Table 4.  Tear degree of rear panel under different energy

      图  5  不同冲击能量下的冲击载荷-时间曲线

      Figure 5.  Force-time curves under different impact energy

      图  6  33.0 J冲击能量下夹芯板的破坏模式

      Figure 6.  Failure mode of sandwich panel under 33.0 J impact energy

      图7给出了锤头位移随冲击时间的变化趋势。冲击起始阶段锤头位移增大,在冲击时间2.2 ms附近时,冲击位移达到最大值,之后减小。图8给出了动能随时间变化的曲线。冲击动能在2.2 ms之前一直减小,2.2 ms之后又变大,说明2.2 ms之后冲击过程已经完成,锤头发生回弹,同时冲击力减小,直到冲击力下降为零,冲击过程结束。可见,冲击持续时间与冲击能量并无关系,在不同的冲击能量下,冲击过程完成的时间没有明显差异。

      图  7  不同冲击能量下锤头位移-时间曲线

      Figure 7.  Displacement-time curves of impactor at different impact energy

      图  8  不同冲击能量下锤头动能-时间曲线

      Figure 8.  Kinetic energy-time curves of impactor at different impact energy

      图9给出了夹层结构内能随时间的变化曲线。在整个冲击过程中,夹层结构的内能先增大后减小,最后趋于稳定。内能大小用来表示夹层结构吸收的能量。图10给出了冲击能量与夹芯板后面板的位移的关系。在一定范围内,随着冲击能量的增大,后面板的最终位移也增大。

      图  9  不同冲击能量下夹芯板内能-时间曲线

      Figure 9.  Internal energy-time curves of the sandwich board at different impact energy

      图  10  不同冲击能量下后面板的最大位移

      Figure 10.  Maximum displacement of rear panel at different impact energy

    • 芯层的相对密度和芯层高度对夹芯板抗冲击性能都有影响。为研究芯层不同相对密度与芯层高度的耦合作用对碳纤维-泡沫铝夹芯板抗冲击性能的影响,保持芯层质量不变,设计5种夹芯结构,按照芯层相对密度从小到大(或芯层高度从大到小)对其进行编号,分别记为1#、2#、3#、4#和5#,初始冲击能量为33.0 J,参数如表5所示。

      Structure typePlane size of specimen/
      (mm × mm)
      Stacking sequenceUpper (lower) panel thickness /mmCore relative density/%The thickness of the core layer/mmDiameter of impactor/mmImpact energy/J
      1#100 × 100[45°/0°/−45°/90°]s110.015.012.533.0
      2#12.512.0
      3#15.010.0
      4#17.5 8.6
      5#20.0 7.5

      表 5  5种不同的夹芯结构

      Table 5.  Five different sandwich structures

      图11图12分别显示了5种夹芯结构受冲击时,锤头的冲击载荷-时间曲线和锤头位移-时间曲线。可以看出,在初始冲击能量相同的条件下,芯层相对密度为10.0%(厚度为15 mm)时,冲击力峰值明显小于其他4种相对密度下的峰值,“平台区”时间更长,锤头最大位移也更大。

      图  11  不同夹芯结构的冲击载荷-时间曲线

      Figure 11.  The impact force-time curves for different sandwich structures

      图  12  不同夹芯结构的锤头位移-时间曲线

      Figure 12.  The displacement of impactor-time curves for different sandwich structures

      图13给出了5种夹芯板受冲击时锤头的动能-时间曲线,锤头动能的变化趋势基本相同。因此芯层质量一定时,芯层相对密度和厚度在一定范围内,冲击过程中冲击动能无显著变化。

      图  13  不同夹芯结构的锤头动能-时间曲线

      Figure 13.  Kinetic energy of impactor-time curves of impactor for different sandwich structures

      图14给出了5种不同夹芯板的芯层的塑性耗散能-时间曲线。可以看出,3#结构的塑性耗散能较其他4种结构大。图15给出了不同结构后面板最大位移与芯层厚度的比值(δ/h)。可以看出:不同结构的后面板最大位移与芯层厚度的比值不同;随着芯层相对密度的增大和芯层厚度的减小,夹芯板的后面板最大位移与芯层厚度的比值增大,抗冲击性能变差。

      图  14  不同夹芯结构的芯层塑性耗散能-时间曲线

      Figure 14.  Plastic dissipation energy of core layer-time curves for different sandwich structures

      图  15  不同结构后面板最大位移与芯层厚度的比值曲线

      Figure 15.  The ratio curves of the maximum displacement of back layer to core thickness for different sandwich structures

      图16给出了5种结构的后面板各层纤维的最大应力变化曲线。可以看出,随着芯层相对密度的提高和厚度的减小,各层上的最大应力总体呈现增大趋势。

      图  16  不同夹芯结构后面板各层的最大应力变化曲线

      Figure 16.  Maximum stress of each back layer for different sandwich structures

      5种结构后面板的撕裂程度见表6。不同结构的芯层相对密度增加,芯层高度减小,后面板撕裂层数逐渐增加。

      Structure No.Tear layers
      1#0
      2#3
      3#5
      4#7
      5#8

      表 6  5种不同结构后面板的撕裂程度

      Table 6.  Tear degree of rear panel for five sandwich structures

      在本研究的5种结构中,当芯层相对密度为15.0%时,芯层相对较薄,冲击载荷峰值较小,芯层的塑性耗散能较大,并且后面板的最大位移也相对较小,后面板剩余3层未撕裂,仍然具有一定的承载能力,所以认为在芯层质量保持一定的情况下,结构类型3#为最优的芯层相对密度与厚度耦合设计。

    • 为研究不同边界条件对碳纤维-泡沫铝夹芯板抗冲击性能的影响,使用简支和固支两种边界条件,冲击能量为33.0 J,分别对相同的夹芯板进行冲击,芯层相对密度为10.0%。

      结合图17图18图19可以看出,边界条件为四边简支时,冲击完成所需时间比固支边界更长,锤头位移的峰值更大,后面板的最大应力是固支条件下的1.99倍。简支条件下后面板最大应力大于固支条件,说明四边固支的夹芯板的抗冲击性能更好。

      图  17  冲击载荷-时间曲线

      Figure 17.  Force-time curves

      图  18  锤头位移-时间曲线

      Figure 18.  Displacement-time curves of impactor

      图  19  后面板应力云图

      Figure 19.  Stress plot of rear panel

    • 为研究碳纤维面板-泡沫铝芯层结构抗低速冲击响应,对复合材料夹芯结构受低速冲击载荷进行了有限元仿真,建立了低速冲击T700碳纤维-泡沫铝夹层板的有限元模型,分析了复合材料夹芯结构在落锤冲击下的损伤起始、损伤扩展和最终破坏模式,通过锤头的接触力、位移、夹芯板的内能、后面板的最大位移研究夹层结构的能量吸收情况及抗冲击特性,主要得出如下结论。

      (1) 给出了碳纤维-泡沫铝夹芯板落锤低速冲击下的损伤失效区域,结果显示夹芯结构可以有效地吸收冲击能量。随着冲击能量的增大,夹层结构面板的失效层数增多,冲击载荷的峰值增大,结构塑性变形吸收的能量也增大,后面板的最大位移增大。

      (2) 保持冲击能量不变,设计面板铺层不变、芯层相对密度在10.0%~20.0%之间的5种结构,通过改变芯层的厚度,使芯层质量一定。在这5种芯层相对密度与厚度的耦合关系中,芯层相对密度15.0%、厚度10 mm为最优的耦合方式,使得结构既满足一定的抗冲击性又不至于太厚。为实际工程中平衡防护层的厚度、质量与抗冲击性能的需求提供了设计指导。

      (3) 相同能量冲击下,四周固支的夹芯板的抗冲击性能较四周简支更好。

参考文献 (17)

目录

    /

    返回文章
    返回