劉 洋,王 楠

(1.武漢港灣工程質(zhì)量檢測(cè)有限公司，武漢 430040；2.海工結(jié)構(gòu)新材料及維護(hù)加固技術(shù)湖北省重點(diǎn)實(shí)驗(yàn)室，武漢 430040；3.中交武漢港灣工程設(shè)計(jì)研究院有限公司，武漢 430040)

現(xiàn)代鋼拱橋建設(shè)在我國(guó)起步較晚，理論研究也落后于國(guó)外[1,2]，得益于近幾十年來國(guó)內(nèi)大跨度拱橋的大量修建，其跨徑不斷增加，矢跨比提高，相對(duì)剛度減弱，可見高聳的桁架結(jié)構(gòu)受側(cè)向風(fēng)荷載作用不容忽略[3,4]，尤其在主拱懸拼過程中，側(cè)向剛度較小，對(duì)側(cè)向風(fēng)荷載更為敏感，因此不得不關(guān)注鋼桁架拱橋的抗風(fēng)穩(wěn)定性。目前國(guó)內(nèi)對(duì)鋼桁架拱橋穩(wěn)定性的研究重點(diǎn)主要集中在恒載、施工臨時(shí)荷載和汽車活載等作用下的第一類穩(wěn)定性[5]，對(duì)風(fēng)荷載的分析仍然以影響因素的形式加以考量，較少關(guān)注施工過程中的抗風(fēng)穩(wěn)定性，尤其是極值點(diǎn)失穩(wěn)問題，為此，對(duì)鋼桁架拱橋施工過程中的彈塑性穩(wěn)定分析很有必要。

1 極值點(diǎn)失穩(wěn)問題的有限元分析

極值點(diǎn)失穩(wěn)問題的分析過程可以理解為同時(shí)考慮幾何非線性和材料非線性的極限承載力問題，結(jié)構(gòu)達(dá)到極限承載能力時(shí)，表現(xiàn)為結(jié)構(gòu)完全崩潰，有限元[1,6]中，即為考慮幾何非線性和材料非線性對(duì)結(jié)構(gòu)剛度的影響，通過求解非線性平衡方程，得出結(jié)構(gòu)的極限承載力，求解方法通常采用荷載增量法[6-8]。結(jié)構(gòu)在外荷載不斷增加的過程中，結(jié)構(gòu)剛度也不斷呈非線性變化，當(dāng)材料進(jìn)入塑性狀態(tài)后，材料非線性對(duì)結(jié)構(gòu)剛度也會(huì)產(chǎn)生影響，隨著結(jié)構(gòu)剛度矩陣趨于奇異，結(jié)構(gòu)承載能力達(dá)到極限，在整個(gè)過程中，可用增量方程表示為[6]

(0[K]o+0[K]L+0[K]σ){Δu}i={ΔP}i

(1)

其中，{ΔP}i為第i步荷載增量，{Δu}i為荷載增量{ΔP}i導(dǎo)致的節(jié)點(diǎn)位移增量。在求解極限承載力時(shí)，有如下兩個(gè)問題需要注意。

1.1 荷載增量的選取

在程序計(jì)算過程中，荷載是被分成若干步荷載的形式施加到結(jié)構(gòu)上的，在計(jì)算開始時(shí)，位移近似呈線性變化，隨著荷載的不斷加大，初應(yīng)力剛度影響越來越顯著，切線剛度矩陣和單元?jiǎng)偠染仃囆辛惺降闹抵饾u減小，荷載位移曲線趨于平穩(wěn)，因此在荷載增量選取時(shí)可以在計(jì)算開始時(shí)適當(dāng)取得大點(diǎn)，節(jié)約計(jì)算時(shí)間，同時(shí)，在將要達(dá)到極限承載能力時(shí)，較大的荷載增量使得計(jì)算結(jié)果難以收斂，因此需要逐步減小荷載增量，在程序分析時(shí)，可以采用自動(dòng)荷載步的方式。ANSYS有限元軟件實(shí)現(xiàn)自動(dòng)荷載步的方式主要是打開自動(dòng)步長(zhǎng)，根據(jù)前一荷載步荷載在分析計(jì)算時(shí)的迭代次數(shù)來確定當(dāng)前荷載步的荷載增量，若前一步荷載在設(shè)定的迭代次數(shù)內(nèi)未能收斂，則采用二分的形式重新劃分荷載增量并繼續(xù)計(jì)算。另外選取荷載增量的方式還可以依據(jù)特征剛度法[9]，表示為

