自上世纪80年代开始，结合一系列国家级和省部级重点科研项目，系统开展土石坝坝体应力变形的数值计算分析方法的研究，自主开发了能考虑Duncan模型，模拟施工期分期填筑、蓄水过程中容重变化、蓄水期心墙固结过程以及接触面的堆石坝二、三维仿真有限元应力变形分析程序。
    针对高土石坝高应力和复杂应力特点，引进了修正的剑桥弹塑性模型、沈珠江双屈服面弹塑性模型、殷宗泽双屈服面弹塑性模型等，并能进行堆石长期变形和湿化变形的精确模拟，土石坝性态仿真模拟技术得到不断完善。
    根据随应力增大堆石料会呈现明显的非线性强度的特点，提出了基于非线性强度指标的坝坡稳定分析方法和相应的安全系数标准，并将其纳入自主开发的土质边坡稳定分析程序STAB中。近年来，又结合糯扎渡等工程进行坝坡稳定的可靠度指标分析，从风险分析的角度确定大坝的安全储备，为大坝分区优化设计提供依据，取得了良好的效果。
    上述分析方法已经成功应用小浪底、糯扎渡、双江口、两河口等国内外数十座大中型土石坝工程的坝体结构与断面分区优化设计研究和咨询工作中，为堆石坝的坝体结构形式与断面分区设计提供科研支撑。
    近年来，结合小浪底斜心墙堆石坝工程，在前期完成的大量的研究工作基础上，针对工程运行出现的问题，追踪该坝的运行，持续进行了多项研究和安全评估工作，掌握了丰富的技术资料，为高水平开展高土石坝科学研究奠定了基础。

      糯扎渡砾石土心墙堆石坝应力变形分析计算模型