基于曲率模態(tài)法和柔度矩陣法的斜拉橋損傷識別對比研究

楊子杰　王華

摘要：現(xiàn)有橋梁損傷識別方法均自成一體，缺乏相互間的橫向比較，識別效果的優(yōu)劣難以確定。文章采用ANSYS軟件建立不同位置、不同程度損傷的斜拉橋模型，對斜拉橋主梁的動力特性進行了分析，對比研究了兩種損傷識別方法（曲率模態(tài)和柔度矩陣法）的靈敏性和定位準確性。結(jié)果表明：曲率模態(tài)法和柔度矩陣法對于斜拉橋主梁微小損傷均具有良好的靈敏性;曲率模態(tài)差隨著損傷的嚴重程度而更大，柔度矩陣的變化率也更大，兩種方法均能識別不同程度的損傷;通過觀察曲率模態(tài)差在損傷位置的突變更易定位損傷位置。

關(guān)鍵詞：曲率模態(tài);柔度矩陣;斜拉橋;損傷識別

中國分類號：U448.27文章標識碼：A250954

0 引言

大型橋梁的健康檢測是保障橋梁服役安全的重要手段，對橋梁損傷進行識別、定位是橋梁健康檢測的核心任務(wù)，損傷識別方法極為重要，其決定著橋梁健康檢測的有效性和可靠性。

目前，發(fā)展較為成熟的大型橋梁損傷識別方法主要是基于動力學參數(shù)變化的方法，主要有頻率法、振型、曲率模態(tài)法、柔度法和頻響函數(shù)法等[1]。曲率模態(tài)法和柔度矩陣法被普遍認為是具有較好工程應(yīng)用前景的損傷識別方法[2-3]，但二者對橋梁損傷識別特別是斜拉橋的適用性仍有待進一步探討。此外，現(xiàn)有損傷識別方法的研究均自成一體，缺乏相互間的橫向比較，識別效果的優(yōu)劣難以確定。

為此，本文采用ANSYS建立不同位置、不同程度損傷的雙塔斜拉橋模型，對斜拉橋主梁的動力特性進行了分析，對比分析兩種損傷識別方法（曲率模態(tài)和柔度矩陣法）的靈敏性和定位準確性，深入了解這兩種方法在斜拉橋主梁損傷識別中的適用性。

1 損傷識別理論

1.1 曲率模態(tài)法

由材料力學和微分學原理可知，位移的二階導數(shù)即為曲率，必有固有曲率分布狀態(tài)與每階位移模態(tài)相對應(yīng)，這種固有曲率分布狀態(tài)即為曲率模態(tài)[4-5]。

由以上分析可知道，如果斜拉橋主梁某個單元發(fā)生損傷，這個單元處的柔度矩陣變化率將發(fā)生突變，形成極大值。因此，可以通過柔度矩陣極大值來識別和定位斜拉橋主梁損傷。

2 有限元模型

本文以雙塔斜拉橋作為研究對象，采用ANSYS建立不同位置、不同程度損傷的斜拉橋模型，并對斜拉橋主梁的動力特性進行了分析，對比分析兩種損傷識別方法（曲率模態(tài)和柔度矩陣法）的靈敏性和定位準確性。用梁單元模擬雙塔斜拉橋的主梁，采用等效為剛度降低的方法來模擬斜拉橋主梁的損傷。

橋梁的跨徑為（171+342+171）m。主梁寬度為28 m，高度為1.6 m，塔高為160 m。主梁采用shell63單元，拉索采用link10單元。塔底固定約束，梁端采用豎向約束。以4 m為一個單元，主梁共劃分171個單元，計算該斜拉橋的前4階模態(tài)。斜拉橋的有限元模型如圖1所示。

3 損傷定位的靈敏性分析

主梁是斜拉橋的關(guān)鍵構(gòu)件之一，直接承受著車輛和人群載荷的作用。在實際橋梁中，對主梁的損傷識別尤為不易，因此，本文選取斜拉橋主梁作為損傷識別的對象。由理論分析可知，應(yīng)用曲率模態(tài)差和柔度矩陣變化率識別斜拉橋主梁損傷，則損傷發(fā)生的位置曲率模態(tài)差和柔度矩陣變化率會發(fā)生明顯變化。

為了對比分析曲率模態(tài)差和柔度矩陣變化率在斜拉橋主梁損傷識別中的靈敏性，在斜拉橋主梁單元81位置設(shè)置了微小損傷（剛度降低5%）（工況一），曲率模態(tài)差和柔度矩陣變化率計算結(jié)果如圖2、圖3所示。

