于 虹

（1.中鐵第四勘察設(shè)計(jì)院集團(tuán)有限公司,湖北 武漢 430063; 2.中國鐵建股份有限公司橋梁工程實(shí)驗(yàn)室,湖北 武漢 430063）

0 引言

斜拉橋在服役過程中，由于荷載作用、材料老化和環(huán)境因素的作用，會(huì)出現(xiàn)拉索銹蝕、主梁及主塔開裂、螺栓松動(dòng)等病害，造成結(jié)構(gòu)性能退化[1]。因此，及時(shí)對(duì)斜拉橋結(jié)構(gòu)的損傷進(jìn)行檢測和修復(fù)，對(duì)保障橋梁結(jié)構(gòu)的安全服役具有重要意義[2-3]。

斜拉索是斜拉橋的主要承重構(gòu)件之一，其受力狀態(tài)是衡量斜拉橋是否處于正常狀態(tài)的重要參量。斜拉橋結(jié)構(gòu)主梁損傷會(huì)引起恒載內(nèi)力重分布，從而引起斜拉索索力變化，因此可以利用斜拉索索力變化識(shí)別斜拉橋主梁損傷。孫宗光[4]首次探討了基于斜拉索張力指標(biāo)與人工神經(jīng)網(wǎng)絡(luò)技術(shù)相結(jié)合，對(duì)斜拉橋主梁中的損傷進(jìn)行定位識(shí)別。朱浩等[5]根據(jù)斜拉橋主梁的線形由于斜拉索索力的改變而改變這一規(guī)律，通過對(duì)主梁位移的改變量進(jìn)行反演，得到索力的變化情況，從而對(duì)斜拉橋的內(nèi)力狀態(tài)進(jìn)行識(shí)別。Hua[6]利用斜拉橋主梁損傷前后的索力變化，通過構(gòu)造斜拉橋主梁參數(shù)對(duì)斜拉索索力的靈敏度矩陣，結(jié)合模型修正技術(shù)對(duì)斜拉橋的主梁進(jìn)行損傷識(shí)別。李延強(qiáng)等[7]基于斜拉橋主梁損傷對(duì)索力的敏感性分析，選擇一組四根（空間位置對(duì)稱）對(duì)主梁損傷最敏感的斜拉索索力構(gòu)造損傷指標(biāo)，利用BP 神經(jīng)網(wǎng)絡(luò)對(duì)斜拉橋有限元模型中的主梁損傷進(jìn)行識(shí)別。

該文利用斜拉索索力作為損傷敏感指標(biāo)，并分析計(jì)算斜拉索索力對(duì)結(jié)構(gòu)損傷的靈敏度，據(jù)此對(duì)高速鐵路大跨度斜拉橋結(jié)構(gòu)中的損傷進(jìn)行識(shí)別分析。以裕溪河特大橋?yàn)檠芯繉?duì)象，對(duì)所提出方法的有效性和準(zhǔn)確性進(jìn)行驗(yàn)證。

1 斜拉索索力對(duì)主梁損傷靈敏度分析

1.1 靈敏度分析理論

靈敏度分析主要是研究系統(tǒng)輸入x的變化會(huì)導(dǎo)致系統(tǒng)響應(yīng)量g如何變化的問題。當(dāng)系統(tǒng)的輸入發(fā)生變化時(shí)，系統(tǒng)的響應(yīng)量也將會(huì)隨之發(fā)生改變，這種改變的相對(duì)程度就稱為系統(tǒng)響應(yīng)量g對(duì)于系統(tǒng)輸入變量x的靈敏度，可以用偏微分?g/x? 表示。

有限差分法的原理是利用差商近似微商，可以表示為

1.2 斜拉索索力對(duì)主梁損傷靈敏度計(jì)算

在結(jié)構(gòu)損傷識(shí)別研究中，通常假設(shè)結(jié)構(gòu)的質(zhì)量在損傷發(fā)生時(shí)不發(fā)生變化，通過結(jié)構(gòu)單元?jiǎng)偠鹊淖兓瘉砟M結(jié)構(gòu)損傷。假設(shè)斜拉橋結(jié)構(gòu)主梁第j個(gè)單元的剛度參數(shù)（如抗彎剛度）由折減至損傷狀態(tài)的，則第j個(gè)單元的剛度折減系數(shù)(Stiffness Reduction Factor，SRF)可以定義為

式中，ne——斜拉橋結(jié)構(gòu)主梁中單元的個(gè)數(shù)。剛度折減系數(shù)jγ的取值范圍為-1~0。

該文通過折減有限元模型中結(jié)構(gòu)單元的彈性模量，實(shí)現(xiàn)對(duì)結(jié)構(gòu)損傷的模擬。利用MATLAB 調(diào)用ANSYS，計(jì)算斜拉索索力對(duì)斜拉橋主梁損傷的靈敏度矩陣，具體步驟如下：

（1）計(jì)算健康狀態(tài)下裕溪河特大橋在恒載作用下的斜拉索索力向量F0。

