基于压电振动能量俘获的弯曲结构损伤监测研究1)

2021-12-21李思谊李映辉

力学学报 2021年11期

赵翔李思谊李映辉

* (西南石油大学土木工程于与测绘学院工程安全评估与防护研究院，成都 610500)

† (西南交通大学力学与工程学院，成都 610500)

引言

近年来，压电材料因其具有能实现电能和机械能互相转换的特性，在工程结构中日益受到重视[1-2]，被广泛应用于结构的形状控制、振动和噪声控制、损伤检测等领域[3-4].从机械振动中提取能量的特性极大满足了传感器自给供电的经济需求与生态需求，因此压电材料广泛应用于俘能器中[5-6].到目前为止，最常用的压电俘能器采用悬臂梁式，通常是由一个结构层和压电材料(如压电陶瓷PZT、压电复合材料PVDF)层组成[7].学者们对直梁压电振动能量俘获系统进行了大量且详细的研究.Erturk和Inman[8-9]在Euler-Bernoulli 梁假设下，研究了在基座横向运动时悬臂式压电俘能器横向振动的精确解析解.另外考虑到常用的单自由度谐波激励对预测悬臂梁的运动可能会产生不准确的结果，推导了修正因子以改进单自由度横向和纵向振动激励模型.Zhao 等[10]应用格林函数法对悬臂式Timoshenko直梁压电俘能器进行了动力学分析，得到了强迫振动的解析解，探究了各类因素对电压响应的影响.Zhou 和Lei[11]研究了非对称三稳态压电俘能器的非线性动力学特性，以提高不同激励下的能量收集性能.Huang 等[12]对电阻-电感谐振电路非线性能量采集器的稳态响应机制进行了理论研究.Fang 等[13]根据哈密顿原理提出了利用离心软化效果分析旋转冲击能量俘获的理论模型.

而曲梁形式的压电俘能器的优点在于可用于多方向振动能量收集[14]，在微型电子系统以及复杂边界的工程结构中应用广泛[15].王剑等[16]构建了空间压电曲梁单元形状控制模型并得到了结构的最优形状控制电压.周勇等[17]提出了一种基于面内驱动的层合压电悬臂曲梁形式微位移驱动结构，忽略方程中曲率半径的影响获得了该曲梁在外荷载作用下控制方程的解析解.丁维高和谢进[18]利用哈密顿原理，研究了在水平与数值方向简谐振动激励的作用下压电曲梁俘能器的机电耦合的解析解.Kathpalia 等[19]建立了采用智能陶瓷压电材料的弯曲式压电传感器的分析模型，并获得了封闭式机电频率响应方程.Chen 等[20]采用数值方法分析了具有弹性约束的功能梯度压电曲梁的自由振动特性和瞬时响应.Tabatabaei-Nejhad 等[21]研究了嵌入压电层的层压曲梁的平面外振动特性，考虑了扭转和弯曲变形造成的剪切变形和旋转惯性效应.何燕丽和赵翔[22]利用格林函数法求解Prescott 型曲梁俘能器强迫振动的响应.

曲梁结构如今广泛地应用于各种领域的工程结构中，在外荷载激励下难免造成损伤产生裂纹，从而严重影响材料结构的整体性能和使用寿命[23-24].对损伤曲梁的振动特性研究也一直是国内外研究的热点[24-27].为识别结构中的损伤发展了许多检测方法，分别有超声波法、电涡流法、磁粉法和红外识别法[28]，以及基于结构振动特性的损伤识别法，后者具有明显的非破坏性，方便、快速、便宜的优点被广泛采用[28].另外，还有基于压电效应的主动感应法对结构的健康监测和损伤检测方法.孙威等[29]利用智能压电传感器对混凝土构件的裂缝损伤发展状况进行监测，探索一种基于压电波动法的裂缝损伤监测方法.Hou 等[30]将压电智能材料和压电陶瓷材料应用于传感器，对地震力作用下的混凝土结构整体剪切应力进行监测研究.Markovic 等[31]建立了使用压电智能传感器的主动监测系统和损坏检测的有效有限元模型，对基于压电波动法的埋入式压电陶瓷混凝土梁损伤识别系统进行模拟分析.由此，本文提出可以采用逆方法，即通过振动引起的压电俘能器的信号响应变化分析曲梁的振动问题，从而对结构进行损伤监测.

目前已经有不少学者对裂纹曲梁的振动问题和曲梁式压电俘能器进行了研究，但是有关损伤曲梁式压电俘能器强迫振动问题的研究较少.在上文已经阐述了裂纹存在给结构带来的负面影响，因此对含裂纹损伤的曲梁压电俘能器的研究是有实际意义的.本文不仅研究其强迫振动问题，还采用逆方法，即利用损伤曲梁振动时的俘能器电压响应变化，提出了一种监测结构损伤的方法.

