结构钢试件不同载荷下断铅声发射信号的混沌特征.pdf

科技与创新Science and Technology&Innovation1282023 年 第 17 期文章编号：2095-6835（2023）17-0128-04结构钢试件不同载荷下断铅声发射信号的混沌特征刘 婷，张 堃，郭晓艳（中山火炬职业技术学院，广东 中山 528437）摘要：Q235 钢作为用途最广泛的结构钢材料，对其内部损伤程度进行无损检测具有非常重要的意义。针对其结构内部受损伤后难以检测的问题，拟用结构钢试件承载不同载荷模拟内部损伤，采用断铅声发射信号检测金属材料内部损伤，并用关联维数、最大 Lyapunov 指数、Kolmogorov 熵 3 种混沌特征量分析其混沌特征。研究结果表明，金属材料内部损伤程度与 3 种混沌特征参数的变化规律具有一定的对应关系，断铅声发射信号的混沌特征可以表征金属材料内部损伤。这将为断铅声发射信号检测金属结构内部损伤提供一种新的方法。关键词：金属材料；结构钢；声发射；混沌中图分类号：TG142.15文献标志码：ADOI：10.15913/ki.kjycx.2023.17.038Q235 钢是一种碳素结构钢，由于它含碳适中，综合性能较好，因此无论是在建筑、工程结构中，还是在机械零件中，Q235 钢都有很大的需求体量。也正因Q235 钢应用的广泛性，在实际的工程应用中，对 Q235钢内部损伤程度进行无损检测，无论在事故预防，还是在旧零件再制造等方面都具有极其重要的意义1。现在常用的无损检测方法有很多，如超声检测、电磁监测、声发射检测、红外检测等，而在诸多无损检测方法中，声发射检测技术因其灵敏度高、抗干扰能力强等优点，日益成为无损检测方面的研究热点2。张颖等3对由Q235钢可视化原位拉伸实验所获取的声发射源信号进行了分析，根据声发射源动态观察结果及图像分析，建立了 Q235 的声发射细观损伤模型；陈冰等4通过时频特征分析在声发射信号中提取带锯条全生命周期的磨损信号特征，判断带锯条的初始剧烈磨损时间；毛汉颖等5通过断铅声发射信号，估算金属材料的声阻抗值，表征了金属材料的损伤程度；贺秀丽等6利用声发射技术对镁合金进行监测，分析了在腐蚀及疲劳载荷等不同情况下的声发射波形特征；綦磊等7以声发射技术为基础，采用混沌-BP 神经网络相结合的方法，实现了泄漏检测及对漏孔大小的评估。为了有效分析 Q235 钢金属材料内部损伤程度，基于金属材料损伤非线性的特点，本文首先对不同损伤情况的 Q235 钢进行断铅声发射试验，激励其内部损伤特征，获取断铅声发射信号；然后采用混沌分形理论对获得的断铅声发射信号进行特征分析，通过计算不同损伤程度下的结构钢的声发射信号的关联维数、最大 Lyapunov 指数及 Kolmogorov 熵值，进行综合比较，以期检测出声发射信号中的损伤信息，从而实现对金属材料内部损伤程度的判断。1断铅声发射试验设置试件采用尺寸为 210 mm30 mm1 mm 的 Q235钢材料，为了模拟损伤状态，对所有钢材进行拉伸试验，最大拉力分别设为 0 kN、4 kN、8 kN、12 kN、14 kN。受到 5 种不同级别的最大拉力的试件，分别用于模拟 5 种不同损伤程度的损伤试件。本试验平台采用美国物理声学公司（PAC）生产的 PCI-II 系统进行断铅试验：在试件的一端放置传感器，试件另一端设置断铅点，位置如图 1 所示。试验中，使用 HB 铅笔，保证其与试件表面成 30夹角，均匀用力压断铅芯，采集产生的声发射信号。单位：mm图 1试件示意图分别对先前制备的 5 种不同损伤程度的试件进行断铅试验，断铅产生的声发射信号沿试件传递，经过损伤部位后被传感器接收，记录声发射信号数据。为减少误差，对每组试件分别进行 5 次试验后，取获得的数据的平均值进行后续分析。基金项目 中山市社会公益与基础研究项目（编号：2021SYF08）；校级产学研专项项目（编号：2022CXY29）；校级产学研课题（编号：2022CXY33）；校级高职教育教学改革研究与实践项目（编号：2022JYQT015）断铅点传感器2106060Science and Technology&Innovation科技与创新2023 年 第 17 期1292混沌特征量算法2.1关联维数在机械故障诊断中，经常用关联维数8-10来定量说明系统的复杂程度，这是机械故障诊断中一种重要的方法。关联维数有多种计算方法，本文采用 G-P 算法进行计算：首先对原始时间序列进行相空间重构，取该相空间中的一个相点为中心，计算相点落入超球体之内的概率，即关联积分为：）（）（），（）（，）（0 112rdrHmNmrCmNjimij式中：m 为相空间重构中的嵌入维数的数值；N（m）为重构向量的个数；r 为超球体半径的数值；H 为Heaviside 函数。Heaviside 函数的公式为：0100 xxxH，）（将 r 适当缩小，根据关联积分与 r 的关系式，可得：drrmrC），（0lim式中：d 为关联维数的数值。关联维数 d 的公式为：rmrCdrlnlnlim0），（如果关联维数随着嵌入维数 m 的增加而不断增加，则表示数据为随机序列；如果关联维数随着嵌入维数 m 的增加而增加，最后无限趋于饱和，则表示为混沌序列。根据 G-P 法获得关联维数的计算结果，如图 2 所示。首先对 lnC（r）-lnr 曲线（如图 2（a）所示）中的直线段部分进行最小二乘拟合，即可求得每条直线段的斜率，结果如图 2（b）所示。关联维数 d 即斜率随 m 的增加而趋向饱和的那个值。从图 2（b）的曲线中可以看出，一开始，斜率随m 的增加而增加，当 m 增加到 14 之后，斜率变化变成在一定的区间内波动，并且逐渐平稳，趋于饱和。此时，对平稳区间内的斜率取平均值，即得到时间序列的关联维数为 1.28。同理，可以计算得到其他不同损伤程度的试件的信号的关联维数值，从而得到不同损伤程度下的试件的声发射信号的关联维数的变化趋势，如图 3 所示。-2.5-2-1.5-1-0.5-4-3.5-3-2.5-2-1.5-1-0.5lnrlnCr（a）lnC（r）-lnr 曲线（b）关联维数与嵌入维数关系图 2关联维数计算过程图 3不同损伤程度下的关联维数值2.2最大 Lyapunov 指数Lyapunov 指数是相邻相空间轨道按指数发散或收敛的速度的一种定量描述，本文采用小数据量法计算0481214210损伤程度/kN嵌入维数值lnr科技与创新Science and Technology&Innovation1302023 年 第 17 期各损伤程度下时间序列的最大 Lyapunov 指数值111-12。首先在重构的相空间中计算第 j 个参考点 Xj与最近相邻点jX之间的初始距离 dj（0），即：pjjXXdjjXjjmin0）（式中：p 为时间序列的平均周期的数值。计算出各参考点 Xj与最近相邻点jX的第 i 个离散时间步长后的距离 dj（i），即：）（）（）（），（，）（tijjijijjedidjMjMiXXid10min21 这里假设参考点Xj与最近相邻点jX具有1的指数发散率，根据上式计算可得到：qjjidtqiy1ln1）（）（利用最小二乘法拟合 y（i）曲线中的线性区域即可得到最大 Lyapunov 指数值1。如果最大 Lyapunov指数值大于 0，则判定系统混沌，并且最大 Lyapunov指数值越大，系统越混沌。以无损条件下试件为例，其最大 Lyapunov 指数值的计算过程如图 4 所示。根据小数据量法，用最小二乘法对 y（i）曲线的线性部分进行拟合，得到的斜率就是此时的最大 Lyapunov 指数值，为 0.012 4。同理，可以得到其他几种不同损伤程度下试件的声发射信号的最大 Lyapunov 指数值，其变化趋势如图 5 所示。05010015020025030011.21.41.61.822.22.42.62.8iy(i)图 4最大 Lyapunov 指数值的计算过程2.3Kolmogorov 熵Kolmogorov 熵（以下简称“K 熵”）是表达混沌系统的一个非常重要的量，它反映了系统的混乱度，其数值根据不同类型的动力学系统而不同：若 K 熵的值为 0，则系统为规则系统；若 K 熵的值是无界的，则系统为随机系统；若 K 熵的值大于 0 且有限，则系统为混沌运动系统，并且 K 熵值越大，信息损失速率越大，系统的混沌程度越高13。本文根据最大似然法计算各组不同损伤程度下试件的声发射信号的 K 熵，得到的结果如图 6 所示。图 5不同损伤程度下的最大 Lyapunov 指数值图 6不同损伤程度下的 K 熵值3断铅声发射信号的混沌特征分析由上述计算结果显示，采集到的不同损伤程度的试件的声发射信号时间序列均具有正的最大 Lyapunov指数，并且关联维数都能够达到饱和状态，表明采集得到的声发射信号时间序列具有混沌特征。为了便于分析，将上述计算得到的 3 种混沌特征量的值进行归一化处理，并绘制趋势图，如图 7 所示。从图中可以看到，随着试件损伤程度的增加，3种混沌特征值都呈上升趋势，不难分析，随着试件所损伤程度/kNi损伤程度/kNScience and Technology&Innovation科技与创新2023 年 第 17 期131受拉力的增加，试件的损伤程度不断增大，大量的微裂纹开始发生萌生，并且不断扩展，导致试件的非线性程度不断增加，所以混沌特征值也随之不断上升。图 7归一化特征值同时可以发现，在模拟试件损伤前期（所受拉力为 08 kN 时），关联维数的增幅相较其他 2 种混沌特征值更明显。分析可知，此时的试件处于拉伸变形的弹性变形阶段，微裂纹开始萌生，试件内部复杂程度高，而关联维数正好能够很好地描述系统的复杂程度，因此该值的增幅程度相对较大；在模拟试件损伤中后期（所受拉力为 812 kN 时），试件处于拉伸曲线由屈服阶段到塑性变形阶段的过程中，分布不均的微裂纹大量萌生和扩展，此时可以发现最大 Lyapunov指数的增幅更大，说明 Lyapunov 指数对金属材料内部损伤的发展阶段更加敏感；在模拟试件损伤后期（所受拉力为 1214 kN 时），试件到达塑性变形阶段，发生形变并趋于断裂，此时内部裂纹扩展大小方向趋于统一，混沌特征量也趋于饱和，其变化幅度不再像受损伤前中期那样剧烈，而是逐渐趋于平稳。4结束语为了实现对金属损伤程度的判定，本文通过利用Q235 钢进行声发射断铅试验，激励对其内部损伤进行混沌特征量分析，得到以下结果：试验结果表明，随着试件损伤程度的增大，混沌特征值也增大。因此，通过声发射信号计算得到的关联维数、最大 Lyapunov指数、K 熵的值与结构钢试件的损伤程度呈正相关关系，这 3 种混沌特征值反映了结构钢试件内部的复杂性和无序性，其变化趋势能够在一定程度上对金属材料损伤程度进行表征。不同损伤程度下试件的声发射信号所对应的混沌特征值敏感程度不同。不难发现，试验结果中计算得到的关联维数对结构钢试件内部发生损伤的萌生阶段更敏感，而在结构钢试件内部损伤的发展阶段，最大 Lyapunov 指数的增量更明显，因此，不同损伤阶段的结构钢利用不同的混沌特征值能够更好地实现损伤程度判定。通过试验证明，断铅声发射信号能够激励结构钢内部的损伤信息，并且结合结构钢声发射信号的多个混沌特征值，可以从不同角度来实现对结构钢损伤程度的共同判定，为金属损伤程度判断奠定了基础，并且为金属损伤识别的应用方面提供了一定的试验依据。参考文献：1李海洋，高翠翠，史慧杨，等.Q235 钢疲劳损伤的非线性Rayleigh 波检测技术研究J.中国测试，2018，44（4）：37-41.2郭景超，李安海.刀具磨损状态监测技术研究进展J.工具技术，2019，53（5）：3-13.3张颖，李越，周俊鹏，等.Q235 钢细观损伤声发射源动态观测实验研究J.实验力学，2020，35（2）：259-266.4陈冰，毕升豪，焦浩文，等.基于声发射的带锯条磨损在线监测试验研究J.工具技术，2021，55（4）：32-38.5毛汉颖，刘婷，范健文，等.金属材料声发射信号传播的声阻抗特性研究J.振动与冲击，2019，38（13）：197-201，214.6贺秀丽，梁红玉，闫志峰，等.腐蚀环境与疲劳载荷作用下 AZ31 镁合金声发射波形特征分析J.中国有色金属学报，2021，31（4）：826-835.7綦磊，孙立臣，闫荣鑫，等.基于混沌-BP 神经网络的真空泄漏声发射在线检测技术J.真空科学与技术学报，2017，37（6）：649-653.8关山，彭昶.金属切削过程刀具磨损信号的混沌特征J.振动与冲击，2015，34（10）：28-33.9羊自力，胡珊珊，王艳，等.基于小波分析的金刚石锯片锯切混凝土加工状态及表面质量分析J.工具技术，2018，52（4）：45-51.10潘伟，徐刚，熊方杰.基于非线性技术锅炉故障分析J.中国安全生产科学技术，2020，16（12）：176-182.11刘婷，毛汉领，黄振峰，等.金属材料声发射 Kaiser 效应的混沌特性分析J.振动与冲击，2017，36（12）：50-54.12杨永锋，仵敏娟，高喆，等.小数据量法计算最大 Lyapunov指数的参数选择J.振动.测试与诊断，2012，32（3）：371-374，511.13毛汉颖，刘畅，刘永坚，等.混沌特征量识别切削颤振的试验研究J.振动与冲击，2015，34（16）：99-103.作者简介：刘婷，广东中山人，专任教师。（编辑：丁琳）损伤程度/kN