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

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

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

水利工程论文-有限元强度折减系数法计算土坡稳定安全系数的精度研究.doc

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

水利工程论文-有限元强度折减系数法计算土坡稳定安全系数的精度研究.doc

水利工程论文有限元强度折减系数法计算土坡稳定安全系数的精度研究摘要有限元强度折减系数法在边坡稳定分析中的应用正逐渐受到人们的重视。本文较为全面地分析了土体屈服准则的种类、有限元法自身计算精度以及H坡高、β坡角、C粘聚力、Φ摩擦角对折减系数法计算精度的影响,并给出了提高计算精度的具体措施。通过对106个算例的比较分析,表明折减系数法所得稳定安全系数比简化Bishop法平均高出约5.7,且离散度极小,这不仅验证了文中所提措施的有效性,也说明了将折减系数法用于分析土质边坡稳定问题是可行的。关键词强度折减系数边坡稳定屈服准则误差分析自弗伦纽期于1927年提出圆弧滑动法以来,至今已出现数十种土坡稳定分析方法,有极限平衡法、极限分析法、有限元法等。不少研究表明,各种方法所得稳定安全系数都比较接近,可以说,这些方法已经达到了相当高的精度。近年来,由于计算机技术的长足发展,基于有限元的折减系数法在边坡稳定分析中的应用备受重视。与极限平衡法相比,它不需要任何假设,便能够自动地求得任意形状的临界滑移面以及对应的最小安全系数,同时它还可以真实的反映坡体失稳及塑性区的开展过程。到目前为止,已有很多学者对折减系数法进行了较为深入的研究1,2,3,并在一些算例中得到了与极限平衡法十分接近的结果。但总体说来,此法仍未在工程界得到确认和推广,究其原因在于影响该法计算精度的因素很多,除了有限元法引入的误差外,还依赖于所选用的屈服准则。此论文的目的有两点1力图全面分析屈服条件和有限元法本身对折减系数法计算精度的影响,并提出应选用何种屈服准则以及提高有限元法计算精度的具体措施2结合工程实例,分析对边坡稳定安全系数影响最大的4个主要参数H坡高、β坡角、C粘聚力、Φ摩擦角对折减系数法计算精度的影响。从以往的计算结果来看,严格法Spencer所得稳定安全系数比简化Bishop法平均高出约2~3,而通过106个算例的比较分析,表明折减系数法所得稳定安全系数比简化Bishop法平均高出约5.7,且误差离散度极小,可以认为是正确的解答4。这有力地说明了将有限元折减系数法用于分析土坡稳定问题是可行的,但必须合理地选用屈服条件以及严格地控制有限元法的计算精度,同时也表明有限元折减系数法所得安全系数稍微偏高,其原因有待进一步研究。1折减系数法的基本原理Bishop等将土坡稳定安全系数F定义为沿整个滑移面的抗剪强度与实际抗剪强度之比,工程中广为采用的各种极限平衡条分法便是以此来定义坡体稳定安全系数。有限元强度折减系数法的基本思想与此一致,两者均可称之为强度储备安全度。因后者无法直接用公式计算安全系数,而需根据某种破坏判据来判定系统是否进入极限平衡状态,这样不可避免地会带来一定的人为误差。尽管如此,仍发展了一些切实可行的平衡判据,如限定求解迭代次数,当超过限值仍未收敛则认为破坏发生或限定节点不平衡力与外荷载的比值大小或利用可视化技术,当广义剪应变等值线自坡角与坡顶贯通则定义坡体破坏3。文中平衡判据取当节点不平衡力与外荷载的比值大于103时便认为坡体破坏。有限元折减系数法的基本原理是将土体参数C、Φ值同时除以一个折减系数Ftrial,得到一组新的C′、Φ′值,然后作为新的材料参数带入有限元进行试算,当计算正好收敛时,也即Ftrial再稍大一些数量级一般为103,计算便不收敛,对应的Ftrial被称为坡体的最小安全系数,此时土体达到临界状态,发生剪切破坏,具体计算步骤可参考文献2,文中如无特别说明,计算结果均指达到临界状态时的折减系数。122屈服准则的影响用折减系数法求解实际边坡稳定问题时,通常将土体假设成理想弹塑性体,其中本构模型常选用摩尔库仑准则MC、DruckerPrager准则以及摩尔库仑等面积圆5准则。摩尔库仑准则可用不变量I1,J2,θσ表述成如下形式3Druckerprager4式中I1为应力张量第一不变量J2为应力偏量第二不变量θσ是应力洛德角。图1各屈服准则在π平面上的曲线MC准则较为可靠,它的缺点在于三维应力空间中的屈服面存在尖顶和棱角的不连续点,导致数值计算不收敛,所以有时也采用抹圆了的MC修正准则6,它是用光滑连续曲线来逼进摩尔库仑准则,此法虽然方便了数值计算,但不可避免地会引入一定的误差而DP准则在偏平面上是一个圆,更适合数值计算。通常取MC准则的外角点外接圆、内角点外接圆或其内切圆作为屈服准则,以利数值计算。各准则的参数换算关系见表1。由徐干成、郑颖人19905实际上是将MC准则转化成近似等效的DP准则形式。该准则要求偏平面上的摩尔库仑不等边六角形与DP圆面积相等。计算表明它与摩尔库仑准则十分接近。见图1,r1为外角外接圆半径r2为内角外接圆半径r3为内切圆半径摩尔库仑准则构成的六角形面积为5对半径为r的圆面积Sπr2,令SSmorl得67式7与式4对应项相等,可得8表1各准则参数换算编号准则种类αφkDP1外角点外接DP圆DP2内角点外接DP圆DP3内切DP圆DP4等面积DP圆注表中αφ、k是与DP有关的材料参数。表2不同屈服准则所得最小安全系数φ/°0.110253545DP10.5251.0441.7692.2543.051DP20.5250.9301.3321.5301.887DP30.4540.8481.2791.4991.870DP40.4770.8961.3961.6892.182简化Bishop法0.4940.8461.3161.6232.073DP1Bishop/Bishop0.0630.2340.3440.3550.472DP2Bishop/Bishop0.0630.0990.0120.0800.090DP3Bishop/Bishop0.0810.0020.0280.0990.098DP4Bishop/Bishop0.0340.0590.0610.0410.053注H20mmβ45°C42kPa。算例分析表明表2、图2DP4准则与简化Bishop法所得稳定安全系数最为接近。对有效算例Φ≠0的误差进行统计分析可知,当选用DP4准则时,误差的平均值为5.7,且离散度很小图3。而DP1的平均误差为29.5,同时采用DP2、DP3准则所得计算结果的离散度非常大,均不可用。因此在数值分析中可用DP4准则代替摩尔库仑准则。图2Φ~折减系数曲线图3DP4准则的计算误差3不同流动法则的影响有限元计算中,采用关联还是非关联流动法则,取决于ψ值剪胀角ψφ,为关联流动法则ψ≠0,为非关联流动法则。总体说来,采用非关联流动法则所得破坏荷载比同一类型材料而采用关联流动法则所得破坏荷载小,如忽略剪胀角ψ0,将会得到较为保守的结果。值得注意的是当ψ0时,正好与郑颖人等提出的广义塑性力学理论相符7,这时对应的塑性势面与q轴垂直。表3不同流动法则的影响φ10°φ17°φ25°非关联0.8711.1051.363关联0.8871.1371.425相对误差0.0180.0290.045β45°C40kPaH20mDP4准则。表4网格疏密对计算结果的影响节点数57711112250DP40.6610.6180.593简化Dishop法0.5830.5830.583DP4Bishop/Bishop0.1340.0600.017注H20mβ45°φ45°c10000Pa。笔者对采用不同流动法则的算例进行了初步分析,表3的计算结果表明对同一边坡,不论采用关联流动法则还是非关联流动法则,计算结果相差不大。这是因为它们只与坡体的体积变形有关,而在边坡稳定分析中,坡体常常为无约束天然坡体,体积变形对坡体稳定影响并不明显。然而,从破坏时位移大小及塑性区的分布来看,还是会有一些差异,有时并不能简单的忽略这种差异8。文中所有的算例均取ψ0,即满足非关联流动法则,算例结果显示出较好的精度。4有限元法引入的误差如前所述,本构模型的选择合理与否会对有限元折减系数法的计算精度造成较大影响,除此之外,有限元法本身也是误差的主要来源之一。4.1网格的疏密网格疏密对单元精度的影响甚至大于单元类型的影响,对于精度较低的单元,可通过加密网格来达到较高的精度。表4列出了不同疏密的网格对计算结果的影响,由表4可知,对于折减系数法,有限元网格不能太稀,否则结果将不

注意事项

本文(水利工程论文-有限元强度折减系数法计算土坡稳定安全系数的精度研究.doc)为本站会员(奋斗不息)主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网([email protected]),我们立即给予删除!

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

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