朱克兆,易壯鵬

(1.長江勘測規(guī)劃設(shè)計研究有限責(zé)任公司,湖北武漢 430010; 2.長沙理工大學(xué)土木與建筑學(xué)院,湖南長沙 410114)

0 引言

橋梁結(jié)構(gòu)的穩(wěn)定性是關(guān)系到安全和經(jīng)濟(jì)的主要問題之一,它與強(qiáng)度問題有著同等重要的意義,近年來隨著橋梁跨度的日益加大和橋梁結(jié)構(gòu)的日益復(fù)雜,穩(wěn)定問題更為突出。橋梁的失穩(wěn)破壞沒有明顯的征兆,帶有很大的突然性,破壞是災(zāi)難性的,在橋梁工程設(shè)計中應(yīng)該避免。

拱橋是以承受壓力為主的結(jié)構(gòu)體系,在過去的橋梁實踐中也曾多次發(fā)生拱橋施工中失穩(wěn)的事故,造成國家和人民生命財產(chǎn)的重大損失,拱橋的穩(wěn)定安全性[1～4]是設(shè)計和施工中需要考慮的重要因素之一。大跨度拱橋存在幾何缺陷時,即在制造、運輸、安裝以及運營等環(huán)節(jié)中拱肋軸線與理想線形存在偏離的情況下,拱肋在外荷載作用下將產(chǎn)生附加內(nèi)力和變形,初始撓度將對拱橋產(chǎn)生不利影響。大跨拱橋考慮幾何缺陷的穩(wěn)定問題是極值點失穩(wěn),極限承載能力分析必須考慮大位移的影響,同時由于拱肋變形很大,應(yīng)變和應(yīng)力隨之增加,線彈性本構(gòu)關(guān)系不再適用,因此須同時考慮幾何、材料非線性的影響。幾何缺陷從分布的范圍來看可以分為整體缺陷和局部分布[5]。整體缺陷在橋跨范圍內(nèi)沿拱肋均有分布,目前關(guān)于拱橋的整體缺陷的分析方法多采用一致缺陷模式法[3,4],實際工程中的幾何缺陷分布具有很大的隨機(jī)性[6,7],筆者采用隨機(jī)缺陷模式法描述幾何缺陷,并得到相關(guān)的分布方式。局部缺陷是指拱肋的某個部位存在局部的初始撓度,主要由節(jié)段的初始連接不當(dāng)和局部損傷造成,它對穩(wěn)定性的影響在目前已有文獻(xiàn)中研究較少。

大跨拱橋考慮幾何缺陷時其穩(wěn)定性能必將受到影響,不同分布方式和大小的整體將使得拱橋的承載能力有一定程度的下降,局部缺陷所處的位置和大小對拱橋穩(wěn)定性的影響目前尚不清楚?；诖?本章從拱橋的兩類穩(wěn)定問題出發(fā),研究大跨度拱橋中幾何非線性、材料非線性和初始幾何缺陷對穩(wěn)定性能的影響,探討大跨拱橋在考慮各種幾何缺陷時的屈曲特性和承載能力。

1 大跨度拱橋的穩(wěn)定性理論

拱橋的穩(wěn)定性分析有分支點失穩(wěn)和極值點失穩(wěn)兩類受力性質(zhì)根本不同的失穩(wěn)形式,對應(yīng)有線彈性屈曲分析和非線性極限承載能力分析兩種求解方法。前者假定結(jié)構(gòu)失穩(wěn)時處于線彈性范圍之內(nèi),結(jié)構(gòu)內(nèi)力與外荷載成比例關(guān)系,將穩(wěn)定分析轉(zhuǎn)化為求解特征值問題,得出的最小特征值就是臨界穩(wěn)定系數(shù);后者則采用考慮了結(jié)構(gòu)非線性的有限元法來求解大跨拱橋的穩(wěn)定問題,它通過結(jié)構(gòu)從加載開始到失穩(wěn)全過程的響應(yīng),得出結(jié)構(gòu)的荷載位移曲線,求解極限承載能力。線性屈曲分析概念清楚,計算簡便,它可以預(yù)測拱結(jié)構(gòu)屈曲的上限;非線性的極限承載能力分析則可以考慮結(jié)構(gòu)的幾何、材料非線性和幾何缺陷等因素的影響,計算準(zhǔn)確,精度高。由于這兩類分析方法各自的優(yōu)點以及兩類穩(wěn)定之間良好的相關(guān)性,本文采用這兩種方法來分析大跨度拱橋考慮幾何缺陷時的穩(wěn)定性能。

