(武漢大學(xué) 水資源與水電工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室，武漢 430072)

1 研究背景

鈣溶出性侵蝕是水泥基材料耐久性病害的一類，又稱溶蝕[1]，根據(jù)水泥基材料所處的環(huán)境的不同，可分為接觸溶蝕與滲透溶蝕[2]。接觸溶蝕是指水泥基材料在與水接觸過程中，水頭較低，水泥基材料內(nèi)部未形成滲流作用或滲流作用效應(yīng)可以忽略時(shí)，孔隙液中的鈣離子在濃度梯度作用下，不斷向外擴(kuò)散，打破了固液兩相鈣濃度間的平衡關(guān)系，使得固相鈣不斷溶出的過程[3]。

在自然條件下，接觸溶蝕過程十分緩慢，溶蝕鋒面由表及里移動(dòng)緩慢[4-5]。據(jù)報(bào)道，建筑物與水接觸100 a，溶蝕鋒面僅向前移動(dòng)5～10 mm[6]，這一溶蝕特性對(duì)水泥基材料鈣溶蝕性能的研究增加了一定的難度。為了能在短時(shí)間內(nèi)研究溶蝕對(duì)水泥基材料各方面性能的影響，加速溶蝕試驗(yàn)和溶蝕數(shù)值模擬成了眾多學(xué)者研究溶蝕的2種方法。在數(shù)值模擬方面，國內(nèi)外學(xué)者做了大量的研究。Reardon等[7]研究了水泥水化物與孔隙液中離子間的熱動(dòng)力平衡關(guān)系；Berner[8]建立了孔隙液中鈣離子濃度與水基材料中鈣含量之間的固液平衡關(guān)系；Gérard等[9]修正了固液平衡關(guān)系，建立起簡化的一維鈣溶蝕數(shù)學(xué)模型；Nguyen等[10]結(jié)合尺度均質(zhì)化分析建立了混凝土鈣離子浸出模型，引入孔隙曲折率來考慮骨料對(duì)溶蝕的影響；Wan等[11]通過試驗(yàn)得到了波特蘭水泥基材料在6 mol/L的硝酸銨溶液中的固液平衡關(guān)系曲線，并據(jù)此建立了加速水泥凈漿溶蝕的一維數(shù)學(xué)模型[6]； Yu等[12]耦合化學(xué)動(dòng)力學(xué)與化學(xué)熱力學(xué)提出了針對(duì)采用6 mol/L的硝酸銨溶液加速接觸溶蝕的一類新的模擬方法。

縱觀以上研究發(fā)現(xiàn)，水泥基材料接觸溶蝕數(shù)值模擬方法雖然較多，但大多以一維展開，求解方法以有限差分法為主，對(duì)于采用三維有限元方法模擬接觸溶蝕尚未見報(bào)道。對(duì)此，本文針對(duì)水泥基材料接觸溶蝕的三維有限元數(shù)值模擬展開研究，擬通過建立三維接觸溶蝕數(shù)學(xué)模型，采用FORTRAN語言設(shè)計(jì)有限元計(jì)算程序，通過數(shù)值計(jì)算研究各溶蝕特征量的時(shí)空分布規(guī)律。

2 三維鈣溶蝕數(shù)學(xué)模型

相關(guān)研究表明，水泥基材料在溶蝕過程中，鈣、硫、鉀、鋁等元素為主要溶出元素，而氫氧化鈣、水化硅酸鈣凝膠又是水泥的主要水化產(chǎn)物，因此，以鈣元素作為研究對(duì)象，可以較好地反映水泥基材料性能退化情況[13]。

參考現(xiàn)有的研究成果，本文根據(jù)鈣離子傳輸數(shù)學(xué)模型、固液平衡關(guān)系、孔隙率演化模型、有效擴(kuò)散系數(shù)模型建立三維鈣溶蝕數(shù)學(xué)模型，并作如下假設(shè)：①孔隙液中的鈣離子與固相鈣含量之間瞬間達(dá)到平衡狀態(tài)，即忽略固體鈣溶解的時(shí)間效應(yīng);②不考慮孔隙液中鈣離子與固體骨架產(chǎn)生吸附與解吸作用;③不考慮溫度對(duì)溶蝕的影響;④材料各向同性。

