李雪艷， 汪羽凡

(暨南大學(xué) 力學(xué)與建筑工程學(xué)院 重大工程災(zāi)害與控制教育部重點(diǎn)實(shí)驗(yàn)室，廣州 510632)

1 引 言

損傷識(shí)別是結(jié)構(gòu)健康監(jiān)測(cè)系統(tǒng)的核心技術(shù)，是對(duì)結(jié)構(gòu)進(jìn)行維護(hù)和加固的前提[1,2]。隨著結(jié)構(gòu)安全問題愈來愈受到重視，多種損傷識(shí)別方法也隨之提出。基于振動(dòng)參數(shù)的結(jié)構(gòu)損傷識(shí)別[3-5]，是近年來土木工程的熱點(diǎn)研究課題，振動(dòng)參數(shù)包括頻率、振型、頻響函數(shù)、模態(tài)應(yīng)變能、應(yīng)變響應(yīng)和加速度響應(yīng)等[6-8]。李雪艷等[9,10]提出了基于振動(dòng)響應(yīng)協(xié)方差參數(shù)的損傷識(shí)別方法[9,10]，如利用加速度響應(yīng)協(xié)方差參數(shù)CoC識(shí)別損傷[11-15]，由于加速度響應(yīng)為結(jié)構(gòu)全局振動(dòng)參數(shù)，當(dāng)不使用結(jié)構(gòu)分析模型，只基于結(jié)構(gòu)損傷前后的CoC參數(shù)改變來判定損傷發(fā)生和識(shí)別損傷位置時(shí)，對(duì)于單損傷的識(shí)別效果良好[13],但是對(duì)多損傷的識(shí)別效果較差，容易出現(xiàn)偽損傷和真實(shí)損傷漏判，多損傷的情況下，CoC參數(shù)在接近約束位置處的變化較小，因此更容易受其他位置損傷導(dǎo)致的CoC值變化影響。

當(dāng)結(jié)構(gòu)損傷時(shí)，損傷位置附近將產(chǎn)生應(yīng)力重分布，從而引起應(yīng)變的變化，因此對(duì)比損傷前后的應(yīng)變或者應(yīng)變響應(yīng)參數(shù)，也可以用來識(shí)別結(jié)構(gòu)損傷[14-17]。而應(yīng)變跟結(jié)構(gòu)的位移緊密相關(guān)，所以位移響應(yīng)也將具有較好的結(jié)構(gòu)局部特性，當(dāng)工程結(jié)構(gòu)上未測(cè)得應(yīng)變響應(yīng)，而錄得位移響應(yīng)時(shí)，可以利用位移響應(yīng)計(jì)算出等效應(yīng)變響應(yīng)來進(jìn)行結(jié)構(gòu)損傷識(shí)別。

因此提出了一種位移響應(yīng)協(xié)方差參數(shù)方法，推導(dǎo)并證明了該參數(shù)是結(jié)構(gòu)模態(tài)參數(shù)(頻率，振型和阻尼等)的函數(shù)，是典型的局部性能指標(biāo)，可用來表征結(jié)構(gòu)狀態(tài)。本文將推導(dǎo)建立該協(xié)方差參數(shù)跟結(jié)構(gòu)模態(tài)參數(shù)之間的解析關(guān)系，并通過一個(gè)七層框架結(jié)構(gòu)進(jìn)行數(shù)值模擬，演示該參數(shù)對(duì)于結(jié)構(gòu)局部剛度改變的敏感性，最后對(duì)一個(gè)簡(jiǎn)支鋼梁進(jìn)行實(shí)驗(yàn)測(cè)試，對(duì)錄得的加速度響應(yīng)進(jìn)行二次積分得到位移響應(yīng)，計(jì)算位移響應(yīng)協(xié)方差參數(shù)，并進(jìn)行損傷識(shí)別，結(jié)果表明該方法能成功識(shí)別出損傷位置，有較好的噪聲魯棒性，不依賴結(jié)構(gòu)分析模型。

2 位移響應(yīng)協(xié)方差參數(shù)(CoD)

