五强溪船闸裂缝的流形元模拟.pdfVIP

水 利 学 报 2003 年11 月 SHUILI XUEBAO 第11 期 文章编号：0559-9350 (2003) 11-0037-06 五强溪船闸裂缝的流形元模拟 1 2 1 3 张国新 ，周立本 ，朱伯芳 ，徐润明 (1.中国水利水电科学研究院 结构材料所，北京 100038，2.国家电力公司 中南勘测设计研究院，湖南长沙 410014，3.湖 南省五强溪水电厂，湖南常德 415000) 摘要：对流形元法进行了扩展，开发了考虑裂隙渗流和变形耦合作用、用断裂力学法模拟离散裂缝的扩展及用杆单 元模拟钢筋作用的流形元法，用以模拟钢筋混凝土结构在水压荷载作用下的破坏过程。用扩展后的软件模拟分析了 五强溪船闸闸墙的大变形问题。模拟结果表明，五强溪船闸闸墙的超标变形是由于闸室底部大量纵深裂缝所致，缝 内渗透压力在裂缝扩展中起了重要作用。建议船闸提高水位运行前应做止水或加固处理。 关键词：流形元；裂缝；船闸；渗流；大变形 中图分类号：TV222.2 文献标识码：A 传统的基于连续介质力学的有限元法、边界元法等可以很好地计算连续介质内的应力与位移，但对于 裂缝扩展的模拟却遇到很大困难，更难以模拟块体运动等大变形问题。建立在非连续力学基础之上的 DEM 法和DDA 法可以模拟块体运动、接触等不连续问题，但却难以给出块体内部尤其是连续介质内部的合理应 力分布，更难以模拟连续体内裂缝的发生和扩展。石根华提出的流形元法［1］克服了以往方法的不足，既吸 收了有限元等连续力学方法的优点，又吸收了DDA 等不连续方法的长处。其像有限元法那样在连续介质中 划分单元，使该方法与有限元法同等的精度计算连续域内的应力分布，物理网格和数学网格的分离及一套 完善的接触理论的引入，使该法能很好地模拟接触、块体运动等不连续力学问题，Cover 的使用又使该方 法较容易和客观地模拟裂缝的产生和扩展。本文在文献［2］的基础上进一步完善了流形元法模拟裂缝产 生与扩展的计算方法和程序，同时开发了渗流与变形耦合的流形元算法，并由此模拟水力劈裂。用开发的 程序模拟了五强溪船闸在使用过程中随水位的不断提高裂缝的产生与扩展过程，并将船闸变形结果与实测 结果进行了比较。 1 模拟方法 1.1 二阶流行法 模拟裂缝的扩展需要应力场、位移场结果具有较高的精度，一阶流行法是难以胜任的。 适用于流形元法的高阶位移函数的构造有两种方法：第一种是利用线性权函数与高阶Cover 函数；另一种 是用高阶权函数和常数Cover 函数。由于裂缝向任意可能方向扩展后，可能产生很小的单元，第一种方法 在遇到很小单元的时候容易使计算发散，且当变形较大时易使物理单元的积分误差增大，因此本文借鉴有 限元的思路，用6 节点三角形单元构造二阶插值函数，即用二阶权函数和数学网格6 个顶点上的常量位移 构成二阶位移函数，即 6 u x y u  ( , )  i    ∑w (x , y )  (1) i v(x ,y )  i 1 v  i  收稿日期：2003-06-26 作者简介：张国新(1960-)，男，山东人，教授级高级工程师，博士，主要研究方向为水工结构数值计算。 37