圖1 表征體元內(nèi)的質(zhì)量守恒Fig.1 Conservation of mass in representative elementary volume

2.1 鈣離子傳輸數(shù)學(xué)模型

2.1.1 鈣離子傳輸控制方程

在求解域內(nèi)任取一

個(gè)長、寬、高分別為Δx，Δy，Δz的六面體表征體元(REV)，如圖1所示。在域內(nèi)僅考慮鈣離子擴(kuò)散作用，根據(jù)表征體元內(nèi)鈣離子的質(zhì)量守恒關(guān)系，結(jié)合非克(Fick)第二定律可推導(dǎo)出水泥基材料溶蝕過程中鈣離子傳輸控制方程。

如圖1所示的表征體元，單位時(shí)間內(nèi)x方向進(jìn)入的鈣離子通量為

(1)

式中：Jx為x方向進(jìn)入的鈣離子通量；Ax為表征體元側(cè)面面積；φ為孔隙率；Dx為x方向擴(kuò)散系數(shù)；C為孔隙液中鈣離子濃度。

單位時(shí)間內(nèi)x方向從表征體元流出的鈣離子通量為

(Jx+ΔJx)Ax=

(2)

根據(jù)式(1)和式(2)可得，單位時(shí)間內(nèi)在x方向從表征體元進(jìn)入和流出的鈣離子通量的差值即為凈有流入量，如式(3)所示。

同理可得，單位時(shí)間內(nèi)在y方向、z方向的凈有流入量，分別如式(4)和式(5)所示。

式中：Jy為y方向進(jìn)入的鈣離子通量；Ay為表征體元側(cè)面面積；Dy為y方向擴(kuò)散系數(shù)；Jz為z方向進(jìn)入的鈣離子通量；Az為表征體元側(cè)面面積；Dz為z方向擴(kuò)散系數(shù)。

在溶蝕過程中，固相鈣不斷向孔隙液中溶解，假設(shè)固相鈣的摩爾濃度為CS，單位時(shí)間內(nèi)表征體元孔隙液內(nèi)增加的鈣離子物質(zhì)的量ΔQ1為

(6)

式中t為溶蝕時(shí)間。

水泥基材料為多孔介質(zhì)，假設(shè)孔隙內(nèi)充滿孔隙液，則表征體元中孔隙液的體積為φΔxΔyΔz，由此可得表征體元內(nèi)孔隙液中鈣離子含量為φΔxΔyΔzC，單位時(shí)間內(nèi)表征體元中鈣離子的物質(zhì)變化量ΔQ2為

(7)

根據(jù)質(zhì)量守恒定律，單位時(shí)間內(nèi)x，y，z方向的鈣離子凈有流入量與孔隙液中因固相鈣溶解的鈣離子增加量之和等于單位時(shí)間內(nèi)表征體元中鈣離子的物質(zhì)變化量，即

(8)

需要指出的是式中Dx，Dy，Dz為鈣離子在孔隙液中不同方向的擴(kuò)散系數(shù)，經(jīng)體積均勻化后得到有效擴(kuò)散系數(shù)，即

(9)

將式(9)代入式(8)整理可得鈣離子三維傳輸控制方程，即

(10)

2.1.2 定解條件

求解上述控制方程組，給定初始條件與邊界條件如下所述。

(1)初始條件：

C(x,y,z,t0)=C0。

(11)

式中C0為初始時(shí)刻孔隙液中鈣離子濃度。

(2)邊界條件：

①定通量邊界條件為

(12)

式中：q為邊界上已知鈣離子的擴(kuò)散通量；nx，ny，nz為邊界的單位外法線方向余弦分量。

