有限元强度折减法
- 格式:docx
- 大小:1.73 MB
- 文档页数:11
基于有限元强度折减法的滑坡稳定性分析滑坡稳定性分析基于有限元强度折减法是一种用于确定滑坡极限稳定性的重要方法。
它主要是通过在滑坡稳定性分析中应用有限元强度折减法,以折减破坏面的形状,计算滑坡受力情况,以及滑坡自重,物质特性及岩土的摩擦特性的数值计算,最终用分析结果来判断滑坡稳定发展的可能性,以确定滑坡稳定状态。
一、有限元强度折减法1、折减原理:有限元强度折减法是一种直接定位破坏面的方法,其原理是通过折减岩体的强度,来确定破坏开裂的面。
在有限元中,折减的本质就是改变模型的材料参数,找到一个最小的一组有限元强度折减设定,以便确定所需的破坏面。
2、折减边界:有限元强度折减法的折减边界就是要折减的破坏开裂的面。
尽管可以采用自然边界,但是最好采用与实际条件有关的先进边界。
二、滑坡受力情况1、岩土特征:滑坡稳定分析包括对岩土特性的计算,例如土壤材料的屈服强度、弹性模量和泊松比以及岩土体内强度、摩擦以及连接情况等,并结合岩土稳定性理论,评价土坡稳定性。
2、受力、物质特性:另外,还需要考虑滑坡体的受力和物质特性,这些元素包含滑坡自重、坡面上的重力、地形力以及雨水等,它们也是滑坡稳定性分析的重要组成部分。
三、岩土的摩擦特性1、析出摩擦角:在滑坡稳定性分析中,析出岩土的摩擦角是计算极限稳定性的重要标准之一。
通过有限元强度折减法分析,可以精准计算出滑体内岩土摩擦角,从而得到表征滑坡发展可能性的结果。
2、摩擦和静定:岩土的摩擦力可以通过契约定理分析求得,它是由滑体摩擦角和坡度决定的,其大小可以被表达为“摩擦-坡度”系数。
此外,只有当滑体内岩土摩擦角足够大时,滑坡才具有静定发展的可能性。
四、滑坡稳定状态1、岩体状态:滑坡稳定状态可以根据岩体状态来评价,只有当滑坡稳定发展时,才能保证滑坡体状态稳定;2、计算结果:通过有限元强度折减法分析,可以根据折减的结果计算出滑体的受力状况,确定极限稳定性;3、应变计算:此外,还需要通过应变计算和时变分析,来评价滑坡稳定状态的发展趋势。
有限元强度折减法应用的几个问题及拓展有限元强度折减法是一种应用广泛的结构设计方法。
然而,在实际应用中,还存在一些问题需要解决。
本文将探讨有限元强度折减法应用中的以下几个问题:
1. 强度折减系数的确定方法:强度折减系数是有限元强度折减法的核心参数。
不同的结构及不同的工况下,强度折减系数的确定方法也不尽相同。
本文将介绍几种常见的强度折减系数确定方法,并探讨其优缺点。
2. 材料的非线性效应:很多结构在工作过程中会发生非线性变形,这会对有限元强度折减法的应用产生影响。
本文将介绍如何考虑材料的非线性效应,并探讨其影响。
3. 极限状态的确定方法:有限元强度折减法基于极限状态设计理论,而极限状态的确定方法对于结构设计具有至关重要的影响。
本文将介绍几种常见的极限状态确定方法,并探讨其优缺点。
4. 多学科设计优化中的应用:有限元强度折减法在多学科设计优化中的应用也是一个重要的研究方向。
本文将探讨有限元强度折减法在多学科设计优化中的应用,并介绍相关的研究进展。
总之,有限元强度折减法是一种重要的结构设计方法,但在实际应用中还存在一些问题需要解决。
本文将介绍这些问题,并提出一些拓展方向,以期为有限元强度折减法的应用和发展提供帮助。
- 1 -。
1 / 1
土力学
S o i l M e c h a n i c s
第七章 Chapter 7
疑难释义7.3
有限元强度折减法
有限元强度折减法是边坡工程中最常用的数值分析方法之一,它的主要原理就是通过不断降低边坡土体的抗剪强度参数,直到边坡达到极限破坏状态,自动形成滑裂面为止。
抗剪强度折减系数是指在外荷载保持不变的情况下,边坡内土体所发挥的最其定义与边坡稳定性分根据莫尔-=tan c τσϕ+
'c c f =,
tan 'ϕ=此时,土体的抗剪强度满足
tan ='tan '=c c f f ϕτσϕ++
计算时,不断的增加折减系数f ,使土体的抗剪强度不断减小,直到f 达到某一个值时,有限元计算不收敛,土体发生失稳,此时的强度折减系数f 即为土坡的安全稳定系数K 。
有限元强度折减系数法计算土坡稳定安全系数的精度研究[ 06-02-22 13:31:00 ] 作者:张鲁渝1,郑颖人1,编辑:studa9ngns摘要:有限元强度折减系数法在边坡稳定分析中的应用正逐渐受到人们的重视。
本文较为全面地分析了土体屈服准则的种类、有限元法自身计算精度以及H(坡高)、β(坡角)、C(粘聚力)、Φ(摩擦角)对折减系数法计算精度的影响,并给出了提高计算精度的具体措施。
通过对106个算例的比较分析,表明折减系数法所得稳定安全系数比简化Bishop法平均高出约5.7%,且离散度极小,这不仅验证了文中所提措施的有效性,也说明了将折减系数法用于分析土质边坡稳定问题是可行的。
关键词:强度折减系数边坡稳定屈服准则误差分析自弗伦纽期于1927年提出圆弧滑动法以来,至今已出现数十种土坡稳定分析方法,有极限平衡法、极限分析法、有限元法等。
不少研究表明,各种方法所得稳定安全系数都比较接近,可以说,这些方法已经达到了相当高的精度。
近年来,由于计算机技术的长足发展,基于有限元的折减系数法在边坡稳定分析中的应用备受重视。
与极限平衡法相比,它不需要任何假设,便能够自动地求得任意形状的临界滑移面以及对应的最小安全系数,同时它还可以真实的反映坡体失稳及塑性区的开展过程。
到目前为止,已有很多学者对折减系数法进行了较为深入的研究[1,2,3],并在一些算例中得到了与极限平衡法十分接近的结果。
但总体说来,此法仍未在工程界得到确认和推广,究其原因在于影响该法计算精度的因素很多,除了有限元法引入的误差外,还依赖于所选用的屈服准则。
此论文的目的有两点:(1)力图全面分析屈服条件和有限元法本身对折减系数法计算精度的影响,并提出应选用何种屈服准则以及提高有限元法计算精度的具体措施;(2)结合工程实例,分析对边坡稳定安全系数影响最大的4个主要参数(H坡高、β坡角、C粘聚力、Φ摩擦角)对折减系数法计算精度的影响。
从以往的计算结果来看,严格法(Spencer)所得稳定安全系数比简化Bishop法平均高出约2%~3%,而通过106个算例的比较分析,表明:折减系数法所得稳定安全系数比简化Bishop法平均高出约5.7%,且误差离散度极小,可以认为是正确的解答[4]。
有限元强度折减法应用的几个问题及拓展引言有限元强度折减法是一种广泛应用的结构抗震容量评估方法,适用于各种类型的建筑结构。
随着理论和实践研究的深入,该方法的理论和应用更加丰富和复杂。
但是,在使用过程中还存在一些问题和待解决的难题。
本文将探讨并提出几个有限元强度折减法应用的问题及拓展。
问题一:有限元强度折减法参数的选择在有限元强度折减法中,参数的选择对结果的准确性和可靠性有很大的影响。
首先,选择强度折减系数M,它是考虑结构的非线性行为和损伤效应的量化描述。
其值的大小直接影响到结构的强度和刚度。
其次,选择分段因子α,分段因子α是用于分段分析时控制分段长度的系数。
同时,还要对分段分析方法进行选择,如求解器类型和分析方法等等。
正确的参数设定将有助于提高有限元强度折减法的准确性和可靠性。
问题二:有限元强度折减法的基本假设在使用有限元强度折减法时,其基本假设是结构的强度退化是因为材料的破坏和梁数据不可靠而发生的。
然而,在实际结构中,强度退化不仅仅是因为这些原因,还有可能与结构的几何形状、结构组成、结构固有特性、结构的工艺技术等因素有关。
因此,准确评估强度退化的原因,对于有限元强度折减法在实际工程中的应用非常重要。
问题三:有限元强度折减法应用时的计算时间和计算精度有限元强度折减法常常需要进行大量的模拟计算,计算时间十分耗费。
同时,随着结构的复杂程度增加,在有限时间内获得高精度的计算结果也面临着很大的挑战。
因此,开发有效的计算方法和提高计算机性能将有助于提高有限元强度折减法的计算精度和计算效率。
问题四:有限元强度折减法的拓展在有限元强度折减法的应用中,其发展有三个方向:其一是将其用于非线性动力学分析、其二是将其扩展到考虑非线性几何行为的分析(如大位移、大变形、逆摆角度等)以及非线性材料行为的分析(如混凝土开裂、骨料断裂等),其三是增加其他结构的性能评估指标,例如板剪强度、面层弯曲刚度等。
结论有限元强度折减法是一种经过实践验证的有效方法,在实际工程中有广泛的应用前景。
有限元强度折减法在边坡稳定计算中的应用近年来,随着科学技术的发展,边坡安全工程成为当今社会的热点问题。
边坡稳定性是边坡工程安全性的重要水平指标之一。
有效地确定边坡稳定性,可以减少边坡垮塌,破坏性侵蚀及其他地质灾害的发生,更有效地保护人民生命财产安全。
传统的边坡稳定计算方法有很多缺点,很难解决大规模复杂的边坡稳定计算问题。
为了解决这些问题,随着计算机技术的发展,数值计算技术的发展,边坡稳定计算中有限元强度折减法也逐渐得到应用,它为边坡稳定计算提供了一种有效的方法。
有限元强度折减法是由著名有限元数值计算理论家Hans Zienkiewicz 于一九六三年提出的,它把数值计算分解为两个步骤:有限元分析和强度折减。
有限元分析不仅可精确的计算土体的应力和变形,还可以求得边坡的稳定系数。
而强度折减步骤则是对这些应力值和变形按照一定的准则进行折减,从而实现边坡的稳定计算。
有限元强度折减法在边坡稳定计算中的具体应用,有很多研究者提出了不同的算法,表达了不同的稳定准则和折减准则。
其中,以倒角条件准则、惯性假设准则、根据节理、滞回因子、失稳指标等为稳定准则的稳定计算模型。
具体的折减准则有力学强度折减法、抗剪强度折减法、抗压强度折减法、有效应变折减法等。
有限元强度折减法在边坡稳定计算中,可以有效地求解复杂参数边坡的稳定性,它把不同的计算模型、稳定准则和折减准则整合在一起,使边坡稳定计算更精确、更准确、更实用。
有限元强度折减法在边坡稳定计算中的应用已经得到了广泛的应用,它可以有效地求解复杂的边坡的稳定性,可以有效地减少不必要的垮塌、破坏性侵蚀等灾害,并可有效地保护人民生命财产安全。
因此,有限元强度折减法在边坡稳定计算中的应用,具有重要的理论意义和实际意义。
综上所述,有限元强度折减法在边坡稳定计算中的应用具有重要的理论意义和实际意义,它可以有效地求解复杂的边坡的稳定性,可以有效地减少不必要的垮塌、破坏性侵蚀等灾害,并可有效地保护人民的生命财产安全。
有限元强度折减法在岩土工程中的应用有限元强度折减法是一种近年来在岩土工程中广泛应用的方法。
该方法基于有限元理论和强度理论,通过对材料本构关系和强度减退规律的描述,可以进行岩土钻井、基坑支护等工程的稳定性分析。
该方法的优点是具有高精度、高可靠性、高效性等特点,不仅可以判断工程在不同荷载条件下的稳定性,还可以为工程设计提供有效的指导。
有限元强度折减法的基本原理是将材料的强度按一定比例进行折减,在外力的作用下,利用强度折减模型描述材料的强度减退规律,从而得到材料的局部损伤情况,再通过有限元分析对工程进行研究。
在进行有限元强度折减法分析时,首先需要建立材料本构模型和强度减退模型。
在建模的过程中,需要考虑材料的非线性和各向异性,同时由于现代岩土工程通常涉及到大变形和断裂,建模还需要考虑这些因素对稳定性的影响。
除此之外,还需要考虑材料的强度损伤规律,并在模型中体现出来。
有限元强度折减法主要在岩土钻井、基坑支护等工程中得到了广泛的应用。
在岩土钻井中,该方法可以模拟出钻头和钻孔过程中的应变和应力变化,进而分析岩土的破坏机制。
在基坑支护工程中,该方法可以对支撑结构的力学性能进行评估,分析支撑结构的破坏模式和破坏机制。
由于该方法具有高精度、高可靠性、高效性等优点,可以为工程设计和施工提供有效的指导,对于保障工程的安全和稳定性具有重要的作用。
总之,有限元强度折减法是一种应用广泛的岩土工程分析方法。
该方法通过建立合适的材料本构模型和强度减退模型,可以对工程的稳定性进行高精度、高可靠性的分析和评估,为工程设计和施工提供有效的指导。
该方法的应用将有助于提高岩土工程的施工质量和安全水平,推动岩土工程技术的发展和进步。
有限元强度折减法1 背景1974年,Smith & Hobbs[1]使用有限元方法分析了φu=0条件下的边坡稳定性并与Taylar[2]的结果进行对比,得到了很好的一致性;1975年,Zienkiewicz等[3]考虑c’、φ’进行有限元边坡稳定性分析,其结果与圆弧滑面解有较好吻合;1980年Griffiths[4]验证了一系列具有不同材料特性和形状的边坡稳定性并通过与Bishop& Morgenstern[5]的结果进行了对比确定了数据的可靠性;此后也有研究证实了利用有限元方法进行边坡稳定性分析的可靠性[6,7,8,9];在文献[9]中,引入一些案例证明了有限元强度折减法的准确性,并证明了有限元强度折减法在分析非均质边坡时相对于传统方法的优越性。
2001年,郑颖人等[10]把有限元强度折减法引入国内,并对此进行了后续研究[11,12,13,14]。
相较于一些传统的边坡稳定型分析方法,有限元强度折减法有以下几个优点[9]:(1)不必假设滑面的位置和形状,当土体自身强度不足以抵抗剪应力时土体失稳会自然发生。
(2)由于有限元强度折减法中没有条分的概念,因此也不必假设条间力,在整体失稳之前土体都处于整体稳定状态。
(3)使用有限元方法能够查看破坏过程。
2 有限元强度系数折减法1.模型参数边坡模型主要包括六个参数,分别是:膨胀角ψ、内摩擦角φ’、黏聚力c’、弹性模量E’、泊松比υ’、重度γ。
膨胀角影响土体屈服后的体积变形,若ψ<0,则土体屈服后体积减小,若ψ>0则体积增大,ψ=0则体积不变。
ψ=φ的情况被称之为关联流动法则,但是此时ψ值通常高于实验观测值,特别是在侧限条件下会提高土的承载力预测值。
边坡稳定型问题通常是处于无侧限条件下,此时膨胀角的选取不再重要[9],因此文献[9]选取ψ=0条件下的非关联流动法则,并且通过案例分析可以得出此膨胀角的选取可以得出准确的安全系数以及滑动面。
有限元强度折减法在边坡稳定计算中的应用有限元强度折减法是一种求解复杂结构力学问题的新方法,用于分析边坡稳定性也有强大的能力。
最近,有关使用有限元强度折减法的研究取得了巨大的进展,在计算边坡稳定性时取得了良好的结果。
本文将就有限元强度折减法在边坡稳定计算中应用的可行性及效果作一介绍。
一、有限元强度折减法背景1、有限元强度折减法是什么?有限元强度折减法是通过改变单元的材料参数,使得最终近似解与实际T失效状态一致,达到分析结构安全性能的一种计算方法。
这一计算方法能够较好地反映出结构的失效过程,从而改善传统的有限元算法的拟合不足的问题。
2、有限元强度折减法的特点有限元强度折减法不仅考虑结构的失效过程,还可以继而模拟出材料的弹性和变形过程,从而改进传统的有限元算法的拟合不足的问题。
此外,它还能模拟多种类型的变形,以保证结构承受能力及临界状态分析。
二、有限元强度折减法在边坡稳定计算中的应用1.计算边坡稳定时的精确性:使用有限元强度折减法计算边坡稳定性能,可以反映出坡面弹性及变形特性,从而更准确地评估边坡的变形性能。
2.降低计算时间:有限元强度折减法可以快速精确地计算边坡稳定性,因此在减少计算时间的同时又能达到边坡稳定性分析的要求。
3.降低精度:有限元强度折减法是一种新的技术,其计算结果与实际物理量有一定的偏差,而这个偏差一般比传统的有限元算法要小,因此使用有限元强度折减法计算边坡稳定性时,可以保证计算的可靠性。
三、结论有限元强度折减法在计算边坡稳定性方面具有优越的性能,具备计算精确、节省时间、降低精度等优点,因此作为计算边坡稳定性的一种有效工具已经得到广泛应用。
有限元强度折减法在岩土工程中的应用文章分析了有限元强度折减法在岩土工程中的应用优势,探析了有限元强度折减法在岩土工程中的应用,旨在为岩土工程相关工作人员提供一定的参考。
标签:有限元强度折减法岩土工程应用1前言英国科学家Zienkiewicz提出了用降低岩土强度或者增加外荷载,进行岩土工程安全系数的计算,该种方法被看做为岩土工程极限分析有限元法的雏形。
有限元强度折减法是极限分析有限元法最常用的方法之一,并且经过多年的实践应用获得了非常好的应用成果。
因此,文章针对有限元强度折减法在岩土工程中应用的研究具有非常重要的现实意义。
2有限元强度折减法在岩土工程中的应用优势分析(1)考虑周全。
有限元强度折减法在岩土工程中的应用,考虑问题比较周全,尤其是在细节方面,能够更加接近实际状况,在进行计算的过程中,计算数值越接近实际状况,其价值越高,能够有效的降低计算误差。
正是由于有限元折减法的这点优势,接近实际状况的计算方法,越来越被岩土工作者所青睐。
(2)适应性强。
我国幅员辽阔,涉及的地形包括高原、高山、山丘、平原等,一些岩土工程的测量工作会涉及到多种地形,如果采用传统的方法进行计算,不仅会增加工作人员的劳动强度,还会影响计算的精度,通过将有限元强度折减法应用在岩土工程测量中,通过分析相关地形,描述该地区地形的范围,以此评估相关地形的稳定性,无论施工现场的地形如何复杂,气候条件多恶劣,工作人员都可以在实验室中轻松的测量。
有限元强度折减法还有一个好处就是化整体为部分、化大为小,这样能够有效的降低工作人员的工作压力,降低计算过程中误差。
(3)准确性高。
国内岩土工程企业非常多,如果工程测量和计算的准确性较低,将会影响岩土工程企业的市场竞争力。
准确度是岩土工程测量中的必备条件,通过实践证明,有限元强度折减法在岩土工程测量中的应用,其结算结果比其他公式计算的结果更加准确,避免了其他公式就是你误差导致的误判,甚至导致地质变软以及工程坍塌等。
有限元强度折减法应用的几个问题及拓展有限元强度折减法是一种用于计算结构强度的方法,但在实际应用中存在一些问题和需要拓展的方向。
其中包括以下几个方面: 1. 结构材料的非线性特性。
有限元强度折减法通常基于线性弹性材料模型,但实际结构材料的行为往往是非线性的。
因此,在应用有限元强度折减法时需要考虑结构材料的非线性特性对结构强度的影响。
2. 结构的失效机制。
有限元强度折减法通常基于一定的失效准则,如最大应力准则和最大应变准则等。
但实际结构的失效机制可能与这些准则不完全符合,因此需要进一步研究结构的失效机制,以确定适合的失效准则。
3. 不确定性分析。
有限元强度折减法通常假设结构中的材料和几何参数是确定的,但实际结构中存在一定的不确定性,如材料性质、载荷大小和方向、制造误差等。
因此,需要进行不确定性分析,以评估这些因素对结构强度的影响。
4. 拓展到多物理场问题。
有限元强度折减法已经广泛应用于机械领域,但在其他领域如电子、生物等多物理场问题中的应用还需要进一步研究和拓展。
综上所述,有限元强度折减法在实际应用中需要考虑材料的非线性特性、结构的失效机制、不确定性分析和拓展到多物理场问题等方面的问题和拓展。
- 1 -。
有限元强度折减法和极限平衡法适用性研究曾红丽 戚明军陕西省建筑设计研究院有限责任公司 710003 西安摘要:简述了有限元强度折减法原理和两种极限平衡理论进行边坡稳定性分析的方法,分别利用ANSYS 和SLOPE/W 程序对一均质土坡进行了稳定性分析。
建立了同样尺寸的两种计算模型,得出了分别采用有限元强度折减法、Bishop 法和Janbu 法在未考虑地下水的情况下边坡安全系数。
将计算结果进行比较分析,表明基于极限平衡理论的分析方法分析得到的结果偏小,而有限元强度折减法由于考虑了土体内部应力应变关系所得结果更切合工程实际。
关键词:有限元强度折减法;极限平衡法;适用性;稳定性 1引言边坡工程是公路、铁路、水利水电以及矿山工程一个不可或缺的组成部分,因此边坡稳定性问题的研究就成为岩土工程界研究的热点问题。
国内外学者,已经取得了大量研究成果。
目前,研究边坡稳定性的传统方法主要有:极限平衡法,极限分析法,滑移线场法等。
随着计算机技术的发展,有限元强度折减法在边坡稳定性分析中已经得到了较好的应用[1~5]。
本文通过算例分析,比较分析有限元强度折减法和基于极限平衡理论的方法在边坡工程中适用性的差别。
2有限元强度折减法原理[6]所谓强度折减,就是在理想弹塑性有限元计算中将边坡岩土体抗剪强度参数逐渐降低直到其达到破坏状态为止,程序可以自动根据弹塑性计算结果得到破坏滑动面,同时得到边坡的强度储备安全系数K ,也称为强度折减系数。
通常,边坡的稳定性安全系数定义为沿滑动面的抗剪强度与滑动面的实际剪力的比值,公式表示为⎰⎰+=dA dAc K τϕσ)tan ( (1)将式(1)两边同时除以K ,得⎰⎰⎰⎰+=+=dA dA c dA dA KK c τϕστϕσ)'tan '()tan (1 (2)式中,K c c =',)arctan(tan 'K ϕϕ=,c 为粘聚力,ϕ为内摩擦角,σ为滑动面上的法向应力,τ为滑动面上的抗剪强度。
有限元强度折减法的原理、优点与超高边坡失稳的判据一、安全系数的定义两种方法可以导致边坡达到极限破坏状态,即:增量加载和折减强度。
传统边坡稳定分析中的安全系数是一个比值,假定一滑动面,根据力学的平衡来计算边坡安全系数,它等于滑动面以上土体条块的抗滑力与下滑力的比值。
式中K——安全系数;τ——滑动面上各点的实际强度。
将式子(4-1)两边同时除以k,上述公式变为其中:式(4-1)的左边等于I,表示滑坡体达到极限平衡状态,这意味着当代表强度的黏聚力和摩擦角被折减为1/K后,边坡最终到达破坏。
这个系数K就是有限元强度折减法中求解的安全系数,其实也就是强度折减系数。
二、有限元强度折减法的原理有限元强度折减法是在理想的弹塑性有限元计算中将边坡岩土体的抗剪强度参数:黏聚力c和内摩擦角φ按照安全系数的定义同时除以一个系数k,得到一组新的c′、φ′值,然后作为一组新的参数输入,再一次试算,如此循环。
当计算不收敛时,所对应的k被称为坡体的安全系数,此时边坡达到极限状态,将会发生剪切破坏,同时可以得到边坡的滑动面。
其中c′、φ′为三、有限元强度折减法的优点有限元强度折减分析法既具备了数值分析方法适应性广的优点,也具备了极限平衡法简单直观、实用性强的特点,目前被广大岩土工程师们广泛应用。
(1)不需要假定滑面的形状和位置,也无须进行条分。
只需要由程序自动计算出滑坡面与强度贮备安全系数。
(2)能够考虑“应力-应变”关系。
(3)具有数值分析法的各种优点,适应性强。
能够对各种岩土工程进行计算,不受工程的几何形状、边界条件等的约束。
(4)它考虑了土体的非线性弹塑性特点,并考虑了变形对应力的影响。
(5)能够考虑岩土体与支护结构的共同作用,并模拟施工过程和渐进破坏过程。
四、有限元强度折减法中超高边坡失稳的判据采用强度折减有限元方法分析超高边坡稳定性时,如何判断边坡是否达到极限平衡状态,十分关键。
这种有限元失稳判据的选取,没有获得共识,常见的失稳判据主要有下列三种。
有限元强度折减法在边坡稳定性分析中的应用有限元强度折减法在边坡稳定性分析中的应用边坡是山地地形中常见的地质体,其稳定性分析对于保护山体和预防地质灾害具有重要意义。
近年来,随着计算机技术的快速发展,有限元强度折减法逐渐成为边坡稳定性分析中的一种有效工具。
本文将介绍有限元强度折减法的原理和在边坡稳定性分析中的应用,并对其优缺点进行了探讨。
有限元强度折减法是一种基于强度准则的边坡分析方法。
它的基本原理是将岩土体的强度按照某个准则进行折减,并在有限元分析中采用折减后的强度参数进行计算。
其中,常用的强度准则有摩尔-库伦准则、维特曼准则等。
有限元强度折减法综合考虑了岩土体的强度特性和应力分布情况,相比传统的极限平衡法和全面滑动面法,能够更准确地评估边坡的稳定性。
有限元强度折减法在边坡稳定性分析中主要包括以下几个步骤:首先,确定边坡的几何形态和岩土体的材料性质。
其次,选择适当的强度准则,并对岩土体的强度参数进行合理的折减。
然后,根据边坡的几何特征和荷载情况,建立有限元模型,并进行边坡的数值计算。
最后,根据计算结果评估边坡的稳定性,并在必要时采取相应的加固措施。
有限元强度折减法在边坡稳定性分析中的应用具有以下优点:首先,有限元方法可以灵活地建立边坡的复杂几何模型,能够较准确地模拟实际工程中的复杂边坡形态。
而且,有限元强度折减法可以根据实际情况对岩土体的强度参数进行合理的折减,准确地反映岩土体的强度特性。
这些都为边坡稳定性分析提供了可靠的基础。
其次,有限元强度折减法可以综合考虑边坡的多种破坏机制,能够对于复杂的边坡情况进行全面的分析。
相比传统的极限平衡法和全面滑动面法,有限元强度折减法在考虑边坡的多种破坏机制时更加灵活准确。
此外,有限元强度折减法在计算过程中可以得到边坡内部的应力和位移分布情况,为边坡的设计和加固提供了参考依据。
通过对边坡各部分的应力和位移分布情况进行分析,可以找到边坡破坏的薄弱环节,并调整加固措施以提高边坡的稳定性。
有限元强度折减法在土坡与岩坡中的应用
有限元强度折减法当今备受瞩目,因它在土坡与岩坡中的应用受到广大人士的欢迎。
它拥有准确性、节省时间、方便使用等优点,但又并不会降低土坡与岩坡的安全性和可靠性。
有限元强度折减法由受力分析开始,根据力学模型计算出受力的工况。
进一步的,通过建立模型,预测出结构的屈服情况,以及各个构件的受力情况。
有限元强度折减法此时正式开始,强度折减前先对整个模型进行极限状态分析,根据极限状态分析结果,对每一个单元进行强度折减,折减强度及方式在此时确定,最后,通过受力分析计算,计算受力是否已足够说明安全性。
有限元强度折减法在实际应用中均取得较佳的效果,尤以土坡与岩坡的应用为典型。
由于储量的不足,常常会出现一些开采后的构件结构非但不坚固,而且结构不平衡,缺乏内在的可延续性,在此情况下,有限元强度折减法的应用大有裨益。
例如可以快速确定土坡与岩坡的最优构筑方式,减少岩体开采对施工效率的影响,有效改善构筑物稳定性,提高土坡与岩坡的安全性。
有限元强度折减法虽然有效,但也不可推而广之,还是要从实际情况出发,调整它的应用模式,做到结合工程实际,准确把握性和安全性,在开发业务中扮演重要角色。
总之,有限元强度折减法在土坡与岩坡的应用引起了广大科学家的共鸣,它的准确性、节省时间、方便使用以及提高安全性等特点,都得到了广泛认可。
有限元强度折减法1 背景1974年,Smith & Hobbs[1]使用有限元方法分析了φu=0条件下的边坡稳定性并与Taylar[2]的结果进行对比,得到了很好的一致性;1975年,Zienkiewicz等[3]考虑c’、φ’进行有限元边坡稳定性分析,其结果与圆弧滑面解有较好吻合;1980年Griffiths[4]验证了一系列具有不同材料特性和形状的边坡稳定性并通过与Bishop& Morgenstern[5]的结果进行了对比确定了数据的可靠性;此后也有研究证实了利用有限元方法进行边坡稳定性分析的可靠性[6,7,8,9];在文献[9]中,引入一些案例证明了有限元强度折减法的准确性,并证明了有限元强度折减法在分析非均质边坡时相对于传统方法的优越性。
2001年,郑颖人等[10]把有限元强度折减法引入国内,并对此进行了后续研究[11,12,13,14]。
相较于一些传统的边坡稳定型分析方法,有限元强度折减法有以下几个优点[9]:(1)不必假设滑面的位置和形状,当土体自身强度不足以抵抗剪应力时土体失稳会自然发生。
(2)由于有限元强度折减法中没有条分的概念,因此也不必假设条间力,在整体失稳之前土体都处于整体稳定状态。
(3)使用有限元方法能够查看破坏过程。
2 有限元强度系数折减法1.模型参数边坡模型主要包括六个参数,分别是:膨胀角ψ、内摩擦角φ’、黏聚力c’、弹性模量E’、泊松比υ’、重度γ。
膨胀角影响土体屈服后的体积变形,若ψ<0,则土体屈服后体积减小,若ψ>0则体积增大,ψ=0则体积不变。
ψ=φ的情况被称之为关联流动法则,但是此时ψ值通常高于实验观测值,特别是在侧限条件下会提高土的承载力预测值。
边坡稳定型问题通常是处于无侧限条件下,此时膨胀角的选取不再重要[9],因此文献[9]选取ψ=0条件下的非关联流动法则,并且通过案例分析可以得出此膨胀角的选取可以得出准确的安全系数以及滑动面。
有限元强度折减法1 背景1974年,Smith & Hobbs[1]使用有限元方法分析了φu=0条件下的边坡稳定性并与Taylar[2]的结果进行对比,得到了很好的一致性;1975年,Zienkiewicz等[3]考虑c’、φ’进行有限元边坡稳定性分析,其结果与圆弧滑面解有较好吻合;1980年Griffiths[4]验证了一系列具有不同材料特性和形状的边坡稳定性并通过与Bishop& Morgenstern[5]的结果进行了对比确定了数据的可靠性;此后也有研究证实了利用有限元方法进行边坡稳定性分析的可靠性[6,7,8,9];在文献[9]中,引入一些案例证明了有限元强度折减法的准确性,并证明了有限元强度折减法在分析非均质边坡时相对于传统方法的优越性。
2001年,郑颖人等[10]把有限元强度折减法引入国内,并对此进行了后续研究[11,12,13,14]。
相较于一些传统的边坡稳定型分析方法,有限元强度折减法有以下几个优点[9]:(1)不必假设滑面的位置和形状,当土体自身强度不足以抵抗剪应力时土体失稳会自然发生。
(2)由于有限元强度折减法中没有条分的概念,因此也不必假设条间力,在整体失稳之前土体都处于整体稳定状态。
(3)使用有限元方法能够查看破坏过程。
2 有限元强度系数折减法1.模型参数边坡模型主要包括六个参数,分别是:膨胀角ψ、内摩擦角φ’、黏聚力c’、弹性模量E’、泊松比υ’、重度γ。
膨胀角影响土体屈服后的体积变形,若ψ<0,则土体屈服后体积减小,若ψ>0则体积增大,ψ=0则体积不变。
ψ=φ的情况被称之为关联流动法则,但是此时ψ值通常高于实验观测值,特别是在侧限条件下会提高土的承载力预测值。
边坡稳定型问题通常是处于无侧限条件下,此时膨胀角的选取不再重要[9],因此文献[9]选取ψ=0条件下的非关联流动法则,并且通过案例分析可以得出此膨胀角的选取可以得出准确的安全系数以及滑动面。
c’和φ’指Mohr-Coulomb准则中边坡土体的有效黏聚力和内摩擦角;E’和υ’是土体材料的弹性参数,这两个参数对土体稳定性分析的影响较小;γ是土体的重度。
应用有限元方法进行边坡稳定性分析中最重要的三个参数是c’、φ’、和γ。
2.屈服条件(1)Mohr-Coulomb准则Mohr-Coulomb准则用大小主应力表示如式(1)所示:σ1′−σ3′2=σ1′+σ3′2sinφ′−c cosφ′(1)其中, σ1′、σ3′分别指土中一点的大小主应力。
在主应力空间中,如果不考虑σ1、σ2、σ3之间的大小关系,屈服面是一个不等角六棱锥,在π平面上是一个等边不等角六边形。
(2)D-P准则(4) (6) D-P 准则可以写成式(2)形式:−βI 1+√J 2=k f (2)其中I 1为第一应力不变量、J 2为第二偏应力不变量,β和k f 为试验常数。
在主应力空间中其屈服面为一个圆锥,在π平面上是一个圆形。
3)D-P 准则转换为Mohr-Coulomb 准则首先引入参数b ,如式(3)所示: b =σ2−σ3σ1−σ3(3) 则,I 1和√J 2分别可转化为式(4):I 1=3(σ1+σ3)2+(b −12)(σ1−σ3)√J 2=λ√31−σ3) 其中λ=√1−b +b 2将其带入(2),得式(5):σ1−σ3= 1.5β2−bβ+λ√3(σ1+σ3)+k f β2−bβ+λ√3 (5)与式(1)对比可知两个准则之间的转换关系如式(6)所示:sinφb = 1.5β2−bβ+λ√3 2c b cosφb =k f β2−bβ+λ√3 因此,当b=0时,即外角点外接DP 圆的两个试验常数分别如式(7)所示,当b=1时,即内角点外接DP 圆的两个试验常数分别如式(8)所示。
β=√3(3−sinφ)k f =√3(3−sinφ) (7) β=√3(3+sinφ)k f =√3(3+sinφ) (8)21σ2≥≥σ2σ1≥b P 圆内角点外接D3.安全系数的定义(1)Mohr-Coulomb 准则中的安全系数1955年,Bishop [15]首先在边坡稳定性分析中提出了抗剪强度折减的概念,在有限元强度折减法中通过将坡体的强度参数:黏聚力c 和内摩擦角φ同时除一个折减系数F t ,得到一组新的c’和φ’值,作为一个新的强度参数输入进行试算,(10) 当计算不收敛时,对应的F t 即为所求的安全系数,此时坡体达到极限状态,发生剪切破坏。
c’=c/F tφ’=arctan(tanφ/F t )(2)D-P(Drucker-Prager)准则中的安全系数取F t 为D-P 准则中的强度折减系数,则D-P 准则可以表示为式(9),−βF t I 1+√J 2=k f F t (9) (3)不同屈服条件下安全系数转换[13]首先引入Mohr-Coulomb 等面积圆屈服准则,在π平面上,其屈服面是一个圆,并且面积与Mohr-Coulomb 准则的不等角六边形相等,Mohr-Coulomb 等面积圆屈服准则中的试验参数如式(10)所示: −β=k f 式中θδ=arcsin −23Asinφ+[49A 2sin 2φ−4(sin 2φ3+1)(A 23−1)]122(sin 2φ3−1),A =√26√3简称外接圆屈服准则为DP1准则,其试验常数分别为β1,k f1;Mohr-Coulomb 等面积圆屈服准则为DP2准则,其试验常数分别为β2,k f2。
把DP1准则表示为f 1=√J 2=β1I 1+k f1,DP2准则可表示为f 2=√J 2=β2I 1+k f2。
令η=β1\β2=k f1\k f2=f(φ),f 1=β1I 1+k f1=ηβ2I 1+ηk f2,所以f 1f 2=ηβ2I 1+ηk f2β2I 1+k f2=η=f(φ)。
由此可知,η是φ的函数,当φ取不同值时可以得到不同的η值如表1所列:表 1 不同内摩擦角时的η值4.失稳判据目前两个比较主流的失稳判据分别是有限元计算中力不平衡和位移的不收敛以及广义塑性应变或者等效塑性应变从坡脚到坡顶贯通。
Griffiths [9]和郑颖人[11,12,13,14]都使用计算不收敛作为失稳判据。
Griffiths [9]提出,当在用户定义的最大迭代数目下计算仍不收敛时,则没有任何一种应力分布方式可以同时满足Mohr-Coulomb 准则以及整体稳定,这种情况可看做边坡失稳判据。
边坡失稳与数值计算不收敛同时发生,并伴随着极大的节点位移,并以1000作为最大的迭代步数。
郑颖人[14]提出,有限元的计算迭代过程就是寻找外力和内力达到平衡状态的过程,整个迭代过程直到一个合适的收敛标准得到满足才停止。
可见,如果边坡失稳破坏,滑面上将产生没有限制的塑性变形,有限元程序无法从有限元方程组中找到一个既能满足静力平衡又能满足应力-应变关系和强度准则的解,此时不管是从力的收敛标准,还是从位移的收敛标准来判断有限元计算都不收敛。
3 案例分析例一,不含地基的均质边坡[9]该边坡如图1所示,有限元程序采用Mohr-Coulomb失效准则,建立平面应变条件下八节点四边形单元减缩积分计算模型,其强度参数为φ’=20°,c’/γh=0.05。
边坡坡度为26.57°(2:1),坡底水平,其边界条件为坡底约束竖直方向位移与水平方向位移,左侧约束水平方向位移,其余面为自由面。
施加重力荷载后使安全系数从0.8到1.4逐步变化直至计算不收敛图 1 不含地基的均质边坡每一个安全系数对应的迭代次数如表2所列,当真正的安全系数接近时需要更多的迭代次数。
表 2 例一计算结果当安全系数为1.4时,无量纲位移E’δmax/γH2突变,并且此时计算无法收敛,在此情况下有限元计算结果与Bishop & Morgenstern[5]给出的结果吻合良好,如图2所示。
图 2 安全系数与无量纲位移边坡失稳时(FOS=1.4)节点位移矢量和网格变形如图3(a)和图3(b)所示,由此可得到边坡的潜在滑动面。
图3安全系数为1.4计算不收敛时边坡变形(a)节点位移矢量(b)网格变形例二,有软弱层的不排水黏性土边坡在本案例中,使用Tresca准则(φu=0)进行总应力分析。
边坡几何形状如图4所示,地基厚度与边坡高度相同,该边坡有一个软弱层,在有限元计算中,令其抗剪强度(C u2)在一定范围内变化但其周围土体抗剪强度保持C u1/γH=0.25不变。
利用有限元方法计算该边坡的安全系数结果如图5所示,对于均质边坡情况,C u2/C u1=1,有限元计算结果与Taylor[2]的结论很接近,随着软弱层的强度逐渐减小,在C u2/C u1≈0.6时,结果发生了明显的变化。
分别假定圆弧滑面和穿过软弱面的三段线滑面并利用Janbu法计算安全系数,可见在C u2/C u1≈0.6处也发生了滑动机制的转换,当C u2/C u1>0.6时,潜在滑面形状为圆弧,当C u2/C u1<0.6时,潜在滑面为结构软弱面。
图6更加清晰的展示了这一现象,图6(a)为均质边坡(C u2/C u1=1)时的潜在滑面,可见此时的滑面形状为圆弧滑面,与Taylor[2]的预测相同;图6(c)为软弱层强度只有其周围土体20%( C u2/C u1=0.2)时的潜在滑面,此时潜在滑面沿软弱层发展;图6(b)为软弱层强度只有其周围土体60%( C u2/C u1=0.6)时的潜在滑面,此时圆弧滑面和沿软弱层的三段线式滑面都有可能发展,至少存在两种明显的滑动机制。
图4有软弱层的不排水黏性边坡图 5 不同软弱层强度时的安全系数图 6 不同软弱层强度下的网格变形(a) C u2/C u1=1.0 (b) C u2/C u1=0.6 (c) C u2/C u1=0.2例三,不同坡度边坡安全系数计算[13],验证Mohr-Coulomb等面积圆屈服准则均质边坡,坡高H=20m,土容重γ=25kN/m3,黏聚力c=42kPa,内摩擦角φ=17°,求坡角β分别为30°,35°,40°,45°,50°时边坡的安全系数。
计算结果如表3所列。
表 3 安全系数计算结果从表中计算结果可以看出,采用外接圆屈服准则计算的安全系数比传统的方法大许多,采用莫尔-库仑等面积圆屈服准则计算的结果与传统极限平衡方法(Spencer 法)计算的结果十分接近,说明采用莫尔-库仑等面积圆屈服准则来代替莫尔-库仑不等角六边形屈服准则是可行的,这样使计算大为方便。
而采用外接圆屈服准则计算的安全系数要比莫尔-库仑等面积圆屈服准则计算的结果大η(1.21)倍。
例四,存在两组节理面的岩质边坡稳定性分析[12]如图7所示,岩体中存在两组方向不同的软弱结构面,贯通率100%,第一组软弱结构面倾角为30°,平均间距10m ;第二组软弱结构面倾角75°,平均间距10m 。