基于數(shù)字圖像內(nèi)核的再生混凝土氯離子擴散分析

謝財奇， 吳宇清

(同濟大學，上海200092)

基于數(shù)字圖像內(nèi)核的再生混凝土氯離子擴散分析

謝財奇， 吳宇清

(同濟大學，上海200092)

沿海及除冰鹽環(huán)境下，再生混凝土 (Recycled Aggregate Concrete,RAC) 的耐久性易受到氯離子擴散的影響。再生混凝土的粗骨料表面往往附著一層老砂漿，結(jié)構(gòu)相對復雜，具有顯著的微觀不確定性。為了研究微觀不確定性對再生混凝土氯離子擴散的影響，本文采用掃描電子顯微鏡對再生混凝土粗骨料進行觀測，得到了一系列再生混凝土細觀結(jié)構(gòu)的數(shù)字圖像。根據(jù)再生混凝土骨料的特點，提出了結(jié)合隨機均質(zhì)化和多尺度的數(shù)字圖像驅(qū)動理論，建立了再生混凝土五層模型，推導了再生混凝土氯離子擴散的理論公式，并對上述圖像進行了數(shù)字圖像驅(qū)動線下部分的分析，求解數(shù)字圖像內(nèi)核方程后得到的特征值和特征向量十分理想，耗費較小的計算量即可得到誤差較小的結(jié)果，表明該研究方法切實可行。

再生混凝土;多尺度;均質(zhì)化;氯離子擴散

鑒于耐久性對混凝土結(jié)構(gòu)應用上的重要影響，越來越多的研究者開始關(guān)注混凝土結(jié)構(gòu)的耐久性問題。氯離子擴散性對于評價沿海及除冰鹽環(huán)境下混凝土結(jié)構(gòu)的耐久性能十分重要[1]。相對于普通混凝土，再生混凝土在環(huán)境保護和資源循環(huán)利用方面有很大優(yōu)勢[2]。再生混凝土的粗骨料表面往往附著了一層老砂漿，使其結(jié)構(gòu)具有顯著的微觀不確定性[3]，氯離子在其中的擴散行為更為復雜。許多學者對再生混凝土氯離子擴散特性進行了相關(guān)的試驗研究[4-6]，由于再生混凝土結(jié)構(gòu)具有較大的不確定性，試驗結(jié)果往往較為離散，需要采用理論手段對其進行進一步的研究。

為了提高計算效率，本文采用結(jié)合多尺度及均質(zhì)化的數(shù)字圖像驅(qū)動理論研究微觀不確定性對再生混凝土氯離子擴散的影響。多尺度理論早期用于研究一維場中橢圓方程的多尺度系數(shù)[7-8]，而后經(jīng)Hou等推廣到一般性的多相材料[9]。均質(zhì)化理論對于研究多相材料的結(jié)構(gòu)模型十分有效。Kamiński提出了研究具有隨機界面缺陷復合材料的隨機多尺度方法[10]。Sakata提出了基于攝動理論的均質(zhì)化方法來研究非均勻材料的彈性性質(zhì)的微觀不確定性[11]。研究結(jié)果表明，使用多尺度理論和均質(zhì)化理論研究多相非均勻介質(zhì)問題均十分有效，因此有必要將其納入到再生混凝土氯離子擴散問題的研究中。

數(shù)字圖像驅(qū)動理論分線下和線上兩部分，其中線下部分需要基于數(shù)字圖像求解協(xié)方差矩陣。目前國內(nèi)已有學者采用數(shù)字圖像對土體剪切帶[12]和再生混凝土疲勞破壞[13]進行研究，但尚無學者采用數(shù)字圖像對再生混凝土中氯離子擴散行為進行研究。

本文采用掃描電子顯微鏡對再生混凝土粗骨料進行觀測，得到其細觀結(jié)構(gòu)的數(shù)字圖像，基于數(shù)字圖像驅(qū)動理論，推導出了分析再生混凝土氯離子擴散的理論公式，并用其對上述圖像進行了數(shù)字圖像驅(qū)動線下部分的分析。

1 數(shù)字圖像驅(qū)動理論

1.1 五層RAC模型

再生混凝土相對普通混凝土而言更為復雜，其組成成分更多，為了便于分析，將再生混凝土簡化為五層模型，即天然粗骨料、老界面過渡區(qū)、老砂漿、新界面過渡區(qū)和新砂漿。其中，骨料部分結(jié)構(gòu)密實，視為不可滲透。砂漿部分往往含有一些細骨料、空隙、裂縫等，對砂漿部分取平均值，采用統(tǒng)一的氯離子擴散參數(shù)。界面過渡區(qū)在參數(shù)處理上與砂漿類似，同時考慮了微觀不確定性對其的影響。圖 1 即為模型再生混凝土的結(jié)構(gòu)示意圖，此處粗骨料取為圓形，即為圖中的黑色部分，圖中白色部分為界面過渡區(qū)，灰色部分為砂漿。其中界面過渡區(qū)假定為環(huán)繞內(nèi)部材料的薄殼形結(jié)構(gòu)。

