吴建营 陈万昕 黄羽立

∗(华南理工大学亚热带建筑科学国家重点实验室,广州 510641)

†(华南理工大学土木工程系,广州 510641)

∗∗(清华大学土木工程系,北京 100084)

引言

混凝土是全世界范围内应用最为广泛的土木工程材料,在建筑结构、桥梁隧道、能源设施(水电大坝、核电安全壳、液化天然气或石油储罐、风电基础等)等土木工程和基础设施建设中发挥着重要作用.然而,一旦混凝土结构在建造期出现裂缝,在服役期载荷和环境作用下这些早龄期裂缝将逐渐扩展和演化,对结构全生命周期的完整性、耐久性和安全性造成严重影响.因此,准确预测早龄期混凝土的抗裂性能并由此控制混凝土结构建造质量,是预防结构先天病患、保证结构优生健康的重要手段.

混凝土浇筑养护过程中,水泥水化过程会释放热量使得混凝土升温,热传导和热对流等机制将导致混凝土温度达到峰值后逐渐降低.这一过程中,与温度有关的热变形以及与水化度有关的自收缩变形会引发较大的拉应力,而此时混凝土力学性能如强度等仍然处于较低水平,导致早龄期混凝土开裂并逐渐扩展和演化[1].上述过程涉及化学、热学、力学等复杂的多场耦合效应,已有理论模型和数值方法尚不能定量预测早龄期混凝土的抗裂性能并准确描述早龄期裂缝的演化过程,这也成为混凝土结构全寿命性能调控研究的主要难点问题之一.

早龄期混凝土裂缝理论模型和数值模拟方法的相关研究可以追溯至20 世纪90 年代.根据所采用的混凝土裂缝模型,大致可以分为以下4 类:

(1)经验或塑性力学模型.早期研究大多采用塑性模型描述开裂后的早龄期混凝土性能,如文献[2-5]等.

(2) 弥散裂缝模型.这方面的代表性研究工作包括De Borst 等[6-8]提出的早龄期混凝土弥散裂缝模型.

(3)损伤力学模型.随着混凝土损伤力学的兴起,若干学者[9-14]将部分经典的损伤模型推广至早龄期混凝土.

(4)微细观力学模型.为了考虑混凝土细观结构特别是粗骨料对早龄期混凝土抗裂性能的影响,部分学者发展了相应的微细观力学模型[15-17].

上述工作为研究早龄期混凝土力学性能奠定了良好的基础.然而,这些研究大多采用较为简单的裂缝模型,也未能考虑裂缝演化与水化反应、热量传输等过程之间的相互耦合,自然也很难描述早龄期混凝土裂缝的演化过程并量化开裂后材料的力学性能.

近年来,法国学者Nguyen 等[18-20]采用脆性断裂相场模型描述开裂混凝土的力学行为,并考虑裂缝与水化反应、热量传输之间的耦合效应,能够较好地描述早龄期的裂缝演化过程.脆性断裂相场模型根植于Francfort 和Marigo[21]提出的线弹性断裂能量变分原理,解决了Griffith 能量方法需要预设裂缝路径的难题;由于将尖锐裂缝正则化为具有一定尺度的裂缝带,并引入在[0,1]之间连续分布的裂缝相场及其梯度描述裂缝状态,脆性断裂相场模型非常便于通过有限元、无网格、物质点等多种数值方法加以实现.可以严格证明[22]:当裂缝尺度趋近于零时,脆性断裂相场模型Γ-收敛于Griffith 线弹性断裂力学.

然而,脆性断裂相场模型仅适用于脆性材料,难以合理反映具有明显软化段的混凝土损伤破坏行为.更糟糕的是,模型的定量分析结果存在严重的裂缝尺度敏感性问题:裂缝尺度越小,结构起裂载荷和峰值载荷越大.因此,若将裂缝尺度视为数值参数,则此类模型与Griffith 线弹性断裂力学类似,同样无法描述完好或仅有弱奇异性固体的裂缝起裂.迫不得已,研究人员只好将裂缝尺度视为材料属性[23-25],这一处理方式仍然无法完全解决裂缝起裂问题[26-27].

