Research on damage performance of steel tube reinforced by CFRP under three-point bending loads based on acoustic emission and bat algorithm
-
摘要:
研究了碳纤维增强树脂复合材料(Carbon fiber reinforced polymer，CFRP)加固Q345钢管在弯曲负荷下的损伤性能。通过三点弯曲试验，采用吸能特性分析方法评估不同加固方式下的抗弯强度和能量吸收性能。采用声发射(Acoustic emission，AE)技术，对比分析了不同CFRP铺层方式对钢管的加固效果及探究结构内部损伤和弯曲破坏的声学特征演化规律。最后提出了蝙蝠算法(Bat algorithm，BA)优化最小二乘支持向量机(Least squares support vector machine，LSSVM)的损伤分类预测模型。研究发现，增加CFRP缠绕层数可以显著提升钢管的抗弯强度和吸能能力，但增大缠绕角度会降低结构性能。通过对比分析不同加固方式下试件的声发射信号，证实了声发射技术在揭示碳纤维复合材料钢管弯曲过程中的损伤模式方面的有效性。能量概率密度的分析和最大似然评估显示，无论加固方式如何，复合管在不同能量级别上均遵循幂律分布，且能量分布指数随着CFRP缠绕层数增加而增大、随着缠绕角度增加而减小。所建立的BA-LSSVM损伤分类模型对试件损伤过程中的损伤程度分类准确性高达98%以上。
Abstract:In this study, the damage performance of carbon fiber reinforced polymer (CFRP) strengthened Q345 steel tubes under bending loads was researched. Through three-point bending tests, the bending strength and energy absorption performance under different reinforcement methods were evaluated using energy characteristic analysis. Additionally, using acoustic emission (AE) techniques, the reinforcement effects of different CFRP layup methods on steel tubes were comparatively analyzed, as well as exploring the evolution of acoustic characteristics of internal damage and bending failure. Finally, a damage classification model optimized by the bat algorithm (BA) for the least squares support vector machine (LSSVM) was proposed. The study finds that CFRP winding layers increasing can significantly enhance the bending strength and energy absorption capacity of the steel tubes, but increasing the winding angle will reduce the structural performance. By comparing the acoustic emission signals of specimens under different reinforcement methods, the effectiveness of acoustic emission technology in revealing the damage modes of carbon fiber reinforced steel tubes during bending was confirmed. Analysis of energy probability density and maximum likelihood estimation shows that, the composite tubes acoustic emission energy follow a power-law distribution at different energy levels, with the energy distribution exponent increasing with the increase of CFRP winding layers and decreasing with the increase of winding angle. The BA-LSSVM model was established to classify the degree of damage during the specimens damage process, with an accuracy of over 98%.
-
柔性压力传感器因其轻便性、柔韧性和生物兼容性等特点在航空航天、智能机器人、汽车、生物医学等领域应用广泛[1-4]。现有的柔性压力传感器按传感机制的不同可以分为压阻式、压电式、压容式和摩擦电式，其中压容式传感器因具有结构简单、稳定性强、响应迅速等优点而被广泛关注。灵敏度[5]作为压力传感单元性能中最为关键的指标之一，反映了传感器整体的力学、电学性能的优异程度。然而，制造具备出色灵敏度的柔性压力传感器也受到传感系统复杂性的严重制约，因此柔性压力传感器的灵敏度性能优化具有重要的研究价值。
触觉传感单元又称传感器触觉敏感单元，作为柔性压力传感器最为核心的组成部件，主要由聚合物纳米复合材料制备而成。聚合物纳米复合材料是以聚合物为基体连续相，以纳米级填充物为分散相的一种复合材料，其多样的基体、填料组合形式及不同的制备方式，赋予了聚合物纳米复合材料与众不同的力学[6]、电学[7]等特性，目前也广泛应用于应力/应变传感器[8]、超级电容[9]、生物医学[10]等工程领域。石墨烯和碳纳米管[11]等碳基材料因其优异的导电性和柔顺性，常作为聚合物填料应用于柔性压力传感单元的制备中。药芳萍等[12]通过设计双层砂纸与石墨烯的微组合结构制备了一种电容式柔性压力传感器，在0~25 kPa的压力载荷下灵敏度最高可达到0.451 kPa−1。Xu等[13]将碳纳米管与MXene结合到多孔微结构的聚二甲基硅氧烷海绵上研制出一种高灵敏度、宽传感范围、快响应速度的应变传感器。吴其皓等[14]以碳纳米管(CNT)与聚二甲基硅氧烷(PDMS)作为填充材料，并引入不同结构的砂纸作为模板对其电极结构进行优化，制备的轻薄型柔性传感器的灵敏度在0~1 kPa下可达到0.602 kPa−1。Ma等[15]将聚氨酯泡沫(PU)制备成“树枝状”骨架并用CNT和BaTiO3进行覆盖，在0~100 kPa的宽测量范围内灵敏度最高可达到2.51 kPa−1。Ke等[16]在热塑性聚氨酯(TPU)基体中以CNT和石墨烯片(MLG)为导电填料，并且通过优化填料比例使传感器在0~1.2 MPa压力下灵敏度最高可以达到2.05 MPa−1。Bi等[17]通过调节镍粉颗粒、碳纳米管、石墨粉混合比例的方式，优化了柔性压力传感器的灵敏度与蠕变性能，使其平均灵敏度在0~16 kPa的压力下达到
0.0589 kPa−1。王苏等[18]则采用商用聚氨酯海绵(PU)为基底，通过多次浸渍涂覆的方法将MXene/多壁碳纳米管(MWCNT)复合材料负载在海绵基底上，制备出一种柔性压力传感器，其在0~50 kPa的宽压力范围内灵敏度达到0.16 kPa−1。根据上述研究可知，设计具有不同微观结构的复合材料和掺杂具有高介电性能的填料颗粒是提高传感器灵敏度的重要手段，如引文[12-15]所述，但其材料结构设计过程复杂，耗时较长，且部分高介电颗粒的成本昂贵，故想要批量化生产此类传感单元十分困难。而基于共混方式制备传感单元的方法操作简单，耗时较少，且适合大批量生产，如引文[16-18]所述，但对于灵敏度性能的提升较微结构设计相对较低。针对以上问题，本文旨在采用操作简易的机械共混方法，以成本较低的碳纳米管(CNT)和石墨烯(MLG)作为导电填料，结合下述的人工智能优化算法，实现无微结构的柔性电容式传感单元的制备与性能优化。在确定所需的填料类型之后，由于填充颗粒之间的协同效应和反应机制的复杂性，这些填料的添加比例仍然难以确定，且在多种因素的交互影响下，不依靠数学模型对制造条件进行估计也是一件耗时漫长、成本高昂的工作。在材料合成和工艺优化领域，传统实验合成方法的分析和优化一直是重要的研究课题。最常用的统计分析模型响应面法(RSM)[19]，由Box和Wilson在1951年开发而来，其主要特点是通过少量试验创建预测模型来验证单个因素或多个因素的相互作用，用于评估和预测多参数实验条件，目前广泛应用于食品和化学领域[20]。但对于复杂的非线性问题，RSM二次多项式有时也无法达到所需的预测精度。支持向量机(SVR)作为机器学习中非线性回归模型，对非线性和非连通数据具有良好的预测能力。通过变换目标函数、改善方程的条件、降低计算的复杂度等操作[21]，可以提高其运算速度，提升预测精度与泛化能力，能在一定程度上克服RSM的缺点。虽然RSM和SVR模型已被用于各类工程领域的过程预测，但鲜有基于此类研究方法对传感器性能进行预测。蜣螂优化算法(DBO)[22]是一种新型群体智能优化算法，相较于传统的遗传算法[23] 和粒子群算法[24]，该算法的求解精度更高、收敛速度更快、全局搜索能力更强。目前该算法在三维路径规划[25]、电力信号预测[26]等问题上已有了初步应用。
综上所述，本文研究了电容式柔性压力传感单元导电机制，建立了相应的解析模型，确定了影响灵敏度性能的主要因素。采用实验设计法(DOE)中的中心组合设计法(CCD)[27]研究了CNT含量、MLG含量、搅拌时间和成型温度4个因素对灵敏度的影响。分别采用RSM和SVR机器学习算法对实验输入条件和输出性能建立预测模型，最后采用改进的DBO算法(IDBO)对该目标模型进行全局迭代寻优，确定传感器的最佳制备条件，最终实现无微结构的柔性电容式传感单元的性能优化。
1. 实验材料及方法
1.1 原材料及流程
本文所需原料包括硅橡胶基体(深圳红叶杰科技有限公司，双组分室温硫化硅橡胶RTV-2)、填充颗粒为多层石墨烯(MLG，苏州碳丰层数5~10层，厚度3.4~8 nm)、碳纳米管(CNT，苏州碳峰，纯度大于95wt%，长度大于3 μm，外径为8~15 nm)，二甲基硅油、偶联剂(KH550)与硅橡胶配套固化剂，皆购于深圳红叶杰科技有限公司。图1为压力传感器的制备流程。
![]()
图 1 触觉敏感单元制备流程图CNT—Carbon nanotubes; MLG—Graphene; KH550—Coupling agentFigure 1. Flow chart for preparation of pressure sensor首先，在室温下(25~30℃)，将一定比例的CNT和MLG加入研磨器研磨搅拌5~10 min，再加入到体积比为3∶1的乙醇水溶液中，超声分散15 min，把预处理后的填料颗粒在高温100℃左右进行干燥处理1 h。随后加入一定比例的偶联剂(KH550)进行表面改性处理，以增强分布宽度，减弱团聚效应[28]。并置于100℃高温干燥箱(DZF-6050B，上海合恒仪器设备有限公司)中进行烘干(1~3 h)，干燥后利用研钵器轻研，即得到改性后的导电填料颗粒。随后将导电填料与45 g RTV基体进行混合，同时加入质量比为5%~8%二甲基硅油进行稀释，并使用搅拌器在
1800 r/min的转速下搅拌一定的时间。搅拌完成后按混合比例1.5wt%加入硅胶固化剂(深圳红叶科技，主要成分为聚硅氧烷)再次搅拌1~2 min，并将搅拌混合后的聚合材料注入模具并放置于RS-1型真空箱(河南沃林仪器设备有限公司)中进行排气处理15 min，然后将模具放入设定的高温干燥箱中加热3 h，最后取出放置于室温下待自然冷却成形。图2为厚度4 mm、直径20 mm的圆形片状体，图3为触觉敏感单元在透射电镜(FEI Tecnai F20，日本电子) 100 nm放大倍数下改性前后的微观结构，从图3(a)可以明显看出，改性前CNT/MLG团聚现象较明显，且分布相对随机。而在改性后，如图3(b)所示，团聚明显减少，导电颗粒联系紧密且在基体中的分布相对更加均匀，保证了其稳定的传感性能。![]()
图 3 触觉敏感单元改性前(a)和改性后(b)对比Figure 3. Comparison between before (a) and after (b) modification of the tactile sensitivity unit1.2 试验方法
1.2.1 灵敏度测量及模型推理
灵敏度是衡量电容式传感器性能的重要指标，它代表了压力传感器电容量与载荷的相对变化关系[29]，式(1)为灵敏度的函数表达式。如图4所示，测量过程中使用微机控制电子万能试验机(FBS200N，弗布斯检测设备有限公司)与LCR数字电桥(UC2836B，优策电子科技有限公司)构成压力传感器电容测试回路，将样品放置在试验机中，使用试验机对样品施加0~30 kPa的单轴压力使样品产生塑性变形，同时通过LCR数字电桥和采集软件监测电容变化数据。最后得到不同压力情况下的压力传感器电容变化：
S=δ(ΔC/C0)δP=(εA/d−εA/d0)/(εA/d0)EΔd/d0 (1) 式中：S表示灵敏度(kPa−1)；C0为不施加压力时的初始电容值(pF)；ΔC=C−C0表示电容变化量，其中C为单轴压力条件下的电容值；P表示施加的压力(kPa)；E为弹性模量；d0是复合材料的初始厚度，Δd=d−d0是压缩过程中的材料厚度变化；复合压敏材料介电常数为ε；电极正对面积为A。
根据渗流理论[30-31]，在以CNT、MLG为填料的渗滤复合材料中，介电常数不再是固定的常数。渗滤复合材料作为介电层的柔性电容式压力传感器的工作原理如图5所示，在压缩过程中，填料浓度局部增加。介电常数随着浓度的增加而增加，当浓度接近渗流阈值时，介电常数增加速率变得更快，最终形成导电通路。
![]()
图 5 渗滤复合材料作为介电层的柔性电容式压力传感器的工作原理图d—Initial thickness; ε—Initial dielectric constant; P—Applied pressure; d–Δd—Material thickness after compression; ε+Δε—Dielectric constant of the composite pressure-sensitive material after compressionFigure 5. Working principle diagram of flexible capacitive pressure sensor with percolation composite as dielectric layer基于如下幂律方程所示：
ε=εd(fC−f)−t/(f<fC) (2) 其中：f和fC分别是导电填料的浓度和渗流阈值；ε和εd分别为复合材料和基体的介电常数；t表示比例常数。
故传感器的电容变化率可表示为下式，其中压敏单元厚度压缩比λ=d/d0：
ΔCC0=(εA/d−εA/d0)εA/d0=ελε0−1=[(fC−f0)/(fC−f)]tλ−1 (3) 综合上述表达式得到最终的传感器灵敏度解析模型如下：
S=[(fC−f0)/(fC−f)]t/λ−1E(λ−1) (4) 式中，f0表示初始状态的导电填料浓度。
由式(4)可知，电容式压力触觉敏感单元的灵敏度与聚合物纳米复合材料的介电性能和弹性模量密切相关。而合适的填料含量及填料在基体中的分布是否均匀是保证传感器优异介电性能的重要因素，同时根据Yang等[32]的研究可知，弹性模量E的屈服强度与温度密切相关。故本文将以填料含量、搅拌时间、成型温度作为主要影响因素，对触觉单元的灵敏度性能进行优化。
1.2.2 单因素实验设计
单因素实验设计是分别研究不同填料配比、搅拌时间、成型温度对传感器灵敏度性能的影响，为后续多因素实验设计提供一个合理的数值范围，如表1所示。试验时，首先固定MLG的含量为0.6 g，搅拌时间15 min，成型温度为 80℃，改变CNT的含量，分别是0.4、0.5、0.6、0.7、0.8、0.9、1.0、1.1 g；同理，固定 CNT的含量为0.7 g，保持搅拌时间与成型温度不变，改变MLG的含量，分别是0.4、0.5、0.6、0.7、0.8、0.9、1.0、1.1 g；而后保持CNT、MLG的含量与成型温度不变，分别为0.7 g、0.6 g、80℃，改变压力传感器制备过程中的搅拌时间分别为10、15、20、25、30 min；最后在保持填料含量与搅拌时间不变的情况下改变固化温度分别为40、60、80、100、120℃。实验重复 3 次，每个指标在相同条件下测定 3 次，最终得到同一条件下制备的传感器的平均灵敏度。
表 1 单因素实验方案Table 1. Single-factor protocolIndex Content of CNT/g Content of MLG/g Mixing time/min Molding temperature/℃ 1 0.4, 0.5, 0.6, 0.7, 0.8, 0.9,
1.0, 1.10.6 15 80 2 0.7 0.4, 0.5, 0.6, 0.7, 0.8, 0.9,
1.0, 1.115 80 3 0.7 0.6 10, 15, 20, 25, 30 80 4 0.7 0.6 15 40, 60, 80, 100, 120 1.3 预测与优化方法
1.3.1 响应面法(RSM)
响应面法(RSM)是一种用以评估一个或多个输出和输入参数之间关系建模方法，其计算量低，采用随机无偏设计，可以减少许多不必要的实验消耗，节约实验的时间与经济成本[33]。RSM中输入参数与输出变量之间的函数关系如下：
y=α0+N∑i=1αixi+N∑i=1αiix2i+∑i<j∑αijxixj (5) 其中：y是响应变量；xi是输入参数；α0表示偏差；αi表示线性效应；αii表示因素的平方效应；αij表示影响因素之间的交互效应。
RSM的实验设计法(DEO)主要分为两类：中心复合设计法(CCD)和Box-Behnken设计法(BBD)。为尽可能多的获取最佳实验条件，本文采用CCD法，进行多因素实验设计，其中每个输入参数有5个编码水平，分别为−2、−1、0、+1、+2。
图6(a)~6(d)分别为在得到每个单因素的实验结果后，根据其实验结果，以影响效果最大的点为中心点，确定了每个影响因素的实验范围，然后在实验范围内应用RSM确定多因素的实验方案，每个输入的具体水平因子显示在表2中，具体的CCD多因素实验方案如表3所示。
![]()
图 6 输入参数的单因素影响：(a) CNT含量；(b) MLG含量；(c) 搅拌时间；(d) 成型温度Figure 6. Single factor effect of input parameter: (a) CNT content; (b) MLG content; (c) Mixing time; (d) Molding temperature表 2 4个因素的水平和编码Table 2. Level and code of four factorsCode Levels factor CNT/g MLG/g Mixing
time/minMolding
temperature/℃−2 0.7 0.6 10 40 −1 0.8 0.7 15 60 0 0.9 0.8 20 80 1 1.0 0.9 25 100 2 1.1 1.0 30 120 表 3 多因素中心组合实验设计(CCD)与结果Table 3. Multifactorial central combinatorial experimental design (CCD) and resultsRun Factor Sensitivity/
(kPa−1)CNT MLG Mixing
timeMolding
temperature1 1 −1 −1 1 0.227 2 0 0 2 0 0.130 3* 0 0 0 0 0.118 4 −1 −1 −1 1 0.058 5 1 1 1 1 0.191 6 1 −1 −1 −1 0.208 7 −1 1 1 1 0.025 8* 0 0 0 0 0.145 9 −1 1 −1 1 0.050 10 −1 −1 1 −1 0.125 11 −2 0 0 0 0.050 12 0 −2 0 0 0.185 13 0 0 −2 0 0.193 14 −1 1 1 −1 0.164 15 1 1 −1 −1 0.375 16 1 −1 1 −1 0.180 17* 0 0 0 0 0.135 18 1 1 1 −1 0.310 19 0 1 0 0 0.198 20 2 0 0 0 0.478 21 −1 1 −1 −1 0.178 22 0 0 0 0 0.154 23 1 1 −1 1 0.326 24 0 0 0 −2 0.124 25* 0 0 0 0 0.117 26* 0 0 0 0 0.112 27 −1 −1 −1 −1 0.078 28 1 −1 1 1 0.248 29 −1 −1 1 1 0.090 30 0 0 0 2 0.073 Note: * indicates the central repeated trials. 1.3.2 支持向量回归(SVR)
支持向量机(SVM)是一种机器学习方法，如图7所示，它通过在高维空间中构建一个或一组超平面来实现分类或回归，其中用于回归的方法称为支持向量回归(SVR)。在SVR模型中核函数用以解决线性不可分的问题，其基本的数学表达式如下：
![]()
图 7 支持向量机回归(SVR)模型示意图Figure 7. Support vector machine regression (SVR) model diagramxn represents the independent variable; wij represents the weight coefficient; (βi, βi*) represents the Lagrange multiplier used for dual space transformation; K(x, xi) represents the kernel function used to solve the quadratic equation; b represents the deviation; y represents the dependent variabley=n∑i=1(βi−β∗i)K(x,xi)+b (6) 其中：(βi, β∗i)表示用于对偶空间变换的拉格朗日乘数；K(x,xi)表示用于求解二次方程的核函数；b表示偏差。
由于核函数是影响SVR模型精度的重要因素之一，因此为SVR模型选择合适的核函数尤为重要。目前，支持向量回归模型的核函数主要分为4种类型，即线性、多项式、径向基函数(RBF)和sigmoid。此外，核函数选取后，算法的超参数如g、C和ε也会影响SVR模型的性能，因此超参数的选取也是SVR建模的关键步骤。在本文中，RBF被用作模型计算的核函数，超参数优化则采用K折交叉验证的网格搜索方法(图8)。
1.3.3 蜣螂优化算法及其改进(IDBO)
蜣螂优化算法(DBO)作为一种新型启发式群智能优化算法，通过模拟蜣螂的滚球、跳舞、觅食、偷窃和繁殖等行为实现种群全局范围内最优值的选择。其适应度值作为衡量优化结果优劣的标准，优化的具体流程如图9所示。
![]()
图 9 蜣螂优化算法(DBO)流程图Figure 9. Dung beetle optimization (DBO) algorithm flowchartt—Number of iterations at the present stage; T—Number of iterations set by the algorithm同其他群智能优化算法类似，DBO通过随机生成位置的方式实现种群的初始化，但这容易出现种群分布不均匀，导致种群多样性减少，种群质量不高，从而影响算法的收敛速度及搜索性能。为了能够让个体在算法开始时有较高的全局搜索能力，需要让种群的位置均匀分布在整个问题的解空间内，混沌作为一种非线性的自然现象，因其混沌序列具有遍历性、随机性等优点，被广泛用于优化搜索问题。利用混沌变量搜索显然比无序随机搜索具有更大的优越性，因此本文采用改进的Tent混沌映射[34]替换原算法的初始化。原算法初始化与改进的Tent混沌初始化表达式如下：
xi,j=lbj+ri,j(ubj−lbj) (7) 式中：ri,j表示随机数列；lbj、ubj分别表示算法种群维度的上下界。
{2xi+rand(0,1)1N, 0⩽ (8) 式中， N 是序列内粒子的个数。引入随机变量 \mathrm{r}\mathrm{a}\mathrm{n}\mathrm{d}\left(0,1\right)/N 不仅仍然保持了Tent混沌映射的随机性、遍历性、规律性，而且能够有效避免迭代落入小周期点和不稳定周期点内。
在对种群初始化进行优化后，为了提升算法在迭代的不同时期全局/局部搜索能力。本文参考了文献[35]采用自适应 t 分布扰动策略，使DBO算法获得更优异的搜索性能。改进后的全局搜索位置更新公式为
\left\{\begin{array}{l}\mathrm{if\ rand} < p \\ x_i^t=x_i+x_it\left(i_{\mathrm{t}\mathrm{e}\mathrm{r}}\right)\end{array}\right. (9) 式中： p=0.5-0.1\left(T-t\right)/T 表示扰动概率； {x}_{{{i}}}^{t} 为突变后的蜣螂位置； {x}_{{{i}}} 为第 i 只蜣螂个体位置； t\left({i}_{\mathrm{t}\mathrm{e}\mathrm{r}}\right) 是以迭代次数 {i}_{\mathrm{t}\mathrm{e}\mathrm{r}} 为自由度的 t 分布函数；T为算法设定的迭代次数。
通过自适应 t 分布突变扰动策略虽然能增强算法全局搜索和跳出局部最优的能力，但是无法确定扰动之后得到的新位置一定比原位置的适应度值要好，因此在进行变异扰动更新后，加入贪婪策略[36]，通过比较新旧两个位置的适应度值，确定是否要更新位置。贪婪策略表达式如下：
{x}_{{{i}}}^{\mathrm{n}\mathrm{e}\mathrm{w}}\left(t\right)=\left\{\begin{array}{l}{x}_{{{i}}}^{}\left(t\right),\left(f\right({x}_{{{i}}}^{}\left(t\right)) < f({x}_{{\mathrm{g}}}\left(t\right)))\\ {x}_{{\mathrm{g}}}\left(t\right),\left(f\right({x}_{{{i}}}^{}\left(t\right)) > f({x}_{{\mathrm{g}}}\left(t\right)))\end{array}\right. (10) 其中： {x}_{{{i}}}^{\mathrm{n}\mathrm{e}\mathrm{w}}\left(t\right) 为经过贪婪规则更新后的个体位置； {x}_{{{i}}}\left(t\right) 为扰动后的个体； {x}_{{\mathrm{g}}}\left(t\right) 为扰动前的个体； f 代表蜣螂适应度函数。
IDBO算法的具体流程如图10所示。
![]()
图 10 改进蜣螂优化算法(IDBO)流程图P— Disturbance probabilityFigure 10. Improved dung beetle optimization (IDBO) algorithm flowchart为了验证所提蜣螂优化算法改进策略的有效性，本文将改进的IDBO与4种被广泛研究的优化算法(粒子群算法(PSO)、遗传算法(GA)、灰狼优化算法(GWO)、DBO)进行了比较。选取多种标准测试函数[37]评价其性能表现。其中，F1~F3为高维单峰测试函数，F4~F6为高维多峰测试函数。为保证实验公平性，所有算法的迭代次数与种群数都设置为100，结果如图11所示。
![]()
图 11 算法性能对比收敛曲线图GA—Genetic algorithm; PSO—Particle swarm optimization; GWO—Grey wolf optimizerFigure 11. Algorithm performance comparison convergence plot由图11可以看出，改进后的蜣螂算法在搜索精度及收敛速度上皆优于其他算法，其性能得到了显著提高。
1.3.4 预测模型精度的评价指标
在本文中，由于RSM、SVR对于传感单元灵敏度的性能预测本质上是一种回归问题，因此采用决定系数(R2)、均方根误差(Rmse)和平均绝对误差(Mae)用于评估二者所建立的预测模型的准确性，具体表达式如下：
{R}^{2}=1-\frac{\displaystyle\sum _{i=1}^{n}{\left({Y}_{{{i}}}-{{Y}_{\mathrm{p}\mathrm{r}\mathrm{e}\mathrm{d}}}_{{{i}}}\right)}^{2}}{\displaystyle\sum _{i=1}^{n}{\left({Y}_{{{i}}}-\overline{Y}\right)}^{2}} (11) {R}_{\mathrm{m}\mathrm{s}\mathrm{e}}=\sqrt{\frac{1}{N}\sum _{i=1}^{N}{\left({Y}_{{{i}}}-{{Y}_{\mathrm{p}\mathrm{r}\mathrm{e}\mathrm{d}}}_{{{i}}}\right)}^{2}} (12) {M}_{\mathrm{a}\mathrm{e}}=\frac{1}{N}\sum _{i=1}^{N}\left|{Y}_{{{i}}}-{{Y}_{\mathrm{p}\mathrm{r}\mathrm{e}\mathrm{d}}}_{{{i}}}\right| (13) 其中： {Y}_{{{i}}} 和 {{Y}_{\mathrm{p}\mathrm{r}\mathrm{e}\mathrm{d}}}_{{{i}}} 分别是通过实验测量的实际灵敏度和通过应用估算模型估算的灵敏度； \overline{Y} 是通过实验所测得灵敏度的平均值；N为本次实验测量的灵敏度样品的总个数。
2. 结果与讨论
2.1 响应面(RSM)模型的建立与分析
本实验采用Design-Expert软件对数据进行处理和分析，并通过方差分析(ANOVA)以检查RSM模型对数据的拟合情况，从而确定自变量对灵敏度的影响程度。表4为灵敏度方差分析结果。可以看出，该回归模型的F-value为32.91，P-value小于
0.0001 ，而失拟项的P-value为0.1422， 大于0.05，表明整体回归模型高度显著。RSM模型的R2为0.9685 ，证明了该模型的拟合程度较高。该回归模型的方差系数(CV)为14.86%，低于15%，CV值较低，表明其置信度较高，模型能很好地反映真实的实验数据。基于F-value的显著性排序为CNT>MLG>成型温度>搅拌时间。根据CCD的ANOVA回归系数，灵敏度与填料配比、搅拌时间、成型温度的响应关系可表示为表 4 灵敏度响应的方差分析Table 4. Analysis of variance for sensitivity responseSource Sum of squares df Mean square F-value P-value Status Model 0.2877 14 0.0205 32.91 < 0.0001 Significant {x}_{1} -CNT 0.1931 1 0.1931 309.36 < 0.0001 {x}_{2} -MLG 0.0117 1 0.0117 18.789 0.0006 {x}_{3} -Mixing time 0.0036 1 0.0036 5.73 0.0302 {x}_{4} -Temperature 0.0106 1 0.0106 17.02 0.0009 {x}_{1}{x}_{2} 0.0047 1 0.0047 7.46 0.0155 {x}_{1}{x}_{3} 0.0038 1 0.0038 6.11 0.0259 {x}_{1}{x}_{4} 0.0036 1 0.0036 5.81 0.0292 {x}_{2}{x}_{3} 0.0060 1 0.0060 9.68 0.0071 {x}_{2}{x}_{4} 0.0136 1 0.0136 21.83 0.0003 {x}_{3}{x}_{4} 0.0001 1 0.0001 0.2211 0.6449 {{x}_{1}}^{2} 0.0267 1 0.0267 42.79 < 0.0001 {{x}_{2}}^{2} 0.0085 1 0.0085 13.65 0.0022 {{x}_{3}}^{2} 0.0008 1 0.0008 1.26 0.2795 {{x}_{4}}^{2} 0.0030 1 0.0030 4.88 0.0432 Residual 0.0094 15 0.0006 Lack of fit 0.0079 10 0.0008 2.70 0.1422 Not significant Pure error 0.0015 5 0.0003 Cor total 0.2970 29 {R}^{2} 0.9685 {R}_{\mathrm{A}\mathrm{d}\mathrm{j}}^{2} 0.9390 {R}_{\mathrm{P}\mathrm{r}\mathrm{e}\mathrm{d}}^{2} 0.8311 CV/% 14.86 Adequate precision 23.8126 Notes: Cor total—Total number of correlation coefficients; {R}_{\mathrm{A}\mathrm{d}\mathrm{j}}^{2} —Adjusted coefficient of determination; {R}_{\mathrm{P}\mathrm{r}\mathrm{e}\mathrm{d}}^{2} —Predicted coefficient of determination; Adequate precision is an important statistical metric for evaluating the predictive ability of response surface models; CV—Coefficient of variance; df—Number of independent variables. \begin{split} Y=&0.1324+0.0897{x}_{1}+0.0250{x}_{2}-0.0122{x}_{3}-\\ &0.0210{x}_{4}+0.0171{x}_{1}{x}_{2}-0.0154{x}_{1}{x}_{3}+0.0151{x}_{1}{x}_{4}- \\ &0.0194{x}_{2}{x}_{3}-0.0292{x}_{2}{x}_{4}-0.0029{x}_{3}{x}_{4}+0.0309{{x}_{1}}^{2}+\\ &0.0228{{x}_{2}}^{2}+0.0053{{x}_{3}}^{2}-0.0104{{x}_{4}}^{2}\\[-1pt]\end{split} (14) 式中： Y 表示响应值，代表预测灵敏度； {x}_{1} 、 {x}_{2} 分别表示CNT和MLG的含量； {x}_{3} 代表搅拌时间； {x}_{4} 表示成型温度。
图12展示了各因素之间交互作用对灵敏度的影响，当研究两个因素之间交互时，其余两个因素处于中心水平状态。
![]()
图 12 灵敏度三维响应面图：(a) CNT/MLG相互作用；(b) CNT/搅拌时间相互作用；(c) CNT/成型温度相互作用；(d) MLG/搅拌时间相互作用；(e) MLG/成型温度相互作用；(f)搅拌/成型温度相互作用Figure 12. Three-dimensional response surface plots of sensitivity: (a) Interaction between CNT/MLG content; (b) Interaction between CNT/mixing time; (c) Interaction between CNT/molding temperature; (d) Interaction between MLG/mixing time; (e) Interaction between MLG/molding temperature; (f) Interaction between mixing time/molding temperature由图12(a)可以看出，随着CNT和MLG含量的提升，灵敏度的提升十分显著，这是由于这类碳基材料具有优异的导电性和柔顺性，两者结合能够实现基于二维石墨烯蜂窝结构与一维碳纳米管的三维导电网络，从而进一步增强其分散性[38]，增加填料与基体间的接触面积[39]，从而提升传感器的介电性能，这也与式(4)的推论相互印证。根据图12(b)、图12(d)可知，在CNT或MLG含量一定时，传感器的灵敏度性能随着搅拌时间的增加改变并不明显，而随着搅拌时间的继续增加，灵敏度有下降的趋势。这是由于长时间的混合剪切破坏了填料内部的三维导电网络结构[40]，这也表明适合的搅拌时间对压力传感器灵敏度的重要性。从图12(c)、图12(e)可知，在固化过程中随着温度的升高灵敏度逐渐下降。这是由于在高温固化的过程中基体持续软化使CNT和MLG之间的平均间距减小，这样便容易发生导电颗粒的团聚现象，从而对压力传感器的力学性能、电学性能产生了负面影响[41]。而从图12(f)可以了解到，在没有导电填料的情况下，仅对硅橡胶(RTV)基体进行搅拌和高温固化并不会对传感器的性能造成实质性的影响。
2.2 支持向量机回归(SVR)模型的建立与优化
基于上述CCD的数据样本，以CNT含量、MLG含量、搅拌时间和成型温度为输入参数，灵敏度为预测的输出变量建立预测模型。如图13所示，利用网格搜索法优化SVR的超参数，当C为
68.5935 ，g为0.7071 时，在5折交叉验证下的均方误差(MSE)达到最小值，平均R2值最大，即模型基于数据样本的预测性能达到最优。其中，基于SVR灵敏度预测模型的R2、Rmse和Mae值如表5所示。表 5 SVR模型的5折交叉验证预测性能评价指标Table 5. Performance evaluation index of SVR model by 5-fold cross-validationIndex {R}^{2} Rmse Mae 1 0.9893 0.0253 0.0193 2 0.9841 0.0328 0.0301 3 0.9735 0.0354 0.0231 4 0.9769 0.0389 0.0321 5 0.9882 0.0188 0.0151 Average 0.9824 0.0320 0.0239 Notes: Rmse—Root mean square error; Mae—Mean error rate. 从表5可以看出，在灵敏度预测性能评估中，平均R2、平均Rmse和平均Mae分别为
0.9824 、0.0320 和0.0239 ，这表明了预测值和实验值之间的高度相关性。此外，R2、Rmse和Mae的范围是0.9735 ~0.9893 、0.0188 ~0.0389 和0.0151 ~0.0321 ，这表明模型预测的波动范围极小，证明了此预测模型的高鲁棒性、高置信度。2.3 RSM与SVR预测模型对比
如表6所示，为比较两种模型的预测性能，基于式(10) RSM回归模型对相同数据进行了5折交叉验证。
表 6 响应面法(RSM)模型的5折交叉验证预测性能评价指标Table 6. Performance evaluation index of response surface methodology (RSM) model by 5-fold cross-validationIndex R2 Rmse Mae 1 0.9563 0.04002 0.03623 2 0.9642 0.02701 0.02443 3 0.9778 0.03103 0.02671 4 0.9669 0.04731 0.03113 5 0.9782 0.03312 0.02608 Average 0.96868 0.035698 0.028916 根据表5、表6的平均比较结果表明，SVR模型基于灵敏度性能的预测平均R2大于RSM模型，表明其适用性更好。SVR模型基于灵敏度预测的平均Rmse和Mae均小于RSM模型，更加表明超参优化后的SVR模型的预测精度优于RSM模型。
图14为SVR和RSM模型的5折交叉验证的雷达图。雷达图中可更加清楚地了解5折交叉验证下两预测模型的性能表现，图14(a)中SVR模型灵敏度的R2面积大于RSM模型，图14(b)、图14(c)中Rmse和Mae的面积均小于RSM模型，充分说明SVR模型在灵敏度的预测上具有更好的准确性。因此，通过此次比较分析表明，RSM和SVR模型在灵敏度性能预测方面都表现出良好的适应性和鲁棒性，而相较于RSM模型，SVR模型具有更好的预测效果。
![]()
图 14 SVR和RSM模型的5折交叉验证的雷达图Figure 14. Radar diagrams of SVR and RSM models with 5 fold cross validation2.4 基于改进蜣螂算法(IDBO)的模型优化
此优化过程是为了获得与灵敏度相适应的最佳材料配比及制作成型的工艺参数。通过上述的研究分析可知，由于通过超参优化后的SVR模型预测性能相优于RSM模型，故选择SVR模型用作优化的目标函数模型。通过IDBO算法不断迭代寻找最优的工艺参数，其中种群数量设为50，变量维度为4，迭代次数100次，迭代过程中的群体最优适应度值如图15所示。可知，随着迭代次数的增加，种群最优的适应度值也在不断上升，直至找到SVR模型的最优解，经过10次迭代后达到最佳适应度，即灵敏度性能达到最优，随后基本保持不变。最优的材料配比及制作工艺为CNT占比2.3wt%、MLG占比1.9wt%、搅拌时间15 min、成型温度78℃，仿真灵敏度为
0.5512 kPa−1。2.5 实验验证与对比
为了验证上述实验分析的合理性，根据上述优化结果的数据进行传感单元样品的制作，图16为所制作传感单元的灵敏度性能图，拟合直线斜率表示压力传感器在0~30 kPa时灵敏度大小。可知，采用优化后数据所制备的传感单元灵敏度S相较未优化时有大幅提升，达到
0.5371 kPa−1，与仿真灵敏度的相对误差为2.625%，也证明了基于SVR-IDBO优化预测模型的可靠性。![]()
图 16 灵敏度测试电容输出响应曲线S—Sensitivity; ΔC/C0—Relative change of the capacitanceFigure 16. Sensitivity test capacitor output response curves表7为本文所制备的柔性电容传感单元与近年来相关柔性压力传感器的灵敏度性能比较情况[12, 42-45]。结果表明，在相同压力载荷下，本文的传感器单元灵敏度性能处于较高的水平。其中，在同等压力下与张鹏等[45]的研究相比，灵敏度性能提升了5倍，与Yan等[44]相比提升了1.4倍，与药芳萍等[12]、赵珂[42]相比较，灵敏度性能分别提升了1.7倍与2.5倍。且与原始实验相比灵敏度最高提升了20倍。进一步表明结合SVR-IDBO方法进行触觉敏感单元性能优化的有效性与优越性。
表 7 传感单元灵敏度性能对比Table 7. Comparison of sensitivity performance of sensing unitsMethod Sensitivity/ \mathrm{kP}\mathrm{a}^{-1} Pressure range/ \mathrm{k}\mathrm{P}\mathrm{a} Comparison/ \mathrm{kP}\mathrm{a}^{-1} Ref. Microstructure design 0.451 0-25 0.7665 [12] Microstructure design 0.69 0-10 1.7531 [42] Microstructure design 1.01 0-25 0.7665 [43] Mechanical blending 0.3850 0-30 0.5371 [44] Microstructure design 0.341 0-5 1.9698 [45] Mechanical blending 0.5371 0-30 – This study 3. 相关性能测试
为满足传感器在实际场景的应用，本文的柔性触觉敏感单元除了在宽传感范围内具备较高的灵敏度之外，还具备其他较理想的特性，包括快速响应、高重复稳定性、良好的耐用性等。
如图17(a)所示，使用万能试验机对触觉敏感单元进行加载和卸载，在10 mm/min的较慢载荷进行速度下，其响应时间和恢复时间均能达到1 s，且每一循环峰值都极为接近，表明其在快速响应的同时也保证了电容信号的稳定传输。
如图17(b)、图17(c)所示，在15 kPa、30 kPa、45 kPa的载荷范围内进行10次循环测试，均能表现出良好的电容输出稳定性。且通过施加0~30 kPa的300次循环载荷也体现出触觉单元在受到载荷的中后期具有良好的重复稳定性。
4. 结 论
采用机器学习与智能优化算法相结合的方式改善了基于机械共混的柔性压力传感单元的制备条件。主要工作如下：
(1)通过推导压容式传感器灵敏度的机制，确定了影响其性能的主要因素有碳纳米管(CNT)含量、石墨烯(MLG)含量、搅拌时间和成型温度。采用中心复合实验设计法(CCD)动态检测了各种因素对输出性能的影响，结果表明，所有的因素均对输出性能有显著影响，并得到了响应面法(RSM)回归模型；
(2)基于CCD的实验数据作为输入，利用超参优化后的支持向量机(SVR)算法对灵敏度进行预测建模，并与RSM预测模型进行了比较，结果表明两种方法对灵敏度性能均表现出极高的预测精度，特别是超参优化后的SVR模型的预测结果与实际结果十分接近，这也表明所提出的方法在预测多因素影响和传感器性能之间的复杂关系方面具有巨大的应用潜力；
(3)利用改进的蜣螂优化算法(IDBO)对SVR模型进行迭代优化，逆向解析出了柔性传感器单元最佳的制备条件，实验结果表明，灵敏度性能大幅提升，且与仿真结果的相对误差较低。
综上所述，本文所采用的研究方法为传感单元性能优化提供了参考，同时为智能优化算法在柔性压力传感器方面的应用提供了理论依据，并且对实现柔性传感单元的大规模、工业化生产具有一定的指导意义。
-
![]()
图 1 碳纤维增强树脂复合材料(CFRP)-钢复合管的制备
θ—Winding angle
Figure 1. Preparation of carbon fiber reinforced polymer (CFRP)-steel composite tube
![]()
图 3 声发射(AE)传感器布置位置
Figure 3. Location of the acoustic emission (AE) sensor arrangement
![]()
图 6 CFRP加固钢管的荷载-位移曲线
Figure 6. Load-displacement curves of CFRP-reinforced steel tubes
![]()
图 7 CFRP加固钢管的吸能特性参数对比
Fmax—Peak load; Fmean—Mean load; CFE—Crushing force efficiency; SEA—Specific energy absorption
Figure 7. Comparison of energy absorption characteristic parameters of CFRP-reinforced steel tubes
![]()
图 8 CFRP加固钢管的声发射信号特征参数图
Figure 8. Acoustic emission signal characteristic parameter diagram of CFRP-reinforced steel tubes
![]()
图 9 CFRP加固钢管的能量释放特征
Figure 9. Energy release characteristics of CFRP reinforced steel tubes
![]()
图 10 CFRP加固钢管的声发射能量概率密度分布图
P/E—Energy probability density; Emin—Lower limit of energy normalisation; k1-k6—Critical indices of AE absolute energy distribution
Figure 10. Probability density distribution of AE energy of CFRP-reinforced steel tubes
![]()
图 11 CFRP加固钢管的声发射能量临界指数最大似然估计曲线
Figure 11. Maximum likelihood estimation curves of AE energy critical index of CFRP-reinforced steel tubes
![]()
图 12 蝙蝠算法(BA)流程图
rj—Transmitting frequency; Aj—Sound intensity
Figure 12. Bat algorithm (BA) flowchart
![]()
图 13 蝙蝠算法优化最小二乘支持向量机(BA-LSSVM)流程图
Figure 13. Bat algorithm for the least squares support vector machine (BA-LSSVM) flowchart
![]()
图 14 CFRP加固钢管的预测损伤类型与实际损伤类型对比图
Figure 14. Comparison chart of predicted damage types and actual damage types of CFRP-reinforced steel tubes
表 1 钢管及碳纤维性能参数
Table 1 Performance parameters of steel pipe and carbon fiber
Steel tube Carbon fiber Tensile strength/
MPaYield strength/
MPaExtension rate Gram weight/
(g·m−2)Tensile strength/
MPaElastic modulus/
GPaElongation rate 670 409 16 300 3870 2.45 1.74 
下载: 导出CSV
表 2 试件名称及参数
Table 2 Specimen name and parameters
Test specimen CFRP winding layer CFRP winding angle/(°) Quality/g ST – – 1651 C2T0 2 0 1799 C4T0 4 0 1953 C2T30 2 30 1798 C2T60 2 60 1801 C2T90 2 90 1800 Note: ST—Steel tube.
下载: 导出CSV
-
[1] 高小育, 段海, 杨正朴, 等. CFRP加固圆形钢管抗压承载力研究[J]. 建筑科学与工程学报, 2020, 37(6): 55-63. GAO Xiaoyu, DUAN Hai, YANG Zhengpu, et al. Research on the compressive capacity of circular steel pipe reinforced by CFRP[J]. Journal of Building Science and Engineering, 2020, 37(6): 55-63(in Chinese).
[2] 阳涛, 胡乔, 陈小兵, 等. 碳纤维布加固带缺陷的三层聚乙烯防腐涂层钢管及聚乙烯管道的试验研究[J]. 工业建筑, 2018, 48(12): 169-173, 180. YANG Tao, HU Qiao, CHEN Xiaobing, et al. Experimental study on carbon fiber cloth reinforced three-layer polyethylene anticorrosive coated steel pipe and polyethylene pipe with defects[J]. Industrial Architecture, 2018, 48(12): 169-173, 180(in Chinese).
[3] 代岩, 赵均海, 张常光. CFRP和角钢复合加固钢管混凝土叠合柱轴心受压承载力分析[J]. 建筑结构, 2018, 48(17): 96-103. DAI Yan, ZHAO Junhai, ZHANG Changguang. Analysis of axial compression capacity of CFRP and angle steel composite reinforced concrete filled steel tube composite column[J]. Architectural Structure, 2018, 48(17): 96-103(in Chinese).
[4] 黄辉, 路四方, 张祥, 等. 碳纤维复材加固钢管的抗内爆性能[J]. 工业建筑, 2022, 52(1): 211-215. HUANG Hui, LU Sifang, ZHANG Xiang, et al. Implosion resistance of carbon fiber reinforced steel pipe[J]. Industrial Building, 2022, 52(1): 211-215(in Chinese).
[5] 颜宇鸿, 卢亦焱, 李杉, 等. 碳纤维编织网增强ECC加固钢管混凝土短柱压弯承载力计算方法[J]. 建筑结构学报, 2023, 44(10): 178-187. YAN Yuhong, LU Yiyan, LI Shan, et al. Calculation method of compressive bending capacity of ECC reinforced concrete filled steel tube short column reinforced by carbon fiber braided mesh[J]. Journal of Building Structures, 2019, 44(10): 178-187(in Chinese).
[6] LI Z B, GAO Y, WANG Y, et al. Failure mechanisms and acoustic emission pattern recognition of all-CFRP cylindrical honeycomb sandwich shell under three-point bending[J]. Composites Science and Technology, 2023, 237: 1-13.
[7] LAVAYEN F D, BUTENEGRO G J A, BOADA M J S, et al. Theoretical and experimental study of the bending collapse of partially reinforced CFRP-steel square tubes[J]. Thin-Walled Structures, 2022, 177: 1-13.
[8] 黄俊杰, 佘艳华, 李猛, 等. 不同碳纤维布加固方式的木构件损伤演化规律[J]. 东北林业大学学报, 2024, 52(1): 115-123. HUANG Junjie, SHE Yanhua, LI Meng, et al. Damage evolution of wood parts with different reinforcement methods of carbon fiber sheet[J]. Journal of Northeast Forestry University, 2024, 52(1): 115-123(in Chinese).
[9] 邵家儒, 刘牛, 杨瑜, 等. CFRP复合材料构件胶接特性及失效规律研究[J]. 应用力学学报, 2023, 40(5): 1058-1067. SHAO Jiaru, LIU Niu, YANG Yu, et al. Study on bonding characteristics and failure rules of CFRP composite components[J]. Chinese Journal of Applied Mechanics, 2023, 40(5): 1058-1067(in Chinese).
[10] 黄书峰, 陈晓周, 刘东, 等. 碳纤维增强复合材料力学性能退化及失效过程的研究进展[J]. 宇航材料工艺, 2022, 52(5): 14-20. DOI: 10.12044/j.issn.1007-2330.2022.05.003 HUANG Shufeng, CHEN Xiaozhou, LIU Dong, et al. Research progress on mechanical degradation and failure processes of carbon fiber reinforced composites[J]. Aerospace Materials & Technology, 2022, 52(5): 14-20(in Chinese). DOI: 10.12044/j.issn.1007-2330.2022.05.003
[11] ALI Z A, TIAN X Y, LIU T F, et al. Mechanical and energy absorption behaviors of 3D printed continuous carbon/Kevlar hybrid thread reinforced PLA composites[J]. Composite Structures, 2023, 303: 1-12.
[12] WANG H, PEI Z, CONG W. A mechanistic cutting force model based on ductile and brittle fracture material removal modes for edge surface grinding of CFRP composites using rotary ultrasonic machining[J]. International Journal of Mechanical Sciences, 2020, 176: 105551.
[13] YU J G, ZHU S Q, FENG X T. Seismic behavior of CFRP-steel composite plate shear wall with edge reinforcement[J]. Journal of Constructional Steel Research, 2023, 203: 1-23.
[14] 陈卓异, 曾剑波, 彭岚, 等. 考虑缺陷位置和CFRP粘贴方式的钢板疲劳性能[J]. 中国公路学报, 2022, 35(2): 212-222. DOI: 10.3969/j.issn.1001-7372.2022.02.019 CHEN Zhuoyi, ZENG Jianbo, PENG Lan, et al. Fatigue properties of steel plate considering defect location and CFRP bonding[J]. China Journal of Highway and Transport, 2022, 35(2): 212-222(in Chinese). DOI: 10.3969/j.issn.1001-7372.2022.02.019
[15] 余海燕, 吴航宇. 碳纤维复合材料/钢的胶铆连接失效机理和选材方法[J]. 上海交通大学学报, 2023, 57(2): 230-240. YU Haiyan, WU Hangyu. Failure mechanism and material selection method of riveted joint of carbon fiber composite/steel[J]. Journal of Shanghai Jiao Tong University, 2023, 57(2): 230-240(in Chinese).
[16] 王庆松, 张玉, 张洪雨, 等. 基于改进Elman神经网络的CFRP补强钢板界面脱粘预测研究[J]. 振动与冲击, 2024, 43(3): 120-127. WANG Qingsong, ZHANG Yu, ZHANG Hongyu, et al. Prediction of interface desticking of CFRP reinforced steel plate based on improved Elman neural network[J]. Journal of Vibration and Shock, 2024, 43(3): 120-127(in Chinese).
[17] ZHANG Y, LIU Z Z, XIN J H, et al. The attenuation mechanism of CFRP repaired corroded marine pipelines based on experiments and FEM[J]. Thin-Walled Structures, 2021, 169: 1-11.
[18] 谢玉强, 佘艳华, 黄俊杰, 等. 碳纤维增强复合材料加固柏木轴压裂纹演化规律的研究[J]. 森林工程, 2024, 40(2): 102-116. DOI: 10.3969/j.issn.1006-8023.2024.02.012 XIE Yuqiang, SHE Yanhua, HUANG Junjie, et al. Study on axial compression crack evolution of cypress wood reinforced by carbon fiber reinforced composites[J]. Forest Engineering, 2024, 40(2): 102-116(in Chinese). DOI: 10.3969/j.issn.1006-8023.2024.02.012
[19] 田宝柱, 徐文涛, 梁鹏, 等. 基于多传感器融合的埋地输水管道泄漏声发射定位方法[J]. 科学技术与工程, 2023, 23(24): 10307-10316. DOI: 10.12404/j.issn.1671-1815.2023.23.24.10307 TIAN Baozhu, XU Wentao, LIANG Peng, et al. Acoustic emission location method for buried water pipeline leakage based on multi-sensor fusion[J]. Science Technology and Engineering, 2023, 23(24): 10307-10316(in Chinese). DOI: 10.12404/j.issn.1671-1815.2023.23.24.10307
[20] 王青原, 许颖, 钱胜. 基于机器视觉和数字图像相关技术的混凝土损伤演化研究[J]. 湖南大学学报(自然科学版), 2023, 50(11): 169-180. WANG Qingyuan, XU Ying, QIAN Sheng. Research on concrete damage evolution based on machine vision and digital image correlation technology[J]. Journal of Hunan University (Natural Science Edition), 2023, 50(11): 169-180(in Chinese).
[21] 许颖, 樊悦, 王青原, 等. 基于DIC的聚丙烯纤维增强混凝土断裂过程分析[J]. 华中科技大学学报(自然科学版), 2024, 52(2): 103-111. XU Ying, FAN Yue, WANG Qingyuan, et al. Fracture process analysis of polypropylene fiber reinforced concrete based on DIC[J]. Journal of Huazhong University of Science and Technology (Natural Science Edition), 2024, 52(2): 103-111(in Chinese).
[22] 袁逸齐, 兰恒星, 刘世杰, 等. 砂岩石窟热诱导裂纹损伤时空特征与分析[J]. 岩石力学与工程学报, 2022, 41(12): 2530-2542. YUAN Yiqi, LAN Hengxing, LIU Shijie, et al. Spatiotemporal characteristics and analysis of heat-induced crack damage in sandstone grottoes[J]. Chinese Journal of Rock Mechanics and Engineering, 2022, 41(12): 2530-2542(in Chinese).
[23] 许鑫浩, 徐福泉, 刘英利, 等. 红外热成像法检测预制混凝土构件外饰面内部缺陷试验研究[J]. 建筑科学, 2021, 37(7): 52-59. XU Xinhao, XU Fuquan, LIU Yingli, et al. Experimental study on detection of internal defects of precast concrete members by infrared thermal imaging[J]. Building Science, 2021, 37(7): 52-59(in Chinese).
[24] LIU Y H, ZHANG L, LI Z X, et al. Investigation on damage evolution of open-hole plain woven composites under tensile load by acoustic emission signal analysis[J]. Composite Structures, 2023, 305: 1-12.
[25] PEI N, XIANG Y X. Torsional damage analysis for pre-delaminated carbon/glass fiber-reinforced hybrid laminates based on acoustic emission[J]. Applied Acoustics, 2023, 202: 1-9.
[26] ZHAO W Z, PEI N, XU C G. Experimental study of carbon/glass fiber-reinforced hybrid laminate composites with torsional loads by using acoustic emission and micro-CT[J]. Composite Structures, 2022, 290: 1-9.
[27] ZHOU J, MATHEWS J V, ADAMS O D. Acoustic emission-based impact location estimation on composite structures[J]. Structural Health Monitoring, 2019, 18(5-6): 1652-1668.
[28] 林俊亭, 王帅. 基于DBN-MPA-LSSVM的无绝缘轨道电路故障诊断研究[J]. 电子测量与仪器学报, 2022, 36(9): 37-44. LIN Junting, WANG Shuai. Research on uninsulated track circuit fault diagnosis based on DBN-MPA-LSSVM[J]. Journal of Electronic Measurement and Instrumentation, 2022, 36(9): 37-44(in Chinese).
[29] 杨兴武, 王江, 孟致丞, 等. 基于电容电流状态估计的MMC多管开路故障诊断方法[J]. 中国电机工程学报, 2023, 43(23): 9297-9310. YANG Xingwu, WANG Jiang, MENG Zhicheng, et al. MMC multi-tube open-circuit fault diagnosis method based on capacitor current state estimation[J]. Proceedings of the CSEE, 2023, 43(23): 9297-9310(in Chinese).
[30] 李云淏, 咸日常, 张海强, 等. 基于改进灰狼算法与最小二乘支持向量机耦合的电力变压器故障诊断方法[J]. 电网技术, 2023, 47(4): 1470-1478. LI Yunhao, XIAN Richang, ZHANG Haiqiang, et al. Power transformer fault diagnosis method based on improved gray wolf algorithm coupled with least squares support vector machine[J]. Power Grid Technology, 2023, 47(4): 1470-1478(in Chinese).
[31] 杨海柱, 石剑, 江昭阳, 等. 基于CEEMD-SSA-LSSVM短期电力负荷预测模型[J]. 武汉大学学报(工学版), 2022, 55(6): 609-616. YANG Haizhu, SHI Jian, JIANG Zhaoyang, et al. Short-term power load forecasting model based on CEEMD-SSA-LSSVM[J]. Journal of Wuhan University (Engineering and Technology Edition), 2022, 55(6): 609-616(in Chinese).
[32] GUAN S Y, HUANG D, GUO S G, et al. An improved fault diagnosis approach using LSSVM for complex industrial systems[J]. Machines, 2022, 10(6): 443.
[33] MENG Q C, CHEN X M, ZHU Y L, et al. Communication signal classification and recognition method based on GA-LSSVM classifier[J]. Journal of Physics: Conference Series, 2019, 1345(2): 1-8.
[34] MA Q, DONG B, ZHA Y, et al. Multi-objective optimization for energy absorption of carbon fiber-reinforced plastic/aluminum hybrid circular tube under both transverse and axial loading[J]. Journal of Materials Engineering and Performance, 2020, 29(9): 1-16.
[35] 中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会. 金属材料弯曲试验方法: GB/T 232—2010[S]. 北京: 中国标准出版社, 2010. General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China, Standardization Administration of the People's Republic of China. Bending test method for metal materials: GB/T 232—2010[S]. Beijing: Standards Press of China, 2010(in Chinese).
[36] LIU Y S, CHEN X H, WU Z Y, et al. Effect of axial yarn distribution on the progressive damage behavior of braided composite tube subjected to three-point bending[J]. Thin-Walled Structures, 2022, 181: 1-12.
[37] 任松, 赵云峰, 张军伟, 等. 煤样巴西劈裂试验声发射能量幂律分布规律[J]. 东北大学学报(自然科学版), 2017, 38(4): 581-585. DOI: 10.3969/j.issn.1005-3026.2017.04.026 REN Song, ZHAO Yunfeng, ZHANG Junwei, et al. Power law distribution of acoustic emission energy in brazil splitting test of coal sample[J]. Journal of Northeastern University (Natural Science Edition), 2017, 38(4): 581-585(in Chinese). DOI: 10.3969/j.issn.1005-3026.2017.04.026
[38] CLAUSET A, SHALIZI C R, NEWMAN M E J. Power-law distributions in empirical data[J]. SIAM Review, 2009, 51(4): 661-703. DOI: 10.1137/070710111
[39] HOSSEIN N, POOYA R. Development of a probability distribution model for the SCFs in tubular X-connections retrofitted with FRP[J]. Structures, 2022, 36: 233-247.
-
目的
碳纤维增强树脂复合材料(Carbon fiber reinforced polymer, CFRP)因其高强度、轻质和优异的耐腐蚀性能而成为薄壁钢管结构的重要加固材料。然而，现有研究大多聚焦于观察CFRP加固钢管后的宏观力学性能和损伤表现，对加固钢管内部裂纹扩展规律与整体损伤的联系缺乏深入探究。为了研究不同CFRP加固方式对钢管内部损伤程度和裂纹发展规律的影响，揭示CFRP加固钢管的损伤机制，本文对6组不同CFRP铺层方式的钢管开展三点弯曲损伤监测试验。
方法试验以缠绕层数(0、2、4)和缠绕角度(0°、30°、60°、90°)为主要变量，设计了6组不同CFRP铺层方式的钢管。试验采用WAW-300B微机伺服控制液压万能试验机对CFRP加固木柱进行三点弯曲加载试验，其加载速率为2 mm/s，加载控制位移设置为60 mm。采用声发射技术(Acoustic emission, AE)实时监测CFRP加固钢管的三点弯曲损伤情况，通过采集的声发射信号，研究了CFRP加固钢管的力学性能、吸能性能和声发射参数特征，并提出了蝙蝠算法(Bat algorithm, BA)优化最小二乘支持向量机(Least squares support vector machine, LSSVM)的损伤分类预测模型。
结果从CFRP加固钢管的三点弯曲-声发射试验中得到了不同CFRP铺层方式下钢管的破坏形式、荷载-位移曲线和声发射特征参数。结果表明：①随着缠绕层数增加，复合管抗弯性能提高，峰值、平均载荷和比吸能均增强。与未加固钢管相比，2层、4层CFRP复合管峰值载荷分别提升12.19%、18.8%。②通过分析声发射信号能量的分布和强度，可以追踪弯曲破坏的累积过程，而通过分析CFRP加固钢管的声发射振铃计数可以评估破坏的程度。③随着CFRP缠绕角度的增加，各试件的AE信号值呈现下降趋势。这一变化表明CFRP的缠绕角度显著影响了AE的分布特征。在不同的缠绕角度下，C2T0试件显示出最高的损伤程度，而C2T90试件则表现出最低的损伤程度和最好的结构完整性。④CFRP加固钢管在三点弯曲载荷作用下，其损伤模式演变主要受CFRP缠绕层数和角度的影响。随着层数的增加，大尺度损伤如纤维断裂、基体剥离等占主导地位；随着缠绕角度的增加，小尺度损伤如界面剥离、基体微裂纹等所占比重逐渐增大。⑤利用直方图法和最大似然估计法得到的声发射能量概率密度临界指数相同，6种加固方式下，临界指数分别为1.31、1.45、1.55、1.69、1.60、1.42。⑥使用蝙蝠算法优化最小二乘支持向量机能够有效防止算法陷入局部最优解，显著加快了算法的收敛速度。建立的BA-LSSVM损伤分类模型对不同加固方式下钢管的损伤程度识别准确率达到98%以上。
结论①CFRP加固钢管在整个弯曲损伤过程中抗弯性能及吸能能力明显优于未加固试件。钢管的吸能性能随着CFRP缠绕层数增大而提升，而随着CFRP缠绕角度的增大而减少。缠绕角度从0°增加至90°时，峰值载荷和平均载荷均表现出下降趋势。②通过分析CFRP加固钢管的声发射振铃计数演化，可将其损伤破坏过程细分为弹性变形、渐进弯曲变形以及结构塌陷三个阶段。③声发射峰值频率的分布区间能够反映CFRP加固钢管的破坏尺度，随着碳纤维缠绕层数增大和角度的减小，CFRP加固钢管的损伤形式从大尺度损伤逐渐转变为小尺度损伤。④不同的缠绕层数和角度对应的复合管试件在声发射能量信号中遵循幂律分布，符合无尺度分布特征。由幂律方法确定的临界指数，能够反映不同缠绕方式下复合管的损伤性能及受弯破坏机理。⑤建立的BA-LSSVM损伤分类预测模型能够地对CFRP不同加固方式下钢管的损伤程度进行识别分类预测。
-
薄壁钢管在航天、公共交通和基础设施等行业中充当重要的缓冲与能量吸收角色。然而，长期的环境暴露和负载作用往往导致钢管结构的损伤和劣化。近年来，碳纤维增强树脂复合材料(CFRP)因其高强度、轻质和优异的耐腐蚀性能而成为薄壁钢管结构中的重要加固材料。然而，现有研究大多聚焦于观察CFRP加固后的宏观力学性能和损伤表现，对加固钢管内部裂纹扩展时声发射信号与整体损伤的联系缺乏深入探究。
本实验通过声发射技术，对六组CFRP加固钢管的弯曲损伤过程进行实时监测，研究了CFRP加固钢管在弯曲负荷下的损伤性能，揭示了CFRP增强钢管的破坏机制，试验结果可为CFRP加固钢管的损伤评估提供一定的依据。通过对比各组试件的压溃效率和比吸能发现，增加CFRP缠绕层数可以显著提升钢管的抗弯强度和吸能能力，但增大缠绕角度会降低结构性能。对声发射信号进行去噪处理提取到的峰值频率可以很好地表征钢管的损伤尺度。随着CFRP缠绕角度增加，低频和高频能量信号差别不大，但中频区的能量信号逐渐增多。这表明缠绕角度增加时，复合管试件的小尺度损伤加重，所占比重增大。不同缠绕层数和角度的复合管试件的声发射能量信号遵循幂律分布，符合无尺度分布特性，同时临界指数k也能反映材料的损伤程度。为了进一步对CFRP加固钢管在弯曲荷载作用下的损伤模式进行识别研究，采用机器学习对声发射数据进行分析，建立了蝙蝠算法优化最小二乘向量机模型，对钢管实际损伤进行分类识别。所建立的BA-LSSVM损伤分类模型对试件损伤过程中的损伤程度分类准确性高达98%以上。
CFRP增强钢管的压溃效率
CFRP增强钢管的比吸能
C2T0预测分类对比图
不同CFRP缠绕层数的声发射能量概率密度分布图(a)和最大似然估计曲线(b)
计量
- 文章访问数: 185
- HTML全文浏览量: 85
- PDF下载量: 23
