随着高性能混凝土在土木工程中的广泛应用,其早龄期开裂现象及其力学性能日益受到重视而成为研究热点,混凝土的早龄期蠕变性能与混凝土的力学性能有着很强的相关性,是预测与减缓混凝土早期开裂时间的重要依据[1]。关于混凝土早期蠕变性能的研究主要分为理论模型和和经验模型两类。理论模型包括有效模量法模型[2]、固化蠕变理论模型[3]、弹性蠕变老化理论模型以及微预应力—固化理论模型等[4, 5, 6]。使用较多的理论模型为Bazant于1995年提出的B3模型[7],该模型是以固结理论为基础,在大量总结徐变试验数据的基础上提出的半经验半理论模型。较其他徐变模型而言,B3模型参数较多,表达式具有较为明确的物理意义。试验拟合模型主要有CEB-FIP(MC90)模型[8]、Muller模型[9]和ACI209系列模型[10]等。CEB-FIP(MC90)模型由欧洲混凝土委员会和国际预应力联合会(Euro-International concrete committee and International Federation for Prestressing)于1990年提出,该模型参数较少,适用于早龄期与成熟混凝土,对于高强度混凝土的徐变特性模拟精度稍差。在此基础上所提出的Muller模型,其形式与CEB-FIP(MC90)模型类似,但增加了3个与混凝土强度有关的修正系数,使其能更好地适于高强混凝土的徐变模拟。ACI209系列模型是由美国混凝土协会基于试验数据提出的经验模型,参数较少,但不适于加载时间早于7 d的混凝土。上述3种经验模型,受限于实验条件以及测量、养护环境的部分不可控因素,且不同模型所考虑的参数类型与数量不尽相同,使得不同模型在针对某一具体实验或混凝土材料时,产生较大的徐变特性预测差异[11]。此外,由于经验公式均由试验数据拟合而得,对于混凝土早龄期蠕变特性的机理解释不尽明确。
基于早龄期混凝土的粘弹性性能和流变性能,采用弹簧和阻尼器串并联的方式,用以模拟早龄期混凝土的粘弹性以及流变性能的元件模型也被用于早龄期混凝土蠕变机理的模拟。 Farah等[12]、Bazant等[13]、Sellier等[14]提出了各自的粘弹性元件模型,其模拟结果与各自实验数据吻合度较高,但这类模型普遍存在使用的弹簧和阻尼器较多,各元件物理意义不明确,,不易确定的不足,且模型的普遍适用性也有待进一步的验证。
对此,笔者提出了变系数四参数Burgers模型[15],该模型为基于早龄期水泥基材料的蠕变松弛机理及其物理化学性能变化特性的元件模型,模型公式简洁,参数较少并均具有明确的物理意义。
本文收集了多组早龄期混凝土徐变特性实验数据[16, 17, 18, 19, 20],对各类模型的蠕变预测效果及其误差原因进行了对比分析,并进一步验证了变系数四参数Burgers模型的适用性以及参数的合理取值范围。
CEB-FIP(MC90)模型为经验模型,主要适用于28 d圆柱体抗压强度为12~80 MPa、养护平均湿度为40%~100%、养护平均温度为5 ℃~30 ℃的早龄期或成熟混凝土:
其中: Ec(t0)为混凝土加载时的弹性模量,如式(2)所示,Ф(t,t0)为徐变系数,如式(3)、(4)所示。
式中: Ac/μ代表混凝土试件等效截面高度。由上述公式可见,CEB-FIP(MC90)模型主要考虑了28 d混凝土强度fcm和弹性模量Ec、养护湿度h、初始加载时间t0与加载持续时间t、水泥品种有关的参数S以及构件尺寸效应Ac/μ等因素的影响。
Muller模型与CEB-FIP(MC90)模型相似,也为经验模型,其适用范围和考虑因素与CEB-FIP(MC90)模型相同,但增加了3个与混凝土强度有关的修正参数αm1、αm2和αm3,改善了对高强混凝土徐变特性模拟的准确性。Muller蠕变模型具体如下
其中,
与混凝土强度有关的修正参数αm1、αm2和αm3计算方式如下所示:αm1=(35/fcm)0.7;αm2=(35/fcm)0.2;αm3=(35/fcm)0.5。
B3模型是Bazant于1995年在BP模型[21]和BP-KX模型[22]基础上进行完善而提出的,它以固结理论为理论基础,认为水化作用引起的混凝土内部微观颗粒物质向水泥胶体的转变是徐变产生的内在原因,并主张将混凝土总应变分为弹性应变、粘弹性应变、流动应变和附加应变(温度、干燥等因素引起的应变)[3]。由此,如式(3)所示,混凝土徐变便也成为基本徐变C0(t,t0)和干燥徐变Cd(t,t0,tc)之和。该模型适用于28 d圆柱体抗压强度为17~70 MPa、水灰比为0.3~0.85、骨料水泥比为2.5~13.5、水泥含量为160~720 kg/cm3的普通混凝土,具体如下
其中:式(9)中第一项为混凝土加载时的瞬时应变率,第二项C0(t,t0)为基本徐变,第三项Cd(t,t0,tc)为与混凝土干燥过程有关的干燥徐变。基本徐变C0(t,t0)包含3部分,如公式(10)所示。第一项为水化后水泥胶体的粘弹性徐变,主要受水泥用量c的影响;第二项为未水化颗粒的粘弹性徐变,主要受水灰比w/c和水泥用量c的影响;第三项为流变,主要受骨料水泥比a/c的影响。
其中: Q(t,t0)为经验拟合公式,如式(11)、(12)所示。m,n为常数,m=0.5,n=0.1。
干燥徐变Cd(t,t0,tc)主要考虑混凝土干燥过程引起的徐变,如式(13)~(17)所示
B3模型中,基本徐变C0(t,t0)主要受水泥用量c、水灰比w/c、骨料水泥比a/c的影响。而干燥徐变Cd(t,t0,tc)主要与混凝土尺寸V/S、形状ks、养护条件α2、混凝土开始干燥时间tc、水泥品种α1、环境湿度h、混凝土强度fcm等有关。总体而言,B3模型所需参数较多,考虑因素较全面,模型较复杂。
如图 1所示,由笔者在文[15]、[16]中所建立的变系数四参数Burgers模型由Maxell模型和Kelvin模型组成。其中,Maxell体对应于水化过程中材料物质形态由流动态向固态的转变过程时其粘性和弹性的变化特性,表征了混凝土徐变中不可恢复的部分;Kelvin体则描述了材料组分中未水化的胶凝颗粒、骨料等的粘弹性性能,表征了混凝土徐变中可部分恢复的部分。
变系数四参数Burgers模型的蠕变模型如下所示:
式(4)中各参数右上标“B”表示Burgers模型,t0为加载时间。其中,弹簧E1B(t)代表早期混合物中骨料及未水化的胶凝颗粒的弹性,随水化进程而变化。如式(5),Ec为混凝土中颗粒材料的初始弹性,Sc(t)为其时间变化函数,ac为对应的与时间相关的参数。
弹簧E2B(t)表示水化固体产物的瞬时弹性,如式(6)。其中,E∞为混凝土固化后的长期弹性模量,一般可取做28 d龄期弹性模量E28,SE(t)为弹性模量随时间的变化函数,aE为对应的与时间相关的参数,可由试验数据拟合得到:
粘壶η1B(t)表示水化胶凝产物的粘性,如式(7)。其中,η1∞为材料长期粘性系数,通常由水泥砂浆材料性质、环境温湿度条件等决定,Sη1(t)为其时间变化函数,aη1为对应的与时间相关的参数
粘壶η2B(t)代表由于胶凝水等的迁移引起的蠕变,具有不可恢复性,如式(8)。其中,η2∞表示长期不可恢复蠕变,Sη2(t)为其时间变化函数,aη2为对应的与时间相关的参数:
为对比各类模型的蠕变预测效果及其误差原因,并进一步验证变系数四参数Burgers模型的适用性以及参数的合理取值范围,本文收集了多组早龄期混凝土徐变特性实验数据[17, 18, 19, 20, 21],其相关实验参数见表 1。
图 2为针对熊维的实验数据采用前述4种蠕变模型的预测结果,其中加载龄期t0=3 d。从图中可以看出,与实验数据相比,B3模型的预测误差最大,CEB-FIP模型和Muller模型的预测结果与试验结果均较吻合,且Muller模型的精度更高,这与Muller模型针对高强混凝土进行了参数修正有关,变系数四参数Burgers模型的预测结果吻合度也较高。
图 3为针对Philippe的实验数据进行的4种蠕变模型的预测结果,其中加载日期t0=7 d。同样可见,CEB-FIP模型和变系数四参数Burgers模型的预测精度均较高,B3模型则高估了较早期混凝土的蠕变,在28 d龄期后与试验数据逐渐吻合。而Muller模型的预测值较原始数据而言偏低,表明其对普通混凝土的适用性较差。
图 4为针对Atrushi的实验数据,采用4种不同模型的蠕变预测对比,其中,加载龄期t0分别为2 d、3 d、6 d,而变系数四参数Burgers模型采用加载龄期为3 d时的数据进行3种不同加载龄期的蠕变预测。由图 4可见,变系数四参数Burgers模型预测曲线与原始数据的吻合度最高,其次分别为CEB-FIP模型和B3模型,而Muller模型偏差最大,且随着初始加载龄期的增长,后3种模型的预测值均逐渐小于实验值,即低估了混凝土的早期蠕变。变系数四参数Burgers模型预测曲线在t0分别为2 d、3 d和6 d时,与原始数据都有很高的贴合度。
图 5为针对Wibke的实验数据,4个不同模型的蠕变函数预测曲线的对比结果。其中,初始加载龄期t0分别为3 d、4 d、7 d,变系数四参数Burgers模型采用加载龄期为7 d时的数据进行模拟。值得注意的是,Wibke采用的是混凝土受拉徐变试验。从图中可以看出,各模型在初始加载龄期较小时(t0=3 d、4 d),其预测曲线与试验数据均偏差较大,随着加载时间增长,如t0=7 d,变系数四参数Burgers模型的预测曲线、CEB-FIP模型以及Muller模型的预测曲线与试验数据逐渐吻合,但B3模型的偏差仍然较大,且整体高估了混凝土的蠕变。
图 6为针对Zheng 等的实验数据,4个不同模型的蠕变函数预测曲线的对比结果,初始加载龄期3 d、7 d、14 d,变系数四参数Burgers模型采用加载日期为7 d时的数据进行参数模拟。从图中可以看出,变系数四参数Burgers模型的预测曲线与原始数据在加载龄期为3 d、7 d时吻合度很高;当加载龄期为14 d时,预测曲线较原始数据偏高。而其他3个模型较原始数据均有不同程度的高估。
综合上述对比,可以看出:CEB-FIP模型和Muller模型作为经验公式,对于不同实验数据,其蠕变预测曲线的准确性有较大差异。这是仅仅依据试验数据的拟合而建立的经验公式的通病,即容易受到实验方法、实验条件或者测量等不可控因素的局限。此外,拟合结果对比还显示,针对强度较高的混凝土(熊维和Wibkle),Muller模型的蠕变预测精度较高。
B3模型的蠕变预测曲线普遍较原始数据偏大。这是由于B3模型以固化理论为基础,对于早龄期混凝土,由于水化作用发展迅速,材料的物理性态变化迅猛,B3模型对此描述不足;其次,B3模型中经验拟合数据均来源于RILEM收缩徐变数据库,该数据库中加载龄期小于7 d的数据非常有限,导致该模型对加载龄期早于7 d的早龄期混凝土蠕变性能模拟较差。
笔者所建立的变系数四参数Burgers模型同样基于固化理论,但结合了早期混凝土的蠕变机制与宏观变化特性,而避免了从细观结构的描述出发导致的模型复杂、参数较多且不易确定的弊端。变系数四参数Burgers模型的参数相对较少,只保留了B3模型中较核心的影响因素,即水化前颗粒材料的粘性和弹性,以及水化后水泥胶体材料的粘弹性和流变性,与各组试验数据较好的吻合度验证了该模型的适用性与合理性。
变系数四参数Burgers模型将早龄期混凝土的蠕变视为胶凝材料的蠕变与固体颗粒材料(包括未水化的水泥、骨料、砂等)的蠕变之和,并分别描述了两个部分弹性和粘性随时间的变化。整体包括四个时变函数,共计8个参数。其中,E28和ac可由实验直接测得,其余6个参数需拟合得到。将针对表 1所列实验数据,分析四参数Burgers模型参数的取值范围。。
将不同试验拟合得到的四参数Burgers模型参数列于表 2。由于Wibke实验数据是针对早龄期混凝土的受拉徐变的,参数拟合结果与其他实验数据的拟合结果相差较大。 熊维、Philippe、Atrushi、Zheng等试验数据的参数拟合结果,其特点如下:
变系数四参数Burgers模型中,Ec表示混凝土中颗粒材料的弹性性能。4组实验数据使用的均为CEMΙ型水泥,骨料与水泥的用料也基本相同,Ec拟合结果的一致性显示,模型中假定混凝土中未水化颗粒材料的弹性模量是合理的。表 2还显示,针对Philippe的试验数据得到的Ec值偏小。这是由于Philippe的实验所采用的混凝土28 d的强度为42 MPa,较其他组实验数据低,据此可以推测,Philippe实验中固体颗粒材料的强度与弹性模量相对较低。此外,Philippe实验采用100%的湿度养护方式,养护时间内颗粒材料的消耗较快,使得在进行加载实验时的实验试件中未水化的颗粒物质相对较少,从而可能引起弹性模量的降低。上述两个原因可能共同导致了Philippe的实验参数Ec值偏低。
参数η1∞表示水化胶凝材料的长期粘性。4组实验数据的拟合结果基本相同,表明模型中引入该项参数的合理性。值得注意的是,根据Atrushi试验数据确定的η1∞值较小,事实上,Atrushi的实验试件尺寸为150 mm×300 mm的圆柱体,其等效截面高度较其他实验试件而言偏大。几何尺寸差异可能是导致η1∞值与其他实验数据有所差异的原因。
参数η2∞表示由于胶凝水等的迁移引起的不可恢复蠕变。由于4组实验数据所使用的水泥品种、骨料与水泥用量、加载时环境湿度和温度都比较一致,因此4组数据的η2∞值也应该相差无几,这从拟合结果中得到了印证。同样可以注意到,由于试件尺寸效应的影响,针对Atrushi实验数据的η2∞值较其他试验数据小。
aE、aη1和aη2这3个与时间变化相关的参数,4组实验数据的拟合结果比较接近。这表明在基本相同的温湿度试验条件下,早龄期混凝土的粘弹性变化规律应该比较接近。
从上述的参数分析结果可以看出,所建立的变系数四参数Burgers模型各参数所对应的物理意义是合理的,且各参数存在一定的取值范围。结合四参数Burgers模型参数拟合结果,可以得出以下的经验取值范围:表征早期混合物中骨料及未水化的胶凝颗粒弹性的弹性模量Ec参考取值为6×104~9×104 MPa,并且应该考虑混凝土28 d强度以及养护方式对Ec的影响;表征水化胶凝产物粘性的参数η1∞参考取值为20×104~30×104 MPa·s,且应考虑混凝土试件尺寸效应;表征水泥基材料早期液态流动性的参数η2∞的参考取值为300×104~600×104MPa·s。表征水化胶凝产物瞬时弹性的参数E2B(t)则可由实验直接测得。
以上模型对比研究表明,针对不同实验数据以及加载龄期,变系数四参数Burgers模型对于早龄期混凝土的蠕变预测准确度较高。对于大部分实验数据,变系数四参数Burgers模型与目前通行的徐变模型也较接近,从而验证了变系数四参数Burgers模型针对于早龄期混凝土的徐变预测的适用性与合理性。
变系数四参数Burgers模型中的8个参数除E28、ac两个参数能够直接从实验室中测量得到外,其余6个参数Ec、η1∞、η2∞、aE、aη1、aη2均需要根据试验数据进行拟合。通过对于不同原始试验数据拟合结果的整理和分析,可以看出当水泥品种、用量以及骨料品种、用量类似的混凝土,其Ec、η1∞、η2∞等参数拟合值较为接近,这表明变系数四参数Burgers模型对于参数的假定及所赋予的物理意义是合理的。当试件养护条件(即环境温度和湿度)相差不大时,Ec、η1∞、η2∞的时间变化规律也较为相似。这说明所建立的四参数Burgers模型反应了早龄期混凝土的粘弹性随时间的变化规律主要受环境温度和湿度影响的机理。
以上变系数四参数Burgers模型的参数分析及其与常用经验徐变模型和B3模型的对比表明,变系数四参数Burgers模型中对混凝土强度、养护方法、尺寸效应、环境温湿度的具体作用机制等因素还缺乏明确反应,仅隐含于参数的具体拟合值中,因此,后续研究将通过对早龄期混凝土水化进程中温度湿度耦合效应的研究,进一步完善Burgers模型。