②離子交換邊界條件[4, 14]。在以往的研究中，鈣溶蝕數(shù)學(xué)模型邊界條件多采用第一類邊界條件，即0濃度邊界條件，這一類邊界條件在模型求解過程中，邊界各點(diǎn)孔隙液濃度始終保持為0。由固液平衡關(guān)系可知，邊界各點(diǎn)的固相鈣濃度也始終保持為0，即邊界各點(diǎn)瞬間完成鈣溶出過程，這顯然與事實(shí)不符；實(shí)質(zhì)上，隨著靠近邊界側(cè)的求解域內(nèi)的濃度梯度變化，邊界各點(diǎn)的孔隙液濃度也在不斷發(fā)生變化。因此，采用物質(zhì)交換邊界條件更為符合實(shí)際情況，這類邊界條件認(rèn)為經(jīng)過邊界的鈣離子通量與邊界孔隙中鈣離子濃度C與環(huán)境水中鈣離子濃度C1之差成正比，其表達(dá)式如式(13)所示。采用這類邊界條件不僅可以反映邊界濃度變化，還可以通過邊界濃度求解鈣離子累積溶出量。本文采用該類邊界條件作為接觸邊界條件,即

式中：kT為邊界上孔隙液與環(huán)境水之間的傳質(zhì)系數(shù)；C1為環(huán)境水中鈣離子濃度。

2.2 固液平衡關(guān)系

控制方程式的固相鈣濃度CS與液相鈣濃度C并非獨(dú)立變量，相關(guān)研究表明，CS與C兩者具有一定的熱動(dòng)力平衡關(guān)系[7-8]。Buil’s等提出了一種便于計(jì)算的等溫平衡關(guān)系模型，Nakarai等進(jìn)一步修正了Buil’s模型得到固液平衡曲線[4]，Berner通過試驗(yàn)數(shù)據(jù)驗(yàn)證了該模型的合理性[8]。本文引用該固液平衡曲線得到C與CS間的關(guān)系，即

式中：x1為C-S-H(水化硅酸鈣凝膠)開始快速溶解時(shí)孔隙液中的鈣離子濃度；x2為Ca(OH)2完全溶解時(shí)孔隙液中的鈣離子濃度；Csatu為孔隙液中鈣離子飽和濃度；CCSH為水泥基材料中C-S-H的摩爾濃度；CCH為水泥基材料中Ca(OH)2的摩爾濃度。

根據(jù)固液平衡關(guān)系，鈣離子傳輸控制方程中固相鈣濃度可轉(zhuǎn)化為孔隙液中鈣離子濃度，進(jìn)一步得到含單一變量的控制方程，即

(15)

2.3 孔隙率演化模型

2.3.1 初始孔隙率計(jì)算

水泥基材料內(nèi)部孔隙分為毛細(xì)孔與凝膠孔，即

φp0=φc0+φg0。

(16)

式中：φp0為水泥漿體的初始孔隙率；φc0為水泥漿體初始毛細(xì)孔隙率；φg0為水泥漿體初始凝膠孔隙率。φc0與φg0根據(jù)POWERS模型計(jì)算可得，即

(17)

式中：W/C為水灰比；α為水化度，可由式(18)計(jì)算。

(18)

根據(jù)水泥漿體在水泥基材料的體積占比，可計(jì)算水泥基材料的初始孔隙率，即

φ0=fp/mfm/cφp0。

(19)

式中：φ0為水泥基材料初始孔隙率;fp/m為水泥漿在砂漿中的體積分?jǐn)?shù)；fm/c為砂漿在混凝土中的體積分?jǐn)?shù)。

2.3.2 孔隙率演化

隨著溶蝕過程的進(jìn)行，Ca(OH)2與C-S-H等水化產(chǎn)物不斷向孔隙液中溶解，使得孔隙體積不斷增大，因而材料孔隙率也不斷增大?？紤]到C-S-H溶解過程較為復(fù)雜，假設(shè)C-S-H以Ca(OH)2形式溶解[4, 6]，據(jù)此可建立孔隙率演化模型，即

(20)

