作为土石坝的主要筑坝材料，堆石料在剪切过程中具有明显的剪胀特性。大量的室内试验结果表明，堆石料在受剪时可能发生剪缩，也可能发生剪胀，其剪胀性取决于堆石料的当前状态，包括堆石料的密度和所处的应力状态。目前在土石坝计算分析过程中，一般采用的堆石料本构模型都无法反映土体密度变化带来的影响，对于具有不同密度的同一种堆石料，往往只能通过试验得到不同的参数，当成不同的材料进行模拟，这给土石坝的应力变形计算带来了一定的问题。 近来一些学者针对砂土提出了状态相关剪胀理论，并建立了基于状态相关剪胀理论的临界状态模型，成功地模拟了砂土的剪胀特性。堆石料与砂土都是颗粒集合而成的粒状材料，其工程性质类似。基于此本文在阅读大量文献的基础上，通过室内试验、理论分析、数值模拟等手段，对堆石料的变性特性与本构理论进行了深入系统的研究，主要工作如下： （1）通过室内试验系统地研究了两河口心墙堆石坝坝体堆石料的变性特性和临界状态，揭示了堆石料变形的密度与应力状态相关性：对具有不同密度和应力状态的堆石料进行大型三轴排水压缩试验，研究了密度和应力水平对堆石料变形特性的影响，得到了在较大剪切变形时该堆石料的临界状态线与临界状态参数。通过试验，还获得了两河口心墙堆石坝坝体堆石料的南水模型参数、邓肯-张模型参数和基于状态相关剪胀理论的临界状态模型参数。 （2）基于堆石料大型三轴排水压缩试验结果，采用基于状态相关剪胀理论的临界状态模型以及获得的一套模型参数，编写了三轴排水试验研究程序（点模型程序），对所有试验结果进行了数值模拟，并与南水模型和邓肯-张模型的模拟结果进行了对比。研究表明，对于处于不同密度不同应力水平的同一种堆石料，基于状态相关剪胀理论的临界状态模型只需一套模型参数就能反映堆石料的各种特性，包括剪胀性与应变软化特性等，而南水模型和邓肯．张模型对于不同密度的同一种堆石料需要多组参数来模拟，即将具有不同状态的同一种堆石料当成不同的材料来处理。这充分显示了基于状态相关剪胀理论的临界状态模型在模拟堆石料变性特性方面的合理性。 （3）对SUMDES2D有限元程序进行了二次开发，嵌入了邓肯，张非线性模型和南水双屈服面弹塑性模型程序，拓展了其应用范围；分别采用基于状态相关剪胀理论的临界状态模型、南水模型和邓肯-张模型，通过SUMDES2D有限元程序对三轴试验结果进行了模拟，并与点模型程序模拟结果进行了对比。此外还对单元的敏感性进行了分析，验证了SUMDES2D软件的可靠性；分别采用基于状态相关剪胀理论的临界状态模型、南水模型和邓肯-张模型，通过SUMDES2D有限元程序对某一均质堆石坝进行了模拟，并与ABAQUS大型商业软件模拟结果进行了比较，进一步验证了SUMDES2D程序在解决岩土工程问题上的适用性。 （4）采用基于状态相关剪胀理论的临界状态模型，通过SUMDES2D有限元程序对长河坝进行了二维静动力有限元计算，得到了坝体施工期和蓄水期的水平位移、沉降以及坝体应力分布情况，揭示了堆石坝在地震作用下变形和破坏的机理。