李庆华，种法澄，张 麒，徐世烺

（浙江大学 高性能建筑结构与材料研究所，浙江 杭州 310058）

1 研究背景

目前我国有众多已建或在建、待建的混凝土水坝工程，工程规模位居世界前列［1］，且多建在西南部多地震区域［2］。由于混凝土自身抗拉性能低而且大体积混凝土在浇筑过程中温度不易控制、振捣质量难以保证，使得混凝土内部容易产生微观裂缝，这些微观裂缝的存在使得结构的宏观承载力下降，存在发生低应力破坏的风险，对结构的安全运行埋下了隐患。因此，采用合适的手段对混凝土裂缝稳定性和安全性进行评估，及时采取加固措施，对于保障结构的安全稳定运行具有重大的意义。

众多混凝土断裂力学试验均表明，由于骨料的咬合桥联作用，混凝土的断裂过程可以分为裂缝起裂、裂缝稳定扩展、裂缝失稳扩展三个阶段。对于抗渗性能要求较高的混凝土水坝的安全性评估，确定混凝土内部裂缝是否起裂十分重要。为此，徐世烺等［3-5］引入起裂韧度作为裂缝起裂时的断裂韧度，并将传统的断裂韧度记为失稳韧度提出了两个应力强度因子表征的双K断裂模型。能量法和应力强度因子法在线弹性断裂力学中是等效的，赵艳华等［6］在双K断裂模型的基础上，提出了能量法表征的双G断裂模型。双K断裂模型和双G断裂模型均可以很好地描述混凝土断裂的全过程，且互为验证，为混凝土断裂性能的评估提供了理论依据。双K断裂模型已多次用作大坝混凝土安全性评估的基本理论，吴瑶等［2］曾基于双K断裂模型对丹江口大坝混凝土的断裂性能进行评估检测，认为上游竖向裂缝尖端应力强度因子已超过起裂韧度需要进行排水加固。

对钻芯样品进行断裂力学试验是大坝混凝土断裂性能评估的常用方法。圆形紧凑拉伸法适用于钻芯样品，因此不断有学者希望将这种常用于沥青材料［7-8］的断裂力学试验方法引入混凝土。2005年，Wagoner等［9］将这种方法应用到混凝土试件上，但是在试验过程中50%左右的试件在预制加载孔处的混凝土发生了提前破坏。为了解决这一问题，Cifuentes等［10］提出了一种改进的圆形紧凑拉伸法，改进了传统方法的加载方式，从而避免了加载孔处混凝土的提前破坏。但是尚未有学者采用该方法测定混凝土双K断裂参数和双G断裂参数，关于该方法能否获得稳定的两种断裂参数有待试验研究。本文采用改进的圆形紧凑拉伸法，用于测定不同尺寸下试件的双G参数和双K参数。一方面通过分析每组数据的离散性，验证这种新型的混凝土断裂力学试验方法能否获得稳定的双K与双G断裂参数，另一方面通过比较由公式计算间接得到的等效断裂韧度和双K模型计算直接得到的断裂韧度，研究双G断裂模型和双K断裂模型在描述断裂性能方面的等效性，最后讨论两种断裂模型参数的尺寸效应。

2 断裂参数的计算方法

2.1 双K断裂参数的计算 由于混凝土断裂过程区的存在，经典线弹性断裂力学公式不能直接运用到混凝土，为此引入线性渐进叠加假定［11］，假定的内容为：（1）P-CMOD曲线的非线性特征是由自由裂缝面前端的虚拟裂缝引起；（2）有效裂缝包括等效弹性自由裂缝和等效弹性虚拟裂缝两部分。

下面对该假定进行简单的推导：图1为某试件的P-CMOD曲线，该试件预制裂缝长度为a0。A点为线弹性点，在A点之前裂缝不起裂，裂缝长度保持为a0不变。将A点的荷载Pini和裂缝长度a0带入线弹性断裂力学公式即可求其起裂韧度。B点为非线性点，其有效裂缝长度用ab表示，则虚拟裂缝长度∆ab=ab-a0。若试件在B点卸载，并忽略残余变形，假设卸载至原点。那么B点就可以看作是预制裂缝为ab试件的线弹性点，线弹性断裂力学公式仍可用。以此类推，非线性的P-CMOD曲线可以看作是一系列线弹性点的组合，因此一个完整的考虑非线性特征的断裂过程可以采用线弹性断裂方法加以描述。