近年来,笔者将断裂力学和经典损伤力学有机结合,建立了同时适用于脆性断裂和准脆性破坏的统一相场理论[28].从该理论出发,不仅经典的脆性断裂相场模型[23,29]可以作为特例给出,还在国际上首次提出了一类相场内聚裂缝模型PF-CZM[30-32].理论分析证明:当裂缝尺度趋近于零时,PF-CZM 收敛为一类混合型破坏的内聚裂缝模型[33-34],能够准确预测线性软化以及指数软化、双曲软化、混凝土科内列森软化[35]等非线性软化曲线.该模型同时考虑了基于强度的裂缝起裂准则、基于能量的裂缝扩展准则以及基于变分原理的裂缝扩展方向判据,仅需少量标准材料参数即可定量预测工程结构的损伤破坏行为;特别是,当裂缝尺度小于某一上限时,其取值对脆性断裂和准脆性破坏的分析结果几乎不产生影响.归功于上述优点,统一相场理论和相场内聚裂缝模型一经提出,迅速得到了国内外学者的广泛认可和跟踪研究[36-41].

基于上述考虑,在统一相场理论框架内,本工作拟进一步考虑裂缝相场与水化反应、热量传输等过程之间的双向耦合影响,建立裂缝相场演化特征与混凝土水化度和温度之间的定量联系,提出混凝土化-热-力多场耦合相场内聚裂缝模型,发展相应的多场有限元数值实现算法,并应用于早龄期混凝土的裂缝模拟和抗裂性能预测.

1 化−热−力多场耦合相场内聚裂缝模型

图1 固体内嵌裂缝/界面及其几何正则化[34]Fig.1 A cracking solid and its geometric regularization[34]

相场模型中,尖锐裂缝/界面S 被弥散为如图1所示的裂缝带B ⊆Ω,其尺度记为b>0,并引入相场变量d(x) :B →[0,1]描述裂缝状态;裂缝带B 的外边界记为∂B,其外法向矢量表示为nB.需要指出的是,在固体损伤破坏过程中,裂缝带B 并非预先设定或保持固定,而是遵循自身特定本构行为发生起裂、扩展和演化等过程.类似于位移场,裂缝相场也可以施加某种强迫边界条件:例如,对于弹性区域有d(x)=0;对于初始预设裂缝有d(x)=1 等.

1.1 混凝土开裂过程

首先介绍描述混凝土开裂行为的相场内聚裂缝模型,包括裂缝相场-位移场耦合方程、应力-应变关系以及本构特征函数等.需要指出的是,早龄期混凝土的变形和力学性能与龄期或水化程度、温度等因素有关,裂缝建模时必须合理地加以考虑.

1.1.1 裂缝相场-位移场耦合方程

统一相场理论[28,32]中,准静态条件下开裂固体的行为由如下耦合方程描述

其中,b∗为体积均布力,裂缝相场通量q和裂缝相场源Q分别表示为

式中,应力σ 和裂缝有效驱动力H 由下节的本构关系给出;尺度b>0 为表征裂缝带宽度的数值参数,其值越小,裂缝带宽越小,越接近尖锐裂缝;裂缝几何函数α(d) ∈[0,1]和能量退化函数ω(d) ∈[0,1]分别为裂缝相场d的递增和递减函数,并满足条件[28,43]

这里,引入了归一化参数cα:=,以保证固体完全破坏时单位裂缝面积的耗能为材料断裂能属性.

需要指出,上述相场模型隶属于断裂力学或梯度损伤力学领域的Bourdin 等[29]变分相场断裂模型,与物理学领域常用的Ginzburg-Landau 相变方程[44]、Cahn-Hillard[45]或Allan-Cahn[46]相分离、晶体生长或多相流方程有明显区别,虽然后者也曾经被美国学者Karma 等[47]用于固体裂缝扩展分析.

