田 豐

(新疆水利水電勘測(cè)設(shè)計(jì)研究院檢測(cè)試驗(yàn)研究中心，新疆 昌吉 831100 )

在水利工程中混凝土材料很難不應(yīng)用到，所建水利設(shè)施壽命在一定程度上取決于所用混凝土材料性能，水利設(shè)施受到凍融循環(huán)、雨水侵蝕等長(zhǎng)此以往的破壞，對(duì)內(nèi)部混凝土性能均是較大的傷害，特別是在我國(guó)西北、東北等氣候寒冷地區(qū)，研究混凝土凍融循環(huán)性能對(duì)開展水利工程選材具有重要幫助[1-3]。國(guó)外已有諸多學(xué)者[4-6]基于室內(nèi)試驗(yàn)與理論推導(dǎo)等手段，開展了混凝土凍融循環(huán)研究。國(guó)內(nèi)一些水利工程師[7-9]基于數(shù)值模擬手段，引入巖石力學(xué)損傷分析理論，探討了混凝土凍融循環(huán)損傷演化特征。還有一些學(xué)者[10-12]利用先進(jìn)的試驗(yàn)系統(tǒng)設(shè)計(jì)凍融循環(huán)實(shí)驗(yàn)方案，開展破壞力學(xué)、細(xì)觀等理論與試驗(yàn)相結(jié)合的研究。利用PFC 3D顆粒流軟件，不同于其他數(shù)值軟件，混凝土這種材料非常契合PFC 3D的研究對(duì)象特征[13-14]，開展多次凍融循環(huán)次數(shù)的損傷破壞實(shí)驗(yàn)研究，為認(rèn)識(shí)水利工程中混凝土材料凍融性能提供一定參考。

1 凍融損傷混凝土模型建立

寒區(qū)水利工程中混凝土常常會(huì)受到凍融侵蝕，導(dǎo)致混凝土從表到里逐漸失去保護(hù)層，強(qiáng)度或其他抗?jié)B特性均會(huì)受到一定損害，為此，本文在利用PFC 3D研究?jī)鋈趽p傷混凝土?xí)r，以混凝土為層狀結(jié)構(gòu)考慮，分為損傷層與無損傷層，如圖1所示構(gòu)造。并依據(jù)前人文獻(xiàn)知，凍融損傷層厚度與凍融循環(huán)次數(shù)關(guān)系如公式(1)：

d=1.742×10-5N3-0.004N2+0.322N+0.179

(1)

式中：d為凍融損傷層厚度，mm；N為凍融循環(huán)次數(shù)。

圖1 凍融損傷混凝土幾何示意圖

利用PFC 3D建立損傷層模型與未損傷層模型，如圖2所示，不論是未損傷層亦或是損傷層，混凝土其所用材料屬性參數(shù)均是一致，例如水泥砂漿膠結(jié)面材料彈性模量均為40 GPa，但是凍融損傷層由于凍融作用，膠結(jié)面上初始彈性模量為40 GPa，在其后凍融循環(huán)次數(shù)中，會(huì)逐漸降低，例如凍融循環(huán)50次，彈性模量為30 GPa。

圖2 損傷層模型與未損傷層模型

圖3(a)為凍融循環(huán)下混凝土應(yīng)力應(yīng)變。從圖中可看出，峰值強(qiáng)度與凍融循環(huán)次數(shù)成反比，循環(huán)次數(shù)為0時(shí)，峰值強(qiáng)度為35 MPa，而凍融循環(huán)次數(shù)達(dá)到150次時(shí)，其峰值應(yīng)力相比前者降低了20%，為28 MPa。另一方面，峰值應(yīng)力點(diǎn)處應(yīng)變值隨凍融循環(huán)次數(shù)增多而逐漸增大，在循環(huán)25次時(shí)，峰值應(yīng)力點(diǎn)應(yīng)變?yōu)?.13%，而經(jīng)過150次循環(huán)后，峰值點(diǎn)應(yīng)變值增大了23.1%，達(dá)0.16%。分析出現(xiàn)該現(xiàn)象是由于凍融循環(huán)次數(shù)增多，試樣內(nèi)部礦物顆粒遭到多次破壞再結(jié)晶，顆粒間黏結(jié)性能大大降低，承受單軸抗壓強(qiáng)度故逐漸減??；但凍融循環(huán)次數(shù)增大，一定程度會(huì)促使礦物顆粒塑性能力增強(qiáng)，混凝土軟化特性加強(qiáng)，反映在變形上則是峰值點(diǎn)應(yīng)變值逐漸增大。

圖3(b)～(d)為凍融循環(huán)次數(shù)與泊松比關(guān)系，從圖中亦可發(fā)現(xiàn)，彈性階段泊松比與凍融循環(huán)次數(shù)亦成反比，循環(huán)次數(shù)25次泊松比值為0.166，而150次循環(huán)次數(shù)下泊松比值相比降低了3.6%，為0.160，但整個(gè)破壞過程中泊松比變化趨勢(shì)不受凍融循環(huán)次數(shù)影響。

