基于纖維模型的梁橋材料非線性研究*

楊國俊　郝憲武　段瑞芳　李子青

(長安大學公路學院1)　西安　710064)　(陜西交通職業技術學院2)　西安　710018)

?

基于纖維模型的梁橋材料非線性研究*

楊國俊1)郝憲武1)段瑞芳2)李子青1)

(長安大學公路學院1)西安710064)(陜西交通職業技術學院2)西安710018)

摘要：為了提高工程結構非線性的計算精度，研究了OpenSEES程序中的纖維模型，以多梁式簡支T梁橋為背景，建立體系的纖維模型，通過偏載、中載、均勻加載三種方式，分析各片主梁的材料非線性變化過程.在偏載加載方式下，數值模擬的各主梁荷載-位移，以及荷載-鋼筋應力曲線分別與破壞性試驗數據進行了對比分析.結果表明，模擬的誤差很小，兩者結果相吻合，說明OpenSEES有限元的纖維模型分析非線性具有更高的精度和優越性.

關鍵詞：橋梁工程；纖維模型；材料非線性；多梁式鋼筋混凝土；本構關系

楊國俊(1988- )：男，博士生，主要研究領域為橋梁橋梁結構及結構分析理論

*中央高校基本科研業務費資助 (批準號：2013G1211010)

0引言

大量的試驗研究表明，鋼筋混凝土結構存在很明顯的非線性特性[1-3]，這是由于組成鋼筋混凝土的材料本身存在非線性.文獻[4-5]指出結構的彈性分析不能較好地反映混凝土和鋼筋之間的變形協調問題，不能合理的反映在長期荷載作用下混凝土構件的內力和變形，所以按照線彈性理論方法求得的結構內力和變形不能反映實際情況.對橋梁結構的材料非線性研究較多，文獻[6-7]研究材料非線性采用了廣義逆力(generalized inverse matrix，GIM)分析方法，該方法在解決材料非線性問題時無需像傳統的基于位移法的逐步增量法那樣逐步求解，但是理論性較強，且在非線性方程組求解過程中穩定性不太好；文獻[8-10]基于有限元軟件分析法，研究了分層有限單元法、分層殼元法以及材料非線性的平面桿系結構有限元分析方法，這些方法各有特點，同時存在一定的不足之處，由于單元劃分復雜，非線性方程組求解過程中不一定收斂，存在穩定性和精確性不夠的缺點；文獻[11]研究材料非線性時采用實測值的反分析方法，該方法需要極高的模型辨識方法，引入進化的思想，利用遺傳算法和遺傳規劃的全局優化搜索技術，進行材料非線性模型識別的進化方法研究，這種方法理論較復雜，在實際工程中很難應用.

綜上所述，很少有文獻研究多梁式T梁橋的材料非線性問題，文中采用軟件OpenSEES[12-13]的纖維模型分析多梁式鋼筋混凝土結構材料非線性，基于OpenSEES這個軟件本身存在的優點能夠克服非線性方程組求解過程中穩定性和精確性不夠的缺點，建立纖維模型可以將構件截面劃分，從而同一截面的不同纖維可以有不同的單軸應力應變關系，這樣就可以采用更加符合構件受力狀態的單軸本構關系，更好的模擬鋼筋和混凝土之間的實際工作性能，用OpenSEES計算結果與試驗結果數據進行了對比分析，結果相吻合，為多梁式鋼筋混凝土橋梁的材料非線性分析提供了一種實用的計算方法.

1纖維模型

纖維模型是OpenSEES程序內部的一種模型，該模型將構件截面劃分成很多纖維，對于鋼筋混凝土結構，截面就包括鋼筋纖維和混凝土纖維，對于每一根纖維考慮它自身的軸向本構關系，不同的纖維就有不同的本構關系.構件的纖維單元劃分見圖1.

圖1　纖維模型劃分示意圖

纖維模型截面分析前，首先應該對構件的橫截面進行劃分.混凝土按其位置及所受約束情況進行劃分.OpenSEES根據構件截面上的保護層混凝土和受約束區核心混凝土的應力應變關系不同將混凝土纖維區分為OpenSEES中的coverconcrete和coreconcrete，鋼筋纖維應該注意使鋼筋單元與原截面鋼筋的形心相重合，劃分見圖2.

圖2　鋼筋混凝土T梁纖維截面圖

2本構關系

