利用Fortran和Python语言联合开发ABAQUS有限元软件用户子程序,建立土坝非稳定渗流的交互式有限元分析模型。模型利用Fortran语言编写用户子程序,实现了坝体和坝基表面非线性动态荷载的施加功能,同时采用Python语言对子程序进行可视化处理,以交互的方式处理库水位降落参数实现降水过程。该模型利用流固耦合原理和强度折减法,对均质土坝库水位降落过程中的渗流场和应力场进行计算,通过与已有研究成果对比水位下降过程中浸润线最高点高程、坝体瞬态安全系数,并结合十里河水库大坝水位降落具体工程,验证模型的合理性。该研究成果可为库水位变动下的土坝渗流特性和稳定性计算提供方法,以及为中小型土石坝遭遇极端情况的风险分析和应急管理提供参考。

An interactive finite element model for analyzing the unsteady seepage flow of earth dam is established by secondary development of ABAQUS user subroutines in Fortran and Python language. The non-linear dynamic loads on the surface of dam body and dam foundation can be applied in the model by programming user subroutines using the Fortran language. Meanwhile, the water level drawdown is realized by interactively processing the water level drawdown parameters by visualizing the subroutine using the Python language. The seepage field and stress field during water level drawdown of homogeneous earth dam can be calculated by the model based on the principle of fluid-structure interaction and strength reduction method. The validity of the model is verified by comparing the largest elevation of phreatic line and the transient safety factor of the dam with previous studies, and also by case study on the Shilihe reservoir project. The research finding serves as a method for calculating the seepage characteristics and stability of dams in the presence of reservoir water level drawdown, and provides a reference for risk analysis and emergency management of small and medium-sized earth and rock dams under extreme working conditions.