基于分形理论的板结构冲击监测方法
【技术领域】
[0001] 本发明属于结构健康监测的冲击监测技术领域，尤其涉及基于分形理论的板结构 冲击监测方法。
【背景技术】
[0002] 现代大型结构如航天飞机、高层结构、新型桥梁、大跨度网架结构等在复杂的服役 环境中将受到设计载荷的作用以及各类突发性外在因素如外来物冲击、振动等影响，而这 些振动或冲击会导致结构不同程度的损伤。这些损伤的隐蔽性强，损伤、失效方面机理复 杂、损伤类型和程度难以判断，但这些损伤如果不及时发现和采取相应的措施，会造成结构 损伤积累，使得其结构力学性能退化，威胁结构安全，造成灾难性的后果。
[0003] 针对冲击的定位识别问题的研宄主要集中于两个方向：一是基于冲击在板结构 中激励出的Lamb波进行研宄，具有代表性的有三角形定位、四边形定位和基于小波变换的 时差定位等，取得了较好的冲击位置识别结果，但由于板结构的各向异性，使得其中传播的 Lamb波的传播波速随着频率和传播方向的不同而变化，直接降低了基于Lamb波定位的精 度。另外这些方法由于其数据采集与传输过程易受到电磁干扰的影响，且其压电监测阵列 需要大量传输电缆线，大大加剧了系统的复杂程度；二是基于低速冲击引起板结构应力变 化进行研宄，具有代表性的如基于HausdorfT距离位置识别和基于支持向量机的位置识别 等，但这些算法受其自身特点的影响很大，如采用支持向量机算法的定位精度取决于核宽 和惩罚因子的选择，但核宽与惩罚因子的选择没有相关理论予以支持，只能依靠实验试凑 或其他方式确定，若参数选择不合适，会降低定位精度。此外，还有基于大量先验知识的神 经网络冲击位置辨识方法，此类方法缺点是需要大量不同的冲击样本来训练网络，泛化性 能不好，易受到温度变化影响，工作量较大，使得其实用性和实时性受到限制。
[0004] 因此，针对目前的冲击监测定位方法的不足，需要研宄无需大量先验知识，能够适 用于采样频率较低的常规光纤光栅解调仪，且监测系统简洁的新方法。为此，本发明提出了 基于分形理论的板结构冲击监测方法。
【发明内容】
[0005] 本发明所要解决的技术问题是提供一种用于工程板结构的光纤冲击定位方法，该 方法采用分布式光纤布拉格光栅传感器网络感知结构中不同位置的冲击响应频谱信号，提 取出冲击信号的特征信息，利用传感器信号频谱分形维数与冲击距离以及位置存在关联， 以此为特征参数实现对冲击位置定位，适用于四边固支板结构，该方法可以简单快捷、实 时、在线监测辨识冲击位置。
[0006] 本发明所采用的技术方案包括以下步骤：
[0007] -种基于分形理论的板结构冲击监测方法，其特征在于，包括以下步骤：
[0008] 步骤一、分布式光纤布拉格光栅传感器布局
[0009] 在四边固支板结构中心部位构建一个正方形监测区域AB⑶，其中点A、B、C、D为逆 时针方向排序的正方形各顶点；建立一个二维直角坐标系，选取板结构待监测区域的左下 角顶点C作为坐标原点，定义X轴平行于AB方向，Y轴平行于BC方向；在板结构正方形监测 区域A、B、C、D四个顶角位置分别布置平行于X轴方向的四个光纤布拉格光栅传感器FBG1、 FBG2、FBG3、FBG4,同时在A、B、C、D四个顶角位置分别布置平行于Y轴方向的四个光纤布拉 格光栅传感器FBG5、FBG6、FBG7、FBG8 ;将这些光纤布拉格光栅传感器粘贴于试件结构的背 面，采用光纤跳线将FBGl和FBG2, FBG3和FBG4, FBG5和FBG8, FBG6和FBG7分别串行连接， 以此构成分布式传感器网络，八个光纤布拉格光栅传感器构成的正方形所覆盖的区域即为 板结构试件的冲击待监测区域；
[0010] 步骤二、光纤布拉格光栅传感器冲击响应信号计盒分形维数计算方法
[0011] 步骤2-1、对不同编号光纤布拉格光栅传感器冲击响应信号利用离散余弦变换分 别计算其频谱，设其中任意一个光纤布拉格光栅传感器冲击响应信号频谱序列表示为=X 1, X2，…，Xi,…，Xn;这里X $频谱中各频率分量，i = 1，2，…，N ;再按以下步骤得到每个 光纤布拉格光栅传感器冲击响应信号频谱的计盒分形维数；
[0012] 步骤2-2、取盒子尺度为r，l彡r << N，将频谱序列划分为N/r个区间段；X1, X2，…，Xr，X r+1，Xr+2···，X2r，…，XN_r+1，X N_r+2，…，Xn;取 N/r 为整数；
[0013] 在每一区间段内，寻找相对离差最大的两样本点XpXr并计算第U段区间最大的 相对离差d u= max(X厂乂』），其中u = 1，2，…，N/r，表示为分段区间序号；
[0014] 计算覆盖第U段区间的盒子数mu:mu= [du/r]+l其中，[.]表示取整运算；则覆 盖整条频谱曲线的盒子数为： Nfr
[0015] = u~\
[0016] 步骤2-3、将尺度r放大为η · r，η = 2, 3，…，ρ，且ρ · r < Ν，重复步骤（2-2)得 到另外P-I个覆盖频谱序列曲线的盒子数N(nr);
[0017] 步骤 2-4、令 X = log(l/nr)，y = log(N(nr))，η = 1，2,3，…，p，其中 p 为迭代 次数；按上述坐标采用描点法描出P个点；在无标度区间对所描出的点使用最小二乘法得 到拟合直线y = kx+h，h为所拟合直线的常数项，此式中斜率k即为计盒分形维数e ;
[0018] 步骤三、样本冲击点对应的光纤布拉格光栅传感器响应信号采集与计盒分形维数 计算
[0019] 作平行于X轴且关于AD边对称的中心线α，作平行于y轴且关于AB边对称的中 心线β，中心线α与⑶边的垂直距离和中心线β与BC边的垂直距离均为L ;在待监测区 域选取a点为样本点，样本点a到中心线α和中心线β的垂直距离均为1 ;采用冲击锤对 样本点a进行冲击，得到八个光纤布拉格光栅传感器对应于样本点a的冲击响应信号；按照 步骤二计算样本点a冲击时八个光纤布拉格光栅传感器信号频谱的计盒分形维数，分别记 为e ai，i =传感器编号1，2，…，8 ;
[0020] 步骤四、待测冲击点对应的光纤布拉格光栅传感器响应信号采集与计盒分形维数 计算
[0021] 采用冲击锤对板结构试件的待监测冲击区域ABCD中的任意点作为待测点进行冲 击加载，记录该冲击下分布式传感网络8个光纤布拉格光栅传感器的冲击响应信号。重复 步骤二，得到相应冲击下的光纤布拉格光栅传感器FBGl、FBG2、"'FBGS对应的计盒分形维 数，分别记为epe；^、…、e8;
[0022] 步骤五、待测冲击点的X轴坐标计算
【主权项】
1. 一种基于分形理论的板结构冲击监测方法，其特征在于，包括以下步骤： 步骤一、分布式光纤布拉格光栅传感器布局 在四边固支板结构中心部位构建一个正方形监测区域AB⑶，其中点A、B、C、D为逆时 针方向排序的正方形各顶点；建立一个二维直角坐标系，选取板结构待监测区域的左下角 顶点C作为坐标原点，定义X轴平行于AB方向，Y轴平行于BC方向；在板结构正方形监测 区域A、B、C、D四个顶角位置分别布置平行于X轴方向的四个光纤布拉格光栅传感器FBG1、 FBG2、FBG3、FBG4,同时在A、B、C、D四个顶角位置分别布置平行于Y轴方向的四个光纤布拉 格光栅传感器FBG5、FBG6、FBG7、FBG8 ;将这些光纤布拉格光栅传感器粘贴于试件结构的背 面，采用光纤跳线将FBG1和FBG2,FBG3和FBG4,FBG5和FBG8,FBG6和FBG7分别串行连接， 以此构成分布式传感器网络，八个光纤布拉格光栅传感器构成的正方形所覆盖的区域即为 板结构试件的冲击待监测区域； 步骤二、光纤布拉格光栅传感器冲击响应信号计盒分形维数计算方法 步骤2-1、对不同编号光纤布拉格光栅传感器冲击响应信号利用离散余弦变换分别 计算其频谱，设其中任意一个光纤布拉格光栅传感器冲击响应信号频谱序列表示为：&， X2，…，Xi，…，XN;这里X#频谱中各频率分量，i= 1，2，…，N;再按以下步骤得到每个 光纤布拉格光栅传感器冲击响应信号频谱的计盒分形维数； 步骤2-2、取盒子尺度为r，1彡r<<N，将频谱序列划分为N/r个区间段；X2，…，Xr，Xr+1，Xr+2...，X2r，…，XN_r+1，XN_r+2，…，XN^N/r为整数； 在每一区间段内，寻找相对离差最大的两样本点X」，并计算第u段区间最大的相对 离差du=max(X厂乂』），其中u= 1，2，…，N/r，表示为分段区间序号； 计算覆盖第u段区间的盒子数mu:mu=[du/r]+l其中，[.]表示取整运算；则覆盖整 条频谱曲线的盒子数为：
步骤2-3、将尺度r放大为n?r，n= 2, 3，…，p，且p?r彡N，重复步骤（2-2)得到另 外P-1个覆盖频谱序列曲线的盒子数N(nr); 步骤 2-4、令x=log(l/nr)，y=log(N(nr))，n= 1，2,3，…，p，其中p为迭代次数； 按上述坐标采用描点法描出P个点；在无标度区间对所描出的点使用最小二乘法得到拟合 直线y=kx+h，h为所拟合直线的常数项，此式中斜率k即为计盒分形维数e; 步骤三、样本冲击点对应的光纤布拉格光栅传感器响应信号采集与计盒分形维数计算 作平行于X轴且关于AD边对称的中心线a，作平行于y轴且关于AB边对称的中心线 0，中心线a与⑶边的垂直距离和中心线0与BC边的垂直距离均为L;在待监测区域选 取a点为样本点，样本点a到中心线a和中心线0的垂直距离均为1 ;采用冲击锤对样本 点a进行冲击，得到八个光纤布拉格光栅传感器对应于样本点a的冲击响应信号；按照步 骤二计算样本点a冲击时八个光纤布拉格光栅传感器信号频谱的计盒分形维数，分别记为 eai，i=传感器编号1，2，…，8 ; 步骤四、待测冲击点对应的光纤布拉格光栅传感器响应信号采集与计盒分形维数计算 采用冲击锤对板结构试件的待监测冲击区域ABCD中的任意点作为待测点进行冲击加 载，记录该冲击下分布式传感网络8个光纤布拉格光栅传感器的冲击响应信号。重复步骤 二，得到相应冲击下的光纤布拉格光栅传感器FBG1、FBG2、…、FBG8对应的计盒分形维数， 分别记Sei、e2、…、e8; 步骤五、待测冲击点的x轴坐标计算
(1) 式中em/en(m、n分别表示传感器编号：1，2,3,4);
;，m、n分别表不传感器编号：1，2,3,4; umn 步骤六、待测冲击点的y轴坐标计算
(2) 式中em/en，m、n分别表不传感器编号：5,6, 7, 8 ;
;=ea ，m、n分别表不传感器编号：5,6,7,8。 um7i
【专利摘要】本发明公开一种基于分形理论的板结构冲击监测方法，属于结构健康监测的冲击监测技术领域。包括以下步骤：步骤一：分布式光纤布拉格光栅传感器布局；步骤二：光纤布拉格光栅传感器冲击响应信号计盒分形维数计算方法；步骤三：样本冲击点对应的光纤布拉格光栅传感器响应信号采集与计盒分形维数计算；步骤四、待测冲击点对应的光纤布拉格光栅传感器响应信号采集与计盒分形维数计算；步骤五：待测冲击点的x轴坐标计算；步骤六：待测冲击点的y轴坐标计算。本发明的定位算法通过计算光纤布拉格光栅传感器冲击响应信号频谱的计盒分形维数来实现冲击点位置的判定，具有无需大量先验知识、简单快速、方法便捷，实用性强等特点。
【IPC分类】G01M7-08, G01B11-00
【公开号】CN104776966
【申请号】CN201510151603
【发明人】曾捷, 熊稚莉, 王博, 潘晓文, 刘苏州, 毛超群, 张先辉
【申请人】南京航空航天大学
【公开日】2015年7月15日
【申请日】2015年4月1日