图1 线性渐进叠加假设

确定起裂荷载Pini是双K断裂模型的关键，常用的测量起裂荷载的方法有激光散斑法、光弹贴片法、扫描电子显微镜法和电阻应变片法等，本文采用电阻应变片法。具体方法为：以初始裂缝为中轴线，在试件两侧平面距预制缝处1 cm的位置对称粘贴一对应变片，应变片与另外两对温度补偿片构成两个全桥回路，用以测定起裂荷载。在裂缝起裂之前，随着荷载的增加，预制裂缝两侧的混凝土不断积聚能量，使得缝端应变片的应变值不断增大，并且基本呈线性增长，直至在某一荷载下发生回缩，这表明观测点间有裂缝出现，混凝土积聚的能量得到释放，导致应变减小。此时应变回缩点所对应的荷载即为起裂荷载Pini。

由于裂缝起裂之前，混凝土可以近似看作是线弹性材料。因此得到了起裂荷载Pini之后，将其和初始裂缝长度a0一起带入线弹性应力强度因子计算公式便可得到起裂韧度的大小，对于本文试件形式，具体的公式为［12］：

其中：

式中：T为试件厚度；W为试件净高；a0为初始裂缝长度；α为缝高比，α=a0/W；Pini为起裂荷载；a为裂缝长度。

直接进行混凝土弹性模量试验测得的弹性模量并不能准确反映每一个试件的弹性模量，每个试件的弹性模量都会稍有不同。因此本文试验采用荷载-裂缝张开口位移公式反算试件弹性模量。对于本文试件形式，荷载-裂缝张开口位移的关系式如下［12］：

其中：

式中：E为混凝土弹性模量。

然后将试验测得的起裂荷载Pini、起裂时裂缝张开口位移CMODini、初始裂缝长度a0、试件厚度T带入到式（3）即可得到该试件的弹性模量E。

双K断裂模型在吸取了两参数模型和等效裂缝模型等经典断裂模型优点的基础上，提出了线性渐进叠加假设。在该假设的前提下，把试验测得的最大荷载Pmax、此时对应的裂缝张开口位移CMODc代替起裂荷载Pini起裂时裂缝张开口位移CMODini和计算得到的试件弹性模量E一并代入到式（3），通过数值计算即可获得临界等效裂缝长度ac。

在线性渐进叠加假设的前提下，将最大荷载Pmax和临界等效裂缝长度ac，代替起裂荷载Pini和初始裂缝长度a0代入式（1），即可得到失稳韧度

根据线性渐进叠加假定，非线性的P-CMOD曲线上可以看作是一系列线弹性点的组合，因此一个完整的考虑非线性特征的断裂过程可以采用线弹性断裂方法加以描述，所以经典的线弹性断裂力学公式是可用的。

2.2 双G断裂参数的计算 根据线弹性断裂力学，对于线弹性材料，其能量释放率为［12］：

式中：C为试件的柔度，对于圆形紧凑拉伸试件而言，C=CMOD/P，P为荷载，CMOD为裂缝张开口位移。

对于圆形紧凑拉伸试件而言，T、W、α均为已知参数，计算G首先要求解dC/dα。

对于圆形紧凑拉伸试件，荷载-裂缝张开口位移的关系式如式（3）、式（4）所示［12］，由C=CMOD/P即：

对式（6）两边同时求导，可得：

将式（7）代入式（5）得基本计算公式。

现将双G断裂模型的计算方法总结概括如下：（1）根据式（3）计算混凝土的弹性模量E；（2）将计算得到的弹性模量E、峰值荷载Pmax、临界裂缝张开口位移CMODc代入式（3），通过数值计算的方式求解临界等效裂缝长度ac；（3）将起裂荷载Pini、初始裂缝长度a0代入式（5）便可得到起裂能量释放率（4）将峰值荷载Pmax、临界等效裂缝长度ac代入式（5）便可得到失稳能量释放率

3 试验过程

3.1 试件制作 浇筑模具外壁为标准尺寸的PVC管材，管材的内径分别为150、200、250和400 mm。管材下部用木模板封口。用黑色马克笔在模具外壁按标准试件尺寸做好相应标记，用电钻在模板外壁标记处钻孔，插入长度50 cm、直径10 mm的螺纹钢筋，并用热熔胶做好钢筋洞口处的密封。同一尺寸试件两两一起浇筑，后期切割成两个试件。