UL列式下總體坐標(biāo)系下結(jié)構(gòu)的增量有限元平衡方程[8]為:

式中:

其中[tK0]和[tKσ]分別為結(jié)構(gòu)的彈性剛度矩陣和幾何剛度矩陣;{t+ΔtQ}為外荷載等效節(jié)點力向量;{tF}為單元應(yīng)力的等效節(jié)點力向量;[tBL]和[tBNL]分別為格林應(yīng)變的線性部分和非線性部分與位移的轉(zhuǎn)換矩陣;[D]是材料矩陣;[tτ]和{tτΛ}分別為Cauchy應(yīng)力矩陣和向量，這些元素都是對應(yīng)時間t位形；tV是t時刻的體積。

1.1 線彈性屈曲分析

在式(1)中如果不考慮荷載項，且以初始荷載作用產(chǎn)生的單元內(nèi)力來計算幾何剛度矩陣，則有：

式中結(jié)構(gòu)坐標(biāo)系下的彈性剛度矩陣和幾何剛度矩陣均在初始未變形的結(jié)構(gòu)位置建立。拱結(jié)構(gòu)存在分支失穩(wěn)時,上式有非零解且結(jié)構(gòu)有分支位移產(chǎn)生,其左端剛度矩陣的行列式必為零,即有:

由上式計算得到的最小特征值λmin即為臨界荷載安全系數(shù),這樣結(jié)構(gòu)的線彈性屈曲分析就轉(zhuǎn)換成了求最小特征值問題。

1.2 非線性穩(wěn)定分析

拱橋的非線性穩(wěn)定分析可采用一定的計算方法直接對式(1)進(jìn)行求解,得到結(jié)構(gòu)從加載直至失穩(wěn)的響應(yīng)全過程,荷載-位移曲線上的最大荷載值就是結(jié)構(gòu)的極限荷載。在求解過程中,若只考慮幾何非線性因素,則得出彈性大位移穩(wěn)定極限荷載,而同時考慮幾何與材料非線性影響,就可求出結(jié)構(gòu)的彈塑性大位移穩(wěn)定極限荷載。大跨度拱橋考慮幾何缺陷的極限荷載求解必須考慮非線性因素,本文分析拱橋的非線性穩(wěn)定時分別計算彈性大位移和彈塑性大位移穩(wěn)定承載能力,并進(jìn)行對比分析。幾何非線性分析通過UL列式基礎(chǔ)上采用增量法和弧長法相結(jié)合進(jìn)行迭代來實現(xiàn),材料非線性則采用高斯積分法實現(xiàn)。

高斯積分法在單元的截面上和軸向取一定數(shù)量的高斯點,對每個高斯點根據(jù)其應(yīng)力和應(yīng)變狀態(tài)進(jìn)行數(shù)值積分,形成單元剛度矩陣,可方便計算結(jié)構(gòu)的彈塑性剛度矩陣和等效節(jié)點力,用較少積分點達(dá)到較高精度。每個積分點的被積函數(shù)[Bi]T[D][Bi]按彈性或彈塑性計算,然后按全單元數(shù)值積分即可得到單元的彈塑性剛度矩陣[Kep],即:

其中ng是高斯積分的點數(shù);分別為形函數(shù)和自然坐標(biāo)與笛卡爾坐標(biāo)之間的亞克比轉(zhuǎn)換矩陣在高斯積分點的取值。

2 結(jié)構(gòu)計算模型及幾何缺陷的分布

2.1 結(jié)構(gòu)計算模型

如圖1-a所示的某下承式鋼管混凝土系桿拱橋,主跨計算跨徑為125 m,計算矢高為25m,橋面凈寬25m;拱軸線為二次拋物線,拱肋采用啞鈴形截面鋼管混凝土,截面總高為2.6m,每個鋼管外徑為110 cm,鋼管及腹板壁厚為2 cm,內(nèi)充50號微膨脹混凝土;橋兩側(cè)各布置有8束27×Φj15(7φ5)高強(qiáng)鋼絲作為系桿;吊桿間距為7m,上端錨固于拱肋頂部,下端錨固于橫梁底;橫梁、縱梁以及蓋梁采用預(yù)應(yīng)力混凝土結(jié)構(gòu),混凝土采用30號,橋面系和拱肋相交處設(shè)置端橫梁;橫撐由拱頂處1道X字形橫撐和順橋向每邊2道橫撐組成;橋面板采用高度為35 cm的預(yù)制空心板。

