聶良學(xué)，許金余，劉志群，羅 鑫

（1．空軍工程大學(xué)機場建筑工程系，陜西西安710038；2．解放軍理工大學(xué)爆炸沖擊防災(zāi)減災(zāi)國家重點實驗室，江蘇南京210007；3．西北工業(yè)大學(xué)力學(xué)與土木建筑學(xué)院，陜西西安710072；4．中央軍委后勤保障部建筑工程規(guī)劃設(shè)計研究院，北京100062）

鹽腐蝕后混凝土的動態(tài)本構(gòu)模型＊

聶良學(xué)1，2，許金余1，3，劉志群1，羅 鑫4

（1．空軍工程大學(xué)機場建筑工程系，陜西西安710038；2．解放軍理工大學(xué)爆炸沖擊防災(zāi)減災(zāi)國家重點實驗室，江蘇南京210007；3．西北工業(yè)大學(xué)力學(xué)與土木建筑學(xué)院，陜西西安710072；4．中央軍委后勤保障部建筑工程規(guī)劃設(shè)計研究院，北京100062）

為探究混凝土受鹽腐蝕后的動態(tài)力學(xué)響應(yīng)，配置了粉煤灰質(zhì)量分?jǐn)?shù)為15%的普通硅酸鹽水泥混凝土，將其置于質(zhì)量分?jǐn)?shù)均為15%的NaCl和Na2SO4溶液中浸泡腐蝕60d后，利用 100mm分離式霍普金森壓桿實驗裝置，測試其受腐蝕后的動態(tài)力學(xué)性能，并結(jié)合宏觀唯象損傷統(tǒng)計理論和Weibull分布思想，建立了混凝土受鹽腐蝕后的動態(tài)統(tǒng)計損傷本構(gòu)模型。結(jié)果表明：受鹽腐蝕后，混凝土試件的動態(tài)抗壓強度均有不同程度的下降，且NaCl溶液腐蝕試件的降幅大于Na2SO4溶液腐蝕試件；模型曲線與實驗曲線的擬合度較高，能夠較準(zhǔn)確地描述混凝土在沖擊荷載作用下的動態(tài)力學(xué)響應(yīng)規(guī)律。

鹽腐蝕；混凝土；動態(tài)抗壓強度；損傷理論；Weibull分布；本構(gòu)模型

混凝土的耐久性直接關(guān)系到建筑結(jié)構(gòu)的壽命及安全性和穩(wěn)定性。在威脅混凝土結(jié)構(gòu)耐久性的諸多因素中，鹽腐蝕尤為突出。鹽腐蝕多發(fā)生于濱海地區(qū)和鹽漬土地區(qū)。地處濱海區(qū)域的建筑或直接建于海水中的混凝土結(jié)構(gòu)，常年經(jīng)受海水的浸泡和沖刷，遭受海水中Cl－和SO2－4的侵蝕；在我國西部鹽漬土區(qū)域，如新疆、甘肅、青海、寧夏等地，土壤中SO2－4的質(zhì)量濃度最高可達(dá)4 200mg／L，高出我國沿海海水中SO2－4約1倍，處于該地區(qū)的建筑，尤其是地下混凝土結(jié)構(gòu)，受鹽離子的侵蝕十分嚴(yán)重。Cl－對混凝土保護層的腐蝕較為嚴(yán)重，可致其脫落，引起內(nèi)部鋼筋銹蝕；硫酸鹽侵入混凝土內(nèi)部引發(fā)石膏腐蝕和鈣礬石腐蝕，導(dǎo)致混凝土內(nèi)部結(jié)構(gòu)發(fā)生膨脹破壞，明顯弱化其各項靜／動力學(xué)性能。雖然目前關(guān)于鹽對混凝土的腐蝕研究較多，但是大多集中于靜態(tài)力學(xué)性能［15］、離子擴散［68］及靜態(tài)損傷模型［911］等，對其在高應(yīng)變率下的性能測試及動態(tài)本構(gòu)模型仍需深入研究。