本文试验材料采用标号为42.5的普通硅酸盐水泥、河沙和直径15 mm以下的碎石，各试件统一按照配合比为：水泥∶石子∶沙子∶水=1∶2.3∶1.19∶0.41的混凝土进行浇筑。按尺寸不同共浇筑4组，每组4个试件，每组试件具体尺寸见表1。同时浇筑3个尺寸为150 mm×150 mm×150 mm的立方体抗压试件，测得28 d抗压强度的平均值为32.07 MPa。

试件在温度为（20±3）℃和相对湿度在90%以上的标准养护室养护28 d后，使用大型精密切割仪器按照图2所示尺寸切割加工。其中初始裂缝由切割机切割制作，初始裂缝长度如表1所示。预埋钢筋随裂缝的切割被一起切断。

3.2 测试及加载过程 以初始裂缝为中轴线，在试件两侧面距离预制缝处约1 cm的位置分别对称粘贴一对应变片，应变片与另外两对温度补偿片构成两个全桥回路，用以测定起裂荷载。除此之外，在断裂过程区两侧依次布置两排应变片，用以监测裂缝的扩展。试验在250 kN的Instron试验机上开展，使用合适的夹头夹住试件的两端钢筋。使用夹式引伸仪测量裂缝张开口处位移CMOD，并以裂缝张开口位移控制加载速率。最小尺寸试件的初始加载速率为2.5 μm/min，荷载达到峰值荷载后可逐步提速到25 μm/min。待CMOD数值增长到4 mm（夹式引伸仪量程）或试件破坏时加载结束。试验基本测得量包括荷载P、裂缝开口位移CMOD以及裂缝尖端处的应变值。

表1 试件尺寸

图2 改进的圆形紧凑拉伸试件的标准几何尺寸

4 结果及分析讨论

4.1 起裂荷载Pini的确定 试件中典型的荷载-应变曲线如图3所示。从图3中可以看出，在开始加载时，应变片应变值基本呈线性增长，当应变达到极值在某一点发生回缩，说明测点间有裂缝产生，混凝土积累的能量得到释放，因此应变片开始回缩时对应的荷载即为起裂荷载Pini。

4.2 峰值荷载Pmax的确定 试件中典型的荷载-裂缝张开口位移曲线如图4所示。从图4可以看出，试件加载初期P-CMOD曲线基本呈线性增长，当裂缝起裂时，曲线呈现非线性增长，达到峰值后平缓下降。图4曲线全程平滑，没有出现峰后荷载陡降，尾部曲线非常平稳，说明试验过程十分稳定。

4.3 试验计算结果 试验按试件直径变化共设计4组16个试件，最终得到15个试件的实测数据，按上述方法计算得到起裂能量释放率和失稳能量释放率由上文提到的线性渐进叠加假设可知，线弹性断裂力学公式是可用的［13］。为了验证双G断裂模型和双K断裂模型在描述断裂性能方面是否具有等效性，首先利用上式将能量释放率转化成相应的应力强度因子表示，记为和同时依据双K断裂模型计算得到双K断裂韧度，记为由于双G参数与双K参数的量纲不同，不可以利用方差直接比较数据离散性的大小，因此计算了每组数据的变异系数，用以比较双G参数与双K参数离散性的差异。将每个试件的计算信息汇总于表2，其中由于试件切割存在误差，初始裂缝长度计算值以实际切割长度为准。

图3 典型试件的荷载-应变曲线（D150-4）

图4 典型试件的荷载-裂缝张开口位移曲线（D300-3）

表2 试验数据及参数计算

对比表2和表3中经修正前后试件的双K和双G断裂参数数据可以发现，修正前后起裂韧度的误差一般在30%以内，最小误差7.4%，最大误差48.6%；修正前后失稳韧度的误差一般在50%以内，最小误差7.4%，最大误差56.4%。这说明混凝土前期的损伤对失稳韧度的影响更大。类似的，起裂能量释放率的误差一般在10%以内，最大误差14.0%；失稳能量释放率的误差一般在15%以内，最大误差19.5%。从上述的对比可以看出，一方面，试件损伤对于失稳断裂参数和的影响更大，另一方面，与双G断裂参数相比，试件损伤对于双K参数的影响更大。同时较多试件发生提前损伤的原因可能是试件在制作的过程中经历过多次的切割，因此如何在试件制作过程中尽量减小损伤，是改进的圆形紧凑拉伸法亟待解决的问题。

