拱坝地基破坏的数值模型与溃坝仿真2行政.docVIP

拱坝地基破坏的数值模型与溃坝仿真2行政论文范文大全 拱坝-地基破坏的数值模型与溃坝仿真 2 玛尔帕塞拱坝失稳机理的研究 2.1 模型的建立 2.1.1 坝肩和坝体的离散和计算条件的设定 法国玛尔帕塞拱坝坝高66m，坝顶高程102.55m，坝顶长222.7m，长高比3.3，坝顶厚1.5m，底厚6.78m，厚高比为0.1.坝址岩体由带状片麻岩组成，含千节理倾向下游右岸，倾角30°～50°.该坝于1954年建成，蓄水历时5年至1959年12月蓄至最高水位约100m时，坝在短瞬间溃决，左岸与中部坝体全部冲走，右岸底部坝体残存.根据玛尔帕塞拱坝的地形和地质资料以及拱坝的体形和横缝设置将坝址处的岸坡简化成如图9所示，坐标系设置和拱坝离散网格如图10.在分析中边界条件设定为除“二面沟”内的岩体外，其他岩体均被固定.坝体、右岸坝肩和中间岩体依照连续介质进行模拟，坝体弹性模量取e=1010pa，坝体材料的抗拉强度取2.0mpa，抗压强度为20mpa，摩擦系数f=1.0；考虑到右岸坝基事后调查仍有坝体残存，表明该区岩体强度较高，故将坝体右岸及中间部分与基础的交界面强度较坝体提高1倍.左岸坝肩依照非连续体计算，重力墩和地基之间的摩擦系数设定为1.0，坝体和左岸坝肩之间以及“二面沟”内的岩体两组节理面与水平层面取相同的摩擦系数f1，根据不同的工况选取，上、下游断层摩擦系数为f2，亦根据不同的工况进行选择.计算中未考虑粘结强度c的作用.主要荷载为库水作用，稳定分析中库水水位均取100m，考虑左坝肩扬压力时，按静水压力作用施加在沿上游面坝基线的裂隙处，作用面见图9中虚线段所示.为了简化计算，所有荷载均一次施加.计算中阻尼参照cundall提供的阻尼模式选取刚度与质量成比例的rayleigh阻尼型式，阻尼系数α、β由试算确定，试算表明：如不出现过阻尼或阻尼过小情况计算结果均收敛于一稳态解的近似值，详见参考文献 . 2.1.2位移分析 计算拱坝位移时左坝肩暂按连续介质考虑，并且不考虑扬压力作用，得到拱坝在1955年9月(水位79.75m)蓄水至1959年7月(水位94.19m)拱坝产生的位移，并与相应实测位移值进行比较，图11列出了顶拱和拱冠梁处的比较结果.可以看出二者最大值基本接近，但在拱冠梁坝基附近，实测的位移值较计算大1cm左右，顶拱计算结果较实测结果略偏大.这一比较结果表明选择的弹簧刚度参数kn、ks以及单元等效模型基本符合实际. 2.2 玛尔帕塞拱坝稳定性分析及失稳过程的仿真 通过改变坝体和左岸坝肩之间以及左岸坝肩“二面沟”内岩体的摩擦系数f1和上下游断层面的摩擦系数f2可以模拟和分析拱坝在各类参数条件下失稳破坏的情形. 图9岸坡离散 图10刚体弹簧元离散网格 图11拱坝位移计算结果与实测结果比较 2.2.1工况1——不考虑扬压力的浅层滑动 不考虑左岸坝肩处的扬压力，将坝体和左岸坝肩之间以及左岸坝肩“二面沟”岩体的摩擦系数f1和上、下游断层面的摩擦系数f2设为相同，从1.0开始折减，逐步计算，直到拱坝发生破坏.计算结果表明当摩擦系数f1和f2折减到0.6时坝肩失稳，拱坝破坏.图12列出了左坝肩附近的块体1(见图10)的位移过程.f=0.7时，块体位移量陡增，接近30cm，此时拱坝处于失稳的临界状态，f=0.6时，块体位移不再收敛，开始沿着坝肩向下游、向上滑动，整个拱坝随之破坏. 图12 块体1的位移过程(工况1) 图13 f=0.6时左坝肩重力墩位移过程(工况1) 分析图13～图15，拱坝的整个破坏过程可以分为两个阶段，第一阶段是大约3s以前拱坝整体位移，第二阶段是3s以后由于过大的变位使拱坝的工作方式发生改变，应力恶化，坝体溃毁.描述如下：在第一阶段，拱坝利用拱的作用将水压力转化为水平作用在坝肩上的推力，沿坝肩向上的分量使拱圈产生向上滑动的趋势，由于左坝肩岩体的摩擦系数较小，不能提供足够的抗力，巨大的拱推力转移到了左岸的重力墩上，使重力墩沿着墩轴线偏向下游移动了近2m(见图13，3s前重力墩主要以平动为主)，同时整个坝体沿建基面向下游、向上滑动，带动坝体以右岸为轴向下游转动(由图14中y方向的位移分布，坝体向下游位移由右到左依次增大，显示以右坝肩为轴向下游转动)和向上旋转(见图14中z方向的位移分布，坝体上抬量由右到左依次增大，显示以右坝肩为轴向上转动).此时，整个拱坝基本上保持整体发生变位，同时有沿右岸转动趋势. (a)y方向位移等值线 (b)z方向位移等值线 图14 f=0.6，3s时拱坝坝体位移分布(工况1，单位：m) 拱坝在经过了第一阶段后，拱圈拱弦拉长，逐渐被压扁，承载能力极度降低，应力发生恶化，坝体破裂，拱坝的破坏进入到第二阶段——破坏阶段.由于中、下部块体受水压力较大，拱坝坝体从该部位偏左岸发生断裂，向下游弯折