功能:求解岩石轴对称弹塑性区域内塑性区半径(井筒)

格式:
RockAxisSymHBr(R0, Rw, E, v, P0, Pw, a, SSL, Stress0, c, h, cr, hr, bL, Error_)
RockAxisSymHBr(R0, Rw, E, v, P0, Pw, a, SSL, Stress0, c, h, cr, hr, bL)
RockAxisSymHBr(R0, Rw, E, v, P0, Pw, a, SSL, Stress0, c, h, cr, hr)
RockAxisSymHBr(R0, Rw, E, v, P0, Pw, a, SSL, Stress0, c, h, cr)
RockAxisSymHBr(R0, Rw, E, v, P0, Pw, a, SSL, Stress0, c, h)
RockAxisSymHBr(R0, Rw, E, v, P0, Pw, a, SSL, Stress0, c)
RockAxisSymHBr(R0, Rw, E, v, P0, Pw, a, SSL, Stress0)
RockAxisSymHBr(R0, Rw, E, v, P0, Pw, a, SSL)
RockAxisSymHBr(R0, Rw, E, v, P0, Pw, a)
RockAxisSymHBr(R0, Rw, E, v, P0, Pw)
RockAxisSymHBr(R0, Rw, E, v, P0)
RockAxisSymHBr(R0, Rw, E, v)
RockAxisSymHBr(R0, Rw, E)
RockAxisSymHBr(R0, Rw)

R0(m):外边界半径
Rw(m):井筒内半径
E(MPa):杨氏弹性模量,默认3000
v:泊松比,默认0.3
P0(MPa):初始孔隙压力,默认10.0
Pw(MPa):井筒内壁压力,默认1.0
a:Biot系数,默认为1.0
SSL:最大膨胀应变,默认为0.01
Stress0(MPa):R0外边界上的径向应力值,默认为22
c(MPa):未失效区域材料的内聚力,默认为4
h(°):未失效区域材料的摩擦角,默认为33.5
cr(MPa):未失效区域材料的内聚力,默认为2
hr(°):未失效区域材料的摩擦角,默认为30
bL(1/MPa):Langmuir等温常数,默认为0.42
Error_:求解半径的控制误差,默认1e-12

原理:
如下图,求解塑性区与弹性区域内指定点的场变量.其中塑性区失效准则采用非线性的Hoek-Brown准则.
{<http://img1.ph.126.net/cUwBz--onDxvHjSl4vhj9Q==/6632240944605398494.png>}

例子:
r=RockAxisSymHBr(100,0.1)//回车得到如下结果,这里假设R0=100,Rw=0.1,其它材料常数采用默认,回车得到塑性区半径r
r =
[ 0.15027527440868 ]