均匀地应力场下圆形隧道静态弹塑性解析方法

针对圆形隧道的弹塑性解析中许多能更好地反映岩石特性的本构模型由于数学求解的困难而得不到广泛应用的问题，进行新的尝试。首先，基于单一曲线和体积不可压缩假设，总结完善一种更为简便的解析方法。其次，运用该解析方法，针对6种不同的本构模型，计算围岩的弹性区和塑性区应力、塑性区半径、围岩平衡方程、围岩自承地应力上下限等，得到精确或近似的解析表达式。针对线性软化本构模型，引入Lambert函数给出塑性区半径R的解析表达式，论证围岩平衡曲线中R的取值范围。针对Nelder非线性软化本构模型，引入幂函数给出拐点横坐标以及拐点处降模量的近似计算式。当峰后为Weibull非线性软化本构模型时，由于whittaker函数无法求反函数和极限，故而提出以Gauss公式拟合整条围岩平衡曲线，再以Gauss反函数的求解得到塑性区半径的新思路。最后，对比分析6种本构模型的解析结果之间的差异后得出，非线弹性情况下得到的围岩自承地应力上限比线弹性情况下高，而应力集中系数低，应力集中位置偏离弹塑区边界位于弹性区内，切向应力再分布曲线上不会出现尖角，这些都与实际更为相符。因此，不论是软岩或是硬岩，使用连续而光滑的本构模型进行弹塑性分析可以得到更为符合实际的结果。