王永寶，趙人達，徐騰飛，占玉林

(西南交通大學　土木工程學院，四川　成都　610031)

?

鋼管混凝土軸壓構件徐變簡化計算方法研究

王永寶，趙人達，徐騰飛，占玉林

(西南交通大學土木工程學院，四川成都610031)

為準確快速計算鋼管混凝土軸壓構件的長期變形，基于按照齡期調整的有效模量法和逐步積分法，提出了鋼管混凝土軸壓構件徐變的簡化分析方法，對比分析了本文方法與傳統(tǒng)方法的優(yōu)缺點，探討了混凝土徐變預測模型及其相對濕度取值對鋼管混凝土軸壓構件長期變形計算精度的影響，并與既有的試驗結果進行對比分析。研究結果表明：本文簡化分析方法簡單且有較高精度，可以用于鋼管混凝土軸壓構件長期變形分析；采用相對濕度取值范圍為90%～98%的CEB90模型和EC2模型，可以較為準確地分析鋼管混凝土軸壓構件長期變形。

橋梁工程；鋼管混凝土；逐步積分法；軸壓構件；徐變模型；濕度

0　引言

高速鐵路上列車的行駛速度較高，對大跨度橋梁結構的平順性要求較高，鋼管混凝土拱橋以其獨特的優(yōu)勢在現代高速鐵路的建設中得到廣泛的應用。鋼管核心混凝土的收縮徐變會引起大跨度鋼管混凝土拱橋產生較大的豎向位移和應力重分布，這將會顯著影響高速鐵路的安全性和舒適性。因此準確分析鋼管混凝土結構的收縮徐變特征，對保證高速鐵路上大跨度鋼管混凝土拱橋的正常運營具有重大理論和實踐意義。

國內外相關學者對鋼管混凝土結構的收縮徐變進行了大量的試驗與理論研究，分析了混凝土時變作用對鋼管混凝土軸心受壓構件長期變形的影響。對鋼管混凝土軸壓構件長期變形計算的理論研究成果有YAMADA[1]和ICHINOSE[2]開爾文鏈模型，HAN L H[3-4]、YANG Y F[5]和UY[6]等的按照齡期調整的有效模量法模型，韓冰和王元豐[7-8]等的繼效流動模型，王玉銀[9-10]等的逐步積分法模型，趙金鋼[11]等的等效溫度荷載法模型，WASSIM[12]提出的鋼管對核心混凝土的緊箍力以及鋼管和核心混凝土屈服的增量計算模型等。韓冰[13]和GENG Y[14]對鋼管混凝土徐變計算方法進行了對比分析，探討了各個預測模型的計算精度。綜合比較以上分析方法可以看出，開爾文鏈模型[1-2]的參數需要試驗擬合，不適合初步設計采用；按照齡期調整的有效模量法和繼效流動模型需要進行較多的迭代計算；等效溫度荷載法需要建立有限元模型，比顯式計算方法復雜。各個學者給出的鋼管混凝土長期變形的分析方法有較大差異[13-14]，在保證精度的前提下，進一步簡化鋼管混凝土的徐變計算方法還需要深入研究。

為了提高計算速度，簡化計算過程，便于設計人員采用，本文根據鋼管混凝土徐變的特點，對按照齡期調整的有效模量法[3]和逐步積分法[10]計算公式進行了簡化，并與傳統(tǒng)分析結果進行了對比分析,探討了修正的ACI209模型[15]、CEB90模型[16]、GL2000模型[17]、EC2模型[18]對預測精度的影響；針對鋼管混凝土的密封性較好，相對濕度的取值與普通混凝土有較大差別的特點，進一步分析了CEB90模型和EC2模型選擇不同相對濕度對預測結果的影響，并給出了這兩種預測模型的相對濕度取值范圍。

1　徐變簡化計算模型

在鋼管混凝土軸壓構件受到外界荷載較小時，鋼管對核心混凝土的緊箍力較小，可以忽略不計[2，7，11，19]。在不考慮鋼管對核心混凝土的緊箍力作用時，任意時刻tk，截面上滿足內力平衡條件和變形協(xié)調條件：