变分相场断裂模型的理论基础来源于线弹性断裂力学中的Griffith 能量原理[48]以及后续Francfort–Marigo 提出的变分原理[21],并借鉴了图像分割理论中Mumford-Shah 泛函[49]的Ambrosio-Tortorelli 椭圆正则化[50].因此,变分相场断裂模型通常满足Γ-收敛性质,即:当正则化尺度参数b趋近于零时,其近似解在能量泛函意义上收敛于原尖锐裂缝问题的精确解,具体可参见Braide[51]和Bourdin 等[22]中的证明和论述.

除应用范围不同外,以上变分相场断裂模型与Ginzburg-Landau 方程主要有以下区别.

(1)Ginzburg-Landau 方程中的裂缝几何函数通常采用双势阱函数α(d)=d2(1-d)2,该函数满足性质α(0)=α(1)=0,即裂缝相场(或序参数)处于中间过渡状态d∈(0,1)时的裂缝耗散比纯相d=0 或d=1时更大,这一特点与材料相变、多相流或结构拓扑优化等问题较为契合.然而,虽然双势阱函数也满足Γ-收敛特性,但该函数在[0,1]之间不具有单调递增性,也无法区分d=0 和d=1 这两种截然不同的裂缝状态.因此,若将其作为裂缝几何函数,将会导致裂缝面积和能量耗散随裂缝扩展呈现先增大后减小的趋势,这与固体断裂时Griffith 能量原理(即应变能持续减小、裂缝表面能单调增大、最终总能量达到最小值)相矛盾.在断裂力学领域被广泛认可的变分相场断裂模型中,其裂缝几何函数必须为裂缝相场(或序参数)d的单调递增函数,且满足条件(2a),这与裂缝扩展的不可逆特征一致,其具体形式将在后文给出.

(2)Ginzburg-Landau 方程通常会引入与裂缝相场(或序参数)一阶时间导数有关的松弛项或黏性项,当松弛模量趋近于零(或黏性系数趋近无穷大) 时,方程退化为准静态式(1b)的形式.因此,此类模型先天具有一定的动力阻尼效应.类似地,对于动态裂缝扩展问题,变分相场断裂模型中的控制方程(1b)1(下标1 指式(1b)中的第1 式,下同),其等号右边也可以引入与裂缝相场一阶时间导数(裂缝扩展速率)有关的粘性阻尼项或(和) 与裂缝相场二阶时间导数(裂缝扩展加速度)有关的惯性项,以描述裂缝动态扩展时的应变率效应,具体参见作者的相关工作[52-53].这里仅考虑早龄期混凝土裂缝的准静态扩展过程,因此并未引入阻尼项和惯性项.

需要指出的是,这里所采用的相场内聚裂缝模型基于作者2017 年提出的统一相场理论[28],虽也属于变分相场断裂范畴,但与脆性断裂变分相场模型有本质不同:归功于后文给出的本构特征函数,相场内聚裂缝模型Γ-收敛于Barenblatt 内聚裂缝模型[54],而非传统的Griffith 线弹性断裂理论[48].因此,相场内聚裂缝模型不仅同时适用于脆性断裂,也适用于粘聚断裂,且模型预测给出的整体响应(包括裂缝路径和破坏模式等定性结果以及峰值载荷和软化段等定量结果)与裂缝正则化尺度无关.

1.1.2 应力-应变本构关系

为了考虑水化过程对混凝土变形的影响,通常可以将总应变张量ϵ 分解为力学应变ϵm、热膨胀应变ϵt和自收缩应变ϵa三部分之和的形式,即

为简单起见,热膨胀应变ϵt和自收缩应变ϵa均假定为二阶各向同性张量,分别表示为