表3 修正后的计算结果

4.4 分析讨论 由表3中的数据可以看出，试验获得的数据较为稳定，离散性较小。这说明改进的圆形紧凑拉伸法可行，能够获得稳定可靠的试验数据。其中断裂应力强度因子的变异系数一般在0.3以下，断裂能量释放率的变异系数一般在0.5以下。与双K参数相比，双G参数的离散性稍大，这可能是由于能量释放率是应力强度因子和弹性模量两个参数的函数，任何一个参数的变化都会对能量释放率的数值产生的影响。

此外，从图5还可以看出，当试件直径从150 mm增长到200 mm时，起裂韧度有较为明显增长，失稳韧度略有下降，但幅度并不明显。当试件直径从150 mm增大到200 mm时，起裂韧度的平均值增长了43.1%、的平均值增长了42.9%；失稳韧度则下降了4.6%、下降了4.3%。当试件直径在200 mm以上时，在本文试验的范围内，起裂韧度没有明显的尺寸效应，但是失稳韧度随着尺寸的增长而增长，但增长幅度并不明显。以每组平均值作为参考，当试件直径从200 mm增长到250 mm时，失稳韧度增长了18.0%，增长了17.8%；当试件直径从250 mm增长到300 mm时，失稳韧度增长了4.2%，增长了4.1%。

图5 断裂韧度散点图

图6 断裂能散点图

本文试验数据也符合以往研究的结论［14-15］。即当试件高度大于200 mm时，起裂韧度基本保持一致，没有明显的尺寸效应，而失稳韧度会随着试件尺寸的增大呈现微小的增大。张秀芳等［16］的研究结果显示，失稳韧度随着试件尺寸的增大而增大，但是增长幅度较小。徐世烺等曾引用文献［17-18］的试验数据计算双K断裂韧度，结果均显示，失稳断裂韧度表现出了一定的尺寸效应，其数值会随着试件尺寸的增大而增大。Ruiz等［19］提出了一种简单的方法计算起裂荷载，无需借助应变片等工具。他们的计算结果表明，起裂韧度随着试件尺寸的变化基本保持一致，失稳韧度随着试件尺寸的增加稍稍增大，但是幅度并不明显。文献［20］采用楔入劈拉试件对碾压混凝土进行了断裂力学试验，同样发现在试验范围内，失稳韧度随着试件尺寸的增大而增大。

为了研究断裂参数的尺寸效应是否是由于加载状态和附加弯矩引起的，徐世烺［14］等采用一种改进的楔入劈拉试件进行试验，改良后试件上各竖向力完全共线，在裂缝尖端不产生附加弯矩。试验结果表明，当试件高度大于200 mm时，双K断裂参数没有明显的尺寸效应。本文试验由于混凝土自身的重力在裂缝尖端同样会产生附加弯矩，影响裂缝尖端的应力应变场，因此本文试验失稳断裂韧度呈现出的尺寸效应可能是由于附加弯矩造成的。同时，从以往的研究可知，对于某一种试件形式，只有当试件高度大于某一范围时，断裂韧度才没有明显的尺寸效应。这个范围对于三点弯曲梁试件一般为200 mm；对于楔入劈拉试件一般为400 mm，本文试验失稳韧度表现出的尺寸效应有可能是因为试件高度尚未达到临界高度。因此的尺寸效应问题，有待进一步的试验研究。

5 结论

改进的圆形紧凑拉伸试件具有体积小、操作简便等优点，而且为既有结构的断裂性能评估提供可能。本文运用改进的圆形紧凑拉伸法进行了断裂力学试验，测定了双G断裂参数和双K断裂参数同时计算了由双G转化的双K参数通过分析可得出以下结论：（1）由改进的圆形紧凑拉伸试验得到的双G参数和双K断裂参数变异系数较小，数据离散性较小，说明改进的圆形紧凑拉伸法稳定可靠。（2）在数值上基本一致，误差在15%以内。这验证了双G断裂模型，说明两者在描述断裂性能方面具有很高的等效性。双G断裂模型是双K断裂模型在能量法方面的延伸和发展。（3）当试件直径不超过200 mm时，均随着试件尺寸的增大而增大。几乎保持不变。当试件直径大于200 mm时，在本次试验范围内，均没有明显的尺寸效应；随着试件尺寸的增大而增大，但是幅度较小。（4）当试件直径不超过200 mm时，随着试件尺寸的增大而增大，几乎保持不变。当试件直径大于200 mm时，没有明显的尺寸效应；随着试件尺寸的增大而增大。（5）的尺寸效应问题，有待进一步的试验研究。

