本文首先扼要介绍了土石坝,特别是心墙坝的发展历程和目前存在的主要问题,对目前土石坝有限元应力变形计算分析常用的三种本构模型—邓肯E-υ模型、修正剑桥模型和南水模型及其特点作了简要介绍和分析。然后,利用上述三种模型,对已取得较为完整原型观测资料的新疆“635”水利枢纽砂砾石粘土心墙坝和云南鲁布革堆石粘土心墙坝的应力变形性状进行了较为详细的有限元计算分析。有限元计算分析结果与原型观测资料的比较表明: 1.三种模型能较好的模拟心墙坝施工期的应力变形性状,只是邓肯模型计算得出的沉降和水平位移与实测值相比偏大,修正剑桥模型计算得出的沉降和水平位移偏小,南水模型的计算结果介于上述两个模型之间,但较之于实测值仍偏小。 2.三种模型在模拟心墙坝蓄水期应力变形性状时存在严重不足,特别是心墙和上游坝壳变形的计算分析结果与实测值定性不一致,即计算得出心墙和上游坝壳在蓄水期明显上抬,邓肯模型上抬量最大,修正剑桥模型次之,南水模型最小。而原型观测资料表明,蓄水期心墙和上游坝壳则明显下沉。究其原因,一是计算时未考虑坝壳料的湿化变形;二是模型的加卸荷准则存在不足。邓肯模型属非线性弹性模型,修正剑桥模型为单屈服面弹塑性,对于蓄水期上游坝壳及心墙中广泛区域内围压或平均应力减小的情况,上述两模型均判定为卸荷,从而不产生塑性变形。而南水模型为双屈服面模型,在较大程度上克服了上述不足,因此计算得出的心墙和上游坝壳的上抬量明显小于邓