拱坝1Word文件下载.docx

1、关键词三维刚体弹簧元 拱坝地基破坏 拱坝 仿真近年来，用非连续介质力学模型研究高坝地基、高边坡与地下工程的变形稳定有长足进展，离散元14、DDA、刚体弹簧元等的发展提供了求解非连续介质力学问题的有效工具.关于高坝-地基系统的破坏机理与仿真分析，由于问题的极端复杂性，迄今仍未见诸报导，本文利用改进的刚体弹簧元研究玛尔帕塞拱坝溃坏的机理与过程，是试图用连续-非连续介质统一的数值模型研究高坝-地基系统破坏过程的尝试. 法国玛尔帕塞拱坝1959年12月2日的溃坝事故是坝工史上的重大事件7,8，引起了坝工界的极大重视.世界各国专家和学者从未停止过对玛尔帕塞拱坝失事原因的探究，并提出了各自见解，归纳起来有

2、以下3种看法：（1）坝基变形，认为左岸坝基修建在弱风化顶层上，岩体弹性模量低，仅1000MPa左右，地基的过大变形使拱的推力向上向下转移，其中向上转移的推力使左岸重力墩超载而产生位移，导致拱坝失去支撑而破坏.（2）坝体上滑失稳，认为岩体浅层的裂隙在不利的条件下不能提供足够的抗滑力，致使坝体沿建基面向上和向下游滑动，拱圈拱弦被拉长，使坝沿一岸或两岸拱端转动，导致拱冠断裂，坝体崩毁.（3）坝肩岩体滑动，认为左岸坝基存在一个潜在的由下游断层、发育的节理面，以及下游临空面组成的滑动楔形体.由于地质条件在楔形岩体上游侧形成很高的静水扬压力使岩体滑移失稳.本文利用改进的三维刚体弹簧元对玛尔帕塞拱坝的失稳破

3、坏过程进行仿真，对该坝的失稳机理作进一步的探讨.1 拱坝破坏仿真模型 三维刚体弹簧元模型 假定三维刚体弹簧元的单元块体为不可变形的刚性体，在引入了刚度、强度等效准则后，将块体变形与强度指标统一集中于接触弹簧上.块体的运动符合刚体的平动和转动定律，块体与块体依靠相互之间连接的切向和法向弹簧和阻尼器来传递相互之间的作用，如图1.弹簧的作用力由块体间的相对位移确定.如图2，块体P对块体Q的相对位移为u，则他们之间的弹簧力变化量为式中：Ks、Kn分别为切向和法向弹簧的刚度；s、n分别为相对位移增量的切向和法向分量；s、n分别为切向和法向弹簧力的变化量.图1 接触模型图2 相对位移增量示意使用三维刚体弹

4、簧元模拟非连续介质 三维刚体弹簧元模型基于离散元方法提出，是非连续介质模型的1种，由于块体之间可以错动、分离，因此可以方便地模拟非连续体如裂隙发育岩体的大变形、大位移问题.离散元方法中接触关系一般可分为角面、角边、角角和边边四类，关于各类接触关系的检索和其它离散元的详细内容已在Cundall和鲁军等人的文章中有详细叙述，参见文献14. 使用三维刚体弹簧元模拟连续介质图3 刚体弹簧元模拟连续介质原理示意图4 中心受拉杆件基本原理 如图3，通过将连续介质的变形和应力凝聚到单元之间连接的法向和切向弹簧上，可以用刚体弹簧元模型来模拟连续介质.即由弹簧的变形刚度（Kn、Ks）来等效连续介质的变形模量（E

5、、G），由弹簧力系统（Fn、Fs）来推求与之静力等效的连续介质的应力（、）分布，从而将刚体弹簧元系统和连续介质联系起来.（1）刚度等效.如图4所示，一连续的ll1l2方杆，弹性模量为E，泊松比为.将杆沿长度方向等分成n个l3l2的块体，其中l3=ln，每个块体间的面面接触被视为四个角点处的点面接触，当n足够大，取弹簧的法向刚度Kn=l1l2E/4l3时，使用三维刚体弹簧元计算它在中心作用力P下的变形与使用材料力学方法的结果是一致的.弹簧的切向刚度借用弹性模量和剪切模量的关系可取Ks=Kn/2（1+）.单元应尽可能采用长方体，对于不规则的块体可以将其均化成长方体使用上述估算式，或另外推求相应的计