圖1 再生混凝土圓形結(jié)構(gòu)模型示意圖

基于Zheng等的研究成果[14]，砂漿部分的氯離子擴散系數(shù)為：

(1)

式中：Dp為空隙的氯離子擴散系數(shù)；fp為空隙率。

1.2 氯離子擴散的多尺度分析

現(xiàn)有的有限元分析方法在研究再生混凝土這樣的多相結(jié)構(gòu)時，過程十分復雜且耗時較多。而工程師往往更關(guān)注混凝土的宏觀性能，故采用多尺度的分析方法十分必要。首先考慮一個代表體元Y?R2，即x3=0作為Ω。采用y=(y1,y2)作為再生混凝土微觀尺度下的坐標系，x=(x1,x2)為宏觀尺度下的坐標系。任意定義在Y上的周期性的函數(shù)F可以表示為

(2)

相應地，氯離子擴散系數(shù)張量可以表示為

Dε(x)=D(y)。

(3)

氯離子擴散方程為

·(Dεμ)+Q=0，x∈Ω。

(4)

式中：Q為邊界的氯離子濃度，為向量；μ為對應的方程解向量。

1.3 氯離子擴散的均質(zhì)化過程

上一節(jié)將微觀的不確定性引入到了宏觀行為的分析中，然而在各相交界處相關(guān)函數(shù)存在一定的跳躍性，故采用均質(zhì)化理論對其進行處理。

在Hilbert空間

下引入一個二維函數(shù)aε(μ,ξ)，其滿足

(5)

引入一個一維函數(shù)L(ξ) ，其滿足

(6)

氯離子擴散問題的變分等價形式就是找到με∈Λ使式(7)對于任意取值范圍內(nèi)的氯離子濃度ξ成立，

aε(με,ξ)=L(ξ)。

(7)

由文獻[15]進一步推導可得有效擴散系數(shù)為

(8)

該均質(zhì)化問題的控制方程有如下形式：

(9)

1.4 數(shù)字圖像驅(qū)動隨機方法

假定材料的擴散系數(shù)為非高斯分布隨機場，從而擴散系數(shù)矩陣可表示為

(10)

數(shù)字圖像驅(qū)動隨機方法包括線下和線上兩個部分。這里介紹線下部分。首先基于圖像數(shù)據(jù)得到協(xié)方差函數(shù)R(y, y′)作為數(shù)字圖像內(nèi)核，然后建立矩陣求解隨機基函數(shù)，根據(jù)Karhunen-Loeve(KL)展開用隨機基函數(shù)表示控制方程，公式推導如下：

H(y,ω)可以表示為

(11)

即將H(y,ω)中的隨機性集中到Ai(ω)上，其中Hi(y)與H(y,ω)相對應，為一參數(shù)。

為求解Ai(ω)，先考慮一個一維問題

(12)

(13)

(14)

式中：η(y)為特征向量；λ為特征值；di為計算中引入的參數(shù)。

將式(14)代入式(13)可得到誤差公式：

(15)

考慮到正交性，有

(16)

寫成矩陣形式

RD=ΓBD，

(17)

(18)

(19)

(20)

Γij=δijλi。

(21)

采用KL展開的截斷面形式，誤差可以表示為

(22)

(23)

2 算例分析

2.1 掃描電鏡觀測

采用S-3400SEM*掃描電鏡，放大倍數(shù)為10 000，對再生混凝土粗骨料預處理后進行觀測，得再生混凝土粗骨料的SEM照片如圖2所示。

圖2 再生粗骨料SEM照片

2.2 協(xié)方差函數(shù)特征方程求解

圖3 M=3時的特征值

圖4 M=3時的特征向量

圖5 M=4時的特征值

圖6 M=4時的特征向量

圖7 M=5時的特征值

圖8 M=5時的特征向量

由圖3—8可知，特征值衰減十分迅速，特征向量呈三角函數(shù)分布，結(jié)果符合預期。

由式(22)可計算相應誤差，限于篇幅，僅列出以下結(jié)果，如圖9—11所示。

圖9 M=3、4、5、6時的誤差

圖10 M=10時的誤差

圖11 M=15時的誤差

從圖9中可以看出，M值較小時，計算誤差總體而言處于較低水平，最大值出現(xiàn)在M=3時，為0.027?？紤]到M=3時，采用KL展開法選取項數(shù)少，誤差確實會偏大，隨著項數(shù)增加，誤差逐漸減小。由圖10—11可知，當M超過10時，誤差顯著增大。綜上所述，M取值為6時能取得十分理想的結(jié)果，同時計算量相對較小。

3 結(jié) 語