参 考 文 献：

［1］管俊峰，李庆斌，吴智敏，等.现场浇筑大坝混凝土断裂参数与等效成熟度关系研究［J］.水利学报，2015，46（8）：951-959.

［2］吴瑶，徐世烺，吴建营，等.双K断裂准则在丹江口大坝安全性评定中的应用［J］.水利学报，2015，46（3）：366-372.

［3］XU S L，REINHARDT H W.Determination of double-K criterion for crack propagation in quasi-brittle fracture，Part I：Experimental investigation of crack propagation［J］.International Journal of Fracture，1999，98（2）：111-149.

［4］XU S L，REINHARDT H W.Determination of double-K criterion for crack propagation in quasi-brittle fracture，Part II：Analytical evaluating and practical measuring methods for three-point bending notched beams［J］.Inter⁃national Journal of Fracture，1999，98（2）：151-177.

［5］XU S L，REINHARDT H W.Determination of double-K criterion for crack propagation in quasi-brittle fracture，Part III：Compact tension specimens and wedge splitting specimens［J］.International Journal of Fracture，1999，98（2）：179-193.

［6］赵艳华，徐世烺，吴智敏.混凝土结构裂缝扩展的双G准则［J］.土木工程学报，2004，37（10）：13-18.

［7］KIM M，BUTTLAR W G，BAEK J，et al.Field and laboratory evaluation of fracture resistance of illinois hot-mix asphalt overlay mixtures［J］.Transportation Research Record Journal of the Transportation Research Board，2009，2127（2127）：146-154.

［8］ZOFKA A ，BRAHAM A.Comparison of low-temperature field performance and laboratory testing of 10 Test sec⁃tions in the Midwestern United States［J］.Transportation Research Record Journal of the Transportation Research Board，2009，2127（2127）：107-114.

［9］WAGNONER M P，BUTTLAR W G，PAULINO G H.Disk-shaped compact tension test for asphalt concrete frac⁃ture［J］.Experimental Mechanics，2005，45（3）：270-277.

［10］CIFUENTES H，LOZANO M，HOLUSOVA T，et al.Modified disk-shaped compact tension test for measuring concrete fracture properties［J］.International Journal of Concrete Structures and Materials，2017，11（2）：215-228.

［11］徐世烺.混凝土断裂力学［M］.北京：科学出版社，2011.

［12］FISHER D M，BUZZARD R J.Comparison tests and experimental compliance calibration of the proposed stan⁃dard round compact plane strain fracture toughness specimen［R］.NASA，1979.

［13］XU S L，ZHANG X F.Determination of fracture parameters for crack propagation in concrete using an energy ap⁃proach［J］.Engineering Fracture Mechanics，2008，75：4293-4308.

［14］徐世烺，卜丹，张秀芳.不同尺寸楔入式紧凑拉伸试件双K断裂参数的试验测定［J］.土木工程学报，2008，42（2）：70-76.

［15］吴智敏，徐世烺，王金来，等.三点弯曲梁法研究砼双K断裂参数及其尺寸效应［J］.水力发电学报，2000（4）：16-24.

［16］ZHANG X F，XU S L.A comparative study on five approaches to evaluate double-K fracture toughness parame⁃ters of concrete and size effect analysis［J］.Engineering Fracture Mechanics，2011，78：15-38.

［17］ZHAO G，JIAO H，XU S L.Study on fracture behavior with wedge splitting test method［C］//Fracture Processes in Concrete Rock Andceramics.London：E＆F.N.Spon.，1991.

［18］BRUHWILER E，WITTMANN F H.The wedge splitting test，a method of performing stable fracture tests［J］.En⁃gineering Fracture Mechanics，1990，35：17-26.

［19］RUIZ G，ORTEGA P J，YU RC，et al.Effect of size and cohesive assumptions on the double-K fracture parame⁃ters of concrete［J］.Engineering Fracture Mechanics，2016，166：198-217.

［20］CUI X，ZHANG J，HOU F，et al.Experimental study of the effect of RCC specimen size and crack depth ratio on double-K fracture parameters［J］.Journal of Testing&Evaluation，2015，43（2）：336-343.