[发明专利]基于FBG传感系统和二阶统计量的声发射定位系统及方法

摘要

本发明公开了基于FBG传感系统和二阶统计量的声发射定位系统及方法，未经平坦的ASE光源发出宽带光，依次经过边缘滤波器、第一耦合器、若干并列布置的第二耦合器传输给各自对应的光纤光栅传感器，光纤光栅传感器将设定波长的光反射回来经第二耦合器传输至光电转换电路，经放大器传至数据采集系统，数据采集系统记录信号变化；声发射源产生高频动态应力波，光纤光栅传感器的反射波长发生变化，通过光纤光栅传感器反射的光功率的变化实现声发射信号的解调；解调得到的声发射信号经过香农小波转换，提取窄带信号；根据二阶统计量计算出声发射源相对于参考传感器的方位角和距离，从而定位出声发射源的位置。本方法具有较好的声发射定位精度和实时性。

主权项

基于FBG传感系统和二阶统计量的声发射定位方法，其特征是，步骤如下：步骤(1)：光纤光栅传感器共同组成线性传感阵列；所述步骤(1)：设定中间的光纤光栅传感器为参考传感器，第l个传感器接收到的信号xl(t)表示为：其中，ωm＝‑2πd/λsinθm，其中，xl(t)表示第l个传感器接收到的信号，t表示时间，bm表示信号幅值，sm(t)表示信号源信号，Nl(t)表示第l个传感器的噪声信号，j表示复数形式，l表示传感器数，l的数量要大于声发射源数量，ωm＝‑2πd/λsinθm，d表示传感器间距，λ表示信号波长，θm表示声发射源相对于参考传感器的方向角；rm表示声发射源到参考传感器的距离；步骤(2)：声发射源产生高频动态应力波，在高频动态应力波的作用下，光纤光栅传感器的反射波长发生变化，通过光纤光栅传感器反射的光功率的变化实现声发射信号的解调；所述步骤(2)中光纤光栅传感器的反射波长发生变化由数据采集系统记录；未经平坦的ASE光源发出宽带光，依次经过第一耦合器、第二耦合器进入光纤光栅传感器，光纤光栅传感器将设定波长的光反射回来经第二耦合器传输至光电转换电路，光电转换电路将光信号转换为电信号，经放大器传至数据采集系统，数据采集系统记录全部信号变化；步骤(3)：对检测到的声发射信号进行香农Shannon小波转换，提取线性传感阵列的窄带信号；计算线性传感阵列的窄带信号的自相关函数；所述步骤(3)中自相关函数r的计算方法为：r=12p+1XXH---(2)]]>其中，X为光纤光栅传感器线性传感阵列接收的信号矩阵，X＝[x1(t),x2(t),…,x2p+1(t)]T；2p+1表示传感器数量；步骤(4)：定义二阶统计量；所述步骤(4)中定义二阶统计量r1(‑l,l)和r2(l+1,l)：r1(-l,l)=E(x-l(t)xl*(t))=Σm=1Mrsme-2jωml---(3)]]>r2(l+1,l)=E(xl+1(t)xl*(t))=Σm=1Mrsmej(ωm+φm)e2jφml---(4)]]>其中，rsm为声发射源信号的自相关函数，xl(t)为第l个传感器接收到的信号，M为声发射源数目；步骤(5)：根据步骤(4)中的二阶统计量定义，得到两个自相关函数矩阵R1和R2；所述步骤(5)：根据步骤(4)中的二阶统计量定义，得到两个p×p自相关函数矩阵R1和R2；R1和R2分别包含声发射源相对于参考传感器位置的方位角θ和声发射源相对于参考传感器位置的距离r；在矩阵R1和R2中，第(k，l)元素分别为：R1(k,l)＝r1(k‑l,l‑k)  (5)R2(k,l)＝r2(k‑l+1,k‑l)  (6)其中，r1和r2为步骤(4)所定义的二阶统计量；其中，R1和R2的向量形式为：R1＝A(ω)RsAH(ω)+σ2I  (7)其中，A(ω)＝[a(ω1) a(ω2) … a(ωM)]a(ωm)=1e2jωm...e2j(p-1)ωmT]]>Rs=diagrs1...rsm...rsM;]]>步骤(6)：对步骤(5)中的两个自相关函数矩阵R1和R2分别作特征值分解，得到两个自相关函数矩阵中的信号子空间ES1和ES2；所述步骤(6)：对R1和R2分别进行特征值分解得：R1＝[ES1 EN1]Λ1[ES1 EN1]H  (9)R2＝[ES2 EN2]Λ2[ES2 EN2]H  (10)其中，Λ1和Λ2分别是R1和R2的特征值矩阵，ES1和ES2分别是R1和R2的最大的M个特征值对应的信号子空间；EN1和EN2分别是R1和R2的L‑M个最小特征值对应的特征向量张成的噪声子空间；步骤(7)：进一步解析步骤(6)中的信号子空间ES1和ES2；所述步骤(7)：令α＝e2jω，和由此得两个多项式：F(α)=1α...αpδ1α...αpT=Σk=-p+1p-1δk′αk;---(11)]]>F(β)=1β...βpγ1β...βpT=Σk=-p+1p-1γk′βk;---(12)]]>利用求根公式对公式(11)和公式(12)进行求根，得到α＝α1,α2,…,αM和β＝β1,β2,…,βM；α和β中包含方位角θ和距离r信息；其中，α＝e2jω，ω＝‑2πd/λsinθ，ES1和ES2分别是R1和R2的最大的M个特征值对应的信号子空间；p表示传感器数量，M声发射信号源数；步骤(8)：根据参数配对和步骤(7)的解析结果，通过root‑music算法解出声发射源相对于参考传感器的方位角和距离；步骤(8)：根据参数配对和步骤(7)中的α和β，通过root‑music算法：θi＝sin‑1(angle(αi)λ/(4πd))和ri＝2πd2cos2θi/(λangle(βi))解出声发射源相对于参考传感器的方位角和距离；其中，θi表示第i个声发射源相对于参考传感器的方位角，angle(αi)表示αi的相位参数，λ表示信号波长，d表示传感器间距，ri表示第i个声发射源相对于参考传感器的距离，αi为步骤(7)中F(α)的第i个根，angle(βi)表示βi的相位参数，βi为步骤(7)中F(β)的第i个根；所述方法所采用的系统，包括：未经平坦的ASE光源，未经平坦的ASE光源发出宽带光，依次经过边缘滤波器和第一耦合器，第一耦合器将光信号分别传输给若干并列布置的第二耦合器，每个第二耦合器将信号传输给各自对应的光纤光栅传感器，光纤光栅传感器等间距排列并粘贴在被检测机构上；光纤光栅传感器将设定波长的光反射回来经第二耦合器传输至光电转换电路，光电转换电路将光信号转换为电信号，经放大器传至数据采集系统，数据采集系统记录全部信号变化；声发射源产生高频动态应力波，在高频动态应力波的作用下，光纤光栅传感器的反射波长发生变化，通过光纤光栅传感器反射的光功率的变化实现声发射信号的解调；解调得到的声发射信号经过香农Shannon小波转换，提取窄带信号；根据二阶统计量的root‑music算法计算出声发射源相对于参考传感器的方位角和距离，从而定位出声发射源的位置。