摘要：广义塑性梯度模型在双屈服函数中分别考虑塑性剪切应变和塑性体积应变的拉普拉斯项,用虚功方程和屈服函数的一致性条件建立控制微分方程。它有望在较全面地反映岩土介质基本力学性质的同时,也合理解释岩土介质的应变局部化现象。为得到上述控制方程的数值解答,同时对位移和塑性乘子进行离散,得到以节点位移和节点塑性乘子为基本未知量的非线性方程组,建立塑性乘子的边界条件,给出求解上述方程组的增量算法。数值算例表明：（1）发生应变局部化时,塑性剪切应变和塑性体积应变主要集中在应变局部化带中;（2）网格密度变化时,网格敏感性问题不再出现;（3）当模型的局部化参数增大时,应变局部化带的带宽也明显增宽。

关键词：数值分析 软化 岩土介质 应变局部化 广义塑性梯度 

单位：长沙理工大学桥梁与结构工程学院; 湖南长沙410076; 西安理工大学岩土工程研究所; 陕西西安710048; 后勤工程学院; 重庆400041