設(shè)計荷載標(biāo)準(zhǔn)為城-A級,人群3.5 kN/m2。在對下承式鋼管混凝土系桿拱橋進(jìn)行穩(wěn)定性分析時,結(jié)構(gòu)的自重荷載保持不變,橫梁、縱梁以及橫撐按實際截面輸入;橋面系不按實際情況考慮,而是將橋面板、橋面鋪裝和橋面的活載視為均布線荷載施加在橫梁上面。穩(wěn)定性分析考慮的工況為恒載+橋面系的活載布置,非線性穩(wěn)定分析通過不斷加大橋面的活載布置使結(jié)構(gòu)拱橋的主要承重結(jié)構(gòu)達(dá)到極限荷載。

圖1 下承式鋼管混凝土系桿拱橋及幾何缺陷

2.2 材料本構(gòu)關(guān)系的選取

鋼管混凝土的本構(gòu)關(guān)系采用統(tǒng)一理論[9]考慮,將鋼管混凝土視為鋼管和混凝土組合而成的單一材料,其多線性隨動強(qiáng)化的軸向應(yīng)力—應(yīng)變關(guān)系如圖2所示。ξ和ξ0分別為鋼管混凝土約束效應(yīng)系數(shù)及其設(shè)計值,ξ=Asfv/Acfc,其中As,Ac分別為截面中鋼材與混凝土面積;fv和fc分別為鋼材屈服強(qiáng)度和混凝土抗壓強(qiáng)度設(shè)計值。當(dāng)ξ≥ξ0時曲線分為彈性(OA)、彈塑性(AB)與強(qiáng)化(BC)3個階段;而在ξ＜ξ0時,曲線可分為彈性(OA)、彈塑性(AB)、強(qiáng)化(BC′)與下降(C′D)4個階段。圖2中彈性模量、強(qiáng)度及應(yīng)變指標(biāo)由以下各式確定。

鋼管混凝土軸壓組合強(qiáng)度標(biāo)準(zhǔn)值fscv:

式中:A=1.212;B=0.974+0.175 9?fv/235;C=0.030 9-0.103 8?fck/20;fv為鋼材屈服點;fck為混凝土軸心抗壓強(qiáng)度標(biāo)準(zhǔn)值;ξ為套箍系數(shù)。軸壓組合比例極限fscp和比例應(yīng)變εscp分別為:

本文鋼管混凝土的屈服強(qiáng)度fscv為63.04MPa,fscp為48.53MPa,抗壓彈性模量Esc為:

在彈塑性階段,假定切線模量Esct按二次拋物線變化,且有:

鋼材應(yīng)力-應(yīng)變曲線采用理想彈塑性模型,鋼材的屈服應(yīng)力和屈服應(yīng)變分別為σv和εv。

圖2 鋼管混凝土的本構(gòu)關(guān)系

2.3 整體缺陷的分布

拱橋的整體幾何缺陷在拱跨范圍內(nèi)均有分布,它既可以采用一致缺陷進(jìn)行描述,也可以采用隨機(jī)缺陷進(jìn)行描述。前者通過某種事先假定的幾何外形來模擬結(jié)構(gòu)的初始缺陷,后者則基于隨機(jī)理論得到結(jié)構(gòu)的隨機(jī)幾何缺陷分布方式。

已有文獻(xiàn)[3,4]大多采用屈曲模態(tài)模擬缺陷的一致描述方法來考慮大跨拱橋中幾何缺陷的分布,這種方法將大跨拱橋特征屈曲分析得到的屈曲模態(tài)作為幾何缺陷模式,施加到結(jié)構(gòu)上用以模擬缺陷。這是由于特征矢量屈曲形狀是最接近于實際屈曲模態(tài)的預(yù)測值,而且特征屈曲荷載作為非線性屈曲荷載的上限,可作為施加初始缺陷或擾動載荷的根據(jù)。同時由于屈曲模式是結(jié)構(gòu)潛在的位移趨勢,如果拱橋結(jié)構(gòu)的初始缺陷分布恰好與其屈曲模態(tài)吻合,將對結(jié)構(gòu)的穩(wěn)定性產(chǎn)生最不利的影響。

