基于FLAC3D的三维边坡稳定性强度折减法计算效率改进算法及其应用

强度折减数值分析方法在边坡工程中应用广泛，其中FLAC^3D已成为三维边坡稳定性分析中的主流计算工具。FLAC^3D自带强度折减计算模块，因其计算过程复杂且耗时，一些学者提出了基于"二分法"的强度折减法，但依然不能令人满意。本文利用Python编程语言对FLAC^3D进行二次开发，实现了"二分法"强度折减法的全面优化，包括初始应力场载入优化和计算时步上限优化，建议计算时步上限取初始应力场计算时步的1.50倍，此时既可以保证计算结果的准确性又可以节省计算时间。利用优化后方法对某公路复杂三维填方路基边坡进行了三维强度折减法计算，结果表明，在保证计算准确性的同时，本文方法计算效率较其他改进的强度折减法至少提升50%，较程序自带模块提升了300%以上。

Optimization and Application of FLAC3D Strength-Reduction Computation in Three-Dimension Slope Stability Analysis

The strength-reduction method (SRM) is widely applied in slope engineering projects, and FLAC^3D has become the main tool for 3D slope stability analysis. Although FLAC^3D has SRM intrinsic modular to complete this task, the computing process is complicated and time-consuming, hence the bisectional method is then involved to simplify the process by some researchers, which, however,is defactive in some aspects. An optimized bisectional method is supplied in this paper based on Python, including in-situ stress field reloading and selection of the upper limit of time-steps, which obviously improves computing efficiency. By using this optimized method, a fill slope in road engineering project was analyzed, showing that the strength reduced factor is consistent with those of the other methods and with much higher computing efficiency (ca. 50% and 300% higher than the other improved methods and the intrinsic modular, respectively).