圖3 試樣破壞力學(xué)特征曲線

2 單軸細(xì)觀破壞模擬研究

2.1 凍融損傷單軸破壞分析

為分析方便，本文以混凝土試樣單軸破壞過程中的峰值強(qiáng)度點(diǎn)A、應(yīng)變達(dá)0.3%點(diǎn)B、峰后達(dá)1/2峰值應(yīng)力點(diǎn)C為對(duì)象，分析凍融循環(huán)下試樣的變形以及裂紋演化過程。

圖4給出了不同凍融循環(huán)次數(shù)下特征點(diǎn)處位移值變化曲線，從圖中可知，在各應(yīng)力或應(yīng)變特征點(diǎn)處，對(duì)應(yīng)的位移均隨凍融循環(huán)次數(shù)增多而增大，在應(yīng)變0.3%B點(diǎn)凍融循環(huán)次數(shù)0次下位移為0.724 mm，而在循環(huán)次數(shù)150次時(shí)，位移值增大了3.5%，為0.749 mm，且在凍融循環(huán)次數(shù)75次、100次時(shí)，兩者之間位移值基本一致，表明應(yīng)變0.3%B點(diǎn)處位移隨凍融循環(huán)次數(shù)增大，變化量較低，受凍融循環(huán)次數(shù)影響較弱。不同于B點(diǎn)，峰后1/2峰值應(yīng)力點(diǎn)位移值隨凍融循環(huán)次數(shù)變化顯著，凍融循環(huán)0次時(shí)位移值為0.918 mm，而次數(shù)為25次時(shí)相比增長(zhǎng)了7.8%，為0.990 mm，凍融循環(huán)次數(shù)愈多，位移增長(zhǎng)幅度愈明顯。從位移分布云圖亦可看出(圖5)，點(diǎn)B位移值端部最大位移在凍融循環(huán)次數(shù)為75次，與150次時(shí)并無太顯著差別，且B點(diǎn)時(shí)試樣整體位移分布亦近乎一致。

圖4 不同凍融循環(huán)次數(shù)下特征點(diǎn)處位移值變化曲線

圖5 特征點(diǎn)應(yīng)力分布云圖

圖6為特征點(diǎn)處應(yīng)力值與凍融循環(huán)次數(shù)關(guān)系曲線，圖7為特征點(diǎn)處應(yīng)力分布圖。從圖中可看出，A點(diǎn)與C點(diǎn)應(yīng)力值基本均隨凍融循環(huán)次數(shù)增多而逐漸降低，循環(huán)次數(shù)0次時(shí)C點(diǎn)荷載為16.6 kN，而在循環(huán)次數(shù)為125次、150次時(shí)，分別降低了3.6%、7.2%；但特征點(diǎn)B應(yīng)力值變化并不呈現(xiàn)趨勢(shì)一致性，而是隨著循環(huán)次數(shù)增多出現(xiàn)波動(dòng)性變化。從分布形態(tài)亦可看出，應(yīng)力分布總是趨向于端部處，導(dǎo)致凍融循環(huán)次數(shù)增多，一定程度會(huì)削弱端部晶體礦物顆粒間的整體承載能力，進(jìn)而呈現(xiàn)特征點(diǎn)A處應(yīng)力逐漸下降。

圖6 特征點(diǎn)處應(yīng)力值與凍融循環(huán)次數(shù)關(guān)系曲線

圖7 不同凍融循環(huán)次數(shù)下特征點(diǎn)處應(yīng)力分布

2.2 細(xì)觀特征

圖8為不同凍融循環(huán)次數(shù)下各特征點(diǎn)處產(chǎn)生的裂紋數(shù)。從圖中可看出，在點(diǎn)A與點(diǎn)C處裂紋總數(shù)在各循環(huán)次數(shù)下幾乎并未有顯著一致性變化，而是隨凍融循環(huán)次數(shù)增多，呈波動(dòng)性變化，點(diǎn)A循環(huán)次數(shù)0次時(shí)總裂紋數(shù)量為7587，而在循環(huán)次數(shù)25次與50次時(shí)，總裂紋數(shù)量幾乎一致，另在點(diǎn)C處亦是如此。分析出現(xiàn)這種現(xiàn)象與A點(diǎn)、C點(diǎn)處選取有關(guān)，A點(diǎn)與C點(diǎn)均是由各凍融循環(huán)次數(shù)下的應(yīng)力值控制確定的點(diǎn)，而不是應(yīng)變控制，當(dāng)采用應(yīng)力控制時(shí)，試樣裂紋的產(chǎn)生與應(yīng)力增長(zhǎng)有關(guān)，應(yīng)力點(diǎn)不變，進(jìn)而該點(diǎn)裂紋數(shù)量亦不會(huì)發(fā)生太大變化，此從室內(nèi)巖石力學(xué)試驗(yàn)亦可得知，當(dāng)采用力控加載時(shí)，巖石應(yīng)力增長(zhǎng)穩(wěn)定，但破壞劇烈。當(dāng)由應(yīng)變控制時(shí)，即B點(diǎn)，該點(diǎn)處裂紋數(shù)量與凍融循環(huán)次數(shù)成反比變化特征，混凝土試樣區(qū)域軟化，抑制了裂紋的產(chǎn)生及貫通，故而不論是張拉裂紋亦或是剪切裂紋數(shù)量，均呈下降態(tài)勢(shì)。從混凝土試樣破壞整體裂紋分布亦可看出，不論是特征點(diǎn)A、B亦或是點(diǎn)C，在各凍融循環(huán)次數(shù)下均是剪切裂紋數(shù)量高于張拉裂紋，由此說明，在經(jīng)凍融循環(huán)后混凝土單軸破壞還是以剪切破壞為主，即破壞內(nèi)因由抗剪能力差引起。

