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

收藏 分享

资源预览需要最新版本的Flash Player支持。
您尚未安装或版本过低,建议您

水利工程论文有限元强度折减系数法计算土坡稳定安全系数的精度研究摘要有限元强度折减系数法在边坡稳定分析中的应用正逐渐受到人们的重视。本文较为全面地分析了土体屈服准则的种类、有限元法自身计算精度以及H坡高、Β坡角、C粘聚力、Φ摩擦角对折减系数法计算精度的影响,并给出了提高计算精度的具体措施。通过对106个算例的比较分析,表明折减系数法所得稳定安全系数比简化BISHOP法平均高出约57,且离散度极小,这不仅验证了文中所提措施的有效性,也说明了将折减系数法用于分析土质边坡稳定问题是可行的。关键词强度折减系数边坡稳定屈服准则误差分析自弗伦纽期于1927年提出圆弧滑动法以来,至今已出现数十种土坡稳定分析方法,有极限平衡法、极限分析法、有限元法等。不少研究表明,各种方法所得稳定安全系数都比较接近,可以说,这些方法已经达到了相当高的精度。近年来,由于计算机技术的长足发展,基于有限元的折减系数法在边坡稳定分析中的应用备受重视。与极限平衡法相比,它不需要任何假设,便能够自动地求得任意形状的临界滑移面以及对应的最小安全系数,同时它还可以真实的反映坡体失稳及塑性区的开展过程。到目前为止,已有很多学者对折减系数法进行了较为深入的研究1,2,3,并在一些算例中得到了与极限平衡法十分接近的结果。但总体说来,此法仍未在工程界得到确认和推广,究其原因在于影响该法计算精度的因素很多,除了有限元法引入的误差外,还依赖于所选用的屈服准则。此论文的目的有两点1力图全面分析屈服条件和有限元法本身对折减系数法计算精度的影响,并提出应选用何种屈服准则以及提高有限元法计算精度的具体措施;2结合工程实例,分析对边坡稳定安全系数影响最大的4个主要参数H坡高、Β坡角、C粘聚力、Φ摩擦角对折减系数法计算精度的影响。从以往的计算结果来看,严格法SPENCER所得稳定安全系数比简化BISHOP法平均高出约2~3,而通过106个算例的比较分析,表明折减系数法所得稳定安全系数比简化BISHOP法平均高出约57,且误差离散度极小,可以认为是正确的解答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不同屈服准则所得最小安全系数Φ/0110253545DP105251044176922543051DP205250930133215301887DP304540848127914991870DP404770896139616892182简化BISHOP法04940846131616232073DP1BISHOP/BISHOP00630234034403550472DP2BISHOP/BISHOP00630099001200800090DP3BISHOP/BISHOP00810002002800990098DP4BISHOP/BISHOP00340059006100410053注H20MM;Β45;C42KPA。算例分析表明表2、图2DP4准则与简化BISHOP法所得稳定安全系数最为接近。对有效算例Φ≠0的误差进行统计分析可知,当选用DP4准则时,误差的平均值为57,且离散度很小图3。而DP1的平均误差为295,同时采用DP2、DP3准则所得计算结果的离散度非常大,均不可用。因此在数值分析中可用DP4准则代替摩尔库仑准则。图2Φ~折减系数曲线图3DP4准则的计算误差3不同流动法则的影响有限元计算中,采用关联还是非关联流动法则,取决于Ψ值剪胀角ΨΦ,为关联流动法则;Ψ≠0,为非关联流动法则。总体说来,采用非关联流动法则所得破坏荷载比同一类型材料而采用关联流动法则所得破坏荷载小,如忽略剪胀角Ψ0,将会得到较为保守的结果。值得注意的是当Ψ0时,正好与郑颖人等提出的广义塑性力学理论相符7,这时对应的塑性势面与Q轴垂直。表3不同流动法则的影响Φ10Φ17Φ25非关联087111051363关联088711371425相对误差001800290045Β45;C40KPA;H20M;DP4准则。表4网格疏密对计算结果的影响节点数57711112250DP4066106180593简化DISHOP法058305830583DP4BISHOP/BISHOP013400600017注H20M;Β45;Φ45;C10000PA。笔者对采用不同流动法则的算例进行了初步分析,表3的计算结果表明对同一边坡,不论采用关联流动法则还是非关联流动法则,计算结果相差不大。这是因为它们只与坡体的体积变形有关,而在边坡稳定分析中,坡体常常为无约束天然坡体,体积变形对坡体稳定影响并不明显。然而,从破坏时位移大小及塑性区的分布来看,还是会有一些差异,有时并不能简单的忽略这种差异8。文中所有的算例均取Ψ0,即满足非关联流动法则,算例结果显示出较好的精度。4有限元法引入的误差如前所述,本构模型的选择合理与否会对有限元折减系数法的计算精度造成较大影响,除此之外,有限元法本身也是误差的主要来源之一。41网格的疏密网格疏密对单元精度的影响甚至大于单元类型的影响,对于精度较低的单元,可通过加密网格来达到较高的精度。表4列出了不同疏密的网格对计算结果的影响,由表4可知,对于折减系数法,有限元网格不能太稀,否则结果将不
编号:201312122203085723    类型:共享资源    大小:114.55KB    格式:DOC    上传时间:2013-12-12
  
2
关 键 词:
专业文献 建筑水利 精品文档 水利工程
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
关于本文
本文标题:水利工程论文-有限元强度折减系数法计算土坡稳定安全系数的精度研究.doc
链接地址:http://www.renrendoc.com/p-205723.html

当前资源信息

4.0
 
(2人评价)
浏览:14次
wenku上传于2013-12-12

官方联系方式

客服手机:17625900360   
2:不支持迅雷下载,请使用浏览器下载   
3:不支持QQ浏览器下载,请用其他浏览器   
4:下载后的文档和图纸-无水印   
5:文档经过压缩,下载后原文更清晰   

精品推荐

相关阅读

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

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

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

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