针对圆形隧道的弹塑性解析中许多能更好地反映岩石特性的本构模型由于数学求解的困难而得不到广泛应用的问题,进行新的尝试。首先,基于单一曲线和体积不可压缩假设,总结完善一种更为简便的解析方法。其次,运用该解析方法,针对6种不同的本构模型,计算围岩的弹性区和塑性区应力、塑性区半径、围岩平衡方程、围岩自承地应力上下限等,得到精确或近似的解析表达式。针对线性软化本构模型,引入Lambert函数给出塑性区半径R的解析表达式,论证围岩平衡曲线中R的取值范围。针对Nelder非线性软化本构模型,引入幂函数给出拐点横坐标以及拐点处降模量的近似计算式。当峰后为Weibull非线性软化本构模型时,由于whittaker函数无法求反函数和极限,故而提出以Gauss公式拟合整条围岩平衡曲线,再以Gauss反函数的求解得到塑性区半径的新思路。最后,对比分析6种本构模型的解析结果之间的差异后得出,非线弹性情况下得到的围岩自承地应力上限比线弹性情况下高,而应力集中系数低,应力集中位置偏离弹塑区边界位于弹性区内,切向应力再分布曲线上不会出现尖角,这些都与实际更为相符。因此,不论是软岩或是硬岩,使用连续而光滑的本构模型进行弹塑性分析可以得到更为符合实际的结果。
A lot of constitutive models,representing the properties of rock well,have not been utilized widely because of the difficulty in elastoplastic solution. So,based on the single curve and constant volume hypothesis, an analytical method that can be used more conveniently was summed up and modified firstly. Using this method, the stress of elastic and plastic zones,radius of plastic zone,equilibrium equation and superior and inferior limits of self-support geostress of surrounding rock after excavation were calculated accurately or approximately. In view of linear softening constitutive model,Lambert function was introduced to calculate the radius of plastic zone R and its range in equilibrium curve of surrounding rock proved. For the Nelder nonlinear softening model, an exponential function was used to describe approximately the abscissa of inflection point and softening modulus. When the Weibull model was taken to represent the post-peak mechanical behaviour,since it is difficult to get the inverse function and limit of Whittaker function,a new idea,fitting the whole equilibrium curve with Gauss function and computing the radius of plastic zone by solving the inverse function of Gauss,was proposed. Then, comparison was taken between these results derived from six constitutive models and it is found that the superior limit of self-support geostress of surrounding rock was higher,the stress concentration factor was lower,the position that stress concentrates was in elastic zone,deviating from the interface between elastic and plastic zones,and the wedge angle in the tangential stress redistribution curve disappeared under nonlinear elastic condition,which is consistent with the practice. Therefore,for either soft or hard rock,a more perfect result from elastoplastic solution can be obtained by taking a smooth and continuous constitutive model.