基于結構撓度和小波變換的梁式結構損傷識別

2014-11-09 05:26:36陳靜嫻

山西建筑 2014年7期

陳靜嫻

(赤峰中天建筑工程有限公司，內蒙古赤峰 024076)

0 引言

工程結構在使用過程中會發生不同程度的損傷，為了不影響結構的正常使用，應當對其是否發生損傷和損傷的具體位置進行識別，進而為制定加固措施提供技術依據。分析時通常將結構簡化為以剛度、質量和阻尼為參數的力學系統，一旦出現損傷，結構的剛度及其固有頻率等特性將會發生變化，進而影響結構的動力響應［1］。但是，由于一般的損傷程度都比較小，即使損傷較大時所得到的結構動態響應曲線也基本平滑，所以不易從結構動力響應上直接判斷結構是否發生損傷。小波變換是一個時間和頻域的局域變換，不僅可看到信號的全貌，又可以表征信號的局部特征，對于那些在結構響應中難以發現的局部損傷信息可以采用小波變換加以分析和識別。

本文利用小波變換的優點，將結構撓度曲線和小波變換有機結合起來，對結構撓度曲線進行小波變換，利用變換系數識別結構的損傷狀態和損傷的具體位置。

1 小波變換及其損傷檢測基本原理

小波變換是通過對Fourier變換的改進而發展起來的一個信號處理工具，其時頻窗口隨著窗口中心的改變，窗口的時寬和頻寬也隨之變換，是一種時—頻局部化分析方法［2-4］。

小波變換的含義是:將某一個稱為基本小波(或者稱為母小波——Mother wavelet)的函數作位移b后，再在不同尺度a下與待分析信號f(t)作內積。一維連續小波函數的定義為:若任意空間函數f(t)是一個可測的，平方可積的函數，L2(R)為f(t)的矢量空間，將f(t)在小波基下展開，稱這種展開為函數f(t)的連續小波變換，簡稱為CWT(Continue Wavelet Transform)，其表達式為:

其中，a為伸縮因子;b為平移因子。所以函數經過小波變換后，就是將一個時間函數投影到二維的時間—尺度相平面上，這樣便于獲取函數信號的某些本質特征。在對信號進行多尺度分析時，如果所用的小波函數是某一光滑函數的一、二階導數，那么在信號突變點處，小波變換后的系數具有過零點或模極大值點的特點，因而可以通過對零點或者模極大值點的檢測來確定損傷的存在及位置。

在工程應用中一般采用離散小波變換，即對小波變換進行二進離散化，取a=2j，b=2j·k。通常將離散化的小波變換稱為二進小波變換。對于任意空間函數f(x)∈L2(－∞，∞)相應的離散小波變換為:

如果一個信號f(x)，它在j尺度下的細節小波系數為:

概貌小波系數為:

則由式(3)，式(4)得到離散小波重構公式為:

其中j尺度的細節函數為:

j尺度的概貌函數為:

信號突變點對應于小波變換系數模的極大值點或過零點，當a=2j，b=k2j時，則:

小波變換的尺度特征j和Lipschitz指數通過式(8)中ja這一項聯系在一起。如果x0點是信號f(x)的突變點，那么在各個尺度上x0點附近的W2jf(x)就會產生一個局部極大值點，而且隨著尺度的減小，這些模的局部極大值點收斂于x0。基于上述理論分析，我們可以利用小波變換的模極大值點來檢測信號f(x)上的突變點。

結構在使用期間一旦發生損傷，就會引起結構剛度的變化，這時結構的空間域信息相應地就會出現十分細微的間斷點。如果對結構的空間域信息進行小波分解，那么所獲得的細節信號中出現的奇異點(小波變換模極大值點)，就是結構的損傷點，從而可以識別結構發生損傷的位置，所以我們可以用小波變換對信號進行奇異性檢測。

當梁損傷后，會造成梁的某一段剛度的降低，損傷識別的目的就是判斷梁體是否存在剛度減少段。在剛度變化截面處，結構滿足變形協調條件和內力平衡條件［5］:

豎向位移:

轉角:

彎矩:

剪力:

由于 EI(d+)≠EI(d－)，所以，即在損傷位置處，位移函數的二階導數不連續，存在第二類間斷點，即為損傷位置。

2 數值模擬

利用Ansys有限元軟件模擬一根長1 m，截面為150 mm×150 mm的懸臂梁，如圖1所示，彈性模量 E=7e10 GPa，泊松比0.3，密度2 700 kg/m3。在數值模擬中，梁損傷可以通過降低單元剛度或者減小截面慣性矩等方法實現，本文采用降低單元彈性模量的方法模擬損傷，損傷程度見表1。將梁分為100個單元，得到梁在各損傷情況下的撓度曲線(見圖2)，然后利用db4小波對其進行4層分解，并分析變換后曲線變化情況。