小湾高拱坝坝踵开裂的有限单元法分析

1、小湾高拱坝坝踵开裂的有限单元法分析摘要：在开裂单元中引入不连续形函数，然后基于虚功原理推导出开裂单元的广义刚度矩阵和广义荷载的详细表达式，以及有限元平衡方程，从而建立了一种三维裂缝扩展的不变网格有限元分析方法。运用该方法对已建的二滩拱坝和待建的小湾工程拱坝在一样的工况下进展了比照计算。计算结果说明：在目前所考虑的影响因素下，小湾拱坝和二滩拱坝的开裂范围、开裂程度及对网格疏密、混凝土抗拉强度和水力劈裂等影响因素的敏感程度都非常相似。考虑到二滩工程已建成，且运行正常，可以初步断定小湾工程的拱坝坝踵开裂危险性不大。关键词：拱坝开裂有限单元法坝踵开裂是高拱坝设计中需要考虑的一个重要问题。奥地利的科恩布

2、赖拱坝建于1977年，坝高200.2，在蓄水前坝下游底部程度施工缝张开，1978年蓄水接近满库时上游坝踵开裂穿通廊道，引起廊道内漏水，后来进展了十几年的漫长修补加固，并在坝的下游坝趾处建造了高约70的混凝土支撑体。俄罗斯的萨扬舒申斯克拱坝建于1987年，坝高242.0，施工期产生大量裂缝，当1990年水位首次到达正常水位时，河床坝段上游坝面靠近坝基处产生程度裂缝，引起廊道内漏水。这些高拱坝的损坏事故引起了世界各国坝工界的极大关注。研究拱坝开裂分析的主要方法可以分为构造模型试验和数值计算两大类。可用于拱坝开裂分析的数值计算方法很多，如拱梁分载法、边界元法、流形元法、无单元法、有限单元法等17，其

3、中有限元法是较为成熟且广为应用的数值方法。但有限元法在本质上是一种连续介质力学方法，必须加以改良才能被运用来分析裂缝扩展这样的不连续问题。改良的方法可以概括为两类：变网格法和固定网格法。变网格法随着裂缝的扩展而重新划分网格8，9，这种方法的优点是不需要开发新的有限元软件，但是存在计算量大和前处理难等问题。固定网格法那么是保持有限元网格不变，通过修改开裂单元的插值关系和本构关系来反映裂缝的影响1012。相比拟而言，固定网格法的应用更为方便。本文在开裂单元中引入不连续形函数，然后基于虚功原理推导出开裂单元的广义刚度矩阵和广义荷载的详细表达式，以及有限元平衡方程，从而建立了一种三维裂缝扩展的不变网格

4、有限元分析方法。本文方法在小湾工程拱坝的分析中获得了应用。由于拱坝开裂的影响因素复杂，计算方法也多种多样，很难做出完全定量的分析评判，且目前也没有相应的开裂深度平安准那么，故研究中采用“校准法的思路：取一个与小湾工程相类似的已建且运转正常的工程，在一样的工况下，采用同一计算方法进展比照分析，通过开裂深度与坝底宽度比值的比拟，评价小湾高拱坝的开裂危险性和平安度。根据设计部门的建议，作为校准的已建工程选为二滩工程。1三维裂缝扩展的不变网格有限单元分析方法1.1开裂单元的位移形式12常规有限单元法中，单元位移场u可以表示为：u=Nue(1)图1开裂单元示意式中：N和ue分别为单元形函数矩阵和单元结点

5、位移向量。式(1)表示的位移场在单元内是连续的，不能描绘单元开裂以后沿开裂面的不连续位移形式。为求改良，需构造不连续位移场插值函数。设单元e被开裂面分为e+和e-两局部(图1)，它们的位移场分别记为u+和u-，可以利用原单元e的形函数矩阵表示为：u+=Naeu-=Nbe(2)式中：ae和be是单元结点的广义位移向量。令：(3)那么：ae=e+de；be=e-de(4)将式(4)代入式(2)可得：u+=Ne+Nde；u-=Ne-Nde(5)或：u=Ne+NH(x)de(6)式中：NH(x)称为不连续形函数矩阵，其中H(x)是一个不连续函数，定义如下：(7)式(6)右端的第一项与式(1)的右端项在

6、形式上是一致的，表示的是单元中的连续位移场，而第二项表示的那么是由于裂缝的出现而附加的不连续位移常根据式(6)，开裂单元的位移形式是在连续单元位移形式的根底上，增加了反映不连续位移的插值项。这些附加的插值项与单元结点对应，但其形函数不连续，故称之为开裂单元的广义结点。开裂单元的应变可以表示为：=Be+BH(x)de(8)式中：B是应变矩阵。应力可以表示为：=D=DBe+DBH(x)de(9)式中：D是弹性矩阵，它是单元材料的弹性模量和泊松比的函数。引入不连续形函数和广义结点的概念后，开裂单元的位移形式和应变形式都可以看作是在连续单元的形式上增加了广义结点的影响项，从而使得固定网格下的裂缝扩展分

