秦偉亮， 孫宗光， 曲翔雨， 司琦， 倫卓

(大連海事大學(xué) 交通運(yùn)輸工程學(xué)院， 遼寧 大連 116023)

大跨橋梁結(jié)構(gòu)的使用期長(zhǎng)達(dá)幾十年、甚至上百年，環(huán)境侵蝕、材料老化、荷載的長(zhǎng)期效應(yīng)、疲勞效應(yīng)與突變效應(yīng)等災(zāi)害因素的耦合作用將不可避免地導(dǎo)致結(jié)構(gòu)和系統(tǒng)的損傷積累和抗力衰減, 從而抵抗自然災(zāi)害、甚至正常環(huán)境作用的能力下降，極端情況下引發(fā)災(zāi)難性的突發(fā)事故。僅憑以往的間歇性的檢測(cè)和評(píng)價(jià)，難以滿(mǎn)足確保結(jié)構(gòu)安全運(yùn)營(yíng)的要求。因此，近20年來(lái)結(jié)構(gòu)長(zhǎng)期健康監(jiān)測(cè)系統(tǒng)，得以迅速發(fā)展。如何利用結(jié)構(gòu)健康監(jiān)測(cè)系統(tǒng)采集到的響應(yīng)信息反演出橋梁結(jié)構(gòu)的工作狀態(tài)，這是實(shí)現(xiàn)大跨橋梁結(jié)構(gòu)損傷識(shí)別與狀態(tài)評(píng)估的關(guān)鍵問(wèn)題。為此，眾多學(xué)者對(duì)其進(jìn)行了大量的研究。斜拉索作為斜拉橋主要的承力構(gòu)件，連接著主梁和主塔并且分布廣泛，本身含有橋面和主塔結(jié)構(gòu)的空間位置信息，而對(duì)于斜拉索和索力的相關(guān)研究大多停留在斜拉索發(fā)生損傷的識(shí)別上。文獻(xiàn)[6]首次提出了索力損傷識(shí)別指標(biāo)，并探討了對(duì)不同損傷工況的敏感性和穩(wěn)定性，顯示了較好的基于該指標(biāo)進(jìn)行斜拉橋損傷識(shí)別的可行性；文獻(xiàn)[10、11]采用了索力損傷指標(biāo)，但損傷識(shí)別僅是對(duì)主梁的損傷識(shí)別。該文以青島海灣大橋滄口斜拉橋?yàn)榉治鰧?duì)象，面向?qū)嶋H監(jiān)測(cè)索力，提出了基于統(tǒng)計(jì)特性而構(gòu)建的索力損傷向量。基于經(jīng)修正后的有限元模型，研究了基于該索力損傷向量進(jìn)行單損傷和組合損傷識(shí)別分析方法。該方法可通過(guò)對(duì)全部或部分斜拉索的索力監(jiān)測(cè)，實(shí)現(xiàn)對(duì)斜拉橋結(jié)構(gòu)的一部分損傷進(jìn)行識(shí)別和健康狀態(tài)評(píng)價(jià)。

1 模型修正與監(jiān)測(cè)索力模擬

1.1 模型修正

滄口橋?qū)嵕皥D如圖1所示，1/2跨立面布置如圖2所示。首先，運(yùn)用Ansys有限元分析軟件創(chuàng)建初始有限元模型；其中，主梁采用Beam4單元，主塔和墩采用Beam188單元，斜拉索采用Link10單元。然后，以荷載試驗(yàn)和部分監(jiān)測(cè)結(jié)果為目標(biāo)，采用響應(yīng)面法進(jìn)行靜、動(dòng)力模型修正。修正后模型計(jì)算的前幾階主要模態(tài)頻率和靜態(tài)索力值與實(shí)測(cè)值的比較如表1、2所示。

圖1 滄口橋?qū)嵕皥D

圖2 滄口橋1/2跨立面圖及索的編號(hào)(單位：m)

表1 修正后模型與實(shí)測(cè)振型頻率對(duì)比

