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

水利工程论文-小湾高拱坝坝踵开裂的有限单元法分析.doc

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

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

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

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

水利工程论文-小湾高拱坝坝踵开裂的有限单元法分析.doc

水利工程论文小湾高拱坝坝踵开裂的有限单元法分析摘要在开裂单元中引入不连续形函数,然后基于虚功原理推导出开裂单元的广义刚度矩阵和广义荷载的具体表达式,以及有限元平衡方程,从而建立了一种三维裂缝扩展的不变网格有限元分析方法。运用该方法对已建的二滩拱坝和待建的小湾工程拱坝在相同的工况下进行了对比计算。计算结果表明在目前所考虑的影响因素下,小湾拱坝和二滩拱坝的开裂范围、开裂程度及对网格疏密、混凝土抗拉强度和水力劈裂等影响因素的敏感程度都非常相似。考虑到二滩工程已建成,且运行正常,可以初步判定小湾工程的拱坝坝踵开裂危险性不大。关键词拱坝开裂有限单元法坝踵开裂是高拱坝设计中需要考虑的一个重要问题。奥地利的科恩布赖拱坝建于1977年,坝高2002M,在蓄水前坝下游底部水平施工缝张开,1978年蓄水接近满库时上游坝踵开裂穿通廊道,引起廊道内漏水,后来进行了十几年的漫长修补加固,并在坝的下游坝趾处建造了高约70M的混凝土支撑体。俄罗斯的萨扬舒申斯克拱坝建于1987年,坝高2420M,施工期产生大量裂缝,当1990年水位首次达到正常水位时,河床坝段上游坝面靠近坝基处产生水平裂缝,引起廊道内漏水。这些高拱坝的损坏事故引起了世界各国坝工界的极大关注。研究拱坝开裂分析的主要方法可以分为结构模型试验和数值计算两大类。可用于拱坝开裂分析的数值计算方法很多,如拱梁分载法、边界元法、流形元法、无单元法、有限单元法等1~7,其中有限元法是较为成熟且广为应用的数值方法。但有限元法在本质上是一种连续介质力学方法,必须加以改进才能被运用来分析裂缝扩展这样的不连续问题。改进的方法可以概括为两类变网格法和固定网格法。变网格法随着裂缝的扩展而重新划分网格8,9,这种方法的优点是不需要开发新的有限元软件,但是存在计算量大和前处理难等问题。固定网格法则是保持有限元网格不变,通过修改开裂单元的插值关系和本构关系来反映裂缝的影响10~12。相比较而言,固定网格法的应用更为方便。本文在开裂单元中引入不连续形函数,然后基于虚功原理推导出开裂单元的广义刚度矩阵和广义荷载的具体表达式,以及有限元平衡方程,从而建立了一种三维裂缝扩展的不变网格有限元分析方法。本文方法在小湾工程拱坝的分析中获得了应用。由于拱坝开裂的影响因素复杂,计算方法也多种多样,很难做出完全定量的分析评判,且目前也没有相应的开裂深度安全准则,故研究中采用“校准法”的思路取一个与小湾工程相类似的已建且运转正常的工程,在相同的工况下,采用同一计算方法进行对比分析,通过开裂深度与坝底宽度比值的比较,评价小湾高拱坝的开裂危险性和安全度。根据设计部门的建议,作为校准的已建工程选为二滩工程。1三维裂缝扩展的不变网格有限单元分析方法11开裂单元的位移模式12常规有限单元法中,单元位移场{U}可以表示为{U}N{U}E1图1开裂单元示意式中N和{U}E分别为单元形函数矩阵和单元结点位移向量。式1表示的位移场在单元内是连续的,不能描述单元开裂以后沿开裂面的不连续位移模式。为求改进,需构造不连续位移场插值函数。设单元E被开裂面分为E和E两部分图1,它们的位移场分别记为{U}和{U},可以利用原单元E的形函数矩阵表示为{U}N{A}E{U}N{B}E2式中{A}E和{B}E是单元结点的广义位移向量。令3则{A}E{C}E{D}E;{B}E{C}E{D}E4将式4代入式2可得{U}N{C}EN{D}E;{U}N{C}EN{D}E5或{U}N{C}ENHX{D}E6式中NHX称为不连续形函数矩阵,其中HX是一个不连续函数,定义如下7式6右端的第一项与式1的右端项在形式上是一致的,表示的是单元中的连续位移场,而第二项表示的则是由于裂缝的出现而附加的不连续位移场。根据式6,开裂单元的位移模式是在连续单元位移模式的基础上,增加了反映不连续位移的插值项。这些附加的插值项与单元结点对应,但其形函数不连续,故称之为开裂单元的广义结点。开裂单元的应变可以表示为{Ε}B{C}EBHX{D}E8式中B是应变矩阵。应力可以表示为{Σ}D{Ε}DB{C}EDBHX{D}E9式中D是弹性矩阵,它是单元材料的弹性模量和泊松比的函数。引入不连续形函数和广义结点的概念后,开裂单元的位移模式和应变模式都可以看作是在连续单元的模式上增加了广义结点的影响项,从而使得固定网格下的裂缝扩展分析变得简便易行。12有限元平衡方程不连续单元的虚功原理可以表示为10式中V和S分别表示单元的体积和边界包括裂缝面;{QV}和{QS}分别是单元的体积力和边界面力包括裂缝面上的分布力。将式6、式8和式9代入式10,整理可得K{U}E{F}E11121314式中K称为单元广义刚度矩阵,{U}E称为单元结点广义位移向量,{F}E称为单元结点广义荷载向量。K由四种类型的子矩阵组成,它们的表达式分别为1516{F}E由两部分组成,分别为1718其中式15和式17在形式上与常规有限元法是一致的。考虑Ⅰ型裂缝的情况,且直接采用裂缝前沿附近单元高斯点的应力,根据最大拉应力准则判断裂缝是否扩展,并且假定开裂单元被裂缝面贯穿。开裂分析中裂缝的扩展将产生应力释放,可以按下式计算其等效荷载19式中的积分沿裂缝面进行。13水力劈裂作用的模拟在拱坝开裂计算中究竟如何考虑水力劈裂的作用是一个复杂的问题。随着裂缝的扩展,水流渗入裂缝,在裂缝面上产生劈裂作用,同时也引起渗流场的改变。如果考虑渗流场与裂缝扩展的耦合作用,将会使计算过程变得非常复杂,甚至可能无法实现。因此,宜对水力劈裂的作用进行适当的简化。对于高拱坝坝踵的张裂缝,可以不考虑水流渗入裂缝引起的渗流场改变,而只考虑水流在裂缝面的劈裂作用,作用力大小近似取上游全水头压强。这样简化是偏于安全的。以上算法已在有限单元法程序CORE3实现,并且通过算例考核证实了其合理性12。2小湾和二滩拱坝应用实例21工程概况小湾水电站位于云南省大理州的南涧县和临沧地区凤庆县的交界,是澜沧江中下游河段规划8个梯级电站中的第二级,电站装机4200MW,水库库容15132108M3。大坝为混凝土抛物线双曲拱坝,最大坝高292M,坝顶高程1245M,是即将开工的世界第一高拱坝。小湾拱坝拱冠梁顶部宽度12M,拱冠梁底部宽度7291M,厚高比为0250。二滩水电站位于四川省攀枝花市境内,雅砻江下游,电站装机3300MW,年发电170亿度。二滩大坝为混凝土抛物线双曲拱坝,最大坝高240M,坝顶高程1025M,是已建的中国第一高拱坝。二滩拱坝拱冠梁顶部宽度11M,拱冠梁底部宽度5574M,厚高比为0232。22计算条件高拱坝在坝踵、坝趾和基础高应力区会出现应力集中,并在坝踵和部分基础内产生裂缝。随着裂缝的开展,库水会进入裂缝,产生水力劈裂作用,从而促使裂缝进一步扩展。如果坝体与地基连接很好,接触面具有很高的强度,开裂可能会深入地基,但是这种裂缝往往属于稳定扩展,对大坝不会构成太大的危害。一般情况下,坝体与地基的交接面附近为一弱面,此面在靠近上游坝踵处往往受一法向拉应力,从而导致裂缝沿坝体与地基的交接面或稍高处扩展,此种裂缝对拱坝安全的危胁最大。本研究将以沿坝体与地基交接面或稍高部位的裂缝扩展作为重点。221计算模型计算模型的建立对三维裂缝问题的求解有显著的影响,网格划分密实,单元形状相对正态,能够充分提高求解精度。为了分析不同计算模型对求解精度的影响,并对小湾拱坝和二滩拱坝进行分析比较,共建立了四套有限元模型,其中小湾拱坝粗网格模型和精细网格模型各一套,二滩拱坝粗网格模型和精细网格模型各一套。图2小湾拱坝精细网格32024单元图3二滩拱坝精细网格20876单元222计算工况分别对考虑水力劈裂作用和不考虑水力劈裂作用两种情况进行了计算。考虑水力劈裂作用时,在已开裂缝面上施加上游全水头。此外,对坝体混凝土的单轴抗拉强度分别取0MPA、1MPA和2MPA进行敏感性分析。以上情况共组合成24种工况。23计算结论图4和图5分别为小湾拱坝和二滩拱坝建基面裂缝扩展范围示意图混凝土抗拉强度1MPA,考虑水力劈裂,图6和图7分别为在相同工况下小湾拱坝和二滩拱坝拱冠梁断面开裂后应力矢量图。图4小湾拱坝建基面裂缝扩展示意混凝土抗拉强度1MPA,考虑水力劈裂图5二滩拱坝建基面裂缝扩展示意混凝土抗拉强度1MPA,考虑水力劈裂

注意事项

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

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

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

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

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

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