本文建立了含裂纹损伤的曲梁压电俘能器的力电耦合模型，得到其强迫振动的解析解并基于压电俘能器的响应监测曲梁结构的损伤情况.基于文献[22]中Prescott 型曲梁压电俘能器的解析解，引入局部刚度模型描述裂纹损伤截面的力学特性，应用格林数法和传递矩阵法推导含多裂纹曲梁压电俘能器强迫振动下的解析解.根据线性系统叠加原理，对力电耦合系统解耦，推导了该裂纹损伤模型的输出电压的解析式.数值模拟了压电俘能器在曲梁结构处于振动状态下的响应信号，提出了对结构进行损伤监测的方法，以期为曲梁式压电俘能器的振动研究以及损伤结构的监测和优化设计提供理论参考.

1 损伤曲梁压电俘能器的力电耦合模型

建立如图1所示的含裂纹损伤的悬臂式曲梁压电俘能器模型.俘能器中的压电层材料比结构层材料的柔性好，因此结构层上更容易产生裂纹.如图1(a)所示，该模型结构层上存在n个开口裂纹，假设受到外力作用时，裂纹截面处的局部刚度发生变化，曲梁的挠度出现不连续性，曲梁的未损伤部分没有受到影响[32].因此，在裂纹损伤截面处建立局部柔度模型用来描述截面的力学特性，将裂纹处视为一个刚度为Keq的等效无质量扭转弹簧，该弹簧连接着相邻的两个完整部分，可表示为

图1 损伤曲梁压电俘能器模型Fig.1 The model of curved piezoelectric energy harvester with n cracks

式中EI为压电曲梁的有效刚度，hs是结构层曲梁的高度，局部柔度的无量纲常数，在本模型中采用

式中h′=hc/hs是裂纹深度与曲梁高度的比值.

如图1(b)所示，压电曲梁被n个裂纹分成n+1段完整部分，每段长度为Li(i=1，2，···，n+1) 并且由一个等效无质量扭转弹簧连接.为了得到裂纹压电曲梁的格林函数，需要分别研究每段的变形情况.因此，建立n+1个局部坐标系oixiyi(i=1，2，···，n+1)[33]表示出出每段的格林函数后，转换局部坐标为全局坐标即即可得到压电曲梁的格林函数.

已知文献[22]中建立的Prescott 曲梁模型，没有考虑曲梁轴向力以及轴向惯性项，仅考虑了曲梁的径向位移.因此曲梁压电俘能器的振动控制方程是一个四阶偏微分方程，引用文献[22]中的格林函数，该格林函数是强迫振动下的压电曲梁径向位移W(x)的解

式中x0表示单位力作用的位置，，，是径向位移Wi(i=1，2，3)的各阶导数，表示压电俘能器的电压.G1(x，x0)表示由外力引起的位移解，G2(x，x0)，G3(x，x0) 是由电耦合效应引起的位移解，且有

式(4)和式(5)中

且si(i=1，2，3，4) 是下列代数方程的根

式(6)中，R是曲梁的半径，μ 和 γ 表示单位长度曲梁的质量和转动惯量，c1表示阻尼系数，A表示截面面积. ϑ1和 ϑ2是力电耦合系数，具体表达式见文献[22].

在本文的后续推导中需要用到轴向位移，根据曲梁控制方程[22]可推出轴向位移V(x)

根据式(3)～式(9)得到由外荷载引起的轴向位移的格林函数

在第i(i=1，2，···，n-1) 个裂纹位置处和在第n个裂纹位置处各状态矢量的传递关系[34]如式(12)和式(13)所示，曲梁的轴向位移、弯矩相等，剪力大小相等但在第n处方向相反，截面转角始终不连续

式中，表示截面转角;Ji和Jn表示第i处和第n处的等效扭转弹簧的局部柔度，是刚度Keq的倒数.

将局部格林函数式(11) 带入传递关系式(12)和式(13)中可得到两个矩阵形式的传递方程

根据传递方程式(14a)和式(14b)建立未知常数向量U1和Un+1的关系

考虑到本文采用的是悬臂式压电俘能器，已知悬臂梁的边界条件是

将边界条件(19)带入传递式(18)即可求解未知常数向量U1和Un+1，再通过式(14a)和式(14b)可以定义其余未知常数向量.

因此，可以得到所有局部坐标系下的格林函数，转换坐标后，含n个裂纹的曲梁压电俘能器的格林函数为

2 损伤力电耦合系统解耦

根据线性系统和格林函数的叠加原理，建立裂纹曲梁的力学位移与电力学方程之间的解析关系，含裂纹曲梁压电俘能器的位移W(x) 为

假设在基础加速度作用下压电曲梁上的外力f(x，t)是简谐力，可将外力[35]写作

式中，A0是指位移的振幅.

结合文献[22]中已知的输出电压和本文求解的多裂纹曲梁压电俘能器的格林函数解析式，可以得到本模型的输出电压

因此，压电俘能器的功率为

3 数值结果与讨论