表2 靜態(tài)索力實(shí)測(cè)值與模擬值對(duì)比 kN

1.2 監(jiān)測(cè)索力模擬

采用基于監(jiān)測(cè)索力構(gòu)建的指標(biāo)進(jìn)行損傷識(shí)別，宜選用相對(duì)穩(wěn)定可靠的索力監(jiān)測(cè)數(shù)據(jù)。正常運(yùn)營(yíng)環(huán)境下索力主要受溫度、車(chē)輛和其他隨機(jī)干擾因素的影響。

為此，首先選用夜間時(shí)段數(shù)據(jù)以盡量消除日照溫度的影響。再通過(guò)數(shù)據(jù)分解分離高頻隨機(jī)干擾部分和車(chē)輛荷載導(dǎo)致的索力突變。然后在年周期內(nèi)選取每天同一時(shí)段的數(shù)據(jù)為樣本，分析建立年溫度與索力的統(tǒng)計(jì)關(guān)系。分析表明：年溫度與索力之間呈現(xiàn)顯著的線(xiàn)性相關(guān)性?；诒O(jiān)測(cè)年溫度數(shù)據(jù)樣本及其分布特性進(jìn)行溫度加載，由修正的有限元模型計(jì)算得到模擬索力。年溫度與實(shí)測(cè)索力和模擬索力的關(guān)系對(duì)比如圖3所示。實(shí)測(cè)索力和模擬索力二者在分布規(guī)律及與溫度的相關(guān)關(guān)系方面吻合較好。表明修正的有限元模型能夠較好地反映索力與環(huán)境溫度的關(guān)系。

圖3 實(shí)測(cè)與模擬索力-環(huán)境溫度關(guān)系比較

由于全年索力F與溫度T之間具有較高的線(xiàn)性相關(guān)關(guān)系，可通過(guò)式(1)消除溫度對(duì)索力的影響，把得到的新索力Fs稱(chēng)為標(biāo)準(zhǔn)索力。

Fs=F-k(T-Tb)

(1)

式中：Fs為標(biāo)準(zhǔn)索力；F為索力；k為索力與環(huán)境溫度關(guān)系的斜率；T為環(huán)境溫度；Tb為基準(zhǔn)溫度(此處取15 ℃)。

標(biāo)準(zhǔn)正態(tài)分布檢驗(yàn)結(jié)果表明：標(biāo)準(zhǔn)索力Fs符合正態(tài)分布。2012年W7索全年實(shí)測(cè)索力和模擬索力的標(biāo)準(zhǔn)值分別如圖4、5所示。以下的損傷識(shí)別分析將基于標(biāo)準(zhǔn)索力。

圖4 W7索的標(biāo)準(zhǔn)索力實(shí)測(cè)值

圖5 W7索的標(biāo)準(zhǔn)索力模擬值

2 損傷敏感性分析

通常情況下拉索索力對(duì)該區(qū)域的損傷更敏感，因此根據(jù)結(jié)構(gòu)對(duì)稱(chēng)性將斜拉橋劃分為4個(gè)區(qū)域，針對(duì)各區(qū)域的損傷工況采用該區(qū)域索力進(jìn)行損傷分析和識(shí)別。這里以一個(gè)區(qū)域?yàn)槔齺?lái)闡述損傷分析和識(shí)別方法。首先基于斜拉橋易損性分析的相關(guān)文獻(xiàn)以及根據(jù)該橋靜載情況下的易損性分析，在該區(qū)域選取8種損傷工況。其中工況1、2為索損傷，3、4為主塔損傷，5、6、7為主梁損傷，8為整體老化。前7種工況的損傷位置如圖6所示。

圖6 第1～7種損傷工況的損傷位置(單位：m)

