張力文，何 華，陳慶志，周建庭，寧金成

(1.重慶交通大學(xué)，重慶400074;2.廣西吉泰投資有限公司，廣西南寧530001;3.河南交通職業(yè)技術(shù)學(xué)院，河南鄭州450005)

1 基本原理

在工程實(shí)踐的基礎(chǔ)上，斜拉索多階自振頻率可以更多以動(dòng)態(tài)測試及根據(jù)計(jì)算值的斜拉索設(shè)計(jì)參數(shù)利用有限元法較準(zhǔn)確地確定。設(shè)計(jì)參數(shù)與實(shí)際參數(shù)的不同，所以存在測定值與計(jì)算值的偏差［δ］。參數(shù)識(shí)別基于實(shí)際斜拉索頻率數(shù)據(jù)，就可以獲得斜拉索索力T、抗彎剛度EI、減振器剛度K、減振器位置εL，對實(shí)際狀態(tài)的斜拉橋結(jié)構(gòu)的準(zhǔn)確評價(jià)借鑒意義。這時(shí)的斜拉索參數(shù)識(shí)別問題可以被轉(zhuǎn)化為一個(gè)有約束力的最優(yōu)化問題［1－2］:

式中:E(p)為殘差向量;n為維數(shù);W為權(quán)重對角矩陣，其權(quán)重系數(shù)定義為元素wk，這是各階誤差在相應(yīng)的目標(biāo)函數(shù)中所占權(quán)重;ek(p)為斜拉索第k階自振頻率計(jì)算值和實(shí)測值之差，ek(p)==［T，EI，K，εL］(T 定義為待識(shí)別參數(shù)向量，T，EI，K，εL 參數(shù)分別用 p1，p2，p3，p4表示)。

問題的解空間分別由pL(待識(shí)別參數(shù)向量取值上限值)和pU(待識(shí)別參數(shù)向量取值下限值)構(gòu)成。用傳統(tǒng)的弦理論公式解決實(shí)際問題得到索力值T0，通過斜拉索全截面的計(jì)算得到抗彎剛度EI0，得到lgK的設(shè)計(jì)值lgK0，其中εL取拉索全長的0.01～0.1。本文中 T 和 EI取值范圍分別為［0.9T0，1.1T0］和［0.3EI0，0.45EI0］;lgK 的解空間?。?，7］。當(dāng)式(1)很小的時(shí)候，就可以認(rèn)為識(shí)別結(jié)果與各待識(shí)別參數(shù)的真實(shí)值密切接近。

對式(1)的優(yōu)化問題，可采用遺傳算法、模擬退火法、單純形法、PSO算法等多種多樣的算法，都可以獲得斜拉索的參數(shù)，如物理參數(shù)、力學(xué)參數(shù)。在此基礎(chǔ)上，筆者結(jié)合上述PSO算法對式(1)進(jìn)行優(yōu)化。采用PSO算法的識(shí)別過程如圖1。

圖1 拉索參數(shù)識(shí)別PSO優(yōu)化算法流程Fig.1 PSO optimization algorithm flow chart identified with cable parameter

1)一群體群微粒，其規(guī)模為m，初始化隨機(jī)位置和速度;本文中，第k個(gè)迭代步里，待識(shí)別參數(shù)向量為，所含4個(gè)參數(shù)向量一個(gè)為拉索索力T，一個(gè)為抗彎剛度EI，一個(gè)為減振器剛度K，還有一個(gè)為減振器位置εL。

2)為每個(gè)微粒適應(yīng)度評估;4個(gè)參數(shù)中T的取值范圍［0.9T0，1.1T0］、EI的取值范圍［0.3EI0，0.45EI0］、lgK 的取值范圍［1，7］、εL 的取值范圍［0.01L，0.1L］。

3)將每個(gè)微粒經(jīng)歷過的pbest和適應(yīng)值做比較，將好的那個(gè)值取代原值為pbest。