式中,1为二阶单位张量;ηt为热膨胀系数,表征与当前温度θ 和初始温度θ0之差有关的热膨胀应变;MacAuley 括号定义为〈x〉:=max(x,0);当水化度χ 超过阈值χ0(本文统一取为χ0=0.05) 后,混凝土出现自收缩变形,其大小与系数ηa有关;ξ(χ)为截断线性函数,表示为[55]

其中,混凝土最终水化度χ∞与水灰比w/c之间服从上述指数关系[56].

早龄期混凝土裂缝通常处于受拉张开状态,因此可采用如下各向同性损伤本构关系

式中G0为弹性刚度张量,E0为杨氏模量;为有效应力张量的最大主应力,ft为材料破坏强度.由于龄期的影响,这里杨氏模量和材料破坏强度均为水化度的递增函数,详见1.2.2 节.

1.1.3 本构特征函数

为了反映混凝土开裂引起的应变软化行为,相场内聚裂缝模型采用如下裂缝几何函数和能量退化函数[28,30-31]

对于线性软化和混凝土常用的科内列森软化曲线[35],模型给出的软化曲线如图2 所示;其他如指数软化、双曲线软化等也可以类似给出,见文献[57].

图2 相场内聚裂缝模型PF-CZM 给出的软化曲线[28,32]Fig.2 Softening curves given by the PF-CZM[28,32]

前期研究表明[58]:裂缝尺度b的取值越小,裂缝面积和相应能量耗散的计算精度就越高;极限情况下,裂缝尺度b→0 时,裂缝带宽度也趋近于零,模型退化为一类混合型的内聚裂缝模型[34].另一方面,模型给出的应力-名义位移跳跃曲线σ(w)与裂缝尺度b无关,而仅取决于弹性模量E0、破坏强度ft、断裂能Gf以及软化段参数(即初始斜率k0、裂缝极限位移wc)等材料属性.因此,对于裂缝几何函数和能量退化函数(8),统一相场理论退化为一类相场内聚裂缝模型.于是,相场内聚裂缝模型不仅适用于各类准脆性破坏,还可通过线性软化描述脆性断裂[31],实现了固体材料结构两种典型破坏特征的统一反映;详见文献[28,31-33]中的讨论.

1.2 混凝土水化过程

早龄期混凝土养护过程中,水化反应会释放热量而导致材料升温膨胀;随后,传导和对流等热传输过程导致温度达到最高值后逐渐降低并引发收缩变形.

一方面,裂缝会对混凝土水化反应和热传输过程产生影响;另一方面,混凝土力学性能也会随着龄期或水化度而发生变化,共同构成了一类典型的化-热-力多场耦合问题,如图3 所示.

图3 早龄期混凝土开裂过程中的化-热-力多场耦合Fig.3 Multi-physical coupling in the cracking process of ear-age concrete

1.2.1 水化度场-温度场耦合控制方程

早龄期混凝土水化过程中,化-热耦合问题由如下傅里叶定律和阿伦尼乌斯定律描述

式中,热流(通量)h和水化速率˙χ 分别表示为

其中,ρc为单位体积热容;Θ 为水化反应所释放的热量,p∞为水化反应的潜热;热传导系数假设为二阶各向同性张量且与水化度χ 无关,即κ0=κ01,其中κ0为初始热传导系数;边界∂Ωh上给定的热流记为h∗,边界∂Ωc上与空气温度θa有关的对流/辐射系数记为g∗,剩余边界∂Ωθ:=∂ΩΩhΩc的温度为给定值θ∗;活化能Ea表征热生成速率,R=8.314 J/(K ·mol)为理想气体常数;对于混凝土,水化反应演化函数φ(χ)一般可由如下多项式函数拟合给出[13]

系数bi(i=0,1,2,···,6) 根据等温量热法试验数据拟合确定,如图4 所示.

图4 典型混凝土水化反应演化函数[13]Fig.4 Evolution function of the hydration reaction[13]

由式(11b)可见,早龄期裂缝演化会影响混凝土水化反应和热量传输过程.

1.2.2 混凝土龄期效应