(1)

(2)

當初始荷載施加后k=0時，混凝土材料滿足胡克定律：

(3)

由式(1)～(3)可以計算得到初始加載后t0時刻核心混凝土的應力：

(4)

式中，σc(tk)，εc(tk)為tk時刻核心混凝土的應力和應變；σs(tk)，εs(tk)為tk時刻鋼管的應力和應變；n=Es/Ec(t0)為鋼管和核心混凝土彈性模量比值；N為鋼管混凝土截面上承擔的軸力；Ac，As分別為核心混凝土和鋼管的截面面積；Es，Ec(t0)分別為鋼管和核心混凝土的彈性模量。

混凝土在變應力狀態(tài)下的徐變計算方法有擬合公式法[1-2]、按照齡期調整的有效模量法[3-6]，繼效流動準則方法[7-8]、老化理論方法[9-11]。擬合公式法計算簡單，但是需要試驗數據進行參數擬合，推廣性較差，不適合橋梁結構的初步設計采用；按照齡期調整的有效模量法和逐步積分法簡單，計算精度能夠滿足設計要求；繼效流動準則方法需要考慮鋼管對核心混凝土的緊箍力作用，并進行大量的迭代，計算方法復雜，對設計人員計算能力要求較高。鑒于此，本文分別選取按照齡期調整的有效模量法[20]和老化理論作為變應力狀態(tài)下核心混凝土徐變計算方法。以下分別用兩種方法推導鋼管混凝土徐變計算公式。

1.1按照齡期調整的有效模量法

變應力狀態(tài)下核心混凝土的徐變計算公式采用BAZANT教授提出按照齡期調整的有效模量法[20]，考慮混凝土的徐變后，任意時刻tk核心混凝土的應變可以表示為：

(5)

(6)

式中，φ(tk,t0)為加載齡期為t0，計算齡期為tk時混凝土的徐變系數；εsh(tk)為計算齡期為tk時混凝土的收縮應變；χ(tk,t0)為加載齡期為t0，計算齡期為tk時混凝土的老化系數，按式(6)計算[3]。

核心混凝土發(fā)生徐變以后，由式(1)和(2)可以得到tj時刻核心混凝土的應力和應變關系：

(7)

式中，σ0=N/Ac,α=As/Ac。

將式(7)代入式(5)，可以得到按照齡期調整的有效模量法計算的tk時刻鋼管混凝土的顯式徐變應變計算公式：

(8)

式(8)即為采用按照齡期調整的有效模量法給出的鋼管混凝土徐變的簡化計算公式。

1.2逐步積分法

變應力狀態(tài)下鋼管核心混凝土的徐變計算公式采用老化理論，將時間離散為t0,t1,…,tk時間點，當時間間隔Δti足夠小時，計算公式為：

(9)

式中J(tk,t0)為加載齡期為t0，計算齡期為tk時的混凝土的徐變函數。

任意時刻tk截面上還滿足內力平衡和變形協(xié)調條件，將式(7)代入式(9)中可以得到采用逐步積分法計算的tk時刻鋼管混凝土的顯式徐變應變計算公式：

(10)

(11)

式中Ekj=J(tk,tj-1)(k≥j)

目前大多數混凝土的收縮徐變預測模型，如ACI209模型、CEB90模型、EC2模型和GL2000模型，直接給出的是徐變系數表達式φ(tk,tj)：

(12)

為了進一步簡化計算公式，將式(12)代入式(10)、(11)中，可以得到由徐變系數表示的鋼管混凝土的徐變計算公式：

(13)

(14)

式(13)、(14)即為用徐變系數表示的逐步積分法鋼管混凝土徐變的簡化計算公式。采用本公式可以通過tk時刻以前的核心混凝土的應力和徐變系數，直接計算tk時刻鋼管混凝土的應變值，大大簡化了計算過程。

2　徐變簡化計算模型驗證

