吳子涵，殷 鳴，劉 智

(1.江西省水利科學(xué)院，江西 南昌 330029；2. 南昌工程學(xué)院，江西 南昌 330099；3. 長(zhǎng)江勘測(cè)規(guī)劃設(shè)計(jì)研究有限責(zé)任公司，湖北 武漢 430010)

鋼筋混凝土是水利工程建設(shè)的重要材料，然而長(zhǎng)期處于干濕交替的工作環(huán)境導(dǎo)致鋼筋極易出現(xiàn)銹蝕現(xiàn)象，是當(dāng)前鋼筋混凝土結(jié)構(gòu)耐久性失效的重要原因。鋼筋銹蝕會(huì)造成鋼筋力學(xué)性能的削弱，與此同時(shí)鋼筋與混凝土之間將出現(xiàn)一層銹蝕產(chǎn)物，影響鋼筋混凝土之間的黏結(jié)滑移性能，并產(chǎn)生銹脹力增大混凝土保護(hù)層的內(nèi)拉應(yīng)力，致使結(jié)構(gòu)發(fā)生銹脹破壞。當(dāng)前學(xué)者關(guān)于鋼筋銹蝕的研究主要集中于，試驗(yàn)提取銹蝕對(duì)鋼筋混凝土材料力學(xué)性能的變化或基于非均勻分布理論的非均勻銹脹作用的研究等[1-5]，少有將兩者結(jié)合起來(lái)探究鋼筋銹蝕對(duì)結(jié)構(gòu)的影響，并缺少相應(yīng)的理論仿真方法。為此，本文在已有鋼筋混凝土分離式模型和混凝土損傷模型的基礎(chǔ)上，結(jié)合鋼筋銹蝕對(duì)鋼筋混凝土材料性能的影響，重點(diǎn)分析鋼筋銹蝕產(chǎn)生的銹脹作用以及對(duì)梁構(gòu)件抗彎性能的影響，為后續(xù)進(jìn)一步研究鋼筋銹蝕對(duì)水工鋼筋混凝土結(jié)構(gòu)安全的影響提供理論基礎(chǔ)。

1 鋼筋銹蝕的影響

根據(jù)已有理論及試驗(yàn)研究，當(dāng)前可將鋼筋銹蝕對(duì)鋼筋混凝土材料性能的影響分為3個(gè)方面。

1.1 鋼筋銹蝕對(duì)鋼筋力學(xué)性能的影響

鋼筋銹蝕后，材料表面將逐漸出現(xiàn)大小不一的銹坑，從而降低鋼筋自身強(qiáng)度及其變形能力。當(dāng)前國(guó)內(nèi)外學(xué)者的研究成果中，銹蝕引起的鋼筋力學(xué)參數(shù)變化趨勢(shì)基本一致，但受限于試驗(yàn)本身的復(fù)雜性與結(jié)果的離散性，銹蝕率與力學(xué)性能變化的定量關(guān)系差異較大，尚未出現(xiàn)較為統(tǒng)一的結(jié)論。本文采用吳慶等[1]的試驗(yàn)結(jié)果，將銹蝕率為5%作為模型的界限；該結(jié)論在CECS 220:2007《混凝土結(jié)構(gòu)耐久性評(píng)定標(biāo)準(zhǔn)》[6]中也有體現(xiàn)，標(biāo)準(zhǔn)認(rèn)為當(dāng)鋼筋銹蝕率達(dá)到5%時(shí)，構(gòu)件已經(jīng)屬于嚴(yán)重銹損構(gòu)件，因此本文僅研究銹蝕率ρ小于5%的情況。

fy，c=fy0(1-0.029ρ)
fu，c=fu0(1-0.026ρ)
δs，c=δs0(1-0.0575ρ)
εy，c=εy0(1-0.0575ρ)
Eu，c=Eu0(1-0.052ρ)

(1)

