欢迎来到人人文库网! | 帮助中心 人人文库renrendoc.com美如初恋!
人人文库网
首页 人人文库网 > 资源分类 > DOC文档下载

水利工程论文-拱坝-地基破坏的数值模型与溃坝仿真.doc

  • 资源大小:44.98KB        全文页数:11页
  • 资源格式: DOC        下载权限:游客/注册会员/VIP会员    下载费用:2
游客快捷下载 游客一键下载
会员登录下载
下载资源需要2

邮箱/手机号:
您支付成功后,系统会自动为您创建此邮箱/手机号的账号,密码跟您输入的邮箱/手机号一致,以方便您下次登录下载和查看订单。注:支付完成后需要自己下载文件,并不会自动发送文件哦!

支付方式: 微信支付    支付宝   
验证码:   换一换

友情提示
2、本站资源不支持迅雷下载,请使用浏览器直接下载(不支持QQ浏览器)
3、本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

水利工程论文-拱坝-地基破坏的数值模型与溃坝仿真.doc

水利工程论文拱坝地基破坏的数值模型与溃坝仿真摘要本文使用刚体弹簧元模型来模拟拱坝地基这类非连续介质和连续介质组成的系统,对一个悬臂梁的位移场和应力场的计算并与有限元方法计算结果的对比表明,这种方法具有良好的精度应用这一模型对玛尔帕塞拱坝的溃坏进行了仿真,重现了拱坝的破坏过程,探讨了拱坝坝肩滑动失稳的原因,计算结果表明坝肩岩体裂隙和断层的强度不足是导致整个拱坝溃毁的关键因素关键词三维刚体弹簧元拱坝地基破坏拱坝仿真近年来,用非连续介质力学模型研究高坝地基、高边坡与地下工程的变形稳定有长足进展,离散元1~4、DDA5、刚体弹簧元6等的发展提供了求解非连续介质力学问题的有效工具关于高坝地基系统的破坏机理与仿真分析,由于问题的极端复杂性,迄今仍未见诸报导,本文利用改进的刚体弹簧元研究玛尔帕塞拱坝溃坏的机理与过程,是试图用连续非连续介质统一的数值模型研究高坝地基系统破坏过程的尝试法国玛尔帕塞拱坝1959年12月2日的溃坝事故是坝工史上的重大事件7,8,引起了坝工界的极大重视世界各国专家和学者从未停止过对玛尔帕塞拱坝失事原因的探究,并提出了各自见解9,归纳起来有以下3种看法1坝基变形,认为左岸坝基修建在弱风化顶层上,岩体弹性模量低,仅1000MPA左右,地基的过大变形使拱的推力向上向下转移,其中向上转移的推力使左岸重力墩超载而产生位移,导致拱坝失去支撑而破坏2坝体上滑失稳,认为岩体浅层的裂隙在不利的条件下不能提供足够的抗滑力,致使坝体沿建基面向上和向下游滑动,拱圈拱弦被拉长,使坝沿一岸或两岸拱端转动,导致拱冠断裂,坝体崩毁3坝肩岩体滑动,认为左岸坝基存在一个潜在的由下游断层、发育的节理面,以及下游临空面组成的滑动楔形体由于地质条件在楔形岩体上游侧形成很高的静水扬压力使岩体滑移失稳本文利用改进的三维刚体弹簧元对玛尔帕塞拱坝的失稳破坏过程进行仿真,对该坝的失稳机理作进一步的探讨1拱坝破坏仿真模型11三维刚体弹簧元模型假定三维刚体弹簧元的单元块体为不可变形的刚性体,在引入了刚度、强度等效准则后,将块体变形与强度指标统一集中于接触弹簧上块体的运动符合刚体的平动和转动定律,块体与块体依靠相互之间连接的切向和法向弹簧和阻尼器来传递相互之间的作用,如图1弹簧的作用力由块体间的相对位移确定如图2,块体P对块体Q的相对位移为ΔU,则他们之间的弹簧力变化量为式中KS、KN分别为切向和法向弹簧的刚度;ΔS、ΔN分别为相对位移增量的切向和法向分量;ΔS、ΔN分别为切向和法向弹簧力的变化量图1接触模型图2相对位移增量示意12使用三维刚体弹簧元模拟非连续介质三维刚体弹簧元模型基于离散元方法提出,是非连续介质模型的1种,由于块体之间可以错动、分离,因此可以方便地模拟非连续体如裂隙发育岩体的大变形、大位移问题离散元方法中接触关系一般可分为角面、角边、角角和边边四类,关于各类接触关系的检索和其它离散元的详细内容已在CUNDALL和鲁军等人的文章中有详细叙述,参见文献1~413使用三维刚体弹簧元模拟连续介质图3刚体弹簧元模拟连续介质原理示意图4中心受拉杆件131基本原理如图3,通过将连续介质的变形和应力凝聚到单元之间连接的法向和切向弹簧上,可以用刚体弹簧元模型来模拟连续介质即由弹簧的变形刚度KN、KS来等效连续介质的变形模量E、G,由弹簧力系统FN、FS来推求与之静力等效的连续介质的应力Σ、Τ分布,从而将刚体弹簧元系统和连续介质联系起来1刚度等效如图4所示,一连续的LL1L2方杆,弹性模量为E,泊松比为Υ将杆沿长度方向等分成N个L3L1L2的块体,其中L3L/N,每个块体间的面面接触被视为四个角点处的点面接触,当N足够大,取弹簧的法向刚度KNL1L2E/4L3时,使用三维刚体弹簧元计算它在中心作用力P下的变形与使用材料力学方法的结果是一致的弹簧的切向刚度借用弹性模量和剪切模量的关系可取KSKN/21Υ单元应尽可能采用长方体,对于不规则的块体可以将其均化成长方体使用上述估算式,或另外推求相应的计算式2强度等效如果刚体弹簧元网格划分得比较密,则可以取单元面上每个弹簧力合力除以该面的面积求得其平均应力作为该面上的应力值因为单元网格较密,可以认为当面上的平均应力达到材料的强度时该面上的弹簧全部断开,则单元在该面发生破坏由于刚体弹簧元的计算工作量比较大,当为了节省计算工作量,网格划分得较稀疏时,可以按照材料力学方法假定单元面上的应力分布方式,由弹簧力根据静力等效原则求出该面上的应力分布3荷载等效对于集中荷载,将荷载作用在相应的块体形心上,并将荷载对中心的力矩加入到该块体的合力矩中即与原荷载的作用等效对于分布荷载也要等效成集中荷载加入到相应的块体中图5简单算例132简单算例如图5,一长40M,高4M,宽1M的悬臂梁受均匀荷载P105N/M2作用,支座A固定,不考虑自重,坐标见图中所示悬臂梁弹性模量E1010N/M2,泊松比Υ03,密度Ρ24103KG/M3,离散网格如图5在厚度方向为一层,厚1M使用上述的刚体弹簧元方法计算悬臂梁的变形和应力并与有限元的结果进行对比,图6为悬臂梁沿Y轴的位移对比结果,图7和图8分别为X方向的正应力ΣX和剪应力ΤXY的对比结果从图中可以看出两种方法计算的位移和应力结果除有限元法在固定端能反应局部应力集中外,整体应力具有很好的近似性图6Y方向位移等值线单位M图7X方向正应力ΣX等值线单位MPA图8剪应力ΤXY等值线单位MPA刚体弹簧元通过引入单元表面应力分布假设和变形近似等效来模拟连续介质,利用这些假设,刚体弹簧元通过上述的刚度等效、强度等效和荷载等效,能够把连续介质和非连续介质统一在一个模型之中,为研究坝肩失稳导致拱坝断裂破坏需要仿真连续与非连续耦合介质的大变形与破坏断裂的复杂过程提供了有力的工具2玛尔帕塞拱坝失稳机理的研究21模型的建立211坝肩和坝体的离散和计算条件的设定法国玛尔帕塞拱坝坝高66M,坝顶高程10255M,坝顶长2227M,长高比33,坝顶厚15M,底厚678M,厚高比为01坝址岩体由带状片麻岩组成,含千节理倾向下游右岸,倾角30~50该坝于1954年建成,蓄水历时5年至1959年12月蓄至最高水位约100M时,坝在短瞬间溃决,左岸与中部坝体全部冲走,右岸底部坝体残存根据玛尔帕塞拱坝的地形和地质资料以及拱坝的体形和横缝设置将坝址处的岸坡简化成如图9所示,坐标系设置和拱坝离散网格如图10在分析中边界条件设定为除“二面沟”内的岩体外,其他岩体均被固定坝体、右岸坝肩和中间岩体依照连续介质进行模拟,坝体弹性模量取E1010PA,坝体材料的抗拉强度取20MPA,抗压强度为20MPA,摩擦系数F10;考虑到右岸坝基事后调查仍有坝体残存,表明该区岩体强度较高,故将坝体右岸及中间部分与基础的交界面强度较坝体提高1倍左岸坝肩依照非连续体计算,重力墩和地基之间的摩擦系数设定为10,坝体和左岸坝肩之间以及“二面沟”内的岩体两组节理面与水平层面取相同的摩擦系数F1,根据不同的工况选取,上、下游断层摩擦系数为F2,亦根据不同的工况进行选择计算中未考虑粘结强度C的作用主要荷载为库水作用,稳定分析中库水水位均取100M,考虑左坝肩扬压力时,按静水压力作用施加在沿上游面坝基线的裂隙处,作用面见图9中虚线段所示为了简化计算,所有荷载均一次施加计算中阻尼参照CUNDALL提供的阻尼模式选取刚度与质量成比例的RAYLEIGH阻尼型式,阻尼系数Α、Β由试算确定,试算表明如不出现过阻尼或阻尼过小情况计算结果均收敛于一稳态解的近似值,详见参考文献3212位移分析计算拱坝位移时左坝肩暂按连续介质考虑,并且不考虑扬压力作用,得到拱坝在1955年9月水位7975M蓄水至1959年7月水位9419M拱坝产生的位移,并与相应实测位移值进行比较,图11列出了顶拱和拱冠梁处的比较结果可以看出二者最大值基本接近,但在拱冠梁坝基附近,实测的位移值较计算大1CM左右,顶拱计算结果较实测结果略偏大这一比较结果表明选择的弹簧刚度参数KN、KS以及单元等效模型基本符合实际22玛尔帕塞拱坝稳定性分析及失稳过程的仿真通过改变坝体和左岸坝肩之间以及左岸坝肩“二面沟”内岩体的摩擦系数F1和上下游断层面的摩擦系数F2可以模拟和分析拱坝在各类参数条件下失稳破坏的情形图9岸坡离散

注意事项

本文(水利工程论文-拱坝-地基破坏的数值模型与溃坝仿真.doc)为本站会员(wenku)主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网(发送邮件至[email protected]或直接QQ联系客服),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。

关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们

网站客服QQ:2846424093    人人文库上传用户QQ群:460291265   

[email protected] 2016-2018  renrendoc.com 网站版权所有   南天在线技术支持

经营许可证编号:苏ICP备12009002号-5