基于此，本文中利用 100mm分離式霍普金森壓桿（split Hopkinson pressure bar，SHPB）實驗裝置，對受氯鹽和硫酸鹽腐蝕的混凝土進(jìn)行動態(tài)力學(xué)性能研究，結(jié)合宏觀唯象損傷統(tǒng)計理論和Weibull分布思想，建立混凝土受鹽腐蝕后的動態(tài)統(tǒng)計損傷本構(gòu)模型，并對模型進(jìn)行驗證。

1 實 驗

1．1 材料和試件

制作混凝土試件所用的材料包括：42．5R級普通硅酸鹽水泥，由硅酸鹽水泥熟料、質(zhì)量分?jǐn)?shù)為6%～20%的少量氧化物、適量磨細(xì)石膏等材料組成；石灰?guī)r碎石，堆積密度為1．62kg／L，粒徑范圍為5～20mm；灞河中砂，堆積密度為1．5kg／L，細(xì)度模數(shù)為2．8；自來水；F類（低鈣）Ⅰ級粉煤灰?；炷恋脑O(shè)計強度為C30，試件尺寸為 100mm×50mm，配合比如表1所示。

表1 混凝土配合比Table 1 Mix proportions of concrete kg／m3

試件澆筑完畢后在室內(nèi)靜置24h后脫模，然后立即置于養(yǎng)護箱內(nèi)，進(jìn)行為期28d的標(biāo)準(zhǔn)養(yǎng)護。

1．2 動態(tài)壓縮實驗

本實驗設(shè)置兩組腐蝕環(huán)境：氯鹽和硫酸鹽。將養(yǎng)護后的試件分別置于質(zhì)量分?jǐn)?shù)均為15%的NaCl和Na2SO4溶液中浸泡，腐蝕60d后，利用 100mm的SHPB實驗裝置測試其在不同應(yīng)變率下的動態(tài)力學(xué)性能。同時，設(shè)置一組空白對照組，該組試件不浸泡任何溶液，其余外部環(huán)境均相同。為便于分析，將浸泡于NaCl溶液和Na2SO4溶液中的試件分別標(biāo)記為S1組、S2組，對照組標(biāo)記為N組。

2 實驗結(jié)果

2．1 應(yīng)力－應(yīng)變曲線

在不同應(yīng)變率下，各組試件的動態(tài)應(yīng)力－應(yīng)變（σ－ε）曲線如圖1所示。由圖1可知：隨著應(yīng)變率的升高，峰值應(yīng)力逐漸增加，具有顯著的應(yīng)變率效應(yīng)；同時應(yīng)力－應(yīng)變曲線的下降段總體上隨著應(yīng)變率的升高而逐漸向高應(yīng)變方向移動，說明試件的變形程度不斷增大。

圖1 3組混凝土試樣的動態(tài)應(yīng)力－應(yīng)變曲線Fig．1 Dynamic stress－strain curves of three groups of concrete specimens

2．2 動態(tài)抗壓強度

動態(tài)抗壓強度（即峰值應(yīng)力）為試件破壞時的應(yīng)力，是衡量材料在動荷載作用下的強度規(guī)律的重要力學(xué)性能指標(biāo)，直接反映材料的動力特性。圖2為各組混凝土試件在不同沖擊荷載作用下動態(tài)抗壓強度fc，d與平均應(yīng)變率的關(guān)系圖。觀察可知，fc，d與滿足線性規(guī)律，對其進(jìn)行線性擬合，可以得到如下關(guān)系式：

由圖2以及式（1）～式（3）可知：隨著應(yīng)變率的升高，各組試件的fc，d明顯增大，且N組試件的fc，d的增長率（擬合直線斜率）最大，即其應(yīng)變率敏感性最強，S1組試件次之，S2組試件最?。幌啾扔贜組試件，S1組和S2組試件的fc，d均有不同程度的下降，且S1組試件的降幅大于S2組。

測得N組、S1組和S2組混凝土試件的靜態(tài)抗壓強度（fc，s）分別為42．65、32．51和35．33MPa。定義動態(tài)強度增長因子I：動態(tài)抗壓強度fc，d與靜態(tài)抗壓強度fc，s的比值，即I＝fc，d／fc，s。圖3為I與應(yīng)變率對數(shù)的關(guān)系圖，對其進(jìn)行線性擬合，得到：