式中，Eu，c、fy，c、fu，c、δs，c、εy，c分別為銹蝕鋼筋的名義彈性模量、名義屈服強(qiáng)度、名義極限強(qiáng)度、名義伸長(zhǎng)率和極限應(yīng)變；Eu0、fy0、fu0、δs0、εy0分別為未銹蝕鋼筋的彈性模量、屈服強(qiáng)度、極限強(qiáng)度、名義伸長(zhǎng)率和極限應(yīng)變。

1.2 鋼筋銹蝕對(duì)鋼筋混凝土粘結(jié)滑移性能的影響

鋼筋銹蝕后將在銹蝕表面產(chǎn)生一層零碎的銹蝕產(chǎn)物，此時(shí)受損混凝土與鋼筋銹蝕產(chǎn)物將在材料的膠結(jié)部位形成新的粘結(jié)界面，該界面的力學(xué)性能受混凝土強(qiáng)度、保護(hù)層厚度以及銹蝕產(chǎn)物厚度等多方面因素的影響。本文鋼筋銹蝕后鋼筋混凝土極限粘結(jié)強(qiáng)度采用曾宇[2]論文中的試驗(yàn)擬合結(jié)果，其研究發(fā)現(xiàn)，常規(guī)分析模型與試驗(yàn)結(jié)果相差較大的原因是未考慮鋼筋銹蝕率與各影響因素間的耦合效應(yīng)，因此其建立的公式包含了與銹蝕產(chǎn)物厚度、保護(hù)層厚度和抗拉強(qiáng)度之間的關(guān)系，即

τcu=βuτou=K1K2K3τou

(2)

式中，τcu為鋼筋銹蝕后材料極限粘結(jié)強(qiáng)度；βu為材料極限粘結(jié)強(qiáng)度折減系數(shù)；τou為材料原始極限粘結(jié)強(qiáng)度；K1為鋼筋銹蝕產(chǎn)物厚度t的影響參數(shù)；K2為相對(duì)保護(hù)層厚度c/d的影響參數(shù)；K3為混凝土抗拉強(qiáng)度的影響參數(shù)。

根據(jù)其試驗(yàn)擬合，可得

K1=-0.058 7t+0.011 3

K2=-0.404 1(c/d)2+2.673 4(c/d)+0.252 6

K3=0.352 6ft+2.984 7

(3)

式中，t可根據(jù)GB/T 51355—2019《既有混凝土結(jié)構(gòu)耐久性評(píng)定標(biāo)準(zhǔn)》[7]5.2.1節(jié)中公式求得。

鋼筋銹蝕后的鋼筋混凝土粘結(jié)強(qiáng)度折減系數(shù)則采用徐善華[8]論文中根據(jù)滕智明試驗(yàn)結(jié)果回歸擬合的公式

β=1+0.562 5ρ-0.337 5ρ2+0.055 625ρ3-0.003ρ4

(4)

式中，ρ為鋼筋銹蝕率；β為銹后鋼筋與混凝土之間的粘結(jié)強(qiáng)度降低系數(shù)。

1.3 銹脹作用

通過(guò)數(shù)值模擬分析鋼筋銹脹作用能夠得到混凝土應(yīng)力分布狀態(tài)和裂縫潛在發(fā)展路徑等試驗(yàn)難以獲取的成果。當(dāng)前學(xué)者們采用有限元法對(duì)銹脹作用的模擬主要有施加徑向力，施加徑向位移和設(shè)置溫度膨脹環(huán)等3種方式[9]。無(wú)論采用何種方法模擬，均需要充分考慮銹蝕非均勻的基本特性，即實(shí)際工程中鋼筋毗鄰保護(hù)層的一側(cè)總是先發(fā)生銹蝕，且銹蝕速度也遠(yuǎn)高于鋼筋內(nèi)側(cè)。其中，施加徑向位移的方法相對(duì)而言使用較為廣泛[3]，本文采用殷鳴[9]論文中引用的銹蝕鋼筋周?chē)炷翉较蛭灰品蔷鶆蚍植寄Ｐ?，如圖1所示。

圖1 混凝土徑向位移分布

非角區(qū)鋼筋周?chē)炷翉较蛭灰品植家?guī)律為

(5)