假定結(jié)構(gòu)第e單元的長度為le,[uiviθiujvjθj]T是平面梁?jiǎn)卧猠的兩個(gè)節(jié)點(diǎn)i和j的位移向量,單元的局部坐標(biāo)跟整體坐標(biāo)之間的夾角為αe,單元e上坐標(biāo)(x,y)處的軸向應(yīng)變可表示為

(1)

因?yàn)檗D(zhuǎn)角位移不容易測(cè)得，所以令式(1)中θi=θj=0,如果用單元中點(diǎn)處(le/2,y)的應(yīng)變響應(yīng)來指示該單元的狀態(tài)，則有

djy(t)sin(αe)-dj x(t)cos(αe)]

(2)

式中di x(t),di y(t),dj x(t)和dj y(t)為節(jié)點(diǎn)i和j處測(cè)得的整體坐標(biāo)x和y方向的位移響應(yīng)，de(t)即為單元e的等效應(yīng)變響應(yīng)。

當(dāng)結(jié)構(gòu)承受沖擊荷載時(shí)，即載荷作用時(shí)間很短，該等效應(yīng)變響應(yīng)的單位脈沖響應(yīng)函數(shù)可以通過式(3)近似得到[14]，以避免傅里葉變換而簡(jiǎn)化計(jì)算。

(3)

單元位移響應(yīng)協(xié)方差參數(shù)(CoD)定義為

(4)

(5)

qi(t)在廣義坐標(biāo)下的表達(dá)式為

qi(t)=(Φf,i/ωd i)e- ξiωitsin(ωd it)

(6)

式中ωi，ωd i和ξi分別為系統(tǒng)的第i階無阻尼頻率、有阻尼頻率和阻尼比，所以CoDe可進(jìn)一步表示為

(7)

(8)

代入式(7)有

(9)

所以CoD的表達(dá)式可寫為

(10)

由式(10)可以看出，CoD僅與振型、頻率和阻尼比等模態(tài)參數(shù)有關(guān)，結(jié)構(gòu)發(fā)生損傷時(shí)會(huì)引起結(jié)構(gòu)物理參數(shù)變化，進(jìn)而影響結(jié)構(gòu)模態(tài)參數(shù)的變化，最終CoD參數(shù)也將變化，所以可通過觀察結(jié)構(gòu)各單元或各測(cè)點(diǎn)的CoD參數(shù)變化，來監(jiān)測(cè)結(jié)構(gòu)的健康狀態(tài)，實(shí)現(xiàn)損傷判定和損傷位置識(shí)別等功能。

3 七層鋼框架結(jié)構(gòu)數(shù)值模擬

對(duì)如圖1所示的七層鋼框架結(jié)構(gòu)進(jìn)行數(shù)值模擬分析，來研究CoD參數(shù)的性質(zhì)和損傷識(shí)別效果。每層豎向柱子分成兩個(gè)等長的梁?jiǎn)卧?，每層橫向梁分成四個(gè)等長的梁?jiǎn)卧?，結(jié)構(gòu)有限元模型的單元和節(jié)點(diǎn)編號(hào)系統(tǒng)如圖2所示，總共56個(gè)平面梁?jiǎn)卧?1個(gè)節(jié)點(diǎn)，每個(gè)節(jié)點(diǎn)3個(gè)自由度，總記153個(gè)自由度，結(jié)構(gòu)每層柱子高為0.3 m，總共2.1 m，柱橫截面為高50 mm、寬8.92 mm的矩形，橫向梁長為0.5 m，截面為高50 mm、寬4.85 mm的的矩形，材料彈性模量為E=206 GPa,結(jié)構(gòu)中柱的質(zhì)量密度和梁的質(zhì)量密度分別為7850 kg/m3和7746 kg/m3,為了模擬樓板的質(zhì)量，每層分別加兩對(duì)質(zhì)量塊，每一對(duì)質(zhì)量塊的質(zhì)量為3.9 kg；框架的底部固支，橫向、豎向和轉(zhuǎn)動(dòng)方向的約束由大剛度1.0×1010kN/m,1.0×1010kN/m和1.0×109kNm/rad來近似，采用瑞利阻尼，前兩階阻尼比為ξ1=ξ2=0.01,其他階阻尼比由ΦTCΦ/ (2ω)計(jì)算得到，Φ為振型矩陣，C為阻尼矩陣，ω為模態(tài)頻率向量，結(jié)構(gòu)在第40節(jié)點(diǎn)承受水平向的三角激勵(lì)，該激勵(lì)的峰值為320.4 N，持續(xù)0.005 s，采樣頻率為2000 Hz，結(jié)構(gòu)前十個(gè)頻率為2.495 Hz,7.503 Hz,12.529 Hz,17.458 Hz,22.049 Hz,25.897 Hz,28.507 Hz,39.999 Hz,42.066 Hz和43.747 Hz，結(jié)構(gòu)的前幾階振型主要為剪切振型。