大跨度拱橋的拱肋在各個環(huán)節(jié)中由于各種原因產(chǎn)生的初始撓度具有很大的隨機(jī)性,基于隨機(jī)理論考慮這些幾何缺陷影響的方法中,可以將實測得到的初撓度分布進(jìn)行調(diào)和分析,然后計算極值點屈曲荷載;或者將隨機(jī)理論的方法用于分析任意分布的初撓度對屈曲的影響,建立屈曲荷載與隨機(jī)初撓度譜密度間的關(guān)系。但前者只能針對已建成的拱橋且實測工作量很大,后者公式推導(dǎo)繁雜,最后得到的也不是最不利的失穩(wěn)形式。本文通過一種新的隨機(jī)方法分析拱的初始幾何偏差,此方法將幾何缺陷導(dǎo)致的節(jié)點坐標(biāo)偏差視為隨機(jī)變量,通過節(jié)點之間一定的相關(guān)性確定結(jié)構(gòu)的條件相關(guān)矩陣,并分解得到按可能性大小排列的隨機(jī)缺陷分布方式。拱結(jié)構(gòu)中隨機(jī)分布的節(jié)點坐標(biāo)偏差用平均值、方差以及一定程度的相關(guān)性描述,為分析方便坐標(biāo)偏差假定為正態(tài)分布且有相同的數(shù)學(xué)期望和方差,節(jié)點之間的聯(lián)系以含相關(guān)長度參數(shù)γ的指數(shù)相關(guān)函數(shù)考慮。令u=[z1,z2,Λ,zn]T為n維隨機(jī)變量來模擬結(jié)構(gòu)的節(jié)點坐標(biāo)偏離向量,設(shè)0和1為拱的邊界條件,z和z′表示拱結(jié)構(gòu)(0,1)中的兩個隨機(jī)點,隨機(jī)變量均為正態(tài)分布且方差均為σ2,z和z′之間的協(xié)方差可表示為:

其中ρ(z,z′)為z和z′的相關(guān)系數(shù),假定z和z′為指數(shù)相關(guān),且其相關(guān)系數(shù)函數(shù)為:

式中γ為相關(guān)長度參數(shù),支撐邊界確定的拱結(jié)構(gòu)中,z和z′的條件協(xié)方差可寫為:

由上述條件得到的條件相關(guān)矩陣Czz可以表示為:

其中,特征向量Ψ即為幾何缺陷的分布方式,在實際分析時只選取前幾階影響較大的缺陷。圖3給出了圖1-a所示系桿拱橋通過條件相關(guān)矩陣分解得到的前三階隨機(jī)幾何缺陷中拱肋軸線在平面內(nèi)的偏差示意圖。

圖3 隨機(jī)整體缺陷下拱肋軸線偏差示意圖

本文分析大跨度拱橋在考慮幾何缺陷的穩(wěn)定性時,采用一致和隨機(jī)兩種方法模擬幾何缺陷的整體分布,分布方式包括特征矢量屈曲形狀和按隨機(jī)理論得到的前三階隨機(jī)分布方式。將這兩種描述所得結(jié)論進(jìn)行對比分析,并與無缺陷拱的結(jié)果進(jìn)行比較,確定缺陷的最不利方式以及對結(jié)構(gòu)穩(wěn)定性能造成非常不利影響的幅值大小。

2.4 局部缺陷的分布

大跨度拱橋的拱肋存在局部幾何缺陷時,必將對吊桿間拱肋的局部穩(wěn)定產(chǎn)生影響。目前研究局部缺陷較多的是柱[10]和板[11]結(jié)構(gòu),現(xiàn)行的規(guī)范均對這些構(gòu)件的整體及局部初始撓度作了限制以避免出現(xiàn)缺陷降低構(gòu)件的承載能力。而拱肋存在局部缺陷時結(jié)構(gòu)穩(wěn)定性的研究現(xiàn)在鮮有文章涉及,拱肋初始局部撓度的大小和位置與承載能力之間的具體關(guān)系尚不清楚。針對這一問題本文假定拱肋存在一定的局部缺陷,研究它對拱橋結(jié)構(gòu)穩(wěn)定性能的影響,這種初始撓度假定沿拱肋呈正弦形式分布,如圖1所示,考察它在拱肋不同位置和大小時拱橋的屈曲特性和承載能力。