為驗證本文提出的簡化分析方法的正確性，以文獻[9]給出的鋼管混凝土軸心受壓構件為對象，進行了鋼管混凝土長期變形計算?；炷恋男熳冾A測模型采用CEB90模型，取密封狀態(tài)下混凝土的相對濕度為98%，構件理論厚度按照實際的計算方法確定，不考慮核心混凝土的收縮。分別用式(8)、(13)、(14)，文獻[3，10-11]給出的公式計算鋼管混凝土軸壓構件的長期變形，并將結果進行對比分析。不同公式給出的結果對比情況見表1，具體截面尺寸見表2的VП試件。由表可知，式(8)與文獻[3]的計算結果相差較小，兩組試件的徐變應變差小于2 με，完全可以滿足初步設計的要求，但是本文提出的顯式計算方法簡單，計算速度快，不需要大量的迭代。式(13)～(14)和文獻[10]的徐變應變相差小于4 με，計算結果相差不大，滿足工程要求，但式(10)、(11)用單個徐變函數代替了平均徐變函數，簡化了文獻[10]中的Ec1k和Ec2kj的計算。由于目前大多數徐變模型直接給出的是徐變系數的表達式，將式(10)、(11)簡化為可以直接代用徐變系數進行計算的式(13)、(14)，進一步簡化了計算過程。

表1　不同計算方法的鋼管混凝土試件徐變計算值對比(單位：με)

圖1給出了采用不同方法計算所得的鋼管混凝土軸心受壓構件的徐變應變與試驗結果的對比情況，由圖可知，上述5種方法的計算結果與試驗結果相差不大，這與文獻[14]中的結論是一致的。按照齡期調整的有效模量法(式(8)和文獻[3])的計算結果偏高，但是此方法可以直接計算任意時刻tk的徐變應變，不需要計算tk時刻以前的徐變應變；逐步積分法(式(13)、(14)和文獻[10])的計算結果比按照齡期調整的有效模量法偏小，比等效溫度荷載法偏大，但是需要計算tk時刻以前的徐變應變；等效溫度荷載法(文獻[11])的計算結果偏低，除需要計算tk時刻以前的應變以外，還需要建立有限元模型，計算方法較前兩種方法復雜。式(13)、(14)計算結果介于按照齡期調整的有效模量法和等效溫度荷載法之間，計算方法相對簡單，不需要將徐變模型計算的徐變系數轉化為徐變函數，計算精度與其他方法一致，鑒于以上優(yōu)點，本文選擇式(13)、(14)作為鋼管混凝土徐變計算公式。

圖1　不同計算方法的試件徐變計算值對比圖Fig.1　Comparison of calculated creep curves of CFST members obtained by different calculation methods

3　混凝土徐變模型選取

以上分析結果表明：用徐變系數表示的逐步積分法簡單，且精度較高?；谠囼灲Y果，采用CEB90模型、EC2模型、GL2000模型和修正的ACI209模型分析鋼管混凝土軸心受壓構件的徐變。CEB90模型、EC2模型、GL2000模型相對濕度取98%；ACI209模型的極限徐變系數用1.5代替文獻[15]給出的2.35[6]，相對濕度取98%，一般情況下內填混凝土添加減水劑后坍落度較大，可以取200 mm。計算試件取表2的VП試件。

表2　典型鋼管混凝土構件的幾何和材料參數

注：D(B)為圓鋼管截面外徑或方形、矩形鋼管截面長和寬；tt為鋼管壁厚度；L為鋼管混凝土試件高度；Es和Ec分別為鋼管和核心混凝土的彈性模量；fcm為圓柱體混凝土28 d的軸心抗壓強度；N為鋼管混凝土施加的軸壓荷載；t0為鋼管混凝土的加載齡期；α為鋼管混凝土的含鋼率。