[Ki]={ΔP}i/Δui

(2)

其中[Ki]表示第i步結(jié)構(gòu)的特征剛度，{ΔP}i為第i步荷載增量，Δui為第i步結(jié)構(gòu)位移中最顯著的一個(gè)節(jié)點(diǎn)增量位移分量，第i+1步荷載增量的計(jì)算依據(jù)第i步結(jié)構(gòu)的特征剛度而選取，表示為

{ΔP}i+1/{ΔP}i=[Ki+1]/[Ki]

(3)

由此可以確保荷載增量開始較大，隨著結(jié)構(gòu)趨于失穩(wěn)時(shí)，荷載增量逐漸減小，由于[Ki+1]是未知的，則下一荷載步的荷載增量可以改寫為

{ΔP}i+1={ΔP}i·[Ki]/[Ki-1]

(4)

1.2 位移增量法

荷載增量法僅適用于極值點(diǎn)前的屈曲分析，在極值點(diǎn)處，結(jié)構(gòu)的切線剛度矩陣的行列式值為零，因此無論荷載增量多小，結(jié)構(gòu)都不可能發(fā)生相應(yīng)的位移增量以達(dá)到平衡，也就是說計(jì)算不能收斂，為此荷載增量法無法逾越極值點(diǎn)，為進(jìn)行后屈曲分析，我們可以在極值點(diǎn)處，采用位移增量法以替代荷載增量法。一般而言，可以將剛度矩陣重新排列，使得要控制的位移列到最后一項(xiàng)，對(duì)原剛度矩陣進(jìn)行分塊，增量平衡方程可表示為

(5)

式中,Δλ為控制荷載增量步長(zhǎng)的系數(shù)，(R1,R2)T為求解迭代過程中的不平衡力向量。上式可進(jìn)一步改寫為

(6)

在極值點(diǎn)處迭代時(shí)，通過給定位移增量Δu2的值，可以求出相應(yīng)的位移增量Δu1和荷載增量步長(zhǎng)的系數(shù)Δλ，由于位移增量影響結(jié)構(gòu)的剛度矩陣，求解過程中需要反復(fù)迭代運(yùn)算，當(dāng)(R1,R2)T的值在收斂容差內(nèi)，則計(jì)算收斂。

2 工程背景

2.1 項(xiàng)目簡(jiǎn)介

此項(xiàng)目北起于規(guī)劃永安南路，自北向南布置，經(jīng)規(guī)劃潘臺(tái)路、航空路、規(guī)劃路、唐白河、濱河路，后止于規(guī)劃金沙路，主線全長(zhǎng)1.945 km。主線橋起訖點(diǎn)樁號(hào)為K0+200～K2+145，其中橋梁起訖點(diǎn)樁號(hào)為K0+340.5～K1+913.5，全長(zhǎng)1.573 km，主橋?yàn)榭鐝浇M合(70+240+70) m的下承式鋼桁架拱橋，橋面全寬43.9 m。

全橋橋型布置圖如圖1所示。

2.2 建立計(jì)算模型

