本文对大型通用软件ANSYS进行二次开发，利用APDL语言将符合混凝土损伤破坏的本构关系以材料弹性模量的变化来体现，最终实现重力坝的动力非线性分析。在ANSYS求解过程中，利用重启动技术实现了时程分析条件下求解与后处理的相互进行，通过循环和嵌套循环来实现地震作用下损伤单元的应变软化以及开裂单元的不抗拉状态，最终根据损伤变量D值和弹性模量降低值可直观量化地了解单元在地震动作用下的损伤发展过程和开裂状态，特别是能够模拟出实验中宏观裂缝出现前坝体内外混凝土的损伤部位和损伤程度。利用本文所述的损伤本构关系，计算得出的加速度最值与试验测值比较接近，说明所选用的本构模型可以很好地模拟高混凝土重力坝在强震下的损伤破坏情况。计算与试验的坝体损伤开裂状况都是：在坝颈处，由坝体下游开始出现裂纹，逐渐扩展至上游，最终汇合贯通成裂缝。因此重力坝的薄弱部位是坝体的折坡段儿，此处单元的应力应变最大，产生的损伤断裂最为明显，故在设计和加固时要对折坡部位予以重视。