圖2 動態(tài)抗壓強度－平均應(yīng)變率曲線Fig．2 Dynamic compressive strength vs．a(chǎn)verage strain rate

圖3 動態(tài)強度增長因子－平均應(yīng)變率對數(shù)曲線Fig．3 Dynamic increase factor vs．the logarithm of average strain rate

觀察圖3可知：對于3組試件，其動態(tài)強度增長因子I的線性相關(guān)性均較強；N組試件的I最??；S1組試件的I在較低的應(yīng)變率下略低于S2組，而在較高的應(yīng)變率下高于S2組。

試件被NaCl溶液腐蝕后，Cl－與混凝土內(nèi)部晶體或膠凝體發(fā)生一系列化學(xué)反應(yīng)，破壞混凝土內(nèi)部物質(zhì)組成和結(jié)構(gòu)，改變了原有水化產(chǎn)物的性質(zhì)，導(dǎo)致水化硅酸鈣膠凝體與水泥石的黏結(jié)性能下降，致使混凝土強度降低。Na2SO4溶液滲入混凝土內(nèi)部后，同樣產(chǎn)生結(jié)晶體，填充于混凝土內(nèi)部孔隙，當(dāng)結(jié)晶體累積到一定量后，對混凝土內(nèi)部進(jìn)行擠壓，產(chǎn)生內(nèi)應(yīng)力；同時，當(dāng)SO2－4進(jìn)入混凝土內(nèi)部后，與氫氧化鈣和水化鋁酸鈣發(fā)生反應(yīng)，生成鈣礬石和石膏，反應(yīng)后體積明顯增大，亦產(chǎn)生內(nèi)應(yīng)力，嚴(yán)重時引發(fā)混凝土產(chǎn)生由內(nèi)向外的裂縫，破壞骨料與水泥石的黏結(jié)，降低混凝土強度。圖4顯示了通過掃描電鏡（scanning electron microscope，SEM）觀測的各組試件微觀形貌。從圖4可以看出：N組試件因未受任何腐蝕，隨著時間的增長，水化產(chǎn)物C─S─H凝膠逐漸增多，致密性較大，微觀缺陷和損傷較小，強度持續(xù)增加；對于被鹽腐蝕的試件（S1組和S2組），其表面的絮狀結(jié)晶體較多，且排列緊密，形成一層致密的晶體層，使裂縫內(nèi)部產(chǎn)生一定擠壓應(yīng)力，致使試件的抗壓強度較N組試件有所下降。

圖4 3組混凝土試樣的掃描電鏡圖像Fig．4 SEM images of three groups of concrete specimens

3 動態(tài)統(tǒng)計損傷本構(gòu)模型

3．1 本構(gòu)模型的構(gòu)建

由于混凝土內(nèi)部存在微裂紋、微孔洞等初始損傷，且裂紋密度和應(yīng)變均服從Weibull分布，因此將混凝土視為連續(xù)介質(zhì)，結(jié)合連續(xù)損傷力學(xué)，基于統(tǒng)計理論，建立統(tǒng)計損傷本構(gòu)模型，即：

式中：E為混凝土無損狀態(tài)下的彈性模量，D為損傷變量。

混凝土內(nèi)部微裂紋密度和應(yīng)變均服從Weibull分布，概率密度函數(shù)φ（ε）的表達(dá)式為［12］：

式中：a、m是與材料物理力學(xué)性質(zhì)相關(guān)的系數(shù)。

3．2 模型參數(shù)的確定

根據(jù)連續(xù)損傷力學(xué)理論［13］，可將混凝土視為多個“微元體”的組合，同時假定每個“微元體”所受的應(yīng)力和損傷參量連續(xù)且均勻分布。設(shè)混凝土內(nèi)部所有“微元體”的數(shù)量為N，受荷載作用后被破壞的“微元體”數(shù)為n，則統(tǒng)計損傷變量D定義為：

