一种时空窗口下的改进型PCA损伤检测方法
技术领域
[0001] 本
发明涉及结构损伤识别的技术领域,尤其是指一种时空窗口下的改进型PCA损伤检测方法。
背景技术
[0002] 结构健康监测的目标是通过实时获取结构的响应信息进而正确评估结构的安全状态。在长期的环境因素、交通荷载和
自然灾害等因素的作用下,结构可能会产生不同程度的损伤,从而引起结构
力学响应发生变化,如固有
频率、应变和位移等。然而,当损伤程度小或者受损
位置远离
传感器时,直接通过监测系统获取的响应信息有时并没有发生显著变
化,因此不能直接利用原始数据正确诊断结构的
健康状态。在这种情况下,需要对原始数据进行深度的学习和训练或者挖掘数据的时
频域特征,从而识别结构的损伤。目前,通过对长期健康监测数据进行
数据处理从而识别损伤的方法有基于神经网络与遗传
算法的方法,小
波分析,ARMA,PCA和MPCA等方法。神经网络与
遗传算法结合的方法需要对监测系统的数据进行
机器学习训练,是一种非线性的损伤检测方法。在实际服役环境中,结构受到噪声、
风、重车荷载等多种随机荷载同时作用,其响应信息是多种随机因素共同作用产生的结果。而
在该方法分析过程中往往不会同时考虑多种荷载作用,因此损伤识别有一定局限性,且鲁
棒性较差。ARMA方法是对时序
信号建立自回归平均滑动模型识别异常信号的数据处理手
法。该方法所研究的时间序列必须是满足同一函数分布的平稳信号,非平稳信号需要通过
差分变成平稳信号再建模,对信号要求较高。PCA(主元素分析)方法由于具备良好的去噪能力而被广泛应用,特点是可以在保证数据信息丢失最少的原则下,将多维
数据压缩到少数
几维。但该方法在实际损伤识别中会因为数据量随时间的增加而导致计算效率下降,且在
降维的过程中会将微小损伤当做噪声处理。随后出现的MPCA方法,是在
选定长度的移动的
时间窗口上对结构所有测点的数据统一使用PCA,求出该时间窗内的
特征向量,并通过对比特征向量在时域上的变化以识别损伤。但是,该方法将所有测点的数据统一使用PCA方法进行分析,忽略了远离损伤区域的测点的响应对损伤不敏感这一特点。因此,在数据分析过程中延迟了发现损伤的时间,降低了测损
精度。
发明内容
[0003] 本发明的目的在于克服
现有技术的缺点与不足,提出了一种时空窗口下的改进型PCA损伤检测方法,该方法是在限定的时间和空间窗口内对结构响应进行PCA分析,从而提
高了损伤识别的精度和灵敏度,可利用长期结构监测数据对大跨度实体
桥梁或其它结构进
行损伤识别。
[0004] 为实现上述目的,本发明所提供的技术方案为:一种时空窗口下的改进型PCA损伤检测方法,首先,利用数值仿真或者实验手段,在损伤处于不同空间位置的多种工况下,确定不同位置损伤所对应的空间影响范围,进而确定对损伤敏感的结构空间区域,即空间窗
口;其次,依据PCA原理确定与外荷载周期有关的时间窗口;最后在确定时间窗口和空间窗口后,对空间窗口内所有测点数据使用PCA方法,以特征向量作为监测指标,实时监测损伤并实现损伤
定位。其包括以下步骤:
[0005] 1)通过无损和不同位置损伤多种工况下的响应信息,确定损伤的空间影响范围。根据空间影响范围,进一步在空间上对结构进行区域划分,实现对测点的分组,从而确定空间窗口;
[0006] 2)基于外荷载周期确定PCA分析所采用的时间窗口;
[0007] 3)在限定的空间窗口和时间窗口内,对结构空间窗口范围内测点所采集的响应数据,使用PCA分析,得到对应的特征向量;
[0008] 4)移动时间窗口,计算不同时间窗口下的特征向量,对比分析特征向量随时间的变化,确定有无损伤及损伤出现的时间;
[0009] 5)计算不同空间窗口所对应的特征向量,对比分析不同空间窗口下的特征向量随时间的变化规律,确定损伤出现的位置,从而实现损伤定位。
[0010] 在步骤1)中,通过无损和不同损伤工况下的响应,确定PCA分析中使用的空间窗口;其中,在同一荷载条件下获取不同损伤工况及结构无损条件下测点的响应数据,分析中所需的样本数据的时间总尺度至少为2个荷载周期;Udi(t)表示有损情况下结构的i号单元
在t时刻的响应,Ui(t)代表无损情况下i号单元在t时刻的响应,当|(Udi(t)-Ui(t))/Ui(t)|>η,判定单元i在t时刻为受到损伤影响的单元,其中η是与传感器的精度有关的常数;假设N为某一损伤工况下一个荷载周期内的测点数据的总
采样数,ni为i号单元在一个荷载周期
内受损伤影响的累积次数,依据ni/N>α,判定该单元为该工况下的受损伤影响单元,其中α为依据结构特点所确定的某一常数,且0<α<1,以此类推,能确定在该损伤工况下结构中所有受损伤影响的单元与非受损伤影响的单元;同理,能够得到每一种损伤工况下,由受损伤影响的单元组成的影响区及对应的非影响区,将影响区内的所有测点组成的集合作为最
优空间窗口。
[0011] 在步骤2)中,所选的时间窗口应为外荷载周期的整数倍。所选的时间窗口与外荷载周期相等时,能提高求解特征向量的计算速度。
[0012] 在步骤3)中所述的在限定时间和空间窗口内使用的PCA数据处理方法,包括以下步骤:
[0013] 3.1)确定限定时空窗口内的响应矩阵
[0014] 假设选定的空间窗口内有m个测点,编号依次为1,2,…m,选定的时间窗口内,每个测点采集的数据个数为Nn,则有
[0015]
[0016] 其中,t为对应的监测时间,时间窗口区间为[t1,tN],ui(tj)为i号传感器在tj时刻的响应,每一列代表一个传感器在选定的时间窗口内测得的所有响应数据;
[0017] 3.2)矩阵零均值化
[0018] 对于i号传感器所采集的数据,其均值 依据下式计算:
[0019]
[0020] 将矩阵中的每一列减去该列的平均值,得到零均值化的矩阵为
[0021]
[0022] 3.3)求对应的协方差矩阵
[0023] 协方差矩阵为: 其中m为列数;
[0024] 3.4)求协方差矩阵的特征值和特征向量
[0025] (C-λjI)Ψ=0
[0026] 依据上述公式,求对应的特征向量Ψ,其中I为对应的单位矩阵,λj为特征值,并有λ1>λ2>...>λm,j取1,2,Ψji为对应的特征值λj下的特征向量分量;
[0027] 在步骤4)中,移动时间窗口求特征向量,包括以下步骤:
[0028] 4.1)当时间向前移动k步时,响应矩阵为:
[0029]
[0030] 假设每个传感器的总采样数为Ns,令k的取值依次为1,2,…,NS-Nn,并重复步骤3.2)至步骤3.4),求得第k个时间步下的特征向量分量Ψji(k);
[0031] 4.2)首先,计算结构无损情况下在1个荷载周期内计算所得的特征向量各分量的均值 即
[0032]
[0033] 依据以下判定准则进行损伤识别:
[0034]
[0035] 其中,σji为特征向量分量Ψji的标准差。当上式成立时,表示结构不存在损伤;当上式不成立时,结构有损,此时依据k并得到对应的损伤出现时刻。
[0036] 在步骤5)中,移动空间窗口求特征向量,包括以下步骤:
[0037] 5.1)根据所得到的最优空间窗口,对结构从空间上进行区域划分和测点分组,计算每一个空间窗口内测点数据所对应的特征向量。首先将每个窗口内测点采集的响应
整理成矩阵,然后重复步骤4)求Ψji(k);
[0038] 5.2)依据以下判定准则判定损伤发生的位置:
[0039]
[0040] 其中,β为依据结构特点所确定的某一常数。当上式成立时,则判定损伤发生在该空间窗口的i号传感器附近。
[0041] 本发明与现有技术相比,具有如下优点与有益效果:
[0042] 1、较现有技术,能识别更小程度的损伤,测损精度更高。
[0043] 2、测损灵敏度方面,较其它的基于监测时序数据的损伤识别方法,更快确定有无损伤及损伤出现的时间。
[0044] 3、本发明方法简单高效,克服了不同程度的噪声干扰,较前人的方法具有较强的鲁棒性,不受环境噪声的干扰。
[0045] 4、所需数据来自与桥梁健康监测系统,不需增加其它检测设备成本,仅需要成桥时桥梁健康状态的数据及现状态下的监测数据。
[0046] 5、本发明方法对于桥梁结构没有限制,基于大量实测或模拟数据,从数据本身即包含桥梁的内在特征,因此可适用于大跨径和中小跨径桥梁。
附图说明
[0047] 图1为本发明逻辑流程示意图。
[0048] 图2为本发明使用的模型传感器位置及损伤工况说明图。
[0049] 图3为本发明检测损伤的灵敏度示意图。
具体实施方式
[0050] 下面结合具体
实施例对本发明作进一步说明。
[0051] 如图1所示,本实例所提供的时空窗口下的改进型PCA损伤检测方法,包括以下步骤:
[0052] 1)确定PCA分析所采用的空间窗口
[0053] 1.1)确定结构的受到损伤影响区域。本实例中的模型为简支梁模型。模型尺寸长宽高为:4m*0.2m*0.4m。长度方向划分50个网格,高度方向划分10个网格。单元尺寸为:
0.04m*0.04m。模型及传感器位置如图2所示。
温度荷载中,
底板Tb服从Tb=20+10sin(πn/
730),n=1,2,...,5840。一天内采集6:00,12:00,18:00,24:00时刻温度,总采样数5840对
应采样时长为四年。顶板温度Tt=Tb+10,内部温度场沿梁高方向呈线性分布。采样时长内,前两年为无损,后两年有损。损伤位置如图2所示。损伤程度分为模量减少20%,40%,60%,
80%。损伤工况如下表所示。
[0054]
[0055] εdi(t)表示i号单元t时刻应变,εi(t)代表无损情况下i号单元t时刻应变。|(εdi(t)-εi(t))/εi(t)|>5%,判定单元i为在t时刻受到损伤影响的单元。N为某一损伤工况下一个周期内的总采样数,ni为i号单元一个荷载周期内受损伤影响的累计采样次数,本实例中N=1460。依据ni/N>60%,判定该单元为该工况下的受到损伤影响的单元。以此类推,能确定在该损伤工况下结构中所有受到损伤影响的单元与非受损伤影响的单元,及受损伤影
响的单元组成的影响区及对应的非影响区。同理,可以得到16种工况下,由受损伤影响的单元组成的影响区及对应的非影响区。
[0056] 1.2)在16种工况中发现,损伤影响范围在自损伤位置起上下左右四个传感器范围内(损伤程度在40%以上时,受损区会超出该位置),以及靠近底板的部位。以损伤工况1为例。实例中,依据损伤影响范围对简支梁从空间上分域,再对测点进行分组。分组后为:(1)[sensor1,sensor2,sensor6,sensor7];(2)[sensor2,sensor 3,sensor 7,sensor 8];(3)[sensor 3,sensor 4,sensor 8,sensor 9];(4)[sensor 4,sensor 5,sensor 9,sensor
10];(5)[sensor 1,sensor 2,sensor 3,sensor 4,sensor 5]。以便于后文阐述本发明中特征向量的计算和与传统的MPCA方法对比,仅例举具有代表性的空间窗口所对应的特征向
量的计算思路,其他窗口计算步骤与之同理。空间窗口A)[sn1,sn2,sn6,sn7],窗口B)[sn1,sn2,sn3,sn4,sn5],传统的MPCA使用的窗口C为:[sn1,sn2,sn3,sn4,sn5,sn6,sn7,sn8,
sn9,sn10],对照组窗口D为:[sn1,sn2,sn3,sn6,sn7,sn8]。
[0057] 2)确定PCA分析所采用的时间窗口。时间窗口与温度荷载的周期相等。本实例中温度荷载周期为1年,时间窗口亦为1年,对应1460个采样数。
[0058] 3)计算限定的时间窗口和空间窗口内所对应的特征向量
[0059] 3.1)将空间窗口内所有测点,在限定的时间窗下所采集的应变整理成矩阵形式,如窗口A整理后的矩阵为1460*4的应变矩阵,窗口B为1460*5的应变矩阵。
[0060] 或
[0061]
[0062] 3.2)矩阵零均值化
[0063] 对于i号传感器所采集的数据,其均值 有:
[0064]
[0065] 将矩阵中的每一列减去该列的平均值,得到零均值化的矩阵。如窗口A整理后形式如下:
[0066]
[0067] 3.3)求对应的协方差矩阵
[0068]
[0069] 协方差矩阵为: m为列数。在此测点组合中,m=4。
[0070] 3.4)求协方差矩阵的特征值和特征向量
[0071] (C-λjI)Ψ=0
[0072] 依据上述公式,求对应的特征向量Ψ。其中I为对应的单位矩阵,λ为特征值,并有λ1>λ2>...>λm,j取1。Ψ11,Ψ21,Ψ31,Ψ41为对应的特征值λ1下,1号、2号、3号、4号传感器对应的特征向量分量。
[0073] 4)移动时间窗口求不同时间窗口所对应特征向量,确定有无损伤及损伤出现的时间。
[0074] 4.1)当时间向前移动k步时,响应矩阵为:
[0075] 或
[0076] 每个传感器的总采样数为5840,令k的取值依次为1,2,…,4380。重复步骤3),求第k个时间步下的Ψ11(k),Ψ21(k),Ψ31(k),Ψ41(k)。
[0077] 4.2)首先,计算简支梁无损情况下在1个荷载周期内计算所得的特征向量各分量的均值 即
[0078]
[0079] 依据以下判定准则进行损伤识别:
[0080]
[0081] 其中,σji为特征向量分量Ψji的标准差。当上式成立时,表示结构不存在损伤;当上式不成立时,结构有损,此时依据k并得到对应的损伤出现时刻。
[0082] 在本实例中,计算可得,窗口A识别到损伤时,k为1735。窗口B识别到损伤时,k为1519。使用传统的MPCA算法,窗口C识别到损伤时,k为1773。对照组窗口D,k为1768。本发明能最快识别损伤,如图3所示。
[0083] 5)确定损伤出现的位置,从而实现损伤定位。
[0084] 5.1)计算测点分组后的各个空间窗口所对应的特征向量Ψji(k)。窗口1)[ε1,ε2,ε6,ε7],窗口2)[ε2,ε3,ε7,ε8],窗口3)[ε3,ε4,ε8,ε9],窗口4)[ε4,ε5,ε9,ε10],窗口5)[ε1,ε2,ε3,ε4,ε5]重复步骤3和步骤4求Ψji(k)。
[0085] 5.2)依据以下判定准则判定损伤发生的位置:
[0086]
[0087] 当上式成立时,则判定损伤发生在该空间窗口的i号传感器附近。
[0088] 在本实例中,计算可得, 能有效地识别损伤处于1号传感器附近,如图3所示。
[0089] 以上所述实施例只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。