模擬損傷方法很多，該文采用降低構(gòu)件彈性模量的方法模擬損傷。為得到標(biāo)準(zhǔn)索力分布特性與損傷程度之間的關(guān)系以及索力均值隨損傷程度的變化趨勢(shì)和趨勢(shì)方程，分別模擬了各損傷工況0～70%的損傷程度。計(jì)算結(jié)果表明：索力均值對(duì)這幾種損傷工況都較敏感，而索力標(biāo)準(zhǔn)差敏感性較差。限于篇幅，這里僅列出工況2的W10索標(biāo)準(zhǔn)索力分布特性變化情況(如圖7所示)，以及相應(yīng)的索力均值改變量與損傷程度的趨勢(shì)(如圖8所示)。

圖7 W10索標(biāo)準(zhǔn)索力分布特性(工況2)

圖8 W10索力均值改變量與損傷程度趨勢(shì)線(xiàn)(工況2)

由圖8可見(jiàn)索力均值隨損傷程度呈顯著變化，而標(biāo)準(zhǔn)差最大改變量?jī)H為為-0.559 4 kN；故選用索力均值作為損傷識(shí)別指標(biāo)。

由損傷敏感性分析可得：斜拉索不僅平衡主梁自重更使主梁產(chǎn)生拱度，而主梁具有一定的抗彎剛度；當(dāng)主梁發(fā)生損傷，索力必然發(fā)生變化。主塔兩側(cè)索力也不完全平衡，從而在橋塔內(nèi)產(chǎn)生可觀(guān)的彎矩，故主塔損傷也必然引起索力變化。對(duì)于斜拉索，當(dāng)拉索發(fā)生損傷必然引起索力重分布。隨著時(shí)間推移，斜拉橋結(jié)構(gòu)構(gòu)件會(huì)發(fā)生老化致使剛度降低(包括梁、塔、索)，索力必然發(fā)生明顯降低。

3 損傷識(shí)別

采用標(biāo)準(zhǔn)索力個(gè)別值抽樣例行監(jiān)測(cè)評(píng)價(jià)，該方法較為簡(jiǎn)單，但有很大的局限性。如果結(jié)構(gòu)發(fā)生顯著改變，必然會(huì)導(dǎo)致索力分布特性發(fā)生顯著改變，但是索力個(gè)別值仍然可能處于正常值域，這時(shí)，上述評(píng)價(jià)方法就會(huì)失效。因此，基于分布特性的索力狀態(tài)評(píng)價(jià)才是更為可靠和科學(xué)的方法。但根據(jù)一根索索力分布特性的改變來(lái)識(shí)別損傷是很難做到的，因此，采用一組索索力分布特性的改變模式來(lái)識(shí)別損傷才有實(shí)現(xiàn)的可能。

以6根索的標(biāo)準(zhǔn)索力均值的改變量作為一組列向量α={w7w8w9w10w11w12}T。根據(jù)向量與向量之間的相關(guān)性來(lái)研究不同損傷工況及相同損傷工況不同損傷程度的損傷向量之間的關(guān)系。

3.1 單損傷識(shí)別

通過(guò)損傷敏感性分析，可以得到前文所述8種損傷工況的損傷向量。“模態(tài)置信準(zhǔn)則”(MAC-Modal Assurance Criterion) 通?？捎糜跈z驗(yàn)兩個(gè)向量之間的相關(guān)性：

(2)

MAC值總是在[0,1]區(qū)間內(nèi)變化，越靠近1表示α1、α2的相關(guān)性越好。但要謹(jǐn)慎區(qū)分“虛假”的MAC值。對(duì)同一損傷工況不同損傷程度之間的損傷向量進(jìn)行MAC檢驗(yàn)，其結(jié)果如表3所示。而對(duì)同一損傷程度下的不同損傷工況之間損傷向量的MAC檢驗(yàn)結(jié)果如表4所示。

表3 同一損傷工況不同損傷程度損傷向量之間的MAC值

表4 不同損傷工況相同損傷程度損傷向量之間的MAC值

表3、4表明：① 只要是同一損傷工況，不同損傷程度的損傷向量之間都具有高度相關(guān)性；② 不同損傷工況之間，即使損傷程度相同，其損傷向量之間相關(guān)性也較差。以上兩點(diǎn)顯示了通過(guò)由索力構(gòu)造的損傷向量識(shí)別損傷的可行性。對(duì)于表4中兩工況的MAC值較大“虛假”的情況，可暫把此損傷工況作為一種損傷可能。在計(jì)算損傷程度的時(shí)候，可將此“虛假”損傷排除。