（2）保持其他參數(shù)不變，利用MATLAB 調(diào)用ANSYS，依次令第j個(gè)主梁結(jié)構(gòu)單元發(fā)生1%剛度折減（γj=-0.01），計(jì)算此時(shí)恒載作用下的斜拉索索力向量。

（3）基于有限差分理論，利用公式（1）計(jì)算斜拉索索力對(duì)第j個(gè)主梁結(jié)構(gòu)單元損傷的靈敏度。

（4）重復(fù)步驟（2）和（3），計(jì)算斜拉索索力對(duì)所有主梁結(jié)構(gòu)單元損傷的靈敏度，得到斜拉索索力對(duì)主梁損傷的靈敏度矩陣。

2 基于斜拉索索力的損傷識(shí)別方程

基于斜拉索索力變化的損傷識(shí)別方程可以表示為

式中，[S]——斜拉索索力對(duì)結(jié)構(gòu)損傷的靈敏度矩陣；{FD}——通過測量得到的結(jié)構(gòu)損傷狀態(tài)的斜拉索索力；{F0}——通過測量得到的未損傷狀態(tài)的斜拉索索力。解向量γ中的非零元素對(duì)應(yīng)的單元編號(hào)反映損傷的位置，非零元素的大小反映損傷的程度。

3 基于稀疏正則化的結(jié)構(gòu)損傷識(shí)別

在實(shí)際工程中，通過測量獲取的斜拉索索力向量中元素的個(gè)數(shù)通常要遠(yuǎn)小于未知結(jié)構(gòu)單元?jiǎng)偠葏?shù)的個(gè)數(shù)。因此，公式（3）通常為欠定方程組，有無窮多個(gè)可能的解，不滿足適定性問題解的唯一性，即公式（3）的求解是一種不適定問題。正則化算法是求解不適定問題的有效方法，它利用具體問題的某些附加信息，通過一個(gè)近似的適定性問題去尋找原不適定問題的近似解。稀疏正則化是指在優(yōu)化問題的求解過程中，將l1范數(shù)引入目標(biāo)函數(shù)，使解向量具有稀疏性。

結(jié)構(gòu)損傷通常只發(fā)生在整個(gè)結(jié)構(gòu)的一小部分構(gòu)件中，即損傷向量具有稀疏性?；谶@個(gè)先驗(yàn)知識(shí)，可以在公式（3）的求解過程中引入對(duì)解向量γ的稀疏約束，將公式（3）轉(zhuǎn)換為一種約束優(yōu)化問題P1：

問題P1等價(jià)于如下的稀疏正則化問題：

式中，J——稀疏正則化問題的目標(biāo)函數(shù)；τ——正則化參數(shù)。

對(duì)于稀疏正則化問題，當(dāng)正則化參數(shù)的取值達(dá)到一定限值時(shí)，問題的解收斂于零，即解中元素全為零。正則化參數(shù)τ的最大值可以表示為

該研究提出的大跨斜拉橋基于索力的結(jié)構(gòu)損傷識(shí)別流程如圖1 所示。首先，選取斜拉索的索力作為損傷敏感指標(biāo)，并利用有限差分法計(jì)算斜拉索索力對(duì)結(jié)構(gòu)損傷的靈敏度；然后，建立一個(gè)由結(jié)構(gòu)剛度折減系數(shù)和損傷前后斜拉索索力組成的欠定方程組；最后，使用稀疏正則化算法對(duì)欠定方程組進(jìn)行求解，解向量中的非零元素表示結(jié)構(gòu)損傷的位置和程度。

圖1 大跨斜拉橋基于索力的結(jié)構(gòu)損傷識(shí)別流程

4 數(shù)值模擬算例

4.1 裕溪河特大橋簡介

商丘至合肥至杭州鐵路裕溪河特大橋主橋采用（61+120+324+120+61）m 雙塔鋼箱桁梁斜拉橋方案，全長686 m。以裕溪河特大橋?yàn)檠芯繉?duì)象，對(duì)所提出損傷識(shí)別方法的精度和有效性進(jìn)行驗(yàn)證。裕溪河特大橋全橋立面布置如圖2 所示。

圖2 裕溪河特大橋全橋立面布置示意圖（m）

大橋主梁為鋼箱桁結(jié)構(gòu)，桁架采用不帶豎桿的華倫式桁架，下弦采用整體式鋼箱梁，如圖3 所示。

圖3 裕溪河特大橋主梁結(jié)構(gòu)示意圖

4.2 斜拉索索力對(duì)主梁損傷靈敏度矩陣

裕溪河特大橋主梁為組合鋼箱桁梁，該文以箱梁發(fā)生損傷的情況為例，對(duì)基于索力的大跨斜拉橋結(jié)構(gòu)損傷識(shí)別過程進(jìn)行說明。裕溪河特大橋斜拉索為空間雙索面，立面上每塔兩側(cè)共13 對(duì)索，全橋共104 根斜拉索。對(duì)每根斜拉索單獨(dú)進(jìn)行編號(hào)，如圖4 所示。按照在商丘至杭州方向線路中心線的左側(cè)或右側(cè)，將全橋斜拉索分為左側(cè)索面和右側(cè)索面。對(duì)于左側(cè)索面，斜拉索自商丘至杭州方向依次編號(hào)為1~52；對(duì)于右側(cè)索面，斜拉索自商丘至杭州方向依次編號(hào)為53~104。