式中：Δφ為孔隙率變化量；ΔCS為固相鈣濃度變化量；MCH為Ca(OH)2的摩爾質(zhì)量；ρCH為Ca(OH)2的密度；CS0為初始固相鈣濃度。

假設(shè)Ca(OH)2的溶解主要造成毛細(xì)孔的增加，而C-S-H的溶解在引起凝膠孔增大的同時(shí)，也會(huì)造成毛細(xì)孔的增大[5]。引入分配系數(shù)γ來考慮C-S-H的溶解對(duì)毛細(xì)孔隙率的影響，因此，可用式(21)、式(22)分別描述毛細(xì)孔隙率與凝膠孔隙率的變化。

(21)

(22)

2.4 擴(kuò)散系數(shù)演化模型

目前，有效擴(kuò)散系數(shù)模型均建立在不含粗骨料的水泥凈漿或砂漿的基礎(chǔ)上，對(duì)于含有粗骨料的水泥基材料無法適用。相關(guān)研究表明[15]，粗骨料對(duì)于擴(kuò)散系數(shù)的影響主要有2方面：一方面，粗骨料所在區(qū)域無擴(kuò)散現(xiàn)象發(fā)生；另一方面，粗骨料與水泥漿間存在界面過渡區(qū)，會(huì)造成孔隙增大，進(jìn)而加劇擴(kuò)散。為此，相關(guān)研究引入孔隙曲折度來考慮粗骨料對(duì)擴(kuò)散系數(shù)的綜合影響，建立水泥基材料的通用擴(kuò)散系數(shù)模型，如式(23)所示。

(23)

3 有限元計(jì)算方法

3.1 計(jì)算格式

采用三維8節(jié)點(diǎn)等參單元離散空間域，試探函數(shù)為式(24)所示。

(24)

式中：C代表空間域上任意一點(diǎn)的濃度；Nm為等參單元的形函數(shù)；Cm為等參單元的第m個(gè)節(jié)點(diǎn)濃度；Hm為等參單元的第m個(gè)節(jié)點(diǎn)總水頭；M為8。

根據(jù)Galerkin有限元計(jì)算方法與隱式差分方法分別對(duì)空間域、時(shí)間域進(jìn)行離散，得到有限元計(jì)算格式，如式(25)所示。

(25)

其中：

[DC]=

3.2 計(jì)算方法

3.2.1 協(xié)調(diào)質(zhì)量矩陣處理

對(duì)于時(shí)間域離散，筆者曾采用過顯式差分法、中心差分法、隱式差分法3種方法，計(jì)算過程中發(fā)現(xiàn)：

(1)采用中心差分格式與顯式差分格式時(shí)，只能選取較小的時(shí)間步長，才能保證數(shù)值計(jì)算結(jié)果的穩(wěn)定性，若時(shí)間步長過大，單元部分節(jié)點(diǎn)將出現(xiàn)無意義的負(fù)濃度。另外，這2種格式均不能有效地減少數(shù)值振蕩問題。

(2)采用隱式差分格式，對(duì)減小數(shù)值振蕩具有較為顯著的效果，而且不會(huì)因時(shí)間步長的變大節(jié)點(diǎn)出現(xiàn)無意義的負(fù)濃度，但時(shí)間步長較小時(shí)，也會(huì)出現(xiàn)較大的數(shù)值振蕩現(xiàn)象。

綜合考慮，為避免在小時(shí)間步長下出現(xiàn)數(shù)值振蕩，筆者在隱式差分格式離散時(shí)間域的基礎(chǔ)上，采用集中質(zhì)量矩陣代替協(xié)調(diào)質(zhì)量矩陣，即對(duì)式(25)中的[DC]進(jìn)行處理，具體處理方法詳見文獻(xiàn)[16]。

3.2.2 程序?qū)崿F(xiàn)

采用FORTRAN語言進(jìn)行程序設(shè)計(jì)，矩陣元素統(tǒng)一采用一維數(shù)組存儲(chǔ)，對(duì)矩陣應(yīng)用稀疏矩陣行壓縮存儲(chǔ)方案，線性方程組求解采用超松弛預(yù)處理共軛梯度法(SSOR-PCG)。