基于以上損傷向量的相關(guān)性分析可知，有的損傷工況易發(fā)生混淆。為此，下面模擬工況1損傷程度5%的情況作為“未知的損傷”來(lái)闡述該文的損傷識(shí)別方法。

首先，模擬“工況1損傷程度5%”，得到損傷向量α={-1.3 -4.4 -13.5 19.9 7.5 2.4}。然后，將上述損傷向量與全部8種損傷工況的某同一損傷程度(不妨選取30%)的損傷向量求MAC值，計(jì)算結(jié)果如表5所示。

表5 未知損傷向量與已知損傷向量的MAC值

表5表明：未知損傷應(yīng)該是工況1或工況3當(dāng)中的某一種工況。下面通過(guò)識(shí)別損傷程度來(lái)排除“虛假”損傷。由損傷敏感性分析可以得出各單損傷工況的均值改變量與損傷程度趨勢(shì)線(xiàn)，并求出如圖8所示的趨勢(shì)方程。將未知損傷工況的均值改變量分別代入與之相應(yīng)的工況1與工況3的趨勢(shì)方程中，求解損傷程度。如果求解的損傷程度為0～100%且數(shù)值接近，就可以確定損傷工況及損傷程度。計(jì)算結(jié)果見(jiàn)表6。

表6 求解的損傷程度

由表6，可以將工況3排除，確定未知損傷為工況1，損傷程度為5%。用該方法進(jìn)行損傷識(shí)別，部分損傷識(shí)別結(jié)果如表7所示。

表7 損傷識(shí)別結(jié)果

3.2 組合損傷識(shí)別

組合損傷是指兩種損傷工況的組合損傷。通過(guò)數(shù)值模擬得出，組合損傷向量可近似看成這兩種單損傷向量的線(xiàn)性疊加。假設(shè)待識(shí)別的未知損傷工況是上述8種損傷工況中某兩種單損傷工況的組合。

首先，選取每種單損傷工況任意兩種已模擬的損傷程度的損傷向量(該文選取10%、30%)，并將8種單損傷工況兩兩組合(共28組)，每組一共4個(gè)損傷向量。

其次，選取第i組的4個(gè)單損傷的損傷向量[αi1αi2αi3αi4]T，i=1,2,3，…，28。假設(shè)存在正交向量βi={xi1xi2xi3xi4xi5xi6}T,使得:

[αi1αi2αi3αi4]T·βi=0

(3)

求正交向量βi(由6個(gè)未知數(shù)4個(gè)方程可知該方程有兩個(gè)線(xiàn)性無(wú)關(guān)的正交向量βi1、βi2)。以同樣的方式求解出所有組的正交向量。

最后，把未知損傷的損傷向量α的轉(zhuǎn)置αT與這28組損傷工況的正交向量βi1、βi2相乘，當(dāng)?shù)玫降臄?shù)值ai1、ai2(i=1,2,3,…,28)越接近于0，則未知損傷工況為該損傷工況的概率越大；當(dāng)?shù)玫降闹礱i1、ai2大于1時(shí)，則未知損傷工況為該損傷工況的概率基本為0。計(jì)算方程式如下：

(4)

其中：αT為未知損傷的損傷向量；βi1、βi2為第i組的正交向量。限于篇幅，該文僅提供一例加以佐證。

模擬“工況1損傷5%+工況2損傷7%”的組合損傷作為待識(shí)別的組合損傷工況，得到損傷向量。將α={-151.4 18.1 -8.1 21.6 -5.1 -20.7}代入式(4)得到的結(jié)果取絕對(duì)值，由于數(shù)值偏差較大，該文只給出[0,1]區(qū)間的結(jié)果展示，如圖9所示(圖中解1表示ai1，解2表示ai2)。

