[发明专利]基于EMD的结构刚度损伤监测方法及系统

摘要

本发明公开了结构刚度损伤的监测方法及系统，方法包括步骤监测结构振动响应并获取结构加速度响应信号；对加速度响应信号进行EMD分解，确定各个时刻第一个IMF分量斜率，并计算结构不同位置的监测因子；利用约束条件剔除虚假监测因子，监测因子与结构刚度损伤程度呈线性关系；根据监测因子随时间的变化特征确定结构刚度损伤发生的时刻，比较结构不同位置的监测因子的分布确定损伤发生的位置；根据刚度损伤发生时刻对应的监测因子的幅值确定结构的刚度损伤程度。本发明特别适用于冲击荷载作用下的结构微小刚度损伤监测判定，能准确的识别微小程度的刚度损伤，具有较强的抗噪能力，通过计算监测因子的大小可以判定损伤程度，具有较高的应用价值。

主权项

一种基于EMD的结构刚度损伤监测方法，其特征是包括如下步骤：步骤S1，在结构不同位置上安装多个加速度传感器实时监测所述结构的动力响应，获取结构不同位置的加速度响应信号；步骤S2，对所述结构不同位置的加速度响应信号进行EMD分解，并获取具有最高频信号成分的第一个IMF分量；步骤S3，计算各个时刻IMF分量的斜率，并计算结构不同位置的一种用以判断结构刚度损伤的监测因子；步骤S4，利用约束条件剔除虚假的监测因子，所述监测因子与所述结构的刚度损伤程度呈线性关系；步骤S5，根据所述监测因子随时间的变化特征确定结构刚度损伤发生的时刻，通过比较结构不同位置的监测因子的分布确定损伤发生的位置；步骤S6，根据刚度损伤发生时刻对应的监测因子的幅值确定结构的刚度损伤程度；所述的步骤S2具体为：通过以下公式对所述结构不同位置的加速度响应信号进行EMD分解，并确定各个不同的IMF分量：x··(t)=Σi=1nfEMDi(x··(t))+rn(t);]]>其中：为加速度响应信号；为加速度信号经过EMD分解后的第i个IMF分量；rn(t)为加速度信号经过EMD分解后的残余趋势项；所述的步骤S3计算结构不同位置的监测因子具体通过以下公式进行：MIi＝|(Di‑Di‑1)+(Di‑Di+1)|＝|2Di‑Di‑1‑Di+1|(i＝2,3,...,tmax‑1)；其中，Di表示i时刻的结构加速度响应第一个IMF分量的变化率；Di‑1和Di+1表示i‑1和i+1时刻的加速度响应第一个IMF分量的变化率；tmax为加速度响应信号的最大时间长度；所述的步骤S3计算结构不同位置的监测因子的步骤如下：S3‑1，建立发生和不发生刚度损伤结构的频率之间的相互关系若等效单自由度结构体系由于杆件屈曲或失稳，其结构的刚度在时刻ti发生了损伤和减小，则结构刚度由正常值K0减小了ΔK，变为Ks：K=K0(0≤t≤ti)Ks(ti<t);]]>ΔK＝K0‑Ks；则无刚度损伤结构的频率f0和有刚度损伤结构的频率fs分别表示为：f0=12πK0M;]]>fs=12πKsM;]]>式中，M为结构的质量；结构的刚度变化用无刚度损伤结构和有刚度损伤结构的频率表示：ΔK=K0-Ks=4π2M(f02-fs2);]]>S3‑2，建立无失稳损伤的原始结构的动力响应计算方法没有发生失稳损伤的结构不存在刚度损伤，因此其等效单自由度体系的运动方程表示为：x··+4πξf0x·+4π2f02x=0;]]>式中，ξ为结构体系的阻尼比；结构在外荷载引起的脉冲作用下将发生具有初始速度为v0的振动，则结构的位移、速度和加速度响应分别为：x(t)=v0sin(2πf0ξdt)2πf0ξd·e-2πf0ξt;]]>x·(t)=v0e-2πf0ξt(cos(2πf0ξdt)-ξξdsin(2πf0ξdt));]]>x··(t)=-2πf0v0e-2πf0ξt·[sin(2πf0ξdt)(2ξd2-1)+2ξξdcos(2πf0ξdt)]ξd;]]>式中：ξd=1-ξ2;]]>S3‑3，建立发生失稳及刚度损伤的结构的动力响应计算方法在结构在外荷载作用下发生振动中，假设在ti时刻发生了构件失稳，由于构件失稳过程具有突然性，因此发生过程很短，这将导致构件的刚度在很短的时间内发生减小；为描述该刚度变化过程，采用一个新的时间坐标轴t1＝t‑ti来描述发生失稳的结构的振动状况；因此，有刚度损伤结构的运动方程表示为：x··s+4πξfsx·s+4π2fs2xs=0,(t>ti);]]>则，有损伤结构振动的初始条件由无损结构在时刻ti的位移和速度响应确定：xs(0)=x(ti)=sin(2πf0ξdti)2πf0ξd·v0e-2πf0ξti;]]>x·s(0)x·(ti)=v0e-2πf0ξti(cos(2πf0ξdti)-ξsin(2πf0ξdti)ξd);]]>由此计算得到在t1时刻有损伤结构的加速度响应为：x··s(t1)=-2πv0fs(E1+E2+E3)ξd2f0·e-2πξ(fst1+f0ti);]]>式中：E1＝fs sin(2πf0ξdti)[ξd cos(2πfsξdt1)‑ξsin(2πfsξdt1)]；E2=f0sin(2πfsξdt1)(1-2ξd2)(ξsin(2πf0ξdti)-ξdcos(2πf0ξdti));]]>E3=-2f0ξcos(2πfsξdt1)[ξξdsin(2πf0ξdti)-ξd2cos(2πf0ξdti)];]]>由于工程结构的阻尼比往往很小，因此构件发生失稳破坏而导致的刚度损伤发生的时间很短，则有：t1＝ti+1‑ti＝Δt≈0；sin(2πfsξdt1)≈0；cos(2πfsξdt1)≈1；由此得结构在发生刚度损伤后后ti+1的时刻的加速度响应为：x··(ti+1)=x··s(Δt)=-2πfsv0e-2πf0ξtif0ξd{fssin(2πf0ξdti)-2f0ξ[ξsin(2πf0ξdti)-ξdcos(2πf0ξdti)]};]]>S3‑4，进行结构振动响应的EMD分解为了建立损伤监测因子，则需对结构的加速度响应进行经验模态分解；具体处理方法是：首先，确定加速度响应的多个局部极大值和局部极小值：采用多次样条函数将的局部极大值点与局部极小值点分别拟和得到其峰值的上包络曲线与下包络曲线然后计算两包络线的均值m1(t)：m1(t)=x··(t)envelopeu+x··(t)envelopel2;]]>将原加速度序列减去该平均包络m1(t)后即得一个去掉低频的新加速度时程序列h1(t)：h1(t)=x··(t)-m1(t);]]>对得到的h1(t)重复以上数据过程，重复k次直至所得到的平均包络趋于零为止：h1k(t)＝h1(k‑1)(t)‑m1k(t)；其中：h1k(t)为第k次处理所得加速度数据；h1(k‑1)(t)为第k‑1次处理所得加速度数据；m1k(t)为h1(k‑1)(t)上下包络线的均值；由此得到该加速度响应的第一个内敛模函数分量c1(t)：c1(t)＝h1k(t)；第一个IMF分量c1(t)代表了原始加速度信号中的最高频成分；将原始加速度响应减去第一个IMF分量c1(t)，得去掉高频成分的加速度响应时程r1(t)；将r1(t)再作为要分解的信号重复上述过程，直至所剩余信号r1(t)已是一单调函数时停止此分解过程；此时的参与量rn(t)代表原始加速度响应的低频趋势项；由此确定加速度响应的一组IMF分量c1(t),c2(t)…cn(t)；原始的加速度响应由全部IMF分量和一个趋势项的叠加来表示：x··(t)=Σi=1nci(t)+rn(t);]]>将结构振动信号的EMD分解过程采用一个隐函数来表示，则无损伤结构加速度相应的EMD分解信号表示为：x··(t)=Σi=1nci(t)+rn(t)=Σi=1nfEMDi(x··(t))+rn(t)=Σi=1nfEMDi(-2πf0v0e-2πf0ξt·[sin(2πf0ξdt)(2ξd2-1)+2ξξdcos(2πf0ξdt)]ξd)+rn(t)]]>其中：第i个IMF分量ci(t)表示为：ci(t)=fEMDi(-2πf0v0e-2πf0ξt·[sin(2πf0ξdt)(2ξd2-1)+2ξξdcos(2πf0ξdt)]ξd);]]>同理，发生刚度损伤后结构的加速度响应表示为：x··s(t)=Σi=1ncis(t)+rns(t)=Σi=1nfEMDi(x··(t))+rns(t)=Σi=1nfEMDi(-2πfsv0e-2πf0ξtif0ξd{fssin(2πf0ξdti)-2f0ξ[ξsin(2πf0ξdti)-ξdcos(2πf0ξdti)]})+rns(t)]]>；其中：第i个IMF分量表示为：cis(t)=fEMDi(-2πfsv0e-2πf0ξtif0ξd{fssin(2πf0ξdti)-2f0ξ[ξsin(2πf0ξdti)-ξdcos(2πf0ξdti)]});]]>S3‑5，确定损伤监测因子EMD分解基于加速度响应局部特征时间尺度，从原加速度时程中提取固有模态函数，其本质是将加速度信号中不同频率和尺度的波动或趋势逐级分解开来；所分解出的各IMF分量分别包含了原加速度信号的不同时间尺度和频率特征的局部特征信息；突变损伤时刻的信号不连续具有两个明显的特点：(1)信号的幅值在损伤时刻ti到时刻ti+1发生了很大的跳跃；(2)在时刻ti‑1和时刻ti+1信号的斜率远小于损伤时刻ti的信号斜率；事实上研究表明，发生刚度突然损伤时的结构加速度响应的第一个IMF分量也具有上述的两个相同的特点；显然可知，发生瞬时失稳事故时，结构构件刚度突然减小，结构的加速度响应出现了一个突然的跳跃；这个突然的跳跃信号具有明显的高频特征和大振幅特点；由于突变刚度损伤具有高频特性，因此，其加速度响应的突变信号只保留在具有最高频成分的第一个IMF分量中；损伤前后的结构加速度响应第一个IMF的变化率Di表示为：Di=Δc1Δt=c1,i+1-c1,iΔt=fEMD1(x··i+1-x··i)=·fEMD1(-2πξdv0(fs2-f02)f0·e-2πf0ξti·sin(2πf0ξdti)),(i=1,2,...,tmax-1);]]>式中：Δt为加速度响应信号的时刻间距，tmax为加速度响应信号的最大时间长度；由于：fs2-f02=ΔK4π2M;]]>则损伤前后的结构加速度响应第一个IMF的变化率Di表示为：Di=·fEMD1(-ΔKξdv02πf0M·e-2πf0ξti·sin(2πf0ξdti)),(i=1,2,...,tmax-1);]]>由于EMD分解过程是一个线性过程，原始加速度信号表示为所有IMF分量和残余分量的线性叠加，因此存在如下关系：|Di|∝|ΔK|；前述的加速度响应信号不连续的第二个特点在数学表示为：|Di|>>|Di-j||Di|>>|Di+j|,(j=1,2,...,tmax);]]>由此，得到一种结构由于瞬时失稳事故引起的刚度损伤的监测因子(Monitoring Index)MIi：MIi＝|(Di‑Di‑1)+(Di‑Di+1)|＝|2Di‑Di‑1‑Di+1|(i＝2,3,...,tmax‑1)；由于存在如下关系：|Di|>>|Di-1||Di|>>|Di+1|,(i=2,3,...,tmax-1);]]>因此有：MIi=|2Di-Di-1-Di+1|=·2|Di|;]]>由上述推导过程可知，所采用的基于EMD的监测因子与结构刚度损伤的程度呈正比关系，即：MIi∝|ΔK|；基于该损伤指标，对应于时刻ti‑1和ti+1的损伤指标MIi‑1和MIi+1表示为：MIi‑1＝|2Di‑1‑Di‑2‑Di|；MIi+1＝|2Di+1‑Di‑Di+2|；考虑i时刻前后的第一个IMF分量斜率均小于Di,则时刻ti‑1和ti+1的损伤指标MIi‑1和MIi+1之和近似等于损伤时刻ti,的损伤指标值MIi，即得出监测因子MIi的约束条件：MIi-1+MIi+1=·MIi;]]>S3‑6，确定损伤监测因子与损伤程度的对应关系建立线性模型：MI＝α*S+β；式中：MI为损伤指标的幅值；S为损伤程度大小；α、β均为线性模型中的常数参数；实际应用过程中首先确定待监测结构的基本信息：质量、刚度和阻尼比；在预先设定的荷载作用下，采用数值分析方法或模型试验模拟结构发生多种不同程度的损伤，然后建立损伤程度与损伤指标的之间的数值关系，并采用数值回归的方法确定参数α和β的数值；由此得到监测因子与损伤程度之间的定量关系。