基于应力-渗流-损伤耦合模型的重力坝三维水力劈裂数值模拟.pdf 14页

'中国科技论文在线http://www.paper.edu.cn基于应力-渗流-损伤耦合模型的重力坝三维水力劈裂数值模拟1,22**沙莎，张国新5（1.清华大学水利水电工程系，北京100084；2.中国水利水电科学研究院流域水循环模拟与调控国家重点实验室，北京100038）摘要：裂缝的高压水力劈裂是混凝土高坝安全评估的重要部分。目前，重力坝水力劈裂的数值模拟绝大多数是二维的，对于坝体上游面经常出现的竖直裂缝的三维水力劈裂的数值研究几乎为零。本文提出一种应力-渗流-损伤耦合模型，可以模拟重力坝三维水力劈裂现象。采10用该耦合模型，模拟了一个内置裂缝的圆柱体试件的水力劈裂，数值结果与试验结果吻合很好，表明该耦合模型可以很好的模拟水力劈裂。在此基础上，采用该耦合模型，模拟了观音岩重力坝劈头裂缝的形成，数值结果表明，随着上游水位的抬升，裂缝逐渐向坝体上部、深部扩展，正常蓄水位时，损伤区域显著扩大，裂缝面呈圆弧状，与设计院预测的范围基本吻合。本文耦合模型可以有效地模拟重力坝的三维水力劈裂现象，可以预测施工期形成的表面15裂缝在蓄水期是否扩展、扩展过程及最终扩展范围，评估表面裂缝的危险程度以及时采取有效的处理措施，为高重力坝的安全运行提供技术保障。关键词：重力坝；水力劈裂；三维；应力-渗流-损伤；耦合中图分类号：TV31320Simulationof3Dhydraulicfracturingofconcretegravitydamsbasedonstress-seepage-damagecouplingmodel1,22SHASha,ZHANGGuoxin(1.TsinghuaUniversity,DepartmentofHydraulicEngineering,Beijing,100084;2.ChinaInstituteofWaterResourcesandHydropowerResearch,StateKeyLaboratoryof25SimulationandRegulationofWaterCycleinRiverBasin,Beijing,100038)Abstract:Highpressurehydraulicfracturingisanimportantpartofsafetyassessmentinhighconcretedams.Atpresent,mostofthenumericalsimulationsofhydraulicfracturinginconcretegravitydamsarelimitedto2D.Thereislackofsimulationof3Dhydraulicfracturingofverticalcrackswhichoftenoccurontheupstreamofthedam.Inthisstudy,astress-seepage-damagecouplingmodelispresented30forthesimulationof3Dhydraulicfracturing.Bythecouplingmodel,hydraulicfracturingofaninternalcrackinacylinderspecimenissimulatedandtheresultsareingoodagreementwiththeexperimentalresults,showingthatthecoupledmodelcansimulatehydraulicfracturingwell.Thenthecoupledmodelisusedtosimulate3DhydraulicfracturingofGuanyinyanconcretegravitydam.Thenumericalresultsshowthatwiththerisingofwaterlevel,thecrackextendstotheupperanddeeper35partofthedamgradually.Atthenormalwaterlevel,thedamageareasignificantlyexpandsandisinarcshape,whichisconsistentwiththeareaforecastedbythedesigninstitute.Thecoupledmodelcansimulate3Dhydraulicfracturingintheconcretegravitydamseffectively,predictwhethersurfacecracksoccurringintheconstructionperiodwouldextendduringthestorageperiod,theextensionprocessandthecrackedarea,assessthedegreeofriskofsurfacecrackstotakeeffectivemeasuresin40timeandprovidetechnicalsupportforthesafeoperationofhighdams.Keywords:Gravitydam;Hydraulicfracture;3D;Stress-Seepage-Damage;Coupling作者简介：沙莎（1989-），女，博士研究生，水工结构数值计算通信联系人：张国新(1960-)，男，博士，教授级高级工程师，博士生导师，主要从事水工结构.E-mail:zhanggx@iwhr.com-1- 中国科技论文在线http://www.paper.edu.cn0引言45随着我国水资源开发需求的增长以及筑坝技术的发展，修建200~300m的超高坝已成为可能，混凝土高坝的安全问题成为各国坝工界一个十分关注的重要课题。由于混凝土材料抗拉强度较低，大体积混凝土坝施工期极容易产生裂缝，尤其是碾压混凝土坝，还存在薄弱层面，这些裂缝和薄弱面在库水压力作用下是否会发生高压水劈裂是一个亟需研究的问题。国内外已有几座重力坝出现了高压水劈裂的现象，例如：加拿大的雷威尔斯托克坝、美国的德50沃歇克坝和卢塞尔坝，施工过程中虽然采用了预冷骨料、水管冷却、表面保温等综合温度控制措施，坝体上游面仍出现了表面裂缝，水库蓄水后经过一段时间，表面裂缝突然大范围扩[1]大，发展为劈头裂缝，深入坝内几十米，有的甚至将整个坝段一分为二，引起严重漏水。[2][3][4]我国柘溪、桓仁大头坝，玉石、观音岩碾压混凝土重力坝也出现了同样的问题。因此，对混凝土重力高坝或超高坝进行安全评估时，考虑高压水劈裂作用的影响是非常重要的。55在混凝土高压水劈裂的试验研究方面，已有学者进行了相关研究。Brühwiler和Saouma[5-6]研究了不同级配的混凝土裂缝内的水压分布变化规律及裂缝中静水压力对混凝土表观断[7]裂能、表观断裂韧度的影响。Slowik和Saouma研究了裂缝边缘张开速度不同时裂缝内水[8]压分布随时间的变化以及裂缝突然闭合时对水压分布的影响。王建敏采用大型尺寸楔入式紧凑拉伸试件进行了水力劈裂试验，研究了裂缝内的静水压力对断裂性能的影响，指出随着[9]60水压力的增加，试件所能承受的最大荷载减小。贾金生等提出了全级配混凝土试件单轴拉、压应力作用下高压水劈裂模拟试验新方法，并推导了判定重力坝坝踵是否会发生水力劈裂的分析公式。通过试验可以揭示水力劈裂的机理，但是由于试验中采取的试件均是立方体、圆柱体试件，所以试验不能真实反映混凝土坝的形态及其受力条件，也无法预测裂缝扩展范围及危险程度。[10]65鉴于试验的不足，很多学者致力于混凝土坝高压水劈裂的数值研究。陈胜宏等采用有限元法研究了高压水劈裂对小湾高拱坝坝踵开裂的影响，指出考虑高压水劈裂坝踵裂缝扩[11]展范围比不考虑时大。Barpi和Valente对混凝土重力坝建基面处裂纹的水力劈裂进行了模[12]拟研究，讨论了断裂过程区的发展对裂缝内水压的影响。董玉文等采用扩展有限元法进[13]行了向家坝坝踵水力劈裂的数值分析，将裂缝面上的水压力视为均布荷载。王克峰等采[14]70用扩展有限元法研究了考虑流固耦合效应的重力坝的水力劈裂问题。SecchiandSchrefler基于分离裂缝模型和三维自动网格重构技术提出一种模拟饱和孔隙介质三维水力劈裂的方法，并且用于混凝土重力坝在水头超载作用下坝踵的开裂分析。水力劈裂是个复杂问题，涉及到四个耦合过程：(1)裂纹面上水压力引发的介质变形；(2)裂纹内及周围介质中流体流动；(3)裂纹的扩展；(4)裂纹内的水渗透到周围介质中[15]。上述混凝土坝高压水劈裂的数值研[14][10][12]75究中，只有文献考虑了这四个耦合过程。文献和将裂缝面上的水压力视为均布荷载，[11][13]没有考虑裂缝扩展过程中水体流动引起的裂缝内的水压重分布。文献和则假定裂缝面是不透水的，且没有考虑周围介质中流体流动。而且已有的重力坝水力劈裂的数值研究大多-2- 中国科技论文在线http://www.paper.edu.cn是二维的，三维方面的十分罕见，特别是重力坝劈头裂缝方面的。[16][17][18]目前模拟水力劈裂的数值方法较多，如相场法，近场动力学，元胞自动机法，[19][20][21][14][22]80离散元法，流形元法，无单元法，有限元法，边界元法。有限元法是最成熟，应用最广的数值方法。有限元法本质上是一种连续介质力学方法，必须加以改进才能被用来分析裂缝扩展这样的不连续问题，改进的方法分为两类：变网格法和固定网格法。变网格法随着裂缝的产生，需要不断地重新划分网格，裂纹面必须与单元边保持一致，在裂尖区域需要细化单元或改用特殊形式的缝端奇异元，产生极大的前处理工作量。固定网格法则是保持85网格不变，通过修改开裂单元插值关系、本构关系反映裂缝的存在，比如扩展有限元(XFEM)、连续损伤模型。相比较而言，固定网格法的应用更为方便。据作者所知，扩展有限元法目前仅限于二维水力劈裂的模拟，采用扩展有限元法进行三维水力劈裂模拟的研究还没有。本文基于有限元采用一种各向同性损伤模型进行水力劈裂的模拟，将混凝土视为饱和的孔隙介质，根据孔隙介质有效应力原理，提出一种应力-渗流-损伤耦合模型，考虑水力劈裂90过程中的耦合效应。采用该耦合模型，本文进行了观音岩重力坝上游面垂直裂缝水力劈裂的数值模拟。1应力-渗流-损伤耦合模型本文假定大坝混凝土是饱和的孔隙介质。文中应力-渗流-损伤耦合模型包括了以下非线性行为：(1)单元损伤后的应力-应变关系；(2)损伤引起的孔隙水压影响系数的改变；(3)应力、95损伤引起的单元渗透系数的改变。计算中假定受压时，应力应变关系是线弹性的，因为重力坝压应力一般不会超过其抗压强度。以下介绍了这种耦合模型的特性。1.1损伤模型混凝土的断裂特性表现为在真实的缝端前面存在一个断裂过程区，在断裂过程区中，通过骨料的咬合和界面间的摩擦力，相互之间仍有应力作用，断裂过程区的存在使得混凝土呈[23]100现出应变软化特性。由Hillerborg等提出的虚拟裂缝模型可以很好地模拟断裂过程区的力学特性。目前，关于混凝土拉伸曲线的数学表达式，大多数学者主张上升段采用直线，主要区别在下降段，包括：单直线下降、分段下降、曲线下降。但无论采用何种形式，均应保持[24][24]拉伸曲线的断裂能相同。本文下降段采用由江见鲸提出的指数下降型表达式。因此，单轴受拉条件下，应力-应变关系表达为：Eεε,≤ε00tt105σ=−−αεε()(1)fet0,εε>tt00式中：E0为材料无损时的弹性模量；ft0为抗拉强度；εt为单元的单轴拉伸应变；ε0为开裂时的应变，ε0=ft0/E0；α为控制下降段的软化系数。断裂能Gf的含义如图1所示，为应力-裂缝宽度曲线与坐标轴围成的面积，可表达为式(2)；单位开裂宽度的断裂能gf的含义如图2中阴影所示，为应力-应变曲线下降段与ε=ε0、横坐标110轴ε围成的面积，可表达为式(3)。Gdf=σw，wd=ε(2)h-3- 中国科技论文在线http://www.paper.edu.cn式中：w为过程区内所有微裂缝张开的位移量之和，h为混凝土微裂缝分布区宽度。∞g=σεd(3)ftε0∞∞−−αεε()f0根据式(1)gdf==σεedt0ε=t(4)fttt0εε00ασft0σGfgf115wε0ε图1应力-裂缝宽度曲线和断裂能GF图2单位开裂宽度的断裂能gfFig.1Stressversuscrackwidthcurve,andFig.2FractureenergyperunitcrackwidthgffractureenergyGF断裂能GF与单位开裂宽度的断裂能gf存在如下关系：120Gg=l(5)Fft式中:lt为混凝土单元特征尺寸。对于平面单元，取积分点区域面积的平方根；对于实体单元，取积分点区域体积的立方根。由式(3)、式(4)，本构曲线下降段的软化系数应满足fltt0α=(6)GF[25]125根据Lemaitre提出的等效应变原理，假定Cauchy应力σ作用在受损材料上的应变与有效应力σ作用在无损材料上的应变等价，即σσσε===(7)EE(1−dE)00得到σε=(1−dE)(8)0式中：E0为材料无损时的弹性模量；E为材料受损后的弹性模量；d为损伤变量，d=0对应130于无损伤状态，d=1对应于完全损伤状态，0(9)t0εt-4- 中国科技论文在线http://www.paper.edu.cn1.2孔隙介质有效应力原理有效应力概念应用于被单相流体浸润的岩石、混凝土材料，可以看成太沙基有效应力原140理应用于土体的推广。基于有效应力原理，Biot首次引进一个标量参数，即Biot系数，反映[27-29][28]孔隙水压对有效应力的影响(Biot,1941,1955,1977)。在文献中，孔隙率bi被定义为有效孔隙面积与横截面Ai的比值。有效孔隙面积则被定义为在垂直于横截面Ai的方向，单位长度所有微小孔隙面积的总和。孔隙率bi同样可以代表孔隙体积(Vp)与代表体体积(Vb)的比值。在混凝土坝弹性分析中，通常假定含有孔隙的大坝混凝土是均质、各向同性的，因而经常采用145bi的一个各向同性值bx=by=bz=b0进行混凝土坝的弹性分析。对于土体的应力计算，通常取b0≈1。然而，对于弹性范围内混凝土的应力计算，取b0≈1则是不合理的，因为混凝土即使在即将破坏的状态下，b0仍小于1。tσbpσyyytσxbpxσx1pp1tt总应力{σ}自平衡应力{b}p有效应力{σ}={σ}-{b}p150图3饱和孔隙介质的应力分解Fig.3Decompositionofstressesinsaturatedporousmedia文中应力符号以拉应力为正。饱和孔隙介质的应力分解如图3所示，作用在单元表面的总应力可以分解成两部分：与内部水压力p相互平衡的外力bp和有效应力，单元的平均变形只与有效应力有关，可表示为：t155{σε}={D}{}，{σσ}=-{}{bp}(10)T式中：{σ}为有效应力向量，{σ}={σx,σy,σz,τxy,τyz,τzx}；{D}为弹性矩阵；{ε}={εx,εy,εz,γxy,Ttγyz,γzx}；{σ}为总应力向量；p为孔隙水压力；{b}为孔隙水压影响向量，三维分析中被定义为{b}={bx,by,bz,0,0,0}。1601.3孔隙水压影响系数与损伤关系单元未损伤时假定bx=by=bz=b0，单元损伤后，垂直于裂缝方向的孔隙水压影响系数与平行于裂缝方向孔隙水压影响系数理应是不一样的。图4定性的展示了开裂对孔隙水压影响系数的影响，1，2，3分别代表第一、二、三主应力的方向，建立以主应力方向为坐标轴的局部坐标系。开裂对垂直于裂缝面方向的b1影响较大，b2和b3在一个范围内，即b0＜b2=b3＜[30][31]165b1。BaryandBournazel通过实验得出了一个各向异性的Biot张量，其与各个方向的损伤值和孔隙水压有关。由于缺乏开裂对平行于裂缝方向的孔隙水压影响系数影响的试验资料，本文假定混凝土损伤之后，孔隙水压影响系数仍各向同性，即b1=b2=b3。弹性状态下，b1=b2=b3=b0，b0为初始Biot系数；完全损伤后，b1=b2=b3=1，孔隙水压影响系数随损伤的演化方程-5- 中国科技论文在线http://www.paper.edu.cn可表示为：170bbbbdb===+−(1)(11)12300孔孔孔孔孔孔孔孔孔孔1渗渗渗孔孔32开开开开图4开裂对孔隙水压影响系数及渗透性的影响Fig.4Effectsofcrackingonpore-pressure-influencecoefficientsandpermeability1.4渗透系数与应力及损伤关系175混凝土在微观上是由骨架和孔隙组成的，这种构造特征使得其在受荷载或扰动作用后，其微观几何形态、孔隙的结构发生改变，从而导致孔隙率和渗透性发生改变。孔隙率的变化主要由两部分组成：一是由于结构形变引起孔隙体积变化；二是微裂纹等缺陷的萌生、扩展及贯通，使得材料内部的孔隙结构和大小发生改变，即材料发生损伤而引起的孔隙结构变化[32]。180当单元未发生拉损伤，处于弹性状态，即损伤变量D=0时，渗透系数与有效应力成指数关系，表示如下：kk=ea01σ110kk=ea02σ(12)220aσkk=e03330式中：k0为初始渗透系数；a0为耦合系数；σ1,σ2,σ3为有效主应力，以受拉为正；k11,k22,k33为方向与有效主应力一致的主渗透系数。185当单元第一主应力达到抗拉强度，发生损伤，即损伤变量D>0时，单元中产生裂缝，图4定性地展示了开裂对单元渗透性的影响。混凝土中的水主要沿着裂缝流动，方向2，方向3的渗透系数大大提高，此时近似按照裂隙渗流考虑，采用在裂隙渗流研究中运用得最为广泛的单裂隙平行板水力模型，则渗透系数表示如下：=a01σkke110gu2(13)kk==223312υ190式中：u为裂缝张开的宽度，近似按u=(ε1-εt)lt计算，εt、lt的定义见2.1节；υ为水的运动粘滞系数，g为重力加速度。由k11,k22,k33形成局部坐标系下的渗透矩阵[k’]，经坐标转换将其转换到整体坐标系中，则整体坐标系下的渗透矩阵[k]可表示为：T[][]["][]kk=ββ(14)kkkxxxyxzk1100lmn111式中：[]kkkk=，["]kk=00，[]β=lmn。[β]为两种yxyyyz22222kkk00klmnzxzyzz33333-6- 中国科技论文在线http://www.paper.edu.cn195坐标系间的转换矩阵，li、mi、ni(i=1,2,3)为主应力σi的方向余弦。1.5渗流场基本微分方程假定水不可压缩，根据Dacy定律，可得无源非稳定情况下渗流连续方程：∂∂HH∂kS=(15)ij∂∂xxt∂ij式中：S为饱和度，H为水头势，kij为渗透系数张量；200对于非恒定渗流，其定解条件包括如下边界条件和初始条件：水头边界：Hxyzt(,,,)=>fxyzt(,,,)t0(16)Γ1∂H流量边界：kq=>(,,,)xyztt0(17)∂ηΓ2∂H隔水边界：=0(18)∂η初始条件：Hxyzt(,,,)=Hxyz(,,)(19)t=00205式中：f为已知水头边界值，q为已知流量边界，H0为已知函数。如果按照恒定渗流计算，即式(14)的右侧为0，没有初始条件式(18)1.6程序实现SAPTIS软件是中国水利水电科学研究院结构材料所独立开发的大型结构多场仿真与[33-35]非线性分析软件，该软件具备混凝土坝温度、渗流、变形、应力等多场耦合仿真分析功能。210本文在SAPTIS原有功能基础上加入了应力-渗流-损伤耦合非线性功能，进行水力劈裂的数值模拟计算。2试验验证[9]为了验证本文应力-渗流-损伤耦合模型的正确性，对贾金生等做的圆柱体混凝土试件水力劈裂试验进行了数值模拟，试件尺寸及受力条件如图5所示，本文计算模型如图6所示。[9][9]215选取文献中3个试件进行数值模拟，试件力学参数、受力状态与文献一致，见表1。水压[9]加载过程与文献一致，如图7所示，即采用梯级加载方式，水压小于1.0MPa时，每小时增加0.5MPa；水压1MPa至2MPa时，每小时增加0.2MPa，水压大于2MPa时，每小时增加-90.1MPa。渗流场计算按照非恒定渗流考虑，混凝土渗透系数k0=5×10m/s，耦合系数a0=0.01。-7- 中国科技论文在线http://www.paper.edu.cn预制裂缝220注：1、2钢管；3混凝土；4预设裂缝；5单轴压力图6计算模型图5高压水劈裂试验（单位：mm）Fig.6ThecalculationmodelFig.5Highpressurehydraulicfracturingtest(unit:mm)表1圆柱体混凝土试件力学参数及受力状态Tab.1Themechanicalparametersandloadingstateofcylindricalconcretespecimens试件混凝土单轴弹性模量断裂能Gf试件应力施加在试件编号龄期抗拉强度E(GPa)(N/m)状态两端的应力(d)ft(MPa)(MPa)Ⅲ904.143.797.9无应力0Ⅳ283.133.494.2单轴压1.0Ⅴ904.443.797.9单轴拉-0.71225图7水压加载过程图8P=3.4MPa，试件损伤分布Fig.7TheloadingprocessofwaterpressureFig.8ThedamagedistributionofspecimenⅣwhenP=3.4MPa230(a)P=1.4MPa(b)P=2.0MPa(c)P=2.5MPa-8- 中国科技论文在线http://www.paper.edu.cn(d)P=3.0MPa(e)P=3.1MPa(f)P=3.4MPa图9不同水压作用下，预制缝所在截面损伤分布Fig.9Thedamagedistributionofthesectionwheretheprefabricatednotchislocatedunderdifferentpressure235(a)P=1.4MPa(b)P=2.0MPa(c)P=2.5MPa(d)P=3.0MPa(e)P=3.1MPa(f)P=3.4MPa图10不同水压作用下，预制缝所在截面孔隙水压分布240Fig.10Theporepressuredistributionofthesectionwheretheprefabricatednotchislocatedunderdifferentpressure试件Ⅳ发生水力劈裂破坏时的损伤分布如图8所示，由图8可以看出，试件损伤区域几乎和预制缝平行。不同水压作用下，试件Ⅳ预制缝所在截面损伤分布如图9所示，孔隙水压分布如图10所示。由图9可知，水压1.4MPa，缝端开始损伤，随着水压的增加，损伤区域越来[9]245越大，水压3.4MPa，损伤区域面积占总面积55%，与文献描述的劈裂面积十分接近，水压超过3.4MPa，截面全部损伤，本文认为水压3.4MPa时，试件发生了水力劈裂破坏。由图10可以看出，随着水压增加，最大孔隙水压分布区域逐渐增大，其边界与图9中损伤区域的边界十分吻合，这是由混凝土单元损伤之后，渗透系数显著增加，水不断渗入引起的。不同试件水力劈裂水压的试验值、计算值见表2，计算值与试验值吻合很好，可以看出，本文耦合250模型可以很好地模拟水力劈裂现象。表2计算与试验对比Tab.2Thecomparisonofthesimulationandthetest试件混凝土单轴试件应力状施加在试件水力劈裂水压(MPa)编号龄期抗拉强度态两端的应力试验数值模拟(d)ft(MPa)(MPa)Ⅲ904.1无应力02.72.8Ⅳ283.1单轴压1.03.23.4Ⅴ904.4单轴拉-0.712.32.2-9- 中国科技论文在线http://www.paper.edu.cn3观音岩重力坝水力劈裂数值模拟255观音岩水电站为金沙江水电基地中游河段“一库八级”水电开发方案的最后一个梯级水电站，位于云南省华坪县与四川省攀枝花市的交界处，上游接鲁地拉水电站，下游距攀枝花市27公里，电站水库正常蓄水位1134m，库容约20.72亿立方米。施工过程中，16号坝段上游面出现了表面裂缝，蓄水过程中这些表面裂缝扩展为劈头裂缝。本文以16号坝段为研究对象，在其上游面设置了一个竖直的初始缝面，模拟了其在蓄水过程中的水力劈裂现象。该坝段坝260顶高程1139m，建基面高程990m，坝体及坝体-地基模型如图11所示，模型中考虑了坝体排水、廊道、地基排水、帷幕。由于缺乏施工期形成的表面裂缝的形状、范围，本文假定初始裂缝面形状为三角形，底高程1005m，顶高程1025m，最大深度为5m。计算模型中坝基的计算范围为向上游、下游和基岩深度方向各取1.5倍坝高，坝基底边界为竖向约束，四个侧边界均为水平约束。265计算中考虑的荷载为自重、上下游水荷载，先施加自重荷载，再施加水压力，上游水位按每荷载步增加10m的方式从1020m抬升至1120m，最后一个荷载步水位从1120m抬升至正常蓄水位1134m，下游水位1010m保持不变。坝体混凝土、地基力学特性见表3。渗流场计算-9按照恒定渗流考虑，混凝土渗透系数k0=5×10m/s，耦合系数a0=0.01。裂缝面270图11坝体、坝体-地基模型Fig.11Thedammodelandthedam-foundationmodel表3观音岩重力坝材料参数Tab.3ThemechanicalparametersofGuanyingravitydam材料弹性模量E泊松比ν密度ρ抗拉强度ft断裂能Gf3(GPa)(kg/m)(MPa)(N/m)混凝土32.20.16724502.5200地基25.70.22026001.0200-10- 中国科技论文在线http://www.paper.edu.cn275(a)水位1070m(b)水位1090m(c)水位1100m(d)水位1110m280(e)水位1120m(f)水位1134m图12不同水位，坝体损伤分布Fig.12Thedamagedistributionofthedamatdifferentwaterlevel计算中，蓄水至1070m时，坝体出现损伤，不同水位坝体损伤分布如图13所示。由图13可以看出，随着水位增加，损伤区域从裂缝面尖端逐步向坝体上部、内部、下部扩展。水位285从1070m上升至1110m过程中，损伤区域缓慢扩展，从1110m上升至1120m，损伤区域明显增大，扩展至坝体排水管处，水位由1120m上升至正常蓄水位1134m时，损伤区域显著增大，上部达到高程1060m处，下部达到高程992m处，最大裂缝深度33m，最终裂缝面呈圆弧状。-11- 中国科技论文在线http://www.paper.edu.cn290图13裂缝预估范围示意图Fig.13Schematicdiagramofthepredictedfracturerange设计院预测的裂缝范围如图13所示，本文预测的正常水位下裂缝范围与设计院预测的范围基本吻合，基本上预测出了表面裂缝在正常蓄水位下的扩展范围，可以为评估表面裂缝的危险程度提供参考。然而仍有一定的差异，设计院预测的裂缝范围比本文向内部扩展的更深295一点，向上扩展的更高一点。经过初步分析差异主要与初始裂缝面形状、范围，坝体力学特性参数，蓄水过程等有关，计算中当初始裂缝面的最大深度增加时，正常水位下裂缝扩展范围增大，坝体抗拉强度降低时，正常水位下裂缝扩展范围也增大，所以掌握真实的初始裂缝面形状、范围，坝体力学特性参数，蓄水过程等资料是采用本文耦合模型准确预测裂缝是否扩展、扩展过程、最终扩展范围的关键。3004结论本文提出一种应力-渗流-损伤耦合模型模拟重力坝的三维水力劈裂现象。该耦合模型具有以下特点：(1)考虑了混凝土断裂过程区的应变软化特性；(2)基于断裂能守恒原理将损伤模型与断裂力学相结合，使得断裂能的消散不会受到网格大小的影响；(3)考虑了水力劈裂过程中的四个耦合过程。305采用该耦合模型，模拟了一个内置裂缝的圆柱体试件的水力劈裂，数值结果与试验结果吻合很好，表明该耦合模型可以很好的模拟水力劈裂。在此基础上，采用该耦合模型，模拟了观音岩重力坝表面裂缝蓄水后扩展成劈头裂缝的过程，数值结果表明，随着上游水位的抬升，裂缝逐渐向坝体上部、内部、下部扩展，水位从1120m抬升至正常蓄水位时，损伤区域显著扩大，最终裂缝面呈圆弧状。本文的模拟结果与设计院预估的裂缝范围基本吻合，然310而仍有一定的差异，主要是因为本文预设的初始裂缝面形状、范围，采用的坝体力学特性参数，蓄水过程等与实际情况存在差异。由于损伤模型通过调整材料本构关系来模拟裂缝，不需要重新划分网格，采用本文耦合模型，可以十分便捷有效地模拟重力坝的三维水力劈裂现象，可以预测施工期形成的表面裂缝在蓄水期是否扩展、扩展路径、最终的扩展范围，评估表面裂缝的危险程度以及时采取有315效的处理措施，为高重力坝的安全运行提供技术保障。-12- 中国科技论文在线http://www.paper.edu.cn[参考文献](References)[1]朱伯芳.重力坝的劈头裂缝[J].水力发电学报,1997(4):85-93.[2]梁文浩,谯常忻,涂传林.柘溪大头坝裂缝问题研究述评[J].水利学报,1982(6):3-12.[3]曹泽生,陈重喜.桓仁大头坝上游坝面劈头裂缝稳定性分析[J].水力发电,1991(3):9-14.320[4]曹刚,伞小雨,王宇.玉石水库6号坝段劈头裂缝原因分析[J].水利建设与管理,2008,28(5):11-14.[5]BruhwilerE,SaoumaVE.Waterfractureinteractioninconcrete-PartⅠ:Fractureproperties[J].ACIMaterialsJournal,1995,92(3):296-303.[6]BruhwilerE,SaoumaVE.Waterfractureinteractioninconcrete-PartII:Hydrostaticpressureincracks[J].ACIMaterialsJournal,1995,92(4):383-390.325[7]SlowikV,SaoumaVE.Waterpressureinpropagatingconcretecracks[J].JournalofStructuralEngineering,2000,126(2):235-242.[8]王建敏.静水压力环境下混凝土裂缝扩展与双K断裂参数试验研究[D].大连:大连理工大学,2008.[9]贾金生,汪洋,冯炜,等.重力坝高压水劈裂模拟方法与特高重力坝设计准则初步探讨[J].水利学报,2013,39(2):127-133.330[10]陈胜宏,邹丽春.小湾高拱坝坝踵开裂的有限单元法分析[J].水利学报,2003(1):66-71.[11]BarpiF,ValenteS.Modelingwaterpenetrationatdam-foundationjoint[J].EngineeringFractureMechanics,2008,75(3-4):629-642.[12]董玉文,任青文.重力坝水力劈裂分析的扩展有限元法[J].水利学报,2011,42(11):1361-1367.[13]王克峰,章青,夏晓舟.考虑流固耦合效应的重力坝水力劈裂模拟[J].应用数学和力学,2015,36(9):335970-980.[14]SecchiS,SchreflerBA.Amethodfor3-Dhydraulicfracturingsimulation[J].InternationalJournalofFracture,2012,178:245-258.[15]MohammadnejadT,KhoeiAR.Anextendedfiniteelementmethodforhydraulicfracturepropagationindeformableporousmediawiththecohesivecrackmodel[J].FiniteElementsinAnalysisandDesign,2013,34073(15):77-95.[16]EhlersW,LuoC.Aphase-fieldapproachembeddedintheTheoryofPorousMediaforthedescriptionofdynamichydraulicfracturing[J].ComputerMethodsinAppliedMechanicsandEngineering,2016,315:348-368.[17]OuchiH,KatiyarA,YorkJ,etal.Afullycoupledporousflowandgeomechanicsmodelforfluiddrivencracks:aperidynamicsapproach[J].ComputationalMechanics,2015,55(3):561-576.345[18]PanPZ,RutqvistJ,FengXT,YanF,andJiangQ.Adiscontinuouscellularautomatonmethodformodelingrockfracturepropagationandcoalescenceunderfluidpressurizationwithoutremeshing[J].RockMechanicsandRockEngineering,2014,47(6):2183-2198.[19]ShimizuH,MurataS,IshidaT.Thedistinctelementanalysisforhydraulicfracturinginhardrockconsideringfluidviscosityandparticlesizedistribution[J].InternationalJournalofRockMechanicsandMiningSciences,3502011,48(5):712-727.[20]ZhangGX,LiX,LiHF.Simulationofhydraulicfractureutilizingnumericalmanifoldmethod[J].ScienceChinaTechnologicalSciences,2015,58(9):1542-1557.[21]OliaeiMN,PakA,SogaK.Acoupledhydro-mechanicalanalysisforpredictionofhydraulicfracturepropagationinsaturatedporousmediausingEFGmesh-lessmethod[J].ComputersandGeotechnics,2014,35555(2):254-266.[22]ChernyS,LapinV,EsipovD,etal.Simulatingfully3Dnon-planarevolutionofhydraulicfractures[J].InternationalJournalofFracture,2016,201:181-211.[23]HillerborgA,ModéerM,PeterssonPE.Analysisofcrackformationandcrackgrowthinconcretebymeansoffracturemechanicsandfiniteelements[J].Cementandconcreteresearch,1976,6(6):773-781.360[24]江见鲸,陆新征,叶列平.混凝土结构有限元分析[M].北京:清华大学出版社,2005.[25]J.Lemaitre.Howtousedamagemechanics[J].NuclearEngineeringandDesign,1984,80(2):233-245.[26]张国新.SAPTIS:结构多场仿真与非线性分析软件开发及应用(之一)[J].水利水电技术,2013,44(1):52-55.[27]任青文,刘爽,陈俊鹏,等.基于损伤理论的重力坝坝基岩体渐进破坏数值模拟研究[J].水利学报,3652014,45(1):1-9.[28]BiotMA.Generaltheoryofthree-dimensionalconsolidation[J].JournalofAppliedPhysics,1941,12(2):155-164.[29]BiotMA.Theoryofelasticityandconsolidationforaporousanisotropicsolid[J].JournalofAppliedPhysics,1955,26(2):182-185.370[30]BiotMA.VariationalLagrangian-thermodynamicsofnonisothermalfinitestrainmechanicsofporoussolidsandthermomoleculardiffusion[J].InternationalJournalofSolidsandStructures,1977,13(6):579-597.[31]BhattacharjeeSSandLégerP.Fractureresponseofgravitydamsduetoriseofreservoirelevation[J].JournalofStructuralEngineering,1995,121(9):1298-1305.[32]BaryB,BournazelJP,andBourdarotE.Poro-damageapproachappliedtohydro-fractureanalysisofconcrete375[J].JournalofEngineeringMechanics,2000,126(9):937-943.[33]林凯生,李宗利.高孔隙水压作用下混凝土渗流-损伤耦合模型[J].水利水运工程学报,2010,(2):51-55.[34]周秋景,张国新.SAPTIS:结构多场仿真与非线性分析软件开发及应用(之二)[J].水利水电技术,2013,44(9):39-43.-13- 中国科技论文在线http://www.paper.edu.cn[35]张国新,沙莎.混凝土坝全过程多场耦合仿真分析[J].水利水电技术,2015,46(6):87-93.380-14-'