3 計算分析

3.1 整體缺陷的影響

拱結(jié)構(gòu)整體幾何缺陷的分布包括特征屈曲形狀缺陷和前三階隨機(jī)分布缺陷,這些初始幾何缺陷的有限元考慮方法是對完善拱結(jié)構(gòu)進(jìn)行坐標(biāo)重新修正。將每種缺陷的幅值取相同值用以比較分布方式對拱橋結(jié)構(gòu)穩(wěn)定性的影響;另外對于某種確定的缺陷不斷變化其缺陷幅值觀察它對結(jié)構(gòu)穩(wěn)定性的影響,并判定使得結(jié)構(gòu)失穩(wěn)的缺陷幅值。

表1給出了缺陷幅值取為L/1 000情況下,即取缺陷幅值為12.5 cm時,拱橋結(jié)構(gòu)含有這4種缺陷時的穩(wěn)定系數(shù)計算結(jié)果。從中可以看出在各種不同分布方式幾何缺陷下,彈性大變形即僅考慮幾何非線性的穩(wěn)定系數(shù)與線彈性情形相差不大,而同時考慮幾何、材料非線性的彈塑性大變形穩(wěn)定系數(shù)則明顯下降,這說明了在拱橋的穩(wěn)定性計算過程中材料非線性是主要的因素,它使結(jié)構(gòu)的承載能力明顯下降。另一方面,雖然各種缺陷的幅值均取為L/1 000,但是不同的分布方式對穩(wěn)定系數(shù)影響不同,其中特征屈曲形狀缺陷和二階隨機(jī)缺陷對拱結(jié)構(gòu)的穩(wěn)定性影響最為不利,在彈塑性大變形分析中,二階隨機(jī)缺陷可使得穩(wěn)定系數(shù)相對無缺陷情形下降15.5%。為了比較這4種缺陷在其幅值發(fā)生變化時拱橋的穩(wěn)定性能,圖4給出了缺陷幅值在變大的情況下拱橋彈性大位移穩(wěn)定系數(shù)的分布規(guī)律,從中可以看出缺陷幅值達(dá)到一定的程度時,安全系數(shù)下降明顯,其中尤以特征屈曲形狀缺陷和二階隨機(jī)缺陷影響最大,當(dāng)幅值為L/500時這兩種缺陷使得安全系數(shù)下降了近50%,而一階隨機(jī)缺陷和三階隨機(jī)缺陷對穩(wěn)定系數(shù)的影響則比這兩種缺陷小得多。特征屈曲形狀和二階隨機(jī)缺陷外形比較相似,與結(jié)構(gòu)的潛在位移趨勢接近,因而對結(jié)構(gòu)穩(wěn)定性能影響非常不利,計算結(jié)果也很好地說明了這一點。特征屈曲形狀缺陷是按照屈曲模態(tài)與結(jié)構(gòu)位移趨勢相似得到的,而二階隨機(jī)缺陷則是基于隨機(jī)理論得到的,與它們類似的整體幾何缺陷外形是拱結(jié)構(gòu)的最不利缺陷分布方式。

表1 不同缺陷大跨拱橋穩(wěn)定系數(shù)

從圖4還可知各種缺陷在其幅值處于一定的范圍之內(nèi)時安全系數(shù)與無缺陷拱相比變化均較小,而當(dāng)缺陷幅值進(jìn)一步增大時缺陷的影響明顯加大。計算結(jié)果表明各種缺陷在其幅值小于L/1 000時可以不考慮它對穩(wěn)定性能的影響。為了說明缺陷增大時拱橋考慮缺陷影響的穩(wěn)定性,圖5給出了二階隨機(jī)缺陷在不同幅值下拱頂?shù)膹椥源笪灰坪奢d-位移曲線,從中可看出缺陷增大時承載能力有一定的下降,缺陷幅值越大下降越明顯,因而將缺陷大小控制一定的范圍內(nèi)對于保證穩(wěn)定性極為必要,計算顯示缺陷幅值小于L/1 000時含缺陷拱橋的承載能力與無缺陷拱差別不大。

