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

格式:
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw, E, v, P0, Pw, a, SSL, Stress0, c, h, cr, hr, bL, Angle)
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw, E, v, P0, Pw, a, SSL, Stress0, c, h, cr, hr, bL)
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw, E, v, P0, Pw, a, SSL, Stress0, c, h, cr, hr)
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw, E, v, P0, Pw, a, SSL, Stress0, c, h, cr)
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw, E, v, P0, Pw, a, SSL, Stress0, c, h)
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw, E, v, P0, Pw, a, SSL, Stress0, c)
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw, E, v, P0, Pw, a, SSL, Stress0)
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw, E, v, P0, Pw, a, SSL)
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw, E, v, P0, Pw, a)
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw, E, v, P0, Pw)
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw, E, v, P0)
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw, E, v)
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw, E)
[z,Rep0]=RockAxisSymHB(r,Rep,R0, Rw)

r(m):需要求解的指定处的半径值,如果为矩阵变量,则对每个半径进行求对应场变量.
Rep(m):塑性区半径,如果此值为空,程序将自动计算塑性区半径
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
Angle(°):剪胀角,默认12

z:返回的一个m*5的矩阵变量,其中m和r当中半径个数一致,每一行对应r的值。其中每一列分别代表:StressRR,StressSS,StrainRR,StrainSS,U
Rep0:返回的塑性区半径大小

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

例子:

r=0.1:0.01:0.2;
[z,rep]=RockAxisSymHB(r,,100,0.1)//回车可得到半径0.1-0.2m之间各点的物理变量值