圖2給出了不同混凝土徐變預測模型計算的鋼管混凝土長期變形隨時間的變化曲線，由圖可知：GL2000模型早期徐變發(fā)展較快，后期變化緩慢，與試驗結果偏差較大，這與文獻[11]的計算結果一致；ACI209模型經過修正后，計算結果與實測結果吻合程度有所提高，但是對不同的試件，都有明顯高估或低估徐變應變的趨勢，且ACI209模型對濕度的敏感性較差[24]，即使進行濕度或者體積系數的修正，其計算結果與濕度假定為98%的計算結果相差不大； EC2模型和CEB90模型計算結果與試驗結果吻合最好。以上分析表明EC2模型和CEB90模型可以作為鋼管混凝土徐變計算的預測模型。

圖2　不同徐變模型的試件徐變計算值對比Fig.2　Comparison of calculated creeps of CFST members obtained by different creep models

4　徐變簡化方法相對濕度取值分析

圖3　鋼管混凝土構件長期荷載作用下徐變應變計算值與試驗值對比Fig.3　Comparison of calculated and tested creep strains of CFST members under long-term loading

EC2模型和CEB90模型的相對濕度顯著影響預測結果，且目前對鋼管混凝土徐變模型相對濕度取值范圍的研究較為欠缺，為了分析EC2模型和CEB90模型計算中相對濕度的取值問題，本文以既有的試件為分析對象，選取相對濕度h為70%，90%，95%和98%作為計算參數，分析相對濕度取值對鋼管混凝土徐變計算預測精度的影響。表2給出了不同加載齡期，承受不同荷載大小的方鋼管、矩形鋼管混凝土和圓鋼管混凝土試件。本次選取的構件包括普通混凝土和膨脹混凝土，混凝土強度變化范圍13～52 MPa，含鋼率為6.6%～20%，加載齡期為5～34 d。

圖3給出了4種不同濕度的EC2模型和CEB90的分析結果與試驗結果對比情況。由圖可知，CEB90模型和EC2模型對相對濕度敏感性較大，當相對濕度取大氣平均濕度70%時，計算結果明顯高估試驗值，證明鋼管核心混凝土受到鋼管的密閉作用，其濕度較高，明顯大于外界大氣相對濕度；當相對濕度取90%時，短齡期加載試件與試驗結果吻合較好,見圖3(f)和(g)；當相對濕度取95%時，28 d 加載試件與試驗結果吻合較好,見圖3(c),(h)～(k)；當相對濕度取98%時，28 d加載試件與試驗結果吻合較好，見圖3(d)。而28 d加載試件在相對濕度取95%時，短期預測結果與試驗結果吻合較好，見圖3(b)、(c)。而長期試驗結果的徐變值突然增加，究其原因可能是調整荷載時荷載增加導致的，其試驗值明顯高于計算值。14 d加載試件CFST受外界溫度影響較大，其長期變形呈現較大的增長速度。文獻[10-11] 認為CEB90模型和EC2模型計算時h應為無窮，但是經計算顯示h取為無窮計算得到徐變系數與相對濕度取98%計算的結果相差不大。從圖3可以看出，如果只取構件理論厚度為無窮(即相對濕度取98%)，大多數試件的徐變計算值小于試驗值，該方法明顯低估了短齡期加載試件混凝土的徐變。因此，本文建議構件理論厚度按照原有的計算方法計算，但是相對濕度取值范圍取90%～ 98%。

當濕度取90%時，短齡期加載的試驗結果與分析結果吻合較好；當濕度取98%時，長齡期加載的試驗結果與分析結果吻合較好。而實際的鋼管混凝土的濕度應該隨著時間的推移逐漸減小，以上結論與實際不符，究其原因，可能是短齡期加載的鋼管混凝土徐變試件較少，計算結果不具有代表性，也可能還有其他因素能夠顯著影響短齡期和長齡期加載的鋼管混凝土的徐變，具體原因還需要進一步探討與分析。

5　結論