4 數(shù)值驗(yàn)證

本文分別對(duì)文獻(xiàn)[6]中一維加速溶蝕試驗(yàn)、文獻(xiàn)[17]中二維加速溶蝕試驗(yàn)開展三維數(shù)值模擬，通過數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果、有限差分法計(jì)算結(jié)果進(jìn)行對(duì)比，驗(yàn)證三維模型及有限元計(jì)算程序的正確性。

4.1 一維溶蝕模擬

引用文獻(xiàn)[6]中的加速水泥凈漿接觸溶蝕試驗(yàn)資料，進(jìn)行數(shù)值模擬。試件尺寸為40 mm×40 mm×20 mm，水灰比為0.35，在標(biāo)準(zhǔn)條件養(yǎng)護(hù)10個(gè)月后，采用6 mol/L的NH4NO3溶液進(jìn)行加速水泥凈漿試塊溶蝕試驗(yàn)。

假設(shè)初始時(shí)刻的試塊孔隙中充滿孔隙液，且孔隙液中鈣離子濃度為飽和濃度，如式(26)所示?；诒疚乃⒌臄?shù)學(xué)模型，根據(jù)文獻(xiàn)中的試驗(yàn)條件給出邊界條件，如式(27)所示，相應(yīng)的邊界條件示意圖如圖2所示。

C(x,y,z,t0)=Csatu。

(26)

(27)

圖2 邊界條件示意圖Fig.2 Schematic diagram of boundary conditions

通過水泥水化模型計(jì)算可得到試塊中C-S-H與Ca(OH)2的摩爾含量，加速溶蝕過程中的固液平衡關(guān)系曲線其他參數(shù)取值與文獻(xiàn)[6]相同，如表1所示，模型計(jì)算其他參數(shù)如表2所示。

表1 固液平衡曲線參數(shù)Table 1 Parameters of solid-liquid equilibrium curve mol/m3

表2 模型計(jì)算基本輸入?yún)?shù)Table 2 The basic input parameters of FEM program

注:kT取值參考文獻(xiàn)[18]。

圖3(a)為本文模擬的溶蝕深度與文獻(xiàn)[6]試驗(yàn)數(shù)據(jù)及有限差分法模擬的結(jié)果對(duì)比，從圖3(a)可知，本文模擬(FEM)結(jié)果與有限差分法模擬(FDM)結(jié)果具有較好的一致性，兩者與試驗(yàn)結(jié)果在試驗(yàn)前期吻合較好，后期具有一定的誤差；經(jīng)分析，這可能是由于經(jīng)過一段浸蝕時(shí)間后，溶蝕介質(zhì)NH4NO3溶液濃度有所降低，減緩了試件溶蝕進(jìn)程[6]。

圖3(b)—圖3(d)分別為本文計(jì)算結(jié)果與文獻(xiàn)采用有限差分法計(jì)算結(jié)果的對(duì)比情況。從孔隙液鈣離子濃度分布、固相鈣濃度分布及孔隙率分布3個(gè)方面對(duì)比可知，本文模型計(jì)算結(jié)果與文獻(xiàn)[6]模擬結(jié)果具有較好的一致性。

圖3 模擬溶蝕特征量與文獻(xiàn)[6]數(shù)據(jù)對(duì)比Fig.3 Comparison between modeling the spatial distributions of characteristic variables and the data from Ref.[6]

4.2 二維溶蝕模擬

引用文獻(xiàn)[17]中的接觸溶蝕試驗(yàn)資料，對(duì)試驗(yàn)中水泥凈漿試塊進(jìn)行加速溶蝕數(shù)值模擬。試驗(yàn)中水泥凈漿試塊尺寸為10 mm×10 mm×2 mm，水灰比為0.5，采用6 mol/L的NH4NO3溶液進(jìn)行加速水泥凈漿溶蝕試驗(yàn)。由于試件上下側(cè)與骨料接觸，本文假設(shè)試件這兩側(cè)不與NH4NO3溶液接觸，即認(rèn)為該試驗(yàn)為二維溶蝕試驗(yàn)。邊界條件示意圖如圖4所示。