2.1混凝土本構關系

纖維模型中混凝土受壓時的單調加載骨架曲線(見圖3)采用Scott等改進后的Kent-Park模型[14].該模型綜合考慮到橫向箍筋對核心混凝土的約束影響，所以通過用修正參數修正圖中混凝土骨架曲線的峰值應力應變和軟化段的斜率來實現目標，其中受壓骨架曲線表達式如下.

(1)

(2)

式中：K為混凝土強度增加系數，0.002K是對應的峰值應變，

(3)

Zm為混凝土應變軟化段斜率，

針對混凝土受箍筋加拉筋約束的情形，Scott給出了混凝土極限壓應變偏保守的取值方法，按下式計算.

(5)

圖3　混凝土單調加載骨架曲線

2.2鋼筋本構關系

該模型為考慮等向應變硬化影響的本構關系模型，見圖4.

圖4　Menegotto-Pinto鋼筋模型

該本構關系模型計算上非常有效率是由于采用了應變的顯函數表達形式，同時又保持了與鋼筋重復加載試驗所得的結果具有很好的一致性，可以反映包興格(bauschinger)效應，模型數學表達式如下.

(6)

式中：

(7)

其中：(ε0，σ0)和(εr，σr)的具體含義見圖4；b為應變硬化率；R為影響過渡曲線形狀的參數，它反映了包興格效應，即

(8)

式中：將隨每次應變的反向而更新數值；R0為首次加載時R的初始值；a1和a2是同R0由試驗確定的參數值.為了解決鋼筋等向硬化的問題，Filippou等指出將線性的屈服漸近線進行應力平移，塑性應變最大值的大小決定的平移大小，表達式為

(9)

式中：εmax為應變反向時應變的最大絕對值；σy，εy分別為屈服應力和屈服應變；a3與a4為由試驗確定的參數值.

3工程實例

3.1工程概況

以某多梁式T梁橋為背景橋，鋼筋混凝土多梁式的上部結構單跨橫向布置4片T梁，T梁截面及鋼筋布置見圖5、圖6，該橋跨徑為9.2 m.分別在跨中和支座處設置橫隔板.該橋使用了40 a后，在1993年廢除使用，廢除前進行了現場破壞性試驗，現場試驗對T梁進行混凝土鉆心取樣，實際測得混凝土的彈性模量Ec=23 530 MPa，混凝土的抗拉強度ft=2.21 MPa，抗壓強度fc=23.5 MPa，對鋼筋試樣進行測試得到鋼筋的屈服強度fy=254.8 MPa.

圖5　橋梁橫斷面布置圖(單位：mm)

圖6　單片T梁截面布置圖(單位：mm)

3.2加載方式

采取3種加載方式，即中載、偏載,以及均勻加載，對該橋做了偏載的破壞試驗，并做了相關的多梁式梁橋全過程破壞分析，編制組合分層殼單元的非線性有限元程序，對該橋進行了非線性有限元的模擬分析[16-17].

1) 偏載加載方式為了模擬實際的車輛作用，采用4個加載點施加集中荷載，每級施加3t的力，逐級加載直到梁破壞以測試多梁式T梁的極限承載能力，在破壞性試驗中，實際加載的配備裝備及加載方式見圖7、圖8，其中G1，G2，G3，G4分別表示第一片至第四片梁.

圖7　偏載加載方式(單位：mm)

圖8　偏載加載設備

在偏載加載方式下，每片主梁的試驗數值與有限元數值的荷載-位移對比曲線見圖9，從圖中分析可知計算值和試驗值比較吻合，說明OpenSEES纖維模型模擬具有較高準確度的特點，由于偏載方式加載，把荷載加到主梁G1和G2的跨中處，使得G1和G2主梁跨中位置成為撓度敏感區，所以其試驗數據與OpenSEES有限元數據計算結果相當接近，吻合程度高，而相對于主梁G3和G4而言，不是撓度敏感區，尤其是主梁G4，試驗數據本身也不光滑而且存在誤差，所以試驗數據和OpenSEES有限元計算結果不是太吻合，存在一定的誤差.

圖9　各片主梁跨中的荷載-撓度曲線(試驗數據對比)

