Loading [MathJax]/jax/output/HTML-CSS/jax.js
  • ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST、EBSCO、DOAJ收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

基于小波的蜂窝板面超高速撞击声发射信号损伤特征提取

刘源 庞宝君 迟润强 曹武雄 张志远

陈昌鑫, 靳鸿, 马铁华. 冲击加速度存储测试的变频采样策略分析[J]. 爆炸与冲击, 2015, 35(4): 501-506. doi: 10.11883/1001-1455(2015)04-0501-06
引用本文: 刘源, 庞宝君, 迟润强, 曹武雄, 张志远. 基于小波的蜂窝板面超高速撞击声发射信号损伤特征提取[J]. 爆炸与冲击, 2017, 37(5): 785-792. doi: 10.11883/1001-1455(2017)05-0785-08
Chen Chang-xin, Jin Hong, Ma Tie-hua. Analysis of frequency-change sampling strategy for impact acceleration storage measurement[J]. Explosion And Shock Waves, 2015, 35(4): 501-506. doi: 10.11883/1001-1455(2015)04-0501-06
Citation: Liu Yuan, Pang Baojun, Chi Runqiang, Cao Wuxiong, Zhang Zhiyuan. Wavelet transformation based damage feature extraction ofhypervelocity impact acoustic emission signalon honeycomb core sandwich[J]. Explosion And Shock Waves, 2017, 37(5): 785-792. doi: 10.11883/1001-1455(2017)05-0785-08

基于小波的蜂窝板面超高速撞击声发射信号损伤特征提取

doi: 10.11883/1001-1455(2017)05-0785-08
基金项目: 

空间碎片专项十二五项目; 中央高校基本科研业务费专项项目 HIT.NSRIF.2015029

详细信息
    作者简介:

    刘源(1987—),男,博士研究生

    通讯作者:

    迟润强, chirq@hit.edu.cn

  • 中图分类号: O384

