碎片云超高速撞击声发射信号特征分析

庞宝君 张凯 林敏 刘源

武晓东, 张海广, 王瑜, 孟祥生. 冲击载荷下仿贝壳珍珠层Voronoi结构的动态力学响应[J]. 高压物理学报, 2020, 34(6): 064201. doi: 10.11858/gywlxb.20200559
引用本文: 庞宝君, 张凯, 林敏, 刘源. 碎片云超高速撞击声发射信号特征分析[J]. 高压物理学报, 2014, 28(6): 664-670. doi: 10.11858/gywlxb.2014.06.004
WU Xiaodong, ZHANG Haiguang, WANG Yu, MENG Xiangsheng. Dynamic Responses of Nare-Like Voronoi Structure under Impact Loading[J]. Chinese Journal of High Pressure Physics, 2020, 34(6): 064201. doi: 10.11858/gywlxb.20200559
Citation: PANG Bao-Jun, ZHANG Kai, LIN Min, LIU Yuan. Characteristic Analysis of Acoustic Emission Signals Caused by Debris Cloud Impact[J]. Chinese Journal of High Pressure Physics, 2014, 28(6): 664-670. doi: 10.11858/gywlxb.2014.06.004

碎片云超高速撞击声发射信号特征分析

doi: 10.11858/gywlxb.2014.06.004
基金项目: 国家空间碎片专题研究项目(kjsp06211)
详细信息
    作者简介:

    庞宝君(1963—), 男,教授,博士生导师,主要从事空间碎片超高速撞击及航天器防护技术研究. E-mail:pangbj@hit.edu.cn

    通讯作者:

    张凯(1984-), 男, 博士研究生, 主要从事空间碎片环境下航天器在轨健康监测研究.E-mail:bluerobbin@163.com

  • 中图分类号: O348.8