鋼管混凝土的徐變受初始應力、徐變預測模型、變應力狀態(tài)下核心混凝土計算公式以及基本假定的影響較大，本文分析了鋼管混凝土徐變的特點，基于按照齡期調整的有效模量法和逐步積分法簡化了鋼管混凝土軸壓構件的徐變計算公式，可以直接利用顯式計算公式計算鋼管混凝土的徐變，極大地簡化了計算過程，運用本文計算公式對鋼管混凝土長期變形進行了分析，將計算結果與試驗結果進行了對比分析，得到如下結論：

(1)基于按照齡期調整的有效模量法和逐步積分法推導的鋼管混凝土軸壓構件徐變簡化計算公式，其計算結果與試驗結果吻合較好，可供工程應用參考。

(2)鋼管混凝土的徐變受混凝土徐變模型相對濕度取值影響較大，大量事實表明：相對濕度取值范圍為90%～98%的CEB90模型和EC2模型在分析低應力狀態(tài)下的鋼管混凝土軸壓構件的徐變時，具有較高的精度。

[1]YAMADA C,MORINO S,KAWAGUCHI J,et al.Creep Behavior of Concrete-Filled Steel Tubular Members[R].Mie,Japan:Mie University,1995:83-98.

[2]ICHINOSE L H,WATANABLE E,NAKAI H.An Experimental Study on Creep of Concrete Filled Steel Pipes[J].Journal of Constructional Steel Research,2001,57(4):453-466.

[3]HAN L H,TAO Z,LIU W.Effects of Sustained Load on Concrete-filled Hollow Structural Steel Columns[J].Journal of Structural Engineering,2004,130(2):1392-1404.

[4]HAN L H,YANG Y F.Analysis of Thin-walled Steel RHS Columns Filled with Concrete under Long-term Sustained Loads[J].Thin-Walled Structures,2003,41(9):849-870.

[5]YANG Y F.Behaviour of Recycled Aggregate Concrete Filled Steel Tubular Columns under Long-Term Sustained Loads[J].Advances in Structural Engineering,2011,14(2):189-207.

[6]UY B,DAS S.Time Effects in Concrete-Filled Steel Box Columns in Tall Buildings[J].Structural Design of Tall Buildings,1997,6(1):1-22.

[7]韓冰，王元豐.鋼管混凝土軸心受壓短柱的徐變分析[J].鐵道學報,1999，21(6):87-90.

HAN Bing,WANG Yuan-feng.Study on Creep of Concrete Filled in Steel Tubular Members under Axial Loading[J].Journal of the China Railway Society,1999,21(6):87-90.

[8]王元豐，韓冰.徐變對鋼管混凝土軸心受壓短柱緊箍應力影響分析[J].鐵道學報,2000,22(增1):92-94.WANG Yuan-feng,HAN Bing.Influence of Creep on Confining Stress of Axially Compressed Concrete Filled Steel Tubular Short Columns[J].Journal of the China Railway Society,2000,22(S1):92-94.

[9]WANG Y Y,GENG Y,RANZI G,et al.Time-dependent Behaviour of Expansive Concrete-filled Steel Tubular Columns[J].Journal of Constructional Steel Research,2011,67(3):471-483.

[10]王玉銀,耿悅,張素梅.鋼管微膨脹混凝土軸壓短柱長期變形研究[J].中國公路學報,2011,24(6):57-63.WANG Yu-yin,GENG Yue,ZHANG Su-mei.Research on Long-term Deformation of Concrete-filled Steel Tubular Stubs with Expansive Additive under Axial Loading[J].China Journal of Highway and Transport,2011,24(6):57-63.

[11] 趙金鋼，趙人達，占玉林.鋼管混凝土軸心受壓構件徐變計算方法及徐變模型對比分析[J].公路交通科技,2013,30(4):46-52.

ZHAO Jin-gang,ZHAO Ren-da,ZHAN Yu-lin.Comparative Analysis of Creep Calculation Methods and Creep Models for Axially Compressed CFST Members[J].Journal of Highway and Transportation Research and Development,2013,30(4):46-52.

[12]NAGUIB W,MIRMIRAN A.Creep Modeling for Concrete-filled Steel Tubes[J].Journal of Constructional Steel Research,2003,59 (11):1327-1344.