3.2 局部缺陷的影響

對于在拱肋上正弦形式分布的局部幾何缺陷,為研究它對局部失穩(wěn)特性的影響,圖6給出了局部缺陷分布在不同的位置時拱橋的彈性大位移穩(wěn)定系數(shù),其中缺陷幅值為吊桿間距的1/1 000,吊桿間的局部缺陷考慮分布于拱肋上和拱肋下兩種情形,其計算結(jié)果與無缺陷時穩(wěn)定系數(shù)進(jìn)行了對比。計算結(jié)果顯示局部缺陷的位置對穩(wěn)定性有一定程度的影響:較為敏感的位置位于拱腳和跨中,當(dāng)局部缺陷分布在這些位置時考慮缺陷的穩(wěn)定系數(shù)與無缺陷情形的差別較大;而當(dāng)局部缺陷位于拱跨徑的1/4附近時缺陷的影響很小。另外由圖6還可看出局部缺陷在拱肋同一位置而分別位于拱肋上方和下方時,局部缺陷的影響程度基本相同,但是它們使得穩(wěn)定系數(shù)變化趨勢恰好相反。不過從整體上來說,不論局部缺陷的位置如何,穩(wěn)定系數(shù)相對無缺陷情形變化程度還是相當(dāng)小的,缺陷的影響均在1%內(nèi)。

圖4 不同缺陷幅值下的非線性穩(wěn)定系數(shù)

圖5 二階隨機(jī)缺陷下的拱肋荷載-位移曲線

為分析局部缺陷大小對穩(wěn)定性能的影響,圖7給出了局部缺陷取不同幅值時拱肋的彈性大位移荷載-位移曲線,其中局部缺陷的分布位置取為圖1-c中編號為1位置的拱肋下方和編號為7位置的拱肋上方。圖7的計算結(jié)果表明不同缺陷幅值下拱頂?shù)暮奢d-位移曲線幾乎完全重合,局部缺陷的大小基本上不對拱橋的承載能力產(chǎn)生影響,因此可得出在計算大跨拱橋的穩(wěn)定性和承載能力時可以不用考慮拱肋局部缺陷影響的結(jié)論。

圖6 不同局部缺陷位置的非線性穩(wěn)定系數(shù)

圖7 局部缺陷幅值變化時拱肋荷載-位移曲線

3.3 拱橋的失穩(wěn)特性

在分析大跨拱橋的穩(wěn)定性時,結(jié)構(gòu)的非線性是必需考慮的因素。為此圖8給出了拱頂?shù)暮奢d-位移曲線分別按照線彈性、彈性大位移和彈塑性大位移3種情況的計算結(jié)果,其中整體幾何缺陷按照隨機(jī)分布方式2施加,缺陷幅值為L/1 000,局部缺陷不予考慮。計算結(jié)果顯示在3種分析方法中,線彈性分析方法所得極限荷載最大,僅考慮幾何非線性的彈性大位移分析方法次之,同時考慮幾何、材料非線性影響的彈塑性大位移分析方法最小,所以在實際大跨拱橋中線性方法和僅考慮結(jié)構(gòu)幾何非線性的方法均過高估計結(jié)構(gòu)的承載能力,穩(wěn)定性計算偏于不安全,在計算拱橋的極限承載能力時必須考慮材料非線性。

圖8 特征屈曲缺陷下拱肋荷載-位移曲線

圖9則給出了在彈塑性大位移分析過程中,塑性區(qū)域的出現(xiàn)順序及分布位置。荷載從0增加到43.5 kPa的過程,即圖9所示的加載階段1,整個結(jié)構(gòu)處于彈性階段,拱肋沒有塑性區(qū)域的分布;荷載從43.5 kPa變大至62.1 kPa時,這一階段結(jié)構(gòu)的靜力非線性分析應(yīng)考慮幾何非線性;當(dāng)荷載超過62.1 kPa時,左端拱腳上緣拱肋進(jìn)入彈塑性狀態(tài),對應(yīng)圖9所示加載階段2;而當(dāng)荷載增至83.2 kPa時,左端拱腳處拱肋的上緣已進(jìn)入塑性屈服狀態(tài),下緣則進(jìn)入了彈塑性狀態(tài),同時1/4處拱肋上緣也開始出現(xiàn)塑性區(qū)域分布,如圖9所示加載階段3;當(dāng)施加荷載達(dá)到95.7 kPa后,拱腳處上、下緣的塑性區(qū)已經(jīng)貫通,形成一個塑性鉸,拱肋1/4處的塑性區(qū)擴(kuò)大;當(dāng)荷載進(jìn)一步增大時,拱肋的塑性區(qū)域繼續(xù)增加,結(jié)構(gòu)位移迅速增大,直至結(jié)構(gòu)發(fā)生失穩(wěn)破壞。