因此，在某一應(yīng)變區(qū)間［ε，ε＋dε］內(nèi)，已破壞的“微元體”數(shù)為Nφ（ε）dε。當(dāng)混凝土應(yīng)變?yōu)棣艜r，已破壞的“微元體”數(shù)為［14］：

將式（10）代入式（8），可得混凝土的統(tǒng)計損傷演化方程：

將式（11）代入式（7），則本構(gòu)模型變?yōu)椋?/p>

在混凝土靜態(tài)應(yīng)力－應(yīng)變曲線的上升段選取兩點A、B，所對應(yīng)的應(yīng)力σA、σB分別為峰值強度的40%和60%，則此兩點的斜率即為所需的彈性模量：

式中：εA、εB分別表示A、B兩點的應(yīng)變。根據(jù)式（13），可確定N組、S1組和S2組的彈性模量分別為43．52、35．79和36．62GPa。

因此，只需再確定材料的物理力學(xué)性能參數(shù)a和m，即可完全確定該本構(gòu)模型。由式（12）可得：

對式（14）等號兩邊取兩次對數(shù)，得：

令x＝lnε，y＝ln［－ln（σ／Eε）］，則式（15）變?yōu)椋?/p>

基于各組混凝土試件的動態(tài)應(yīng)力－應(yīng)變曲線，可得本構(gòu)模型參數(shù)a和m，見表2。對a和m分別進(jìn)行非線性擬合，得到其與平均應(yīng)變率的關(guān)系：

將式（17）～式（19）分別代入式（12），即可得到最終的動態(tài)統(tǒng)計損傷本構(gòu)模型。

表2 不同應(yīng)變率下模型參數(shù)Table 2 Model parameters at different strain rates

3．3 應(yīng)變率效應(yīng)修正

表3 動態(tài)強度增長因子與平均應(yīng)變率對數(shù)的關(guān)系Table 3 Relation of dynamic compressive strength increase factor with the logarithm of average strain rate

式中：Kε為應(yīng)變率強化因子，

綜上所述，考慮應(yīng)變率強化效應(yīng)的動態(tài)統(tǒng)計損傷本構(gòu)模型為：最終的動態(tài)統(tǒng)計損傷本構(gòu)模型的表達(dá)式為：

3．4 模型驗證

至此，動態(tài)統(tǒng)計損傷本構(gòu)模型已構(gòu)建完成，將實際的動態(tài)應(yīng)力－應(yīng)變曲線（部分）與計算曲線進(jìn)行對比，如圖5所示?？梢钥闯?，各組試件的應(yīng)力－應(yīng)變實驗曲線與相應(yīng)的模型曲線吻合較好，雖然存在一定的偏差，但整體上能夠較準(zhǔn)確地描述各組混凝土試件在沖擊荷載作用下的動態(tài)力學(xué)響應(yīng)規(guī)律。

圖5 實驗曲線與模型曲線的對比Fig．5 Comparison of experimental and model curves

模型中，將用于確定參數(shù)a、m及C值的動態(tài)應(yīng)力－應(yīng)變曲線記為基準(zhǔn)曲線，因此圖5反映的是基準(zhǔn)實驗曲線與本構(gòu)模型曲線的擬合程度。為進(jìn)一步對模型進(jìn)行驗證，將模型曲線與各組非基準(zhǔn)曲線作進(jìn)一步的對比分析。分別取N組、S1組和S2組的非基準(zhǔn)曲線，對應(yīng)的應(yīng)變率分別為95．79、67．15和82．64s－1，根據(jù)式（17）～式（19）計算得到a、m值，分別為：N組，a＝1．207×10－4，m＝0．183；S1組，a＝7．367×10－4，m＝0．360；S2組，a＝0．774×10－4，m＝0．198。再由表3及式（21）計算得到各應(yīng)變率下的模型曲線，并與實際曲線進(jìn)行對比，如圖6所示。由圖6可知，各組試樣的應(yīng)力－應(yīng)變實驗曲線與模型曲線同樣具有較高的吻合度。由于該本構(gòu)模型是基于統(tǒng)計理論、Weibull分布函數(shù)及部分假定構(gòu)建而成，而試件的應(yīng)力－應(yīng)變行為受試件本身尺寸、內(nèi)部結(jié)構(gòu)及外界環(huán)境等影響，因此模型曲線與實驗曲線之間的偏差在可接受范圍之內(nèi)。

