会员注册 | 登录 | 微信快捷登录 支付宝快捷登录 QQ登录 微博登录 | 帮助中心 人人文库renrendoc.com美如初恋!
站内搜索 百度文库

热门搜索: 直缝焊接机 矿井提升机 循环球式转向器图纸 机器人手爪发展史 管道机器人dwg 动平衡试验台设计

   首页 人人文库网 > 资源分类 > DOC文档下载

水利工程论文-含节理单元的三维P型自适应有限元解法.doc

  • 资源星级:
  • 资源大小:102.21KB   全文页数:15页
  • 资源格式: DOC        下载权限:注册会员/VIP会员
您还没有登陆,请先登录。登陆后即可下载此文档。
  合作网站登录: 微信快捷登录 支付宝快捷登录   QQ登录   微博登录
友情提示
2:本站资源不支持迅雷下载,请使用浏览器直接下载(不支持QQ浏览器)
3:本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

水利工程论文-含节理单元的三维P型自适应有限元解法.doc

水利工程论文含节理单元的三维P型自适应有限元解法摘要本文提出了含三维无厚度节理单元、等厚度节理单元和变厚度节理单元的p型自适应有限元模型,给出了三维节理单元升阶谱有限元法的解题步骤,通过具体算例,验证了p型升阶谱有限元法在求解含三维节理单元的有限元问题时的可行性及优越性。关键词节理单元升阶谱有限元p型有限元有限元法是求解微分方程数值解的一种重要方法,对于一个给定的问题,为改善其有限元解的精度,可以采用以下3种方法。1h型有限元法1,这种方法通过减小单元尺寸来提高有限元解的精度。2p型有限元法2,这种方法通过增加基底函数的阶次来提高有限元解的精度。3hp型有限元法3,这种方法是以上两种方法的综合,它既减小单元尺寸,又增加基底函数的阶次。作者所在的研究小组从1995年开始研究水工结构的h型弹粘塑性有限单元方法,目前已建立了实用的二维分析软件体系4,5,并在三维分析方面取得了进展6。从1999年开始,在水工结构的p型自适应分析方面也有所突破,1999年,程昭7等人针对水工结构分析问题提出了三维升阶谱有限元分析方法。2001年,陈胜宏8等人进一步提出了二维问题的p型自适应分析策略,并将自适应有限元方法归类为全域升阶方法、单元升阶方法和自由度升阶方法等三类。之后,费文平9等人将p型自适应有限元分析方法推广到三维弹粘塑性领域。但是,以上有关p型有限元的研究成果中均未涉及到断层、节理这一类特殊单元。大坝坝基、坝肩和岩石高边坡等部位总是存在断层、节理和软弱夹层等大规模的不连续面,且对结构的变形和稳定影响巨大,故在有限元分析中应给予高度的重视。古德曼Goodman最初运用有限元技术模拟岩体工程中的非线性不连续面问题,并提出了无厚度的节理元的概念10。随后,朱伯芳于1979年提出了等厚度节理元模型,并将其与无厚度节理元模型形成统一的计算公式11,在此基础上,王鸿儒等人提出了变厚度的节理单元的弹塑性模型并将其应用到工程实践中12。目前,国内外在有关p型自适应有限元分析的研究中,尚未涉及到这类特殊单元的处理问题,从而使研究成果的工程应用受到一定程度的限制。对于常规块体单元的三维p型有限元模型,作者在文献9中已有详尽的论述,本文主要给出三维无厚度节理单元、等厚度节理单元和变厚度节理单元的p型有限元模型,并给了升阶谱的计算格式。实例分析结果表明,用p型有限元法来求解含三维节理单元的有限元问题具有收敛速度快、计算精度高的优点。1三维p型无厚度节理单元和等厚度节理单元模型对厚度很小和厚度变化不大的节理,可以分别采用无厚度的节理单元和等厚度的节理单元进行模拟,可取如图1所示的六面体节理单元,进行单元的网格划分。1.1形函数及位移函数对节理单元上下面的相应点、棱和面可取相同的点基函数、棱基函数和面基函数,对无厚度节理单元或等厚度节理单元,不存在连续上、下两面的棱基函数和面基函数,也不存在体基函数。基底函数的具体形式如下13点基函数p≥11式中,。棱基函数p≥2,2面基函数p≥4(ijp,i,j≥2)3式中而为Legender多项式。令位移函数为456同理可以写出V下,V上,W下,W上及Δv,Δw的具体表达式。将基函数Ni,pEi,pF统一记为φi,位移差uN,i4uNi,uE,i4uEi,uF2uF1ui,设单元基底函数个数为fep,上式简化为,同理有,1.2三维等厚度节理单元或无厚度节理单元升阶过程当p1时N1=φ1Iφ2Iφ3Iφ4Iφ1Iφ2Iφ3Iφ4I,I为33阶的单位阵。当p2时,N2可在N1的基础上进行扩充,扩充矩阵ΔNp2为ΔNp2φ5Iφ6Iφ7Iφ8Iφ5Iφ6Iφ7Iφ8I。依此类推,可得最终的形函数矩阵为71.3坐标插值及坐标变换对于节理上、下面坐标的插值仍采用各面上的四个节点进行插值,即89式中xi,yi,zii1~8为节理间面体单元8个顶点的整体坐标。定义整体坐标系的x轴朝北,y轴朝西,z轴朝上,定义等厚度的节理单元或无厚度的节理的局部坐标系的z′为中面的法线朝上方向,y′指向节理面的倾向,x′轴由右手法则确定,并设等厚度节理的倾角为α,倾向为β。三维等厚节理或无厚节理单元局部坐标与整体坐标的转换矩阵为101.4三维等厚度节理的单元刚度矩阵1112式中弹性矩阵单元应变矩阵B=LB=1/eLNp根据虚功原理,单元刚度矩阵为131.5三维无厚度节理的单元刚度矩阵当等厚度节理单元的厚度e→0时,即形成无厚度的节理单元。此时,可假定单元内应力分量与位移差成正比,同理可得单元刚度矩阵为14式中λ′为单元劲度矩阵。2三维p型变厚度的节理单元模型当节理单元的厚度变化较大时,应将等厚度节理单元推广得到变厚度的节理单元,如图2所示建立局部坐标系。2.1形函数及位移函数六面体变厚度的节理单元的基底函数,由点基函数、棱基函数、面基函数和体基函数组成,各基底函数的具体形式如下点基函数p≥1Niξ,η,ζ1/81+ξiξ1ηiη1ζiζi1,2,,815式中ξi1i,ηi1i/20.5,ζi1i/40.75。棱基函数p≥2pEi1/41ξ1iξ1η1iη1ζ1iζξ2iΦpξη2iΦpηζ2iΦpζ,i,1,2,1216式中ξ2i1i/120.65,η2i=1211i1/4ζ2ii18,将ξ1i,η1i,ζ1i用向量的形式表达为{ξ1}{000011111111}T,{η1}={11110000111}T,{ζ1}={111111110000}T。面基函数p≥4171819

注意事项

本文(水利工程论文-含节理单元的三维P型自适应有限元解法.doc)为本站会员(奋斗不息)主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网([email protected]),我们立即给予删除!

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

copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5