圖9 極限加載過程的拱肋塑性區(qū)域位置

4 結(jié)論

本文在進(jìn)行大跨度拱橋穩(wěn)定性分析時,考慮了幾何缺陷、幾何非線性和材料非線性等因素的影響,幾何缺陷的分布從整體和局部上分別考慮,整體缺陷包括由隨機(jī)和一致兩種分布方式,局部缺陷通過吊桿間拱肋的局部初撓度考慮,得出了以下結(jié)論:

1)沿拱肋分布的各種整體幾何缺陷,當(dāng)其幅值小于跨徑的1/1 000時,可以不用考慮整體缺陷的影響。特征矢量屈曲形狀和二階隨機(jī)分布方式的缺陷對結(jié)構(gòu)的穩(wěn)定性能影響最大,在其幅值達(dá)到一定程度時結(jié)構(gòu)承載能力有明顯的下降,它們是穩(wěn)定性分析必須考慮的影響因素。

2)對于局部缺陷,通過計算發(fā)現(xiàn)存在缺陷時拱橋的穩(wěn)定系數(shù)雖然和無缺陷拱橋存在差異,但是缺陷的影響很小,局部缺陷位于拱肋不同位置時非線性穩(wěn)定系數(shù)相差不到1%,而且局部缺陷的幅值增大時,缺陷對拱肋的荷載-位移曲線的影響基本上可以忽略。

3)通過計算發(fā)現(xiàn)線彈性和僅考慮幾何非線性的分析方法會過高估計結(jié)構(gòu)的承載能力,大跨拱橋的穩(wěn)定性分析必須考慮材料的非線性因素?？紤]材料非線性的彈塑性大位移分析時,拱肋的塑性區(qū)首先出現(xiàn)在拱腳位置上緣,達(dá)到極限荷載時分布區(qū)域包括整個拱腳與1/4拱肋的上緣。

[1]李國豪.橋梁結(jié)構(gòu)的穩(wěn)定與振動[M].北京:中國鐵道出版社,1996.

[2]程 進(jìn),江見鯨,肖汝誠,等.大跨度鋼拱橋結(jié)構(gòu)極限承載力分析[J].工程力學(xué),2003,20(2):7-10.

[3]謝 旭,李 輝,黃劍源.大跨度兩鉸鋼拱橋面內(nèi)穩(wěn)定分析[J].土木工程學(xué)報,2004,37(8):43-49.

[4]徐文平,裝銘海,戴 捷,等.大跨簡支系桿拱橋非線性穩(wěn)定性[J].交通運輸工程學(xué)報,2005,5(4):53-57.

[5]Palamas.J,Coussy.O,Bamberger.Y.Effects of surface irregularitiesupon the dynamic response of bridges under suspended moving loads[J].Journal of Sound and Vibration,1985,99(2):235-245.

[6]Most.T,Bucher.C,Schorling.Y.Dynamic stability analysis of non-linear structures with geometrical imperfections under random loading[J].Journal of Sound and Vibration,2004,276(1-2):381-400.

[7]Schorling.Y,Bucher.C.Dynamic stability analysis for structures with geometrical imperfections[J].In:Structural Safetyand Reliability,Rotterdam,1998:771-777.

[8]王勖成,邵 敏.有限單元法基本原理和數(shù)值分析[M].北京:清華大學(xué)出版社,1996.

[9]韓林海.鋼管混凝土結(jié)構(gòu)[M].北京:科學(xué)出版社,2000.

[10]沈祖炎,張其林.受壓方管鋼柱的屈曲后極限承載力[J].土木工程學(xué)報,1991,24(3):15-26.

[11]陳國棟,郭彥林.非加勁板抗剪極限承載力[J].工程力學(xué),2003,20(2):49-54.