圖6 模型驗證Fig．6 Model verification

4 結(jié) 論

利用 100mm SHPB實驗裝置，對經(jīng)質(zhì)量分?jǐn)?shù)均為15%的NaCl和Na2SO4溶液浸泡腐蝕的混凝土的動態(tài)力學(xué)性能進(jìn)行了研究，建立并驗證了混凝土腐蝕后的動態(tài)統(tǒng)計損傷本構(gòu)模型，主要結(jié)論如下：

（1）各組試件的動態(tài)抗壓強度和動態(tài)強度增長因子均具有較強的應(yīng)變率相關(guān)性，而受鹽腐蝕后，各組試件的動態(tài)抗壓強度均有不同程度的下降，且NaCl腐蝕試件的降幅大于Na2SO4腐蝕試件；

（2）基于宏觀唯象損傷統(tǒng)計理論和Weibull分布思想，將混凝土視為連續(xù)介質(zhì)，以混凝土內(nèi)部“微元體”定義損傷，構(gòu)建動態(tài)統(tǒng)計損傷本構(gòu)模型，各參數(shù)的物理意義明確，參數(shù)a和m精確控制著擬合曲線的形狀、彎曲度和幅值等，同時計入應(yīng)變率強化效應(yīng)修正系數(shù)C；

（3）動態(tài)統(tǒng)計損傷本構(gòu)模型曲線與實驗曲線的擬合度較高，能夠較準(zhǔn)確地描述各組混凝土試件在沖擊荷載作用下的動態(tài)力學(xué)響應(yīng)規(guī)律，模型曲線與實驗曲線之間存在的偏差在可接受范圍之內(nèi)。

［1］ Uysal M，Yilmaz K，Ipek M．The effect of mineral admixtures on mechanical properties，chloride ion permeability and impermeability of self－compacting concrete［J］．Construction ＆Building Materials，2012，27（1）：263－270．

［2］ 袁承斌，張德峰，劉榮桂，等．不同應(yīng)力狀態(tài)下混凝土抗氯離子侵蝕的研究［J］．河海大學(xué)學(xué)報（自然科學(xué)版），2003，31（1）：50－54．Yuan Chengbin，Zhang Defeng，Liu Ronggui，et al．Diffusivity of chloride in concrete in different stress states［J］．Journal of Hehai University（Natural Sciences），2003，31（1）：50－54．

［3］ Hossain K M A，Lachemi M．Performance of volcanic ash and pumice based blended cement concrete in mixed sulfate environment［J］．Cement ＆Concrete Research，2006，36（6）：1123－1133．

［4］ 范穎芳，黃振國，郭樂工，等．硫酸鹽腐蝕后混凝土力學(xué)性能研究［J］．鄭州工業(yè)大學(xué)學(xué)報，1999，20（1）：91－93．Fan Yingfang，Huang Zhenguo，Guo Legong，et al．Research on mechanics property of concrete after corrosion of sulphate［J］．Journal of Zhengzhou University of Technology，1999，20（1）：91－93．

［5］ Chindaprasirt P，Chalee W．Effect of sodium hydroxide concentration on chloride penetration and steel corrosion of fly ash－based geopolymer concrete under marine site［J］．Construction ＆Building Materials，2014，63（8）：303－310．

［6］ Kwon S J，Na U J，Sang S P，et al．Service life prediction of concrete wharves with early－aged crack：Probabilisticapproach for chloride diffusion［J］．Structural Safety，2009，31（1）：75－83．

［7］ 史美倫，張雄，吳科如．混凝土中氯離子滲透性測定的電化學(xué)方法［J］．硅酸鹽通報，1998（6）：55－59．Shi Meilun，Zhang Xiong，Wu Keru．Electrochemical methods for the measurement of chloride permeability in concrete［J］．Bulletin of the Chinese Ceramic Society，1998（6）：55－59．

