混凝土的徐变是指在持续应力作用下,其应变随时间而持续增长的特性[1]。徐变是混凝土固有的材料特性。一般将混凝土徐变据结构受拉/压不同的受力状态定义为拉伸徐变和压缩徐变;据不同的加载龄期又可分为早龄期徐变和成熟后的徐变[2]。一般认为养护龄期在7 d以内的混凝土为早龄混凝土[2-4]。
已有研究表明,成熟后的混凝土其压缩、拉伸徐变规律大致相同,但早龄混凝土的压缩、拉伸徐变却有较大差别[1-2]。既有的混凝土徐变研究主要针对于成熟混凝土的压缩徐变,早龄混凝土的压缩徐变并未得到关注,其规律并不明确[5]。此外,混凝土在早龄期硬化期间的抗拉性能试验不易实现,故目前国内外相关研究较少。由于混凝土早龄期拉伸、压缩徐变预测模型不足,在对结构进行受力分析时,常以成熟混凝土的压缩徐变代替其早龄期拉伸、压缩徐变。此举将导致结构应力求解精度低, 甚至出现不可信的结果[4-6]。
混凝土材料以其良好的经济性而广泛应用于现代土木工程。但许多混凝土结构在施工养护期间即有不同程度的开裂现象。混凝土早期裂缝对结构的安全性和耐久性造成了严重影响,已成为结构最主要的病害之一[7]。混凝土早龄期拉伸与压缩徐变规律尚不明确,以及结构徐变应力理论计算方法的不完善,是混凝土结构早期开裂无法进行准确的数值模拟、预测并控制的主要原因[3-8]。
早龄混凝土徐变特性比成熟混凝土的压缩徐变更为复杂。其主要特点为:1) 加载龄期早:大体积混凝土在水化期间由于自身水化放热会导致较大的温度应力,该温致应力甚至会致使结构开裂[9];另一方由于其他因素如较大的收缩、支座沉降等不利因素均会导致结构过早受载[10],因此,混凝土在养护期间即较早地参与结构受力;2) 影响因素众多:混凝土配有钢筋,且一般外掺粉煤灰,这些因素均会影响早龄混凝土的徐变[11]。3) 计算理论不完善:一方面早龄混凝土拉伸和压缩徐变发展规律不同,计算时需同时考虑二者不同的徐变效应[5];另一方面,早龄混凝土结构徐变应力计算应考虑不可恢复徐变的影响,需采用非线性徐变理论计算方法[12]。
本文对早龄混凝土的压缩和拉伸徐变研究成果、测试方法及其徐变应力计算方法进行了详细综述。基于既有研究成果,对早龄混凝土的压缩和拉伸徐变测试方法、预测模型以及非线性徐变应力计算提出建议,可为早龄混凝土的徐变研究提供参考。
自Hatt于1907年发现混凝土徐变现象以来,众多学者对这一问题进行了长期研究,形成了诸多徐变理论、徐变预测模型和计算理论,并已反映到相关设计规范中。这些研究成果基本上是针对于成熟混凝土的压缩徐变[13-14]。研究表明,成熟混凝土的压缩徐变与混凝土强度等级、加载龄期、持荷时间、体表比、养护温度、湿度、配筋和粉煤灰掺量等因素相关[15-17]。
目前,混凝土早龄期压缩徐变并未引起足够重视,其关注度甚至不及早龄拉伸徐变,相关研究成果极少。早龄混凝土的压缩徐变与成熟混凝土徐变影响参数基本相同,有可能对某些影响参数更敏感[18]。Li等[19]基于3 d加载的压缩徐变试验发现混凝土早龄期压缩徐变与其强度发展速率存在一定的相关性,而与自身的设计强度等级关系不大。祁广星[20]对加载龄期为1、3、7 d的密封混凝土早龄徐变进行了试验研究,结果表明其压缩徐变较成熟混凝土明显偏大。祝昌暾等[21]研究了掺硅粉混凝土在常温常湿条件下2 d内加载的压缩徐变,基于试验结果对MC-90徐变预测模型进行了修正。谢楠等[22]对喷射混凝土进行了2、5、7、15 d加载龄期的压缩徐变试验研究,基于试验结果修正了GL2000徐变预测模型。试验概况详见文献[22]的相关试验结果如表 1所示,由表 1可知,在徐变发展规律方面,配合比1和配合比2的徐变试件5 d加载持荷30 d时徐变度为2 d加载试件的57.8%和55.4%;7 d加载分别2 d加载试件的46.4%和47.9%;16 d加载分别2 d加载试件的30.3%和28.6%。这说明早龄混凝土的徐变度要远大于成熟混凝土。在徐变预测方面,配合比1和配合比2的徐变试件2 d加载时,其徐变度实测值与GL2000模型预测值误差分别为49.3%和27.4%,误差不可接受;5 d加载时误差值分别为35.6%和13.7%;16d加载时误差值分别为16.2%和3.3%。这说明既有的徐变预测模型不适用于早龄混凝土的压缩徐变预测。
同时,上述两组试验均发现,相对于成熟混凝土压缩徐变而言,混凝土早龄期压缩徐变对加载龄期更敏感。有些参数对早龄期和成熟混凝土压缩徐变的影响甚至呈现相反的规律,这一点可由已有的研究成果可以看出。Alexander[23]研究了粉煤灰对成熟混凝土徐变的影响,发现10 d加载持荷两年后25%掺量的粉煤灰混凝土徐变度下降约33%。而惠荣炎等[24]的研究表明,掺20%粉煤灰大坝混凝土2 d加载,持荷180 d的徐变比未掺试件增大23%。这说明粉煤灰对成熟混凝土的压缩徐变有抑制作用,但增大了早龄期压缩徐变。
目前,成熟混凝土压缩徐变研究已取得了较为丰硕的成果,但混凝土早龄期压缩徐变却有待进一步研究。具体表现在:混凝土强度、粉煤灰掺量和加载龄期对混凝土早龄期压缩徐变影响相关的试验数据不足;其他参数如构件尺寸效应(体表比)、温度和配筋率的影响规律亦少见相关研究;相关试验数据少,预测模型较为缺乏。
在国外,Illston[25]和Domorte[26]较早地报道了混凝土早龄期拉伸徐变并展开了基础性研究。国外相关研究也只是近20年以来才逐渐活跃。Bissonnette等[27-28]的试验研究表明,混凝土拉伸徐变随加载龄期增大而逐渐减小,随环境湿度升高而降低。Nielsen[29]认为应力强度比小于0.6时拉伸徐变与应力呈线性关系, 应力强度比在0.6~0.8之间时,拉伸徐变与应力呈非线性关系。Ya等[30]对掺30%矿渣的混凝土在23℃、33℃和约束作用下早龄期拉伸徐变进行了试验研究,发现掺入矿渣和养护温度升高均加剧了混凝土早龄期拉伸徐变。在混凝土早龄期拉伸徐变预测模型方面,Forth指出既有成熟混凝土的压缩徐变模型不能直接应用于早龄混凝土拉伸徐变预测[31]。Muller等[32]基于试验结果在CEB-FIP 1990模型修正的基础上提出了高性能混凝土早龄期拉伸徐变预测模型。
在中国,惠荣炎等[24]于1986年较早地报道了混凝土的拉伸徐变,对大坝混凝土进行了早龄期压缩(2 d和7 d加载)和拉伸(7 d加载)徐变试验,研究表明,压缩徐变随粉煤灰掺量增加而增大, 而拉伸徐变随粉煤灰掺量增加而减小,粉煤灰掺量对混凝土早龄期拉伸和压缩徐变影响呈现相反的规律。中国相关研究近10年以来才逐渐展开。马新伟[33]对高性能混凝土拉伸徐变特性进行了试验研究,发现混凝土在早龄期尤其是1 d以内其徐变系数要远大于成熟混凝土。汪伦焰等[34]、李司晨[35]等利用自制的拉伸徐变仪较系统的研究了水灰比、外掺物(如粉煤灰、矿渣等)掺量等因素对早龄混凝土拉伸徐变的影响。叶德艳[2]的研究发现,0.5、0.75、1 d加载的混凝土表现出较高的徐变能力,且具备很强的非线性特征。梁思明等进行了早龄混凝土的1 d、3 d和7 d试件拉伸徐变测试,基于B3徐变模型的修正提出了早龄混凝土拉伸徐变模型[36]。杨杨等[37]研究了不同水灰比、加载龄期、加载应力水平和养护温度等参数条件下的早龄期高性能混凝土拉伸徐变特性,并基于流变学原理修正了经典的Burgers模型,提出了评价混凝土早龄期拉伸徐变的ZC模型。以下对Burgers模型(图 1) 和ZC模型(图 2) 进行阐述。
Burgers模型由Voigt体(Kelvin体)与Maxwell体串联而成。混凝土的变形是Maxwell体和Voigt体变形之和。当拉伸荷载σ在t0时刻作用于模型时,将产生弹性应变εe;持荷到t1时刻时,徐变变形在增长,该变形包括Maxwell体中的粘性元件产生粘性流动εd;Voigt体中粘性元件的变形在弹性元件的制约下产生可恢复的延迟弹性变形εv。当荷载t1时刻卸除时,εe立即恢复,随后便是徐变恢复。徐变成分中粘性流动εd在卸载后是不可恢复的。模型与混凝土在加载及卸载过程中所产生变形响应十分类似,但Burgers模型未被赋予明确的物理化学意义。在ZC模型中,引入了一个Hooke弹性元件[H],与Voigt体并联后,再与Maxwell体串联。在混凝土中,水和水泥形成的水泥石用Voigt体表示,水泥石中的结晶体部分是弹簧,凝胶体则是粘壶。当荷载持续不变的情况下,Voigt体开始分担荷载,并随着时间而变形,这时结构内的应力将从Voigt体粘壶向弹簧转移,并发生应力重分布,同时发生的还有整个Kelvin体向单独的弹簧体[H] (代表集料)上的应力重分布[37]。杨杨等利用ZC模型模型预测了早龄混凝土的拉伸徐变,但并未用于早龄结构的徐变分析。
目前,早龄混凝土的拉伸徐变已逐渐成为研究热点,相关研究已取得了一定的成果。拉伸徐变主要研究方法为采用自制仪器在室内进行小试件试验研究,所关注的影响参数基本与压缩徐变相同。由于混凝土早龄期拉伸徐变研究起步较晚,尤其是国内在近10年才展开相关研究,试验数据和预测模型仍较缺乏,亦少见针对早龄混凝土材料力学特性建立的早龄期拉伸徐变预测模型。
由上述文献综述可知,早龄混凝土和成熟混凝土的徐变影响参数基本相同。因此,可借鉴成熟混凝土的徐变研究方法设计试验对混凝土早龄期拉伸、压缩徐变进行研究。成熟混凝土的徐变试验有标准试验方法,试验手段成熟[38]。早龄混凝土的徐变试验尚无规范可循。此外,混凝土早龄期强度较低,拉伸及压缩徐变试验均不易进行,试验的关键在于早龄徐变仪。早龄混凝土的徐变试验装置应具备以下功能特点:为消除试件重力对测试结果的影响,徐变试件应卧置,且应采用有效手段减小试件和接触面间的摩阻力;能实现不同的应力水平加载,并具备良好的持荷能力;可自动、连续地采集试验数据。
目前,混凝土徐变试验采用的试验仪器主要有以下5种:1) 弹簧立拉式徐变装置[38];2) 杠杆立拉式徐变装置[38];3) 杠杆卧拉式徐变装置[39];4) 液压卧拉式徐变装置[34, 37];5) 清华大学土木工程系建筑材料研究所研发的TSTM装置[40]。上述各种仪器各有特点,其中,第1、2种装置只能用于晚龄期拉伸徐变试验。第3~5种装置可用于混凝土早龄徐变试验。
目前,早龄混凝土徐变测试方法较为缺乏,主要体现在:无规范可循、无标准试件尺寸、自制实验仪器昂贵、无法进行批量试验。
混凝土结构在其施工期间的早期开裂现象较为普遍,这些裂缝对结构安全性和耐久性产生极为不利的影响。而早龄期结构开裂的有效预测并采取针对性控制措施依赖于对其应力的准确预测。随着计算机的普及和有限元法的广泛应用,复杂结构的徐变应力求解成为可能。由于混凝土早龄期徐变预测模型的缺乏,以成熟混凝土的压缩徐变代替早龄期徐变进行混凝土结构受力分析,其误差较大[3-5, 41]。Kwak等[42]和Cusson等[43]建议采用拉伸徐变模型对结构早龄期非荷载裂缝进行数值模拟分析。苏安双[44]采用混凝土拉伸徐变模型对高性能混凝土早期收缩应力进行预测,获得了较好的结果。巨玉文等[40]分别采用混凝土早龄期拉伸和压缩徐变预测模型对地下连续墙养护期间的水化热温致应力进行了对比分析,计算结果表明采用拉伸与压缩徐变预测模型其结果相差较大,在早龄期受拉应力状态下使用压缩徐变规律会得到偏于不安全的计算结果。
除了能反映混凝土材料性能的早龄期徐变预测模型外,适用的徐变应力理论分析方法亦尤为重要。目前,早龄混凝土结构徐变应力分析一般采用增量初应变法,即将混凝土徐变历程分为若干时步,每一计算时步采用上一时步的徐变应变增量计算徐变应力增量。Zienkiewicz等[45]提出了等时步条件下徐变增量的递推算法,对增量初应变法进行了简化。朱伯芳院士[46]又提出了混凝土结构徐变应力分析的隐式解法。秦煜等[47]采用初应变增量有限元法建立了混凝土箱梁水化热温致应力的弹性徐变隐式解法数值模型,据此得到的仿真分析结果与实测结果吻合较好。大型计算机的应用可方便的实现上述算法。
而徐变应力分析结果决定性因素为徐变计算理论。经典的徐变计算理论主要有老化理论、有效模量法、弹性老化理论、弹性徐变理论及继效流动理论等。上述算法及理论均假定混凝土徐变与应力呈线性关系,并服从叠加原理[48]。线性徐变理论未考虑混凝土的不可恢复徐变在不同应力水平下的非线性性质[49]。如大体积混凝土在水化时经历先升温后降温的过程, 对徐变应力的求解需要考虑内力的加载和卸载时程, 应采用相应的非线性徐变理论计算。由Bazant等[50-51]提出的固化徐变理论可实现对早龄期结构的非线性徐变应力求解。混凝土固化徐变理论将弹性理论、粘弹性理论和流变理论结合起来,模拟混凝土宏观物理力学性质。Bazart模型有非常清晰的物理意义,认为混凝土的基本徐变由老化粘弹性项徐变Ca(t, τ)、非老化粘弹性项徐变Cna(t, τ)、粘性流动项徐变(不可复徐变)Cf(t, τ)组成[48],即
采用混凝土固化徐变理论建立混凝土徐变度函数,可据此建立混凝土非线性徐变应力求解的理论模型。张涛等[12]在此方面做了一些研究工作,采用固化徐变理论对混凝土拱坝水化热温致应力场进行非线性徐变应力分析,取得了较好的效果。
综上所述:混凝土早龄期徐变理论计算方法是实现对徐变应力准确求解的关键。目前国内外相关研究已取得了一定的成果。但对于早龄混凝土结构,在低应力水平条件下其徐变即表现出非线性性质。如何在早龄期徐变应力理论分析中,既反映其早龄期在受拉和受压时两种不同的松弛特性,又考虑其非线性特征,是有待深入研究的问题。
1) 混凝土早龄期拉伸、压缩徐变测试方法无规范可循,试验数据较为缺乏。一方面早龄混凝土徐变测试仪较少,早龄期拉伸、压缩徐变试验研究较少,已有研究主要探讨了强度等级、加载龄期、应力水平等参数对混凝土早龄期徐变的影响,构件的尺寸效应(体表比)、粉煤灰掺量和配筋率等参数对混凝土早龄期徐变影响规律尚未明确。另一方面,已有研究多数只是单独进行拉伸或压缩徐变试验,混凝土早龄期拉伸、压缩徐变与成熟混凝土压缩徐变在数值上存在多大差别,发展规律有何不同,缺乏定量的报道,亦未引起足够重视。
2) 混凝土早龄期拉伸、压缩徐变预测模型不足。目前混凝土早龄期徐变预测模型还少有比较明确的规定。少数几个修正模型基本上源于成熟混凝土压缩徐变模型的修正,未能反映混凝土早龄期物理力学性能的发展和徐变的不可恢复性。因此现有徐变规律并不能准确反映混凝土早龄期徐变性能。
3) 早龄混凝土结构徐变应力理论分析方法不尽完善。早龄混凝土结构在受拉和受压时其应力具有不同松弛特性。此外,早龄混凝土结构的徐变应力理论计算方法应考虑其非线性特征。基于固化徐变理论,建立早龄混凝土结构同时考虑受拉和受压不同应力松弛特性的非线性徐变应力理论计算方法,应可满足该要求。