為了研究微觀不確定性對再生混凝土氯離子擴散行為的影響，本文采用掃描電子顯微鏡觀測了再生粗骨料的細觀結(jié)構(gòu)，基于數(shù)字圖像驅(qū)動理論推導出了再生混凝土氯離子擴散行為的理論分析方程，并對觀測結(jié)果進行了數(shù)字圖像驅(qū)動線下部分的分析。從數(shù)值計算的結(jié)果可知，該方法是切實可行的，求解出的數(shù)字圖像內(nèi)核特征方程的特征值和特征向量十分理想，耗費較小的計算量即可得到誤差較小的結(jié)果。限于篇幅，本文僅對數(shù)字圖像驅(qū)動的線下部分進行了研究，并未對再生混凝土氯離子擴散行為的過程進行分析，該部分將納入后續(xù)研究中，本文的成果將為再生混凝土氯離子擴散行為的分析提供一定的參考。

[1]Zheng Jianjun,Zhou Xinzhu,Wu Yufei,et al.A numerical method for the chloride diffusivity in concrete with aggregate shape effect[J].Construction and Building Materials,2012,31:151-156.

[2]Ying J,Xiao J,Tam V W Y.On the variability of chloride diffusion in modelled recycled aggregate concrete[J].Construction and Building Materials,2013,41:732-741.

[3]Xiao J,Ying J.Meso-level numerical simulation on two-dimensional chloride diffusivity in modeled recycled aggregate concrete[J].Journal of Tongji University(Natural Science),2012,40(7):1051-1057.

[4]Wang H L,Wang J J,Sun X Y,et al.Chloride diffusivity characteristics of the new interface transition zone in recycled aggregate concrete[J].Advanced Materials Research,2011,261:104-110.

[5]Xiao J,Ying J,Shen L.FEM simulation of chloride diffusivity in modeled recycled aggregate concrete[J].Construction and Building Materials,2012,29:12-23.

[6]Ying Jingwei,Xiao Jianzhuang,Shen Luming,et al.Five-phase composite sphere model for chloride diffusivity prediction of recycled aggregate concrete[J].Magazine of Concrete Research,2013,65(9):573-588.

[7]Babuska I,Osbom E.Generalized finite element methods:their performance and their relation to mixed methods [J].SIAM Journal on Numerical Analysis,1983,20:510-536.

[8]Babuska I,Caloz G,Osbom E.Special finite element methods for a class of second order elliptic problems with rough coefficients[J].SIAM Journal on Numerical Analysis,1994,31:945-981.

[9]Hou T Y,Wu Xiaohui.A multiscale finite element method for elliptic problems in composite materials and porous media[J].Journal of Computational Physics,1997,134(1):169-189.

[10]Kamiński M.Sensitivity and randomness in homogenization of periodic fiber-reinforced composites via the response function method[J].International Journal of Solids and Structures,2009,46(3-4):923-937.

[11]Sakata S,Ashida F,Kojima T.Stochastic homogenization analysis on elastic properties of fiber reinforced composites using the equivalent inclusion method and perturbation method[J].International Journal of Solids and Structures,2008,45(25-26):6553-6565.

[12]方慶軍,王篤波,洪寶寧．基于數(shù)字圖像的土體剪切帶細觀結(jié)構(gòu)試驗研究[J]．華北水利水電大學學報(自然科學版)，2014，35(2)：60-63．

[13]肖建莊,李宏,袁俊強．數(shù)字圖像技術(shù)在再生混凝土性能分析中的應用[J]．建筑材料學報，2014，17(3)：459-464．

[14]Zheng Jianjun,Zhou Xinzhu.Analytical solution for the chloride diffusivity of hardened cement paste[J].Journal of Materials in Civil Engineering,2000,47(5):1001-1027.

(責任編輯:杜明俠)

Analysis of Chloride Diffusivity in Recycled Aggregate Concrete Based on Digital Image Kernel

XIE Caiqi, WU Yuching

(Tongji University, Shanghai 200092, China)

The durability of recycled aggregate concrete (RAC) is easy to be affected by the chloride diffusivity in coastal or deicing salt environment. The coarse aggregates of RAC are always attached by old mortar so that RAC has more complex structures and impressive microscope uncertainties. In order to study the effects of the microscope uncertainties on the chloride diffusivity in recycled aggregate concrete, this paper used scanning electron microscope to observe the meso-scale structure of the recycled coarse aggregates and obtained a series of digital pictures. A digital image driven method was developed, which combined stochastic homogenization method and multi-scale method. A five-phase model of RAC was built and the formulas for studying the chloride diffusivity were derived, the digital pictures were analyzed for the offline part by the digital image driven method. The results of the offline part show that this method is practicable, and the eigenvalues and eigenvectors obtained from the digital image kernel are very desirable with low cost and slight error.

recycled aggregate concrete; multi-scale; homogenization; chloride diffusivity

2015-12-29

國家自然科學基金項目(10972162)。

謝財奇(1991—),男,江西宜春人,碩士研究生,主要從事再生混凝土結(jié)構(gòu)分析方面的研究。E-mail:09qsh@#edu.cn。

10.3969/j.issn.1002-5634.2016.02.011

TV431+.3;TU528

A

1002-5634(2016)02-0062-05