本文于1999年11月4日收到。 众多的研究结果表明，试件高度是影响混凝土断裂韧度KIc的主要因素，而造成KIc尺寸效应的主要原因是忽略了混凝土裂缝失稳断裂前亚临界扩展阶段。为此，文献[1，2]采用激光散斑及光弹贴片等测试技术研究了混凝土裂缝失稳断裂前的亚临界扩展过程，结果表明，混凝土裂缝是先由起裂、稳定扩展然后失稳破坏，当试件高度h>1.0m时，其扩展量Δαc达210mm.根据试验中测得的起裂荷载Pini、Δαc及最大荷载Pmax获得了与试件高度无关的起裂断裂韧度KiniIc及失稳断裂韧度KunIc，据此提出了混凝土结构裂缝扩展的双K断裂准则，即当KI=KiniIc时，裂缝起裂；当KiniIc<KI≤KunIc时，裂缝稳定扩展；当KI>KunIc时，裂缝失稳扩展。然而，KiniIc、KunIc是根据大型试件断裂试验中测得的，它要求测得混凝土的Δαc、Pini，这在普通实验室是难以做到的。为此，文献［3］采用虚拟裂缝模型并结合线弹性断裂理论通过迭代计算了混凝土Δαc，并得到了当试件高度h>400mm时，KiniIc、KunIc及临界裂缝尖端张开位移CTODc值均与试件高度无关的结论。然而在以上的研究中，均未考虑混凝土初始缝高比α0/h的影响。因此本文采用尺寸400mm×400mm×２00mm，α0/h比分别为0.2、0.4、0.5、0.6的楔入劈拉试件研究α0/h对双K断裂参数KiniIc、KunIc及CTODc的影响。
1 混凝土有效裂缝长度的确定
由于混凝土在失稳断裂前存在着裂缝的稳定扩展阶段，使得失稳断裂前试件实际裂缝长度α大于预制缝长α0。若将失稳断裂前裂缝稳定扩展长度记为Δαc，则
|
但在实验中要精确测量Δαc值较困难，需要先进的测试技术。因此本文根据文献[4]的线性渐近叠加原理，采用线弹性断裂力学公式计算αc的值。
对于如图1所示的楔入劈拉试件，当外荷载P达最大值Pmax时，其对应的有效裂缝长度αc可由下式确定： | CMODc＝Pmax/tE[11.56(1-αc/h)-2－9.397)] | (2) | | 
图1 楔入劈拉试件
|
式中：t、h、E分别为试件厚度、高度及弹性模量。
2 混凝土双K断裂参数的确定
2.1 闭合力产生的应力强度因子KcI的计算 由于混凝土裂缝失稳扩展前存在着主裂缝的稳定扩展阶段，根据虚拟裂缝模型，当裂缝张开位移ω小于ω0时，尚能传递应力，这个应力称为闭合力。因此，试件除了受到外荷载P作用外还存在着阻止裂缝扩展的闭合力的作用(如图2所示).根据叠加原理，可将图2(a)分解为图2(b)及图2(c).对于图2(a)，其裂缝尖端处的应力强度因子为：
|
其中，KIp为由集中力P所产生的应力强度因子，可直接由线弹性断裂力学公式计算：
| 
| (4) | |
图2 试件受力特性
|
而闭合力σ(x)产生的应力强度因子(负值)可由文献［5］所给的公式计算：
| 
| (5) |
式中
| 
| (6) |
|
如图3所示，当裂缝处于临界状态时，P=Pmax，其裂缝尖端张开位移CTOD也达其临界值CTODc，则其相应的闭合力σ(x)在虚拟裂缝上的分布为：
| 
| (7) | | [ft-σs(CTODc)] (α0≤x≤αc) | |
图3 闭合力的分布 |
(7)式中:σs(CTODc)为当裂缝尖端张开位移达临界值CTODc时的应力，它可由图4的混凝土应力软化曲线得出[6]：
| 
| (8) |
式中：c1、c2为材料常数。当用CTODc代替式中的ω时即可得到σs值。对于普通混凝土，c1=3，c2=7，ω0=160μm.为简化计算，也可采用双线性模型，其表达式为：
| | { | σ= ft-(ft-σs)ω/ωs (0≤ω≤ωs) | |
| | σ= σs（ω0-ω/ω0－ωs） (ωs<ω≤ω0) | (9) |
| | σ= 0(ω>ω0) | |
| 其中，ωs＝0.8Gf/ft,σs＝ft/3，ω0＝3.6Gf/ft,ft和Gf分别为混凝土的抗拉强度及断裂能。将CTODc代入式(9)中的w即可得到σs(CTODc)值。而CTODc值则可由下式计算： | | 
| | 图4 混凝土应力软化曲线 | |
| CTODc＝（αc－α0）/αc ·CMODc | (10) |
因此，在临界状态时，在虚拟裂缝区内由闭合力σ(x)产生的断裂韧度KcIc为：
| 
| (11) |
2.2 混凝土双K断裂参数的确定 由式(3)可知，当裂缝处于临界状态时，
式中：KiniIc、KunIc称为起裂断裂韧度及失稳断裂韧度。其中失稳断裂韧度KunIc可将试验中测得的最大荷载Pmax及由式(2)计算的αc代入式(4)获得，而起裂断裂韧度KiniIc定义为混凝土裂缝起裂时的荷载Pini对应的应力强度因子，它表示材料抵抗裂缝扩展的能力。但要精确测量Pini比较困难。
对于理想弹塑性材料，根据DM模型，由于裂缝尖端附近出现较大的塑性区，使得裂缝尖端的应力奇异性消失，即KI＝KIp－KcI＝0，也就是说由闭合力σ(x)产生的负的应力强度因子完全抵消由拉应力产生的正的应力强度因子，而混凝土为半脆性材料，由σ(x)产生的负的应力强度因子不能完全抵消由拉应力产生的正的应力强度因子，因此KI≠0。这样，若在试验中测得混凝土试件的抗拉强度ft、最大荷载Pmax及对应的裂缝口张开位移CMODc、弹性模量E等参数，根据以上各式便可获得混凝土的双K断裂参数。
3 试验结果
试件尺寸均为400mm×400mm×200mm，初始缝高比α0/h分别为0.2、0.4、0.5、0.6，每组试件为3个。试件所用材料为大连产河砂，石灰岩碎石，最大骨料粒径为20mm，水泥为大连水泥厂生产的425＃普通硅酸盐水泥，试件配合比为水泥∶砂子∶石子∶水=1∶1.73∶3.01∶0.52.试验时测得150×150×150mm3立方体抗压强度fcu为38.8MPa，抗拉强度ft为3.0MPa，弹性模量E为32.2GPa.试验测得的各试件Pmax、CMODc、Gf等参数见表1.
| | 表1 试验结果 | |
|
|
| 试件编号 | α0/h | Pmax/N | CMODc/mm | Gf/(N·m-1) |
|
|
| ωs1
ωs2
ωs3
ωs4 | 0.2
0.4
0.5
0.6 | 34244.7
25206.0
16042.1
10779.3 | 0.1445
0.1712
0.1673
0.2364 | 488.2
437.8
376.7
291.5 |
|
|
4 混凝土双K断裂参数计算结果
根据试验测得的Pmax、ft、E、CMODc、Gf等参数即可计算混凝土的αc、KiniIc、KunIc及CTODc等值。详细结果见表2.从表2中可以发现，试件初始缝高比的不同，裂缝亚临界扩展量Δαc也有所不同，而KunIc值则与试件的α0/h无关；当试件缝高比α0/h≥0.4时，所计算的KiniIc、CTODc也基本上与α0/h无关，但采用不同的闭合力方程所计算的KiniIc却有所不同，可见断裂能Gf及σ-ω曲线形状对KiniIc有明显影响。
| | | 表2 双K断裂参数计算结果 | | |
|
|
| | | | KcIc/(MPa·m1/2) | | KiniIc/(MPa·m1/2) | KunIc/(MPa·m1/2) | CTODc/mm |
| 试件编号 | αc/mm | Δαc/mm |
| |
|
| | | | 方程(8) | 方程(9) | | 方程(8) | 方程(9) | | |
|
|
| ωs1
ωs2
ωs3
ωs4 | 175.1
213.4
344.6
289.2 | 95.1
53.4
44.6
49.3 | 1.3538
0.9164
0.8378 0.9804 | 2.0000 1.3969
1.2865
1.4836 | | 1.0063
1.3825
1.0863
1.1675 | 0.3601
0.9021
0.6376
0.6642 | 2.3601
2.2990
1.9241
2.1479 | 0.07848
0.04286
0.03049
0.04021 |
|
|
5 结 论
根据本文的计算及分析，可以得到如下结论：(1)混凝土结构的裂缝在起裂后，由于闭合力的作用使之在失稳断裂前存在较大的亚临界扩展阶段，而不象脆性材料一经起裂便失稳扩展。(2)混凝土失稳断裂韧度KunIc与试件初始缝高比无关，可以作为混凝土材料的断裂参数。(3)当试件初始缝高比α0/h≥0.4时，根据虚拟裂缝模型计算的混凝土KiniIc及临界裂缝尖端张开位移CTODc值也与试件α0/h无关。(4)采用不同的闭合力方程所计算的KiniIc相差较大，因此在计算KiniIc时应考虑断裂能Gf及σ-ω曲线形状的影响。(5)当α0/h≥0.4时，KiniIc、KunIc及CTODc可以作为混凝土的材料常数，因此在测定KiniIc、KunIc及CTODc时α0/h宜取0.4～0.6。
参 考 文 献：
［1］ 徐世烺，赵国藩。混凝土裂缝的稳定扩展过程与临界裂缝尖端张开位移[J]。水利学报，1989，(4).
［2］ 徐世烺，赵国藩。光弹贴片法研究混凝土裂缝扩展过程[J]。水力发电学报，1991，(3).
［3］ 吴智敏，徐世烺，王金来。基于虚拟裂缝模型的混凝土双K断裂参数[J]。水利学报，1999，(7).
［4］ Shilang Xu，Reinhardt H W. Determination ofdouble K criterion for crack propagation in quasi brittle fracture, part Ⅱ:Analytical evaluating andpractical measuring methods for three point bendingnotched beams[J]。InternationalJournal of Fracture,1999,98(2).
［5］ 中国航空研究院主编。应力强度因子手册[M]。北京：科学出版社，1981.
［6］ Reinhardt H W, Cornelissen H A W, Hordijk D A.Tensile tests andfailure analysis of concrete［J］。Journal of Structural Engineering, 1986,112(11).