圖10為各主梁跨中的荷載-撓度.由圖10可見，由于偏載方式加載，每片梁進入屈服階段時的屈服點不同，在相同荷載作用下，主梁G1的跨中撓度最大，剛度折減現象最嚴重，而主梁G4還幾乎處于彈性階段，沒有達到屈服狀態，剛度折減現象不太明顯.由于是偏載加載方式，當數值模擬的實際加載量在105 t時，G1梁完全破壞，整個結構失效，但是G4梁還仍然存在一定的剛度，還能承受荷載.然而，實際試驗測得實際加載量為102 t時，才達到上述狀態，說明數值模擬和試驗數據誤差很小，誤差僅是2.9%.

圖10　各主梁跨中的荷載-撓度

圖11　各梁受拉區縱向鋼筋應力隨荷載的變化

圖11為各梁受拉區縱向鋼筋應力隨荷載的變化圖.由圖11可見，由于偏載加載方式，G1主梁的受拉區鋼筋應力最大而且應力發展相對最快，荷載加載到24級左右時，即加載量達到72 t時，G1主梁的受拉區縱向鋼筋首先開始屈服，進入屈服階段，之后鋼筋進入硬化階段，應力發展較緩慢.主梁G2的受拉區縱向鋼筋加載量達到78 t時，進入屈服階段，主梁G3的受拉區縱向鋼筋加載量達到102 t時，進入屈服階段，此時，整個結構完全破壞.然而，只有G4主梁的受拉區縱向鋼筋還尚未屈服，仍然處于彈性階段.

2) 中載加載圖12為各主梁跨中的位移-荷載等級.由圖12可見，由于中載方式加載，荷載施加在主梁G2和G3上，G1主梁和G4主梁荷載-位移曲線相同，G2和G3主梁荷載-位移曲線相同.計算分析說明，荷載加載至25級左右時，即加載量為75 t跨中位移在15 mm時，該結構進入非線性階段，當繼續加大荷載，當34.5級時，即加載量為103.5 t時，跨中位移18 mm左右時，結構發生破壞.從圖13分析可知，在中載加載方式下，G1主梁和G4主梁的受拉區縱向鋼筋應力變化是相同的，而G2主梁和G3主梁的受拉區縱向鋼筋應力變化是相同的，荷載加至28級時，主梁G2和G3的受拉區縱向鋼筋開始屈服，繼續增加荷載，當荷載增加到34.5級時，主梁G1和G4的受拉區縱向鋼筋開始屈服.與偏載加載方式相比較，中載加載時荷載橫向分布系數比偏載時較均勻，所以中載加載時鋼筋屈服明顯的滯后于偏載加載.在偏載加載直到結構破壞時，主梁G4的受拉區鋼筋還處于彈性階段，沒有達到屈服狀態，而中載作用下，主梁G1和G4的受拉區鋼筋在荷載加至34.5級時剛開始屈服.

圖12　各主梁跨中的位移-荷載等級

圖13　受拉區鋼筋應力隨加載等級的變化曲線

3) 均勻加載圖14為各主梁在均勻加載下跨中位移-荷載曲線.由圖14知，由于均勻加載，各主梁的荷載-撓度曲線基本是一樣的，即G1，G2，G3和 G4主梁荷載-位移曲線重合在一起.計算分析說明，當荷載加至48 t，跨中位移為16 mm時，結構進入非線性階段，繼續加載至破壞時上部結構各片主梁的跨中位移約為28 mm，此時該橋上部結構的極限承載能力約為54 t.

圖14　各主梁在均勻加載下跨中位移-荷載曲線

圖15為各梁受拉區鋼筋應力隨加載等級的變化曲線.由圖15可見，在均勻加載方式下，各主梁受拉區縱向鋼筋應力基本相同的.當荷載加至16.5級，即實際加載量為49.5 t時，所有梁的梁底受拉區縱向鋼筋基本同時達到屈服階段，從此開始進入非線性階段，由于鋼筋的硬化作用，鋼筋應力隨加載增加，但增加的相當緩慢，之后隨著荷載的繼續增加，非線性越來越明顯，當加載到54 t時，結構開始破壞.

圖15　各梁受拉區鋼筋應力隨加載等級的變化曲線

縱向比較3種加載方式，偏載加載方式下，主梁G1的受拉區縱向鋼筋加至24級時進入屈服階段，在中載下，荷載加至28時，主梁G2和G3的受拉區縱向鋼筋開始屈服，而在勻載下，荷載加至16.5級時，所有梁的受拉區縱向鋼筋同時屈服，綜上所述，在均勻加載方式下主梁鋼筋屈服明顯早于偏載方式和中載方式.

