桁架結構損傷預測方法初探

陳詩超　陳建飛

(陽光學院土木工程系　福建福州　350015)

?

桁架結構損傷預測方法初探

陳詩超陳建飛

(陽光學院土木工程系福建福州350015)

綜合運用模態曲率變化量與模態應變能變化率，提出了桁架結構損傷預測的方法。首先，應用ANSYS有限元軟件建立梁式桁架的模型，分析其位移模態并模擬其損傷，并通過分析桁架上下弦的模態曲率變化量預測損傷區域。其次，計算損傷區域內桿件的模態應變能變化率，并將應變能變化率較大的桿件認定為損傷單元。數值模擬結果表明，該法可準確定位梁式桁架的損傷部位。

損傷預測；桁架結構；模態曲率變化量；模態應變能變化率

0　引言

相對于其他結構型式，桁架結構可以極大地節約材料。該結構可以在相同的豎向撓度限制下獲得更大的跨度，并且可以方便、迅速地安裝，因此在大跨度鐵路橋中大量運用了該結構型式。隨著對結構健康監測需求的增長，越來越多的學者開始關注結構損傷的識別，也因此提出了一系列的方法[1-2]。總體來說，模態曲率變化量法廣泛用于梁式結構的損傷鑒別[3-5]，模態應變能變化率法則用于桁架結構[6-8]。

梁式結構中桿件眾多，下弦桿的模態曲率變化量僅能甄別可能的破壞區域，而不能定位破壞桿件。雖然模態應變能變化率可以鑒別破壞的桿件，但需要分析的桿件數量較大，同時分析結果并非十分準確[6]。因此，本文將模態曲率變化量法與單元應變能變化率法相結合，提出了結構損傷預測的兩步法。首先，按照第一階模態的曲率變化定位可能的損傷區域，隨后計算可能損傷區域內桿件的應變能變化率，以最終確定真正的損傷桿件。所進行的數值模擬表明：本文提出的損傷預測法適用于梁式桁架結構損傷的預測。

1　基本原理

1.1模態曲率(MC)

學術界多采用結構的位移描述其振動模態，本文所采用的模態曲率，其作為受彎結構中性層變形的描述，則顯得更為精確。直梁軸線的曲率可以通過下列彎曲變形近似微分方程計算：

(1)

式中的x是沿梁長度方向的坐標，y則是梁的撓度。式(1)在節點i處的微分形式則可寫成：

(2)

在本文中，將較大模態曲率變化量節點所圍成的區域作為可能的損傷區域。

1.2模態應變能變化率(MSECR)

隨著結構剛度的降低，損傷單元的模態應變能將隨之改變。由于第一步中已經定位了損傷區域，因此對其中的桿件進行應變能變化量的分析是可行的。模型第j單元i階模態的應變能可以表示為如下形式：

MSEij={φi}T[Kj]{φi}

(3)

式(3)中Kj為第j單元的結構剛度矩陣，φi為第i階模態，另外d表示損傷的結構。

由于Kj、φi均過于龐大，為了節約計算資源，可以采用單元剛度矩陣ki代替Ki，同時采用單元節點向量vi代替φi：

MSEij={vij}T[kj]{vij}

(4)

因此應變能改變率(MSECRij)可以按照下式計算

(5)

實際應用的過程中，可以采用多階模態的復合應變能改變率(MSECRj)減小噪音：

(6)

損傷單元即為最大復合應變能改變率(MSECRj)所在的單元。

2　數值模擬

本文采用ANSYS，建立了一個兩端鉸支的梁式桁架的數值模型，模型具體尺寸如圖1所示，桿件的材料屬性詳見表1。該數值模型采用35個Link1單元模擬桿件，計算所得未損傷結構的一二階模態如圖2所示。研究過程中通過降低桿件的剛度模擬結構的損傷。

2.1下弦桿損傷的模擬

假設3號桿件受到損傷后，剛度分別降低5%、10%、 20% 及 50%，上下弦桿第一階模態曲率變化量如圖3所示。由圖3可知，3號、4號節點數值較大，這就意味著下弦桿的3、4號節點位于破壞區域。對于上弦桿，顯然節點13位于破壞區域。因此破壞區域由3號、4號、13號節點圍成，即破壞桿件應該在3號、14號、15號桿件中。

接著計算3號、14號、15號桿件的模態應變能變化率，并將結果列于表2。對于不同的損傷工況，均可以按照最大的模態應變能變化率確定3號桿為損傷桿件。

2.2邊界桿件的損傷模擬

