毕业设计（论文）-混凝土面板堆石坝面板开裂机理及效应研究.doc

摘要 I 论文题目：混凝土面板堆石坝面板开裂机理及效应研究论文题目：混凝土面板堆石坝面板开裂机理及效应研究 学科学科名名称称：水工结构工程：水工结构工程 研研 究究 生：生： 签签 名：名： 指导教师：指导教师： 签签 名：名： 摘摘 要要 混凝土面板堆石坝具有对不同地基条件适应性强、施工快、工期短、投资省等优 点，无论是大型还是中小型水电工程，都成为首选比较坝型，与碾压混凝土坝和高拱 坝并列成为最有发展前景的三类坝型。对混凝土面板堆石坝而言，面板 及其接缝止水 组成的防渗系统是大坝的主要防渗结构，面板的防裂问题是关系 到坝体安全和运行性 状的一个重要课题。我国在面板堆石坝建设 20 多年的时间里，取得了一定的成绩，但 也有一些失败的教训。西北口面板堆石坝的大量结构性裂缝、株树桥面板堆石坝面板 的塌陷以及沟后面板砂砾石坝的垮坝等事故都昭示着我们：进行混凝土面板堆石坝面 板破坏机理及效应研究具有重要意义。 目前，关于混凝土面板堆石坝面板应力变形的分析研究，主要集中在计算其在水 荷载和自重作用下的应力变形，而认为面板尺寸相对较小，将温度应力忽略不计，这 显然与混凝土面板堆石坝尤其是处于日温差和年温差较大的寒冷地区的面板坝面板的 实际应力变形情况不符。国内外关于面板堆石坝面板裂缝成因的研究结果表明：面板 裂缝尤其是贯穿性裂缝的产生主要是由于坝体的不均匀沉降及面板的温度应力共同作 用的结果。针对这方面问题，本文模拟混凝土面板堆石坝的施工过程，建立考虑面板 温度应力、水荷载及自重等情况的应力变形有限元计算模型，全面系统的研究混凝土 面板的应力变形特性，基于断裂力学的基本原理，分析混凝土面板裂缝的产生机理及 开裂过程。 在充分掌握混凝土面板裂缝分布规律的基础上， 预测分析了混凝土面板开裂情况 下堆石坝的渗流及应力变形特性。首先， 针对面板缝隙的物理特性，以等宽缝隙水流 运动的规律为基础，研究建立了面板缝隙渗流 的计算模型。然后，以堆石体的孔隙率 为桥梁，对多孔岩土介质渗透系数与体积应变之间的数学表达式进行了推导，研究建立 了混凝土面板堆石坝渗流场与应力场耦合分析的数学模型、有限元计算模型。最后，应 用本文建立的计算模型，对公伯峡面板堆石坝进行了 渗流场、渗流场与应力场耦合分 析有限元仿真计算，预测分析了假定的面板开裂情况下公伯峡混凝土面板堆石坝的渗 西安理工大学博士学位论文 II 流场及应力变形分布规律，计算结果与实测结果的对比分析表明，本文的计算方法 与 计算模型是合理准确的。 本论文的研究成果，可以方便、快速、准确地计算面板的应力变形分布规律，研究 分析面板裂缝产生的机理及破坏过程，预测分析面板开裂情况下，堆石坝的渗流与应力 变形特性，是进一步研究混凝土面板堆石坝的理论基础。 关键词：关键词：混凝土面板堆石坝；面板裂缝；开裂机理；渗流特性；应力变形特性 Abstract III Title: Research on the Mechanism and Effect of Concrete Face Rockfill Dam Cracking Major: Hydraulic Structure Engineering Name: Li Yanlong Signature： Supervisor: Prof. Li Shouyi Signature： Prof. Wang Ruijun Signature： Abstract CFRD has a lot of advantages such as strong adaptability to the different soil conditions, fast construction, short construction period, less investment and so on, whether large or small and medium hydropower projects, CFRD has become the first choice of dam type comparison. It is parallel to RCC and high arch dam which are the most promising three types of dam. For CFRD, face slab and joint sealing constitute the seepage system of main dam. Face slab cracking problem is an important topic related to the dam safety and operation characteristics. In our country,there are some achievements are obtained in the construction of CFRD in the last 20 years, but also some failures. Many accidents such as large number of structural concrete cracks in the XIBEIKOU concrete face rockfill dam ,the collapse of ZHUSHUQIAO concrete face rockfill dam,and the dam-failure of GOUHOU concrete face rockfill dam are clear to tell us that it has great significances to research on the mechanism and effect of the face-failure in CFRD. At present, the analysis of the stress and deformation in CFRD are mainly in the calculation of the stress and deformation under the loads of its weight and water press, in which case the size of face slab is relatively small, the temperature stress is negligible, that is obviously difffreren from the actual stress and deformation of concrete face rockfill dam, especially in the cold regions where the daily temperature difference and year temperature difference are very large. The results of research on the resons of the face rockfill dam craks at home and abroad showed that face slab cracks, especially penetrating cracks are mainly due to both uneven settlement of the dam and the temperature stress of face slab. Focusing on this problem, simulating the construction process of concrete face rockfill dam in this paper,the stress deformation finite element model is established which has considered temperature stress weight and water press etc. to comprehensively and systematically analysis the stress and deformation properties of concrete face slab. The mechanism of concrete slab cracks and the cracking process are analysed based on the basic principles of fracture mechanics. 西安理工大学博士学位论文 IV The stress and deformation properties and seepage characteristic of rockfill dam are analysed and predicted based on fully grasping the distribution of cracks in concrete slab. First, according to the physical characteristics of face slab gap ,the calculation model of face slab gap flow is established based on the movement law of water flow in the equal-width cracks.Second, use the porosity of rockfill as a bridge and deduce the mathematical expression between permeability coefficient of porous soil medium and volumetric strain. And then the coupling analysis model of seepage/stress and finite element model of CFRD are established. Finally, use the model established to calculate the seepage field and coupling analyse of seepage/stress. The stress deformation distribution law of Gongboxia CFRD is predicted in the case of hypothetical face slab cracks. The comparative analysis of calculated and measured results show that the calculation method and calculation model is reasonable and accurate. With the result of this paper, the stress deformation distribution law can be calculated easily, quickly and accurately. Researching and analyzing the mechanism of the slab cracks and damage process of production, predicting and analyzing the seepage and stress and deformation properties in the case of face slab cracks,which are the theory foundation to the further study in CFRD. Key words: concrete face rockfill dam; face slab crack; cracking mechanism; seepage characteristic; stress and deformation properties 目录 V 目目 录录 1 1 绪论绪论 1 1 1.1 混凝土面板堆石坝发展概况.1 1.2 混凝土面板开裂机理研究进展.2 1.2.1 面板开裂影响因素研究现状.2 1.2.2 混凝土面板堆石坝面板温度应力的研究现状.3 1.2.3 混凝土面板堆石坝面板开裂过程研究现状.4 1.3 混凝土面板堆石坝的渗流及应力变形特性研究进展.6 1.3.1 混凝土面板堆石坝渗流特性研究现状.6 1.3.2 混凝土面板堆石坝应力变形特性研究现状.6 1.3.3 混凝土面板堆石坝渗流场与应力场耦合分析研究现状.7 1.4 本文研究内容和方法.9 1.4.1 本文研究内容.9 1.4.2 预期目标10 1.4.3 研究方法和技术路线10 1.5 研究意义及创新之处11 1.5.1 研究意义11 1.5.2 本论文的创新之处11 2 2 混凝土面板堆石坝温度应力仿真计算研究混凝土面板堆石坝温度应力仿真计算研究 1313 2.1 混凝土水化热温升变化规律13 2.1.1 水泥水化热变化规律13 2.1.2 混凝土的绝热温升变化规律13 2.2 热传导理论14 2.2.1 热传导的基本微分方程14 2.2.2 热传导方程的定解条件15 2.3 混凝土面板温度应力分析的弹性基础梁法16 2.4 混凝土面板温度应力分析的有限元法17 2.4.1 稳定温度场有限元分析17 2.4.2 非稳定温度场有限元分析18 2.4.5 温度应力有限元分析21 2.5 混凝土面板的温度徐变应力23 2.5.1 混凝土材料的徐变23 西安理工大学博士学位论文 VI 2.5.2 混凝土徐变应力分析25 2.6 公伯峡混凝土面板堆石坝面板温度应力仿真计算27 2.6.1 公伯峡混凝土面板堆石坝工程概况27 2.6.2 计算基本条件28 2.6.3 计算工况30 2.6.4 有限元计算模型31 2.6.5 计算结果分析 31 2.6.6 计算结果与观测结果的对比41 2.7 本章小结42 3 3 混凝土面板应力变形特性有限元分析混凝土面板应力变形特性有限元分析 4343 3.1 面板堆石坝应力变形数值计算的材料本构模型43 3.1.1 混凝土面板的本构模型43 3.1.2 散粒体材料的本构模型44 3.1.3 中点增量法45 3.2 混凝土面板与垫层之间的接触面45 3.3 基于 ADINA 软件的面板堆石坝应力变形仿真计算47 3.3.1 ADINA 软件的二次开发 47 3.3.2 ADINA 软件对几个关键问题的模拟 49 3.3.3 ADINA 软件计算流程 50 3.4 公伯峡混凝土面板堆石坝应力变形有限元分析52 3.4.1 有限元计算模型52 3.4.2 计算参数52 3.4.3 计算工况53 3.4.4 计算结果分析53 3.4.5 计算结果与观测结果的对比分析56 3.5 本章小结58 4 4 混凝土面板开裂机理及破坏过程研究混凝土面板开裂机理及破坏过程研究 5959 4.1 混凝土面板结构特征59 4.2 混凝土面板初始微裂缝的形成机理61 4.2.1 面板与垫层结合关系的描述61 4.2.2 面板微裂缝形成机理的定性分析63 4.2.3 面板微裂缝形成机理的定量分析64 4.3 基于断裂力学原理的混凝土面板裂缝破坏过程研究66 4.3.1 断裂力学的基本原理66 目录 VII 4.3.2 应力强度因子的有限元计算方法68 4.3.3 应力强度因子的线性拟合求解方法70 4.3.4 混凝土面板裂缝断裂力学分析模型71 4.3.5 混凝土面板裂缝扩展过程分析71 4.4 本章小结72 5 5 面板开裂情况下堆石坝的渗流计算模型研究面板开裂情况下堆石坝的渗流计算模型研究 7373 5.1 稳定渗流计算的基本理论73 5.1.1 达西定律73 5.1.2 稳定渗流的连续性方程74 5.1.3 稳定渗流基本微分方程75 5.1.4 稳定渗流基本微分方程的定解条件75 5.2 三维稳定渗流有限元计算原理76 5.3 混凝土面板堆石坝面板缝隙渗流计算模型77 5.3.1 面板缝隙稳定流的运动规律77 5.3.2 面板接缝的渗流计算模型78 5.3.3 面板裂缝的渗流计算模型80 5.4 混凝土面板堆石坝渗流有限元计算模型84 5.5 应用工程实例85 5.5.1 计算工况85 5.5.2 有限元计算模型 86 5.5.3 计算参数87 5.5.4 计算结果分析87 5.5.5 与实测结果的对比90 5.6 本章小结91 6 6 面板开裂情况下堆石坝的渗流场与应力场耦合分析研究面板开裂情况下堆石坝的渗流场与应力场耦合分析研究 9292 6.1 岩土多孔介质饱和渗流的基本规律92 6.1.1 非稳定渗流的微分方程92 6.1.2 非稳定渗流基本微分方程的定解条件96 6.1.3 有限元方程97 6.2 混凝土面板堆石坝渗流场与应力场的耦合原理98 6.2.1 渗流场对应力场的影响98 6.2.2 应力场对渗流场的影响.100 6.3 面板堆石坝渗流场与应力场耦合分析的有限元方程.100 6.3.1 应力场有限元方程.100 西安理工大学博士学位论文 VIII 6.3.2 渗流场有限元方程.101 6.3.3 渗流场与应力场耦合分析的有限元方程.102 6.4 面板堆石坝渗流场与应力场耦合分析的步骤.102 6.5 应用实例.103 6.5.1 计算模型.103 6.5.2 计算参数.104 6.5.3 计算结果分析.105 6.5.3 计算结果分析.105 6.6 本章小结.108 7 7 结论与展望结论与展望 109109 7.1 结论 .109 7.2 展望.110 致致 谢谢 111111 参考文献参考文献 112112 附附 录录 118118 1 绪论 1 1 1 绪论绪论 1.11.1 混凝土面板堆石坝发展概况混凝土面板堆石坝发展概况 混凝土面板堆石坝（Concrete Face Rockfill Dam）是以堆石分层碾压填筑成坝体， 在上游面布置混凝土面板作为防渗体的一种土石坝【1】。混凝土面板堆石坝具有安全性高、 造价较低、适用性强、施工方便等优点，因此，在实际工程中得到了广泛应用和发展。 美国的莫拉(Morena)坝是世界上第一座面板堆石坝，坝高 54m，修建于 1895 年，距 今已有 100 多年的历史。著名的混凝土面板堆石坝工程专家库克（J.B.Cook） 【2，3】认为， 混凝土面板堆石坝的发展主要经历了三个阶段，分别为：1850 年1940 年的早期抛填堆 石阶段，1940 年1965 年的抛填堆石向碾压堆石过渡阶段以及 1965 年以后的堆石薄层 碾压为特征的现代阶段【4】。早期抛填堆石阶段：1850 年1900 年正值美国的淘金热， 在美国的加利福尼亚州由于淘金的需要修建蓄水库，最初的面板堆石坝是用抛填块石作 为坝体，表面覆盖木板作为防渗面板，直到十八世纪末，随着筑坝技术的进步，才开始 修建混凝土面板堆石坝。1931 年建成的美国盐泉坝（高 100m）及 1955 年建成的葡萄牙 帕拉德拉坝(高 112m)，混凝土面板难以适应堆石体较大的变形，水库蓄水后，面板出现 严重的开裂和大量的漏水，因此，这一坝型的发展在 20 世纪 4050 年代处于停滞状态。 过渡阶段：1960 年美国土木工程师协会发表了美国一次堆石坝学术会议的论文集，在 C.M.Roberts 的讨论文章【5】中，介绍了最早采用薄层碾压堆石坝的工程实例即 1958 年 完成的 Quoich 坝。土力学家太沙基提出，堆石体经过碾压后，变形就会很小，混凝土面 板的应力变形状况就会大大改善，也就可以修建更高的面板堆石坝。同时，施工机械的 发展，尤其是大型振动碾的发明，为混凝土面板堆石坝的再次兴起起到了推进作用。从 此，堆石体的填筑施工均采用薄层碾压堆石代替抛填堆石【68】。以堆石薄层碾压为特征 的现代混凝土面板堆石坝阶段：1965 年以后，由于堆石体采用分层填筑碾压的施工方法， 使得密实度和变形模量大大提高，因此坝体的变形也大为减小。1971 年建成的澳大利亚 Cethena 坝（高 110m） ，奠定了现代混凝土面板堆石坝的技术基础。此后，在学者、专家 及工程技术人员的共同努力下，混凝土面板堆石坝在世界范围内得到了迅速发展，筑坝 技术也越来越成熟。从目前已建的工程看，大多数混凝土面板堆石坝运行状况良好，有 的工程还成为了坝工史上的经典力作。 我国以现代技术修建混凝土面板堆石坝始于 20 世纪 80 年代初期【1】。1985 年，我 国第一座开工建设的混凝土面板堆石坝是坝高为 95m 的湖北西北口大坝，第一座完工的 是辽宁关门山水库大坝，坝高为 58.5m。中国的混凝土面板堆石坝建设与国外相比，虽 然起步较晚，但起点高，发展快。目前，我国的混凝土面板堆石坝无论是在建坝高度、 建坝数量，还是在筑坝技术方面，都处于世界领先地位。据统计【9】，截止 2004 年我国 的混凝土面板堆石坝有 150 多座，其中，坝高超过 100m 的有 37 座。我国的清江水布娅 面板坝，是目前世界上最高的混凝土面板堆石坝，坝高为 233.2m。 西安理工大学博士学位论文 2 与其它坝型相比，混凝土面板堆石坝在筑坝技术方面、适用性方面以及经济方面均 拥有独特的优势【10】，在国内外获得了广泛的推广和发展，设计理论和施工技术水平也 不断得提高，因此，工程规模越来越大，坝越建越高是混凝土面板堆石坝发展的基本趋 势【11】。 1.21.2 混凝土面板开裂机理研究进展混凝土面板开裂机理研究进展 1.2.11.2.1 面板开裂影响因素研究现状面板开裂影响因素研究现状 混凝土面板是以堆石体为支撑结构，它承担着向下游堆石体传递水压力及大坝防渗 的重任【12】，对大坝的安全与工程的运行起着决定性的作用。从空间结构看，面板的长 度较长，宽度较小，而厚度相对来说很小，是一块长条形的薄板，因此，受环境温度的 影响较大，尤其是高坝及处于高寒地区的混凝土面板堆石坝【13】，在气温骤变条件下， 面板内将产生拉应变，再加上坝体的沉降及其它不确定性因素，面板很容易出现裂缝。 从国内外修建的面板堆石坝看，面板开裂的情况比较普遍，个别的坝还相当严重【1417】， 因此，面板的抗裂问题是面板堆石坝的重要课题。国内外研究表明：混凝面板的温度应 力、外界荷载（主要为水荷载）以及堆石体变形是面板开裂的主导因素。 a.a. 混凝土力学特性及荷载对面板开裂的影响混凝土力学特性及荷载对面板开裂的影响 混凝土面板开裂问题，实质上是面板中的破坏力与抵抗力之间的较量。破坏力主要 指的是由外因及内因引起的面板应力，抵抗力则指的是面板混凝土的强度。面板应力主 要受混凝土的力学性能、面板的自重、外部荷载以及坝高、河谷形状等因素的影响。混 凝土力学特性及荷载对面板开裂的影响主要体现在以下几个方面： （1）混凝土材料内部物理化学作用引起的面板应力。混凝土中的碱骨料反应、钢筋 锈蚀、混凝土的徐变以及混凝土遭受冻融循环等因素引起的面板应力。 （2）混凝土面板的温度应力、干缩应力。面板的温度、湿度变化会产生引起结构变 形，由于面板不能自由变形，则会产生应力。面板的温度、干缩应力主要受环境气温条 件、面板的保温保湿条件、极限拉伸及线膨胀系数等影响。 （3）水荷载及自重荷载引起的结构应力。施工期，面板的挠度变形主要由自重引起， 因此，其值较小。运行期，面板挠度变形的主要影响因素为库水压力，水压力越大，混 凝土面板挠度相应越大。同时，坝高、混凝土面板的结构尺寸、堆石的材料特性及河谷 形状对面板的应力变形也有影响。 （4）作用于面板上的库水结冰荷载、地震作用等也会引起面板的开裂破坏。 （5）施工过程及质量。工程实践表明【13】，面板混凝土的施工质量与养护，坝体材 料分区、填筑次序与压实质量、水库蓄水与泄水过程、垫层坡面的保护方式等是影响面 板开裂的因素。 此外，在周边缝位置，如果面板相对于趾板的变形（张开及剪切变形）过大会破坏 止水结构，也会引起面板的开裂。周边缝的变形主要受水荷载、面板变形、坝体变形、 1 绪论 3 趾板的几何形状以及河谷坝肩的形状等因素的影响。 b.b. 堆石体变形对面板开裂的影响堆石体变形对面板开裂的影响 混凝土面板以堆石体为支撑结构，因此，面板的应力变形除了受混凝土力学特性及 荷载对面板开裂的影响外，还受堆石体变形的影响。尤其是运行期，在水压力的作用下， 面板开裂与否主要由堆石体变形决定，因此，堆石体变形的影响因素亦是混凝土面板变 形特性的影响因素。大量监测数据及数值计算结果表明【18,19】，在堆石体的流变、湿化 变形以及库水位升降的影响下，堆石体的变形呈现出随时间而变化的趋势，则混凝土面 板的应力变形也随之变化。 混凝土面板堆石坝的坝体是分层碾压填筑的，在自重及碾压力的作用下，堆石体会 产生竖直和水平位移，浇筑在垫层上的面板变形如果不能适应坝体的变形，就会出现垫 层料脱空，此时，面板失去了坝体的支撑作用，在水压力的作用下，其应力变形特性会 发生变化，主要表现为：脱空部位的面板被压向垫层料，在水压力及垫层反力的作用下， 内部产生较大的应力，当应力超过混凝土的强度则会产生裂缝，发生开裂破坏。 综上所述，混凝土面板发生裂缝的原因主要是堆石体变形和面板的温度应力两个方 面。由温度变化引起的混凝土收缩，加之基础约束将导致面板出现大量密集的有规律的 水平裂缝。由坝体的不均匀变形引起的拉应力，将导致面板出现较大的裂缝，甚至出现 断裂。 1.2.21.2.2 混凝土面板堆石坝面板温度应力的研究现状混凝土面板堆石坝面板温度应力的研究现状 国内外关于混凝土面板堆石坝面板温度应力的计算方法主要有解析计算法（简化条 件下）和有限元数值计算法【20】，相比较而言，有限元法可以考虑混凝土的水化热温升、 徐变、弹性模量、自生体积变形等参数随龄期的变化，考虑气温、水温、自重、水压、 寒潮、等诸多因素的影响，计算结果更为准确可靠。 有限元数值计算方法分为不考虑接触面影响的有限元法和考虑接触面影响的有限元 法。不考虑接触面影响的有限元法，是将面板与垫层之间的接触面按弹性约束处理，不 考虑接触面非线性变形特性对面板温度应力的影响，在混凝土面板堆石坝温度场的计算 结果上，计算面板的温度应力。面板与垫层之间的关系按“弹性约束”处理，夸大了垫 层对面板的约束作用，因此计算结果的准确性不够准确。考虑接触面影响的有限元法， 是将面板与垫层之间的接触面按摩擦约束处理【21】，假定面板与垫层之间为点面接 触，在法向面板与垫层之间用弹簧联接，切向则为摩擦接触。文献【21】根据摩尔库 仑准则，采用嵌接处罚法建立了面板与垫层之间接触面的附加刚度和附加摩擦力荷载矩 阵，模拟了接触面锁定和摩擦两种接触状态。显然，文献【21】研究建立的力学模型考 虑了接触面非线性变形特性及不连续变形特性对面板温度应力的影响，但也存在明显的 不足，如：该计算模型没有考虑面板与垫层之间脱开的自由状态；对于面板与垫层之间 的导热关系，该模型没有明确表述等。混凝土面板堆石坝面板温度应力的计算模型及方 法的研究还处于初步探索阶段，仍有大量科研难题值得深入研究。 西安理工大学博士学位论文 4 面板接缝（周边缝、垂直缝与水平缝的统称）止水结构的模拟是混凝土面板温度应 力研究的关键技术问题之一，尤其是对于高坝，面板的变形会严重影响接缝缝间止水的 工作性态。因此，在面板堆石坝的有限元数值仿真计算中，必须考虑面板接缝的非连续 变形。目前，模拟面板接缝一般采用以下几种模型： 节理材料模型：受压时材料的模量与面板混凝土的模量相同或接近，接缝一定程度 上可以传递压力和剪力、有应变和位移；受拉时，将模量取为一个极小值，近似认为材 料不承受拉应力。 复合模型：将面板与接缝简化为复合材料，在复合材料的刚度中考虑接缝的特性， 将其均化在复合材料之中，把面板处理为各向异性板。 经验模型：按照试验获得的数据提出数学拟合表达式，将其作为接缝的本构模式进 行分析。 上述几种模型虽然有所不同，但是均不同程度地考虑了接缝可以受压难以受拉的变 形非连续特性。 1.2.31.2.3 混凝土面板堆石坝面板开裂过程研究现状混凝土面板堆石坝面板开裂过程研究现状 断裂力学是利用力学的成就，研究含缺陷材料和结构的破坏问题。和其他科学一样， 断裂力学是在生产实践中产生和发展的，是从 20 世纪 70 年代发展起来的一门新兴学科。 尽管起步较晚，但实验与理论均发展迅速，在实际工程中得到了广泛应用。 a.a. 混凝土的断裂过程混凝土的断裂过程 混凝土材料的拉伸、压缩和弯曲等力学性能是混凝土结构设计的重要依据，而试验 则是揭示混凝土力学性能的主要方法。从国内外的试验资料可知：混凝土是典型的非均 匀材料，其中含有夹碴、气泡、孔穴等，其变形、强度以及破坏过程均与裂纹的扩展有 关。混凝土的裂纹可以分为两大类【22】：一类是随机分布的微裂纹，微裂纹的存在决定 了混凝土的抗压与抗拉强度；另一类是宏观裂纹，它使得混凝土的力学性质呈现各向异 性。混凝土的破坏是由混凝土的微裂纹与宏观裂纹共同作用，一方面微裂纹的存在影响 宏观裂纹的萌生，另一方面微裂纹对主裂纹产生屏蔽和劣化的双重作用【23】。混凝土的 破坏过程实质上是其中微裂纹萌生、扩展、贯通，直至微裂纹变成宏观裂纹，进而导致 混凝土失稳破裂的过程。 b.b. 混凝土断裂力学的研究现状混凝土断裂力学的研究现状 1920 年，Griffith 在研究玻璃的强度时首次提出了断裂力学理论，此后，经过 Irin 和 Owen 的补充和完善，在 20 世纪 50 年代适用于高强度钢材的线弹性断裂力学体系基 本形成。1961 年，Kaplan 首先将断裂力学的概念引用到混凝土中，用断裂力学的理论研 究混凝土的破坏机理，并进行了混凝土的断裂韧度实验。此后，国内外诸多学者专家在 混凝土断裂试验方面做了许多研究，积累了大量的实测资料，并提出了一系列应力强度 因子的计算方法和经验断裂判据。 根据受力情况的不同，裂缝可以分为三种基本类型：（1）张开型裂缝，作用力与裂 1 绪论 5 缝面及裂缝前缘均垂直；（2）滑开型裂缝，作用力与裂缝面平行，而与裂缝前缘垂直； （3）撕开型裂缝，作用力与裂缝面及裂缝前缘均平行。图 1-1 为三种类型的裂缝示意图。 （1）张开型裂缝 （2）滑开型裂缝 （3）撕开型裂缝 （1) openning mode crack （2）Slide mode crack （3）Tear mode crack 图 1-1 三种类型的裂缝 Fig.1-1 The three types of crack 断裂力学中，能量释放率 G 是衡量裂缝扩展的重要指标，是指扩展单位裂缝表面积 上裂缝体释放的应变能。能量释放率 G 的表达式为： （1.1）lUllU lB G l 1 lim 0 式中：B 平行裂缝前缘且平行于裂缝面的裂缝宽度； l 裂缝扩展的长度； U 裂缝体释放的应变能； l 垂直于裂缝前缘且平行于裂缝面的长度。 在断裂力学中，还有一个表示抗断裂性能的指标，即裂缝尖端的应力强度因子 K， 它代表了裂缝的应力场强度。假设 为裂缝位置上按无裂缝计算的应力，a 为裂缝的深 度，引入一个表征裂缝大小、位置等的形状系数 Y，应力强度因子可以用下式表示： （1.2）aYK 对于张开型裂缝，裂缝尖端应力强度因子 K 与能量释放率 G 之间存在如下关系： （1.3） E K G 1 1 其中，E = E（平面应力情况)；（平面应变情况)。 2 1 E E 滑开型裂缝同张开型裂缝是相同的。 对于撕开型裂缝： （1.4） E K G 2 1 式中，v 为泊松比；E 为杨氏弹性模量。 c.c. 混凝土面板断裂损伤研究现状混凝土面板断裂损伤研究现状 目前，关于混凝土结构的断裂损伤研究主要集中在大体积混凝土结构、碾压混凝土 西安理工大学博士学位论文 6 拱坝、混凝土重力坝方面，而对于面板堆石坝混凝土面板断裂方面的研究甚少。混凝土 面板对环境温度的变化较敏感，再加上坝体的沉降，面板很容易出现裂缝。因此，基于 断裂力学研究混凝土面板微裂纹萌生、扩展、贯通，直到产生宏观裂纹的破坏机理，对 于面板的抗裂，大坝的安全运行具有重要意义。 1.31.3 混凝土面板堆石坝的渗流及应力变形特性研究进展混凝土面板堆石坝的渗流及应力变形特性研究进展 1.3.11.3.1 混凝土面板堆石坝渗流特性研究现状混凝土面板堆石坝渗流特性研究现状 混凝土面板堆石坝以堆石体做为支撑结构，安全性较高，许多学者认为混凝土面板 堆石坝不会发生渗透破坏，因此，对面板坝的渗流分析研究不够重视，甚至在设计时不 作为计算的主要内容。混凝土面板是大坝防渗的主要结构，一旦发生开裂破坏，将直接 危及到工程的安全与经济效益，应当引起足够的重视【24】。混凝土面板堆石坝的渗流分 析对于工程的设计、运行及管理具有重要的意义【25】。 关于混凝土面板堆石坝的渗流计算模型或计算方法，国内外研究甚少，澳大利亚的 伦吉卡希纳得通过工程实例渗漏观测数据的分析，利用面板的结构尺寸、面板的坡率 及水深与面板裂缝渗漏量的关系，获得了半经验性的化引表面渗透系数，即表面渗透性 法【1】 【25】，该方法将面板作为各向同性的均质体，且面板裂缝中的渗流符合达西定律。 表面渗透性法忽略了垫层、过渡层等面板基础的渗流，只能估算面板裂缝及接缝的渗漏 量，不能确定渗透坡降、渗流速度等其他要素，更不能用于整个大坝的渗流仿真计算。 张嘎，张建民等通过对混凝土面板堆石坝渗流形态的分析【26】，假定裂缝渗流为层流、 面板缝隙为等宽缝隙，基于达西渗流定律进行了面板出现裂缝工况下堆石坝的渗流有限 元计算，该方法没有考虑面板裂缝宽度沿厚度方向变化的特性及裂缝面的粗糙度，因此 计算结果也不尽如人意。正如文献【1】所分析指出的，现有的混凝土面板堆石坝渗流分 析的计算模型都不成熟，仍需要进行更深入的研究，尤其是面板开裂情况下堆石坝的渗 流机理研究。 1.3.21.3.2 混凝土面板堆石坝应力变形特性研究现状混凝土面板堆石坝应力变形特性研究现状 现代混凝土面板堆石坝采用分层碾压的施工方法，堆石体的密实度与变形模量与早 期面板堆石坝相比有了较大的提高，但对于面板坝来说，大坝的应力变形性状仍是制约 其发展的瓶颈问题【2729】。自 20 世纪 60 年代后期，有限单元法开始应用于土石坝的应 力变形分析已来，至今已有五十多年的时间，在此期间，国内外诸多学者专家在这一方 面做了大量的科研工作，取得了丰硕的研究成果【3033】。 1970 年，N.O.Broughton 以 Wilmot 坝（坝高 36m）为例，建立了土石坝应力变形有 限元计算模型，采用的是等应变三角形单元，基于弹性有限元法进行了混凝土面板堆石 坝应力变形数值计算，由于他没考虑面板和堆石体界面的滑动，因此，其计算结果与实 际情况相差较大。此后，J.K.Wilkings 等人对 Alto Anchicaya, Foz do Areia 和 Centhana 坝均进行了应力变形有限元计算，数值计算结果与实测结果的对比表明：计算 1 绪论 7 的下游坝体大部分范围内的竖向位移小于水平位移，与实测结果不一致；下游坝脚附近 坝体位移数值计算结果与实测结果基本相符；数值计算的面板拉应力区比实测结果大。 这种现象主要是因为坝坡位移受材料的各向异性、泊松比的取值等因素影响。Inaizumi H.和 SardinhaA E.等人研究发现：当材料的泊松比取值趋近于零时，计算值与实测值符 合程度较高。 混凝土面板堆石坝应力变形数值分析的研究主要集中在堆石材料、面板与垫层接触 面、接缝(周边缝、面板缝)等关键因素本构模型的建立，以及施工、运行过程的仿真模 拟方面。堆石体的变形特性主要表现在四个方面【34，35】：压硬性，堆石体的变形模量随 着围压的增加而增加；非线性，堆石体的应力应变本构关系呈现典型的材料非线性特性； 应力引起的各向异性，不同方向应力差异引起变形刚度不同；剪缩性，堆石体受荷时， 颗粒会产生破碎和滑移，引起体积收缩。国内外许多学者致力于坝体材料的数学、力学 模型研究，目前，面板堆石坝应力变形数值分析中堆石料采用的的本构模型主要有清华 K-G 模型、剑桥模型、 “南水”双屈服面模型、邓肯-张 E-B 本构模型及邓肯-张 E- 本 构模型。堆石体的本构模型可以分为弹塑性模型和非线性模型两大类，弹塑性模型以 “南水”双屈服面模型为代表；非线性模型以邓肯-张 E- 模型和邓肯-张 E-B 模型为代 表。由于邓肯-张 E-B 本构模型能反映堆石体的主要变形特性，而且参数测定简单，应用 方便，所以得到较为广泛的应用。但邓肯-张 E-B 本构模型只能反映堆石料的非线性和压 硬性，不能考虑堆石体的剪缩性，因此计算所得的水平位移比实际值偏大【36】。 近年来，随着电子计算机的发展，混凝土面板堆石坝应力变形有限元仿真计算得到 了较快发展。国内外的公司及个人研究开发出了大量有限元通用软件与专用计算程序， 考虑的因素越来越多，仿真程度也越来越高，为混凝土面板堆石坝的发展增添了新的活 力。利用计算机强大的计算功能结有限单元法的优势，是当今混凝土面板堆石坝的应力 变形数值仿真计算的主要技术手段。 1.3.31.3.3 混凝土面板堆石坝渗流场与应力场耦合分析研究现状混凝土面板堆石坝渗流场与应力场耦合分析研究现状 目前，对面板堆石坝的渗流场与应力场耦合研究不够重视，渗流场与应力场的计算 是单独进行的，两场分开考虑。应力场方面，已充分考虑了堆石体材料的非线性、混凝 土面板与垫层之间的接触、面板接缝单元等因素，但基本不考虑考虑渗透体积力对坝体 应力场的影响；渗流场方面，基于达西定律，运用有限元法计算确定面板堆石坝坝体浸 润线，进而确定渗流量和渗流水力坡降，没有考虑应力场对堆石体孔隙率变化的影响。 对混凝土面板堆石坝而言，渗流场与应力场的耦合现象是客观存在的，二者之间相 互影响、相互作用。一方面，渗流场的变化引起了渗透体积力的改变，即作用于面板堆 石坝的外荷载发生了变化，进而改变面板堆石坝的应力场；另一方面，运行期的水荷载 会使堆石体变得更加密实，孔隙率减少，改变堆石体的渗透系数，进而改变面板堆石坝 的渗流场。 关于渗流场与应力场的耦合分析的研究成果较多地偏重于裂隙岩体的渗流场与应力 西安理工大学博士学位论文 8 的场耦合问题以及土体渗流场与应力场耦合问题，关于面板堆石坝的渗流场与应力场耦 合问题的研究鲜见报导。 1941年，Biot【37】基于达西渗流定律，假定材料为线弹性且各向同性，把孔隙流体 力学应用到线弹性的应力变形计算中，建立了固结计算的基本理论。 Sandhu和Wilson【38】将土体与孔隙水作为不可压缩固体，应用变分原理，推导建立 了Biot方程。 O.C.Zienkeiwicz和Shiomi【39】基于Biot建立的三维固结理论，研究提出了广义Biot 公式，该公式考虑了材料的非线性与几何非线性。 Savage和Bradock【40】不考虑多孔介质的流固耦合效应，对于多孔介质的不同应力应 变本构关系，建立了各向同性弹性多孔介质的三维固结计算模型。 我国诸多学者在渗流场与应力场耦合分析的研究方面做了大量的科研工作，取得了 一批有价值的成果，如： 柴军瑞【41】从土体的渗透性出发，研究建立了均质土坝渗流场与应力场耦合分析的 连续介质计算模型。 陈晓平等【42】以下式作为渗透系数的计算表达式： （1.5） v v n n aKK 0 0 0 1 exp 采用间接耦合分析的方法，实现了非均质土坝渗流场和应力场耦合分析。 李锡夔，范益群等【43】采用下式计算考虑饱和度及孔隙比影响的各向同性渗透率： （1.6） e w r rw ww S SS bK 3 1 研究建立了非饱和土应力变形和渗流场的有限元计算模型。 杨林德【44】假定饱和土体为均质、连续、各向同性的弹塑性多孔介质，有效应力系 数为二阶张量，根据虚位移原理，建立了各向异性饱和土体渗流场和应力场耦合分析的 数学模型。 李培超【45】以孔隙度和渗透率之间的关系作为耦合分析的动态模型，将多孔介质和 流体作为可压缩体，基于渗流连续性方程，建立了饱和多孔介质流固耦合渗流计算模型。 王媛【46】基于固结理论，将结点位移和孔隙水压力作为未知量，实现了渗流场与应 力场的耦合分析。 李筱艳等【47】通过抽水试验观测数据的分析，建立了土体渗透系数与有效应力相关 的动态模型： （1.7）exp 0 kk 对深基坑土体渗流场与应力场的耦合分析问题进行了研究。 柳厚祥，李宁等【48】根据尾矿坝的勘察试验数据，提出了尾矿坝的渗透系数与初始 1 绪论 9 渗透系数、应力及中值粒径的经验关系式，进行了尾矿坝应力场与渗流场耦合的非稳定 渗流分析。 曾海容、宋惠珍等【49】基于多相渗流理论及连续介质力学理论，建立了解耦形式的 固相变形方程，采用改进的初应力增量加载方法对非线性固体变形问题进行求解。 平扬【50】通过比奥固结理论在弹塑性分析的扩展应用研究，建立了深基坑开挖过程 中渗流场与应力场耦合分析的有限元计算模型。 徐则民【51】对斜坡稳定性中的渗流场与应力场耦合分析问题进行了研究，提出了斜 坡稳定性评价中的耦合分析的技术路线。 罗晓辉【52】通过对深基坑开挖过程中渗流场及应力场耦合分析研究，提出了稳定渗 流场与非稳定渗流渗流场有限元计算的方法。 陈波、李宁等【53】基于多孔介质温度场、应力场、渗流场三场耦合分析数学模型， 从微分控制方程出发，推导建立了固液两相介质的三场耦合分析的有限元计算方程。 陈庆中、冯星梅等【54】在Sandhu和钱伟长的研究成果基础上，提出了渗流场、应力 场及流场耦合分析的方法。 1.41.4 本文研究内容和方法本文研究内容和方法 1.4.11.4.1 本文研究内容本文研究内容 本文的研究工作主要是针对混凝土面板堆石坝面板开裂机理及效应开展的，其主要 研究内容如下： （1）混凝土面板堆石坝温度应力仿真计算研究 混凝土面板的温度应力是面板早期初始裂缝产生的主导因素之一，本文基于公伯峡 混凝土面板堆石温度场、温度应力的仿真计算结果，研究分析了混凝土面板施工期、运 行期及遭遇寒潮 3 种工况下的温度场、温度应力分布规律。 （2）混凝土面板应力变形特性有限元分析 混凝土面板是以堆石体为基础的，其应力变形特性除了受自身因素（自重、温度应 力、徐变）和