圖4 裕溪河特大橋斜拉索編號(hào)簡圖

將裕溪河大橋主梁箱梁按照節(jié)間劃分為57 個(gè)梁段，對(duì)每個(gè)梁段單獨(dú)進(jìn)行編號(hào)，自商丘至杭州方向依次編號(hào)為1~57，如圖5 所示。

圖5 裕溪河特大橋箱梁編號(hào)簡圖

將箱梁的每一個(gè)梁段看作一個(gè)損傷位置，依次計(jì)算每一個(gè)梁段發(fā)生1%剛度折減時(shí)的斜拉索索力，基于有限差分理論可以得到斜拉索索力對(duì)箱梁損傷的靈敏度矩陣，如圖6 所示（為方便展示，取絕對(duì)值）。

圖6 斜拉索索力對(duì)箱梁損傷的靈敏度矩陣（絕對(duì)值）

4.3 損傷工況

為了研究所提出的損傷識(shí)別方法是否對(duì)結(jié)構(gòu)中不同位置的損傷敏感，從箱梁結(jié)構(gòu)中選取具有代表性的4 個(gè)箱梁單元進(jìn)行損傷識(shí)別分析，所選取的4 個(gè)箱梁單元的編號(hào)分別為3、9、22 和28，箱梁編號(hào)如圖5 所示。一共模擬6 種箱梁結(jié)構(gòu)損傷工況，包含4 種單損傷工況和2種多損傷工況。各損傷工況具體配置如表1 所示。

表1 箱梁結(jié)構(gòu)損傷工況

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

利用損傷前后斜拉索索力變化和斜拉索索力對(duì)箱梁損傷的靈敏度，可以建立如公式（3）所示的求解方程。使用稀疏正則化算法對(duì)公式（3）進(jìn)行求解，得到各損傷工況的剛度折減系數(shù)SRF。

箱梁結(jié)構(gòu)單損傷工況的損傷識(shí)別結(jié)果如圖7 所示。從圖中可以看出，在實(shí)際損傷位置處，剛度折減系數(shù)SRF 的識(shí)別值十分明顯，而其他位置處剛度折減系數(shù)SRF 的識(shí)別值為零或者接近零。并且，3、9、22 和28號(hào)腹桿處剛度折減系數(shù)SRF 的識(shí)別值分別為-0.312、-0.307、-0.336 和-0.319，與實(shí)際值-0.3 吻合較好。以上結(jié)果表明，所提出的損傷識(shí)別方法能夠準(zhǔn)確識(shí)別裕溪河特大橋箱梁結(jié)構(gòu)中不同位置的損傷，并且識(shí)別出的損傷程度接近實(shí)際損傷程度。

圖7 單損傷工況識(shí)別結(jié)果

箱梁結(jié)構(gòu)多位置損傷的識(shí)別結(jié)果如圖8 所示。從圖中可以看出，在實(shí)際損傷位置處，剛度折減系數(shù)SRF 的識(shí)別值十分明顯，而其他位置處剛度折減系數(shù)SRF 的識(shí)別值為零或者接近零。對(duì)于工況5，3 號(hào)箱梁單元處剛度折減系數(shù)SRF 的識(shí)別值為-0.414，與實(shí)際值-0.3一致；22 號(hào)箱梁單元處剛度折減系數(shù)SRF 的識(shí)別值為-0.347，與實(shí)際值-0.3 接近。對(duì)于工況6，9 號(hào)箱梁單元處剛度折減系數(shù)SRF 的識(shí)別值為-0.336，與實(shí)際值-0.3 接近；28 號(hào)箱梁單元處剛度折減系數(shù)SRF 的識(shí)別值為-0.408，與實(shí)際值-0.3 一致。以上結(jié)果表明，所提出的損傷識(shí)別方法可以正確識(shí)別裕溪河特大橋箱梁結(jié)構(gòu)中多個(gè)位置同時(shí)發(fā)生的損傷，并且識(shí)別出的損傷程度與實(shí)際損傷程度一致。

圖8 多損傷工況識(shí)別結(jié)果

5 結(jié)論

為了識(shí)別斜拉橋結(jié)構(gòu)中的損傷，該文提出了一種基于索力的大跨斜拉橋結(jié)構(gòu)損傷識(shí)別方法。針對(duì)裕溪河特大橋數(shù)值模型的計(jì)算結(jié)果表明，所提出的損傷識(shí)別方法可以正確識(shí)別裕溪河特大橋桁架和箱梁結(jié)構(gòu)中的單位置損傷和多位置損傷，并且識(shí)別出的損傷程度與實(shí)際損傷程度一致。