本文探究了弧长L=0.1 m，宽b=0.01 m 受简谐激励的悬臂曲梁压电俘能器，曲梁由上层的压电层和下层的结构层组合，其厚度分别为0.005 m 和0.005 m.曲梁几何和材料参数、压电参数取值如表1 所示.其次，从裂纹几何参数(裂纹深度和位置)对裂纹压电曲梁进行了模拟，裂纹深度比h′分别为0.02，0.04 和0.06，用裂纹位置L1表示离固定端的距离.在以下算例中，没有考虑阻尼的影响.为了方便起见，引入阻尼效应的无量纲化

表1 裂纹压电曲梁的几何参数和压电参数取值Table 1 Geometrical and electromechanical parameters of the cracked beam

3.1 解的有效性验证

令裂纹深度为零，含裂纹的曲梁可退化为完整曲梁模型.利用这一点将含裂纹压电曲梁的解析解与文献[22]得到的完整压电曲梁模型的解作对比，验证结果是否一致.如图2 所示，本文的位移频率响应结果与文献中的结果基本吻合，从而验证了本文解的有效性.

图2 位移-频率响应的对比图Fig.2 The displacement for comparison of present solutions and the solutions from reference

3.2 裂纹损伤对压电俘能器电压-频率响应的影响

本节探究了压电俘能器电压响应在裂纹几何参数影响下的变化.一方面为相关设计提供了理论依据，另一方面介绍了如何通过电压响应的变化对结构进行损伤监测.

图3 绘制了不同裂纹深度时压电俘能器的电压-频率响应图.从图中可以看出，频率在到5 kHz的时候，各裂纹深度的电压响应之间有重叠部分.当振动频率大于此重叠频率振动时，损伤压电曲梁的响应效果与完整曲梁有了显著差异.第一阶频率9088.6 Hz (虚线处)的峰值随着裂纹深度加深成比例的减小，很快激励出二阶频率9144.4 Hz，峰值的波动较小并且同样成比例减小深度加深至0.7 左右后，一阶频率甚至出现在9088.6 Hz 之前，且峰值极低.特别说明，此时的二阶频率9144.4 Hz 时开路条件下压电曲梁的一阶频率值.

图3 不同裂纹深度的电压-频率响应(L1=0.06 m)Fig.3 The voltage frequency responses with different crack depths(L1=0.06 m)

对于浅裂纹(0＜h′＜0.2)，激励出正向峰值的三阶频率;而在较深裂纹(h′＞0.2)的情况下，三阶频率的峰值远远小于浅裂纹时三阶频率的峰值，且随着裂纹加深频率向一阶频率方向偏移.因此，在不同裂纹深度时，压电曲梁的电压响应变化情况的趋势一致的.

裂纹损伤可能会出现在曲梁结构的任何位置，由图3 可知，裂纹深度的深浅引起的电压-频率响应变化的差别较大.因此在研究俘能器在裂纹位置变化下的电压响应时，分别考虑了浅裂纹和深裂纹两种情况.

如图4 所示，随着裂纹位置不断从曲梁固定端靠近曲梁中部，一阶固有频率逐渐增加.但在图4(a)中，一阶频率的峰值会不断增大，二阶频率9088.6 Hz和三阶频率9144.4 Hz 的峰值变大，且四阶频率也向右偏移.而在图4(b)中，压电俘能器在9088.6 Hz 前就被激励出了一阶和二阶频率，且频率皆向右偏移.当裂纹位置从曲梁中部不断靠近曲梁的自由端时，一阶频率和二阶频率仍是9088.6 Hz 和9144.4 Hz，且峰值变大，而三阶频率同样发生向右偏移.

图4 不同裂纹位置下的电压-频率响应Fig.4 The voltage frequency responses with different crack locations

总的来说，不同裂纹深度和位置的损伤曲梁压电俘能器的电压响应的趋势是一致的，不同的是被激励出的频率以及峰值.

因此想要检测结构的损伤情况，就要根据结构振动响应引起的俘能器的电压响应变化.一旦有裂纹出现，健康曲梁的一阶频率就不再是损伤曲梁的一阶频率，且随着裂纹加深该频率处峰值变小，由此可以判断出出现损伤时裂纹的深度情况.如若损伤曲梁被激励出的一阶频率小于健康曲梁的频率，说明曲梁的固定端到中部这一段结构中出现了裂纹，此时的裂纹深度也可以根据损伤曲梁是否在健康曲梁的一阶频率值前出现反向峰值的二阶频率.相反，如若损伤曲梁被激励出的一阶频率等于健康曲梁的一阶频率，则损伤出现的范围在曲梁的中部至自由端.至此，本文提出的根据结构振动引起的压电俘能器信号响应变化监测结构损伤的逆方法得以验证.

从上述裂纹几何参数对压电俘能器电压响应的分析中可以看出，在很多情况下固有频率都会发生偏移.因此，利用固有频率是否发生偏移这一条件仅能够监测结构中是否有损伤产生，并不能准确检测出结构的损伤情况.