水利工程论文-试件初始缝长对混凝土双K断裂参数的影响.doc

水利工程论文-试件初始缝长对混凝土双K断裂参数的影响摘要：本文采用尺寸为400mm400mm200mm，初始缝高比分别为0.2、0.4、0.5、0.6的楔入劈拉试件，应用线性渐近叠加原理，利用试验中测得的最大荷载PmaxCMODc计算了混凝土裂缝的亚临界扩展量c，在此基础上利用虚拟裂缝模型并结合线弹性断裂理论，采用不同的闭合力方程计算了混凝土起裂断裂韧度KiniIc、失稳断裂韧度KunIc及临界裂缝尖端张开位移CTODc.结果表明，采用不同闭合力方程所计算的混凝土KiniIc值有所差异，但当试件缝高比大于或等于0.4时KiniIc及CTODc均与试件初始缝高比无关，而KunIc值则与初始缝高比无关，这表明双K断裂参数及CTODc可以作为描述混凝土裂缝起裂、稳定扩展及失稳破坏全过程的材料常数。关键词：混凝土裂缝断裂参数虚拟裂缝模型众多的研究结果表明，试件高度是影响混凝土断裂韧度KIc的主要因素，而造成KIc尺寸效应的主要原因是忽略了混凝土裂缝失稳断裂前亚临界扩展阶段。为此，文献1，2采用激光散斑及光弹贴片等测试技术研究了混凝土裂缝失稳断裂前的亚临界扩展过程，结果表明，混凝土裂缝是先由起裂、稳定扩展然后失稳破坏，当试件高度h1.0m时，其扩展量c达210mm.根据试验中测得的起裂荷载Pini、c及最大荷载Pmax获得了与试件高度无关的起裂断裂韧度KiniIc及失稳断裂韧度KunIc，据此提出了混凝土结构裂缝扩展的双K断裂准则，即当KI=KiniIc时，裂缝起裂；当KiniIcKunIc时，裂缝失稳扩展。然而，KiniIc、KunIc是根据大型试件断裂试验中测得的，它要求测得混凝土的c、Pini，这在普通实验室是难以做到的。为此，文献3采用虚拟裂缝模型并结合线弹性断裂理论通过迭代计算了混凝土c，并得到了当试件高度h400mm时，KiniIc、KunIc及临界裂缝尖端张开位移CTODc值均与试件高度无关的结论。然而在以上的研究中，均未考虑混凝土初始缝高比0/h的影响。因此本文采用尺寸400mm400mm00mm，0/h比分别为0.2、0.4、0.5、0.6的楔入劈拉试件研究0/h对双K断裂参数KiniIc、KunIc及CTODc的影响。1混凝土有效裂缝长度的确定由于混凝土在失稳断裂前存在着裂缝的稳定扩展阶段，使得失稳断裂前试件实际裂缝长度大于预制缝长0。若将失稳断裂前裂缝稳定扩展长度记为c，则c=0ac(1)但在实验中要精确测量c值较困难，需要先进的测试技术。因此本文根据文献4的线性渐近叠加原理，采用线弹性断裂力学公式计算c的值。对于如图1所示的楔入劈拉试件，当外荷载P达最大值Pmax时，其对应的有效裂缝长度c可由下式确定：CMODcPmax/tE11.56(1-c/h)-29.397)(2)图1楔入劈拉试件式中：t、h、E分别为试件厚度、高度及弹性模量。2混凝土双K断裂参数的确定2.1闭合力产生的应力强度因子KcI的计算由于混凝土裂缝失稳扩展前存在着主裂缝的稳定扩展阶段，根据虚拟裂缝模型，当裂缝张开位移小于0时，尚能传递应力，这个应力称为闭合力。因此，试件除了受到外荷载P作用外还存在着阻止裂缝扩展的闭合力的作用(如图2所示).根据叠加原理，可将图2(a)分解为图2(b)及图2(c).对于图2(a)，其裂缝尖端处的应力强度因子为：KIKIpKcI(3)其中，KIp为由集中力P所产生的应力强度因子，可直接由线弹性断裂力学公式计算：(4)图2试件受力特性而闭合力(x)产生的应力强度因子(负值)可由文献5所给的公式计算：(5)式中(6)如图3所示，当裂缝处于临界状态时，P=Pmax，其裂缝尖端张开位移CTOD也达其临界值CTODc，则其相应的闭合力(x)在虚拟裂缝上的分布为：(7)ft-s(CTODc)(0xc)图3闭合力的分布(7)式中:s(CTODc)为当裂缝尖端张开位移达临界值CTODc时的应力，它可由图4的混凝土应力软化曲线得出6(8)式中：c1、c2为材料常数。当用CTODc代替式中的时即可得到s值。对于普通混凝土，c1=3，c2=7，0=160m.为简化计算，也可采用双线性模型，其表达式为：=ft-(ft-s)/s(0s)=s（0-/0s）(s0)其中，s0.8Gf/ft,sft/3，03.6Gf/ft,ft和Gf分别为混凝土的抗拉强度及断裂能。将CTODc代入式(9)中的w即可得到s(CTODc)值。而CTODc值则可由下式计算：图4混凝土应力软化曲线CTODc（c0）/cCMODc(10)因此，在临界状态时，在虚拟裂缝区内由闭合力(x)产生的断裂韧度KcIc为：(11)2.2混凝土双K断裂参数的确定由式(3)KiniIcKunIcKcIc(12)式中：KiniIc、KunIc称为起裂断裂韧度及失稳断裂韧度。其中失稳断裂韧度KunIc可将试验中测得的最大荷载Pmax及由式(2)计算的c代入式(4)获得，而起裂断裂韧度KiniIc定义为混凝土裂缝起裂时的荷载Pini对应的应力强度因子，它表示材料抵抗裂缝扩展的能力。但要精确测量Pini比较困难。对于理想弹塑性材料，根据DM模型，由于裂缝尖端附近出现较大的塑性区，使得裂缝尖端的应力奇异性消失，即KIKIpKcI0，也就是说由闭合力(x)产生的负的应力强度因子完全抵消由拉应力产生的正的应力强度因子，而混凝土为半脆性材料，由(x)产生的负的应力强度因子不能完全抵消由拉应力产生的正的应力强度因子，因此KI0。这样，若在试验中测得混凝土试件的抗拉强度ft、最大荷载Pmax及对应的裂缝口张开位移CMODc、弹性模量E等参数，根据以上各式便可获得混凝土的双K断裂参数。3试验结果试件尺寸均为400mm400mm200mm，初始缝高比0/h分别为0.2、0.4、0.5、0.6，每组试件为3个。试件所用材料为大连产河砂，石灰岩碎石，最大骨料粒径为20mm，水泥为大连水泥厂生产的425普通硅酸盐水泥，试件配合比为水泥砂子石子水=11.733.010.52.试验时测得150150150mm3立方体抗压强度fcu为38.8MPa，抗拉强度ft为3.0MPa，弹性模量E为32.2GPa.试验测得的各试件Pmax、CMODc、Gf等参数见表1.表1试验结果试件编号0/hPmax/NCMODc/mmGf/(Nm-1)s1s2s3s40.20.40.50.634244.725206.016042.110779.30.14450.17120.16730.2364488.2437.8376.7291.54混凝土双K断裂参数计算结果根据试验测得的Pmax、ft、E、CMODc、Gf等参数即可计算混凝土的c、KiniIc、KunIc及CTODc等值。详细结果见表2.从表2中可以发现，试件初始缝高比的不同，裂缝亚临界扩展量c也有所不同，而KunIc值则与试件的0/h无关；当试件缝高比0/h0.4时，所计算的KiniIc、CTODc也基本上与0/h无关，但采用不同的闭合力方程所计算的KiniIc却有所不同，可见断裂能Gf及-曲线形状对KiniIc有明显影响。表2双K断裂参数计算结果KcIc/(MPam1/2)KiniIc/(MPam1/2)KunIc/(MPam1/2)CTODc/mm试件编号c/mmc/mm方程(8)方程(9)方程(8)方程(9)s1s2s3s4175.1213.4344.6289.295.153.444.649.31.35380.91640.83780.98042.00001.39691.28651.48361.00631.38251.08631.16750.36010.90210.63760.66422.36012.29901.92412.14790.078480.042860.030490.040215结论根据本文的计算及分析，可以得到如下结论：(1)混凝土结构的裂缝在起裂后，由于闭合力的作用使之在失稳断裂前存在较大的亚临界扩展阶段，而不象脆性材料一经起裂便失稳扩展。(2)混凝土失稳断裂韧度KunIc与试件初始缝高比无关，可以作为混凝土材料的断裂参数。(3)当试件初始缝高比0/h0.4时，根据虚拟裂缝模型计算的混凝土KiniIc及临界裂缝尖端张开位移CTODc值也与试件0/h无关。(4)采用不同的闭合力方程所计算的KiniIc相差较大，因此在计算KiniIc时应考虑断裂能Gf及-曲线形状的影响。(5)当0/h0.4时，KiniIc、KunIc及CTODc可以作为混凝土的材料常数，因此在测定KiniIc、KunIc及CTODc时0/h宜取0.40.6。参考文献：1徐世烺，赵国藩。混凝土裂缝的稳定扩展过程与临界裂缝尖端张开位移J。水利学报，1989，(4).2徐世烺，赵国藩。光弹贴片法研究混凝土裂缝扩展过程J。水力发电学报，1991，(3).3吴智敏，徐世烺，王金来。基于虚拟裂缝模型的混凝土双K断裂参数J。水利学报，1999，(7).4ShilangXu，ReinhardtHW.DeterminationofdoubleKcriterionforcrackpropagationinquasibrittlefracture,part:Analyticalevaluatingandpracticalmeasuringmethodsforthreepointbendingnotchedbeam