角區(qū)鋼筋周?chē)炷翉较蛭灰品植家?guī)律為

(6)

式中，r為鋼筋半徑；u1為外側(cè)最大位移；u2為內(nèi)側(cè)最小位移；θ為各部位與x軸的夾角。

2 鋼筋混凝土材料仿真分析模型

2.1 鋼筋混凝土分離式模型

當(dāng)前主流的鋼筋混凝土材料有限元仿真分析模型主要可分為嵌入式模型、分布式模型和分離式模型，其中嵌入式模型和分布式模型假定鋼筋與混凝土之間不產(chǎn)生相對(duì)滑移，主觀消除了鋼筋位移自由度的運(yùn)算，因此雖然計(jì)算效率較高，但無(wú)法準(zhǔn)確反映鋼筋混凝土的相互作用。本文基于趙蘭浩等[10]建立的單彈簧聯(lián)結(jié)單元法，結(jié)合GB 50010—2010《混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范》[11]所采用的徐有鄰五段線(xiàn)模型構(gòu)建鋼筋混凝土分離式模型來(lái)模擬鋼筋與混凝土之間的相互作用[12]。

單彈簧聯(lián)結(jié)單元法如圖2所示。該方法基于混合坐標(biāo)系建立，采用實(shí)體單元模擬混凝土，使用桿單元模擬鋼筋，并在混凝土與鋼筋接觸部位設(shè)置重復(fù)節(jié)點(diǎn)分別作為混凝土與鋼筋插值點(diǎn)，建立彈簧單元作為兩者間的聯(lián)結(jié)單元。通過(guò)法向約束強(qiáng)制混凝土單元與鋼筋單元間的法向自由度協(xié)同，通過(guò)設(shè)置切向彈簧模擬兩者之間的切向作用，切向剛度則通過(guò)粘結(jié)-滑移本構(gòu)關(guān)系確定，其有限元平衡方程為

圖2 單彈簧聯(lián)結(jié)單元法求解示意

(7)

正確的鋼筋混凝土粘結(jié)-滑移本構(gòu)模型是準(zhǔn)確模擬鋼筋混凝土結(jié)構(gòu)性態(tài)的重要前提，本文采用的徐有鄰五段線(xiàn)模型可分為膠結(jié)階段、強(qiáng)度上升階段、強(qiáng)度峰值階段、強(qiáng)度退化階段和摩擦穩(wěn)定階段(見(jiàn)圖3)，與鋼筋混凝土相互作用機(jī)理吻合。其相應(yīng)計(jì)算公式為

圖3 鋼筋混凝土粘結(jié)-滑移本構(gòu)關(guān)系示意

(8)

式中，τ，τs，τcr，τu，τr分別為粘結(jié)應(yīng)力及各階段粘結(jié)強(qiáng)度特征值；s，ss，scr，su，sr分別為相對(duì)滑移量及各階段滑移量特征值；k1，k2，k3，k4分別為膠著段、上升段、峰值段、退化段的斜率。

模型骨架線(xiàn)各階段粘結(jié)強(qiáng)度與滑移量特征值見(jiàn)表1。

表1 鋼筋混凝土粘結(jié)-滑移骨架線(xiàn)粘結(jié)強(qiáng)度與滑移量特征值

2.2 混凝土損傷模型

本文采用韋未等[13]基于四參數(shù)破壞準(zhǔn)則的基本思路提出的四參數(shù)等效應(yīng)變計(jì)算方法，通過(guò)式(10)的等效應(yīng)變求解可將真實(shí)空間中復(fù)雜的多軸問(wèn)題轉(zhuǎn)化為等效空間中簡(jiǎn)單的單軸問(wèn)題。

(9)

(10)

本文損傷模型選取過(guò)鎮(zhèn)海提出的受拉與受壓應(yīng)力應(yīng)變?nèi)€(xiàn)，該曲線(xiàn)已得到國(guó)內(nèi)外科研工作者的認(rèn)可，并納入我國(guó)混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范[11]。

混凝土單軸受拉應(yīng)力應(yīng)變曲線(xiàn)表達(dá)式為