圖1 七層框架結(jié)構(gòu)及尺寸

3.1 結(jié)構(gòu)剛度改變對(duì)CoD的影響

假定第12單元?jiǎng)偠纫?%的梯度從100%逐漸減少到50%，總計(jì)有11個(gè)損傷狀態(tài)，即剛度減少0%，5%，10%，…，50%。計(jì)算結(jié)構(gòu)每個(gè)狀態(tài)各節(jié)點(diǎn)處的位移響應(yīng)，并用式(4)計(jì)算位移單位脈沖響應(yīng)函數(shù)IRF，進(jìn)而計(jì)算得到每個(gè)狀態(tài)各單元的CoD，11個(gè)狀態(tài)各單元的CoD分布如圖3所示，除了第12單元和第45單元外，11個(gè)狀態(tài)的CoD分布曲線在其他單元處幾乎重合，第12單元為損傷單元，第45單元為其對(duì)稱單元，每個(gè)單元的CoD與損傷程度的關(guān)系曲線如圖4所示，可以看出只有第12、11、13、21、41、42、44、45和46等單元的

圖2 七層框架結(jié)構(gòu)有限元模型單元和節(jié)點(diǎn)編號(hào)

圖3 第12單元?jiǎng)偠葴p少時(shí)CoD的分布

CoD變化顯著，其他單元的CoD曲線呈水平走向，表示隨著損傷程度的增加，變化較少；變化較明顯的9個(gè)單元的CoD曲線重新集中如圖5所示，可以看出，最陡的曲線為第12單元的CoD曲線，而且是單調(diào)遞增的，該研究表明，損傷單元的CoD變化最大，能夠準(zhǔn)確識(shí)別出損傷位置，但是在損傷單元附近和對(duì)稱單元的CoD變化也較明顯。這是由于損傷單元和鄰近單元有共同節(jié)點(diǎn)，損傷單元跟其對(duì)稱單元有相同的振型分量，所以導(dǎo)致鄰近單元和對(duì)稱單元的CoD的改變也較大，造成誤判。

對(duì)第6單元也做同樣的分析研究，即以5%的遞度,剛度減少從0%增加到50%，總共11個(gè)結(jié)構(gòu)損傷狀態(tài)，11個(gè)狀態(tài)各單元CoD的分布如圖6所示，可以看出只有第6、51和52單元，由于損傷程度不同，而有CoD的顯著變化，不同損傷狀態(tài)時(shí)CoD分布曲線在其他單元幾乎重合，九個(gè)變化顯著單元的CoD與損傷程度的關(guān)系曲線如圖7所示，可以看出，第6單元的CoD有最顯著的變化，其次是第51和52單元，它們是單元6的對(duì)稱單元；說明單元?jiǎng)偠葴p少會(huì)使該單元的CoD產(chǎn)生最大的改變，所以利用CoD參數(shù)可以進(jìn)行結(jié)構(gòu)損傷識(shí)別。

圖4 第12單元?jiǎng)偠葴p少時(shí)各單元CoD與損傷程度的關(guān)系曲線

圖5 第12單元?jiǎng)偠葴p少時(shí)九個(gè)變化顯著單元的CoD與損傷程度的關(guān)系曲線

圖6 第6單元?jiǎng)偠葴p少時(shí)CoD的分布

3.2 損傷識(shí)別

利用CoD進(jìn)行損傷識(shí)別，結(jié)構(gòu)仍然在第40節(jié)點(diǎn)處承受水平方向的三角沖擊激勵(lì)。