4)將每個(gè)微粒經(jīng)歷過的gbest和適應(yīng)值做比較，將好的那個(gè)值取代原值為gbest。

5)由式(1)改變微粒位置和速度。

6)最后得到minj(p)，與一期望足夠小的數(shù)ε做比較，當(dāng)ε較大時(shí)，輸出優(yōu)化參數(shù)的結(jié)果;否則，返回第2)重復(fù)循環(huán)。

2 算法應(yīng)用

通過對比實(shí)際建模時(shí)由有限元法測得的和基于PSO優(yōu)化算法得到的拉索參數(shù)，驗(yàn)證本文所結(jié)合的PSO算法的正確性。索所需參數(shù)示意圖如圖2，參數(shù)如表 1。對于拉索 B1，取 T=2 981 kN，EI=5.811 78 ×105N·m2，K=3 ×105N/m，εL=2.374 m，拉索前9階頻率和拉索B2的各階頻率可以由有限元方法建立的拉索實(shí)際模型得到(表2)。

圖2 拉索參數(shù)示意Fig.2 Cable parameter diagram

表1 拉索基本參數(shù)Table 1 Basic cable parameters

表2 拉索的各階頻率Table 2 Each frequency order of cable/Hz

通過MATLAB的程序計(jì)算，結(jié)合PSO優(yōu)化得到的參數(shù)識(shí)別結(jié)果如圖3。

表2中的頻率現(xiàn)假設(shè)成初始獲得的拉索參數(shù)，可清楚的看到對比B1、B2號(hào)拉索每一個(gè)參數(shù)初始值與識(shí)別結(jié)果(表3)。

表3 各參數(shù)識(shí)別結(jié)果和設(shè)計(jì)初值的對比Table 3 The contrast of the recognition results and the initial value of design of each parameter

利用前2，3，4，…，10階頻率，結(jié)合PSO優(yōu)化算法與有限元方法對比頻率階數(shù)對于識(shí)別結(jié)果產(chǎn)生的影響，由此得到了索力T、抗彎剛度EI、減振器等效剛度lgK、減振器位置εL識(shí)別結(jié)果的相對誤差如圖4。

圖3 PSO優(yōu)化計(jì)算結(jié)果Fig.3 PSO optimization calculation results

圖4 不同頻率階次對識(shí)別結(jié)果的影響Fig.4 The recognition results affected by different frequency order

由表2及圖4可知:兩種算法所求得的值偏差是比較小的，T，EI，lgK，εL的最大識(shí)別誤差分別為－0.932%，－2.347%，1.521%，4.701%。與此同時(shí)可得到以下結(jié)論:識(shí)別結(jié)果與所采用的頻率階次正相關(guān)，為提高振動(dòng)法索力測試的準(zhǔn)確性，可以在實(shí)際運(yùn)用過程中測出斜拉索的前10階自振頻率，這可以通過環(huán)境激勵(lì)等方法達(dá)到。

3 實(shí)例求證

干溪溝1號(hào)特大橋?yàn)橹鳂驗(yàn)?55 m+360 m+155 m雙塔雙索面PC斜拉橋，位于渝湘高速公路黔彭段上，橋梁跨越深切割的干溪溝峽谷，為全漂浮體系。做為渝湘高速公路上的“三最”斜拉橋(索塔的高度最高、主跨跨度最大、橋面距溝底最深)，干溪溝1號(hào)特大橋全橋布有180根斜拉索，拉索采用112套橡膠組合減振器，為采購定型產(chǎn)品，其中A9～A22、B9～B22號(hào)斜拉索采用體外減振器。

分別對 3 號(hào)塔斜拉索 A4，A6，A7，A8，A9 進(jìn)行現(xiàn)場測試，其基本參數(shù)如表4。對此5根拉索進(jìn)行動(dòng)力測試，其中A4號(hào)索(帶有橡膠減振器)所測試的動(dòng)力響應(yīng)如圖5、圖6。

