本构理论是岩土工程中的关键部分,至今仍存在一些问题：基于虎克定律的弹性及非线性弹性模型不能反映土的剪胀这一重要特性；非关联模型不符合Drucker塑性公设；塑性应变增量方向不具有唯一性；传统弹塑性模型本构矩阵非正定性等,因而传统理论不足以表述岩土的本构特性。杨光华根据岩土的变形特性提出了广义弹塑性本构理论的思想,传统的弹性位势理论和塑性位势理论都可归结为数学上的势函数问题,因而可以用广义位势理论来表达和统一,这给岩土材料的本构关系研究指出了更开阔的数学背景。从广义位势理论出发来建立土的本构模型突破了传统理论中以塑性公设为前提的限制,可以无需推求塑性势函数和屈服函数,且其所能表述的本构关系比通常的塑性位势理论更为广泛,因而具有更广阔的前景。本文主要围绕广义位势理论进行进一步深入的研究。剖析了各模型族本构矩阵的特性,通过算例展示了广义位势理论多重势面模型简便优越的特点；在广义位势理论的基础上,研究了塑性应变增量的分解准则,提出了考虑拟弹性塑性变形的土体弹塑性本构模型,把传统不可恢复的塑性应变增量分解为具有弹性应变特性的拟弹性塑性变形部分和符合传统塑性理论的纯塑性部分；提出和建立了带拟弹性参数的数值弹塑性模型,并进行了相关验证；在多重势面模型基础上研究了碎石桩的机理问题。本文主要完成的具体工作如下：首先,对土的本构理论相关文献资料进行全面回顾,针对国内外学者研究岩土本构模型的方法进行了简单归纳和分类概述,分析了岩土材料本构理论的发展和一些具体模型的特点,阐述了岩土非线性弹性本构理论和传统弹塑性本构理论以及最新的广义位势理的研究成果及其发展现状,并对各种理论的优缺点进行了比较,分析它们之间的联系与区别,找出存在的问题和有待进一步深入的方向,指出统一的数学基础和具体土类在特定假设下的本构模型是以后发展趋向。第二,系统阐述了广义位势理论的提出和发展,从理论基础的角度讨论了广义位势理论的内在特性及重要意义。针对土的变形的特殊性,主要介绍了应力空间多重势面模型和应变空间多重势面模型,对模型特性进行了全面展示,分析了其对于邓肯-张模型和传统弹塑性模型理论上的优越性。由于不再用屈服函数,即不再区分明显的弹、塑性阶段,这也更符合工程土的实际变形特性,有更好的适应性。通过试验验证表明多重势面模型参数简单,计算可靠,计算结果与试验符合较好,能充分利用试验成果,又能反映邓肯-张模型反映不了的剪胀特性,证明其不论理论基础还是工程实用性上都比较优越。第三,进行了弹性矩阵与弹塑性矩阵的对比剖析,从数学本质上探讨了非线性弹性模型和传统弹塑性模型以及多重势面模型的的数学基础和假设及其联系；阐明了其主要是刚度(柔度)矩阵中角占优以及正定的问题,指出邓肯模型和传统弹塑性模型以及多重势面模型的本构矩阵内在差异为：有没有反映球应力对剪应变和偏应力对体应变交叉影响的那些元素,以及这些元素之间不同的等量关系,进而建议了对几种主要土的模型使用原则。第四,指出了邓肯-张模型和传统塑性理论关于应力计算存在的问题,即本构矩阵的不同将在反算应力时导致较大差异,邓肯模型使应力计算结果偏大,传统关联流动弹塑性模型因为假设交叉项相等也会导致与实际不符,只有满秩且上下三角元素独立的多重势面模型才能得出更符合实际的应力计算结果。通过算例比较了用不同本构模型进行应力计算所得结果的差异,分析了多重势面模型的合理性。第五,进一步深入分析了Ottowa砂土的变形规律。为解决传统理论对塑性应变增量方向描述的不足,本文在广义位势理论的基础上,研究了塑性应变增量的分解准则,提出了考虑拟弹性塑性变形的土体弹塑性本构模型,把传统不可恢复的塑性应变增量分解为具有弹性应变特性的拟弹性塑性变形部分和纯塑性部分,这在理论上是合理的。拟弹性部分遵从弹性法则,并与应力增量有相同的方向,采用弹性模型表示；纯塑性部分遵从传统塑性理论,方向具有唯一性,采用关联流动法则进行建模。这相当于在塑性矩阵中拿出拟弹性这一部分合并到弹性矩阵部分,从理论上阐明了这一措施可改善总的弹塑性矩阵性质,也更符合土的实际变形机理。这样分解后建立的模型将更为合理和简便,又可以解决土塑性应变增量方向不具有唯一性的问题。第六,推导了拟弹性分解方法并给出相应参数的确定准则；在广义位势理论的框架下再增加两个拟弹性参数,从而在计算中一开始就包含两项塑性应变增量：拟弹性和纯塑性,两项各自随着应力水平的变化而变化；当拟弹性系数为零时即退化为传统的关联流动模型,因而涵盖更广。通过新模型对水坠坝冲填土的计算证明了进行分解的可行性和优越性。在此基础上研究了关联和非关联流动法则产生的差异,评估了塑性应变中符合关联法则的部分和非关联法则的部分的比例；给出关联与非关联流动法则的使用建议。第七,建立了数值弹塑性模型。首先通过直接从试验曲线求导得到多重势面本构模型增量计算所需参数,对得到的塑性系数按所提出的分解准则进行拟弹性分解；然后把多组围压下的常三轴试验结果综合来建立数据库,通过立方插值得到其他应力状态下单元的模型参数,即得到数值弹塑性模型。包含三种模型：做关联流动假设的两参数简化模型、不做关联流动假设并分解出拟弹性应变的六参数模型、直接将试验与简化模型结果差值部分归为拟弹性的四参数模型。因为真实应力路径是由σ1、σ2、σ3控制的,不同的p值也反映了三个主应力的不同变化路径,数值模型充分利用试验成果所蕴含的力学信息,可在一定程度上反映应力路径的影响；数值模型参数纳入力学基本方程,有着明确的力学概念；处理后分别用弹性准则和关联流动法则进行建模,计算稳定可靠,且能更简便、更真实地反映土的变形特性。第八,用编入了多重势面模型的NM2001程序计算振冲碎石桩,给出了荷载作用下土中碎石桩的体积膨胀与其承载能力提高程度的关联因子；研究了碎石芯体膨胀后的受力规律；多重势面模型能更好地反映碎石桩的作用机理,为这项工程实践提供了更深入的理论依据。