还剩4页未读,继续阅读
文本内容:
对非线性破坏准则下边坡稳定性分析的线性简化马崇武 武生智 苗天德摘 要基于极限分析的上限理论对非线性破坏准则的边坡稳定性问题提出了一个线性简化方法.该方法把复杂的非线性问题化成线性问题从而在实际工程应用时可以利用已有的关于线性问题的分析结果或根据本文论述的方法使问题简化.数值算例表明了该方法的简单实用性.关键词边坡稳定;非线性破坏准则;极限分析中图分类号P
642.22 文献标识码A 文章编号0455-2059199901-0049-04LinearizationonSlopeStability____ysiswithNonlinearFailureCriterion__Chongwu,WuShengzhi,MiaoTiandeDepartmentofMechanicsLanzhouUniversityLanzhou730000ChinaAbstract Basedontheupperboundmethodoflimit____ysisamethodissuggestedwhichconvertstheslopestabilityproblemwithanonlinearfailurecriteriontoonewithalinearcriterion.Comparisonofthenumericalresultsbythelinearmethodsuggestedherewiththoseobtainedbyusinganonlinearfailurecriterion__ailableintheliteratureshowsthatthis______methodissatisfactoryforengineeringpracti__.Keywords slopestability;nonlinearfailurecriterion;limit____ysis 在边坡稳定性分析中极限平衡方法已得到广泛应用但是一般无法确保它的解答是精确解答的上限还是下限.极限分析的上限分析方法所得的解答确为精确解答的上限,因此,该方法越来越得到广泛的应用.对于线性破坏准则的边坡稳定性问题,已有许多人利用极限分析理论进行了研究和探讨[1~3]. 在许多实际工程问题中,已有充分的实验数据表明,其破坏准则具有极高的非线性,且往往因表征材料特性的资料难以得到而无法直接用有限元等方法解决.由于在分析边坡的稳定性问题时,我们关心的是边坡的稳定性,并非详尽的应力应变历史,为此,Zhang等[4]利用极限分析的上限理论对非线性屈服条件的边坡稳定性问题进行了详细分析. 本文利用极限分析理论,对非线性屈服条件的边坡稳定性问题提出了一个把复杂的非线性问题简化为线性问题的方法,以下简称线性简化法.附图 屈服准则和流动法则Fig.Ageneralyieldandflowrule1 基本假定
(1)土体发生破坏的瞬间,几何形状的改变较小,虚功原理可以应用;
(2)土体是刚塑性的,在破坏面上满足形式为τ=fσ1的非线性破坏准则,其中τ,σ分别为破坏面上的剪应力和正压力;
(3)土体材料满足相关联的流动法则. 若土体满足上述假定,在(σ,τ)空间,如附图所示,定义φt为切线摩擦角,22 极限分析的上限方法 按照极限分析的上限理论,对于任何假定的破坏机制,如果边坡体在外力作用下所作的功率等于耗散的内功率,便可得边坡发生破坏的极限荷载(或安全因子)的上限.该方法在边坡稳定性分析中已广泛应用[1,4].对于二维问题,表达式为3该等式简称为能量平衡方程.对于刚塑性体,式
(3)可简化为4 对仅受自重作用的均质边坡本文仅分析旋转式破坏机制.若记外力即重力所作的总功率为,耗散的内功率为D,则有56这里γ为容重,vγ为区域A中点沿γ方向的速度分量,区域A由边坡边界S与速度间断面线L围成,[v]为速度间断标量. 若材料满足线性Mohr-Coulomb屈服条件τ=c+σtgφ,7则方程6简化为D=∫Lccosφ[v]dL.8这里cφ分别为土的粘聚力和内摩擦角.3 极限荷载
3.1 Zhang-Chen法 Zhang和Chen依照极限分析的上限理论,对非线性破坏准则的边坡稳定性问题,用变分原理得到了极限荷载最小上限必须的条件——两个微分方程和横截条件.为简化计算,他们把复杂的微分方程组转化为较简单的初值问题,并把该数值过程称为“逆方法”,详见文[4].
3.2 线性简化法 极限分析理论表明[5]“在结构的任何部分提高材料的屈服极限,不会降低结构的极限载荷.”因此,若给定一屈服条件超过真实屈服条件,则对应该屈服条件的任一运动学解均为极限载荷对应真实屈服条件的上限. 考察附图中由直线l确定的线性屈服条件,l为经过非线性屈服面线上点B并与其相切的直线,若非线性屈服面线的斜率随应力σ单调递减,则直线l完全落在该屈服面线的外侧.于是,对应线性屈服条件由l确定的极限荷载不比对应非线性屈服条件的极限荷载小.极限荷载的最小上限值由相应线性屈服条件的极限荷载该值随B的位置变化而改变的最小值确定.直线l确定的线性屈服条件为τ=ct+σtgφt.9若非线性屈服条件取1式,即τ=fσ则φt=arctgf′σB,10ct=fσB-σBtgφt=fσB-σBf′σB.11其中σB为切点的σ值. 以上分析表明,对非线性破坏准则的边坡稳定性问题可根据相对应的线性破坏准则估算极限荷载的最小上限.不过在计算时,须用φtct代替8式的φc.4 举 例 Zhang等在文[4]中对具有12非线性屈服条件的边坡稳定性问题进行了详细分析,他们定义稳定因子13这里h1为边坡的允许高度.为方便计,我们记线性简化法中对应非线性屈服条件的边坡稳定因子为Nn,对应线性屈服条件的边坡稳定因子为Nl,则有14联立式9~12和1415其中 关于Nl在这里采用Chen[1]的结果,由于Chen的结果是关于φ每隔5°给出的,我们可用抛物线插值法将Nl表示成φt的函数.求极限荷载的最小上限成为求Nn关于φt的最小值.线性简化法与Zhang-Chen法的数值结果比较见附表.结果表明,线性简化法与Zhang-Chen法的结果相差很小,前者比后者仅高不到
1.5%.附表 稳定因子C=90kPaσt=
247.3kPaTable ThestabilityfactorsC=90kPaσt=
247.3kPa系数m坡角β/°
907560451.
05.
5107.
48010.
39016.
1805.
5107.
48010.
39016.
1801.
25.
1306.
7708.
95012.
5505.
1576.
7958.
99812.
6021.
44.__
06.
3308.
13010.
8204.
9256.
3668.
17710.
9081.
64.
7306.
0407.
6109.
7004.
7636.
0737.
6579.
8531.
84.
6005.
8207.
2409.
1004.
6445.
8607.
2879.
1902.
04.
5205.
6606.
9708.
7804.
5505.
7017.
0178.690*
2.
54.
3505.
4006.
5407.
9504.___
5.
4336.
5937.918*注上行是Zhang-Chen的结果;下行是线性简化法的结果.5 结果与讨论 a 本文提出的线性简化方法,对于仅受自重作用的均质边坡,采用与文[4]相同的屈服条件,其稳定因子仅比文[4]的结果高估了不到
1.5%,但计算过程要简单得多. b 附表中有两个注的数值(对应低坡角,高非线性指标)与理论分析不符,它们应该比Zhang-Chen的计算值大,现在反而小,估计差别是由Zhang-Chen的计算方法使其值偏大所引起. c 对于非线性屈服条件的稳定性问题,可以利用已有的关于线性屈服条件的结果或根据本文提出的方法化为线性问题,从而使问题得以简化. d 只要能找到超出相应非线性屈服条件的线性屈服条件.线性简化法对诸如考虑孔隙水压力作用的稳定性问题以及地基承载力的计算等均适用.__简介马崇武1965-,男,__生.__单位兰州大学力学系, 甘肃兰州 730000____[1] ChenWF.Limit____ysisandsoilplasticity[M].Amsterdam:Elsevier,1975.[2] MichalowskiRL.Slopestability____ysis:akine__ticalapproach[J].Geotechnique1995452:283~293.[3] DrescherADetournayE.Limitloadintranslationalfailuremechani__sforassociativeandnon-associative__terials[J].Geotechnique1993433:443~456.[4] ZhangXJ,ChenWF.Stability____ysisofslopeswithgeneralnonlinearfailurecriterion[J].IntJNumer____MethodsGeomech1987111:33~50.[5] 余同希.塑性力学[M].北京:高等教育出版社1991.收稿日期1998-04-06.。