拱坝-地基破坏的数值模型与溃坝仿真.doc

拱坝-地基破坏的数值模型与溃坝仿真 摘要：本文使用刚体弹簧元模型来模拟拱坝-地基这类非连续介质和连续介质组成的系统，对一个悬臂梁的位移场和应力场的计算并与有限元方法计算结果的对比表明，这种方法具有良好的精度.应用这一模型对玛尔帕塞拱坝的溃坏进行了仿真，重现了拱坝的破坏过程，探讨了拱坝坝肩滑动失稳的原因，计算结果表明坝肩岩体裂隙和断层的强度不足是导致整个拱坝溃毁的关键因素. 关键词：三维刚体弹簧元 拱坝地基破坏 拱坝 仿真 　　近年来，用非连续介质力学模型研究高坝地基、高边坡与地下工程的变形稳定有长足进展，离散元[1～4]、dda[5]、刚体弹簧元[6]等的发展提供了求解非连续介质力学问题的有效工具.关于高坝-地基系统的破坏机理与仿真分析，由于问题的极端复杂性，迄今仍未见诸报导，本文利用改进的刚体弹簧元研究玛尔帕塞拱坝溃坏的机理与过程，是试图用连续-非连续介质统一的数值模型研究高坝-地基系统破坏过程的尝试. 法国玛尔帕塞拱坝1959年12月2日的溃坝事故是坝工史上的重大事件[7,8]，引起了坝工界的极大重视.世界各国专家和学者从未停止过对玛尔帕塞拱坝失事原因的探究，并提出了各自见解[9]，归纳起来有以下3种看法：(1)坝基变形，认为左岸坝基修建在弱风化顶层上，岩体弹性模量低，仅1000mpa左右，地基的过大变形使拱的推力向上向下转移，其中向上转移的推力使左岸重力墩超载而产生位移，导致拱坝失去支撑而破坏.(2)坝体上滑失稳，认为岩体浅层的裂隙在不利的条件下不能提供足够的抗滑力，致使坝体沿建基面向上和向下游滑动，拱圈拱弦被拉长，使坝沿一岸或两岸拱端转动，导致拱冠断裂，坝体崩毁.(3)坝肩岩体滑动，认为左岸坝基存在一个潜在的由下游断层、发育的节理面，以及下游临空面组成的滑动楔形体.由于地质条件在楔形岩体上游侧形成很高的静水扬压力使岩体滑移失稳.本文利用改进的三维刚体弹簧元对玛尔帕塞拱坝的失稳破坏过程进行仿真，对该坝的失稳机理作进一步的探讨. 1 拱坝破坏仿真模型 1.1 三维刚体弹簧元模型 假定三维刚体弹簧元的单元块体为不可变形的刚性体，在引入了刚度、强度等效准则后，将块体变形与强度指标统一集中于接触弹簧上.块体的运动符合刚体的平动和转动定律，块体与块体依靠相互之间连接的切向和法向弹簧和阻尼器来传递相互之间的作用，如图1.弹簧的作用力由块体间的相对位移确定.如图2，块体p对块体q的相对位移为δu，则他们之间的弹簧力变化量为 ? 式中：ks、kn分别为切向和法向弹簧的刚度；δs、δn分别为相对位移增量的切向和法向分量；δs、δn分别为切向和法向弹簧力的变化量. ? 图1 接触模型 图2 相对位移增量示意 1.2 使用三维刚体弹簧元模拟非连续介质 三维刚体弹簧元模型基于离散元方法提出，是非连续介质模型的1种，由于块体之间可以错动、分离，因此可以方便地模拟非连续体如裂隙发育岩体的大变形、大位移问题.离散元方法中接触关系一般可分为角面、角边、角角和边边四类，关于各类接触关系的检索和其它离散元的详细内容已在cundall和鲁军等人的文章中有详细叙述，参见文献[1～4]. 1.3 使用三维刚体弹簧元模拟连续介质 ? 图3 刚体弹簧元模拟连续介质原理示意 ? 图4 中心受拉杆件 1.3.1基本原理 如图3，通过将连续介质的变形和应力凝聚到单元之间连接的法向和切向弹簧上，可以用刚体弹簧元模型来模拟连续介质.即由弹簧的变形刚度(kn、ks)来等效连续介质的变形模量(e、g)，由弹簧力系统(fn、fs)来推求与之静力等效的连续介质的应力(σ、τ)分布，从而将刚体弹簧元系统和连续介质联系起来.(1)刚度等效.如图4所示，一连续的l×l1×l2方杆，弹性模量为e，泊松比为υ.将杆沿长度方向等分成n个l3×l1×l2的块体，其中l3=l／n，每个块体间的面面接触被视为四个角点处的点面接触，当n足够大，取弹簧的法向刚度kn=l1×l2×e/4l3时，使用三维刚体弹簧元计算它在中心作用力p下的变形与使用材料力学方法的结果是一致的.弹簧的切向刚度借用弹性模量和剪切模量的关系可取ks=kn/[2(1+υ)].单元应尽可能采用长方体，对于不规则的块体可以将其均化成长方体使用上述估算式，或另外推求相应的计算式. (2)强度等效.如果刚体弹 簧元网格划分得比较密，则可以取单元面上每个弹簧力合力除以该面的面积求得其平均应力作为该面上的应力值.因为单元网格较密，可以认为当面上的平均应力达到材料的强度时该面上的弹簧全部断开，则单元在该面发生破坏.由于刚体弹簧元的计算工作量比较大，当为了节省计算工作量，网格划分得较稀疏时，可以按照材料力学方法假定单元面上的应力分布方式，由弹簧力根据静力等效原则求出该面上的应力分布. 1 / 1 1 /