σ=(1-dt)Ecε

(11)

(12)

(13)

式中，σ和ε分別為混凝土的應(yīng)力與應(yīng)變；ft，r為混凝土單軸抗拉強(qiáng)度；εt，r為ft，r對(duì)應(yīng)的應(yīng)變；dt為單軸受拉損傷變量；Ec為混凝土彈性模量；at為混凝土受拉應(yīng)力應(yīng)變曲線(xiàn)軟化段參數(shù)。

混凝土單軸受壓應(yīng)力應(yīng)變曲線(xiàn)表達(dá)式為

σ=(1-dc)Ecε

(14)

(15)

(16)

(17)

式中，fc，r為混凝土單軸抗壓強(qiáng)度；εc，r為fc，r對(duì)應(yīng)的應(yīng)變；dc為單軸受壓損傷變量；ac為混凝土受壓應(yīng)力應(yīng)變曲線(xiàn)軟化段參數(shù)。

當(dāng)結(jié)構(gòu)處于多軸狀態(tài)時(shí)，需用等效應(yīng)力應(yīng)變?chǔ)?和ε*替代真實(shí)應(yīng)力應(yīng)變和ε，通過(guò)應(yīng)變張量第一不變量I′1判斷單元拉壓狀態(tài)，該模型正確性已在文獻(xiàn)[14]中得以驗(yàn)證。

3 算例驗(yàn)證

3.1 銹脹作用算例

截取文獻(xiàn)[15]中梁試件的橫截面建立二維有限元模型，如圖4所示。模型尺寸為100 mm×160 mm，在布置鋼筋處設(shè)立與鋼筋截面相對(duì)應(yīng)的孔洞，保護(hù)層厚15 mm。本算例僅考慮梁底部的受拉鋼筋和頂部的受壓鋼筋的銹脹作用。

圖4 梁試件二維模型示意(單位：mm)

本算例材料參數(shù)如表2所示，推求損傷模型四參數(shù)取值為：A=0.013 47，B=0.131 49，C=0.742，D=0.245 19。本文銹脹加載取角區(qū)最大徑向位移u1為15 μm，u1/u2=30，將模型鋼筋孔洞處各節(jié)點(diǎn)位移線(xiàn)性逐步增大，以模擬銹脹的過(guò)程。

表2 算例材料參數(shù)

根據(jù)混凝土損傷演化發(fā)展，整個(gè)銹脹過(guò)程大致可分為4個(gè)階段，如圖5所示。從圖5可知：①內(nèi)裂期。銹蝕發(fā)生后將首先在面向角區(qū)方向出現(xiàn)開(kāi)裂，形成3條明顯的損傷演化路徑，此時(shí)鋼筋附近混凝土將不斷發(fā)生應(yīng)力重分布，內(nèi)部隨鋼筋銹脹發(fā)展沿?fù)p傷路徑持續(xù)演化。②擴(kuò)展期。隨著銹蝕產(chǎn)物不斷增加，混凝土損傷沿初始路徑向外部發(fā)展，此時(shí)4個(gè)角區(qū)部位的混凝土損傷沿兩條相互垂直路徑迅速演化，主拉應(yīng)力區(qū)不斷向外部轉(zhuǎn)移，混凝土保護(hù)層外表面同時(shí)出現(xiàn)受拉破壞并向內(nèi)部延申。③貫通期。隨著銹蝕程度進(jìn)一步增加，混凝土內(nèi)部裂縫與外部裂縫形成貫通，此時(shí)裂縫貫通路徑上的應(yīng)力完全釋放轉(zhuǎn)移，混凝土表層出現(xiàn)脫落。由于上部鋼筋直徑較下部鋼筋偏小，因此上部區(qū)域的混凝土保護(hù)層出現(xiàn)損傷貫通較下部更遲。④深化期。雖然混凝土保護(hù)層已形成貫穿裂縫，但銹脹作用仍隨銹蝕程度加深而繼續(xù)，此時(shí)混凝土損傷程度不斷加深，混凝土表層脫落現(xiàn)象加劇，同時(shí)裂縫逐步向內(nèi)部發(fā)展。