随着龄期增加,混凝土表现出明显的龄期效应,即弹性模量E0、破坏强度ft和断裂能Gf均逐渐增大,而泊松比ν0基本保持不变.

为简化起见,这里假设早龄期混凝土开始形成弹性模量、破坏强度和断裂能的水化度阈值与自收缩应变的水化度阈值χ0取为相同值,即考虑如下混凝土龄期效应[4,59-60]

式中,函数ξ(χ)由式(5)给出;E∞,Gf∞和ft∞分别为水化度达到饱和χ∞时的弹性模量、断裂能和抗拉强度.

上述混凝土力学性能的龄期效应也可以采用不同的水化度阈值或其他函数形式,如文献[19].此时,混凝土的特征长度lch将与水化度χ 有关;相应地,相场内聚裂缝模型中,参数a1也将随之发生变化.另一方面,参数a2和a3仅与软化曲线类型有关,不受材料力学参数变化的影响.

2 数值实现和求解算法

为便于有限元数值实现,通常将化-热-力多场耦合控制方程(1)和式(11)表示为如下弱形式

式中,δθ,δχ,δu,δd分别为任意容许的节点温度、水化度、位移和裂缝相场.与温度和位移条件类似,裂缝相场d(x) 也可以施加诸如d(x∈ΩB)=0 或d(x∈S0)=1 等强制边界条件,这里S0为预制裂缝.

上述弱形式的化-热-力多场耦合问题可以采用多场有限元方法进行空间离散,并选用合理的数值算法进行求解.

2.1 多场有限元数值实现

在进行有限元空间离散时,可以将整个计算区域内所有单元节点均赋予温度、水化度、位移和裂缝相场自由度.为提高计算效率,也可以(不一定必要)将整个计算域Ω 分为两部分[58],即:可能出现裂缝的相场子区域B 和剩余部分ΩB.对于前者,单元节点同时具有温度、水化度、位移和裂缝相场等自由度,而后者的单元节点仅考虑温度、水化度和位移等自由度.已有研究表明[28,30-31,57,61],对于三角形常应变单元,为保证计算精度,相场子区域B 内的单元大小一般取为he≤b/5;四节点双线性单元,其大小可适当放宽至he≤b/3.

有限元空间离散后,温度、水化度、裂缝相场和位移等场变量及其梯度分别表示为(数值实现时,张量采用Voigt 标记)

方程组(16) 通常采用增量迭代法进行求解,即将时间步划分为N个子步.对于某典型子步[Tn,Tn+1],其时间增量为ΔT:=Tn+1-Tn且Tn时刻所有的状态变量已知,在此基础上采用迭代方法求解Tn+1时刻的状态变量.

2.2 控制方程的数值求解

考虑早龄期混凝土开裂过程的时间尺度较水化作用小得多,上述多场耦合问题可以视为化-热耦合以及裂缝相场-位移耦合两个串行过程,分别进行求解后判断是否收敛;如不收敛,则重新进行迭代.上述数值求解过程如流程图5 所示.

图5 化-热-力多场耦合相场内聚裂缝模型数值求解流程图Fig.5 Flow chart in solving chemo-thermo-mechanical coupled problems

2.2.1 水化过程

对于水化过程,采用后退欧拉法进行时间离散

下标n标识上一增量步结束后的状态变量.

水化过程涉及的裂缝相场均采用上一增量步结束后的dn.由此,化-热耦合方程(16)1和(16)2可采用牛顿迭代方法进行求解,其线性化形式表示为

式中,子刚度矩阵分别表示为

上述刚度矩阵不对称,这是耦合问题的特点之一.

2.2.2 开裂过程

对于早龄期混凝土开裂和演化过程,有效应力¯σ需考虑热膨胀应变ϵt和自收缩应变ϵa,这些状态变量可通过水化过程求解得到的温度和水化度计算给出.相应的,裂缝相场-位移场耦合控制方程的求解与纯力学问题基本相同.

