剪胀角对分层边坡稳定性的影响.docVIP

PAGE PAGE 1 剪胀角对分层边坡稳定性的影响 　　摘要：利用有限元软件ABAQUS，采用强度折减法，对不同层数、不同材料、不同高度的分层边坡剪胀角进行折减，分析了剪胀角对分层边坡稳定性的影响。结果表明，对分层边坡，任意土层的剪胀角进行折减时，安全系数Fs随着剪胀角的增加而增大，增大的幅度逐渐减小，其变化规律不受土层材料、高度、层数的影响。与边坡整体的剪胀性相比，单层土的剪胀性对Fs的影响要小，且边坡上层土层的剪胀性对边坡稳定性的影响较大。随着剪胀角的增大，边坡位移趋势越加明显，且塑性区的范围逐渐扩大，滑弧半径增大，安全系数随之增大。在实际工程中，宜对分层边坡整体的剪胀角进行折减并慎重取值。 　　关键词：剪胀系数；分层土；软弱土层；强度折减；安全系数 　　中图分类号：O319文献标识码：A文章编号2013　　岩土在剪切过程中会产生体积变形，正常固结的黏土和低密度的砂土一般会产生剪缩，而超固结土及密实度较高的砂土则产生剪胀[1]。由于对岩土剪胀性的认识还不充分，大多数岩土工程分析中选择了两种极端情况：关联流动法则（剪胀角ψ=φ）和非关联流动法则（剪胀角ψ=0）。关联流动法则意味着土体在剪切的过程中可以产生无限制的体积膨胀现象[2]，非关联流动法则ψ=0则认为土体是不可膨胀的，由此获取极限荷载。 　　关于剪胀角对岩土工程稳定性的影响问题，不少学者进行过探讨：康亚明[3]等通过剪胀角对边坡潜在滑动面影响的研究，得出的结论认为，随着剪胀角的减小，边坡的塑性区从坡顶向坡脚有贯通的趋势；蒋青青等[4]研究了坡顶超载情况下的边坡剪胀效应；吴春秋[5]提出边界荷载会使剪胀性更显著。至于工程实践中剪胀角的取值问题，赵星光[7]研究了岩石峰值内摩擦角和剪胀角的关系，认为岩石在零围压时的峰值剪胀角小于并近似等于峰值内摩擦角；张培文[7]提出剪胀角对地基极限承载力的影响是不可忽略的，并认为只有考虑了材料的剪胀性、正确选用材料的屈服准则，才能得到符合理论与实际的边坡稳定安全系数[8]；孔位学等[9]通过对速度矢量和破坏面的夹角分析，得出剪胀角取0.5倍的摩擦角的结论，但却不一定适应具有软弱层的分层边坡。 　　因此，本文在前人的研究成果基础上，利用有限元软件ABAQUS分析剪胀角对分层边坡以及有软弱层的边坡稳定性的影响，以期为工程实践提供指导。 　　1有限元强度折减法失稳破坏判据 　　目前，判断边坡失稳破坏的标准通常包括：有限元数值计算的不收敛、塑性区的贯通、广义剪应变的贯通等[10]。吕庆[11]认为当采用小应变模式假设时，应用数值计算不收敛作为判据。然而，计算不收敛的原因比较多，如荷载过大或者计算单元奇异等，因此以不收敛为判据有一定的局限性。栾茂田[12]建议采用塑性应变贯通作为判据。不过尽管塑性应变在一定程度上可以描述塑性区的发展过程，但作为评判依据时因为受主观因素影响大，并且未排除弹性塑性应变的影响，所以破坏界限比较模糊。 　　分析边坡失稳破坏的主要特征，可以发现不管岩土内部变形机理是广义的剪应变还是塑性应变，最终结果都是边坡产生位移，所以位移是边坡内部作用的外在表现，而滑动主要是由剪应变和位移造成的。有研究发现[13]，随着强度参数的不断折减，边坡上的位移矢量和剪应变不断向坡脚增大，因此以典型的坡顶特征点位移突变为失稳判据，意义明确、界限清晰。由此，笔者选择了边坡特征点的位移矢量与强度折减系数关系曲线中的突变点作为边坡的破坏判据。因为在该点，边坡处于极限平衡状态，塑性区刚刚贯通但不是非常明显，若以塑性区贯通为判据，在感官上判别并不准确。该点过后特征点的位移矢量会突然增大，意味着进入了破坏阶段。 　　2MohrCoulomb模型的剪胀性及屈服准则 　　2.1MohrCoulomb模型的屈服面 　　ABAQUS提供了许多岩土体材料的本构模型，大多数的模型参数都可以与温度场变量结合，从而实现强度折减。MohrCoulomb模型在单调加载条件下的适用性、收敛性较好，本文拟采用该模型，其屈服面函数为： 　　F（p，q，Θ）=Rmcq-ptanφ-c=0（1） 　　其中：c为粘聚力；φ为摩擦角；q为广义剪应力；p为球应力；Θ是极偏角；Rmc控制了屈服面在π平面的形状，按下式计算： 　　Rmc=113cosφsin（Θ+π13）+113cos（Θ+π13）tanφ（2） 　　其中，cos（3Θ）=r31q3，r是第三偏应力不变量J3；Θ是极偏角。 　　2.2MohrCoulomb模型的剪胀性 　　MohrCoulomb屈服面存在尖角，如采用相关流动法则，即塑性势面与屈服面相同，将会在尖角处出现塑性流动方向不是唯一的现象，为了避免这一问题，采用如下形式的连续光滑的椭圆函数作为塑性势面： 　　G=（εc|0t