基于Teager能量算子和样本熵的板结构冲击监测方法技术

本发明专利技术提出一种基于Teager能量算子和样本熵的板结构冲击监测方法。本发明专利技术方法包括以下步骤：步骤一：分布式光纤布拉格光栅传感网络布置；步骤二：待测冲击的光纤布拉格光栅传感器响应信号采集及Teager能量算子处理；步骤三：经处理冲击响应信号的样本熵计算；步骤四、确定待测点所在区域，并计算区域内样本冲击点响应信号的样本熵；步骤五：冲击点到传感器的距离与相应样本熵差值之间的关系模型构建；步骤六：基于三角区域法，计算待测冲击点的具体坐标位置。本发明专利技术通过Teager能量算子提高冲击响应信号质量，并通过样本熵定量估计冲击响应信号的复杂度，可应用于航空航天领域四边固支板结构的冲击定位，无需大量先验知识并可以在线监测辨识冲击位置。

【技术实现步骤摘要】
基于Teager能量算子和样本熵的板结构冲击监测方法
本专利技术属于结构健康监测的冲击监测
，尤其涉及能抑制干扰和提高时间分辨率的Teager能量算子和定量估计信号复杂度的样本熵的板结构冲击监测方法。
技术介绍
现代诸多大型结构如航空航天飞行器、高层结构、新型桥梁、大跨度网架结构等在复杂服役环境中将受到设计载荷作用以及各类突发性外在因素如太空垃圾、飞鸟撞击、冰雹冲击、雷击和维修过程中人为低能量冲击，而这些冲击载荷的施加有可能对机械结构造成基体断裂、分层和穿孔等一些不能目检的损伤，给航空航天飞行器带来严重安全隐患，减少了飞行器使用寿命。因此有必要对航天航空飞行器进行安全检测，而常规检测主要是无损检测方法如超声波、X射线、敲击、电涡流射线、热应力场以及电位测量等，相关检测设备对使用工况要求较高，人员费用高、不适于实现在线监测。目前，国内外对冲击载荷位置识别方法主要有频域识别法、时域识别法、人工神经网络法和时间差定位等。不少学者将光纤光栅传感技术应用到结构冲击载荷位置识别的研究，并取得一定成果。由于光纤光栅解调仪采样频率较低，基于光纤光栅的冲击监测方法大都采用神经网络、支持向量机等模式辨识算法，但此类方法缺点是需要事先采集高密度不同冲击样本进行网络训练，这不仅会对被测结构造成损伤，而且还存在泛化性能不好，易受温度变化影响，工作量较大等问题，使得其实用性和实时性受到限制。针对实际冲击监测定位的要求，需要研究无需大量先验知识，能够适用于采样频率较低的光纤光栅解调仪，且监测原理简洁的新方法。为此，本专利技术提出了一种基于Teager能量算子与样本熵的冲击位置监测方法。
技术实现思路
本专利技术所要解决的技术问题是提供一种用于工程板结构的光纤冲击定位方法，该方法采用分布式光纤布拉格光栅传感器网络感知结构中不同位置的冲击响应信号，首先用Teager能量算子对响应信号进行处理，再提取冲击响应信号的样本熵。利用该样本熵与冲击距离存在的关系，实现冲击定位。该方法通过Teager能量算子提高冲击响应信号质量，并通过样本熵定量估计冲击响应信号的复杂度，可应用于航空航天领域四边固支板结构的冲击定位，无需大量先验知识并可以在线监测辨识冲击位置。本专利技术的技术方案是：基于Teager能量算子和样本熵的板结构冲击监测方法，其采用分布式光纤布拉格光栅传感器网络感知结构中不同位置的冲击响应信号，首先用Teager能量算子对响应信号进行处理，再提取冲击响应信号的样本熵，通过样本熵定量估计冲击响应信号的复杂度，利用该样本熵与冲击距离存在的关系，实现冲击定位。所述的基于Teager能量算子和样本熵的板结构冲击监测方法，其特征在于：包括以下步骤：步骤一：分布式光纤布拉格光栅传感网络布置；步骤二：待测冲击的光纤布拉格光栅传感器响应信号采集及Teager能量算子处理；步骤三：经处理冲击响应信号的样本熵计算；步骤四、确定待测点所在区域，并计算区域内样本冲击点响应信号的样本熵；步骤五：冲击点到传感器的距离与相应样本熵差值之间的关系模型构建；步骤六：基于三角区域法，计算待测冲击点的具体坐标位置。步骤一、分布式光纤布拉格光栅传感网络布置；在四边固支板结构中心位置选取一个正方形监测区域ABCD，其中点A位于监测区域左上角顶点，A、B、C、D按逆时针方向分布于正方形监测区域各顶点；选取板结构待监测区域的中心O作为坐标原点，建立一个二维直角坐标系，定义X轴平行于AD方向，Y轴平行于AB方向；在监测区域顶点位置A、C布置轴向平行于正方形监测区域对角线BD的光纤布拉格光栅传感器，分别计作FBG1、FBG3，在监测区域顶点位置B、D布置轴向平行于正方形监测区域对角线AC的光纤布拉格光栅传感器，分别计作FBG2、FBG4，同时在正方形待监测区域中心位置分别布置垂直对角线AC、BD方向的光纤布拉格光栅传感器，分别计作FBG5、FBG6；将这些光纤布拉格光栅传感器连成一串并粘贴于板结构背面以此构成分布式传感器网络。步骤二：待测冲击的光纤布拉格光栅传感器响应信号采集及Teager能量算子处理步骤2-1、采用冲击锤对板结构的待监测冲击区域内的点施加低速冲击载荷，利用光纤光栅解调仪记录该冲击下分布式传感网络中6个光纤布拉格光栅传感器的冲击响应信号x(n)；步骤2-2、Teager能量算子定义为：Ψ[x(n)]＝[x(n)]2-x(n-1)x(n+1)(1)该算子具有抑制信号干扰和提高时间分辨率的效果，将传感器采集到的响应信号x(n)先去除直流分量，再通过Teager能量算子计算得到经处理的信号u(n)。步骤三：经处理冲击响应信号的样本熵计算；步骤3-1、设经处理的冲击响应信号的时间序列为{u(i),i＝1,2,...,N}，可以按连续顺序组成一组m维矢量X(1),X(2),...,X(N-m)，其中X(i)＝{u(i),u(i+1),...,u(i+m-1)},i＝1,2,...,N-m+1；(2)步骤3-2、定义冲击响应信号的m维相空间中矢量X(i)与量X(j)之间的距离d[X(i),X(j)]：d[X(i),X(j)]＝max[u(i+k)-u(j+k)],k＝0,1,...,m-1；(3)步骤3-3、给定一个阈值r，统计d[X(i),X(j)]小于r的数目称之为模板匹配数，并求出该数目与距离总数N-m-1的比为：求其对所有i的平均值为：步骤3-4、维数m增加1，重复步骤3-1到步骤3-4，得到Bm+1(r)；步骤3-5、则样本熵的为：SampEn(m，r，N)＝-ln[Bm+1(r)/Bm(r)](6)按此步骤计算出光纤布拉格光栅传感器FBG1、FBG2、...、FBG6对应的样本熵SampEni，i＝1,2,3,4,5,6为传感器编号。步骤四、确定待测点所在区域，并计算该区域内样本冲击点经Teager能量算子处理响应信号的样本熵；4-1将正方形待监测区域通过位于顶点的4个传感器和中心的2个传感器沿监测区域对角线划分为4个三角形子监测区域，其中定义由FBG1、FBG4、FBG5和FBG6构成的监测区域为区域Ⅰ，由FBG4、FBG3、FBG5和FBG6构成的监测区域记为区域Ⅱ，由FBG3、FBG2、FBG5和FBG6构成的监测区域记为区域Ⅲ，由FBG2、FBG1、FBG5和FBG6构成的监测区域记为区域Ⅳ；4-2由步骤二及步骤三计算得四个顶点处传感器响应信号的样本熵值SampEn1，SampEn2，SampEn3，SampEn4，将位于正方形对角的传感器样本熵值进行两两比较，其中样本熵更大的两个传感器所在的监测区域确定为待测点所在的三角形子监测区域，从而完成待测点的区域定位；4-3在预测待测点所在区域的各边上均匀施加若干与待测冲击点能量相同的样本冲击载荷，并分别记录每次冲击下冲击点所在边两端光纤布拉格光栅传感器的响应信号。对于整个正方形监测区域中心位置的两个传感器FBG5和FBG6所采集冲击响应信号的取舍，选择其轴向与样本冲击点所在边方向垂直的那个传感器的响应信号；按照步骤二所述方法计算各样本点冲击下，冲击点所在边两端光纤布拉格光栅传感器经Teager能量算子处理响应信号的样本熵和i，k＝1,2,3,4,5,6,i≠k为传感器编号；j＝1,2,3,...,n为第n个样本冲击点本文档来自技高网...

