孫 浩

(湖南省高速公路集團(tuán)有限公司, 湖南 長(zhǎng)沙 410026)

近年來(lái)，隨著我國(guó)科技力量不斷發(fā)展，新的橋梁結(jié)構(gòu)體系層出不窮，鋼-混疊合梁組合結(jié)構(gòu)是其中一種具有代表性的橋型，一般通過(guò)剪力連接件使鋼結(jié)構(gòu)主梁與混凝土橋面形成整體。與混凝土結(jié)構(gòu)相比，可大幅減輕結(jié)構(gòu)自重，有效提高橋梁跨越能力；與鋼結(jié)構(gòu)相比，混凝土橋面板可解決鋼結(jié)構(gòu)箱梁頂板鋪設(shè)施工的施工難題。此外鋼-混組合結(jié)構(gòu)箱梁整體性好，抗扭能力強(qiáng)，具有廣闊的發(fā)展前景。目前，我國(guó)已經(jīng)針對(duì)鋼-混組合結(jié)構(gòu)開(kāi)展了多項(xiàng)研究，研究方向主要為結(jié)構(gòu)整體受力及局部剪力連接件強(qiáng)度及抗疲勞工作性能，而對(duì)有限元模型修正方向研究較少。考慮到有限元理論模型與實(shí)際結(jié)構(gòu)之間不可避免地存在誤差，因此很有必要對(duì)其開(kāi)展有限元修正方法研究，本文以某鋼-混疊合梁為研究對(duì)象，基于MPGA算法對(duì)其進(jìn)行有限元模型修正，并通過(guò)現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)予以驗(yàn)證，研究成果可為鋼-混組合體系有限元模擬提供一定參考。

1 工程概況

某兩跨鋼-混疊合梁，橋跨布置為(38+42)m，主梁結(jié)構(gòu)形式為雙箱單室等截面鋼箱梁，材料為Q345D鋼材，橋面采用C40混凝土。梁高1.95 m，頂板寬12 m，頂板混凝土厚度為22 cm。橫截面分為左右兩幅，中間通過(guò)橫梁連接，為加強(qiáng)箱梁橫向剛度，約束箱梁扭轉(zhuǎn)及翹曲變形，在全橋范圍內(nèi)每隔5 m設(shè)置一橫隔板，板厚20 mm，橫向設(shè)置4個(gè)腹板，腹板厚25 mm，混凝土橋面板與鋼箱梁之間通過(guò)剪力連接件(剪力釘)連接，其中墩頂剪力釘布置間距為15 cm，其余位置25 cm，采用φ22 mm×200 mm圓柱頭焊釘。為獲取橋梁結(jié)構(gòu)的應(yīng)力變化及振動(dòng)頻率，在每個(gè)吊裝節(jié)段選取一典型截面布置應(yīng)變測(cè)點(diǎn)及加速度計(jì)，根據(jù)動(dòng)載試驗(yàn)結(jié)果拾取結(jié)構(gòu)加速度響應(yīng)值。橫截面如圖1所示，測(cè)點(diǎn)布置見(jiàn)圖2、圖3。

2 初始有限元模型建立

根據(jù)結(jié)構(gòu)實(shí)際材料特性(見(jiàn)表1)，通過(guò)有限元軟件ANSYS APDL建立該橋初始參數(shù)化有限元模型，考慮到結(jié)構(gòu)的空間受力特性，使用Solid45實(shí)體單元模擬混凝土橋面板，鋼結(jié)構(gòu)主梁采用Shell63殼單元模擬，鋼箱梁與混凝土橋面板的連接使用彈簧單元Combin14模擬，為精準(zhǔn)模擬剪力釘?shù)目辜艉涂估涡阅?，在Combin14單元實(shí)常數(shù)中賦予X、Y、Z這3個(gè)方向的剛度以模擬其三向受力狀態(tài)。邊界條件按照橋梁實(shí)際支座布置施加對(duì)應(yīng)約束，采用六面體映射法劃分單元網(wǎng)格[1-2]，有限元模型見(jiàn)圖4、圖5。

圖1 標(biāo)準(zhǔn)節(jié)段截面示意圖(單位： mm)

圖2 應(yīng)變測(cè)點(diǎn)布置圖

圖3 動(dòng)力試驗(yàn)測(cè)試測(cè)點(diǎn)布置圖

表1 主要材料特性構(gòu)件材料彈性模量/MPa容重/(kN·m-3)混凝土板C403.25×10425.0鋼箱梁Q3452.1×10578.5剪力釘ML15AL2.0×10578.5

圖4 橋梁整體及鋼結(jié)構(gòu)箱梁有限元模型

圖5 剪力連接件有限元模型

3 基于MPGA算法的有限元模型修正方法