4結論

1) OpenSEES有限元軟件能夠克服非線性方程組求解過程中穩定性和精確性不夠的缺點.

2) 在偏載加載方式下，提取了多梁式簡支T梁橋破壞性試驗數據，并與OpenSEES軟件模擬的計算值進行了對比分析，誤差是2.9%，具有較高的精度.

3) 從荷載等級與位移關系曲線及鋼筋應變隨荷載等級的變化曲線可以展現該橋梁非線性發展變化的全過程，從線性階段到最后的破壞階段.

4) 在3種加載方式下，從荷載-位移曲線和鋼筋應力-荷載曲線中可以得出，橋梁結構均勻加載方式比其他兩種方式更容易破壞和主梁受拉區縱向鋼筋更早開始屈服.

參 考 文 獻

[1]過鎮海.混凝土的強度和變形(試驗基礎和本構關系)[M].北京:清華大學出版社,1997..

[2]王心勇,辛全才,宋娟.鋼筋混凝土結構的非線性有限元分析[J].人民黃河,2006,28(8):60-61.

[3]李揚.基于材料非線性的構件計算方法的研究[D].武漢：華中科技大學,2009.

[4]何政,歐進萍.鋼筋混凝土結構非線性分析[M].哈爾濱：哈爾濱工業大學出版社,2006.

[5]劉禹.纖維布加固鋼筋混凝土抗彎構件的非線性有限元分析[D].西安：西安理工大學,2008.

[6]劉西拉,王開健.材料非線性問題的廣義逆力法有限元格式表達[J].巖石力學與工程學報,2004,23(21):3629-3635.

[7]J Z C, L L X. A large increment method for material[J]. Advances in Structural Engineering,1997(1):99-110.

[8]潘家英,張國政.大跨度橋梁極限承載力的幾何與材料非線性耦合分析[J].土木工程學報,2000,33(1):5-8.

[9]張劍,葉見曙,陳艷,等.基于分層殼元法多梁式T梁的極限承載力評估[J].公路交通科技,2007,24(5):81-85.

[10]曾永革.考慮材料非線性的平面桿系有限元分析方法[J].湘潭師范學院學報：自然科學版,2009,31(2):98-102.

[11]楊成祥.材料非線性模型的進化識別方法研究[D].哈爾濱：東北大學,2001.

[12]陳學偉,韓小雷,林生逸.基于宏觀單元的結構非線性分析方法、算例及工程應用[J].工程力學,2010(A01):59-67.

[13]齊虎,孫景江,林淋.OPENSEES中纖維模型的研究[J].世界地震工程,2007,23(4):48-54.

[14]S M, F M, H S. OpenSees users manual[M]. Berkeley: Univ. of California,2006.

[15]凌炯.面向對象開放程序OpenSees在鋼筋混凝土結構非線性分析中的應用與初步開發[D].重慶：重慶大學,2004.

[16]張劍.多梁式混凝土梁橋的極限承載力研究[D].南京：東南大學，2007.

[17]SONG Ha-Won， DY K B. Finite element failure analysis of reinforced concrete T-girderbridges[J]. Engineering Structures,2002,24(2):151-162.

Research on the Material Nonlinear of Reinforced

Concrete Bridge Based on Open SEES Fiber

YANG Guojun1)HAO Xianwu1)DUAN Ruifang2)LI Ziqing1)

(SchoolofHighway,Chang’anUniversity,Xi’an710064,China)1)

(ShaanxiCollegeofCommunicationTechnology，Xi’an710013,China)2)

Abstract：In order to improve the nonlinear accuracy of engineering structure, the fiber model of OpenSEES program is studied. Taking the multi-beam reinforced concrete as an example, fiber model is established to analyze the changing process of material nonlinear with each beam through the partial, middle and uniform load. In the partial load mode, the numerical data of load-midspan displacement and load-steel stress is compared with the destructive test data. The results show the simulation error is very small, and the numerical data is consistent with the test data, which shows that analyzing material nonlinearity with the fiber element of OpenSEES finite program has higher accuracy and superiority.

Key words：bridge engineering; fiber model; material nonlinearity; multi-beam reinforced concrete; constitutive relationship

收稿日期：2015-10-28

doi:10.3963/j.issn.2095-3844.2016.01.017

中圖法分類號：U448.217