一种光纤布拉格光栅传感网络结构多裂纹定位方法
【专利摘要】本发明涉及一种光纤布拉格光栅传感网络结构多裂纹定位方法，属于光纤传感的信号处理及应用技术领域。该方法包括以下步骤：S1：搭建具有温度补偿功能的FBG应变花传感网络检测声发射信号；S2：通过信源估计方法改进分离算法，对声发射信号进行分离，得到定位所需的单个声发射信号；S3：通过滤除非裂纹源算法对分离得到的单个声发射信号集合进行干扰源滤除，得到所需的裂纹声发射信号；S4：通过设定经验阈值对FBG检测到的单个声发射信号设计判定声发射事件开始和结束时刻的策略，进而求取各个声发射事件导致的应变花FBG波长漂移量；S5：建立FBG传感器横向效应补偿因子模型H，引入应变解耦模型，根据应变花定位原理对每个裂纹声发射事件进行定位。本方法考虑了外部干扰源问题，能够用于检测与定位结构多裂纹情况。
【专利说明】
-种光纤布拉格光栅传感网络结构多裂纹定位方法
技术领域
[0001] 本发明属于光纤传感的信号处理及应用技术领域，设及一种光纤布拉格光栅传感 网络结构多裂纹定位方法。
【背景技术】
[0002] 随着我国经济及科学技术的发展，越来越多的大型设备，如航空飞行器、舰船、海 洋平台、轨道、桥梁等逐渐服役用来服务于社会，然而运些设备在服役期间不可避免地受到 各种各样的冲击、疲劳等因素的影响，进而产生结构裂纹。运些结构裂纹如若不及早发现， 便会随着服役发生扩展，不仅会对结构部位产生损坏，极有可能影响设备的安全运行，对人 们的生命财产产生巨大的威胁，因此，设计一种可行的结构裂纹定位传感系统意义重大。
[0003] 光纤布拉格光栅(F;Lber Bragg grating,FBG)传感器具有体积小、灵敏度高、抗电 磁干扰、耐腐蚀、可靠性高、易串接复用和不影响宿主材料结构特性等优点，广泛应用于航 空航天、舰船工业、海洋平台、上木工程及医学工程等领域。
[0004] 目前，FBG传感技术已被国内外很多研究者应用于结构单个裂纹的检测，而对结构 多裂纹定位研究上，由于观测信号是多个未知数目的源信号的混叠，常规的单裂纹定位技 术便不再适用。依据FBG对应变敏感的特性，针对结构多裂纹检测检测方法，根据其原理可 分为谐波定位法和声发射(Acoustic emission,AE)定位法。谐波定位法主要是依据把裂纹 视作造成谐波的来源，通过激励源激励结构材料，对采集数据进行时间反转技术，进而对裂 纹造成的谐波源进行聚焦，该方法局限在数值仿真上，且信号采集困难，而不当的激励方式 也可能会对结构产生一定的影响;声发射定位法主要依据是裂纹在产生和扩展的时候会产 生AE现象，进而通过传感器接收AE应力波，再通过定位算法进行定位，概括起来主要有基于 独立分量分析的多裂纹定定位法和基于最小方差无失真响应的多裂纹定定位法。基于独立 分量分析的定位法存在着默认裂纹源数目已知的缺陷;而基于最小方差无失真响应的定位 法在采集信号上，每个信号对应一个FBG传感器，局限于FBG传感网络上。并且运两种方法都 没有考虑外部干扰源的问题。
【发明内容】
[0005] 有鉴于此，本发明的目的在于提供一种光纤布拉格光栅传感网络结构多裂纹定位 方法，该方法能够用于检测与定位结构多裂纹情况。
[0006] 为达到上述目的，本发明提供如下技术方案：
[0007] -种光纤布拉格光栅传感网络结构多裂纹定位方法，该方法包括W下步骤：
[000引S1:搭建具有溫度补偿功能的FBG应变花传感网络检测声发射信号；
[0009] S2:通过信源估计方法改进分离算法，对声发射信号进行分离，得到定位所需的单 个声发射信号；
[0010] S3:通过滤除非裂纹源算法对分离得到的单个声发射信号集合进行干扰源滤除， 得到所需的裂纹声发射信号；
[0011] S4:通过设定经验阔值对FBG检测到的单个声发射信号设计判定声发射事件开始 和结束时刻的策略，进而求取各个声发射事件导致的应变花FBG波长漂移量；
[0012] S5:建立FBG传感器横向效应补偿因子模型H，引入应变解禪模型，根据应变花定位 原理对每个裂纹声发射事件进行定位。
[0013] 进一步，在步骤S4中，通过设置经验阔值来对FBG应变花结构检测到的裂纹声发射 事件进行判定其开始和结束时刻，再通过设定其最大波长漂移量时刻，进而确定波长漂移 量，具体包括W下步骤：
[0014] S41:开始时刻：通过设置经验阔值，其与FBG检测信号波形的交点即为开始时刻， 且对Ξ个FBG传感信号做同样处理；
[001引S42:最大波长漂移量时刻:在开始时刻后，对各个FBG传感信号进行一次求导，第 一个零点所对应的时刻即为对应传感器的峰值时刻，记为ti(i = 1，2，3为FBG传感器序号）， 则声发射事件的峰值时刻为
[0016]
[0017] S43:结束时刻：当Ξ个传感器波长值，都在我们设置的阔值内时的时刻设置为声 发射事件结束时刻。
[0018] 进一步，在步骤S5中，采用粘合剂(如环氧树脂等)将一不同中屯、波长的FBG传感器 单端固定在被测结构位置上，实现对应变花结构的溫度补偿。
[0019] 进一步，在步骤S5中，通过构建FBG传感器横向效应补偿因子模型Η来减少FBGii量 应变误差，由FBG传感原理推导出Η如下式所示：
[0020]
[0021] 式中：Pe和P^e分别是指在轴向应力和横向应力作用下的FBG有效弹光系数，且
P11、P12分别为单模光纤的弹光系数，μ日为光纤的泊 松比，neff是纤忍折射率。
[0022] 进一步，在步骤S5中，引入解禪模型，使得FBG二维应变之间的禪合作用被消除，W 提高测量精度，减小定位误差。根据解禪原理，推出FBG二维应变解禪模型，如下式所示：
[0023]
[0024] 式中：ε功FBG粘贴方向的应变值，ε/功FBG的测量值，i = 1，2，3为传感器序号，Η为 FBG横向效应补偿因子。
[0025] 本发明的有益效果在于:本发明提供的光纤布拉格光栅传感网络结构多裂纹定位 方法考虑了外部干扰源问题，能够用于检测与定位结构多裂纹情况。
【附图说明】
[0026] 为了使本发明的目的、技术方案和有益效果更加清楚，本发明提供如下附图进行 说明：
[0027] 图1为本发明所述方法的流程示意图；
[00%]图2为FBG直角应变花结构示意图；
[0029]图3为剔除干扰源方法流程图；
[0030]图4为已知源数目的信号分离方法框图；
[0031 ]图5应变解禪模型示意图；
[0032] 图6FBG1-3光谱检测开始、结束时刻图。
【具体实施方式】
[0033] 下面将结合附图，对本发明的优选实施例进行详细的描述。
[0034] 本实施例包括两个部分，第一部分是信号处理，即是对FBG传感信号进行处理，W 得到我们定位所需的FBG波长漂移量数据。第二部分是多裂纹定位，也即是根据裂纹声发射 应力波导致的FBG波长漂移量，采用具有溫度补偿功能的FBG直角应变花模型对裂纹进行定 位。具体的设计方法包括如下步骤：
[0035] 1.信号处理。首先，通过声发射源数目估计方法（如本实施例选用交叉验证思想） 来改进盲分离算法，进而实现多声发射信号的分离;其次，采用合适的权值算法(如本实施 例选用质屯、加权算法)构建滤除非裂纹源算法，得到裂纹声发射信号；最后设计分布式FBG 动态传感信号波长漂移量提取策略，对单个裂纹声发射现象导致的FBG波长漂移量进行提 取。
[0036] 1.1多声发射信号分离。首先，采用交叉验证技术对声发射信号源数目进行估计； 其次，采用盲均衡算法对已知数目的声发射信号进行分离。
[0037] 1.2非裂纹源滤除。采用质屯、加权算法构建滤除非裂纹源算法，根据声发射应力波 特征参数，构建声发射应力波综合指标，当综合指标落在[μ-3σ ,μ+3σ]的区间外时，剔除该 干扰信号。其中σ表示特征参数的标准偏差，μ表示特征参数的平均值。
[0038] 1.3FBG应变花定位模型波长漂移量提取策略。由于本方法是个串接的FBG组 成的应变花为对象，因此需要确定同一声发射事件造成的Ξ个应变FBG信号何时为提取波 长值时刻。通过设置经验阔值，来确定每一个声发射事件造成的Ξ个应变FBG信号的开始时 亥IJ、峰值时刻和结束时刻，实现对动态FBG信号波长漂移量的求取。
[0039] 2.多裂纹定位。采用带溫度补偿功能的FBG应变花结构(如实施例采用直角应变花 结构)对声发射信号进行检测，并构建裂纹感应与传感网络，将多个FBG直角应变花串接实 现分布式传感检测，通过FBG反射光谱中屯、波长漂移量来对结构裂纹实时检测。
[0040] 2.1定位模型。将FBG传感器和力学测量中测量应力波主应变方向的应变花原理相 结合，采用Ξ个FBG传感器构建一组FBG直角应变花结构模型，即可测量一个应力波的主应 变方向角，当布置两组不同位置的FBG直角应变定位结构模型，两条相交的直线即可确定裂 纹的具体位置。
[0041 ] 2.2溫度补偿。在FBG直角应变花结构上增加一单端粘贴且中屯、波长不同的FBG传 感器。
[0042] 2.3构建横向效应补偿因子模型。通过对平面应变条件下的FBG受力分析，构建Η的 理论模型。
[0043] 2.4FBG应变解禪。引入控制理论中的解禪原理对FBG传感器在二维平面条件下的 横纵应变进行解禪。
[0044] 图1为本发明所述方法的流程示意图，如图所示，本方法具体包括：
[0045] 1、获得光谱
[0046] 选取四个中屯、波长不同的FBG传感器，按照图2方式进行布置，即一组FBG应变花。 其中FBG1-FBG3为定位传感器，使用粘合剂（如环氧树脂等），将FBG1-FBG3完全粘贴在相应 位置上，FBG4为溫度补偿传感器，且为单端固定。
[0047] 2、多裂纹声发射信号分离方法 [004引2.1裂纹数目估计
[0049]本发明引入基于交叉验证技术的时变信源数估计方法来进行盲源数估计，基本思 想如下：
[(K)加 ]
(1)
[0051]式中：n(t)为裂纹估计数。i = l，2，L，m，m为FBG检测到的信号个数。trace(g)为矩 阵求迹运算。>//''=曲g(c-0#) Χ = χχΤ为观测信号X的协方差矩阵。#=昭ν，Λι为矩阵C前i个特 征值组成的对角矩阵，Ui的列向量是矩阵C前i个特征值对应的特征向量；而计算…的 S =民,义，.对角矩阵怎的对角线与Am的交叉，且两个矩阵在对角线上的元素排序的相反的。为 了满足估计时变裂纹数目的要求，C采用增量式(2)、递推式(3)、（4)计算：
[0055] 2.2裂纹数目已知的信号分离技术
[0056] 图3中，S化)是裂纹声发射信号，G是未知信道，X化)是观测信号，W化)是分离矩阵， 如3是输出信号y化)的广义能量函数，f(g)是非线性函数，常见的有f(y)=y，f(y)=y3，f(y) =化nh(y)等，由设备所处的特殊环境决定J脚为源信号的期望信号。e化)是分离矩阵的调 节因子。
[0057] 采用最小均方(LMS)算法最小化输出信号y化)的广义能量目标函数
[0化引
(S)
[0059] 可得分离矩阵W化)的迭代公式
[0060]
(.紛
[0061] 式中为迭代系数。
[0062] 而非线性函数f(g)，由设备的服役环境决定，采用盲解卷结果与期望信号的相似 程度，即相似系数来检验，表示如下：
[0063]
巧）
[0064] 相似度系数最大的非线性函数确定为本系统的f(g)。
[0065] 同时采用估计值和输入值之间的误差e(k)来对分离矩阵进行补偿，图4中的盲解 卷滤波器补偿系数e化)的目标函数可表示为：
[0066]
(8)
[0067]通过常模算法(CMA)，最小化上式，可得w(k)的自适应迭代式 [006引
巧）
[0069] 式中:k = 2，y为自适应迭代式的步长，根据实验分离效果确定μ的最优值。
[0070] 3、非裂纹源滤除
[0071] 由于每个应力波可W用少数几个特征参数来表示，本发明采用复杂度较低的质屯、 加权算法求取权重，求取第i个声发射事件中的第j项参数wj，设Ej是每个应力波的第j个特 征值检测值，得应力波的综合指标
[0072]
(10)
[0073] 进而利用30原理对非裂纹源进而滤除。其具体流程如图4所示。
[0074] 4、应力波提取波峰策略
[0075] 经非裂纹源滤除、信号分离等步骤处理过的声发射应力波，由于是动态信号，为避 免手动提取FBG波长漂移量，因此在运用应变花求取声发射应力波主应变方向时需要确定 出输入波长的值，此时需要确定提取波长值的时刻，处理步骤如下：
[0076] 开始时刻：WFBG1为例，通过设置经验干扰波长变化量Δ λ，其与波形的交点时刻 即作为FBG1开始时刻，对Ξ个FBG波形进行同样处理，设置最早时刻作为开始时刻；
[0077] 波长漂移量确定时刻：在开始时刻后，通过一次求导，第一个零点即为所求时刻， 记为ti，同理，对FBG2，FBG3做同样处理，即可得同一次裂纹声发射事件内各自的波长漂移 量最大时刻，即定位模型所需的波长值，则τ时刻对应的FBG1-3的各自波长值，其中
[007引
（21)
[0079] 然后提取τ时刻FBGi (i = 1，2，3)的波长漂移量，求取相应的应变值，实现裂纹定 位。
[0080] 结束时刻：当Ξ个传感器波长值，都在我们设置的阔值内的时刻设置为声发射事 件结束时刻。如图5所示，其中矩形点表示开始标识，圆点表示结束标识。
[0081] 5、多裂纹定位
[0082] 5.1FBG直角应变原理
[0083] 图2所示为带有溫度补偿的FBG直角应变花示意图。传感器粘贴处的Ξ个应变分量 εχ,εγ，丫 xy和FBG测量应变之间的关系
[0084]
山)
[0085] 主应变方向由式（12)给出。
[0086]
(似
[0087] 式中α为主应变的方向和X轴主方向的夹角，ε〇。，645。和690。分别是相应下标角度方 向上的测量应变值。
[0088] 5.2FBG传感器横向效应修正
[0089] 光纤化agg光栅由于具有横向效应，因此需要对式（12)中的ε0。，ε45。和ε90。进行横向 效应补偿，平面应变条件下两个相互垂直的FBG横向效应补偿表达式
[0090]
(14)
[0091] 式中ει、ε2分别为FBGi方向、FBGj方向应变的真实值，ε/ι、ε/2分别为FBGi方向、 FBG巧向的测量值巧中i = 1，j = 3或i = 2，j = 4)，K = (1 -μ地)/(l-tf)，μ。为光纤的泊松比， Η=Κβ/肚为横向效应补偿因子，Kb为其横向灵敏系数，肚为轴向灵敏系数。
[0092] 而横向效应因子Η模型由下式给出
[0093]
牌
[0094] 式中Pe和P^e分别是指在轴向应力和横向应力作用下的FBG有效弹光系数，且
P11、P12分别为单模光纤的弹光系数，μ日为光纤的泊 松比，neff是纤忍折射率。
[00巧]5.3应变解禪模型
[0096] 由于平面应变下，FBG对二维应变（ει、ε2)存在交叉敏感，即禪合关系，采用控制理 论中的解禪原理来实现应变的解禪，其系统解禪模型如图6所示。图中输入信号ε/ι、ε/2分别 为FBGi、FBGj方向应变的测量值，输出信号ει、ε2分别为FBGi、FBGj的方向应变真实值;Gn (S)、Gi2(S)、G21(S)和 G22 ( S )为禪合对象,结合式(3 )可知,Gll(S)二 K、Gi2(S)二-kH、G21 ( S )二- 地及622(3)=10,〇12(3)、〇21(3)为解禪因子刷6横向效应解禪模型由下式给出。
[0097] ε = (1-Η^)κε' (16)
[009引式中ε是经解禪得到的应变值，κ=(1-μ曲)/(1-Η2)，Η为横向效应补偿因子，ε'为测 量应变值。
[0099] 5.4多裂纹定位
[0100] 裂纹产生的应力波经过检测位置时，通过一组应变花可W确定一个应力波的主应 变方向角，当在不同位置放置两个应变花时，可W得到两个方向角，其交点便可确定裂纹位 置。对于多裂纹情况，通过上述信号处理，也即把多裂纹分离成了单个裂纹的集合，再采用 FBG直角应变定位原理实现多裂纹定位。
[0101] 最后说明的是，W上优选实施例仅用W说明本发明的技术方案而非限制，尽管通 过上述优选实施例已经对本发明进行了详细的描述，但本领域技术人员应当理解，可W在 形式上和细节上对其作出各种各样的改变，而不偏离本发明权利要求书所限定的范围。
【主权项】
1. 一种光纤布拉格光栅传感网络结构多裂纹定位方法，其特征在于：该方法包括以下 步骤： Sl:搭建具有温度补偿功能的FBG应变花传感网络检测声发射信号； S2:通过信源估计方法改进分离算法，对声发射信号进行分离，得到定位所需的单个声 发射信号； S3:通过滤除非裂纹源算法对分离得到的单个声发射信号集合进行干扰源滤除，得到 所需的裂纹声发射信号； S4:通过设定经验阈值对FBG检测到的单个声发射信号设计判定声发射事件开始和结 束时刻的策略，进而求取各个声发射事件导致的应变花FBG波长漂移量； S5:建立FBG传感器横向效应补偿因子模型H，引入应变解耦模型，根据应变花定位原理 对每个裂纹声发射事件进行定位。2. 根据权利要求1所述的一种光纤布拉格光栅传感网络结构多裂纹定位方法，其特征 在于:在步骤S4中，通过设置经验阈值来对FBG应变花结构检测到的裂纹声发射事件进行判 定其开始和结束时刻，再通过设定其最大波长漂移量时刻，进而确定波长漂移量，具体包括 以下步骤： S41:开始时刻：通过设置经验阈值，其与FBG检测信号波形的交点即为开始时刻，且对 三个FBG传感信号做同样处理； S42:最大波长漂移量时刻:在开始时刻后，对各个FBG传感信号进行一次求导，第一个 零点所对应的时刻即为对应传感器的峰值时刻，记为ti( i = 1，2，3为FBG传感器序号），则声 发射事件的峰值时刻为S43:结束时刻：当三个传感器波长值，都在我们设置的阈值内时的时刻设置为声发射 事件结束时刻。3. 根据权利要求1所述的一种光纤布拉格光栅传感网络结构多裂纹定位方法，其特征 在于：在步骤S5中，采用粘合剂将一不同中心波长的FBG传感器单端固定在被测结构位置 上，实现对应变花结构的温度补偿。4. 根据权利要求1所述的一种光纤布拉格光栅传感网络结构多裂纹定位方法，其特征 在于:在步骤S5中，通过构建FBG传感器横向效应补偿因子模型H来减少FBG测量应变误差， 由FBG传感原理推导出H如下式所f式中：Pe和P7e分别是指在轴向应力和横向应力作用下的FBG有效弹光系数，且 ',,Π '?,υ] _.Pii、Pi2分别为单模光纤的弹光系数，μO为光纤的泊 松比，neff是纤芯折射率。5. 根据权利要求1所述的一种光纤布拉格光栅传感网络结构多裂纹定位方法，其特征 在于:在步骤S5中，引入解耦模型，使得FBG二维应变之间的耦合作用被消除，以提高测量精 度，减小定位误差。根据解耦原理，推出FBG二维应变解耦模型，如下式所示： Ei= (I-H2)Ke7 i 式中：￡i为FBG粘贴方向的应变值，'为FBG的测量值，i = I，2，3为传感器序号，H为FBG 横向效应补偿因子。
【文档编号】G01N29/14GK105842345SQ201610234657
【公开日】2016年8月10日
【申请日】2016年4月15日
【发明人】陈勇, 刘保林, 刘焕淋
【申请人】重庆邮电大学