Characteristic Analysis of Acoustic Emission Signals Caused by Debris Cloud Impact

  • 摘要: 为了掌握带防护屏的航天器结构受空间碎片超高速撞击时的声发射信号特征,利用二级轻气炮发射球形弹丸撞击铝合金双层板结构,获取了碎片云撞击铝合金板舱壁产生的声发射信号,并利用小波包技术和能量熵理论对信号进行了分析。实验结果表明:弹丸初始速度、防护屏厚度及弹丸直径是决定二次碎片云形态及声发射信号特征的重要因素;在本实验工况范围内,小波包能量熵值能够描述声发射信号频率的复杂程度;当弹丸初始速度处于破碎段(3~7 km/s)时,随着初始速度的增大,二次碎片云进一步细化,撞击产生的声发射信号幅值趋于减小、频率成分趋于复杂化,其小波包能量熵值逐渐增大;防护屏厚度对声发射信号的小波包能量熵值影响较大,弹丸直径对其影响较小。研究结果有助于实现对碎片云撞击舱壁结构的损伤模式识别。

     

  • 仿贝壳珍珠层复合材料作为新型复合材料,有着优异的力学性能,不仅具有很高的强度,而且具有很好的韧性,近年来引起了学术界的不断关注[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)所示,站点位置(x,y)由极坐标控制方程决定

    图  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=x0+rcosθ,y=y0+rsinθ
    (1)

    式中:(x0,y0)为网格内站点的参考位置;r(0,R)之间的随机值,R为圆的半径;θ(0,2π)之间的随机值。通过站点的随机移动,最终形成新的随机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={tn,ts,tt},其中下标n、s和t分别表示一个法向和两个切向分量。这些变量之间满足双线性二次黏聚律[25]

    (tnt0n)2+(tst0s)2+(ttt0t)2=1
    (2)

    式中:t0nt0st0t分别表示变形垂直于界面以及在第一、第二剪切方向上的最大应力。

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

    (GnG0n)2+(GsG0s)2+(GtG0t)2=1
    (3)

    式中:GnGsGt分别为正向和两个切向的断裂能,G0nG0sG0t分别为正向和两个切向引起破坏所需的最大断裂能。

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

    表  1  铝片的材料参数
    Table  1.  Parameters of aluminum plate
    Materialρ/(kg·m−3νE/GPa
    Aluminum27500.372
    下载: 导出CSV 
    | 显示表格
    表  2  Cohesive模型的材料参数
    Table  2.  Parameters of Cohesive model
    t0n/MPat0s/MPat0t/MPaG0n/(kJ·m−2)G0s/(kJ·m−2)G0t/(kJ·m−2)ρ/(kg·m−3)Es/GPaEt/GPa
    8080801111 85041.5
    下载: 导出CSV 
    | 显示表格

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

    σ=(A+Bεn)(1+Cln˙ε˙ε0)[1(TTrTTm)m]
    (4)

    式中:σ为等效应力,A为材料在常温准静态下的屈服强度,B为应变强化系数,n为强化指数,ε为等效塑性应变,˙ε为塑性应变率,˙ε0为初始应变率,C为应变率敏感系数,m为温度软化系数,Tm为材料的熔点,Tr为参考温度(取常温),T为当前温度。表3列出了Johnson-Cook模型参数[25]

    表  3  Johnson-Cook模型参数[25]
    Table  3.  Parameters of Johnson-Cook model[25]
    MaterialA/MPaB/MPanCm
    Aluminum3916840.4360.009592
    下载: 导出CSV 
    | 显示表格

    图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模型的数值模拟工作进行验证,进而开展更深入的研究。

  • 图  碎片云撞击信号获取实验示意图

    Figure  1.  Schematic of signal acquisition experiments for debris cloud impact

    图  实验结果与预测曲线的比较

    Figure  2.  Comparison between experimental results and predicted curve from Eq.(1)

    图  不同破碎程度碎片云撞击损伤结果

    Figure  3.  Impact damage caused by different levels of fragmentation projectiles

    图  弹丸不同程度破碎形成的碎片云撞击声发射信号波形

    Figure  4.  Debris cloud impact signals by projectiles with different fragmentation degrees

    图  小波包能量熵值随弹丸初始速度的变化

    Figure  5.  Wavelet packet energy entropy changing with initial velocity of projectile

    图  距离板中心不同位置处信号的小波包能量熵值

    Figure  6.  Wavelet packet energy entropy versus velocity at different positions

    表  1  实验工况及结果

    Table  1.   Experimental conditions and results

    No. dp/
    (mm)
    tb/
    (mm)
    v/
    (km/s)
    Dmax/
    (mm)
    Damage of
    the rear wall
    1 3.20 0.5 2.78 110 P
    2 3.20 0.5 3.25 120 P
    3 3.20 0.5 3.52 129 P
    4 3.20 0.5 4.24 142 P
    5 3.20 0.5 4.63 150 P
    6 3.20 0.5 4.92 158 P
    7 3.20 1.0 2.50 102 NP
    8 3.20 1.0 3.05 126 P
    9 3.20 1.0 3.65 133 NP
    10 3.20 1.0 4.00 152 NP
    11 3.20 1.0 4.45 159 NP
    12 3.20 1.0 5.05 165 NP
    13 3.20 1.5 2.81 120 NP
    14 3.20 1.5 3.05 130 NP
    15 3.20 1.5 3.73 148 NP
    16 3.20 1.5 4.21 160 NP
    17 3.20 1.5 4.46 165 NP
    18 3.20 1.5 4.81 170 NP
    19 3.97 1.0 2.65 120 P
    20 3.97 1.0 3.25 138 P
    21 3.97 1.0 3.57 148 P
    22 3.97 1.0 4.15 160 NP
    23 3.97 1.0 4.55 170 NP
    24 3.97 1.0 4.85 180 NP
    25 4.76 1.0 2.31 110 P
    26 4.76 1.0 3.01 140 P
    27 4.76 1.0 3.30 145 P
    28 4.76 1.0 3.75 160 P
    29 4.76 1.0 4.64 180 P
    30 4.76 1.0 5.20 195 P
    Note: P-Penetration in the rear wall, NP-No perforated point in the rear wall.
    下载: 导出CSV
  • [1] Miyachi T, Hasebe N, Ito H, et al. Real-time detector for hypervelocity microparticles using piezoelectric material[J]. Adv Space Res, 2004, 34(5): 935-938. doi: 10.1016/j.asr.2003.11.019
    [2] Schäfer F, Janovsky R. Impact sensor network for detection of hypervelocity impacts on spacecraft[J]. Acta Astronaut, 2007, 61(10): 901-911. doi: 10.1016/j.actaastro.2007.02.002
    [3] Prosser W H, Gorman M R, Humes D H. Acoustic emission signals in thin plates produced by impact damage[J]. J Acoust Emiss, 1999, 17(1/2): 29-36.
    [4] 唐颀, 庞宝君, 韩增尧, 等.单层板超高速撞击声发射波的频谱特征分析[J].宇航学报, 2007, 28(4): 1059-1064.

    Tang Q, Pang B J, Han Z Y, et al. Analysis of frequency spectrum character of acoustic emission wave from hypervelocity impact on single-sheet plate[J]. Journal of Astronautics, 2007, 28(4): 1059-1064. (in Chinese)
    [5] 刘武刚, 庞宝君, 韩增尧, 等.基于声发射的单层铝板高速撞击损伤类型识别[J].宇航学报, 2011, 32(3): 671-675.

    Liu W G, Pang B J, Han Z Y, et al. Damage identification of single aluminum plate produced by hypervelocity impact based acoustic emission[J]. Journal of Astronautics, 2011, 32(3): 671-675. (in Chinese)
    [6] 管永红, 胡八一, 黄超.基于小波包的爆炸容器振动分析[J].爆炸与冲击, 2010, 30(5): 551-555.

    Guan Y H, Hu B Y, Huang C. Vibration analysis of an explosion vessel based on wavelet packet transform[J]. Explosion and Shock Waves, 2010, 30(5): 551-555. (in Chinese)
    [7] 凌同华, 廖艳程, 张胜.冲击荷载下岩石声发射信号能量特征的小波包分析[J].振动与冲击, 2010, 29(10): 127-130, 255.

    Ling T H, Liao Y C, Zhang S. Application of wavelet packet method in frequency band energy distribution of rock acoustic emission signals under impact loading[J]. Journal of Vibration and Shock, 2010, 29(10): 127-130, 255. (in Chinese)
    [8] 管公顺, 庞宝君, 哈跃, 等.铝合金Whipple防护结构高速撞击实验研究[J].爆炸与冲击, 2005, 25(5): 461-466.

    Guan G G, Pang B J, Ha Y, et al. Experimental investigation of high-velocity impact on aluminum alloy Whipple shield[J]. Explosion and Shock Waves, 2005, 25(5): 461-466. (in Chinese)
    [9] Christiansen E L. Design and performance equations for advanced meteoroid and debris shields[J]. Int J Impact Eng, 1993, 14(1): 145-156.
    [10] 唐颀.超高速撞击板波特性与声发射空间碎片在轨感知技术[D].哈尔滨: 哈尔滨工业大学, 2008: 56-76.

    Tang Q. Characteristics of plate waves induced by hypervelocity impact and onboard monitoring technique for detection of impact on spacecraft by space debris[D]. Harbin: Harbin Institute of Technology, 2008: 56-76. (in Chinese)
    [11] Rosso O A, Blanco S, Yordanova J, et al. Wavelet entropy: A new tool for analysis of short duration brain electrical signals[J]. J Neurosci Meth, 2001, 105(1): 65-75. doi: 10.1016/S0165-0270(00)00356-3
    [12] 印欣运, 何永勇, 彭志科, 等.小波熵及其在状态趋势分析中的应用[J].振动工程学报, 2004, 17(2): 49-53.

    Yin X Y, He Y Y, Peng Z K, et al. Study on wavelet entropy and its applications in trend analysis[J]. Journal of Vibration Engineering, 2004, 17(2): 49-53. (in Chinese)
  • 加载中
图(6) / 表(1)
计量
  • 文章访问数:  6304
  • HTML全文浏览量:  2181
  • PDF下载量:  353
出版历程
  • 收稿日期:  2012-12-26
  • 修回日期:  2013-03-08

目录

/

返回文章
返回