7、析变得简便易行。1.2有限元平衡方程不连续单元的虚功原理可以表示为：(10)式中：v和s分别表示单元的体积和边界(包括裂缝面)；qv和qs分别是单元的体积力和边界面力(包括裂缝面上的分布力)。将式(6)、式(8)和式(9)代入式(10)，整理可得：Kue=fe(11)(12)(13)(14)式中：K称为单元广义刚度矩阵，ue称为单元结点广义位移向量，fe称为单元结点广义荷载向量。K由四种类型的子矩阵组成，它们的表达式分别为：(15)(16)fe由两局部组成，分别为：(17)(18)其中式(15)和式(17)在形式上与常规有限元法是一致的。考虑型裂缝的情况，且直接采用裂缝前沿附近单元高斯点的应力

8、，根据最大拉应力准那么判断裂缝是否扩展，并且假定开裂单元被裂缝面贯穿。开裂分析中裂缝的扩展将产生应力释放，可以按下式计算其等效荷载：(19)式中的积分沿裂缝面进展。1.3水力劈裂作用的模拟在拱坝开裂计算中终究如何考虑水力劈裂的作用是一个复杂的问题。随着裂缝的扩展，水流渗入裂缝，在裂缝面上产生劈裂作用，同时也引起渗流场的改变。假如考虑渗流场与裂缝扩展的耦合作用，将会使计算过程变得非常复杂，甚至可能无法实现。因此，宜对水力劈裂的作用进展适当的简化。对于高拱坝坝踵的张裂缝，可以不考虑水流渗入裂缝引起的渗流场改变，而只考虑水流在裂缝面的劈裂作用，作用力大小近似取上游全水头压强。这样简化是偏于平安的。以

9、上算法已在有限单元法程序RE3实现，并且通过算例考核证实了其合理性12。转贴于论文联盟.ll.2小湾和二滩拱坝应用实例论文联盟.LL.编辑。2.1工程概况小湾水电站位于云南省大理州的南涧县和临沧地区凤庆县的交界，是澜沧江中下游河段规划8个梯级电站中的第二级，电站装机4200，水库库容151.321083。大坝为混凝土抛物线双曲拱坝，最大坝高292，坝顶高程1245，是即将开工的世界第一高拱坝。小湾拱坝拱冠梁顶部宽度12，拱冠梁底部宽度72.91，厚高比为0.250。二滩水电站位于四川省攀枝花市境内，雅砻江下游，电站装机3300，年发电170亿度。二滩大坝为混凝土抛物线双曲拱坝，最大坝高240，

10、坝顶高程1025，是已建的中国第一高拱坝。二滩拱坝拱冠梁顶部宽度11，拱冠梁底部宽度55.74，厚高比为0.232。2.2计算条件高拱坝在坝踵、坝趾和根底高应力区会出现应力集中，并在坝踵和局部根底内产生裂缝。随着裂缝的开展，库水会进入裂缝，产生水力劈裂作用，从而促使裂缝进一步扩展。假如坝体与地基连接很好，接触面具有很高的强度，开裂可能会深化地基，但是这种裂缝往往属于稳定扩展，对大坝不会构成太大的危害。一般情况下，坝体与地基的交接面附近为一弱面，此面在靠近上游坝踵处往往受一法向拉应力，从而导致裂缝沿坝体与地基的交接面或稍高处扩展，此种裂缝对拱坝平安的危胁最大。本研究将以沿坝体与地基交接面或稍高部

11、位的裂缝扩展作为重点。2.2.1计算模型计算模型的建立对三维裂缝问题的求解有显著的影响，网格划分密实，单元形状相对正态，可以充分进步求解精度。为了分析不同计算模型对求解精度的影响，并对小湾拱坝和二滩拱坝进展分析比拟，共建立了四套有限元模型，其中小湾拱坝粗网格模型和精细网格模型各一套，二滩拱坝粗网格模型和精细网格模型各一套。图2小湾拱坝精细网格(32024单元)图3二滩拱坝精细网格(20876单元)2.2.2计算工况分别对考虑水力劈裂作用和不考虑水力劈裂作用两种情况进展了计算。考虑水力劈裂作用时，在已开裂缝面上施加上游全水头。此外，对坝体混凝土的单轴抗拉强度分别取0Pa、1Pa和2Pa进展敏感性