選取中間截面進(jìn)行溶蝕深度分析，溶蝕深度隨時(shí)間變化關(guān)系如圖5所示，與文獻(xiàn)[17]試驗(yàn)數(shù)據(jù)對(duì)比發(fā)現(xiàn)，本文模擬計(jì)算結(jié)果與試驗(yàn)數(shù)據(jù)吻合良好。

圖4 邊界條件示意圖Fig.4 Schematic diagram of boundary conditions

圖5 模擬溶蝕深度與文獻(xiàn)[17]數(shù)據(jù)對(duì)比Fig.5 Comparison between modeling leaching depth and the experimental results from Ref.[17]

圖6為溶蝕24 h后各溶蝕特征量分布云圖，可以看出，固體鈣濃度由內(nèi)及外逐漸降低，溶蝕程度逐漸增強(qiáng)，與文獻(xiàn)[17]試驗(yàn)結(jié)果對(duì)比發(fā)現(xiàn)，固體鈣濃度分布與試驗(yàn)所到溶蝕分布具有較好的一致性?？紫兑褐锈}離子分布與固體鈣濃度分布規(guī)律基本一致。需要指出的是，孔隙率分布規(guī)律與固相鈣濃度分布規(guī)律正好相反，這是由于接觸溶蝕是一個(gè)由表及里的發(fā)展過程，與內(nèi)部相比，外部固相鈣溶出量更大，孔隙體積增大程度更大。

圖6 模擬溶蝕特征量分布與文獻(xiàn)[17]試驗(yàn)結(jié)果對(duì)比Fig.6 Comparison between modeling the spatial distributions of characteristic variables after 24 hours and the experimental results from Ref.[17]

選取接觸邊界節(jié)點(diǎn)進(jìn)行分析，其溶蝕特征量隨時(shí)間變化關(guān)系如圖7所示。從圖7(a)可以看出，邊界節(jié)點(diǎn)隨著時(shí)間增加孔隙液鈣離子濃度逐漸減小并最終趨于穩(wěn)定；從圖7(b)可以看出，固相鈣含量變化呈現(xiàn)出相同的規(guī)律。這充分表明離子交換邊界條件能夠有效地避免采用“0濃度”邊界條件造成的問題，即接觸邊界節(jié)點(diǎn)孔隙液中鈣離子濃度始終為0，接觸邊界節(jié)點(diǎn)固相鈣溶出過程瞬間完成。

圖7 邊界節(jié)點(diǎn)孔隙液鈣離子濃度、固相鈣濃度 隨時(shí)間變化關(guān)系曲線Fig.7 Time-histories of calcium ion concentration in pore fluid and solid phase

通過數(shù)值模擬結(jié)果與試驗(yàn)數(shù)據(jù)、有限差分法模擬結(jié)果對(duì)比可知，各項(xiàng)溶蝕特征量數(shù)據(jù)吻合良好，本文模型及程序的正確性得到驗(yàn)證；采用離子交換邊界條件能夠有效避免“0濃度”邊界條件造成的問題。

圖8 試件5有限元模型Fig.8 FEM model of specimen 5

5 算 例

5.1 模型簡介

目前溶蝕試驗(yàn)研究大多以一維、二維展開，三維較少涉及，因此，本文給出一個(gè)三維溶蝕的數(shù)值算例，研究三維溶蝕下的各溶蝕特征量的空間分布規(guī)律?？紤]一個(gè)圓柱體水泥凈漿試件，尺寸為Φ20 mm×20 mm，將其置于6 mol/L的NH4NO3溶液中，使試件與溶液充分接觸；水泥化學(xué)組分及其他模型計(jì)算基本輸入?yún)?shù)取值見表1、表2；試件的有限元模型如圖8所示。

圖9 各截面溶蝕特征量模擬計(jì)算結(jié)果Fig.9 Modeling results of characteristic variables at each section

