Comparison of Self-consistent Method and Mori Tanaka Method for Predicting Cement Paste Creep
-
摘要:
为了预测水泥浆体的徐变性能, 并对比不同预测模型的精度。使用微米压痕测试确定了水泥浆体的徐变模量, 使用纳米压痕点阵测试获得微观物相的接触徐变模量。随后基于一系列微观测试获得物相的体积分数。将水泥浆体划分为2个尺度, 分别基于自洽(self-consistent, SC)法和Mori Tanaka法, 采用4种组合方式进行徐变均匀化预测。结果表明, 在凝胶尺度采用Mori Tanaka法, 在浆体尺度采用SC法时, 水泥浆体的长期徐变性能的预测精度较高。该研究为进行水泥基材料徐变均匀化时选择合适的计算方法提供了依据。
Abstract:To predict the creep behavior of cement paste and compare the accuracy of different prediction models, the creep modulus of cement paste was determined using microindentation test and the contact creep modulus of microscopic phases was obtained using grid nanoindentation test. The volume fraction of the phases was subsequently obtained based on a series of microscopic tests. The cement paste was divided into two scales and four combinations were used for the prediction of creep homogenization based on the self-consistent (SC) and Mori Tanaka methods, respectively. Results show that the long-term creep properties of cement paste are predicted with high accuracy when the Mori Tanaka method is used at the gel scale and the SC method is used at the slurry scale. This study provides a basis for selecting an appropriate calculation method when performing the creep homogenization of cement-based materials.
-
Keywords:
- cement paste /
- creep /
- indentation test /
- homogenization /
- self-consistent (SC) method /
- Mori Tanaka method
-
水泥基材料的徐变是影响结构耐久性和安全性的重要因素，自20世纪以来受到学者们的广泛关注[1-2]。国内外学者对于徐变的研究大多是针对宏观徐变试验，然而宏观试验往往持续周期为几年甚至几十年，这大大增加了时间成本[3]。
水泥基材料的宏观徐变起源于微观，在微观层面上研究水泥基材料徐变已经取得一定进展[3-4]。有学者应用纳米压痕技术在纳米尺度上对微观物相进行徐变测试并得到诸多有进步意义的徐变认识，如得到水化硅酸钙凝胶(calcium silicate hydrates, C-S-H)是水泥基材料的主要徐变承担者[5]；C-S-H凝胶会因为自身堆积密度的差异导致不同的徐变性能[6]；微观物相的压痕徐变性质与土体的徐变有一定的相似[7]。在诸多压痕徐变研究中，Vandamme等[6]基于纳米压痕试验推导出可以很好地表征物相徐变行为的接触徐变模量。接触徐变模量可以很好地表达物相以及材料所具有的抵抗徐变变形的属性[8]。随后，在微观物相的接触徐变模量的基础上应用均匀化计算，可以经过多尺度模型得到硬化浆体乃至混凝土的徐变[9]。
目前，在水泥基材料领域应用较为普遍的均匀化方法是自洽(self-consistent, SC)法[10]和Mori Tanaka法[11]。2种均匀化方法对水泥基材料的力学性质的可行性已得到证明[12-13]。同时，有学者将以上均匀化方法应用到水泥基材料的徐变性能的研究中[12]。但是，在当前的徐变均匀化研究中，一方面，由于不同尺度的物相之间具有较大的性能差异，对不同尺度的物相进行模型划分未有一致结论。另一方面，当前大都使用单一的均匀化方法，这可能会造成较大的预测偏差。因此，有必要就均匀化方法在计算水泥基材料徐变行为的应用方面进行更深入的研究。
针对以上问题，本研究采用Vanadamme等[6]提出的徐变压痕测试获得了水泥浆体和微观物相的徐变性能。随后确定了各种物相的体积分数。基于2个尺度的均匀化计算，分析SC法和Mori Tanaka法2种均匀化方法的异同点以及精度问题。
1. 试验材料与方法
1.1 原材料和试样制备
试验使用的胶凝材料为PO 42.5级别的普通硅酸盐水泥，化学组成及主要物理性质见表 1。试验用水为普通自来水，采用0.35的水灰比。将水泥与水按比例混合后搅拌至均匀状态，灌入直径为20 mm的足够长度的塑料管内，两端封闭后放入温度为(20±2)℃的养护室内养28 d后，取出制作微观测试试样。从塑料管上切割厚度为8 mm的薄片放入硅胶套内，倒入配置好的冷凝树脂冷却24 h后制得原始试样。然后依次使用400、1 200、2 000目的金相砂纸对原始试样打磨，随后使用5.00、1.00、0.25 μm的金相抛光膏在抛光机上依次抛光40 min，最后在无水乙醇超声波清洗液中清洗后制得压痕测试和BSE测试的试样。
表 1 PO 42.5水泥的化学组成与物理性质Table 1. PO 42.5 Chemical composition and physical properties of cementw(CaO)/% w(Al2O3)/% w(Fe2O3)/% w(MgO)/% w(Na2O)/% w(SiO2)/% w(SO3)/% ρ/(g·cm-3) 比表面积/(m2·g-1) 64.04 5.2 3.26 1.47 0.48 21.48 0.3 3.12 1.16 1.2 压痕测试
1.2.1 微米压痕
对养护28 d后的浆体进行微米压痕测试以获得硬化浆体的徐变模量。采用美国安捷伦公司生产的G200型号纳米压痕仪进行测试，该仪器压头为Berkovich压头。微米压痕点阵规格设置为5×5点阵。控制压痕间距为500 μm以保证相邻压痕互不影响。微米压痕测试控制加、卸载时间均为10 s，持载180 s。测试的峰值荷载为1.5 N。该测试荷载水平及持载时间均可以有效表征硬化浆体的匀质徐变特性[14]。
徐变模量计算过程如下[5]。首先，对于Berkovich压头，当压头与材料表面发生接触时，根据压头的位移与时间、荷载与时间的变化关系，可以提取出表征材料徐变特性的徐变柔度函数L(t)，且在加载时间t=0时刻，柔度完全由材料的弹性特性决定，即L(0)=1/M，M表示压痕点处测得的弹性模量。徐变柔度函数表达式为
$$L(t)-L(0)=L(t)-\frac{1}{M}=\frac{2 r_{\mathrm{u}} \Delta h(t)}{P_{\mathrm{m}}}$$ (1) 式中：ru为卸载时压头接触投影面积的等效压痕半径；h表示压头位移；Δh(t)为持载阶段压头位移随时间的变化；Pm表示峰值荷载。
应用对数型公式Δh(t)=p1ln(p2t+1)+p3t+p4对压痕测试得到的持载阶段压痕深度变化曲线Δh(t)进行拟合。式中的p1、p2、p3和p4为拟合参数。因为参数p3和p4与材料的徐变性能无关，且在实际数据拟合过程中过小，在计算过程中可以忽略，因此，接触徐变函数可以表示为
$$L(t)-L(0)=\frac{2 r_{\mathrm{u}}\left[p_1 \ln \left(p_2 t+1\right)\right]}{P_{\mathrm{m}}}$$ (2) 同时，徐变函数也有如下表达式
$$L(t)-L(0)=\frac{\ln (t / \tau+1)}{C}$$ (3) 式中：τ表示材料开始展现对数型徐变特性时的特征时间，即材料开始发生偏徐变的时间点；C代表徐变模量，能够表征材料固有徐变属性，徐变模量值越大，材料徐变度越小，即材料抵抗徐变变形的能力越强。最后，结合式(2)(3)得到徐变模量的公式为
$$C=\frac{P_{\mathrm{m}}}{2 r_{\mathrm{u}} p_1}$$ (4) 1.2.2 纳米压痕
为了得到硬化浆体中微观物相的徐变特性，采用与微米压痕测试相同的仪器进行纳米压痕测试。在试样上随机选取位置采用10×10的点阵进行纳米压痕测试。在进行压痕测试时，首先采用线性加载的形式，10 s内达到Pm，本研究中Pm取值为2 mN。随后进行180 s的保载阶段，最后线性卸载10 s后使荷载达到0。Pm取值为2 mN时可以很好地记录物相的徐变行为。而且此荷载下的压痕深度为155~245 nm，按照压痕影响面积半径与压痕深度为4倍关系[15]，控制压痕间距为15 μm。由此形成的压痕点阵不但满足测试物相的深度，而且可以较好涵盖所有微观物相。基于纳米压痕测试，使用1.2.1节给出的徐变模量计算方法，可以求得物相的徐变性能。最后使用数据解卷积处理方法，从总的物相徐变模量集合中分理出具体物相的徐变性质[16]。
1.3 物相体积分数测试
使用压汞法测试获得水泥浆体中孔隙的体积分数，采用可测孔径范围0.003~400.000 μm的AutoPore Ⅳ 9500型号的压汞仪进行测试。使用Q5000IR型热重(thermogravimetry, TG)测试仪测定浆体中的氢氧化钙(calcium hydroxide, CH)体积含量，该仪器灵敏度＜0.1 μg，质量称量的准确度是±0.1%。利用QUANTA 250型号的场发射环境扫描电镜(scanning electron microscope, SEM)仪器获得背散射电子(backscattering electron, BSE)图像。测试时的分辨率为2.5 nm，测试时加速电压为15 kV，工作距离为10~15 mm。
2. 结果和讨论
2.1 硬化浆体的徐变模量
微米压痕测试点阵如图 1(a)所示，可以看到每个压痕点之间距离适中，不会互相影响。从微米压痕测试点阵中随机选取的2个压痕点，分别记为压痕点1和压痕点2，其BSE图像见图 1(b)(c)。可以看到每个压痕点覆盖了足够多的物相(BSE图像中不同颜色表示不同的物相)。
![]()
图 1 微米压痕点阵以及压痕点BSE图像Figure 1. Microindentation dot matrix and BSE images of indentation dot图 2所示为根据微米压痕徐变测试结果得到的徐变柔度函数随着时间的变化曲线。可以看到在40 s之前，已经完成了大部分徐变，表明水泥浆体徐变主要发生在试验早期。基于图 2所示函数，使用1.2.1节给出的徐变模量计算方法，得到水泥浆体试样徐变模量的平均值为166.3 GPa。
![]()
图 2 微米压痕测试在持载阶段所得水泥浆体徐变函数Figure 2. Creep function of cement paste obtained by micrometer indentation test in the holding stage2.2 微观物相的性质
2.2.1 微观物相的接触徐变模量
已有研究表明，在水泥浆体中主要存在的物相为C-S-H凝胶、CH、未水化的水泥颗粒(unhydrated clinkers, UC)以及孔隙，凝胶根据密度的不同分为低密度凝胶(low density C-S-H, LD C-S-H)和高密度凝胶(high density C-S-H, HD C-S-H)[6, 17]。因此本研究主要针对以上物相进行后续研究和分析。图 3(a)给出的纳米压痕测试全过程曲线中，可以看到未水化的水泥颗粒、CH、HD C-S-H和LD C-S-H的压痕深度呈现递增趋势。分析上述物相在持载阶段的接触徐变函数随时间的变化如图 3(b)所示，发现UC、CH、HD C-S-H和LD C-S-H的徐变逐渐增大。
![]()
图 3 纳米压痕测试所得物相的荷载位移曲线以及持载阶段的接触徐变函数Figure 3. Load displacement curves of material phases obtained from nanoindentation tests and contact creep functions in the holding phase至于孔隙相，由于其无法直接与压痕探针接触，因此不能直接获得其徐变性能。但是，在水泥浆体中，孔隙相并非独立存在，其通常会与水泥浆体形成混合体[15]。因此，在考虑孔隙的徐变性能时，可以将压痕探针与浆体接触点处含有孔隙的压痕点记为特征孔隙相的压痕测试，所得测试结果记为孔隙的特征接触徐变模量。而且，根据图 4所示的压痕点在加载阶段的荷载-位移曲线，可以发现孔隙相在加载阶段均会出现突变点。
![]()
图 4 特征孔隙相的纳米压痕全过程曲线Figure 4. Nanoindentation full process curve of characteristic pore phase最后，根据纳米压痕点阵测试结果以及解卷积数据处理，得到水泥浆体中的各种物相的接触徐变模量如表 2所示。
表 2 微观物相的接触徐变模量Table 2. Contact creep modulus of microscopic phasesGPa 孔隙相 LD C-S-H HD C-S-H CH UC 52.6±17.4 97.7±26.8 192.6±45.2 394.5±48 896±197.3 2.2.2 微观物相的体积分数
使用压汞法测试得到水泥浆体的孔隙率为12.1%。TG测试得到的质量随温度的变化曲线如图 5所示。可以看到，在400~500 ℃，质量损失曲线会出现明显的下降峰段，此温度段即为CH的受热分解段，此时的质量损失即为水分蒸发的质量。根据CH受热分解的公式，可以计算得到水泥浆体中CH的质量分数为16.9%。随后，根据质量分数与体积分数的换算关系，结合本研究中水泥浆体的实测密度为1 929 kg/m3以及CH的密度为2 240 kg/m3[18]，得到CH在水泥浆体中的体积分数为14.6%。
![]()
图 5 水泥浆体的质量随温度升高的变化曲线Figure 5. Variation curve of the mass of cement paste with increasing temperature图 6所示为基于BSE图像，使用MATLAB软件进行阈值分割计算未水化的水泥颗粒在水泥浆体中的体积分数的示意图。因为在BSE图像中，颜色最亮的即为未水化的水泥颗粒[19]。统计BSE图像中的所有灰度值得像素点个数，得到图 6(b)所示的结果。随即将像素点分布图中在灰度值较大处的拐点即为未水化的水泥颗粒与其余物相的灰度值分割点。以该分割点的灰度值为阈值分割点，得到图 6(c)所示的阈值分割结果，图中灰色部分即为未水化的水泥颗粒。统计灰色部分占总面积的比例即为未水化的水泥颗粒的体积分数。为了消除实验和计算误差，本研究处理了30张BSE图像作为最终的计算结果。最后得到未水化的水泥颗粒在水泥浆体中的体积分数为13.3%。
![]()
图 6 水泥浆体的BSE图像以及二值化处理示意Figure 6. BSE images of cement paste and the diagram of binarization processing将浆体中所有物相的占比即为100%，根据以上3种物相的计算结果，可以得到C-S-H凝胶的体积分数为60%。为了得到LD C-S-H和HD C-S-H的体积分数，对去掉孔隙相的压痕点阵测试所得的接触徐变模量进行解卷积处理，结果如图 7所示。统计图中2个拟合曲线的面积比，即为2种凝胶的比例。最终得到LD C-S-H凝胶的体积分数为34.2%，HD C-S-H凝胶的体积分数为25.8%。统计5种物相在水泥浆体中的体积分数如表 3所示。
![]()
图 7 解卷积处理获得LD C-S-H和HD C-S-H的比值Figure 7. Data deconvolution process to obtain the ratio of LD C-S-H and HD C-S-H表 3 微观物相的体积分数Table 3. Volume fraction of microscopic phases% 孔隙相 LD C-S-H HD C-S-H CH UC 12.1 34.2 25.8 14.6 13.3 2.3 不同均匀化计算方法及组合的比较
2.3.1 多尺度模型
通常在对水泥基材料进行均匀化预测时会按照尺寸大小划分多个尺度。本文采用当前较为通用的Ulm[20]多尺度模型方法将硬化浆体划分为2个尺度，如图 8所示。
![]()
图 8 水泥浆体的SC法模型和Mori Tanaka法模型的示意Figure 8. Schematic diagram of the self-consistent method model and the Mori Tanaka method model for cement paste2.3.2 均匀化理论与方法
混凝土作为一种黏弹性材料，其徐变均匀化处理是仿照材料弹性问题来解决的。当材料在长期受力时，徐变的表达式为$\varepsilon(t)=\int_0^t J(t-\theta) \sigma^{\prime}(\theta) \mathrm{d} \theta$，式中：θ表示时间的积分变量；J表示徐变柔度；ε(t)为随时间变化的应变；σ表示应力。松弛的表达式为$\sigma(t)=\int_0^t R(t-\theta) \varepsilon^{\prime}(\theta) \mathrm{d} \theta$，式中：R表示松弛模量；σ(t)为随时间变化的应力；ε代表应变。根据拉普拉斯转换定义和卷积定理，对上述2个式子进行拉普拉斯转换得到
$$\left\{\begin{array}{l}\varepsilon^*(p)=J^*(p) \sigma^*(p) \\ \sigma^*(p)=R^*(p) \varepsilon^*(p) \\ p J^*(p)=1 / p R^*(p)\end{array}\right.$$ (5) 式中：p表示拉普拉斯变换的复频域；*表示经过laplace-carson变换后的表达形式。由于式(5)中的表达形式与胡克定律的形式相似，因此仿照弹性问题来处理黏弹性问题在理论上是可行的。
本研究分别采用SC法和Mori Tanaka法在2个尺度上进行均匀化处理，按照图 8所示2种模型在每个尺度上的物相进行均匀化计算。均匀化后的均质体由均质体的体积模量和剪切模量得来，而均质体的体积模量和剪切模量又由2种均匀化方法中的基质与夹杂体的体积模量与剪切模量而来[6]。
在使用SC法计算时, 有效剪切柔度(Jhd)和有效体积柔度(Jhv)为
$$J_{\mathrm{h}}^{\mathrm{d}}=\frac{\sum\limits_{i=1}^n \varphi_i J_i^{\mathrm{d}}\left[1+\alpha\left(J_i^{\mathrm{d}} / J_{\mathrm{h}}^{\mathrm{d}}-1\right)\right]^{-1}}{\sum\limits_{i=1}^n \varphi_i\left[1+\alpha\left(J_i^{\mathrm{d}} / J_{\mathrm{h}}^{\mathrm{d}}-1\right)\right]^{-1}}$$ (6) $$J_{\mathrm{h}}^{\mathrm{v}}=\frac{\sum\limits_{i=1}^n \varphi_i J_i^{\mathrm{v}}\left[1+\beta\left(J_i^{\mathrm{v}} / J_{\mathrm{h}}^{\mathrm{v}}-1\right)\right]^{-1}}{\sum\limits_{i=1}^n \varphi_i\left[1+\beta\left(J_i^{\mathrm{v}} / J_{\mathrm{h}}^{\mathrm{v}}-1\right)\right]^{-1}}$$ (7) 其中
$$\begin{gathered}\alpha=6\left(J_{\mathrm{h}}^{\mathrm{v}}+2 J_{\mathrm{h}}^{\mathrm{d}}\right) / 5\left(3 J_{\mathrm{h}}^{\mathrm{v}}+4 J_{\mathrm{h}}^{\mathrm{d}}\right) \\ \beta=3 J_{\mathrm{h}}^{\mathrm{v}} /\left(3 J_{\mathrm{h}}^{\mathrm{v}}+4 J_{\mathrm{h}}^{\mathrm{d}}\right)\end{gathered}$$ 式中：角标i代表夹杂物相；φi是第i夹杂的体积分数。基于水泥基材料的长期对数型假设的拉普拉斯形式有
$$\begin{aligned} J_i^{\mathrm{d}} & =-\ln p /\left(p C_i^{\mathrm{d}}\right) \\ J_i^{\mathrm{v}} & =-\ln p /\left(p C_i^{\mathrm{v}}\right) \\ J_{\mathrm{h}}^{\mathrm{d}} & =-\ln p /\left(p C_{\mathrm{h}}^{\mathrm{d}}\right) \\ J_{\mathrm{h}}^{\mathrm{v}} & =-\ln p /\left(p C_{\mathrm{h}}^{\mathrm{v}}\right)\end{aligned}$$ 将假设带入式(6)(7)得到SC法预测的均匀化徐变剪切模量(Chd)和均匀化徐变体积模量(Chv)的拉普拉斯变换公式为
$$\begin{gathered}C_{\mathrm{h}}^{\mathrm{d}}=\frac{\sum\limits_{i=1}^n \varphi_i C_i^{\mathrm{d}}\left[1+\alpha\left(C_i^{\mathrm{d}} / C_{\mathrm{h}}^{\mathrm{d}}-1\right)\right]^{-1}}{\sum\limits_{i=1}^n \varphi_i\left[1+\alpha\left(C_i^{\mathrm{d}} / C_{\mathrm{h}}^{\mathrm{d}}-1\right)\right]^{-1}}, \\ \alpha=\frac{6\left(1 / C_{\mathrm{h}}^{\mathrm{d}}+2 / C_{\mathrm{h}}^{\mathrm{v}}\right)}{5\left(3 / C_{\mathrm{h}}^{\mathrm{d}}+4 / C_{\mathrm{h}}^{\mathrm{v}}\right)}\end{gathered}$$ (8) $$\begin{gathered}C_{\mathrm{h}}^{\mathrm{v}}=\frac{\sum\limits_{i=1}^n \varphi_i C_i^{\mathrm{v}}\left[1+\beta\left(C_i^{\mathrm{v}} / C_{\mathrm{h}}^{\mathrm{v}}-1\right)\right]^{-1}}{\sum\limits_{i=1}^n \varphi_i\left[1+\beta\left(C_i^{\mathrm{v}} / C_{\mathrm{h}}^{\mathrm{v}}-1\right)\right]^{-1}}, \\ \beta=\frac{3 / C_{\mathrm{h}}^{\mathrm{d}}}{3 / C_{\mathrm{h}}^{\mathrm{d}}+4 / C_{\mathrm{h}}^{\mathrm{v}}}\end{gathered}$$ (9) 在Mori Tanaka法中，有效剪切柔度和体积柔度为
$$J_{\mathrm{h}}^{\mathrm{d}}=\frac{\varphi_{\mathrm{m}} J_{\mathrm{m}}^{\mathrm{d}}+\sum\limits_{i=1}^n \varphi_i J_i^{\mathrm{d}}\left[1+\alpha\left(J_i^{\mathrm{d}} / J_{\mathrm{m}}^{\mathrm{d}}-1\right)\right]^{-1}}{\varphi_{\mathrm{m}}+\sum\limits_{i=1}^n \varphi_i\left[1+\alpha\left(J_i^{\mathrm{d} /} / J_{\mathrm{m}}^{\mathrm{d}}-1\right)\right]^{-1}}$$ (10) $$J_{\mathrm{h}}^{\mathrm{v}}=\frac{\varphi_{\mathrm{m}} J_{\mathrm{m}}^{\mathrm{v}}+\sum\limits_{i=1}^n \varphi_i J_i^{\mathrm{v}}\left[1+\beta\left(J_i^{\mathrm{v}} / J_{\mathrm{m}}^{\mathrm{v}}-1\right)\right]^{-1}}{\varphi_{\mathrm{m}}+\sum\limits_{i=1}^n \varphi_i\left[1+\beta\left(J_i^{\mathrm{v}} / J_{\mathrm{m}}^{\mathrm{v}}-1\right)\right]^{-1}}$$ (11) 式(10)(11)中α=6(Jmv+2Jmd)/(3Jmv+4Jmd)，β=3Jmv/(3Jmv+4Jmd)，角标m表示基质。依据水泥基材料的长期对数型假设的拉普拉斯形式，将假设带入公式(10)和(11)得到Mori Tanaka法预测水泥基材料徐变剪切徐变模量和体积徐变模量的拉普拉斯转换公式为
$$\begin{gathered}C_{\mathrm{h}}^{\mathrm{d}}=\frac{\varphi_m C_m^{\mathrm{d}}+\sum\limits_i^n \varphi_i C_i^{\mathrm{d}}\left[1+\alpha\left(C_i^{\mathrm{d}} / C_m^{\mathrm{d}}-1\right)\right]^{-1}}{\varphi_m+\sum\limits_i^n \varphi_i\left[1+\alpha\left(C_i^{\mathrm{d}} / C_m^{\mathrm{d}}-1\right)\right]^{-1}}, \\ \alpha=\frac{6\left(C_m^{\mathrm{v}}+2 C_m^{\mathrm{d}}\right)}{3 C_m^{\mathrm{v}}+4 C_m^{\mathrm{d}}}\end{gathered}$$ (12) $$\begin{gathered}C_{\mathrm{h}}^{\mathrm{v}}=\frac{\varphi_m C_m^{\mathrm{v}}+\sum\limits_i^n \varphi_i C_i^{\mathrm{v}}\left[1+\beta\left(C_i^{\mathrm{v}} / C_m^{\mathrm{v}}-1\right)\right]^{-1}}{\varphi_m+\sum\limits_i^n \varphi_i\left[1+\beta\left(C_i^{\mathrm{v}} / C_m^{\mathrm{v}}-1\right)\right]^{-1}}, \\ \beta=\frac{3 C_m^{\mathrm{v}}}{3 C_m^{\mathrm{v}}+4 C_m^{\mathrm{d}}}\end{gathered}$$ (13) 弹性体中体积模量(K)、剪切模量(G)与M的关系式为M=4G(3K+G)/(3K+4G)，对该关系式应用拉普拉斯转化得到SC法和Mori Tanaka法中徐变模量的表达式为C=4Chd(3Chv+Chd)/(3Chv+4Chd)。
2.3.3 均匀化计算结果
将3.2节中得到的各物相的接触徐变模量值和各物相的体积分数分别带入第3.3.2节的均匀化预测公式中，即可得到水泥浆体试样的徐变多尺度预测结果。为了充分分析SC法和Mori Tanaka法对于徐变预测结果的影响，本研究采用4种计算方式进行计算。分别是在2个尺度上均采用SC法或者Mori Tanaka法进行计算，得到水泥浆体的徐变模量预测值分别为185.2、183.6 GPa。在凝胶尺度上采用SC法(Mori Tanaka法)，在浆体尺度上采用Mori Tanaka法(SC法)进行计算，得到水泥浆体的徐变模量预测值分别为176.9、176.5 GPa。
可以发现，4种均匀化计算过程的预测结果均大于微米压痕实测的水泥浆体的徐变模量。其中SC法、Mori Tanaka法、SC法+Mori Tanaka法以及Mori Tanaka法+SC法的预测结果与实测值的误差分别为11.4%、10.4%、6.4%和6.1%。可以看到，SC法和Mori Tanaka法的徐变预测所得结果较为接近。在2个尺度上使用SC法和Mori Tanaka法混合计算时，所得计算结果与实测值更为接近。这可能是因为在不同尺度上，由于物相的存在形式存在差异，而不同的均匀化方法的使用条件也存在一定差别。其中，预测效果最优的均匀化组合是Mori Tanaka法+SC法。这可能是因为，LD C-S-H和HD C-S-H的性能差异较大，即两者的体积分数以及存在形式与基质包裹夹杂较为相近。在浆体尺度上，C-S-H凝胶、CH和UC等物相之间的性能差异较小，可将其所有物相构成的总体视为基质，所有物相视为夹杂相处理。因此，在今后进行水泥基材料的徐变多尺度预测时，应充分结合每个尺度上的物相的具体特征去进行均匀化方法的选择。
3. 结论
1) 采用微米压痕测试可以获得硬化水泥浆体的徐变性能，采用点阵纳米压痕测试可以获得微观物相的徐变性能。在水泥浆体中，可以将孔隙相和水泥浆体的混合材料记为特征孔隙相，纳米压痕测试可以得到特征孔隙相的徐变性能。在水泥浆体中，特征孔隙、LD C-S-H凝胶、HD C-S-H凝胶、CH和未水化水泥颗粒的接触徐变模量值依次增加。
2) 在确定水泥浆体中的LD C-S-H凝胶和HD C-S-H凝胶的体积分数时，可以通过先确定CH、孔隙等物相，然后根据纳米压痕测试结果，采用数据解卷积方法，得到2种密度凝胶的比值，最后根据总凝胶的体积分数确定。
3) 在对水泥浆体进行徐变均匀化时，可以划分为2个尺度。凝胶尺度上包含LD C-S-H和HD C-S-H。在水泥浆体尺度上包含C-S-H凝胶、CH、未水化的水泥颗粒以及孔隙。
4) 在水泥浆体的多尺度徐变预测中，SC法和Mori Tanaka法的预测结果较为接近。在不同尺度上将SC法和Mori Tanaka法结合使用时，会显著提升徐变预测结果的精度。在徐变预测时，应充分考虑不同尺度的物相特征，选择合适的均匀化方法。
-
![]()
图 1 微米压痕点阵以及压痕点BSE图像
Figure 1. Microindentation dot matrix and BSE images of indentation dot
![]()
图 2 微米压痕测试在持载阶段所得水泥浆体徐变函数
Figure 2. Creep function of cement paste obtained by micrometer indentation test in the holding stage
![]()
图 3 纳米压痕测试所得物相的荷载位移曲线以及持载阶段的接触徐变函数
Figure 3. Load displacement curves of material phases obtained from nanoindentation tests and contact creep functions in the holding phase
![]()
图 4 特征孔隙相的纳米压痕全过程曲线
Figure 4. Nanoindentation full process curve of characteristic pore phase
![]()
图 5 水泥浆体的质量随温度升高的变化曲线
Figure 5. Variation curve of the mass of cement paste with increasing temperature
![]()
图 6 水泥浆体的BSE图像以及二值化处理示意
Figure 6. BSE images of cement paste and the diagram of binarization processing
![]()
图 7 解卷积处理获得LD C-S-H和HD C-S-H的比值
Figure 7. Data deconvolution process to obtain the ratio of LD C-S-H and HD C-S-H
![]()
图 8 水泥浆体的SC法模型和Mori Tanaka法模型的示意
Figure 8. Schematic diagram of the self-consistent method model and the Mori Tanaka method model for cement paste
表 1 PO 42.5水泥的化学组成与物理性质
Table 1 PO 42.5 Chemical composition and physical properties of cement
w(CaO)/% w(Al2O3)/% w(Fe2O3)/% w(MgO)/% w(Na2O)/% w(SiO2)/% w(SO3)/% ρ/(g·cm-3) 比表面积/(m2·g-1) 64.04 5.2 3.26 1.47 0.48 21.48 0.3 3.12 1.16 
下载: 导出CSV
表 2 微观物相的接触徐变模量
Table 2 Contact creep modulus of microscopic phases
GPa 孔隙相 LD C-S-H HD C-S-H CH UC 52.6±17.4 97.7±26.8 192.6±45.2 394.5±48 896±197.3
下载: 导出CSV
表 3 微观物相的体积分数
Table 3 Volume fraction of microscopic phases
% 孔隙相 LD C-S-H HD C-S-H CH UC 12.1 34.2 25.8 14.6 13.3
下载: 导出CSV
-
[1] SU L, WANG Y, MEI S, et al. Experimental investigation on the fundamental behavior of concrete creep[J]. Construction and Building Materials, 2017, 152(15): 250-258.
[2] 张震雷, 史美伦, 杜江. 聚合物改性水泥浆体黏弹性的频率响应(Ⅰ): 动态力学分析[J]. 建筑材料学报, 2022, 25(4): 339-343. https://www.cnki.com.cn/Article/CJFDTOTAL-JZCX202204002.htm ZHANG Z L, SHI M L, DU J. Frequency response of viscoelasticity of polymer modified cement paste (Ⅰ): dynamical mechanical analysis[J]. Journal of Building Materials, 2022, 25(4): 339-343. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JZCX202204002.htm
[3] 宋昊, 曾晓辉, 谢友均, 等. 水泥乳化沥青砂浆的徐变模型[J]. 建筑材料学报, 2020, 23(6): 1305-1312. doi: 10.3969/j.issn.1007-9629.2020.06.007 SONG H, ZENG X H, XIE Y J, et al. Creep model of cement emulsified asphalt mortar[J]. Journal of Building Materials, 2020, 23(6): 1305-1312. (in Chinese) doi: 10.3969/j.issn.1007-9629.2020.06.007
[4] 曹国辉, 胡佳兴, 张锴, 等. 混凝土徐变预测模型修正分析[J]. 建筑结构, 2014, 44(5): 45-49. https://www.cnki.com.cn/Article/CJFDTOTAL-JCJG201405013.htm CAO G H, HU J X, ZHANG K, et al. Correction analysis on prediction model of concrete creep[J]. Building Structure, 2014, 44(5): 45-49. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JCJG201405013.htm
[5] VANDAMME M, ULM F J. Nanogranular origin of concrete creep[J]. PNAS, 2009, 106(26): 10552-10557. doi: 10.1073/pnas.0901033106
[6] VANDAMME M, ULM F J. Nanoindentation investigation of creep properties of calcium silicate hydrates[J]. Cement and Concrete Research, 2013, 52: 38-52. doi: 10.1016/j.cemconres.2013.05.006
[7] 魏亚, 梁思明, 高翔. 基于土体isotache模型的硬化水泥净浆微观徐变预测[J]. 复合材料学报, 2018, 35(6): 1609-1618. https://www.cnki.com.cn/Article/CJFDTOTAL-FUHE201806031.htm WEI Y, LIANG S M, GAO X. Prediction of micro creep behavior of hardened cement paste based on the isotache model of clays[J]. Acta Materiae Composite Sinica, 2018, 35(6): 1609-1618. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-FUHE201806031.htm
[8] WEI Y, LIANG S, GAO X. Indentation creep of cementitious materials: experimental investigation from nano to micro length scales[J]. Construction and Building Materials, 2017, 143: 222-233. doi: 10.1016/j.conbuildmat.2017.03.126
[9] YU P, DUAN Y, CHEN E, et al. Microstructure-based homogenization method for early-age creep of cement paste[J]. Construction and Building Materials, 2018, 188: 1193-1206. doi: 10.1016/j.conbuildmat.2018.08.196
[10] KRONER E. Berechnung der elastischen konstanten des vielkristalls aus den Konstanten des Einkristalls[J]. Zeitschrift für Physik, 1958, 151: 504-518. doi: 10.1007/BF01337948
[11] MORI T, TANAKA K. Average stress in matrix and average elastic energy of materials with misfitting inclusions[J]. Acta Metallurgica, 1973, 21(5): 571-574. doi: 10.1016/0001-6160(73)90064-3
[12] ZHAN Q, LEROY R, VANDAMME M, et al. Long-term creep properties of cementitious materials: comparing microindentation testing with macroscopic uniaxial compressive testing[J]. Cement and Concrete Research, 2014, 58: 89-98. doi: 10.1016/j.cemconres.2014.01.004
[13] 张庆华, 刘西拉, 朱振宇. 用自洽方法计算混凝土的弹性模量[J]. 上海交通大学学报, 2001, 35(10): 1503-1506. https://www.cnki.com.cn/Article/CJFDTOTAL-SHJT200110014.htm ZHANG Q H, LIU X L, ZHU Z Y. Applying self-consistent method to calculate concrete modulus[J]. Journal of Shanghai Jiao Tong University, 2001, 35(10): 1503-1506. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SHJT200110014.htm
[14] 魏亚, 梁思明, 高翔. 应变率对硬化水泥净浆微观力学性能及徐变行为的影响[J]. 硅酸盐学报, 2017, 45(11): 1605-1612. https://www.cnki.com.cn/Article/CJFDTOTAL-GXYB201711008.htm WEI Y, LIANG S M, GAO X. Strain rate effect on micro-mechanical properties and creep behavior of hardened cement paste[J]. Journal of the Chinese Ceramic Society, 2017, 45(11): 1605-1612. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GXYB201711008.htm
[15] WEI Y. GAO X, HUANG W. Effect of individual phases on multiscale modeling mechanical properties of hardened cement paste[J]. Construction and Building Materials, 2017, 153(30): 25-35.
[16] 李一凡, 管学茂, 刘松辉, 等. 压痕点数及解卷积法对水泥纳米压痕试验的影响[J]. 建筑材料学报, 2021, 24(2): 291-295. https://www.cnki.com.cn/Article/CJFDTOTAL-JZCX202102010.htm LI Y F, GUAN X M, LIU S H, et al. Effect of number of indentation points and deconvolution method on cement nanoindentation Test[J]. Journal of Building Materials, 2021, 24(2): 291-295. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JZCX202102010.htm
[17] CONSTANTINIDES G, ULM F J. The effect of two types of C-S-H on the elasticity of cement-based materials: Results from nanoindentation and micromechanical modeling[J]. Cement and Concrete Research, 2004, 34(1): 67-80.
[18] TENNISAP D, JENNINGS H M. A model for two types of calcium silicate hydrate in the microstructure of Portland cement pastes[J]. Cement and Concrete Research, 2000, 30(6): 855-863.
[19] 王培铭, 丰曙霞, 刘贤萍. 背散射电子图像分析在水泥基材料微观结构研究中的应用[J]. 硅酸盐学报, 2011, 39(10): 1659-1665. https://www.cnki.com.cn/Article/CJFDTOTAL-GXYB201110024.htm WANG P M, FENG S X, LIU X P. Application of backscattered electron imaging and image analysis in microstructure research on cement-based materials[J]. Journal of the Chinese Ceramic Society, 2011, 39(10): 1659-1665. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GXYB201110024.htm
[20] BERNARD O, ULM F J, LEMARCHAND E. A multiscale micromechanics-hydration model for the early-age elastic properties of cement-based materials[J]. Cement and Concrete Research, 2003, 33(9): 1293-1309.
计量
- 文章访问数: 96
- HTML全文浏览量: 17
- PDF下载量: 32