圖5 銹脹作用演化發(fā)展過(guò)程

3.2 四點(diǎn)彎曲梁算例

選取Walraven四點(diǎn)彎曲梁試驗(yàn)作為研究對(duì)象，研究不同銹蝕率對(duì)受彎梁試件的影響。模型如圖6所示。模型尺寸為4.1 m×0.2 m×0.45 m，在底部40 mm處布置3根不同直徑鋼筋；模型頂部施加對(duì)稱(chēng)的分級(jí)荷載，在右側(cè)支點(diǎn)施加法向約束，左側(cè)支

圖6 四點(diǎn)彎曲梁試驗(yàn)?zāi)Ｐ褪疽?單位：mm)

點(diǎn)施加雙向約束。模型混凝土材料參數(shù)如表3所示。

表3 算例材料參數(shù)

對(duì)受彎梁構(gòu)件處于銹蝕率0～5%之間的6種情況進(jìn)行計(jì)算，計(jì)算工況如表4所示。

表4 算例材料參數(shù)

計(jì)算所得梁試件的損傷演化過(guò)程如圖7所示，由于各工況下?lián)p傷演化過(guò)程相似，因此僅列演化過(guò)程的示意圖。由圖7可見(jiàn)，其過(guò)程與常規(guī)試驗(yàn)規(guī)律一致。不同銹蝕率下梁試件的荷載-擾度曲線(xiàn)如圖8所示。由圖8可見(jiàn)，在各種銹蝕率下，混凝土梁開(kāi)裂前的荷載-撓度關(guān)系均呈線(xiàn)性變化，各工況下數(shù)值基本一致。當(dāng)梁在跨中底部出現(xiàn)開(kāi)裂現(xiàn)象后，荷載-撓度曲線(xiàn)斜率較開(kāi)裂前減小，此時(shí)不同銹蝕率的荷載-撓度曲線(xiàn)存在較為明顯的區(qū)別，底部鋼筋銹蝕導(dǎo)致粘結(jié)滑移曲線(xiàn)中的化學(xué)膠著段消失，極限粘結(jié)強(qiáng)度也有所降低，因此梁的抗彎能力將隨著銹蝕率增大而明顯降低。

圖7 四點(diǎn)彎曲梁試件損傷演化過(guò)程

圖8 不同銹蝕率下梁跨中荷載-撓度曲線(xiàn)對(duì)比

對(duì)于梁受彎構(gòu)件數(shù)值仿真而言，通常計(jì)算的極限荷載會(huì)略低于試驗(yàn)數(shù)據(jù)，主要原因有：①試驗(yàn)結(jié)果受試驗(yàn)條件和試驗(yàn)過(guò)程影響，存在一定的離散性；②試驗(yàn)推求的材料本構(gòu)模型未能充分考慮鋼筋的強(qiáng)化階段；③數(shù)值計(jì)算中所考慮的不均勻銹蝕與真實(shí)情況存在偏差，全截面銹蝕相較真實(shí)情況更為嚴(yán)重，導(dǎo)致結(jié)果偏大。

4 結(jié) 語(yǔ)

本文針對(duì)鋼筋混凝土結(jié)構(gòu)無(wú)法規(guī)避的鋼筋銹蝕問(wèn)題，梳理了鋼筋銹蝕對(duì)鋼筋混凝土材料性能的影響因素，并根據(jù)相關(guān)研究成果篩選出可靠的影響因素?cái)?shù)值計(jì)算方法。結(jié)合已有的鋼筋混凝土分離式模型和混凝土損傷模型，通過(guò)算例分析鋼筋銹蝕后的銹脹作用以及對(duì)混凝土受彎構(gòu)件抗彎性能的影響。計(jì)算結(jié)果能夠較為準(zhǔn)確地描述鋼筋銹蝕后鋼筋混凝土材料性能的變化，為進(jìn)一步研究鋼筋銹蝕對(duì)水工鋼筋混凝土結(jié)構(gòu)安全的影響奠定了良好基礎(chǔ)。