該文采用ANSYS有限元軟件進(jìn)行分析計(jì)算，計(jì)算模型中鋼桁桿件、橫縱梁和臨時(shí)吊索塔架構(gòu)件均采用BEAM189單元進(jìn)行模擬，吊桿、扣索、背索和臨時(shí)系桿均采用LINK10單元模擬，橋面板采用SHELL63單元模擬，實(shí)常數(shù)根據(jù)實(shí)際橋面板厚度設(shè)置，吊索塔架與主體結(jié)構(gòu)進(jìn)行鉸接，通過不耦合相應(yīng)兩節(jié)點(diǎn)繞橫橋向的轉(zhuǎn)動(dòng)自由度加以實(shí)現(xiàn)。主體結(jié)構(gòu)及臨時(shí)施工構(gòu)件的幾何參數(shù)按照設(shè)計(jì)圖紙取值，鋼材的彈模取為206 GPa，拉索的彈模取為195 GPa，為實(shí)現(xiàn)材料非線性，鋼材采用雙線性隨動(dòng)強(qiáng)化(BKIN)的理想彈塑性材料，屈服準(zhǔn)則采用Von Mises屈服準(zhǔn)則。主橋鋼結(jié)構(gòu)除部分中支點(diǎn)處下弦桿采用Q390D鋼材外其余均采用Q345D鋼材；橋面板采用C50的混凝土；吊桿為強(qiáng)度等級(jí)為1 770 MPa的φ7成品高強(qiáng)鋼絲，規(guī)格為PES7-109；柔性系桿采用強(qiáng)度等級(jí)為1 860 MPa的56φ15.24鋼絞線?？鬯⒅蜋M聯(lián)采用Q345B的鋼材外，其余立柱剪力撐均采用Q235B鋼材；扣索和背索采用φ15.24低松弛鋼絞線，每桁均由2組140φ15.24鋼絞線組成。

鋼桁拱橋考慮施工階段的有限元模型如圖2所示，共有個(gè)13 567節(jié)點(diǎn)，共10 732個(gè)單元。另外，成橋狀態(tài)下拱橋縱梁邊界條件布置圖如圖3所示。

2.3 結(jié)果分析

穩(wěn)定性分析也即屈曲分析，屈曲分析是一種用于確定結(jié)構(gòu)在外荷載作用下開始變得不穩(wěn)定時(shí)的臨界荷載值和屈曲結(jié)構(gòu)發(fā)生屈曲響應(yīng)時(shí)的模態(tài)形狀的技術(shù)。屈曲分析包括線彈性屈曲分析和非線性屈曲分析，線彈性屈曲分析在ANSYS軟件中對(duì)應(yīng)特征值屈曲分析，非線性屈曲分析也即考慮幾何非線性和材料非線性雙重非線性的極限承載力分析。特征值屈曲分析的理論前提是材料始終處于線彈性階段，結(jié)構(gòu)變形處于小撓度范圍內(nèi)，只適用于理想結(jié)構(gòu)，求解結(jié)果偏安全，但也有其工程意義，此外，其求解結(jié)果為極限承載力分析提供了荷載理論上限，對(duì)后續(xù)非線性屈曲分析的荷載取值具有指導(dǎo)意義，因此對(duì)結(jié)構(gòu)進(jìn)行線彈性穩(wěn)定分析也很有必要。

2.3.1 一階穩(wěn)定性在ANSYS中分析方法

特征值屈曲分析求解結(jié)果包括λ和{u}，即為荷載系數(shù)和屈曲模態(tài)，在軟件中實(shí)施的具體步驟如下：

1)靜力分析

進(jìn)入靜力分析求解類型，打開預(yù)應(yīng)力效應(yīng)，PSTRES,ON，施加荷載進(jìn)行求解；

2)模態(tài)分析

重新進(jìn)入求解模塊，選擇模態(tài)分析類型，指定屈曲分析選項(xiàng)，求解完成，進(jìn)入后處理模式。

2.3.2 恒載及安裝荷載下的一階穩(wěn)定性