圖8 不同凍融循環(huán)次數(shù)下特征點(diǎn)處產(chǎn)生裂紋數(shù)

圖9 混凝土內(nèi)部破裂形態(tài)

針對(duì)混凝土此類材料，其裂紋的貫通主要在峰值應(yīng)力后出現(xiàn)，故本文給出點(diǎn)B、點(diǎn)C處混凝土內(nèi)部破裂形態(tài)，如圖9所示。從圖中可看出，點(diǎn)B處凍融循環(huán)次數(shù)0次、75次試樣內(nèi)部出現(xiàn)顯著貫通性裂紋，而在凍融循環(huán)150次時(shí)，試樣僅局部出現(xiàn)裂紋，整體處于穩(wěn)定狀態(tài)。分析此是由于點(diǎn)B為由應(yīng)變控制確定，當(dāng)凍融循環(huán)次數(shù)較少時(shí)，峰值點(diǎn)對(duì)應(yīng)的應(yīng)變值亦較小，而為達(dá)到點(diǎn)B的應(yīng)變值，試樣內(nèi)部需發(fā)生更大的變形，故造成試樣裂紋貫通上下端部，而凍融循環(huán)次數(shù)多時(shí)，峰值點(diǎn)對(duì)應(yīng)的應(yīng)變較大，而點(diǎn)B確定的應(yīng)變值靠近峰值點(diǎn)，故而出現(xiàn)的裂紋還未及貫通。點(diǎn)C處不論是凍融循環(huán)次數(shù)0亦或是150，內(nèi)部裂紋發(fā)育充分，均產(chǎn)生了貫通性破壞裂紋，不同于點(diǎn)B處，點(diǎn)C處為峰后階段1/2的峰值應(yīng)力，即裂紋均是在峰值應(yīng)力后期產(chǎn)生，內(nèi)部本身就已充滿較多裂隙，峰后逐漸聯(lián)通，并在點(diǎn)C處貫通成大裂紋，顆粒破碎性顯著，進(jìn)而導(dǎo)致試樣失穩(wěn)破壞。

3 結(jié) 論

針對(duì)某水利工程中所用混凝土材料，引入PFC 3D顆粒離散元分析方法，建立凍融損傷混凝土模型，研究了不同凍融循環(huán)次數(shù)下單軸破壞特征，得到了以下結(jié)論：

(1)研究了峰值強(qiáng)度、彈性階段泊松比與凍融循環(huán)次數(shù)成反比，但峰值點(diǎn)處應(yīng)變值成正比，凍融循環(huán)次數(shù)150次相比0次峰值應(yīng)力降低了20%，峰值點(diǎn)處應(yīng)變值增大了23.1%，達(dá)0.16%。

(2)獲得了峰值應(yīng)力點(diǎn)A、應(yīng)變0.3%點(diǎn)B、峰后1/2峰值應(yīng)力點(diǎn)C處應(yīng)力與變形特征，點(diǎn)B處位移受凍融循環(huán)次數(shù)影響較弱，點(diǎn)C處位移隨循環(huán)次數(shù)增長(zhǎng)較明顯；點(diǎn)A與點(diǎn)C應(yīng)力隨循環(huán)次數(shù)逐漸降低，點(diǎn)B應(yīng)力變化并無趨勢(shì)一致性。

(3)分析了典型特征點(diǎn)裂紋變化，點(diǎn)A與點(diǎn)C產(chǎn)生裂紋總數(shù)隨循環(huán)次數(shù)無一致性規(guī)律，點(diǎn)B處裂紋總數(shù)隨凍融循環(huán)次數(shù)逐漸減少，三個(gè)特征點(diǎn)均是剪切裂紋數(shù)量高于張拉裂紋；點(diǎn)B處循環(huán)0次、75次出現(xiàn)貫通性裂紋，循環(huán)150次時(shí)，僅局部出現(xiàn)裂紋，點(diǎn)C為峰后階段，在各循環(huán)次數(shù)下均出現(xiàn)有貫通性破壞裂紋。