【技术保护点】
基于Teager能量算子和样本熵的板结构冲击监测方法，其特征在于：采用分布式光纤布拉格光栅传感器网络感知结构中不同位置的冲击响应信号，首先用Teager能量算子对响应信号进行处理，再提取冲击响应信号的样本熵，通过样本熵定量估计冲击响应信号的复杂度，利用该样本熵与冲击距离存在的关系，实现冲击定位。

【技术特征摘要】
1.基于Teager能量算子和样本熵的板结构冲击监测方法，其特征在于：采用分布式光纤布拉格光栅传感器网络感知结构中不同位置的冲击响应信号，首先用Teager能量算子对响应信号进行处理，再提取冲击响应信号的样本熵，通过样本熵定量估计冲击响应信号的复杂度，利用该样本熵与冲击距离存在的关系，实现冲击定位。2.根据权利要求1所述的基于Teager能量算子和样本熵的板结构冲击监测方法，其特征在于：包括以下步骤：步骤一：分布式光纤布拉格光栅传感网络布置；步骤二：待测冲击的光纤布拉格光栅传感器响应信号采集及Teager能量算子处理；步骤三：经处理冲击响应信号的样本熵计算；步骤四、确定待测点所在区域，并计算区域内样本冲击点响应信号的样本熵；步骤五：冲击点到传感器的距离与相应样本熵差值之间的关系模型构建；步骤六：基于三角区域法，计算待测冲击点的具体坐标位置。3.根据权利要求2所述的基于Teager能量算子和样本熵的板结构冲击监测方法，其特征在于：步骤一中分布式光纤布拉格光栅传感网络布置的具体过程如下：在四边固支板结构中心位置选取一个正方形监测区域ABCD，其中点A位于监测区域左上角顶点，A、B、C、D按逆时针方向分布于正方形监测区域各顶点；选取板结构待监测区域的中心O作为坐标原点，建立一个二维直角坐标系，定义X轴平行于AD方向，Y轴平行于AB方向；在监测区域顶点位置A、C布置轴向平行于正方形监测区域对角线BD的光纤布拉格光栅传感器，分别计作FBG1、FBG3，在监测区域顶点位置B、D布置轴向平行于正方形监测区域对角线AC的光纤布拉格光栅传感器，分别计作FBG2、FBG4，同时在正方形待监测区域中心位置分别布置垂直对角线AC、BD方向的光纤布拉格光栅传感器，分别计作FBG5、FBG6；将这些光纤布拉格光栅传感器连成一串并粘贴于板结构背面以此构成分布式传感器网络。4.根据权利要求2所述的基于Teager能量算子和样本熵的板结构冲击监测方法，其特征在于：步骤二：待测冲击的光纤布拉格光栅传感器响应信号采集及Teager能量算子处理的具体过程如下：步骤2-1、采用冲击锤对板结构的待监测冲击区域内的点施加低速冲击载荷，利用光纤光栅解调仪记录该冲击下分布式传感网络中6个光纤布拉格光栅传感器的冲击响应信号x(n)；步骤2-2、Teager能量算子定义为：Ψ[x(n)]＝[x(n)]2-x(n-1)x(n+1)(1)该算子具有抑制信号干扰和提高时间分辨率的效果，将传感器采集到的响应信号x(n)先去除直流分量，再通过Teager能量算子计算得到经处理的信号u(n)。5.根据权利要求2所述的基于Teager能量算子和样本熵的板结构冲击监测方法，其特征在于：步骤三：经处理冲击响应信号的样本熵计算具体过程如下：步骤3-1、设经处理的冲击响应信号的时间序列为{u(i)，i＝1，2，...，N}，可以按连续顺序组成一组m维矢量X(1)，X(2)，...，X(N-m)，其中X(i)＝{u(i)，u(i+1)，...，u(i+m-1)}，i＝1，2，...，N-m+1；(2)步骤3-2、定义冲击响应信号的m维相空间中矢量X(i)与量X(j)之间的距离d[X(i)，X(j)]：d[X(i)，X(j)]＝max[|u(i+k)-u(j+k)|]，k＝0，1，...，m-1；(3)步骤3-3、给定一个阈值r，统计d[X(i)，X(j)]小于r的数目称之为模板匹配数，并求出该数目与距离总数N-m-1的比为：求其对所有i的平均值为：步骤3-4、维数m增加1，重复步骤3-1到步骤3-4，得到Bm+1(r)；步骤3-5、则样本熵的为：SampEn(m，r，N)＝-ln[Bm+1(r)/Bm(r)](6)按此步骤计算出光纤布拉格光栅传感器FBG1、FBG2、…、FBG6对应的样本熵SampEni，i＝1，2，3，4，5，6为传感器编号。6.根据权利要求2所述的基于Teager能量算子和样本熵的板结构冲击监测方法，其特征在于：步...

【专利技术属性】
技术研发人员：李彤韡，
申请(专利权)人：中国航空工业集团公司金城南京机电液压工程研究中心，
类型：发明
国别省市：江苏,32