MPGA算法原理為首先設(shè)置多個(gè)初始子種群，對(duì)各個(gè)控制參數(shù)予以賦值，確定種群的進(jìn)化方向，同時(shí)引入遷移算子以保證變異后子種群數(shù)據(jù)交換的便捷性和精確性，選取第(N-1)個(gè)母體變異產(chǎn)生下代種群的(N-1)個(gè)子體，引入精英策略篩選子種群適應(yīng)度最高的個(gè)體作為修正結(jié)果，由該法選出的最優(yōu)個(gè)體可有效防止優(yōu)秀基因流失，提升優(yōu)化計(jì)算精度。MPGA算法修正結(jié)果能較好地滿足適應(yīng)度的單調(diào)性，確保修正過(guò)程始終處于全局最優(yōu)收斂狀態(tài)[3-5]。

1) 目標(biāo)函數(shù)：通過(guò)對(duì)敏感度較高的參數(shù)進(jìn)行修正，使有限元理論計(jì)算結(jié)果與實(shí)際結(jié)構(gòu)受力狀態(tài)更為接近，以結(jié)構(gòu)應(yīng)力指標(biāo)作為控制指標(biāo)，目標(biāo)函數(shù)擬取為：

(1)

式中：f(σ)為應(yīng)力目標(biāo)函數(shù)；σ1為理論計(jì)算應(yīng)力值；σ2為現(xiàn)場(chǎng)實(shí)測(cè)應(yīng)力值。

2) 設(shè)計(jì)變量：選取上部結(jié)構(gòu)鋼箱梁彈性模量、鋼箱梁質(zhì)量密度、混凝土橋面彈性模量、混凝土質(zhì)量密度、剪力釘彈性模量等5個(gè)參數(shù)作為待修正參數(shù)，即設(shè)計(jì)變量。

根據(jù)以上參數(shù)設(shè)置，初始種群設(shè)置為20×10，以MAC≥0.9為約束條件，初始種群長(zhǎng)度設(shè)置為20，交叉率取0.6～0.9，變異率取0.001～0.005，子種群最優(yōu)個(gè)體保持代數(shù)設(shè)置為30，收斂容差設(shè)置為e-3。使用Matlab編制數(shù)學(xué)優(yōu)化程序，將程序指令轉(zhuǎn)換為inp輸入文件，將待修正參數(shù)以APDL語(yǔ)言的模式輸入至相應(yīng)求解模塊中，可極大節(jié)約計(jì)算時(shí)間，提升計(jì)算效率，避免復(fù)雜的前后處理操作[6]。具體算法步驟如下：

1) 在Matlab中設(shè)置相應(yīng)參數(shù)生成初始種群，對(duì)初始種群實(shí)施變異生成下一代子種群，賦予目標(biāo)函數(shù)、指定待修正參數(shù)、設(shè)置收斂準(zhǔn)則后進(jìn)行全局最優(yōu)搜索求解。

2) 對(duì)于每次迭代結(jié)果進(jìn)行收斂性驗(yàn)證，若滿足收斂條件，則認(rèn)為該解即為全局最優(yōu)結(jié)果，若不滿足收斂性準(zhǔn)則，則繼續(xù)變異至下一代子種群重新進(jìn)行全局最優(yōu)求解，每次變異均應(yīng)對(duì)種群適應(yīng)度進(jìn)行檢驗(yàn)，若適應(yīng)度滿足要求，則終止變異。

3) 重復(fù)以上步驟，直至收斂，輸出計(jì)算結(jié)果。

4 修正前后對(duì)比結(jié)果

通過(guò)Matlab與ANSYS APDL結(jié)合，在一定可行域范圍內(nèi)對(duì)待修正參數(shù)進(jìn)行修正，多次迭代優(yōu)化后，可獲得目標(biāo)函數(shù)最小時(shí)的設(shè)計(jì)變量最優(yōu)解，具體如表2所示。

表2 基于MPGA算法修正后設(shè)計(jì)參量對(duì)比結(jié)果參數(shù)容重彈性模量初始值/(e4MPa)修正值/(e4MPa)修正幅度/%初始值/(kN·m-3)修正值/(kN·m-3)修正幅度/%鋼箱梁78.579.961.862.102.299.05橋面混凝土25.026.084.323.253.9421.23剪力釘———2.002.063.00