從圖2可以看出，斜拉橋主梁81單元位置發(fā)生微小損傷（5%），1～4階曲率模態(tài)差在損傷單元位置發(fā)生明顯突變，通過曲率模態(tài)法可靈敏地識別微小損傷且能夠進行定位。從圖3可以看出，柔度矩陣變化率在微小損傷位置出現(xiàn)明顯峰值。兩種方法對于斜拉橋主梁微小損傷都具有較高的靈敏性。

4 損傷定位準確性對比分析

為驗證損傷定位的準確性，可通過觀察模型曲率模態(tài)差和柔度矩陣變化率與損傷位置是否一致來對比分析。本文通過設(shè)置主梁單損傷和多損傷兩種工況來對比分析兩種損傷識別方法對于損傷定位的準確性。

工況二：單元81（81、82兩節(jié)點間）發(fā)生30%的損傷。

工況三：單元81發(fā)生10%的損傷，單元131（131、132兩節(jié)點間）發(fā)生30%的損傷。

工況二曲率模態(tài)法差計算如圖4所示，柔度矩陣變化率計算如圖5所示。從圖4可以看出，主梁出現(xiàn)損傷，通過第一至第四階曲率模態(tài)差均能夠明顯看出在節(jié)點78～83范圍內(nèi)出現(xiàn)損傷，且隨著階數(shù)越大，曲率模態(tài)差絕對值越大。從圖5可以看出，節(jié)點81單元柔度矩陣變化率遠大于其他單元的變化率。對于單損傷，曲率模態(tài)差和柔度矩陣變化率均能夠較好識別。

工況三曲率模態(tài)法差計算如下頁圖6所示，柔度矩陣變化率計算如下頁圖7所示。從圖6可以看出，曲率模態(tài)法的第一階曲率模態(tài)能夠較好識別兩處損傷，但隨著階數(shù)增大，逐漸出現(xiàn)了一些偽損傷位置。從圖7可以看出，在81單元和131單元處均形成了兩個明顯峰值，柔度矩陣變化率能夠清晰地識別這兩處損傷。

因此，針對損傷定位準確性，對比兩種方法可知：

（1）斜拉橋主梁發(fā)生多處損傷情況下，1～4階曲率模態(tài)差在多處損傷單元位置均突變明顯。

損傷越大，對應(yīng)的損傷單元位置突變的值越大，在對損傷的識別方面，兩種方法對一處或多處不同程度的損傷均能識別。

（2）柔度矩陣變化率在損傷單元位置出現(xiàn)極大值。

（3）通過觀察曲率模態(tài)差在損傷位置的突變更易定位損傷位置。

5 結(jié)語

本文采用ANSYS建立斜拉橋模型，設(shè)置主梁不同位置損傷，對其進行動力特性分析，通過曲率模態(tài)差和柔度矩陣變化率識別主梁損傷，對比兩種方法對損傷定位的準確性和靈敏性進行分析。具體分析結(jié)論如下：

（1）在發(fā)生微小損傷情況下，1～4階曲率模態(tài)差在損傷單元位置突變明顯，且隨著階數(shù)增加，突變更加明顯。柔度矩陣變化率在損傷單元位置達到峰值，曲率模態(tài)差和柔度矩陣變化率對于斜拉橋主梁損傷均具有良好的靈敏性。

（2）斜拉橋主梁發(fā)生多處損傷情況下，1～4階曲率模態(tài)差在多處損傷單元位置均突變明顯。損傷越大，對應(yīng)的損傷單元位置突變的值越大，在對損傷的識別方面，兩種方法對一處或多處不同程度的損傷均能識別。柔度矩陣變化率在損傷單元位置出現(xiàn)極大值。通過觀察曲率模態(tài)差在損傷位置的突變更易定位損傷位置。

參考文獻：

[1]李帛書. 基于優(yōu)化智能算法與低階曲率模態(tài)的斜拉橋損傷識別[D].重慶：重慶大學，2019.

[2]李忠獻，齊懷展，朱勁松.基于模態(tài)曲率法的大跨度斜拉橋損傷識別[J].地震工程與工程振動，2007（4）：122-126.

[3]荊龍江，項貽強.基于柔度矩陣法的大跨斜拉橋主梁的損傷識別[J].浙江大學學報（工學版），2008（1）：164-169.

[4]范小寧，杜晨輝，羅志宏.基于小波變換和模態(tài)曲率差的起重機金屬結(jié)構(gòu)損傷識別[J].安全與環(huán)境學報，2021，21（1）：147-153.

[5]肖儀清，李成濤.基于曲率模態(tài)和神經(jīng)網(wǎng)絡(luò)的斜拉橋損傷識別[J].武漢理工大學學報，2010，32（9）：275-279.

[6]趙 玲，李愛群.基于模態(tài)柔度矩陣變化指標的結(jié)構(gòu)損傷預(yù)警方法[J].東南大學學報（自然科學版），2009，39（5）：1 049-1 053.