在求解耦合方程(16)3和(16)4时,整体牛顿迭代或子问题交错迭代等算法存在收敛性差、计算速度慢等问题.近年来,笔者[62-63]首次将整体BFGS 拟牛顿迭代算法[64]应用于这一耦合方程的求解,在保证算法收敛性的同时,收敛速度和计算效率大幅度提升.这里对该算法做简要介绍.

在整体BFGS 拟牛顿算法中,对于非线性方程组f(z)=0,其状态变量z按如下格式进行更新

整体BFGS 拟牛顿算法可以根据需要设置一定迭代次数(5～10 次)间隔后才对刚度矩阵˜K进行更新计算.因此,虽然所需迭代次数可能比牛顿迭代算法更多(后者不一定能收敛),但BFGS 算法计算过程较为简单且无需每次迭代都更新刚度矩阵,计算效率通常更高.较之子问题交错迭代算法,整体BFGS拟牛顿算法的求解效率提高5～8 倍[62-63,66-67].

3 应用算例

针对上述化-热-力多场耦合相场内聚裂缝模型,本节给出其在早龄期混凝土裂缝演化模拟和抗裂性能预测的验证算例.对于混凝土,采用科内列森软化曲线.这里,二维问题均假定为平面应力状态,网格划分采用四节点双线性插值单元.

由于水化反应过程的控制方程(11)涉及时间导数,在求解过程中按照式(18) 进行时间离散并需要给定节点温度和节点水化度的初始值.算例1 中,节点温度的初始值取为室温20°C,算例2 中,节点温度的初始值取为室温21°C;算例1 和算例2 中,节点水化度的初始值均取为0.

3.1 带方形孔早龄期混凝土圆板收缩开裂[18]

首先考虑如图6 所示的早龄期混凝土圆形板,其外直径为0.5 m,方形孔的尺寸为0.4 m.材料参数由文献[13,18]给出(详见表1);所采用混凝土水化反应演化函数如图4 所示.由于模型具有对称性,建模时仅考虑右上角1/4 模型,并分别在模型左侧和底部施加X向和Y向位移约束.

图6 带方形孔圆板收缩开裂:试件几何尺寸、边界和环境条件[18]Fig.6 A circular plate with a square hole:Specimen geometry,boundary and environment conditions[18]

表1 带方形孔早龄期混凝土圆形板模型参数取值Table 1 Model parameters for the cracking analysis of an early-age circular concrete?plate with a square hole

首先考虑裂缝尺度对数值模拟结果的影响.模拟中,采用b=5 mm 和b=10 mm 两种裂缝尺度,相应的单元大小取为he=b/5.模拟给出典型时刻的裂缝相场云图如图7 所示.可以看出,裂缝尺度b的取值大小仅改变损伤带的宽度,而对裂缝路径和破坏模式几乎没有任何影响.

图7 带方形孔圆板收缩开裂:不同时刻的裂缝相场云图(上排:b=10 mm;下排: b=5 mm)Fig.7 A circular plate with a square hole:Predicted damage profiles at various instants(top: b=10 mm;bottom: b=5 mm)

图8 给出了两种尺度参数情况下圆形板内、外边缘和中间位置处混凝土的温度、水化度、弹性模量、破坏强度以及断裂能随养护龄期的演化曲线.同样,相场尺度参数b不影响上述定量计算结果.

图8 带方形孔圆板收缩开裂:不同位置处温度、水化度和典型力学参数随时间的演化曲线Fig.8 A circular plate with a square hole:Predicted evolution curves of the temperatures,hydration degrees and mechanical properties

对于b=5 mm 和he=1 mm 情况,图9 中给出了不同时刻的裂缝相场、温度场和水化度场云图.可以看出,早龄期混凝土的裂缝演化过程主要分为以下3个阶段:

图9 带方形孔圆板收缩开裂:不同时刻的裂缝相场、温度场和水化度场云图Fig.9 A circular plate with a square hole:Predicted evolution curves of the temperatures and hydration degrees at various instants

