专利名称：：一种转子碰摩声发射源定位方法
技术领域：
：本发明涉及一种声源定位方法，特别涉及一种转子碰摩声发射源定位方法。
背景技术：
：在基于声发射技术的转子碰摩故障诊断中，不仅可以由声发射信号判断碰摩的发生，而且可通过声发射源定位技术快速找到碰摩发生的位置，为分析故障原因进而排除故障提供十分重要的信息。但在转子系统结构中，从碰摩源到传感器之间常常是一段非连续非单一介质的复杂体传播路径，碰摩激励的多模态声发射波传播速度不同，并且在传播过程中受散射、边界条件、频散效应、模式转换等多方面影响，信号畸变严重，如果采用传统的线或面时差定位法，必然带来严重偏差。对于复杂体的声源定位，主要方法有基于时延估计的定位算法和基于空间谱估计的定位算法。基于时延估计的算法通过计算不同传感器信号的时域互相关获得样本时延，再根据声信号传播理论和阵元位置获得声源的空间位置，它要求获得的传感器时延必须十分精确，否则对算法性能影响很大。基于空间谱估计的声源定位算法以MUSIC算法为典型代表，相比基于时延估计的声源定位方法，它的定位精度更高，但是其算法复杂度也相应提高。
发明内容本发明要解决的技术问题是针对现有技术存在的缺陷提出一种转子碰摩声发射源定位方法。本发明一种转子碰摩声发射源定位方法，其特征在于包括如下步骤-(1)采用碰摩声发射试验装置获得声发射信号；(2)建立基于步骤(1)所述的能量衰减的声发射信号的传播模型由N个传感器组成的传感器阵列中，第/个传感器和第y个传感器的能量比如下当0<~#1，所有满足上式的声源坐标r(O都位于中心为c,"半径为^的d维超球上<formula>formulaseeoriginaldocumentpage5</formula>n个传感器组成w(iv—"1对能量比传感器，对于任意超球数m《^^—I使用最22小方差标准构造m维超球的估计目标代价函数<formula>formulaseeoriginaldocumentpage5</formula>声发射源的位置argminJ"(r)=argmin"||r-cm||_/m|，即目标代价函数JO0最小时所对应的声源位置r就是声发射源的位置^，其中r(O为声源在时刻的坐标，乂o)为第f个传感器在,时刻检测得到的声发射信号的能量，1(0为第7'个传感器在r时刻检测得到的声发射信号的能量，A为第/个传感器模型累积误差及观测噪声引入误差的均值，/^为第y个传感器模型累积误差及观测噪声引入误差的均值，r,为第/个传感器的坐标，g,为第/个传感器的增益系数，。为第J'个传感器的坐标，^为第J个传感器的增益系数，"为能量衰减因子，w维超球的中心和半径分别为c^、pm，n、/、y和附都为自然数，下同；(3)利用自适应次梯度投影方法更新声发射源的位置^;:基于凸函数g(r)的半空间/T(rj:={rer":(r-rjft+)S0}，那么声发射源的位置^对半空间/T(rJ的投影为p0;):=其中Vg,O;)为声发射源的位置iv的梯度算子，得到声发射源的位胃l^的&代更新为:r*+1=R+4(;(rt)(rJ-0，ir为w对能量比传感器构成的声源位置r的集合，r为声发射源的发射周期，々为松弛系数，其满足^e。本发明优点和效果在于1.根据声发射信号传播特性建立的信号能量衰减模型，将目标声源位置限定于圆心和半径都为传感器对的能量比函数的超球，并将定位问题转化为估计问题，使得定位问题的处理变得简单。2.采用自适应次梯度投影系统估计方法逼近声源位置，该方法收敛性能好，并且具有很高的定位精度。图1——转子碰摩试验台结构图2——位于(-l,O)和(l，O)的传感器对声源定位超球图3—传感器接收到的实测声发射信号图4——次梯度投影定位方法估计声源坐标图5——广义相关时延估计法估计时间延迟图6^互功率谱相位法估计时间延迟图7—传感器接收到的实测声发射信号图8——次梯度投影定位方法估计声源坐标图；图9——广义相关时延估计法估计时间延迟图；图10——互功率谱相位法估计时间延迟图；图11——三种定位方法的性能比较图。具体实施例方式1.获得声发射信号如图1所示的转子碰摩试验台提取声发射信号。该试验台通过电动机输入电压来调节转速，电动机与转轴之间为半挠性连轴节，转子支承为具有滑动轴承的轴承座。转子系统动静部件的碰摩通过一安装在转子台底座上可移动的碰摩装置来模拟实现。碰摩装置安装在轴承座1、2之间，装置侧面的螺孔上安装一个可伸縮的螺栓，沿转轴径向对着转轴中心，6通过调节螺栓产生转子碰摩。图中标号对应装置如下l-电机；2-增速箱；3-联轴器；4-轴承；5-轴承座1;6-碰摩装置；7-底座；8-转盘；9-轴；10-轴承座22.建立基于能量衰减的声发射信号的传播模型，将定位问题转化为系统估计问题在声发射信号传播过程中，其能量与声源距离按反比关系衰减。假设空间有一点声源，其坐标表示为r，由N个传感器组成传感器阵列，其中第/个传感器位置为V",2,…，iV。在？时刻第/个传感器测量到的声源能量表示为,、<卜0"(1)lr(卜()-r,l其中(表示从声源传播到传感器的时延，s(0表示时间间隔^发射的能量，r(O为声源在f时刻的坐标，r,为第/个传感器的坐标，g,为第/个传感器的增益系数，"(*2)为能量衰减因子，。为第J'个传感器的坐标，^为第y个传感器的增益系数，"为能量衰减因子，m维超球的中心和半径分别为^、/7m，s,(O为模型累积误差及观测噪声引入的误差的共同结果。假定^(0为独立同分布归一化随机变量，均值;",，方差为W，则X(O的概率密度函数为iV(,g气其似然函数(,)，r(o)=/(y。(0,…,;v—々)IJ2,W0,KW)，-y"「WV|r(0-r,rOC6Xp1w-丄》^=0(2)最大似然估计的目的是找到声源能量和位置^(0,r(OM吏似然函数最大，等价于最小化对数似然函数少,(0H，|r(0_r,(3)在不知道声源能量s(O的情况下，弓l入能量比以消去W0的影响。在噪声《(0以其均值估计A代替的情况下，第/个传感器和第j'个传感器的能量比计算如下:(4)(W)-A)化(0-A)这里对O<^*1的情况下，所有满足上式的声源坐标r(O都位于中心为cy.、半径为&的d维超球上，该超球的定义公式为|r(0_%|=《f一^Sr'_r)(5)1-《"y该超球被称为目标定位超球。因此，使用一对传感器的能量比，可以将目标声源位置限制于圆心和半径都为能量比的函数的超球上。图2表示了在不同能量比情况下的目标定位超球示意图。如果使用更多的传感器，则可以决定更多的超球。如果所有传感器都接收来自同一个声源的信号，则相应的目标定位超球一定在特定点相交。此即能量比声源定位的基本原理。由于该定位方法基于能量比而非能量，故即使在声源能量显著变化时其定位精度也不会受到显著影响。假定N个声源传感器可以组成—0=n对能量比传感器，对于任意M《翠—"，22可以使用最小方差标准构造估计目标代价函数/(r)=X||r-c』-&|，(6)其中^为超球数。那么声源位置r的定位问题可以转化为如下式所表示的估计问题M2argmjnJ(r)=argri^inZ||r—cm|->om|(7)m=l即，目标代价函数最小时所对应的r就是声发射源的位置。3.利用自适应次梯度投影方法估计出声发射源的位置自适应滤波器法是最常见的系统估计方案。这种方案的核心是设计自适应滤波器，使其系数逐渐逼近被估计系统。由于声发射源的位置可能变化，故自适应滤波器应能快速跟踪这种变化。该方法利用次梯度取代了正交投影运算，从而简化了算法并减小了计算量。由式(7)，定义集合A^卜R":Ur-cgi、/^，因此，声源位置的估计等价于寻找交门A定义凸函数:g(r):=r—c』—pM,VreR"(8)和满足下列约束的凸集:={reR":g(r)20}(9)由于凸集&以较大的概率包含了实际声源位置，那么对Z的估计问题就转变成了向凸集G的投影问题。定义梯度算子:t:=Vg(r)=2(r—O,VreR"，(10)那么可以定义基于该凸函数的半空间zt(rJ:={reR":(r-rtft+g(rjS。}，(11)满足当rtgQ时有Q(/Jczt(rJ和rg/t(rJ，因此对Q(户)的投影可以扩展为对半空间/T(rJ的投影，投影公式如下式所示9尸""(叫^;+:，)2V^)(12)那么，R的迭代更新公式为(ri)(rJ-rj，(13)i"为w对能量比传感器构成的声源位置r的集合，r为声发射源的发射周期，这里松弛系数4应满足^e。为了验证本发明提出的定位方法的优势，选取两种最常用的声源定位方法——广义相关时延估计法和互功率谱相位法进行比较。实施例l:传感器1在轴承座1，传感器2在轴承座2，声发射源在转轴上，距离传感器l=20cm。图3为传感器接收到转子摩擦时的声源信号，采样频率2MHz，点数32768。为确定摩擦声源的具体位置距离，两个传感器分别位于(xl，yl)=(0，0)和(x2，y2)=(43，0)。采用采用本发明提出的定位方法、广义相关时延估计法和互功率谱相位法3种方法分别仿真的结果如图4、图5和图6所示。图4估计声源在20.3684cm处，图5中互相关函数在采样点31956处取得峰值，而图6中互相关函数取得峰值在采样点31824处。实施例2:传感器l-轴承座l，传感器S2-底座，AE源在转轴上，距离轴承座147cm。图7为传感器接收到的转子摩擦时的声源信号，采样频率lMHz，点数16384。为确定摩擦声源的具体位置距离，2个传感器分别位于(xl，yl)=(0，0)和(x2，y2)=(35，5)。同样采用本发明提出的定位方法、广义相关时延估计法和互功率谱相位法3种方法分别仿真的结果如图8、图9和图10所示。图8估计声源在17.0799cm处，图9中互相关函数在采样点325处取得峰值，而图10中互相关函数取得峰值在采样点451处。由于广义相关时延估计法和互功率谱相位法属于时间延迟方法，利用延迟值可以定位出声源的坐标。表1列出了三种定位算法的性能比较。表l<table>tableseeoriginaldocumentpage11</column></row><table>从表1中可以看出次梯度投影算法比时延估计算法有更高的定位精度。比较图4和图8可以发现，经过5次迭代以后，本发明提出的方法就已收敛，因而同时具有较低的复杂度。权利要求1.一种转子碰摩声发射源定位方法，其特征在于包括如下步骤(1)采用碰摩声发射试验装置获得声发射信号；(2)建立基于步骤(1)所述的能量衰减的声发射信号的传播模型由N个传感器组成的传感器阵列中，第i个传感器和第j个传感器的能量比如下当0<kij≠1，所有满足上式的声源坐标r(t)都位于中心为cij、半径为ρij的d维超球上N个传感器组成对能量比传感器，对于任意超球数使用最小方差标准构造m维超球的估计目标代价函数声发射源的位置即目标代价函数J(r)最小时所对应的声源位置r就是声发射源的位置rk，其中r(t)为声源在t时刻的坐标，yi(t)为第i个传感器在t时刻检测得到的声发射信号的能量，yj(t)为第j个传感器在t时刻检测得到的声发射信号的能量，μi为第i个传感器模型累积误差及观测噪声引入误差的均值，μj为第j个传感器模型累积误差及观测噪声引入误差的均值，ri为第i个传感器的坐标，gi为第i个传感器的增益系数，rj为第j个传感器的坐标，gj为第j个传感器的增益系数，α为能量衰减因子，m维超球的中心和半径分别为cm、ρm，N、i、j和m都为自然数，下同；(3)利用自适应次梯度投影方法更新声发射源的位置rk基于凸函数g(r)的半空间H-(rk)＝{r∈Rn(r-rk)Tt+g(rk)≤0}，声发射源的位置rk对半空间H-(rk)的投影为其中为声发射源的位置rk的梯度算子，得到声发射源的位置rk的迭代更新为Rn为n对能量比传感器构成的声源位置r的集合，T为声发射源的发射周期，λk为松弛系数，其满足λk∈。全文摘要本发明公布了一种转子碰摩声发射源定位方法。本发明通过转子碰摩试验台获得声发射信号，建立基于能量衰减的声发射信号的传播模型，并且将定位问题转化为系统估计问题，利用自适应次梯度投影方法进行系统估计，从而估计出声发射源的位置。本发明具有计算处理简单，收敛性能好，定位精度高等优点，从而有效地应用于声发射源定位。文档编号G01S5/00GK101477194SQ20091002508公开日2009年7月8日申请日期2009年2月17日优先权日2009年2月17日发明者包永强,力赵,邓艾东申请人:东南大学