由于采用了微分算法，實際上我們無法得到支撐節點的模態曲率變化量。如果邊界位置的桿件發生破壞，上下弦桿間的模態曲率變化量峰值點將在一條直線上，而不是圍成一個區域。現假設10號桿件受到損傷，且損傷程度為5%、10%、20% 及 50%。上、下弦桿的第一階模態曲率變化量如圖4所示。由圖4可見：模態曲率變化量的峰值出現在2號及12號節點，因此破壞區域應處于節點的連線左側，即可能的破壞桿件為1號、10號、11號、28號桿。

1號、10號、11號、28號桿的模態應變能變化率的計算如表3所示，表中數據可知：損傷單元應該為具有最大值的10號桿件。

2.3多位置損傷

假設2號、33號桿件受到損傷，損傷程度為10%。上下弦桿的第一階模態曲率變化量如圖5所示。由圖5可知：上弦桿的模態曲率變化量峰值出現于12號、16號、17號節點，下弦桿的峰值這出現于2號、3號、7號節點。此時，由這些節點所圍成的損傷區域有兩個，可能的損傷桿件分別為為2號、12號、13號和21號、22號、23號桿件。將這些桿件的模態應變能變化率計算結果列于表4，由計算結果可知：損傷桿件為具有較大變化率的2號、3號桿件。

3　結論

本文通過將模態曲率變化量法和模態應變能變化率法結合，提出了梁式桁架損傷的預測方法。方法將損傷預測分為兩部分：其一，通過模態曲率變化量定位可能的損傷區域；其二，依據模態應變能變化率準確確定損傷桿件。本文的數值模擬結果表明：無論桁架損傷的位置如何，本文所提供的方法都可以準確定位損傷的桿件。

[1]YU Ling, ZHU Junhua, YU Lili. Structural Damage Detection in a Truss Bridge Model Using Fuzzy Clustering and Measured FRF Data Reduced by Principal Component Projection [J]. Advances in Structural Engineering, 2013, 16:207-217.

[2]Sazonov E, Klinkhachorn P. Optimal Spatial Sampling Interval for Damage Detection by Curvature or Strain energy Mode Shapes [J]. Advances in Structural Engineering, 2013, 16:207-217. Journal of Sound and Vibration, 2005, 285:783-801.

[3]Dawari V B, Vesmawala G R. Modal Curvature and Modal Flexibility Methods for Honeycomb Damage Identification in Reinforced Concrete Beams [J]. Procedia Engineering, 2013, 51:119-124.

[4]Baghiee N, Esfahani M R, Moslem K. Studies on Damage and FRP Strengthening of Reinforced Concrete Beams by Vibration Monitoring[J]. Engineering Structures, 2009, 31: 875-893.

[5]葉黔元, 翟立祥, 陳錦輝. 曲率模態及其在桁架橋梁損傷識別中應用 [J].力學季刊, 2005, 26:286-292.

[6]Seyedpoor S M. A Two Stage Method for Structural Damage Detection Using a Modal Strain Energy Based Index and Particle Swarm Optimization [J]. International Journal of Non-Linear Mechanics, 2012, 47:1-8.

[7]Yan W J, Ren W X, Huang T L. Statistic Structural Damage Detection Based On the Closed-form of Element Modal Strain Energy Sensitivity [J]. Mechanical Systems and Signal Processing, 2012, 28:183-194.

[8]史治宇, 羅紹湘, 張令彌. 結構破損定位的單元模態應變能變化率法 [J].振動工程學報, 1998, 11:356-360.

陳詩超(1988.03-)，男，助教，主要從事沖擊動力學、結構健康監測與檢測方向的研究。

陳建飛(1987.12-)，男，助教，主要從事土體注漿加固方向的研究。

A Prieliminary Study on Damage Forecast Method for Truss Structure

CHENShichaoCHENJianfei

(Department of Civil Engineering,Yango College,Fuzhou 350015)

By integrated use of Variation of modal curvature and Rate of modal strain energy changing，a damage forecast method for truss structure is proposed. Firstly, a finite element model of truss structure is established with ANSYS，after which displacement modes are got and damage is simulated，the damage region is estimated via analyzing the variation of modal curvature of top and bottom chords of truss structure. Subsequently, the modal strain energy changing rate of members in the damage region is calculated and the element with the maximum modal strain energy changing rate is deemed as the damaged element. The simulation results show that this method can accurately locate the damage in the beam-like truss structure.

Damage forecast；Truss structure；Variation of modal curvature; Rate of modal strain energy changing

陳詩超(1988.03-) ，男，助教。

E-mail:ligong0601@163.com

2015-12-01

TU3

A

1004-6135(2016)02-0043-05