(1)时间段[0 h,16 h]:水化加速和裂缝起裂阶段.该阶段混凝土水化度快速增大,水化反应释放的热量超过周边环境对流和热传导引起的热量损失,导致混凝土温度升高,方孔边缘处混凝土的最高温度可达35.3°C.热膨胀变形引起方形孔顶点处应力集中,裂缝起裂首先出现此处.

(2)时间段[22 h,80 h]:水化减速和裂缝缓慢发展阶段.该阶段水化速率降低,释放的热量小于环境对流和热传导引起的热量损失,导致混凝土温度降低,材料逐渐从热膨胀变形转变为收缩变形.由于水化作用,混凝土获得了70%的破坏强度,足以抵抗热膨胀和自收缩引起的裂缝扩展,该阶段裂缝发展缓慢.

(3) 时间段[80 h,+∞]:水化减速和裂缝迅速扩展并贯穿阶段.该阶段,自收缩应变进一步增大,方形孔顶点位置处的拉应力逐渐增加并超过破坏强度,致使裂缝迅速扩展并最终贯穿整个试件,应力得到释放.最后,混凝土获得了80%左右的破坏强度,其他部位不再有裂缝进一步发展.

可以看出,在不同阶段,混凝土热膨胀变形和自收缩变形两种机制相互竞争,并对早龄期混凝土的应力发展和抗裂性能产生显著影响.

与纯力学载荷作用下的分析结果类似,化-热-力多场耦合相场内聚裂缝模型的数值模拟结果不受裂缝尺度和网格大小影响,能够合理描述早龄期混凝土水化反应和热量传输耦合作用下复杂的开裂过程,特别是热膨胀变形和自收缩变形这两类相互竞争的物理机制得到了准确反映.

3.2 受约束圆环试验收缩开裂试验

考虑如图10 所示的受约束圆环收缩开裂试验[68].

图10 约束混凝土圆环试验:试件几何尺寸和环境条件[68]Fig.10 Restrained concrete ring test:Specimen geometry and environment conditions[68]

混凝土圆环试件的厚度为75 mm,内直径为150 mm,外直径为300 mm.试验中考虑了3.1 mm,9.5 mm 和19.0 mm 3 种厚度的约束钢环(其中,3.1 mm 厚度圆环试件由于刚度较小、约束效应较低,混凝土试件并未开裂,故这里仅考虑后两种厚度情况),并通过贴在钢环中部位置的4 个应变片监测钢环应变.所采用的材料(水灰比w/c=0.3 的砂浆)和模型参数见表2;混凝土水化反应演化函数如图11 所示.

表2 受约束圆环试验收缩开裂试验参数取值Table 2 Model parameters for the cracking analysis of a constrained concrete ring

图11 约束混凝土圆环试验:混凝土水化反应演化函数Fig.11 Restrained concrete ring test:Evolution function of the hydration reaction

类似地,为考察裂缝尺度和网格大小的影响,数值模拟采用两种裂缝尺度b=2.5 mm 和b=5.0 mm,相应的单元网格大小取为he=b/5.

图12 和图13 分别给出了9.5 mm 和19.0 mm 两种厚度试件不同时刻的裂缝相场云图.可以看出,在开始阶段,由于试件的环向应力均匀,微裂缝沿内边缘均匀分布;随着龄期增长,收缩应变增大,由于随机性的影响(这里体现为网格划分引起的数值误差),部分微裂缝较其他裂缝扩展更为明显,但混凝土强度也逐渐提高,导致该阶段总体上裂缝扩展缓慢;随着龄期进一步增长,收缩应变显著增大,导致环向某一位置出现控制性宏观裂缝并迅速扩展至圆环试件外边缘.

图12 约束混凝土圆环试验:约束混凝土圆环试验:两种相场尺度参数条件下不同时刻的裂缝相场云图模拟结果(圆环厚度19 mm)Fig.12 Restrained concrete ring test:Predicted contours of the crack phase-field at various time instants for different length scale parameters(thickness of the steel ring:19 mm)