損傷工況1： 第9單元?jiǎng)偠葴p少10%。測(cè)得該損傷狀態(tài)下結(jié)構(gòu)各節(jié)點(diǎn)的位移，并在位移響應(yīng)中添加10%的白噪聲來模擬測(cè)試噪聲，然后計(jì)算位移脈沖響應(yīng)函數(shù)IRF，再得到各單元的CoD，把該損傷狀態(tài)時(shí)各單元的CoD減去結(jié)構(gòu)未損傷時(shí)各單元的CoD，得到各單元CoD的變化值，并進(jìn)行歸一化得到如圖8所示的損傷向量，可以看出第9單元的CoD有最大的改變，因此有最大的損傷概率，而其對(duì)稱單元(第48單元和第47單元)有第二和第三大的CoD變化和損傷概率，因此通過CoD損傷前后的改變，可以成功識(shí)別出該損傷位置。

損傷工況2： 第3和第48單元均有10%的剛度減少。同上得到位移響應(yīng)，并添加白噪聲和計(jì)算IRF，得到各單元的CoD值，最后得到各單元CoD變化值，并進(jìn)行歸一化得到損傷概率向量，如圖9所示，可以看出第3單元和第48單元有最顯著的CoD變化和最大的損傷概率，而對(duì)稱單元(第9、10、53和54單元)有較大的損傷概率，其他單元的損傷概率很小。表明該方法也能成功識(shí)別出兩個(gè)損傷時(shí)的損傷位置。

圖7 第6單元?jiǎng)偠葴p少時(shí)九個(gè)變化顯著單元的CoD與損傷程度的關(guān)系曲線

損傷工況3： 第4單元和第8單元均有15%的剛度減少。進(jìn)行同上的分析，可以得到如圖10所示的損傷概率向量，仍然可以看出，損傷單元第4和第8單元有最大的損傷概率，對(duì)稱單元有第二大的損傷概率，而其他單元的損傷概率非常小，再一次表明該方法能成功識(shí)別出多損傷時(shí)的損傷位置。

圖8 第9單元?jiǎng)偠葴p少時(shí)由CoD改變得到的損傷概率向量

圖9 第3和48單元?jiǎng)偠葴p少時(shí)由CoD改變得到的損傷概率向量

圖10 第4和8單元?jiǎng)偠葴p少時(shí)由CoD改變得到的損傷概率向量

通過上述損傷工況研究發(fā)現(xiàn)，雖然可能會(huì)對(duì)損傷單元的對(duì)稱單元產(chǎn)生誤判，但是CoD參數(shù)對(duì)結(jié)構(gòu)剛度減少敏感，對(duì)噪聲魯棒，不依賴結(jié)構(gòu)分析模型，是很好的結(jié)構(gòu)健康狀態(tài)監(jiān)測(cè)指標(biāo)。

4 實(shí)驗(yàn)驗(yàn)證

對(duì)圖11所示的簡(jiǎn)支鋼梁進(jìn)行實(shí)驗(yàn)室測(cè)試，來進(jìn)一步驗(yàn)證本文提出的損傷識(shí)別方法，鋼梁長為1996 mm，截面是寬為50.75 mm,高為9.69 mm的矩形，楊氏模量為191.1 GPa，密度為7790.6 kg/m3,鋼梁兩端簡(jiǎn)支，支座間跨度為1920 mm，由于加速度響應(yīng)更容易測(cè)量，而且數(shù)據(jù)質(zhì)量更高，所以本實(shí)驗(yàn)通過測(cè)量加速度響應(yīng)，再通過數(shù)值二次積分得到位移響應(yīng)來進(jìn)行損傷識(shí)別，七個(gè)加速度傳感器等間距安裝在梁的上表面,如圖11所示，采樣頻率為2000 Hz，在距梁右端638 mm處的上表面，用錘子進(jìn)行敲擊產(chǎn)生振動(dòng)，使用商用數(shù)據(jù)采集系統(tǒng)INV303和數(shù)據(jù)分析系統(tǒng)DASP2003測(cè)試并記錄加速度響應(yīng)，進(jìn)行多次重復(fù)測(cè)試以后，再在距離梁右端158 mm左邊處，梁的前后表面鋸掉長9 mm、高9.69 mm、寬0.9 mm的缺口，如圖12所示，來制造損傷，并測(cè)試記錄加速度響應(yīng)。

