边坡及挡土墙变形局部化分析.pdfVIP

水 利 学 报 2004 年12 月 SHUILI XUEBAO 第12 期 文章编号: 0559-9350(2004)12-0021-06 边坡及挡土墙变形局部化分析 1 2 2 1 冯吉利，孙东亚，丁留谦，隋允康 (1.北京工业大学，北京 100022；2.中国水利水电科学研究院 防洪减灾研究所，北京 100044) 摘要：本文将局部和非局部软化损伤线性结合在一起表示非局部软化损伤内变量，给出变形体边界上积分点非局 部内变量的复合映射方式。用非局部连续莫尔-库仑模型和次塑性模型(Hypoplasticity)，分别对边坡和挡土墙的 变形局部化破坏过程进行了模拟分析。二维、三维有限元计算结果表明，本文采用的软化损伤内变量非局部方式 和非局部边界的处理是行之有效的，反映材料微观性态的特征尺度以统计平均的方式表达，在材料本构关系中加 以体现，所得计算结果对有限元网格的划分具有不敏感性，因而用非局部连续模型模拟土工建筑物的变形破坏过 程具有可行性。 关键词：变形局部化；软化损伤；边坡稳定；次塑性模型；有限元 中图分类号：O347 文献标识码：A 土工建筑物如边坡、挡土墙后的土体在外载作用下，常沿某一部位形成剪切变形带。由于传统弹塑性 理论未涉及材料的微观尺度和微观缺陷，因而不能用其合理地描述剪切带现象。近几十年来，人们不断尝 试使用各种手段、寻找各种可能的途径，以揭示变形局部化问题的物理实质，使工程设计更加安全合理[1]， 如二阶梯度弹塑性模型、增强应变强不连续弹塑性模型，基于统计物理的非局部连续模型等。二阶梯度模 型要求在变形梯度表达式中引进变形梯度二次项，即在应力应变关系中自然包含变形的高阶效应，所以其 本构表达式是高度非线性的[2]。增强应变强不连续模型目前仅限于二维问题，向三维问题扩展尚有许多困 难，比如剪切带空间面发展方向如何确定，如何合理选用跳跃间断函数[3]。诸如岩石、土等材料，在变形 过程中其局部缺陷和微观裂纹一般跨越多种尺度，而目前适用于金属材料的细观尺度晶格原子模型[4]未经 改进，尚无法满足岩土工程分析计算的要求。非局部模型从统计物理的观点来考虑材料的微观缺陷，以积 分形式考虑材料特征尺度对材料本构关系的影响，从而揭示了材料相互作用的集体行为。本文应用非局部 连续模型，尝试用局部化和非局部化软化损伤内变量的线性组合来表示非局部软化量，其目的是增强计算 迭代过程的稳定性。变形体边界上的软化损伤变量采用映射复制方式形成非局部量，探讨用非连续局部模 型来揭示边坡和挡土墙变形破坏特性，二维、三维有限元计算结果显示非连续局部模型在解决岩土工程变 形破坏问题中是有效的。 1 非局部连续模型的基本形式及其边界条件 一般而言，非局部场变量可引入任意材料的本构模型中。在具有各向同性软化损伤内变量的标准塑性 流动理论中[5]，屈服函数为以下形式 f( σ, κ)=F( σ)- σ(κ) (1) Y 式中：σ为应力；κ为软化损伤量；F( σ)为等效应力；σY为当前屈服应力。 收稿日期：2003-10-14 作者简介：冯吉利(1963-)，男，博士，副教授，主要研究方向为应用力学、结构材料及岩土工程。 1 水 利 学 报 2004 年12 月 SHUILI XUEBAO 第12 期 屈服应力的演化由下列软化损伤定律表示 σ(κ)= σ+h( κ)