修正結(jié)果表明：通過(guò)MPGA算法修正后，5個(gè)待修正參數(shù)值均有一定幅度增大，其中橋面混凝土彈性模量變化幅度最大，達(dá)21.23%，鋼箱梁彈性模量變化次之，達(dá)9.05%，這是因?yàn)榛炷翉椥阅Ａ勘旧砭哂休^強(qiáng)離散性，加上橋面板中布設(shè)有大量鋼筋，有限元建模時(shí)并未考慮普通鋼筋的影響，故修正后混凝土折合彈性模量大幅提高。同時(shí)由于剪力釘?shù)倪B接作用，使混凝土橋面板與鋼箱梁形成整體，二者在剪力釘連接位置局部剛度主要受剪力釘本身彈性模量變化影響，故修正后鋼箱梁彈性模量計(jì)入了因剪力釘彈性模量變化造成的影響。

自振頻率能直接反映結(jié)構(gòu)剛度分布，因此可利用自振頻率驗(yàn)證參數(shù)修正后有限元模型與實(shí)際結(jié)構(gòu)的匹配程度，模態(tài)階次越高，越能驗(yàn)證有限元模型的正確性?？紤]到理論計(jì)算振型與實(shí)測(cè)振型對(duì)應(yīng)階次不一定完全吻合，引入MAC評(píng)價(jià)標(biāo)準(zhǔn)對(duì)計(jì)算與實(shí)測(cè)振動(dòng)頻率進(jìn)行振型配對(duì)[7]，具體描述為以下形式：

(2)

其中，φn,φt分別表示計(jì)算與實(shí)測(cè)振型向量，MAC值一般介于0～1之間，其值越接近于1，表明二者匹配程度越高，相關(guān)性越好。本文以0.9作為匹配程度界限值，當(dāng)MAC≥0.9時(shí)，即認(rèn)為二者對(duì)應(yīng)統(tǒng)一振動(dòng)頻率。

提取修正前后結(jié)構(gòu)自振頻率結(jié)果，并與實(shí)際測(cè)試結(jié)果進(jìn)行對(duì)比，如表3所示。

基于MPGA算法進(jìn)行參數(shù)修正后，理論計(jì)算頻率與實(shí)測(cè)頻率更為吻合，誤差均在10%以內(nèi)，MAC值均大于0.9，特別是5～10階高階振動(dòng)頻率，與實(shí)測(cè)頻率誤差均在5%以內(nèi)，MAC值較修正前更趨于1?？烧J(rèn)為有限元模型修正取得了良好效果，修正后有限元模型能較為精確地反映實(shí)際結(jié)構(gòu)的受力狀態(tài)及動(dòng)力特性，可作為基準(zhǔn)有限元模型，為運(yùn)營(yíng)階段疊合梁長(zhǎng)期健康監(jiān)測(cè)提供理論依據(jù)。

表3 修正前后結(jié)構(gòu)動(dòng)力特性對(duì)比結(jié)果振型階數(shù)振型特征實(shí)測(cè)頻率/Hz計(jì)算值/Hz誤差/%MAC值修正前修正后修正前修正后修正前修正后1一階豎彎1.131.011.06-10.62-6.190.830.912二階豎彎2.262.092.17-7.52-3.980.850.933三階豎彎3.843.363.52-12.50-8.330.830.904一階橫彎4.634.314.49-6.91-3.020.850.945二階橫彎5.515.045.28-8.53-4.170.840.936四階豎彎5.985.295.73-11.54-4.180.820.927三階橫彎7.236.847.02-5.39-2.900.870.968五階豎彎7.867.137.55-9.29-3.940.820.949六階豎彎8.147.657.88-6.02-3.190.880.9310一階扭轉(zhuǎn)8.928.088.49-9.42-4.820.820.94

5 結(jié)語(yǔ)

基于MPGA算法原理，以某鋼-混疊合梁為研究對(duì)象，聯(lián)合Matlab與ANSYS對(duì)其進(jìn)行有限元模型修正，可得到以下結(jié)論：

1)通過(guò)構(gòu)建合理修正參數(shù)，獲得的有限元理論計(jì)算結(jié)果能較為全面地反映橋梁結(jié)構(gòu)的動(dòng)力特性，修正效果比基于單一靜載或動(dòng)載試驗(yàn)的模型修正效果更好。

2)通過(guò)有限元模型修正，選取的待修正參數(shù)數(shù)值均有不同幅度增大，其中參數(shù)值變化最大的為橋面板混凝土彈性模量、其次為鋼箱梁彈性模量，修正幅度分別為21.23%和9.05%。

3)有限元模型修正后，結(jié)構(gòu)前10階自振頻率與實(shí)測(cè)值更為吻合，最大僅為8.33%，相比于修正前的12.5%有一定程度降低；對(duì)于5～10階高階頻率，修正后理論值與實(shí)測(cè)值更為接近，MAC值也更趨于1，說(shuō)明修正后模型能更好地模擬實(shí)際結(jié)構(gòu)的剛度分布，較好地反映實(shí)際結(jié)構(gòu)的受力狀態(tài)。