Hoek-Brown准则奇异屈服面的圆化方法及其强度折减技术与应用

 为克服基于Hoek-Brown强度准则的屈服面数值奇异问题，提出新的棱角处屈服面圆化方法，在此基础上，借助ABAQUS二次开发平台，采用FORTRAN语言编制理想弹塑性模型子程序，应力更新算法为完全隐式的向后欧拉积分算法； 提出基于Hoek-Brown准则的有限元强度折减技术，并将其嵌入到圆化后的屈服函数中，探讨了强度折减技术应用于Mohr-Coulomb强度准则和Hoek-Brown非线性强度准则的异同点。对某圆形隧洞的模拟分析表明，围岩应力和塑性区的数值计算结果和解析结果吻合很好；对某岩质边坡的模拟分析表明，与等效Mohr-Coulomb模型的线性强度折减技术相比，基于Hoek-Brown模型的非线性强度折减技术获得的潜在滑动面位置较浅、形状较陡，更符合岩质边坡滑动面的特征，两个算例验证了该圆化方法和非线性强度折减技术的正确性和适用性。

A rounding approach the singular surface of Hoek-Brown criterion and its strength reduction technique

To overcome the numerical singularity of the yield surface based on the Hoek-Brown strength criterion， a new approach was presented to smooth the edges of the yield surface. The ideal elasto-plastic constitutive model was implemented in ABAQUS software with the fully implicit backward Euler integral algorithm employed in the stress updating. The strength reduction method in finite element method based on the Hoek-Brown criterion was implemented with the smoothed yield function. The similarities and differences of the strength reduction techniques applied with the Mohr-Coulomb strength criterion and the Hoek-Brown nonlinear strength criterion were discussed. The modeling on a circular tunnel shows that the numerical results of stresses and plastic zones in the surrounding rock agree well with the analytical results. The modeling on a rock slope shows that in comparison with the linear strength reduction technique of the equivalent Mohr-Coulomb model，the nonlinear strength reduction technique based on the Hoek-Brown model produced the slip surface which agree better with the characteristics of the slip surface of a rock slope，which are shallower in location and steeper in shape.