本文针对直立重力坝面提出了一种新的高阶精度时域动水压力公式，将此公式嵌入了ABAQUS有限元软件中，并利用开发的计算软件分析了Koyna混凝土重力坝非线性地震反应。主要研究工作和成果如下:　　(1)提出了一种计算直立坝面上地震动水压力的新的高阶精度时域公式。首先在频域内基于库水区域分离变量解析解和坝-水交界面条件推导出精确的坝面动水压力-坝面绝对水平运动关系，进一步利用时间卷积局部化方法将该关系转化为含有辅助变量的高阶精度时域公式。刚性坝面地震动水压力分析表明，提出的公式计算结果稳定，并与解析解吻合较好。　　(2)将提出的动水压力公式的有限元方程转化为物理模型系统，然后通过ABAQUS用户单元子程序(User Element)结合自编Fortran代码开发了基于该物理模型的有限元计算软件。该计算软件在动力计算过程中无需重启动等任何特殊处理，整个实现过程与结构有限元模型无缝结合。现有的实现方式是在商用有限元软件中基于重启动技术进行有限元模型和高阶精度时域动水压力公式交替求解，每一时步均需重启动计算，计算效率低，本文实现方式避免了该问题。　　(3)应用开发的计算软件分析了Koyna混凝土重力坝的非线性地震反应。在分析模型中，采用了ABAQUS内置的塑性损伤本构模型模拟坝体混凝土材料非线性，并考虑了坝体自重与静水压力作用。在线弹性和非线性坝体分析中，将高阶精度时域动水压力公式的计算结果与Westergaard附加质量公式和空库计算结果进行了对比分析。