由圖9得出：未知損傷為工況1+工況2組合損傷的概率高達(dá)90%以上，其他組合損傷的概率基本為0，但不能排除未知損傷為工況1或工況2的單損傷的可能性。為識(shí)別損傷，采用3.1節(jié)中單損傷識(shí)別的方法進(jìn)行識(shí)別，結(jié)果為與工況1的MAC值為0，與工況2的MAC值為0.95。

圖9 未知損傷向量與解向量的乘積

由此可得：未知損傷與工況2的單損傷相關(guān)性較高，與工況1的相關(guān)性很差，故把工況1的單損傷情況排除。用3.1節(jié)中求單損傷程度的方法求解損傷程度，求解結(jié)果見(jiàn)表8。

表8 工況2求解的損傷程度

通過(guò)表8，可以將工況2的單損傷可能性排除，并確定損傷工況為工況1+工況2的組合損傷工況。用該方法進(jìn)行損傷識(shí)別，部分組合損傷識(shí)別結(jié)果見(jiàn)表9。

表9 組合損傷識(shí)別結(jié)果

3.3 損傷個(gè)數(shù)識(shí)別

眾所周知，在實(shí)際工程中損傷個(gè)數(shù)是無(wú)法預(yù)知的(此處損傷個(gè)數(shù)為1個(gè)或2個(gè))，損傷個(gè)數(shù)未知就無(wú)法判斷用單損傷識(shí)別還是組合損傷識(shí)別。然而，如果用組合損傷識(shí)別的方法判斷單損傷，那么包含此單損傷的所有組的正交向量βi1、βi2與未知單損傷的損傷向量轉(zhuǎn)置αT的乘積ai1、ai2都應(yīng)該接近于0。

下面以3.1節(jié)中的單損傷例子加以證明，“工況1發(fā)生5%”的損傷，得到標(biāo)準(zhǔn)索力均值改變量向量α={-1.3 -4.4 -13.5 19.9 7.5 2.4}，用組合損傷識(shí)別的方法進(jìn)行識(shí)別得到結(jié)果，如圖10所示。

圖10 未知損傷向量與正交向量的乘積

由圖10可知：包含單損傷工況1的7種組合損傷的正交向量βi1、βi2與未知損傷向量轉(zhuǎn)置αT的乘積ai1、ai2都趨近于0，這與預(yù)期一致。由此得出一套未知損傷識(shí)別的方法：在得到未知損傷向量α之后，先用組合損傷識(shí)別的方法進(jìn)行驗(yàn)證，如果只有一組正交向量βi1、βi2與未知損傷向量轉(zhuǎn)置αT的乘積ai1、ai2都趨近于0，那么就可以確定損傷工況為該組合損傷；如果有包含某一單損傷工況的所有組的正交向量βi1、βi2與未知損傷向量轉(zhuǎn)置αT乘積ai1、ai2都趨近于0，那么就可以確定損傷工況為此單損傷，并可用單損傷識(shí)別的方法加以確定，并計(jì)算出損傷程度。

4 結(jié)論

以索力構(gòu)造損傷向量是基于索力相對(duì)穩(wěn)定，易于監(jiān)測(cè)、便于應(yīng)用的特點(diǎn)，并且斜拉索本身包含了橋面和主塔結(jié)構(gòu)的空間位置信息，因此以索力作為損傷識(shí)別指標(biāo)具有可行性。通過(guò)對(duì)索力損傷向量的研究，得出以下結(jié)論：

(1) 只要是同一損傷工況，不同損傷程度的損傷向量之間都具有高度相關(guān)性；不同損傷工況之間，即使損傷程度相同，其損傷向量之間相關(guān)性也較差。由此可顯示基于索力構(gòu)建損傷向量的損傷識(shí)別具有顯著可行性。

(2) 該文方法不但可以識(shí)別單損傷情況，結(jié)合損傷程度的分析，也可以較好地識(shí)別多損傷情況；并且可以有效地區(qū)分損傷個(gè)數(shù)。