12、分析。以上情况共组合成24种工况。2.3计算结论图4和图5分别为小湾拱坝和二滩拱坝建基面裂缝扩展范围示意图(混凝土抗拉强度1Pa，考虑水力劈裂)，图6和图7分别为在一样工况下小湾拱坝和二滩拱坝拱冠梁断面开裂后应力矢量图。图4小湾拱坝建基面裂缝扩展示意(混凝土抗拉强度1Pa，考虑水力劈裂)图5二滩拱坝建基面裂缝扩展示意(混凝土抗拉强度1Pa，考虑水力劈裂)(1)计算裂缝的扩展范围大小与网格形式有较大的关系。网格加密后，开裂范围普遍有所减小，小湾模型的开裂深度占底厚的百分比平均减小约8.68个百分点，二滩模型平均减小约7.14个百分点。(2)随着混凝土抗拉强度的增加，裂缝的扩展范围逐渐缩校混凝土抗

13、拉强度由0Pa进步到1Pa时，小湾模型的开裂深度占底厚的百分比平均减小约10.42个百分点，二滩模型平均减小约9.52个百分点。混凝土抗拉强度由1Pa进步到2Pa时，小湾模型的开裂深度占底厚的百分比平均减小约13.02个百分点，二滩模型平均减小约9.52个百分点。图6小湾拱坝拱冠梁断面应力矢量图7二滩拱坝拱冠梁断面应力矢量(3)在其它条件一样情况下，考虑水力劈裂后的裂缝扩展范围比不考虑水力劈裂时大，小湾模型的开裂深度占底厚的百分比平均增加约7.99个百分点，二滩模型平均增加约10.32个百分点。(4)混凝土抗拉强度取1Pa时，小湾拱坝的最大绝对开裂深度为31.90，最大相对开裂深度(开裂深度:

14、坝基宽度)为0.4375；二滩拱坝的最大绝对开裂深度为25.73，最大相对开裂深度(开裂深度:坝基宽度)为0.4615。(5)在目前所考虑的影响因素下，小湾拱坝和二滩拱坝的开裂范围、开裂程度及对网格疏密、混凝土抗拉强度和水力劈裂等影响因素的敏感程度都非常相似。考虑到二滩工程已建成，且运行正常，可以初步断定小湾工程的拱坝坝踵开裂危险性不大。3结语本文介绍了三维裂缝扩展的不变网格有限元分析方法。拱坝开裂分析结果的影响因素很多，其中荷载、边界条件、开裂力学模型与参数、计算方法等都很重要却又很难准确把握。而且，作为一种新方法，配套的平安准那么取值也需进一步研究。基于拱坝开裂问题的复杂性，在对小湾工程拱

15、坝进展应用研究时，取已建的二滩工程为校准工程，在一样的工况下进展比照计算，比拟其裂缝扩展范围与坝基宽度的相比照值，从而可对小湾高拱坝的开裂危险性和平安度进展评价。能对拱坝裂缝扩展进展三维分析的方法很多。笔者体会，包括本文介绍的方法在内，各种模型和方法都尚不够成熟以可靠地模拟裂缝发生位置、开展过程和构造的真实平安度。以下几个方向的研究可能会有助于问题的解决：(1)考虑真实荷载与边界条件的施工和运行过程仿真模拟；(2)在拉剪和压剪状态下混凝土及岩体的开裂机理和准那么(包括起裂条件、开裂长度与方向)以及对应的参数；(3)任意形状与分布的三维裂缝的定位与追踪技术；(4)开裂深度平安准那么。参考文献：1

16、朱伯芳，栾丰.拱与梁产生裂缝后的失效角J.水力发电学报，1997，(3)：55-60.2杨海霞，杜成斌，王德信.拱坝非线性开裂分析的分载位移法J.水力发电，1997，(7)：19-22.3LinsbaurerHN,IngraffeaAR.Siulatinfarakinginlargearhda.Part1J.J.fStrut.Eng.,1989,115(7):1599-1615.4LinsbaurerHN,IngraffeaAR.Siulatinfarakinginlargearhda.Part2.J.J.fStrut.Eng.,1989,115(7):1616-1630.5石根华.数值流形方法与非连续变形分析.北京：清华大学出版社，1997.6BelytshT,LuYY,GuL.Eleent-freeGalerkinethdsJ.Int.J.frNu.eth.inEng.,1994,37:229-256.7冠晓东，周维垣.应用无单元法近似计算拱坝开裂J.水利学报，2000，(10)：28-35.8DhndtG.Autati3DdeIrakprpagatinalulatinsithfiniteeleentsJ.Int.J.frNu.eth.inEing.,1998,41:739-757.9王光纶，张楚汉，王少敏.混凝土重力坝的非线性断裂分析J.水利学报，1997，(增