圖5 A4號(hào)斜拉索加速度時(shí)程響應(yīng)Fig.5 Acceleration time-history response figure of No.A4 stayed-cables

圖6 A4號(hào)斜拉索頻譜分析Fig.6 Spectrum analysis figure of No.A4 stayed-cables

表4 斜拉索基本參數(shù)Table 4 Basic stayed-cables parameters

表5 各拉索前5階頻率值Table 5 Top-five frequency order value for each cable/Hz

表5為各拉索前5階頻率值，利用本文PSO優(yōu)化算法對拉索索力T、抗彎剛度EI、減振器剛度K、減振器的位置εL進(jìn)行精細(xì)識(shí)別，識(shí)別結(jié)果及誤差如表6。

4 結(jié)論

表6 拉索參數(shù)識(shí)別結(jié)果及誤差Table 6 Recognition results and errors of cable parameter

對于加減振器的斜拉索，利用PSO識(shí)別算法所得的索力與應(yīng)用有限元法所求得值偏差是很小的，因此筆者基于拉索實(shí)測頻率之間復(fù)雜的非線性關(guān)系，針對拉索參數(shù)進(jìn)行識(shí)別采用PSO優(yōu)化算法能夠得到準(zhǔn)確的和實(shí)用的參數(shù)。本研究對今后的斜拉索參數(shù)識(shí)別有一定的借鑒意義。

［1］張清華，卜一之，李喬.基于非線性振動(dòng)的斜拉索參數(shù)識(shí)別無模型方法［J］.振動(dòng)工程學(xué)報(bào)，2008，21(6):572 －575.Zhang Qinghua，Bu Yizhi，Li Qiao.Based on parameter identification of nonlinear vibration of cable-stayed model-free method ［J］.Vibration Engineering，2008，21(6):572 －575.

［2］張清華，冉再紅，卜一之，等.拉索非線性振動(dòng)問題求解及參數(shù)識(shí)別方法研究［J］.土木工程學(xué)報(bào)，2009，42(6):86 －91.Zhang Qinghua，Ran Zaihong，Bu Yizhi，et al.Nonlinear vibration problem solving and parameter identification method ［J］.Journal of Civil Engineering，2009，42(6):86－91.

［3］陳杰.MATLAB寶典［M］.北京:電子工業(yè)出版社，2010:22－35.

［4］Chen Qingzhi，Yan Banfu，Wang Guangwu.Vibration-based cable force measurement of inclined cable with an attached vibration absorber［C］.Changsha，P.R.China:The 4th International Symposium on Lifetime Engineering of Civil Infrastructure，2009:26 －28.

［5］曾建潮，崔志華.一種保證全局收斂的PSO算法［J］.計(jì)算機(jī)研究與發(fā)展，2004，41(8):1333 －1338.Zeng Jianchao，Cui Zhihua.A guarantee of global convergence of the PSO algorithm ［J］.Computer Research and Development，2004，41(8):1333－1338.

［6］孫越泓，魏建香，夏德深.一種基于粒子對稱分布多樣性的PSO算法［J］.模式識(shí)別與人工智能，2010，23(2):137 －142.Sun Yuehong，Wei Jianxiang，Xia Deshen.Distribution based on particle symmetry diversity of PSO algorithm ［J］.Pattern Recognition and Artificial Intelligence，2010，23(2):137 －142.

［7］劉雨青.橋梁結(jié)構(gòu)模態(tài)參數(shù)識(shí)別與應(yīng)用研究［D］.武漢:武漢理工大學(xué)，2005.

［8］周云，易偉建.帶減振圈斜拉索的參數(shù)識(shí)別及損傷診斷研究［J］.中外公路，2007，27(1):99 －103.Zhou Yun，Yi Weijian.Damping ring with a cable-stayed parameter identification and damage diagnosis［J］.Journal of China ＆ Foreign Highway，2007，27(1):99－103.