6、算式.（2）强度等效.如果刚体弹簧元网格划分得比较密，则可以取单元面上每个弹簧力合力除以该面的面积求得其平均应力作为该面上的应力值.因为单元网格较密，可以认为当面上的平均应力达到材料的强度时该面上的弹簧全部断开，则单元在该面发生破坏.由于刚体弹簧元的计算工作量比较大，当为了节省计算工作量，网格划分得较稀疏时，可以按照材料力学方法假定单元面上的应力分布方式，由弹簧力根据静力等效原则求出该面上的应力分布.（3）荷载等效.对于集中荷载，将荷载作用在相应的块体形心上，并将荷载对中心的力矩加入到该块体的合力矩中即与原荷载的作用等效.对于分布荷载也要等效成集中荷载加入到相应的块体中.图5简单算例简单算例

7、如图5，一长40m，高4m，宽1m的悬臂梁受均匀荷载P=105N/m2作用，支座A固定，不考虑自重，坐标见图中所示.悬臂梁弹性模量E=1010N/m2，泊松比=03，密度=103kg/m3，离散网格如图5（在厚度方向为一层，厚1m）.使用上述的刚体弹簧元方法计算悬臂梁的变形和应力并与有限元的结果进行对比，图6为悬臂梁沿y轴的位移对比结果，图7和图8分别为x方向的正应力x和剪应力xy的对比结果.从图中可以看出两种方法计算的位移和应力结果除有限元法在固定端能反应局部应力集中外，整体应力具有很好的近似性.图6 y方向位移等值线（单位：m） 图7 x方向正应力x等值线（单位：MPa） 图8 剪应力xy

8、等值线（单位：MPa）刚体弹簧元通过引入单元表面应力分布假设和变形近似等效来模拟连续介质，利用这些假设，刚体弹簧元通过上述的刚度等效、强度等效和荷载等效，能够把连续介质和非连续介质统一在一个模型之中，为研究坝肩失稳导致拱坝断裂破坏需要仿真连续与非连续耦合介质的大变形与破坏断裂的复杂过程提供了有力的工具.2 玛尔帕塞拱坝失稳机理的研究 模型的建立 坝肩和坝体的离散和计算条件的设定 法国玛尔帕塞拱坝坝高66m，坝顶高程，坝顶长，长高比，坝顶厚，底厚，厚高比为坝址岩体由带状片麻岩组成，含千节理倾向下游右岸，倾角3050.该坝于1954年建成，蓄水历时5年至1959年12月蓄至最高水位约100m时，坝

9、在短瞬间溃决，左岸与中部坝体全部冲走，右岸底部坝体残存.根据玛尔帕塞拱坝的地形和地质资料以及拱坝的体形和横缝设置将坝址处的岸坡简化成如图9所示，坐标系设置和拱坝离散网格如图10.在分析中边界条件设定为除“二面沟”内的岩体外，其他岩体均被固定.坝体、右岸坝肩和中间岩体依照连续介质进行模拟，坝体弹性模量取E=1010Pa，坝体材料的抗拉强度取，抗压强度为20MPa，摩擦系数f=；考虑到右岸坝基事后调查仍有坝体残存，表明该区岩体强度较高，故将坝体右岸及中间部分与基础的交界面强度较坝体提高1倍.左岸坝肩依照非连续体计算，重力墩和地基之间的摩擦系数设定为，坝体和左岸坝肩之间以及“二面沟”内的岩体两组节理

10、面与水平层面取相同的摩擦系数f1，根据不同的工况选取，上、下游断层摩擦系数为f2，亦根据不同的工况进行选择.计算中未考虑粘结强度c的作用.主要荷载为库水作用，稳定分析中库水水位均取100m，考虑左坝肩扬压力时，按静水压力作用施加在沿上游面坝基线的裂隙处，作用面见图9中虚线段所示.为了简化计算，所有荷载均一次施加.计算中阻尼参照Cundall提供的阻尼模式选取刚度与质量成比例的Rayleigh阻尼型式，阻尼系数、由试算确定，试算表明：如不出现过阻尼或阻尼过小情况计算结果均收敛于一稳态解的近似值，详见参考文献.位移分析 计算拱坝位移时左坝肩暂按连续介质考虑，并且不考虑扬压力作用，得到拱坝在1955

11、年9月（水位）蓄水至1959年7月（水位）拱坝产生的位移，并与相应实测位移值进行比较，图11列出了顶拱和拱冠梁处的比较结果.可以看出二者最大值基本接近，但在拱冠梁坝基附近，实测的位移值较计算大1cm左右，顶拱计算结果较实测结果略偏大.这一比较结果表明选择的弹簧刚度参数Kn、Ks以及单元等效模型基本符合实际. 玛尔帕塞拱坝稳定性分析及失稳过程的仿真 通过改变坝体和左岸坝肩之间以及左岸坝肩“二面沟”内岩体的摩擦系数f1和上下游断层面的摩擦系数f2可以模拟和分析拱坝在各类参数条件下失稳破坏的情形.图9岸坡离散图10刚体弹簧元离散网格图11拱坝位移计算结果与实测结果比较工况1不考虑扬压力的浅层滑动 不