來自第一個(gè)傳感器的未損傷和損傷狀態(tài)下的 2 s 內(nèi)的加速度響應(yīng)如圖13所示，對(duì)測(cè)得的加速響應(yīng)進(jìn)行實(shí)驗(yàn)?zāi)B(tài)分析，得到前10個(gè)頻率和阻尼比列入表1，可知兩個(gè)狀態(tài)的頻率變化很小，但是第1、2、7和第8階阻尼比有較大的改變。

把未損傷狀態(tài)和損傷狀態(tài)下的各個(gè)傳感器的加速度響應(yīng)對(duì)時(shí)間進(jìn)行數(shù)值二次積分，得到位移響應(yīng)，來自第一個(gè)傳感器的未損傷和損傷狀態(tài)下的位移響應(yīng)如圖14所示，可以看出兩個(gè)狀態(tài)下的位移響應(yīng)，頻率幾乎重合，但是幅值有一定的改變，再利用式(2，4)得到各個(gè)傳感器的CoD的改變，并進(jìn)行歸一化，如圖15所示的損傷向量，可以看出第7個(gè)傳感器有最大的CoD改變和損傷概率，表示損傷發(fā)生在第7個(gè)傳感器附近，與實(shí)際損傷位置一致，表明該方法能成功判定損傷發(fā)生和識(shí)別損傷 位置。

圖11 實(shí)驗(yàn)室測(cè)試的簡(jiǎn)支鋼梁

圖12 簡(jiǎn)支鋼梁損傷切口

圖13 第一個(gè)傳感器測(cè)得的結(jié)構(gòu)未損傷和損傷后的加速度響應(yīng)

表1 實(shí)驗(yàn)測(cè)試得到的頻率和阻尼

圖14 由第一個(gè)傳感器測(cè)得的結(jié)構(gòu)未損傷和損傷后的加速度響應(yīng)二次積分得到的位移響應(yīng)

圖15 損傷向量

5 結(jié) 論

本文通過理論推導(dǎo)得到位移響應(yīng)協(xié)方差參數(shù)CoD的解析公式，表明它是結(jié)構(gòu)固有參數(shù)(頻率、振型和阻尼)的函數(shù)，結(jié)構(gòu)物理參數(shù)的改變會(huì)導(dǎo)致CoD參數(shù)改變，所以可利用CoD參數(shù)來進(jìn)行結(jié)構(gòu)損傷識(shí)別。

通過對(duì)七層框架結(jié)構(gòu)數(shù)值模擬分析，發(fā)現(xiàn)損傷單元的剛度減少對(duì)該單元的CoD改變最大，其次是對(duì)稱和附近單元，而且隨著損傷程度的增加，損傷單元的CoD變化也增加，通過對(duì)單損傷和多損傷等三種損傷工況的分析，表明該方法能成功判定損傷發(fā)生和準(zhǔn)確識(shí)別損傷位置，只需比較結(jié)構(gòu)損傷前后單元的CoD改變，無需結(jié)構(gòu)分析模型，而且具有較好的噪聲魯棒性。

最后對(duì)實(shí)驗(yàn)室測(cè)試的簡(jiǎn)支鋼梁進(jìn)行分析研究，人為引入損傷，比較結(jié)構(gòu)損傷前后的頻率和阻尼比，發(fā)現(xiàn)其變化較小，難以判定損傷發(fā)生，而利用本文提出方法，對(duì)錘擊激勵(lì)下的加速度響應(yīng)數(shù)據(jù)進(jìn)行數(shù)值二次積分，得到位移響應(yīng)，比較結(jié)構(gòu)損傷前后各傳感器CoD數(shù)據(jù)的改變，得到損傷概率向量，損傷位置附近傳感器的CoD有最大的改變，表明該方法能成功識(shí)別出損傷位置。

本文提出方法無需結(jié)構(gòu)分析模型，簡(jiǎn)單易算，具有較好的工程應(yīng)用前景。