［8］ Deby F，Carcassès M，Sellier A．Probabilistic approach for durability design of reinforced concrete in marine environment［J］．Cement ＆Concrete Research，2009，39（5）：466－471．

［9］ Sun C，Chen J，Zhu J，et al．A new diffusion model of sulfate ions in concrete［J］．Construction ＆Building Materials，2013，39（1）：39－45．

［10］ Thomas M，Bamforth P B．Modeling chloride diffusion in concrete：Effect of fly ash and slag［J］．Cement ＆Concrete Research，1999，29（4）：487－495．

［11］ 金祖權(quán)，孫偉，張云升，等．混凝土在硫酸鹽、氯鹽溶液中的損傷過程［J］．硅酸鹽學(xué)報，2006，34（5）：630－635．Jin Zuquan，Sun Wei，Zhang Yunsheng，et al．Damage of concrete in sulfate and chloride solution［J］．Journal of the Chinese Ceramic Society，2006，34（5）：630－635．

［12］ 吳政，張承娟．單向荷載作用下巖石損傷模型及其力學(xué)特性研究［J］．巖石力學(xué)與工程學(xué)報，1996，15（1）：55－61．Wu Zheng，Zhang Chengjuan．Investigation of rock damage model，and its mechanical behaviour［J］．Chinese Journal of Rock Mechanics and Engineering，1996，15（1）：55－61．

［13］ 李兆霞．損傷力學(xué)及其應(yīng)用［M］．北京：科學(xué)出版社，2002．

［14］ 徐衛(wèi)亞，韋立德．巖石損傷統(tǒng)計本構(gòu)模型的研究［J］．巖石力學(xué)與工程學(xué)報，2002，21（6）：787－791．Xu Weiya，Wei Lide．Study on statistical damage constitutive model of rock［J］．Chinese Journal of Rock Mechanics and Engineering，2002，21（6）：787－791．

Dynamic constitutive model of concrete after salt corrosion

Nie Liangxue1，2，Xu Jinyu1，3，Liu Zhiqun1，Luo Xin4

（1．Department of Airfield and Building Engineering，Air Force Engineering University，Xian 710038，Shaanxi，China；2．State Key Laboratory of Disaster Prevention ＆Mitigation of Explosion ＆Impact，PLA University of Science and Technology，Nanjing210007，Jiangsu，China；3．College of Mechanics and Civil Architecture，Northwest Polytechnic University，Xian 710072，Shaanxi，China；4．Construction Engineering Planning ＆Design Institute，Logistic Support Department，Central Military Comm ission，Peoples Republic of China，Beijing100062，China）

To explore the dynamic mechanical behaviors of the concrete corroded by salt solution，we fabricated concrete specimens of ordinary Portland cement with fly ash of 15%mass fraction and，having them immersed for 60din NaCl and Na2SO4solutions with a 15%mass fraction，studied their dynamic mechanical properties using a 100mm split Hopkinson pressure bar apparatus．Then，based on the macro phase－only statistical damage theory and the Weibull distribution theory，we built a dynamic statistical damage constitutive model for the specimens mechanical behaviors．The results indicate that the dynamic compressive strength of the corroded specimens experiences a significant decrease，and the declines of the specimens corroded by NaCl are bigger than those corroded by Na2SO4．Based on the closely fitted curves of the model and the experiment，this established model can accurately estimate the dynamic mechanical behaviors of concrete under impact loading．

salt corrosion；concrete；dynamic compressive strength；damage theory；Weibull distribution；constitutive model

O347．3國標(biāo)學(xué)科代碼：13015

A

10．11883／1001－1455（2017）04－0712－07

（責(zé)任編輯 王 影）

2015－12－12；

2016－04－18

國家自然科學(xué)基金項目（51078350，51208507）；爆炸沖擊防災(zāi)減災(zāi)國家重點實驗室開放課題（DPMEIKF201406）；陜西省青年科技新星計劃項目（2013KJXX－81）

聶良學(xué)（1990－ ），男，博士研究生，nieliangxue＿kgd＠126．com。