图13 约束混凝土圆环试验:两种相场尺度参数条件下不同时刻的裂缝相场云图模拟结果(圆环厚度为9.5 mm)Fig.13 Restrained concrete ring test:Predicted contours of the crack phase-field at various time instants for different length scale parameters(thickness of the steel ring:9.5 mm)

数值模拟给出的9.5 mm 和19.0 mm 两种厚度圆钢环应变演化曲线与试验结果的对比如图14 所示.可以看出,考虑化-热-力多场耦合效应的相场内聚裂缝模型能够反映两种不同厚度钢环试件的应变演化规律,特别是能够准确地预测混凝土开裂后、环向拉应力得到释放导致钢环所受压应力迅速降低、压应变出现陡降等典型特征.

图14 约束混凝土圆环试验:圆环应变演化曲线Fig.14 Restrained concrete ring test:predicted evolution curves of the compressive strains of the steel ring at various time instants

同样,与纯力学载荷作用下的结果类似,无论是定性的裂缝演化过程、亦或是定量的圆环压应变演化曲线,裂缝尺度和网格大小对预测结果几乎不产生影响,这一结果对模型应用于实际混凝土结构早龄期裂缝演化过程的准确模拟、抗裂性能的定量预测至关重要.

4 结论与展望

针对早龄期混凝土的抗裂性能预测和裂缝扩展全过程模拟,本工作在统一相场理论的理论框架内,采用相场内聚裂缝模型描述裂缝相场-位移场耦合效应,进一步考虑了温度引起的混凝土热膨胀变形、水化反应导致的混凝土收缩和随龄期变化的力学性能(弹性模量、抗拉强度和断裂能等),以及裂缝对水化过程和传热过程的影响,建立了化-热-力多场耦合相场内聚裂缝模型.

与纯力学载荷作用下的相场内聚裂缝模型一样,该模型同时融入了基于强度的裂缝起裂准则、基于能量的裂缝扩展准则以及基于变分原理的裂缝扩展方向判据,因此能够描述早龄期混凝土水化作用和热量传输引起的复杂裂缝起裂和演化过程.不同的是,刻画裂缝起裂的混凝土抗拉强度并非固定值,而是随龄期增长而逐渐增大并趋于最终强度;类似地,随龄期变化的混凝土弹性模量和断裂能,也影响着早龄期混凝土的力学性能和裂缝演化行为.因此,不同阶段的水化作用和热量传输导致混凝土变形在热膨胀和自收缩之间相互竞争,并对混凝土的早龄期抗裂性能产生显著影响.

上述化-热-力多场耦合相场内聚裂缝模型可以采用有限元等数值方法加以实现,相应单元节点的自由度除位移和裂缝相场外,还包括温度和水化度.考虑到水化过程和传热过程的时间尺度远比开裂过程大,将化-热-力多场耦合作用作用下早龄期混凝土的裂缝演化过程视为化学-热学耦合以及位移-裂缝相场耦合两个串行过程并分别进行求解.前者采用标准的牛顿迭代方法进行求解,而后者则利用笔者近期发展的整体BFGS 拟牛顿算法,实现收敛性和计算效率的双提升.

数值算例表明,所提出的化-热-力多场耦合相场内聚裂缝模型不仅能够描述早龄期混凝土水化反应和热量传输条件下复杂的裂缝演化过程,特别是热膨胀变形引起的裂缝起裂以及自收缩变形导致的裂缝扩展这两类相互竞争的裂缝演化机制,而且数值模拟结果不受裂缝尺度和网格大小的影响,实现了早龄期混凝土抗裂性能的准确预测.

需要指出的是,上述工作忽略了徐变变形(包括基本徐变和瞬态徐变)[20,69]以及水分挥发引起的干缩裂缝[1,70]对混凝土早龄期力学性能的影响,这将在后续研究中进一步加以考虑.