下载文档原格式
综述CHINA SYNTHETIC RESIN AND PLASTICS合 成 树 脂 及 塑 料 ， 2022， 39（4）： 70随着现代科技的飞速发展，对高性能材料的需求日益增加，聚四氟乙烯（PTFE）作为一种性能优良的工程塑料，在许多领域具有广泛的应用［1-2］。
PTFE是由单体四氟乙烯聚合而成［3］，分子结构为一种螺旋构象，即C—C骨架全部被周围的F原子包裹。
同时由于C—F的键能很高不易断裂，使PTFE可以抵抗强酸、强碱、油脂、纯氧化剂和有机溶剂等的腐蚀，但缺点是强度较低，不利于成型加工，机械磨损率高，特别是在受外力作用下会产生严重的蠕变现象，极大地限制了PTFE 的应用。
因此对PTFE的改性显得尤为重要［2］。
目前，PTFE的改性方法主要有表面改性、填充改性和共混改性。
本文详细阐述了PTFE改性的几种方法，并研究了改性方法对PTFE复合材料力学性能、摩擦性能和介电性能的影响。
DOI：10.19825/j.issn.1002-1396.2022.04.15 *1 PTFE的改性1.1 表面改性由于PTFE表面结合能较小，不易与其他化合物和小分子反应，同时其他填料也很难附着在PTFE表面。
采用物理化学法对PTFE表面进行处理，可以在PTFE表面产生反应位点同时提高表面的粗糙程度，改善PTFE表面的疏水性、亲核性和防污性能。
常见的处理方法主要有等离子体处理法、电子辐照处理法、偶联剂处理法［4］。
聚四氟乙烯改性现状及研究进展左 程1，肖 伟2*（1. 江苏扬建集团有限公司 扬州华正建筑工程质量检测有限公司，江苏 扬州 202105；2. 上海工程技术大学 数理与统计学院，上海 201620）摘要：综述了近几年国内外聚四氟乙烯（PTFE）改性的研究进展，并总结了表面改性、填充改性和共混改性的优缺点，着重分析了填料对PTFE力学性能、摩擦性能和介电性能的影响。
最后对PTFE改性工艺的发展趋势和前景进行了展望。
聚四氟乙烯结构简式简介聚四氟乙烯（Polytetrafluoroethylene，简称PTFE）是一种常用的高性能聚合物材料，具有许多优异的特性，如耐高温、耐腐蚀、低摩擦系数等。
它在工业领域被广泛应用于润滑剂、密封材料、电子器件等方面。
本文将对聚四氟乙烯的结构简式进行详细介绍。
分子结构聚四氟乙烯是由四氟乙烯单体（C2F4）通过聚合反应形成的线性高分子化合物。
它的化学式为(-CF2-CF2-)n，其中n表示重复单元的个数。
每个重复单元由两个碳原子和四个氟原子组成。
从分子结构上看，聚四氟乙烯具有以下特点：1.长链结构：聚四氟乙烯由大量重复单元组成，形成了长链结构。
这使得其具有良好的柔韧性和拉伸强度。
2.全氟化合物：聚四氟乙烯中的所有氢原子被氟原子取代，使其成为全氟化合物。
这使得聚四氟乙烯具有出色的耐腐蚀性和化学稳定性。
3.线性排列：聚四氟乙烯中的重复单元呈线性排列，没有支链或交联结构。
这使得聚四氟乙烯具有低摩擦系数和良好的电绝缘性能。
物理性质熔点和热稳定性聚四氟乙烯具有很高的熔点，约为327℃。
它在高温下也能保持较好的稳定性，可以长时间使用在260℃以下的温度范围内。
密度聚四氟乙烯的密度相对较低，约为2.2g/cm³。
这使得它成为一种轻质材料，在航空航天等领域具有广泛应用。
机械性能由于聚四氟乙烯分子链之间存在相互作用力较弱，因此其机械强度相对较低。
但是，它具有良好的耐磨性和抗冲击性，能够承受较大的外力。
电绝缘性能聚四氟乙烯是一种优异的电绝缘材料，具有很高的电阻率和低介电常数。
它在高频率和高温下仍能保持良好的绝缘性能。
摩擦系数聚四氟乙烯具有非常低的摩擦系数，约为0.04。
这使得它成为一种理想的润滑材料，在机械设备中广泛应用。
应用领域聚四氟乙烯由于其优异的特性，在许多领域得到了广泛应用，包括但不限于以下几个方面：1.密封材料：由于聚四氟乙烯具有出色的耐腐蚀性和化学稳定性，被广泛应用于密封件制造。
乙烯-四氟乙烯共聚物复合体系的结构与性能研究乙烯-四氟乙烯共聚物(ETFE)是一种由乙烯和四氟乙烯共聚而成的半结晶性含氟聚合物,由于含有四氟乙烯单体,使其具有优异的耐腐蚀、耐候、耐高温、耐老化、耐辐照等性能。
ETFE因其质轻、透明度高、环保等特点,可加工成涂层、薄膜、线缆等,被广泛的应用到航空航天、建筑、农业、化工、井下等诸多领域。
本文主要围绕ETFE材料在加工性能等方面的不足,通过引入同为含氟材料的聚偏氟乙烯(PVDF)对其进行共混改性,制备了 ETFE/PVDF二元复合材料。
通过对该复合材料的相结构、相容性、热性能、流变性能、机械性能等进行表征分析,从而建立起复合体系的微观结构与宏观性能间的对应关系,深入地对ETFE材料研究探索。
此外,通过电子束辐照以及引入无机蒙脱土(MMT)纳米粒子等改性方式对ETFE/PVDF复合体系进行二次改性,逐步建立起以ETFE为主体材料的综合加工应用体系。
本论文的主要研究内容以及相应的实验结论如下:1.采用熔融共混的方式制备了不同比例的ETFE/PVDF二元复合体系。
在空气环境下,ETFE/PVDF(90/10 wt%)复合体系具有更高的初始热分解温度(374.5℃)和初始热分解活化能(187 kJ/mol),分别比纯ETFE组分高52.7℃和121 kJ/mol。
随着PVDF含量的增加,ETFE/PVDF复合体系的储能模量和复合黏度逐渐增加,复合体系的剪切变稀现象增强;PVDF的引入,改善了ETFE的挤出加工性能,ETFE/PVDF复合体系仍具有很好的拉伸强度和断裂伸长率。
2.采用熔融共混的方式制备了PVDF、ETFE与不同敏化剂组成的复合材料,再通过电子束辐照对复合体系进行改性处理。
选用的敏化剂分别是三烯丙基异氰酸酯(Triallyl Isocyanurate,TAIC)、三羟甲基丙烷三甲基丙烯酸酯(Trimethylolpropane Trimethacrylate,TMPTMA)和1、2、4-苯三甲酸三烯丙酯(Trially trimellitate)。
王靖岩，王义智，韩志伟铝‐氟聚物反应性物质制备及研究进展王靖岩1，王义智2，韩志伟1（1.南京理工大学化工学院，江苏南京210094；2.南京长山化工有限责任公司，江苏南京211103）摘要：综述了物理混合法、球磨法、气相沉积法、静电喷雾/纺丝法、溶剂/非溶剂法、3D 打印法六大类制备方法，从产品性能、方法的优缺点等角度对近些年来铝⁃氟聚物反应性物质的研究进展进行简要综述。
介绍了铝⁃氟聚物反应性物质在慢升温速率和快升温速率下的反应过程。
指出今后研究的重点方向为：设计一种能集各种制备方法优点于一身的新方法，加强铝⁃氟聚物反应性物质受热时氟聚物和铝的反应机理的研究。
关键词：铝粉；氟聚物；反应性物质；制备方法；反应过程中图分类号：TJ55文献标志码：ADOI ：10.11943/CJEM20201131引言传统的单分子炸药（TNT 、RDX 、HMX 等）大多数由C 、H 、O 、N 元素组成，实测能量已基本达到其理论能量值，没有太大的提升空间。
为了寻求更高的能量释放效率，研究者们开始注意到单位质量和体积燃烧热都很高的活性金属燃料铝粉。
德国于一战末期时率先将铝粉加入到炸药中，以增加爆炸威力［1］，开创了铝粉在含能材料中的应用，在炸药中添加金属粉的方式改善了炸药的性能，满足了当时的使用需求。
尽管如此，随着对炸药能量密度更高的要求，研究人员发现不论如何减小铝粉的粒径［2］、提高活性铝含量［3］，甚至添加氧化剂［4］，混合炸药的放热速率仍然远低于单质炸药，实际放热量与理论计算值也相差甚远［5］。
与此同时，铝粉的加入导致含铝炸药长贮性存在问题［6］。
为了提高含铝炸药的能量利用率并改善其长贮性能，研究人员采用对铝粉进行包覆改性的方法改善其性能，取得了较好的使用效果［7-9］。
氟聚物因具有优异的耐腐蚀性、耐热性、低表面能以及很好的耐氧化特性而被广泛用于包覆铝粉。
使用氟聚物包覆铝粉，不仅可以有效阻止铝粉的进一步氧化，而且氟聚物的能量特性可提升体系的能量释放。
聚四氟乙烯复合材料的制备及其应用研究近年来，聚四氟乙烯（PTFE）复合材料在工业制造和科技领域中得到了广泛的应用，成为了新型材料领域的研究热点之一。
作为一种具有高强度、高稳定性、耐腐蚀性和生物惰性等优异性能的材料，PTFE不仅可以单独使用，而且还可以与其他材料复合加工，制成更加优质的复合材料，用于制造和生产多种产品。
一、PTFE复合材料的制备PTFE复合材料的制备方法多种多样，常见的有机械混合法、化学修饰法、物理吸附法、化学沉积法等。
在制备复合材料时，首先需要选取可与PTFE相容的材料，然后进行充分的混合和加工。
1.机械混合法机械混合法是指将PTFE和其他材料用机械方式进行混合。
这种方法的优点是简单易操作，生产成本低，但需要耗费大量的能量和时间。
机械混合法常用于制备各种PTFE/复合材料密封材料和弹性材料。
2.化学修饰法化学修饰法是指对PTFE表面进行化学修饰，使其表面具有亲和力能够与其他材料进行复合。
这种方法优点是可以制备出优异的化学和物理性能，一般适用于生产电子、化工和环保等行业的材料。
3.物理吸附法物理吸附法是利用PTFE表面的分子力或静电作用，将材料自然吸附于其表面。
这种方法优点是简单快捷，但存在着吸附量小、不牢固的问题。
物理吸附法常用于制备PTFE的表面性能改良剂。
4.化学沉积法化学沉积法是通过一个或多个反应进行PTFE/复合材料的制备。
这种方法的优点是制备速度快，材料齐全，但存在着制备条件严苛、成本高等问题。
化学沉积法的应用范围很广泛，可以用于制备高级传热材料、高性能材料、电子器件材料等。
二、PTFE复合材料的应用研究随着科技的不断发展，PTFE复合材料的应用领域也在不断拓展。
据统计，PTFE复合材料已广泛用于制造化学、电子、纺织、航空、船舶、汽车、建筑等工业领域。
以下是几种常见的PTFE复合材料的应用研究。
1.PTFE防爆电缆此电缆采取了PTFE与FEP的共混和配合制成的新型复合材料作为绝缘材料，具有高使用温度、强抗拉强度、抗化学腐蚀、不爆燃等特点，非常适用于石油、化工、冶金、纺织、军工等行业的防爆设备。
应材料ＰＴＦＥ/Ａｌꎬ反应破片沿着战斗部端面等距排列ꎬ破片总个数为２１个ꎬ密度为ρꎻ战斗部壳体材料为硬质铝合金ꎮ２㊀破片飞散特性破片的初始抛撒速度是破片飞散特性的主要内容之一ꎬ它是评估杀伤战斗部威力的重要依据[２]ꎮ发散角的大小直接影响到形成的破片束的侵彻能力和打击面积ꎮ本文的战斗部模型为完全轴对称的ꎬ选择破片速度方向与对称轴ｚ轴反向之间的夹角为发散角ꎮ根据破片的速度分析ꎬ可以确定每个破片的发散角α的值ꎬ进而近似求得破片的打击密度ꎮ本文只考虑破片发散角稳定时的状况ꎮ破片的发散角为:α＝ａｒｃｃｏｓｖｚｖæèçöø÷＝ａｒｃｔａｎｌＨæèçöø÷ꎻ(１)Ｈ＝ʃｔ００ｖｚｄｔꎻ(２)ｌ＝ʃｔ００ｖｘｄｔꎮ(３)式中:ｔ０为破片从开始飞散至击中目标的时间ꎻＨ为ｔ＝ｔ０时刻破片在ｚ方向的位移ꎻｌ为破片的发散半径ꎻｖｚ为ｚ轴速度ꎬ即轴向速度ꎻｖｘ为破片的径向速度ꎻｖ为破片初始抛撒速度ꎮ定义打击面积为Ｓꎬ破片总数为Ｎꎬ打击密度为ω(Ｓ)ꎬ打击密度作为打击面积的函数ꎬ即ω(Ｓ)＝ＮＳꎮ(４)其中:打击面积Ｓ＝π(Ｈｔａｎα＋ｌｚ)２ꎬｌｚ为某破片到弹轴线的距离ꎮ３㊀数值模拟３.１㊀计算模型、材料模型及参数㊀㊀采用ＡＮＳＹＳ/ＬＳ￣ＤＹＮＡ进行轴向驱动反应破片战斗部的数值模拟分析ꎬ采用流固耦合算法ꎮ数值模型由战斗部壳体㊁主炸药㊁内衬和反应破片(壳体和反应材料)组成ꎮ其中ꎬ主炸药㊁空气和内衬采用欧拉网格建模ꎬ单元使用多物质ＡＬＥ算法ꎬ战斗部壳体和反应破片均采用拉格朗日网格建模ꎮ采用ＴｒｕｅＧｒｉｄ软件进行有限元模型的建立ꎬ模型简化如图２所示ꎮ㊀㊀㊀㊀图２㊀战斗部有限元模型Ｆｉｇ.２㊀Ｆｉｎｉｔｅｅｌｅｍｅｎｔｍｏｄｅｌｏｆｔｈｅｗａｒｈｅａｄ㊀㊀网格单元采用Ｓｏｌｉｄ１６４八节点六面体单元ꎻ战斗部壳体材料模型选用Ｐｌａｓｔｉｃ＿Ｋｉｎｅｍａｔｉｃ材料模型ꎻ破片壳体和反应材料选用Ｍａｔ＿Ｊｏｈｎｓｏｎ＿Ｃｏｏｋ材料模型和Ｇｒｕｎｅｉｓｅｎ状态方程ꎻ主炸药采用Ｈｉｇｈ＿Ｅｘｐｌｏｓｉｖｅ＿Ｂｕｒｎ材料模型和ＪＷＬ状态方程描述ꎻ内衬模型选用Ｍａｔ＿Ｅｌａｓｔｉｃ＿Ｐｌａｓｔｉｃ＿Ｈｙｄｒｏꎻ空气采用Ｎｕｌｌ材料模型和Ｇｒｕｎｅｉｓｅｎ状态方程描述ꎮ主要材料参数如表１所示ꎬ参考值见文献[６ꎬ９￣１１]ꎮ３.２㊀数值模拟结果分析反应破片在主装药爆轰过程中ꎬ瞬间受到爆轰压力的驱动ꎬ以防止破片在加速飞散过程中不发生破碎或者受热提前反应ꎬ确保反应破片的完整性ꎬ避免毁伤能力削弱ꎮ本文中ꎬ研究了６种方案:无内衬和分别以泡沫铝㊁尼龙㊁泡沫塑料㊁橡胶和ＰＴＦＥ低密度材料作为内衬ꎬ其方案编号为１＃㊁２＃㊁３＃㊁４＃㊁５＃㊁６＃ꎬ数值模拟分析其对反应破片飞散特性和完整性的影响ꎮ３.２.１㊀反应破片飞散数值模拟结果㊀㊀通过数值模拟ꎬ得到了不同内衬时反应破片的初始抛撒速度和飞散角曲线ꎬ如图３和图４所示ꎮ由仿真计算可知ꎬ起爆后大约在２５μｓ时ꎬ破片的抛掷速度趋于稳定ꎬ爆轰压力与空气阻力相等ꎬ爆轰产物对破片束的作用基本消失ꎮ破片的最大速度出现在端面中心处ꎬ破片的最小速度出现在靠近壳体处ꎬ反应破片的最大速度为２４００ｍ/ｓꎬ最小速度为表１㊀主要材料参数Ｔａｂ.１㊀ＭａｉｎｐａｒａｍｅｔｅｒｓｏｆｍａｔｅｒｉａｌｓＢ炸药ρ/(ｇ ｃｍ－３)Ｄ/(ｍ ｓ－１)ｐＣＪ/ＧＰａＡ/ＧＰａＢ/ＧＰａＲ１Ｒ２ω１.７１７９８０２９.５５４２２７６.７４.２１.１０.３４ＰＴＦＥ/ＡｌＡ０/ＭＰａＡ１/ＭＰａＡ２/ＭＰａＢ０/ＭＰａＢ１/ＭＰａＢ２/ＭＰａ－４３７.９２３８.１７－３０.７２１８００.６５－１０９３.５２１７４.４１４５＃钢ρ/(ｇ ｃｍ－３)Ｇ/ＧＰａＡ/ＧＰａＢ/ＧＰａｎＣｍ７.８３７７７９.２５１０.２６０.０１４１.０３ ０５ ㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀爆㊀破㊀器㊀材㊀ＥｘｐｌｏｓｉｖｅＭａｔｅｒｉａｌｓ㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀第４５卷第４期㊀㊀图３㊀破片初始抛撒速度曲线Ｆｉｇ.３㊀Ｉｎｉｔｉａｌｄｉｓｐｅｒｓａｌｖｅｌｏｃｉｔｙｏｆｆｒａｇｍｅｎｔｓ㊀㊀图４㊀破片发散角曲线Ｆｉｇ.４㊀Ｅｊｅｃｔｉｏｎａｎｇｌｅｏｆｆｒａｇｍｅｎｔｓ１５９５ｍ/ｓꎮ由公式知:破片发散角主要集中在－１６ʎ~１６ʎ范围内ꎬ表明破片基本上沿着战斗部端面法线方向飞行ꎮ破片飞散稳定后ꎬ发散角标准差值较大ꎬ说明破片发散角度呈离散化ꎬ导致破片发散角增大ꎬ同时破片束的密度有所降低ꎮ３.２.２㊀内衬对反应破片飞散的影响分析通过数值模拟得到了不同内衬时反应破片的初始抛洒状态ꎮ由表２可以看出ꎬ内衬对破片的飞散特性影响较大ꎮ随着内衬材料密度的增加ꎬ破片的平均抛撒速度增大ꎬ平均发散角减小ꎬ并且破片的破损现象也相对减小ꎮ采用泡沫铝为内衬时ꎬ破片的速度最低ꎬ而采用密度相对较大的ＰＴＦＥ时ꎬ破片的速度最大ꎮ采用ＰＴＦＥ作为内衬时ꎬ其速度比无内衬时提高了２７.８％ꎮ从图５不同方案的破片变形状态也可以看出ꎬ随着内衬密度的增加ꎬ破片的完整性逐渐趋于完好ꎮ这是由于无内衬时ꎬ炸药与反应破片直接接触ꎬ当爆轰波传至端面时ꎬ爆轰产物迅速向周围运动ꎬ压缩空气形成冲击波ꎬ削弱了爆轰产物对破片的有效作功ꎻ同时ꎬ驱动反应破片沿着端面法线方向高速飞行ꎻ破片受到爆轰的直接载荷ꎬ造成外壳的变形ꎮ而采用内衬时ꎬ避免了反应破片与主装药的直接接触ꎬ一方面防止过度的动压加载ꎬ避免反应破片破碎ꎬ另一方面阻止过热的炸药爆轰产物的直接作用ꎬ削弱表２㊀不同内衬材料对反应破片飞散的影响Ｔａｂ.２㊀Ｅｆｆｅｃｔｓｏｆｄｉｆｆｅｒｅｎｔｌｉｎｅｒｍａｔｅｒｉａｌｓｏｎｄｉｓｐｅｒｓｉｏｎｏｆｒｅｓｐｏｎｓｅｆｒａｇｍｅｎｔｓ防护材料材料密度ρ/(ｇ ｃｍ－３)平均抛撒速度/(ｍ ｓ－１)平均发散角α/ʎ破片状态无内衬０１５００１１.８有破碎㊁变形泡沫铝１.０４０１６５０７.６有破碎㊁变形尼龙１.１４０１７８０７.３少量破碎㊁变形泡沫塑料１.２６５１８０４７.１少量破碎㊁变形橡胶１.４００１８５０７.０无破碎㊁变形ＰＴＦＥ２.１６０１９１７６.４无破碎㊁变形㊀㊀㊀㊀㊀㊀(ａ)方案１＃㊀㊀㊀㊀(ｂ)方案２＃㊀㊀㊀㊀㊀㊀(ｃ)方案３＃㊀㊀㊀㊀(ｄ)方案４＃㊀㊀㊀㊀㊀㊀(ｅ)方案５＃㊀㊀㊀㊀(ｆ)方案６＃图５㊀８０μｓ时不同方案破片的变形状态Ｆｉｇ.５㊀Ｄｅｆｏｒｍａｔｉｏｎｓｔａｔｅｏｆｆｒａｇｍｅｎｔｓｆｒｏｍｄｉｆｆｅｒｅｎｔｐｒｏｇｒａｍｓａｔ８０μｓ１５２０１６年８月㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀轴向驱动反应破片的作用过程㊀张世林ꎬ等㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀了爆轰波压力ꎬ避免反应破片过早受高热而发生快速分解反应ꎮ同时ꎬ内衬有效阻止稀疏波向爆轰产物内传入ꎬ起到了约束㊁延缓的作用ꎬ使得对破片的加载时间延长ꎬ提高了爆轰产物对破片的做功能力ꎮ４㊀试验及分析４.１㊀试验装置根据数值模拟分析ꎬ开展试验验证ꎮ如图６试验装置简图所示ꎬ战斗部距目标靶板距离３.０ｍꎬ靶板为１.２ｍˑ１.０ｍˑ３.０ｍｍ的钢板ꎮ靶心距离地面１.０ｍꎮ图７为反应破片战斗部ꎬ战斗部内置反应破片２１枚ꎬ破片密度为ρꎬ内衬层采用ＰＴＦＥꎬ共计４枚试验样机(图６所示)ꎮ试验时ꎬ战斗部横向放置在试验架ꎬ且中心与靶板中心水平对正ꎮ１－战斗部ꎻ２－测速靶ꎻ３－目标靶ꎻ４－高速摄像机ꎮ图６㊀试验装置布置示意图Ｆｉｇ.６㊀Ａｒｒａｎｇｅｍｅｎｔｏｆｔｈｅｔｅｓｔｄｅｖｉｃｅ㊀㊀㊀㊀图７㊀战斗部实物图Ｆｉｇ.７㊀Ｔｈｅｗａｒｈｅａｄ４.２㊀试验结果通过图８摄像可以看出ꎬ燃烧式反应破片在撞击靶板瞬间ꎬ伴有强烈的火焰产生ꎬ火焰长度约半米ꎮ这是由于高速飞行的破片在撞击靶板瞬间ꎬ破片内部压力瞬间变大ꎬ温度升高ꎬ反应材料受到高温高压后ꎬ达到临界爆压值ꎬ发生燃烧效应ꎮ由试验结果(图９)可知:靶板表面均出现穿孔痕迹ꎬ穿孔２１个ꎬ与预制反应破片数目相同ꎬ且弹坑处均有燃烧熏蚀痕迹ꎬ表明反应破片战斗部在炸药爆轰后ꎬ预制破片比较完整ꎮ在穿过靶板后ꎬ受到靶板的挤压摩擦后ꎬ反应材料发生爆炸燃烧反应ꎬ释放大量的热ꎬ大部分反应弹丸在靶板上能够有较完整的穿孔ꎬ穿孔孔径均大于１０ｍｍꎬ破片穿透靶板后仍有较高的剩余速度ꎬ且反应性黑色痕迹较多ꎬ反应破片分布在宽９４０~１２００ｍｍ㊁高８７０~１０００ｍｍ的区域内ꎬ最大飞散角为１８ʎ~２２ʎꎬ其打击密度为１７.５~２５.６ｍ－２ꎮ通过高速摄影图像ꎬ可近似求得４枚样机反应破片㊀㊀图８㊀破片撞击靶板瞬间燃烧图Ｆｉｇ.８㊀Ｂｕｒｎｉｎｇｒｅｓｕｌｔｓｏｆｆｒａｇｍｅｎｔａｔｔｈｅｍｏｍｅｎｔｏｆｉｍｐａｃｔｉｎｇｔｈｅｔａｒｇｅｔｐｌａｔｅ㊀㊀㊀㊀(ａ)１＃靶板㊀㊀(ｂ)２＃靶板㊀㊀(ｂ)３＃靶板图９㊀破片侵彻靶板毁伤效果图Ｆｉｇ.９㊀Ｄａｍａｇｅｒｅｓｕｌｔｓｏｆｆｒａｇｍｅｎｔｓｐｅｎｅｔｒａｔｉｎｇｔｈｅｔａｒｇｅｔｐｌａｔｅ２５ ㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀爆㊀破㊀器㊀材㊀ＥｘｐｌｏｓｉｖｅＭａｔｅｒｉａｌｓ㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀第４５卷第４期的平均速度为１８００~２０００ｍ/ｓꎮ㊀㊀试验表明:设计的轴向驱动反应破片战斗部发射的反应破片完整性较好ꎬ破碎情况较少ꎬ与数值模拟结果相吻合ꎮ５㊀结论１)数值模拟发现:无内衬时ꎬ轴向驱动反应破片抛撒速度最低ꎬ且在爆轰载荷下受损严重ꎻ采用密度相对较大的ＰＴＦＥ作为内衬ꎬ可以有效提高破片的完整性和抛撒速度ꎮ２)轴向驱动反应破片在侵彻靶板后ꎬ具有燃烧等二次毁伤效应ꎬ造成的复合毁伤效果优于传统型单一利用动能打击目标的惰性破片ꎮ３)轴向驱动反应破片战斗部具有定向功能ꎬ其速度提升和打击密度均优于传统型破片战斗部ꎮ参考文献[１]㊀周智超ꎬ吴晓锋ꎬ冷画屏.舰炮近炸引信预制破片弹在反导中的弹丸威力分析[Ｊ].军事运筹与系统工程ꎬ２００５ꎬ１９(２):６７￣７０.[２]㊀隋树元ꎬ王树山.终点效应学[Ｍ].北京:国防工业出版社ꎬ２０００ＳＵＩＳＹꎬＷＡＮＧＳＳ.Ｔｅｒｍｉｎａｌｅｆｆｅｃｔｓ[Ｍ].Ｂｅｉｊｉｎｇ:ＮａｔｉｏｎａｌＤｅｆｅｎｄＩｎｄｕｓｔｒｙＰｒｅｓｓꎬ２０００.[３]㊀邢恩峰ꎬ钱建平ꎬ赵国志.装药结构参数对轴向预制破片抛掷速度的影响[Ｊ].火炸药学报ꎬ２００７ꎬ３０(１):４９￣５３.ＸＩＮＧＥＦꎬＱＩＡＮＪＰꎬＺＨＡＯＧＺ.Ｅｆｆｅｃｔｏｆｓｔｒｕｃｔｕｒａｌｐａｒａｍｅｔｅｒｓｏｆｃｈａｒｇｅｏｎｔｈｅｄｒｉｖｉｎｇｖｅｌｏｃｉｔｙｏｆａｘｉａｌｐｒｅｆｏｒｍｅｄｆｒａｇｍｅｎｔｓ[Ｊ].ＣｈｉｎｅｓｅＪｏｕｒｎａｌｏｆＥｘｐｌｏｓｉｖｅｓ＆Ｐｒｏｐｅｌｌａｎｔｓꎬ２００７ꎬ３０(１):４９￣５３.[４]㊀张辉ꎬ李向豪ꎬ张选明.轴向增强炮弹破片场优化控制及反导效能分析[Ｊ].计算机仿真ꎬ２０１２ꎬ２９(９):２２￣２７.ＺＨＡＮＧＨꎬＬＩＸＨꎬＺＨＡＮＧＸＭ.Ｏｐｔｉｍａｌｃｏｎｔｒｏｌａｂｏｕｔａｘｉａｌｅｎｈａｎｃｅｄａｒｔｉｌｌｅｒｙｓｈｅｌｌｓｅｘｐｌｏｄｉｎｇａｎｄｄａｍａｇｅａｎａｌ￣ｙｓｉｓｏｆｍｉｓｓｉｌｅｓ[Ｊ].ＣｏｍｐｕｔｅｒＳｉｍｕｌａｔｉｏｎꎬ２０１２ꎬ２９(９):２２￣２７.[５]㊀黄亨建ꎬ黄辉ꎬ阳世清ꎬ等.毁伤增强型破片探索研究[Ｊ].含能材料ꎬ２００７ꎬ１５(６):５６６￣５６９.ＨＵＡＮＧＨＪꎬＨＵＡＮＧＨꎬＹＡＮＧＳＱꎬｅｔａｌ.Ｐｒｅｌｉｍｉｎａｒｙｒｅｓｅａｒｃｈｏｎｄａｍａｇｅｅｎｈａｎｃｅｄｆｒａｇｍｅｎｔ[Ｊ].ＣｈｉｎｅｓｅＪｏｕｒａｌｏｆＥｎｅｒｇｅｔｉｃＭａｔｅｒｉａｌｓꎬ２００７ꎬ１５(６):５６６￣５６９. [６]㊀徐松林ꎬ阳世清ꎬ张炜ꎬ等.ＰＴＦＥ/Ａｌ含能复合物的本构关系[Ｊ].爆炸与冲击ꎬ２０１０ꎬ３０(４):４３９￣４４４.ＸＵＳＬꎬＹＡＮＧＳＱꎬＺＨＡＮＧＷꎬｅｔａｌ.Ａｃｏｎｓｔｉｔｕｔｉｖｅｒｅｌａ￣ｔｉｏｎｆｏｒａｐｒｅｓｓｅｄＰＴＦＥ/Ａｌｅｎｅｒｇｅｔｉｃｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌ[Ｊ].ＥｘｐｌｏｓｉｏｎａｎｄＳｈｏｃｋＷａｖｅｓꎬ２０１０ꎬ３０(４):４３９￣４４４. [７]㊀阳世清ꎬ徐松林ꎬ张彤.ＰＴＦＥ/Ａｌ反应材料制备工艺及性能[Ｊ].国防科技大学学报ꎬ２００８ꎬ３０(６):３９￣４２ꎬ６２.ＹＡＮＧＳＱꎬＸＵＳＬꎬＺＨＡＮＧＴ.Ｐｒｅｐａｒａｔｉｏｎａｎｄｐｅｒｆｏｒ￣ｍａｎｃｅｏｆＰＴＦＥ/Ａｌｒｅａｃｔｉｖｅｍａｔｅｒｉａｌｓ[Ｊ].ＪｏｕｒｎａｌｏｆＮａｔｉｏｎａｌＵｎｉｖｅｒｓｉｔｙｏｆＤｅｆｅｎｓｅＴｅｃｈｎｏｌｏｇｙꎬ２００８ꎬ３０(６):３９￣４２ꎬ６２.[８]㊀张世林.轴向预制破片战斗部破片飞散特性影响因素分析[Ｄ].太原:中北大学ꎬ２０１２.ＺＨＡＮＧＳＬ.Ａｎａｌｙｚｅｏｎｔｈｅｆａｃｔｏｒｓｏｆｄｉｓｐｅｒｓｉｏｎｃｈａｒａｃ￣ｔｅｒｉｓｔｉｃｏｆａｘｉａｌｐｒｅｆａｂｒｉｃａｔｅｄｆｒａｇｍｅｎｔｓ[Ｄ].Ｔａｉｙｕａｎ:ＮｏｒｔｈＵｎｉｖｅｒｓｉｔｙｏｆＣｈｉｎａꎬ２０１２.[９]㊀孙文旭ꎬ李尚斌ꎬ黄亨建ꎬ等.防护材料对爆炸驱动反应破片的影响[Ｊ].科技导报ꎬ２０１１ꎬ２９(１６):３０￣３４.ＳＵＮＷＸꎬＬＩＳＢꎬＨＵＡＮＧＨＪꎬｅｔａｌ.Ｉｍｐａｃｔｏｆｔｈｅｐｒｏ￣ｔｅｃｔｉｖｅｍａｔｅｒｉａｌｓｏｎｅｘｐｌｏｓｉｏｎｄｒｉｖｉｎｇｒｅａｃｔｉｖｅｆｒａｇｍｅｎｔｓ[Ｊ].Ｓｃｉｅｎｃｅ＆ＴｅｃｈｎｏｌｏｇｙＲｅｖｉｅｗꎬ２０１１ꎬ２９(１６):３０￣３４.[１０]㊀才鸿年ꎬ王鲁ꎬ李树奎.战斗部材料研究进展[Ｊ].中国工程科学ꎬ２００２ꎬ４(１２):２１￣２７.ＣＡＩＨＮꎬＷＡＮＧＬꎬＬＩＳＫ.Ｒｅｓｅａｒｃｈｐｒｏｇｒｅｓｓｉｎｗａｒ￣ｈｅａｄｍａｔｅｒｉａｌｓ[Ｊ].ＥｎｇｉｎｅｅｒｉｎｇＳｃｉｅｎｃｅꎬ２００２ꎬ４(１２):２１￣２７.[１１]㊀时党勇ꎬ李裕春ꎬ张胜民.基于ＡＮＳＹＳ/ＬＳ￣ＤＹＮＡ８.１进行显示动力分析[Ｍ].北京:清华大学出版社ꎬ２００５.ＲｅｓｅａｒｃｈｏｎｔｈｅＭｅｃｈａｎｉｓｍｏｆＡｘｉａｌＤｒｉｖｉｎｇＲｅｓｐｏｎｓｅＦｒａｇｍｅｎｔｓＺＨＡＮＧＳｈｉｌｉｎ①ꎬＨＵＡＮＧＤｅｙｕ②ꎬＦＡＮＸｕｅｆｅｎｇ①ꎬＷＵＪｉａｎｐｉｎｇ①①７１０ＩｎｓｔｉｔｕｔｅꎬＣｈｉｎａＳｈｉｐｂｕｉｌｄｉｎｇＩｎｄｕｓｔｒｙＣｏｒｐｏｒａｔｉｏｎ(ＨｕＢｅｉＹｉｃｈａｎｇꎬ４４３０００)②Ｕｎｉｔ９５８７４ꎬＰＬＡ(ＪｉａｎｇｓｕＮａｎｊｉｎｇꎬ２１００２２)[ＡＢＳＴＲＡＣＴ]㊀Ｉｎｏｒｄｅｒｔｏｉｍｐｒｏｖｅｎａｖａｌｇｕｎ ｓａｂｉｌｉｔｙｏｆｔｅｒｍｉｎａｌｉｎｔｅｒｃｅｐｔｉｎｇａｎｔｉ￣ｓｈｉｐｍｉｓｓｉｌｅꎬａｎｄｉｍｐｒｏｖｅｔｈｅｄａｍａｇｅｅｆｆｉｃｉｅｎｃｙｔｏｔａｒｇｅｔꎬａｎａｘｉａｌｄｒｉｖｉｎｇｒｅｓｐｏｎｓｅｆｒａｇｍｅｎｔｗａｒｈｅａｄｗａｓｄｅｓｉｇｎｅｄ.ＡＮＳＹＳ/ＬＳ￣ＤＹＮＡꎬｔｈｅｎｏｎ￣ｌｉｎｅａｒｆｉｎｉｔｅｅｌｅｍｅｎｔａｎａｌｙｓｉｓｓｏｆｔｗａｒｅꎬｗａｓｕｓｅｄｔｏｐｅｒｆｏｒｍｔｈｅｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｔｈｅｒｅｌｉａｂｉｌｉｔｙｏｆｔｈｅｆｒａｇｍｅｎｔｗａｒｈｅａｄ.ＲｅｓｕｌｔｓｓｈｏｗｔｈａｔｕｓｉｎｇＰＴＦＥ(ｔｅｆｌｏｎ)ａｓｔｈｅｉｎｎｅｒｌｉｎｅｒｃａｎｅｆｆｅｃｔｉｖｅｌｙｉｍｐｒｏｖｅｔｈｅｉｎｔｅｇｒｉｔｙａｎｄｔｈｅｄｉｓｐｅｒｓａｌｓｐｅｅｄｏｆｆｒａｇｍｅｎｔｓ.Ｄａｍａｇｅｅｆｆｉｃｉｅｎｃｙｏｆｔｈｅｄｅｓｉｇｎｅｄｗａｒｈｅａｄｗａｓｖｅｒｉｆｉｅｄｂｙｅｘｐｅｒｉｍｅｎｔｓ.Ｔｅｓｔｒｅｓｕｌｔｓｓｈｏｗ:ａｆｔｅｒｐｅｎｅｔｒａｔｉｎｇｔｈｅｔａｒｇｅｔｐｌａｔｅꎬｔｈｅｂｕｒｎｉｎｇｆｒａｇｍｅｎｔｓｐｏｓｓｅｓｓｅｓｓｅｃｏｎｄａｒｙｄａｍａｇｅｅｆｆｅｃｔｓｄｕｅｔｏｄｅｆｌａｇｒａｔｉｏｎꎬｗｈｉｃｈｉｓｓｕｐｅｒｉｏｒｔｏｔｒａｄｉｔｉｏｎａｌｆｒａｇ￣ｍｅｎｔꎬｂｙｍｅａｎｏｆｔｈｅｋｉｎｅｔｉｃｅｎｅｒｇｙ.Ｉｔｃｏｕｌｄｇｒｅａｔｌｙｉｍｐｒｏｖｅｔｈｅｅｆｆｉｃｉｅｎｔｏｆａｎｔｉｍｉｓｓｉｌｅ.[ＫＥＹＷＯＲＤＳ]㊀ｒｅｓｐｏｎｓｅｆｒａｇｍｅｎｔꎻｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎꎻｄａｍａｇｅꎻｓｃａｔｔｅｒｉｎｇｃｈａｒａｃｔｅｒｉｓｔｉｃ ３５２０１６年８月㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀轴向驱动反应破片的作用过程㊀张世林ꎬ等㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀ｄｏｉ:１０.３９６９/ｊ.ｉｓｓｎ.１００１￣８３５２.２０１６.０４.０１２隧道掘进爆破振动在地表及上部岩体传播特征❋费鸿禄①㊀国志雨①㊀杨智广①㊀曲广建②㊀高㊀英③①辽宁工程技术大学爆破技术研究院(辽宁阜新ꎬ１２３０００)②广州中爆数字信息科技股份有限公司(广东广州ꎬ５１０６７０)③华南理工大学计算机科学与工程学院(广东广州ꎬ５１０６４１)[摘㊀要]㊀为了解隧道掘进爆破振动信号在地表及上部岩体内传播的特征ꎬ以鹤大高速小沟隧道掘进爆破为工程背景进行研究ꎮ对地表部分ꎬ现场监测后利用小波包分析对监测数据进行信号去噪重构处理ꎬ对得到的纯净信号进行时频分析及能量分布研究ꎬ进而得出振动信号沿地表的频率分布规律及能量衰减规律ꎮ对上部岩体ꎬ利用ＡＮＳＹＳ/ＬＳ￣ＤＹＮＡ进行数值模拟ꎬ提取关键点振速ꎬ分析得出开挖与未开挖两侧振动信号在隧道上部岩体内传播过程中的质点振速与能量衰减规律ꎮ为隧道爆破掘进动力响应研究㊁爆破震害评估和爆破设计的改进提供依据ꎮ[关键词]㊀振动信号ꎻ地表及上部岩体ꎻ小波包分析ꎻＬＳ￣ＤＹＮＡ数值模拟ꎻ传播特征[分类号]㊀ＴＤ２３５.１ꎻＯ３８２＋.２引言作为爆破信息的载体ꎬ振动信号的获取与分析始终是爆破工程的重要研究内容ꎬ尤其在当今数字爆破的大趋势下ꎬ振动信号数字化处理与研究对爆破震害的控制㊁爆破施工数字化的实现皆具有重要意义[１]ꎮ信号去噪是爆破振动信号分析的首要步骤ꎮ近年来ꎬ小波包变换取得了快速发展ꎮ小波包分解对低频部分㊁高频部分均可分解ꎬ克服了小波分解的局限性ꎬ比小波分解更为精细ꎮ研究证明ꎬ小波包及其衍生算法[２￣３]能够更准确地提取爆破地震波中不同频带下的振动分量ꎬ更清晰地研究各频带下爆破地震波的时频特征及能量分布特征[４￣５]ꎬ更清楚地揭示爆炸能量的传播及衰减规律[６￣７]ꎬ为现行的考虑质点振速和频率等多指标的爆破振动安全综合判据[８￣９]研究提供精确的数据信息ꎬ为深入研究隧道掘进爆破地震效应对周围环境的振动响应[１０]提供可靠的依据ꎮ利用现场监测与信号处理手段可准确分析振动信号沿地表的分布规律ꎬ但对于地表以下㊁掌子面以上的岩体内部空间ꎬ显然很难进行实地监测ꎬ而工程中往往遇到相邻隧道㊁管沟爆破施工的复杂情况ꎬ使研究振动信号在地表以下空间的传播规律变得十分必要ꎮ所以ꎬ笔者在分析振动信号沿地表衰减规律的基础上ꎬ利用ＡＮＳＹＳ/ＬＳ￣ＤＹＮＡ软件根据工程实际情况进行爆破掘进数值模拟ꎬ提取关键点振速ꎬ分析开挖与未开挖两侧振动信号在岩体内部的振速与能量衰减规律ꎮ１㊀现场监测１.１㊀工程概况鹤大高速小沟隧道[１１]位于吉林通化境内ꎮ在监测范围内ꎬ围岩主要为Ⅳ㊁Ⅲ级玄武岩ꎮ隧道掘进采用光面爆破技术ꎬ为避免爆破振动叠加导致振动过大[１２]ꎬ使用微差间隔不小于５０ｍｓ的导爆管雷管ꎮ根据现场条件ꎬ对处于Ⅲ级围岩区段的隧道掘进使用全断面一次爆破成型ꎻ对处于Ⅳ㊁Ⅴ级围岩使用台阶法施工ꎮ１.２㊀监测方案设计监测点集中于Ｋ５２３＋９６０~Ｋ５２５＋１１０段ꎮ该段围岩主要为强㊁中风化杏仁状玄武岩ꎬ围岩级别为Ⅳ级ꎬ岩体饱和抗压强度Ｒｃ＝６.２ＭＰａꎬ岩石完整系数Ｋｖ＝０.１１ꎬ岩体无断层ꎬ可基本排除地质因素对试验结果产生的影响ꎮ由于研究的是隧道掘进爆破振动在上部地表的响应ꎬ因此ꎬ主要对上台阶爆破引起的地面振动进行监测ꎮ本次试验共进行５组监测(１＃㊁２＃㊁３＃㊁４＃㊁５＃)ꎬ每组监测分别设置５个监测点(Ａ㊁Ｂ㊁Ｃ㊁Ｄ㊁Ｅ)ꎬ各组测点布设位置相同ꎮ如图１所示ꎬ其中Ｂ点位于掌子４５ ㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀爆㊀破㊀器㊀材㊀ＥｘｐｌｏｓｉｖｅＭａｔｅｒｉａｌｓ㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀第４５卷第４期❋收稿日期:２０１６￣０１￣２８作者简介:费鸿禄(１９６３－)ꎬ男ꎬ教授ꎬ博导ꎬ从事工程爆破㊁地下工程教学和科研工作ꎮＥ￣ｍａｉｌ:ｆｅｉｈｏｎｇｌｕ＠１６３.ｃｏｍ通信作者:国志雨(１９９２－)ꎬ男ꎬ硕士ꎬ主要从事地下工程方面的研究ꎮＥ￣ｍａｉｌ:９７０２２２８０３＠ｑｑ.ｃｏｍ㊀㊀(ａ)监测位置(ｂ)监测点的布置图１㊀监测位置及监测点的布置(单位:ｍ)Ｆｉｇ.１㊀Ｍｏｎｉｔｏｒｉｎｇｐｏｓｉｔｉｏｎｓａｎｄｐｏｉｎｔｓ(ｕｎｉｔ:ｍ)面正上方ꎬ相邻各点水平间距均为１０ｍꎮ１.３㊀监测结果监测共获取５组爆破振动数据ꎬ以第４组为例进行分析(表１)ꎮ该次爆破掘进共使用炸药１０８.０ｋｇꎬ最大一段起爆药量为２３.１ｋｇꎮ将Ｄ点ｙ方向振动信号简记为Ｄ(ｙ)ꎬ以Ｄ(ｙ)为研究对象进行小波包去噪重构处理及时频㊁能量谱分析ꎮ２㊀振动信号的小波包分析２.１㊀小波包去噪本次爆破振动监测仪器采用ＴＣ￣４８５０爆破测振仪ꎬ信号的记录时长为２ｓꎬ采集频率为２０００Ｈｚꎬ其奈奎斯特频率为１０００Ｈｚꎮ用小波包分解法对监测数据进行３层小波包分解ꎬ其二叉树如图２所示ꎮ小波包分解与重构算法[１３]不再赘余ꎮ图２㊀３层小波包分解树Ｆｉｇ.２㊀Ｔｈｒｅｅｌａｙｅｒｗａｖｅｌｅｔｐａｃｋｅｔｄｅｃｏｍｐｏｓｉｔｉｏｎｔｒｅｅ㊀㊀采用ｄｂ８作为基函数ꎬ应用Ｍａｔｌａｂ平台对信号进行小波包分析ꎮ将Ｄ(ｙ)信号进行小波包去噪ꎬ得到的重构信号与原信号进行对比ꎬ如图３所示ꎮ(ａ)Ｄ(ｙ)原信号(ｂ)去噪重构信号图３㊀Ｄ(ｙ)原信号及去噪重构信号Ｆｉｇ.３㊀ＯｒｉｇｉｎａｌｓｉｇｎａｌａｎｄｄｅｎｏｉｓｉｎｇｒｅｃｏｎｓｔｒｕｃｔｉｏｎｓｉｇｎａｌｏｆＤ(ｙ)㊀㊀对比可知ꎬ由于低频部分所含噪声信号很少ꎬ所以去噪后变化微小ꎻ大部分的高频噪声信号被过滤表１㊀质点振动监测结果Ｔａｂ.１㊀Ｍｏｎｉｔｏｒｉｎｇｒｅｓｕｌｔｓｏｆｐａｒｔｉｃｌｅｖｉｂｒａｔｉｏｎ监测点径向(ｘ)振动分量ｖｘ/(ｃｍ ｓ－１)切向(ｙ)振动分量ｖｙ/(ｃｍ ｓ－１)垂向(ｚ)振动分量ｖｚ/(ｃｍ ｓ－１)径向(ｘ)主频ｆｘ/Ｈｚ切向(ｙ)主频ｆｙ/Ｈｚ垂向(ｚ)主频ｆｚ/ＨｚＡ１.０２０.８９１.０７４３.５４７.０７１.４Ｂ１.２９１.３０１.６４８０.０８８.３５６.９Ｃ１.１４１.３０１.２８６７.６６６.７８３.３Ｄ０.６６０.３９０.５５５６.７５７.６６５.５Ｅ０.４００.４２０.４９５６.３５８.０６０.０５５ ２０１６年８月㊀㊀㊀㊀㊀㊀㊀㊀㊀隧道掘进爆破振动在地表及上部岩体传播特征㊀费鸿禄ꎬ等㊀㊀㊀㊀㊀㊀㊀㊀掉ꎬ使高频部分变得相对顺滑ꎮ２.２㊀振动信号时频分析为更加直观地反映振动信号去噪效果ꎬ利用Ｍａｔｌａｂ绘制振动信号Ｄ(ｙ)去噪前㊁后的三维小波包时频图ꎬ图４中ꎬＣ为主振频带能量百分比ꎮ(ａ)去噪前(ｂ)去噪后图４㊀Ｄ(ｙ)信号去噪声前㊁后的三维时频图Ｆｉｇ.４㊀Ｔｈｒｅｅｄｉｍｅｎｓｉｏｎａｌｔｉｍｅ￣ｆｒｅｑｕｅｎｃｙｍａｐｏｆｏｒｉｇｉｎａｌｓｉｇｎａｌａｎｄｄｅｎｏｉｓｉｎｇｓｉｇｎａｌｏｆＤ(ｙ)㊀㊀在去噪后的时频图中ꎬ中㊁高频部分以及０.５ｓ之后区域的信号能量大幅削减ꎬ干扰成分明显减少ꎬ爆破振动信号的能量分布更加集中于低频区域ꎬ可以认为得到了基本纯净的爆破振动信号ꎮ２.３㊀振动信号能量谱分析根据Ｐａｒｓｅｖａｌ定理[１４]ꎬ可得爆破振动信号的能量谱为:Ｐｉꎬｊ＝ʃｆｉꎬｊ２ｄｔ＝ðｎｋ＝１ｘｊꎬｋ２ꎮ(１)式中:ｘｊꎬｋ(ｊ＝０ꎬ１ꎬ２ꎬ ꎬ２ｉ－１ꎻｋ＝１ꎬ２ꎬ ꎬｎ)为ｆｉꎬｊ(ｔｊ)离散点的幅值ꎻｎ为爆破振动信号采集点数ꎮ爆破振动信号Ｓ(ｔ)的总能量Ｐ为:Ｐ＝ð２ｉ￣１ｊ＝０Ｐｉꎬｊꎮ(２)爆破振动信号Ｓ(ｔ)小波包分解到第ｉ层时ꎬ各频带能量占信号总能量的比例Ｃ为:Ｃ(ｎ)＝ＰｉꎬｊＰˑ１００％ꎮ(３)利用Ｍａｔｌａｂ编程求得并绘制振动信号Ｄ(ｙ)去噪前㊁后各频带能量分布规律ꎬ如图５所示ꎮ(ａ)原信号(ｂ)去噪后图５㊀Ｄ(ｙ)信号去噪声前㊁后的能量分布Ｆｉｇ.５㊀ＥｎｅｒｇｙｄｉｓｔｒｉｂｕｔｉｏｎｏｆｏｒｉｇｉｎａｌｓｉｇｎａｌａｎｄｄｅｎｏｉｓｉｎｇｓｉｇｎａｌｏｆＤ(ｙ)㊀㊀由图５可知ꎬ信号Ｄ(ｙ)去噪后ꎬ高频部分能量显著减少ꎬ基本为零ꎬ能量分布更加向低频区域集中ꎬ表明振动信号能量主要集中于２００Ｈｚ以内的低频部分ꎮ２.４㊀振动信号沿地表衰减规律分析对各测点振动分量进行小波包去噪重构ꎬ并进行能量谱分析ꎬ得到各振动分量的主频带㊁总能量及主频带能量比ꎬ如表２所示ꎮ表２㊀质点能量谱分析结果㊀Ｔａｂ.２㊀Ａｎａｌｙｓｉｓｒｅｓｕｌｔｓｏｆｐａｒｔｉｃｌｅｅｎｅｒｇｙｓｐｅｃｔｒｕｍ监测点主频带ｆ/(Ｈｚ)总能量Ｐ/(ｃｍ ｓ－１)２主振频带能量百分比Ｃ/％Ａ(ｘ)７.８~１０５.５９３.６９６６１.９６８Ａ(ｙ)０~１２１.１９２.４５１６９.２３０Ａ(ｚ)２３.４~１４４.５７０.２１７６０.３３７Ｂ(ｘ)７.８~１２５.０１１７.４９２６７.３９６Ｂ(ｙ)７.８~１２１.１１２０.７５６６３.７４０Ｂ(ｚ)０~１１７.２１１２.８９８６１.８８４Ｃ(ｘ)１５.６~１２１.１７４.９６９７０.９５３Ｃ(ｙ)０~５８.６７１.２６４６８.３７４Ｃ(ｚ)０~１２１.１６０.８１８７９.２６５Ｄ(ｘ)７.８~１２５.０４０.９５８７１.７３５Ｄ(ｙ)０~１２５.０４１.８２８８８.４１９Ｄ(ｚ)１５.６~１２１.１２４.１２１８１.９３７Ｅ(ｘ)１５.６~１２１.１３１.３７０８６.９３９Ｅ(ｙ)３.９~１２１.１２８.４３１９０.７５９Ｅ(ｚ)０~１２１.１３２.１８０８９.２０６６５ ㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀爆㊀破㊀器㊀材㊀ＥｘｐｌｏｓｉｖｅＭａｔｅｒｉａｌｓ㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀第４５卷第４期。