施工過程中除了主體結(jié)構(gòu)及臨時(shí)結(jié)構(gòu)自重和拉索內(nèi)力外，中跨拱肋合龍前還存在邊跨壓重，通過ANSYS有限元軟件對(duì)鋼桁拱橋施工過程中4個(gè)特定工況在恒載和安裝荷載作用時(shí)的彈性穩(wěn)定性進(jìn)行分析，得出如下結(jié)果，見表1。

表1 各工況下恒載和安裝荷載一階穩(wěn)定系數(shù)

注：表中穩(wěn)定系數(shù)λ=極限荷載/單倍荷載。

考慮恒載和安裝荷載作用各工況下失穩(wěn)模態(tài)如圖4所示。

從分析結(jié)果可以看出，吊索扣塔安裝前穩(wěn)定安全系數(shù)較高，失穩(wěn)形態(tài)主要表現(xiàn)為桁架腹桿局部失穩(wěn)；安裝吊索扣塔后，拆除邊跨的臨時(shí)支架，隨著中跨節(jié)段的拼裝，荷載不斷加大，此時(shí)結(jié)構(gòu)失穩(wěn)表現(xiàn)為吊索塔架面外壓彎屈曲；隨著中跨荷載不斷加大，穩(wěn)定安全系數(shù)不斷減小，安裝至主梁合龍段時(shí)，結(jié)構(gòu)的失穩(wěn)安全系數(shù)降至10.26；吊索塔架拆除后，卸載邊跨壓重，成橋狀態(tài)在恒載作用下的安全系數(shù)達(dá)到9.81，失穩(wěn)形態(tài)表現(xiàn)為主拱面外反對(duì)稱失穩(wěn)?？傊?，在恒載和安裝荷載作用下，施工過程中及成橋狀態(tài)的穩(wěn)定性系數(shù)均較高。

3 結(jié) 論

a.靜風(fēng)荷載一階穩(wěn)定系數(shù)值均較大，可以看出施工過程及成橋狀態(tài)時(shí)整體抗風(fēng)性能良好，結(jié)構(gòu)抗風(fēng)可靠度較高；

b.各工況下失穩(wěn)模態(tài)均表現(xiàn)為第一道風(fēng)撐平聯(lián)桿件面外的局部失穩(wěn)，對(duì)于高聳的扣塔而言抵抗側(cè)向風(fēng)安全性能良好，拱肋合龍前，隨著懸臂長(zhǎng)度不斷增加，拱肋的橫向風(fēng)荷載穩(wěn)定安全系數(shù)不斷減小，隨著拱肋和主梁的合龍，側(cè)向剛度增大，穩(wěn)定系數(shù)略有增大。

[1] 姜禮尚．有限元方法及其理論基礎(chǔ)[M]．北京：人民交通出版社，1980．

[2] 梅建松．大跨度鋼拱橋穩(wěn)定性及極限承載力的研究[D]．武漢：武漢理工大學(xué)，2005．

[3] 李國(guó)豪．橋梁結(jié)構(gòu)穩(wěn)定與振動(dòng)[M]．北京：中國(guó)鐵道出版社，2003．

[4] 王勖成．有限單元法[M]．北京.清華大學(xué)出版社，2003．

[5] Bathe K J,Bolourchi.Large Displacement Analysis of Three Dimensional Beam Structures[J].International Journal for Numerical Methods in Engineering,1979,14:961～986.

[6] 張洪偉．ANSYS非線性有限元分析方法及范例應(yīng)用[M]．北京：中國(guó)水利水電出版社，2013．

[7] Ever J Barbero.Finite Element Analysis of Composite Materials[M].Boca Raton:CRC Press,2007.

[8] 蓋衛(wèi)明．大跨度鋼桁拱橋的穩(wěn)定與極限承載力研究[D]．長(zhǎng)沙：中南大學(xué)，2009.

[9] 劉 劍，劉永健，楊炳成．關(guān)于橋梁結(jié)構(gòu)線彈性穩(wěn)定安全系數(shù)的談?wù)揫J]．廣西大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，34(6)，719-724.