过大的永久变形往往是导致堆石坝震害的主要原因。准确预测堆石体的地震变形在堆石坝的抗震设计中尤为重要。目前广泛使用的等价粘弹塑性模型对堆石体在地震荷载作用下变形的描述具有很大经验性,只能粗略估算地震荷载引起的残余变形,采用循环弹塑性模型可弥补这一缺陷。本文在前人研究的基础上,着重对堆石料在循环荷载作用下的应力应变规律进行了总结,建立了一个实用化的三维弹塑性循环本构模型,提出了基于该模型的数值化方法,开发了基于开源动力有限元软件OpenSees的面板堆石坝动力计算程序。本文取得的研究成果主要包括以下几个方面:第一,建立了一个可合理地反映堆石料在循环荷载作用下应力-应变特性的三维弹塑性本构模型。与目前已有的弹塑性循环本构模型相比,该模型具有以下几个明显优势:1)提出了一个基于修正空间滑动面强度理论的双参数三维非线性强度准则,使复杂循环加载下堆石料强度的描述更为简便合理;2)可反映偏应力循环、球应力循环及二者共同作用下不可逆体应变的累积规律,使循环荷载作用下堆石料体应变的描述更为准确合理;3)可反映循环剪切荷载作用下剪应力-剪应变关系的滞回特性、硬化特性和应变累积特性,使循环荷载作用下堆石料剪应变的描述更为精细合理;4)共有11个物理意义明确的模型参数,均可由振动三轴试验测得的最大剪切模量、模量衰减曲线和残余变形累积曲线等常规试验结果来确定,使该模型更加实用化。第二,提出了该本构模型的数值实现方法,开发了基于开源动力有限元软件OpenSees的用户材料子程序以及针对面板堆石坝动力计算的辅助性前后处理程序,实现了较为高效的数值建模。通过对二维理想坝模型、三维均质坝模型的系列计算和对比分析,完成了对程序有效性的测试。第三,通过三维动力弹塑性有限元分析,较为合理地再现和揭示了紫坪铺高面板堆石坝的震动响应和变形规律,模拟结果与汶川地震后实测残余变形的分布基本上吻合,进一步检验了本文所建立模型及相应数值方法的适用性。