12、考虑左岸坝肩处的扬压力，将坝体和左岸坝肩之间以及左岸坝肩“二面沟”岩体的摩擦系数f1和上、下游断层面的摩擦系数f2设为相同，从开始折减，逐步计算，直到拱坝发生破坏.计算结果表明当摩擦系数f1和f2折减到时坝肩失稳，拱坝破坏.图12列出了左坝肩附近的块体1（见图10）的位移过程.f=时，块体位移量陡增，接近30cm，此时拱坝处于失稳的临界状态，f=时，块体位移不再收敛，开始沿着坝肩向下游、向上滑动，整个拱坝随之破坏.图12 块体1的位移过程（工况1） 图13 f=时左坝肩重力墩位移过程（工况1）分析图13图15，拱坝的整个破坏过程可以分为两个阶段，第一阶段是大约3s以前拱坝整体位移，第二阶段是3

13、s以后由于过大的变位使拱坝的工作方式发生改变，应力恶化，坝体溃毁.描述在第一阶段，拱坝利用拱的作用将水压力转化为水平作用在坝肩上的推力，沿坝肩向上的分量使拱圈产生向上滑动的趋势，由于左坝肩岩体的摩擦系数较小，不能提供足够的抗力，巨大的拱推力转移到了左岸的重力墩上，使重力墩沿着墩轴线偏向下游移动了近2m（见图13，3s前重力墩主要以平动为主），同时整个坝体沿建基面向下游、向上滑动，带动坝体以右岸为轴向下游转动（由图14中y方向的位移分布，坝体向下游位移由右到左依次增大，显示以右坝肩为轴向下游转动）和向上旋转（见图14中z方向的位移分布，坝体上抬量由右到左依次增大，显示以右坝肩为轴向上转动）.此时

14、，整个拱坝基本上保持整体发生变位，同时有沿右岸转动趋势.（a）y方向位移等值线（b）z方向位移等值线 图14 f=，3s时拱坝坝体位移分布（工况1，单位：m）拱坝在经过了第一阶段后，拱圈拱弦拉长，逐渐被压扁，承载能力极度降低，应力发生恶化，坝体破裂，拱坝的破坏进入到第二阶段破坏阶段.由于中、下部块体受水压力较大，拱坝坝体从该部位偏左岸发生断裂，向下游弯折，首先被冲出，随后其周围坝块跟着折向下游，上部坝块受水压力较小，发生翻转后随在中下部坝块之后冲出，拱坝彻底破坏.与重力墩连接的坝段在被向下游冲出的过程中又带动重力墩发生了向下游的位移和旋转（见图13中3s以后的位移，向下游移动1m左右，在平面上

15、旋转约3）.图15给出了上述过程一些典型时刻拱坝的状态，拱坝破坏后最后残留状态见最后一幅.图15 f=时拱坝破坏过程（工况1）上述拱坝破坏过程的计算结果和拱坝失事后的调查结果基本相符：（1）f值的大小，“实有摩擦角为3035”，本文计算表明f值取时拱坝处于临界状态，对应的摩擦角为35，当拱坝左坝肩的实有平均摩擦角小于该值时，拱坝会发生上滑失稳，与拱坝破坏事实相符；（2）坝体的残留部分，拱坝失事后“右岸部分的坝体和中间部分的坝体底部基础留存了下来”，坝基由于嵌固作用，一般不易破坏，当拱坝上滑失稳时，应力发生恶化，拱坝沿横缝和工作缝等薄弱处断裂，同时底部的坝体得到卸载，从而保留了下来，计算中将坝体右岸及中间部分与基础的交界面强度较坝体提高了1倍，坝体的残留部分如图15，和拱坝破坏后的结果相符；（3）重力墩的位移，“左按重力墩在沿墩轴线方向上已移动约2m，沿垂直于墩轴线方向，向下游移动了约80cm，并略有倾侧”，本文计算的过程解释了上述重力墩先沿轴线方向移动（拱坝整体滑动时的推力所致），又向下游移动并转动（拱坝破坏后，重力墩附近的坝段向下游破坏时的带动）的现象，位移量也基本接近； （4）残留部分的位移规律，“AB坝段的位移很小，各坝段的位移值自右岸向左岸逐渐增大，至JK坝段处增达80cm；位移