Wavelet transformation based damage feature extraction ofhypervelocity impact acoustic emission signalon honeycomb core sandwich

  • 摘要: 为了通过超高速撞击声发射信号识别蜂窝结构受空间碎片撞击后的损伤状态,提出一种基于小波的损伤特征提取方法。采用超高速撞击声发射技术,以铝合金蜂窝板为研究对象,通过超高速撞击实验获取实验信号。分析超高速撞击声发射信号的时频特征及板波模态等特征,采用Daubechies小波变换将信号中模态分离,根据小波系数计算各尺度小波能量分数及小波能量熵特征,分析各特征参数与损伤间的关系,并通过Kruskal-Wallis检验方法验证各特征值对损伤识别的贡献。结果表明:小波能量分数和小波能量熵具有一定的损伤模式分类能力;250 kHz以上的小波能量分数具有良好的损伤模式分类能力;非超声部分的低频信号对损伤识别存在干扰。
  • 舰载设备抗冲击能力是舰艇生命力的重要组成部分,轻级、中级舰载设备用冲击机进行考核,重级舰载设备用浮动冲击平台进行考核[1]。浮动冲击平台作为检验设备抗冲击能力的重要工具,它所能提供的冲击谱谱值大小直接影响了其检验能力的阈值[2-3]。设备安装在平台内底上如图1所示,当受到舷侧方向爆炸载荷作用后,会产生横向、纵向、垂向3个方向的振动。其中纵向振动不在爆炸主要冲击方向,相比横向和垂向振动幅度小很多,由此可见横向和垂向的冲击震动谱值过大是导致设备损坏的主要原因。

    图  1  水下爆炸示意图
    Figure  1.  Underwater explosion sketch

    不同舰载设备有不同的抗冲击要求,检验重级设备抗冲击能力需要使用浮动冲击平台,因为浮动冲击平台能够提供足够严厉的冲击环境。工程上常用设计谱的低频谱位移、中频谱速度和高频谱加速度来量化冲击环境的严厉程度[4],更高的谱值对应着更严厉的冲击环境。对于设计谱的3个谱值,谱速度和谱位移过大往往是导致设备损坏的直接原因[5],因此如何提高平台谱速度、谱位移成为了冲击平台研究领域内的一个核心问题[6]。此外,高频段谱值很大程度上受滤波影响,因此冲击谱高频段不作为研究重点。对于某些设备,不仅要考虑垂向谱速度,横向谱速度也是考核重点,此时应在不会明显影响垂向谱速度的同时尽量提高横向谱速度。改变结构形式是实现目的的一种有效手段,张磊等[7]发现浮动冲击平台的导流装置能够很好的在低频段调节垂向和横向的冲击强度比例。王军等[8]经过四种形式的对比发现外板夹层箱形梁结构形式的平台强度最好,能够提供尽可能高的谱值。通过改变基础的结构形式可以改变结构振动的主振型,也可改变其上冲击谱的整体趋势、峰值频率和峰值高度[9]。本文结合特定中型浮动冲击平台,对平台外形结构进行设计,利用有限元软件研究平台外形结构特征对其冲击环境的影响,重点分析怎样的结构形式能够提高冲击平台横向冲击谱谱速度和谱位移。

    作者参考多篇文献,所用平台模型引用王军[10]的设计,如图23所示,平台外形参数如表1所示。为防止平台进水,需要在上面加一个顶盖,顶盖上覆盖幕布,顶盖质量不可忽略,数值模拟时需考虑顶盖的影响。

    图  2  浮动冲击平台三视图[10]
    Figure  2.  Three views of floating impact platform[10]
    图  3  平台模型图
    Figure  3.  Platform model diagram
    表  1  中型浮动冲击平台外形尺寸
    Table  1.  Size of medium floating impact platform
    长度/m宽度/m本体高度/m顶盖高度/m内底高度/m质量/t吃水深度/m
    12.26.14.223.051.092.271.33
    下载: 导出CSV 
    | 显示表格

    材料选用强度较高的921A钢。平台设备放置于内底上,针对各种设备质量,利用ANSYS建出不同深度流场,在ABAQUS中完成组装,组装完成如图4所示。

    图  4  平台流场装配图
    Figure  4.  Assembly drawing of platform flow field

    水下爆炸研究通常采取arbitrary Lagrange-Euler (ALE)法或声固耦合法[11]。ALE法要求流场包含炸药在内,且对网格要求较高,计算近场爆炸也能得到较准确的结果,但对于远场爆炸,其网格数量过多会导致计算时间变得很长甚至直接计算失败。而声固耦合法通过模拟声学介质和边界无反射技术允许流场不包含炸药,对于远场爆炸可以得到较准确的结果,相比ALE方法花费时间更少、计算更稳定[12]。本文涉及的水下爆炸工况均属远场爆炸,因此选用声固耦合法对平台冲击响应进行分析。

    声音在弹性介质中的传播体现为振动的传播,而振动在流体介质中的传播体现为压缩波的传播,即声音在水中传播表现为纵波,衰减很慢。实际海洋流场是无限大的,不存在声波反射问题,而有限元软件中的流场模型是有限大,流场边界对声波的反射会导致计算结果产生错误,因此声固耦合法利用边界阻抗技术模拟无限流畅,即无反射边界条件。

    本计算方法已在文献[11]中加以验证,结果显示声固耦合法能够比较准确的模拟舰船遭受水下爆炸冲击时产生的响应。本文后续部分将用这种数值方法对一特定中型浮动冲击平台的冲击响应特性进行研究。

    中型浮动冲击平台最大考核设备质量为113 t,据分析计算可知设备越轻,平台冲击响应越大,冲击谱谱值越大,为使冲击谱的变化更明显,运用文献[4]工况中最近的工况,即装药55 kg TNT、水深8.2 m、舷侧水平距离6.71 m(参考图1),安装10 t设备采用一体化分析方法[13]进行数值模拟。仿真时间设为0.25 s,已经包含了整个冲击波作用阶段。

    在结果文件中提取出内底板上如图5所示3个纵横桁材交叉点上的时历加速度响应。所提取的时历响应含有很多高频成分,而实际工程中,设备安装频率一般较低,普遍低于250 Hz。时历响应中的高频段成分对其影响很小,因此需要对初始响应进行滤波操作,把研究重点放在中低频。参考实际舰船冲击环境及其他相关领域研究材料,本文采用250 Hz低通滤波,滤波前后平台内底时历加速度响应如图6所示。

    图  5  内底测点
    Figure  5.  Inner bottom point
    图  6  250 Hz滤波前后平台内底时历加速度响应
    Figure  6.  The acceleration response of the front and rear platform of the 250 Hz filter is obtained

    图中红色曲线是原始时历曲线,蓝色曲线是250 Hz低通滤波之后的时历曲线,用此曲线绘制冲击谱。利用自主开发程序首先根据冲击谱理论完成Duhamel积分得到原始冲击谱,然后将原始谱圆整为设计谱,即可得到各测点处谱位移、谱速度和谱加速度。

    浮动冲击平台通过水面以下的湿表面接收冲击波载荷[14],在质量不变的情况下,受到的力越大,则冲击谱谱值越大,因此考虑如何在不明显改变平台质量的情况下增加横向载荷的接收。平台所受的水下爆炸冲击波载荷大小与湿表面积密切相关,因此考虑在平台舷侧底端迎爆面加装竖直挡板,在不明显改变平台质量的情况下增加横向湿表面面积,从而提高横向载荷的接收,达到提高谱值的目的。竖直挡板位于舷侧下方,厚0.032 m,长1 m,如图7所示。中型浮动冲击平台最大考核设备质量为113 t,据分析计算知设备越重,平台距药包越近,塑性变形越大,现将流场和113 t设备组装好进行安全性分析,检验挡板结构能否满足平台考核设备时的强度要求,结果如图8所示。

    图  7  舷侧竖直挡板
    Figure  7.  Side vertical baffle
    图  8  加舷侧挡板后塑性应变
    Figure  8.  Plastic strain after adding side baffle

    图8可以看出,载荷作用时间结束后等效塑性应变为0.038,满足强度要求,可展开下一步研究。在此结构上施加不变的工况,分析此结构内底上的横向冲击谱,与相同载荷下无挡板结构横向冲击谱进行对比。为保证冲击载荷不变,药包水深随着吃水改变而改变,保证爆距不变,冲击因子不变[15]。取图5中的平台内底上CLR三个测点,横向冲击谱(谱位移u、谱速度v、谱加速度a)如图9所示。由图可以看出,能略微提高横向谱速度值但效果并不明显,下一节分析其主要原因。

    图  9  舷侧竖直挡板与无挡板的各测点水平冲击谱对比
    Figure  9.  Horizontal shock spectrum with and without vertical baffle plates

    图9,竖直挡板受某些因素影响不能发挥预想的效果,查阅相关文献资料,将主要影响因素归结为以下两点。

    (1)冲击波的绕射:当冲击波被一个小于其波长的物体阻挡,就会绕过这个物体继续行进[16]。挡板尺寸远小于冲击波波长,冲击波绕过挡板两侧和底部后会产生指向挡板背面的力,与挡板正面的力相抵消,如图10所示,由于挡板厚度很小,绕射波不会有很大的衰减,即增加舷侧挡板不能明显增加横向载荷的接收。

    图  10  竖直挡板冲击波绕射现象
    Figure  10.  Vertical baffle wave diffraction phenomenon

    (2)阻力的增加:竖直挡板在船体做横向振动时无疑会增加阻力,阻碍振动,降低谱速度。

    综合以上两点,如图9所示,竖直挡板几乎不能提高冲击平台谱速度。冲击波绕射现象广泛存在于水下爆炸作用过程中,与阻力一样,二者都是不可避免的,但可以通过合理的结构设计进行尽量削弱其影响。图10只给出了增加的竖直挡板处的冲击波绕射现象,但实际上平台本体也一样受到绕射波的作用,如图11所示。但由图9所示,平台依然由很大的横向振动,这是由于平台横向尺寸很大绕射波传到平台另一侧后有明显衰减。考虑到这一点,下一节将对平台外形结构进行优化,以尽量削弱冲击波绕射和阻力增加的影响。

    图  11  平台本体冲击波绕射现象
    Figure  11.  Platform body wave diffraction phenomenon.

    针对竖直挡板存在的问题,作者提出流线型挡板优化方案,为方便对比分析,取图12图13所示大小不同两个尺寸,其为参考文献[17]中流线型结构进行设计,装配图如图14图15所示。

    图  12  挡板1截面线型
    Figure  12.  Baffle 1 section type
    图  13  挡板2截面线型
    Figure  13.  Baffle 2 section type
    图  14  挡板1模型图
    Figure  14.  Baffle 1 model diagram
    图  15  挡板2模型图
    Figure  15.  Baffle 2 model diagram

    当冲击波从迎流一侧传来时,会从挡板底部绕到挡板另一面并产生指向挡板的压力,挡板表面受力如图16所示。

    图  16  底部流线挡板冲击波绕射现象
    Figure  16.  The bottom streamline baffle wave diffraction phenomenon

    在同一高度微元dh处冲击波水平方向合力为:

    F=p1dL1cosθ1p2dL2cosθ2=(p1p2)dh (1)

    式中:θ1p1方向与水平方向夹角;θ2p2方向与水平方向夹角;L1为dh对应的迎流面长度;L2为dh对应的背流面长度。由式(1),F仅与冲击压力波衰减有关,而冲击波在流线挡板两侧的衰减要比竖直挡板的正反两面明显得多,流线型挡板所受合力显然更大。并且采用流线造型可以大幅度的降低挡板所受阻力,减少阻力对振动的阻碍作用,流线型挡板设计很好地解决了竖直挡板存在的两个问题。

    流线型挡板重心不在平台中心处,在其重心对称位置沿长度方向一条线上加上若干压载,使重心回到平台中心处,并将爆点位置随平台吃水的改变竖直移动,保证与平台相对位置不变。对流线型挡板进行安全性分析,如图17图18所示。由图可以看出,两种流线挡板等效塑性应变分别为0.029 3和0.035 1,满足强度要求。

    图  17  挡板1塑性变形
    Figure  17.  Plastic deformation of baffle 1
    图  18  挡板2塑性变形
    Figure  18.  Plastic deformation of baffle 2

    在两种结构上施加与之前相同的工况,分析内底上CLR三个测点冲击谱,与相应载荷下无挡板结构垂向及横向冲击谱进行对比,如图19所示。由图可以看出流线型挡板对平台垂向冲击谱几乎无影响,但却可以明显提高平台横向冲击谱中低频段谱值。对于无挡板、挡板1、挡板2分别取CLR三测点设计谱谱速度和谱位移进行对比,如表2所示。挡板2由于尺寸比挡板1大一些,对横向载荷接收较大,相应的平台冲击谱谱值也比较大。

    图  19  流线挡板结构与无挡板结构各测点冲击谱对比
    Figure  19.  Comparison of the impact spectrum between the flow baffle structure and the non-baffle structure
    表  2  有无挡板4测点谱值平均值对比
    Table  2.  comparison of the mean value of the four measuring points
    结构垂向谱速度/(m∙s−1)垂向谱位移/cm横向谱速度/(m∙s−1)横向谱位移/cm
    无挡板2.69811.810.8463.08
    挡板12.61611.941.0783.22
    挡板22.77312.261.2383.65
    下载: 导出CSV 
    | 显示表格

    表2看出,流线型挡板可以在几乎不影响(影响小于4%)平台垂向谱速度和谱位移的同时提高平台横向谱速度和谱位移。挡板2可以将横向谱速度提高46.3%,谱位移提高18.5%,可见效果十分显著。当待考核舰载设备对横向冲击环境有较高要求的时候,应考虑添加流线型挡板。而流线挡板几乎不会影响垂向冲击环境,这就说明在设计好内底隔板数量和决定是否使用甲板模拟器之后,可以再进行挡板的考虑,而不会干扰到内底隔板和甲板模拟器对冲击谱的影响。

    本文对浮动冲击平台外形结构进行了结构设计和优化,研究了外形结构特征对其冲击环境的影响并提出了提高横向冲击谱谱值的方法。由于冲击谱谱加速度很大程度受滤波影响,平台结构特征对冲击环境的影响主要体现在冲击谱谱位移和谱速度上。研究得出以下结论:(1)增加舷侧竖直挡板对垂向和横向冲击谱均没有明显影响;(2)增加底部流线型挡板对对垂向冲击谱影响很小,小于4%;(3)增加底部流线型挡板对横向冲击谱有很大影响,能明显提高谱值。随着流线型挡板尺寸变大,对横向冲击谱的影响也会变大,对于美军标中最近工况,采用图13所示挡板2结构,可使横向冲击谱谱位移提高18.5%,谱速度提高46.3%。而何种尺寸挡板形式最优,高度应如何选取最佳等相关问题将在之后的工作中进一步探索和研究。

  • 图  1  时间窗选取示意图

    Figure  1.  Schematic diagram oftime window selection

    图  2  传感器装置示意图

    Figure  2.  Sketch of sensor positions on target plate

    图  3  蜂窝板弹道极限曲线

    Figure  3.  Ballistic limit curve ofhoneycomb core sandwich

    图  4  蜂窝板后蒙皮面板超高速撞击损伤情况

    Figure  4.  Hypervelocity impact damage of back panelsof honeycomb core sandwich

    图  5  蜂窝板的超高速撞击声发射信号波形

    Figure  5.  Oscillograph of hypervelocity impact acoustic emission signalin honeycomb core sandwich

    图  6  高通滤波信号

    Figure  6.  Signal by high-pass filtering

    图  7  信号模态

    Figure  7.  Modal acoustic emission of signal

    图  8  各频带能量分数与弹丸速度关系

    Figure  8.  Relationship between energy fraction of every frequency band and projectile velocity

    图  9  小波能量熵与弹丸速度关系

    Figure  9.  Relationship between wavelet energy entropyand projectile velocity

    表  1  信号频带划分

    Table  1.   Frequency band division of signal

    尺度 小波系数 频带范围/kHz
    1 D1(n) 0~62.5
    2 D2(n) 62.5~125
    3 D3(n) 125~250
    4 D4(n) 250~500
    5 D5(n) 500~1 000
    下载: 导出CSV

    表  2  Kruskal-Wallis检验结果

    Table  2.   Results of Kruskal-Wallis test

    参数 κKW 显著性水平
    D1(n) 0.115 0.735
    D2(n) 3.206 0.073
    D3(n) 3.304 0.069
    D4(n) 24.815 0
    D5(n) 35.672 0
    H 3.268 0.071
    下载: 导出CSV
  • [1] Upper stage explosion places LEO satellites at risk[J/OL]. Orbital Debris Quarterly News, 2013, 17(1): 8[2016-04-24]. https://orbitaldebris.jsc.nasa.gov/quarterly-news/pdfs/odqnv17i1.pdf.
    [2] Liou J C, Giovane F, Corsaro R, et al. LAD-C: A large area debris collector on the ISS[C]//36th COSPAR Scientific Assembly. Beijing, China, 2006: 36.
    [3] Prosser W H, Madaras E I. Distributed impact detector system (DIDS) health monitoring system evaluation[R]. Hampton: Langley Research Center, 2010.
    [4] Spencer G, Schäfer F, Tanaka M, et al. Design and initial calibration of micrometeoroid space debris detector (MDD)[C]//Proceedings of the 4th European Conference on Space Debris. Darmstadt, Germany: ESA, 2005: 18-20.
    [5] 刘武刚, 庞宝君, 韩增尧, 等.基于小波分析技术的高速撞击声发射源定位[J].高技术通讯, 2009, 19(2):181-187. http://d.old.wanfangdata.com.cn/Periodical/gjstx98200902013

    Liu Wugang, Pang Baojun, Han Zengyao, et al. Acoustic emission detection and location for hypervelocity impacts based on wavelet transform[J]. Chinese High Technology Letters, 2009, 19(2):181-187. http://d.old.wanfangdata.com.cn/Periodical/gjstx98200902013
    [6] Liu Z D, Pang B J. A method based on acoustic emission for locating debris cloud impact[C]//4th International Conference on Experimental Mechanics. 2009: 7522.
    [7] 张凯, 庞宝君, 林敏.碎片云撞击声发射信号能量特征小波包分析[J].振动与冲击, 2012, 31(12):125-128. http://d.old.wanfangdata.com.cn/Periodical/zdycj201212025

    Zhang Kai, Pang Baojun, Lin Min. Wavelet packet analysis for acoustic emission signals caused by debris cloud impact[J]. Journal of Vibration and Shock, 2012, 31(12):125-128. http://d.old.wanfangdata.com.cn/Periodical/zdycj201212025
    [8] 熊秋鹏.基于神经网络技术的空间碎片损伤模式识别研究[D].哈尔滨: 哈尔滨工业大学, 2012.
    [9] Taylar E M, Glanville J P, Clegg R A, et al. Hypervelocity impact on spacecraft honeycomb: Hydrocode simulation and damage laws[J]. International Journal of Impact Engineering, 2003, 29(1):691-702.
    [10] Nia A A, Razavi S B, Majzoobi H H. Ballistic limit determination of aluminum honeycombs: Experimental study[J]. Materials Science and Engineering A, 2008, 488(1):273-280. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ027439541/
    [11] Liu Y, Pang B J, Jia B, et al. Modal acoustic emission based location method in honeycomb core sandwich structure[C]//Sixth European Conference on Space Debris. 2013: 183.
    [12] Chakraborty N, Rathod V T, Mahapatra D R, et al. Guided wave based detection of damage in honeycomb core sandwich structures[J]. NDT & E International, 2012, 47(7):27-33.
    [13] 唐颀.高速撞击板波特性与声发射空间碎片在轨感知技术[D].哈尔滨: 哈尔滨工业大学, 2008. http://cdmd.cnki.com.cn/Article/CDMD-10213-2009223857.htm
    [14] 刘武刚.基于声发射的空间碎片在轨撞击感知技术研究[D].哈尔滨: 哈尔滨工业大学, 2008.
    [15] Rosso O A, Blanco S, Yordanova J, et al. Wavelet entropy: A new tool for analysis of short duration brain electrical signals[J]. Journal of Neurosci Methods, 2001, 105(1):65-75. doi: 10.1016/S0165-0270(00)00356-3
    [16] 向明江.屏舱声发射信号耦合对在轨感知系统定位的影响研究[D].哈尔滨: 哈尔滨工业大学, 2011. http://cdmd.cnki.com.cn/Article/CDMD-10213-1012002331.htm
  • 期刊类型引用(2)

    1. 余孙全,樊程广,张翔,付康佳,陈勇,綦磊. 航天器结构中导波健康监测技术的若干进展. 宇航学报. 2024(04): 487-498 . 百度学术
    2. 吕伟臻,宋燕,黄雪刚,殷春. 一种改进的特征匹配算法在碎片云测量建模中的应用. 振动与冲击. 2021(02): 29-38 . 百度学术

    其他类型引用(4)

  • 加载中