[13] 韓冰,王元豐.圓鋼管混凝土軸心受壓構件徐變分析的比較[J].中國公路學報,2007,20(2):83-86.

HAN Bing,WANG Yuan-feng.Creep Analysis Comparison of Circular Axially Compressed Concrete-filled Steel Tubular Members[J].China Journal of Highway and Transport,2007,20(2):83-86.

[14]GENG Y，RENZIG WANG Y，et al.Time- dependent Behaviour of Concrete-filled Steel Tubular Columns:Analytical and Comparative Study[J].Magazine of Concrete Research,2012,64 (1):55-69.

[15]American Concrete Institute.ACI Committee 209，Prediction of Creep,Shrinkage and Temperature Effects in Concrete Structures[S].Detroit:American Concrete Institute,1992.

[16]Comite Euro-International Du Beton,CEB-FIP Model Code 1990[S].Lausanne:Thomas Thelford,1993.

[17]GARDNER N J,LOCKMAN M J.Design Provisions for Drying Creep and Shrinkage of Normal Strength Concrete[J].ACI Materials Journal,2001,98(2):159- 167.

[18]BS EN 1992-1-1:2004,Eurocode 2:Design of Concrete Structures [S]

[19]UY B.Static Long-Term Effects in Short Concrete-filled Steel Box Columns under Sustained Loading[J].ACI Structural Journal,2001,98(1) :96-104.

[20]BAZANT Z P.Prediction of Concrete Creep Effects Using Age:Adjusted Effective Modulus Method[J].ACI Journal &Proceedings,1972,69:212-217.

[21]王戎令，王起才，馬麗娜，等.膨脹劑摻量和應力比對鋼管混凝土徐變性能的影響[J].中南大學學報：自然科學版,2014,45(7).2416-2423.

WANG Rong-ling,WANG Qi-cai,MA Li-na,et al.Effect of Interaction between Expanding Agent Proportion and Stress Ratio on Creep Characteristics of Concrete Filled Steel Tube[J].Journal of Central South University：Science and Technology Edition,2014,45(7):2416-2423.

[22]KWON S H，KIM T H，KIM J K，et al.Long-term Behaviour of Square Concrete-Filled Steel Tubular Columns under Axial Service Loads[J].Magazine of Concrete Research,2007,59 (1) :53-68．

[23]HOWEL R W,LARK R J,BARR B I G.A Sensitivity Study of Parameters Used in Shrinkage and Creep Prediction Models[J].Magazine of Concrete Research,2005,57(10):589-602.

Study on Simplified Creep Calculation Method of CFST Members under Axial Loading

WANG Yong-bao,ZHAO Ren-da,XU Teng-fei,ZHAN Yu-lin

(School of Civil Engineering，Southwest Jiaotong University，Chengdu Sichuan 610031，China)

In order to calculate long-term deformation of concrete-filled steel tube (CFST) members properly and rapidly,based on the age adjusted effective modulus method and step-by-step integration method,a simplified analytical method to predict the creep of CFST under axial loading is proposed.The advantages and disadvantages of simplified and traditional methods are comparatively analyzed.The effect of creep prediction model of concrete and its relative humidity value on the calculation accuracy of long-term deformation of CFST axially compressed members are discussed,and the analytical result and test result are compared.The research result shows that (1) the simplified method is simple and has high accuracy,it could be used to analyze the long-term deformation of CFST axially compressed members;(2) using the CEB90 and EC2 creep models with relative humidity of 90% to 98%,the long-term behavior of the CFST members can be well predicted.

bridge engineering;CFST;step-by-step integration method;axially compressed member;creep models;humidity

2015-06-04

國家自然科學基金項目(51208431);中央高校基本科研業(yè)務費專項資金(SWJTU12CX064)

王永寶(1989-)，男，山東濰坊人，博士研究生.(yonbaowang@163.com).

10.3969/j.issn.1002-0268.2016.01.009

U448.34

A

1002-0268(2016)01-0057-07