5.2 結(jié)果分析

為了研究溶蝕的空間分布規(guī)律，選取y=0 mm和z=10 mm 2個(gè)截面展開分析。經(jīng)過40 d溶蝕后，各截面的各溶蝕特征量模擬計(jì)算結(jié)果云圖如圖9所示。從圖9可以看出，無論是徑向還是軸向，各溶蝕特征量的空間分布規(guī)律與上文算例呈現(xiàn)出一致的規(guī)律。

為了研究各溶蝕特征量隨時(shí)間的變化規(guī)律，在試件徑向取4個(gè)特征點(diǎn)，分別距離表層為1，2，3，5 mm，編號(hào)依次為①,②,③,④；在試件軸向，以試件頂面圓心為基點(diǎn)，從上至下取4個(gè)特征點(diǎn)，分別距離表層為1，2，3，5 mm，編號(hào)依次為⑤,⑥,⑦,⑧。根據(jù)計(jì)算可得，各溶蝕特征量與時(shí)間的變化關(guān)系，如圖10所示?？梢钥闯?，隨著時(shí)間的變化，特征點(diǎn)的各溶蝕特征量變化速率均呈現(xiàn)出由快變慢的總體規(guī)律，說明各處的溶蝕速率由快變慢；從圖10(b)可以看出，固相鈣濃度在7 000 mol/m3左右時(shí)，溶蝕速率開始變慢，這是由于此時(shí)Ca(OH)2溶解完成，C-S-H開始溶解，而C-S-H前期的溶解速率遠(yuǎn)遠(yuǎn)小于Ca(OH)2的溶解速率。無論是徑向還是軸向，各個(gè)特征點(diǎn)開始溶蝕的時(shí)刻均不相同，呈現(xiàn)出外部早于內(nèi)部的總體規(guī)律，體現(xiàn)出接觸溶蝕是一個(gè)由表及里的過程。需要指出的是，在徑向與軸向上，距表層距離相同的特征點(diǎn)的溶蝕開始時(shí)間并不相同，徑向明顯早于軸向，據(jù)分析，這是由于所選取的徑向特征點(diǎn)所在處的濃度梯度大于軸向特征點(diǎn)，從而造成孔隙液中鈣離子徑向擴(kuò)散速率大于軸向擴(kuò)散速率。

圖10 各截面特征點(diǎn)孔隙液中鈣離子濃度、固相鈣濃度、 孔隙率隨時(shí)間變化的變化關(guān)系曲線Fig.10 Time-histories of calcium ion concentration in pore fluid and solid phase and porosity at feature points

6 結(jié) 論

(1)在隱式差分格式離散時(shí)間域的基礎(chǔ)上，用集中質(zhì)量矩陣代替協(xié)調(diào)質(zhì)量矩陣可以消除對(duì)時(shí)間步長選擇的限制，能夠有效避免數(shù)值振蕩現(xiàn)象。

(2)稀疏矩陣行壓縮存儲(chǔ)方案雖然具有內(nèi)存消耗較小的特點(diǎn)，但在形成總體剛度矩陣時(shí)卻十分耗時(shí)，對(duì)大型有限元計(jì)算較為不利。

(3)本文建立了三維鈣溶蝕模型，通過對(duì)加速接觸溶蝕過程數(shù)值模擬，經(jīng)過結(jié)果對(duì)比驗(yàn)證，表明該模型能夠用于模擬三維狀態(tài)下發(fā)生的各維度的接觸溶蝕，具有一定的通用性。

(4)采用離子交換邊界條件，得到邊界濃度變化規(guī)律，解決了采用“0濃度邊界條件”造成接觸邊界各點(diǎn)瞬間完成溶出過程的問題。

(5)通過開展三維溶蝕數(shù)值試驗(yàn)，結(jié)果表明：在三維狀態(tài)下，接觸溶蝕發(fā)展呈現(xiàn)出由表及里，各處溶蝕速率由快到慢，徑向溶蝕速率大于軸向溶蝕速率的總體特征。