推荐阅读
剪切增稠液填充蜂窝夹芯板的低速冲击响应
李雨薇 等, 爆炸与冲击, 2025
高速冲击下混凝土动力学性质和动态温度研究
黄晨瑞 等, 爆炸与冲击, 2025
超高速撞击条件下混凝土靶体内 应力波的测量和分析
钱秉文 等, 爆炸与冲击, 2025
考虑动态拉压比影响的岩石损伤本构模型
胡学龙 等, 爆炸与冲击, 2025
聚酰亚胺的动静态压缩力学性能及本构模型
于文峰 等, 高压物理学报, 2022
玻璃纤维增强聚碳酸酯复合材料的动态拉伸力学行为特征及其失效机理
管海陆 等, 高压物理学报, 2023
高速颗粒群冲击多层夹芯复合结构的动态力学响应特性研究进展
郑伟 等, 高压物理学报, 2025
Quaternary ammonium salts: insights into synthesis and new directions in antibacterial applications
Zhou, Zhenyang et al., BIOCONJUGATE CHEMISTRY, 2023
Critical dynamic stress and cumulative plastic deformation of calcareous sand filler based on shakedown theory
JOURNAL OF MARINE SCIENCE AND ENGINEERING, 2023
Performance and geocell-soil interaction of sand subgrade reinforced with high-density polyethylene, polyester, and polymer-blend geocells: 3d numerical studies
COMPUTERS AND GEOTECHNICS
Powered by
图(9) / 表(2)
计量
  • 文章访问数:  4668
  • HTML全文浏览量:  1342
  • PDF下载量:  465
  • 被引次数: 6
出版历程
  • 收稿日期:  2016-01-20
  • 修回日期:  2016-04-24
  • 刊出日期:  2017-09-25

目录

/

返回文章
返回