网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

饱和结构性黄土不排水柱孔扩张问题弹塑性解  PDF

  • 李镜培 1,2
  • 周攀 1,2
  • 李亮 1,2
  • 谢峰 1,2
  • 崔纪飞 1,2
1. 同济大学 岩土及地下工程教育部重点实验室,上海 200092; 2. 同济大学 土木工程学院,上海 200092

中图分类号: TU473

最近更新:2021-03-17

DOI:10.11908/j.issn.0253-374x.20237

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
目录contents

摘要

为导出饱和黄土中不排水柱孔扩张问题严格解,采用结构性黄土修正剑桥模型描述黄土力学特性。在弹性区采用小变形假设求解孔周土体应力与变形,同时将塑性区柱孔扩张边值问题转化为3个以应力分量为基本未知量的一阶非线性常微分方程求解问题。以弹塑性界面上应力分量作为初始条件求解控制方程,得到孔周塑性区应力分布,然后对径向平衡方程积分求得孔隙水压力。通过算例分析验证了该方法的正确性。结果表明:退化解与已有数值解完全吻合,并且黄土结构性对孔周应力、超静孔压以及弹塑性半径等影响显著。

西部大开发战略的实施使黄土分布最广的西北地区岩土工程问题研究变得十分紧

1-2。尤其是在流域低阶地区,因地质环境改变、地下水位上升而形成了饱和黄土,此类黄土虽然有着与软土相近的某些特征,但是具有独特的结构性,因此对饱和黄土地区的岩土工程问题研究十分重要。对于许多复杂岩土工程问题,如原位土工试验(旁压测试、静力触探测试等)、桩基工程、隧道和地下工程以及井壁失稳等,圆孔扩张理论是既简单又适用的分析工3-7

在岩土体中的圆孔扩张或收缩实质上是一个一维边值问题,其求解往往需要一个合理描述岩土材料应力‒应变关系的本构模型。理想弹塑性本构模型(Mohr‒Coulomb理论、Tresca准则和统一强度理论等)和临界状态应变硬化或软化本构模型(剑桥模型、修正剑桥模型以及众多用于特殊土体的改进修正剑桥模型)是目前广泛应用于反映土体力学特征的2类本构模型。圆孔扩张理论应用初期,许多学者采用摩尔‒库伦准则求解柱孔和球孔扩张问题。Vesic

8基于土体塑性屈服后的体积改变推导出了柱孔扩张和球孔扩张问题的近似解;Carter9采用非关联流动法则描述土体塑性流动,求得了小应变下孔压与孔扩张关系的显式解,同时求得了有限变形下孔压与圆孔扩张关系的数值解;邹金锋10同样采用非关联流动法则并考虑体变守恒和大变形理论推导了砂性土中极限扩孔压力表达式。同时,其他理想弹塑性模型适用性也较为广泛,丰富了岩土工程中圆孔扩张问题的解答。饶平平11假定土体服从Lade‒Duncan屈服准则,采用大变形和非关联流动法则对圆柱孔扩张问题进行了解答;曹黎娟12、綦春明13在弹性区采用小变形理论,在塑性区采用统一强度准则与大变形理论,推导出了扩孔问题的塑性区半径、极限扩孔压力以及超静孔压的解析式。

理想弹塑性模型无法考虑土体应力历史与超固结比的影响,因此临界状态模型逐渐受到广大学者的青睐而广泛应用于柱孔扩张问题的研究。由于临界状态模型自身的复杂性,因此在使用该类模型时或多或少会进行适当简化,其中包括假定塑性区内剪切模量恒定与偏应力线性分

14、应力路径线性简15-17以及平均塑性应变值估18等,以便求解柱孔和球孔扩张问题解析解。为了反映土体各向异性、强度三维力学特性对扩张问题计算结果的影响,李镜培19-20采用各向异性修正剑桥模型和基于空间滑动面(SMP)准则改进的修正剑桥模型分析了不排水和排水条件下柱孔扩张机制。为了得到柱孔扩张问题的精确解,Chen21-22采用修正剑桥模型,在未对平均有效应力和偏应力做任何简化的前提下求得了柱孔不排水与排水扩张半解析解。随后,郑金辉23在此求解框架基础上采用考虑砂土颗粒破碎的本构模型,分析了砂土圆孔扩张问题。

对于天然的岩土材料,特别是广泛分布于我国西北地区的黄土,其结构性对土体强度和变形特性的影响是土力学最本质的一个因素。然而,以上研究大部分是针对于求解饱和黏土、软土和砂土等的柱孔扩张问题,现有解析解中对考虑土体结构性的扩孔机理分析仍存在不足,尤其是黄土地区的圆孔扩张问题解答较少。因此,基于Chen

21的分析框架,在遵循2个有效应力不变的严格定义前提下,通过采用邓国华24-25提出的黄土结构性修正剑桥模型来描述土体塑性屈服后的材料强度和变形特性,为结构性黄土的圆孔排水扩张边值问题求解推导出精确的半解析解,最终得到孔周围土体精确的有效径向、切向和竖向应力以及超静孔压的分布规律。

1 不排水柱孔扩张问题力学模型

1.1 基本假设与模型建立

图1为无限空间内的饱和黄土柱孔扩张力学模型。土体受到的初始水平应力与初始竖向应力分别为σh0'σv0',且σh0'=K0σv0',其中K0为静止侧压力系数。假设柱孔初始半径为a0,初始扩孔压力为σ0,扩孔稳定后柱孔半径为a,同时扩孔压力增加到σa。当扩张完成后,孔周任意一点土体从原始位置rx0移动到当前位置rx,其位移用Ux表示,孔周塑性区半径从原始位置rp0移动到当前位置rp,其位移用Up表示。当孔压增大到一定值后,孔周土体应力处于3种状态:孔壁附近为临界状态区(arrf,其中rf为破坏区半径),应力不再改变,应变不断发展;向孔外延伸一定厚度的土体为塑性区(rf < r < rp),随着孔压增大应力应变仍在改变;距离孔壁较远处为弹性区(rrp),随着距孔壁距离越大,土体应力状态越接近初始应力值。

图1 结构性黄土中的柱孔扩张力学模型

Fig.1 Mechanical model of cylindrical cavity expansion in structured loess

柱孔扩张问题的求解前提是基于以下假设:

(1)土体为均质、连续且各向不等压的饱和黄土。

(2)孔隙水不可压缩,应力和应变以压为正。

(3)平面应变和不排水条件下柱孔扩张,因此扩孔过程中土体体积不变。

(4)在弹性区土体服从线弹性和小变形假定,在塑性区采用大变形理论和结构性黄土修正剑桥模型描述土体的力学特性。

孔周任意土体微单元总应力满足以下平衡微分方程:

dσrdr+σr-σθr=0 (1)

式中:σrσθ分别为孔周土体总径向应力和切向应力。式(1)还可表示为有效应力形式,如下所示:

dσr'dr+ur+σr'-σθ'r=0 (2)

式中:σr'σθ'分别为孔周土体有效径向应力和切向应力;u为孔隙水压力。

1.2 结构性黄土本构模型

在孔周弹性区采用胡克定律描述土体应力状态与变形特征,应力‒应变关系的增量形式如下所示:

dεij,e=1+νEdσij'-νEdσmm'δij (3)

式中:E为弹性模量,E=2G(1+ν),其中ν为泊松比,G为剪切模量;εij,e为应变分量;σij'为应力分量;σmm'为主应力分量;δij为克罗内克函数。G的计算式如下所示:

G=3(1-2ν)υp'2(1+ν)κ (4)

式中:υ为土体比容;p'为有效平均应力;κ为回弹曲线斜率。

在孔周塑性区采用邓国华

24-26、邵生俊27提出的结构性黄土修正剑桥模型,如下所示:

f=q2-M2p'pC'-p'-CM-MCpC'=0 (5)

式中:pC'为平均有效屈服应力;C为综合结构性参数,反映应力比结构性参数mη对屈服面平移量和形状尺寸的改变。综合结构性参数C、临界状态应力比M、有效平均应力p'以及偏应力q的计算式如下所示:

C=33c(mη)cos φ(mη)3sin θσsin φ(mη)+3cos φ(mη) (6)
M=33sin φ(mη)3sin θσsin φ(mη)+3cos φ(mη) (7)
p'=13(σr'+σθ'+σz') (8)
q=12(σr'-σθ')2+(σr'-σz')2+(σθ'-σz')2 (9)

式中:σz'为孔周土体有效竖向应力;θσ为应力罗德角;cmη)、φmη)分别为黄土中随应力比结构性参数变化的黏聚力和内摩擦角。cmη)和φmη)的计算式分别为

c(mη)=mηa+bmη-Aφ(mη)=φ (10)

式中:a、b、A为试验拟合参数;φ为土体内摩擦角。mη的表达式如下所

24

mη=ηi2/(ηrηs),(prps)/pi2,qi2/(qrqs), (11)

式中:下标i、r、s分别表示原状土、重塑土和饱和土;ηiηrηs分别为原状土、重塑土和饱和土在一般应力路径下的应力比;piprps分别为原状土、重塑土和饱和土在一定体应变条件下的等向压应力或球应力;qiqrqs分别为原状土、重塑土和饱和土在某一剪切变形条件下的剪应力。

考虑应变硬化规律后,式(5)可进一步表示为

f=λ-κ1+e0lnp'p0'+λ-κ1+e0ln1+q2M2p'2-εv,p-λ-κ1+e0ln(Mp'+C)(M2p'2+q2)(M2p'2+MCp'+q2)Mp'=0 (12)

式中:λ为土体压缩曲线斜率;p0'为初始有效平均应力;e0为土体初始孔隙比;εv,p为塑性体应变。当综合结构性参数C = 0时,式(12)退化为修正剑桥模型屈服方程表达式。

2 柱孔扩张问题半解析解

2.1 弹性区分析

由Chen

21-22分析可知,根据小变形理论、胡克定律与孔周土体微单元的应力平衡方程可推得孔周弹性区应力分量、位移与超静孔隙水压力的解析式,如下所示:

σr'=σh0'+(σp'-σh0')rpr2 (13)
σθ'=σh0'-(σp'-σh0')rpr2 (14)
σz'=σv0' (15)
Ur=σp'-σh0'2G0rp2r (16)
Δu=0 (17)

式中:σp'为弹塑性界面上的有效径向应力;G0为初始剪切模量。G0的计算式如下所示:

G0=3(1-2ν)υ0p0'2(1+ν)κ (18)

式中:υ0为土体初始比容。

2.2 弹塑性界面分析

对于初始应力各向异性的结构性黄土,孔周初始有效径向应力σr0'、切向应力σθ0'和竖向应力σz0'之间的关系如下所示:

σr0'=σθ0'=σh0-u0=K0σz0'=K0(σv0-u0) (19)

式中:σh0σv0分别为总初始水平、竖向应力;u0为初始孔隙水压力。式(19)可转化为

σr0'p0'=σθ0'p0'=3K01+2K0,σz0'p0'=31+2K0 (20)

由弹性区应力解析式可知,扩孔过程中孔周土体的有效平均应力p'保持不变,则根据应力连续条件弹塑性界面处有效平均应力pp'可表示为

pp'=p0' (21)

将(pp'q0)代入式(5)屈服方程中,可推得初始加载点所对应的当前屈服面与平均正应力p轴的交点pA'pA'的计算式如下所示:

pA'=q02+M2p0'2+MCp0'M2p0'+MC (22)

式中:q0为初始有效偏应力。

将弹塑性界面上点(p0'qp)代入式(5)屈服方程中,并联立式(21)式(22),可解得弹塑性界面上偏应力

qp=Mp0'Rq02+M2p0'2+MCp0'M2p0'+MC1p0'+CMp0'2-CMp0'-1 (23)

式中:qp为弹塑性界面上的偏应力;R为结构性黄土超固结比。R的计算式如下所示:

R=pC'pA' (24)

弹塑性界面上的有效径向应力σrp'、切向应力 σθp'、竖向应力σzp'和偏应力qp分别表示为如下关系式:

σrp'+σθp'=σr0'+σθ0'=6K01+2K0p0' (25)
σzp'=σz0'=31+2K0p0' (26)
qp=12(σrp'-σθp')2+(σrp'-σzp')2+(σθp'-σzp')2 (27)

结合式(25)~(27),可解得弹塑性界面上应力分量

σrp'=σr0'+σr0'2-13(4σr0'2+σz0'2-2σr0'σz0'-qp2) (28)
σθp'=σr0'-σr0'2-13(4σr0'2+σz0'2-2σr0'σz0'-qp2) (29)
σzp'=11+2K0(σr0'+σθ0'+σz0') (30)

设刚进入塑性状态的微单元位置为rxp,其初始位置为rx0,由式(16)可得该土体单元的径向位移

Urp=rxp-rx0=σrp'-σr0'2G0rxp (31)

对于不排水柱孔扩张,土体单元体积不变,因此孔周任意一点的当前径向坐标rx和初始径向坐标rx0与圆孔当前半径a和初始半径a0之间的关系可表示为

rx2-rx02=a2-a02 (32)

结合式(31)式(32)可推得柱孔不排水过程中孔周任意一点土体进入塑性状态时径向位置

rxp=2G0a2G0-(σrp'-σr0')rxa2+a0a2-1 (33)

rx = rxp = rp,则由式(33)可得柱孔不排水扩张过程中孔周塑性区半径rp与当前扩张半径a之间的关系为

rpa=a0a2-1σrp'-σr0'2G0-12-1 (34)

2.3 塑性区分析

孔周土体塑性屈服后总应变增量由弹性体应变增量和塑性体应变增量两部分组成,如下所示:

dεij=dεij,e+dεij,p (35)

式中:ε为应变。土体屈服后弹性体应变增量任然采用胡克定律计算,由结构性黄土剑桥模型和相关联流动法则,可通过下式计算塑性体应变增量:

dεij,p=Λfσij' (36)

式中:Λ为塑性因子,表示塑性应变增量大小的一个标量;f为屈服函数。

塑性体应变增量dεv,p与主应力方向上的塑性应变增量dεr,p、dεθ,p和dεz,p满足

dεv,p=dεr,p+dεθ,p+dεz,p (37)

联立式(36)式(37),由塑性一致性条件可将结构性黄土修正剑桥模型的塑性因子表示为

Λ=-(f/p')dp'+(f/q)dq(f/εv,p)(f/σij') (38)
fp'=λ-κ1+e01p'C(M2p'2+MCp'+q2)(M2p'2+q2)(Mp'+C)+                      M2p'2-q2M2p'2+q2 (39)
fq=λ-κ1+e02qM2p'2+q2 (40)
fεv,p=-1 (41)
fσij'=fp'p'σij'+fqqσij' (42)
p'σij'=13δij (43)
qσij'=3(σij'-p')2q (44)

有效平均应力增量dp'与偏应力增量dq可表示为

dp'=13(dσr'+dσθ'+dσz') (45)
dq=qσr'dσr'+qσθ'dσθ'+qσz'dσz' (46)

联立式(27)~(35),孔周塑性区应力‒应变关系式采用对称矩阵形式表示,如下所示:

dεr,pdεθ,pdεz,p=yar2araθarazaθaraθ2aθazazarazaθaz2dσr'dσθ'dσz' (47)
y=λ-κ(1+e0)(p'3(M4-η4)+p'(2MCp'+C2)(M2+η2))(1+C/Mp') (48)
ar=p'(M2-η2)+2MC+C2/p'3+3(σr'-                    p')1+CMp' (49)
aθ=p'(M2-η2)+2MC+C2/p'3+3(σθ'-                     p')1+CMp' (50)
az=p'(M2-η2)+2MC+C2/p'3+3(σz'-                      p')1+CMp' (51)

结合孔周弹性区应力‒应变关系式(3)和塑性区应力‒应变关系式(47),能够得到用矩阵表示的弹塑性本构方程,如下所示:

dεrdεθdεz=yarar+1Eyaraθ-νEyaraz-νEyaθar-νEyaθaθ+1Eyaθaz-νEyazar-νEyazaθ-νEyazaz+1Edσr'dσθ'dσz' (52)

式中:εrεθεz分别为径向、切向和竖向应变。

式(52)转化为如下形式:

dσr'dσθ'dσz'=1Δb11b12b13b21b22b23b31b32b33dεrdεθdεz (53)
b11=1E2(1-ν2+Eaθ2y+2Eνaθazy+Eaz2y) (54)
b12=1E2(-Ear(aθ+νaz)y+ν(1+ν-
 Eaθazy+Eaz2y)) (55)
b13=1E2(-Ear(νaθ+az)y+ν(1+ν-
Eaθazy+Eaθ2y)) (56)
b22=1E2(1-ν2+Ear2y+2Eνarazy+Eaz2y) (57)
b23=1E2(ν+ν2+Eνar2y-Eaθazy-
                          Eνar(aθ+az)y) (58)
b33=1E2(1-ν2+Ear2y+2Eνaraθy+Eaθ2y) (59)
b21=b12 (60)
b31=b13 (61)
b32=b23 (62)
Δ=-1+νE3((-1+ν+2ν2)+E(-1+ν)ar2y+E(-1+ν)aθ2y-2Eνaθazy-Eaz2y+
           Eνaz2y-2Eνar(aθ+az)y) (63)

基于大变形理论,孔周塑性区r处的径向应变增量dεr和切向应变增量dεθ的拉格朗日描述形式为

dεr=-(dr)r (64)
dεθ=-drr (65)

柱孔扩张属于典型的平面应变问题,因此竖向应变增量亦为零,即

dεz=0 (66)

根据不排水条件可知,柱孔扩张过程中土体体应变增量为零,即

dεv=dεr+dεθ+dεz (67)

联立式(64)~(67),可得

dεr=-dεθ=drr (68)

式(68)代入式(53),可将柱孔扩张过程的本构方程转化为如下形式的一阶非线性常微分方程组:

dσr'dr-b11-b12Δr=0 (69)
dσθ'dr-b21-b22Δr=0 (70)
dσz'dr-b31-b32Δr=0 (71)

注意到b11b12b21b22b31b32均为关于3个应力分量σr'σθ'σz'的显式函数,因此能够通过3个应力分量在弹塑性界面上的初始值对式(69)~(71)进行数值求解。

孔周塑性区内任意位置处超静孔隙水压力Δurx可通过对式(2)积分得出,如下所示:

Δurx=σrp'-σrx'+rxrpσr'-σθ'rdr (72)

式中:σrx'为当前计算点处的有效径向应力。

3 解析解验证与讨论

上述解析式和求解程序为结构性黄土不排水柱孔扩张问题提供了严格的理论框架,然而对于控制方程式(69)~(71)无法直接得到显式解析解,必须通过数值方法求解。首先,通过Matlab数值软件对Chen

21的算例参数进行编程求解,用于验证本方法的正确性;然后,研究初始结构性与超固结比对柱孔扩张的影响与孔周土体应力和超孔隙压力分布。为了与已有计算结果进行对比以研究土体结构性在扩孔过程中对土体受力和变形的影响规律,采用文献[21]中算例参数,并且分别对超固结比为1.2、3.0和10.0的土体假定综合结构性参数C为0、50、100 kPa,计算参数如表1所示,其中υcs表示土体正应力为1 kPa时的土体比容。

表1 计算参数
Tab.1 Calculated parameters (M=1.2, λ=0.15, κ=0.03, ν=0.278, υcs=2.740)
Rσr0'/kPaσθ0'/kPaσz0'/kPaυ0K0G0/kPaC/kPa
1.2 100 100 160 2.06 0.625 4 302 0
3.0 120 120 120 1.97 1.000 4 113 50
10.0 144 144 72 1.80 2.000 3 756 100

表1中参数代入第2节公式中,通过Matlab数值方法,可得到孔周土体应力和超静孔压分布、孔壁处扩孔压力、超静孔压以及应力路径,计算结果如图2~4所示。图2~4中,虚线表示Chen

21基于修正剑桥模型(未考虑结构性影响)的计算值,实线为本方法在相同土体参数中考虑结构性影响的解。当C = 0时,本方法的解退化为Chen21的解,间接验证本方法的正确性。当C ≠ 0时,由本方法解与Chen21解的差异性来分析土体结构性对沉桩扩孔力学机理的影响。

图2 归一化的扩孔压力与超静孔压随扩孔半径的变化

Fig.2 Variation of normalized internal cavity pressure and excess pore pressure with cavity radius

图3 不同超固结比结构性黄土的有效径向、切向与竖向应力分布(a/a0=2)

Fig.3 Effective radial, tangential and vertical stress distributions for different over-consolidated structured loesses (a/a0=2)

图4 不同超固结比结构性黄土的有效平均应力、偏应力与超静孔压分布(a/a0=2)

Fig.4 Effective mean stress, deviator stress and excess pore pressure distributions for different over-consolidated structured loesses (a/a0=2)

图2为不同先期固结压力与结构性参数条件下孔壁总应力σa、孔壁超静孔压Δua)(通过除以不排水强度su进行归一化处理)与扩孔半径之间的关系。

图2所示,在扩孔过程中随着扩孔半径的增加扩孔压力与孔壁超静孔压均随之增大,当a/a0介于1和2之间时,σa/su与Δua)/su均变化剧烈,直到a/a0超过2时两者逐渐趋于稳定,表明土体已达到临界状态,此时应力不再变化而剪应变持续发展。还可以看出,归一化的扩孔压力和孔壁超静孔压随超固结比增加而减小。由此可知,无论是否考虑土体结构性影响,扩孔压力、孔壁超静孔压随扩孔半径的变化规律都与Chen

21的分析结果一致,当结构性参数C = 0时,两者的计算结果完全相同,因此Chen21的计算结果为本解析解的一种特殊情况(C = 0)。还可知,当考虑土体结构性时,扩孔压力σa/su、孔壁超静孔压Δua)/su均大于未考虑结构性时的计算结果,但是当土体受到的先期固结压力较大时,扩孔压力σa/su、孔壁超静孔压Δua)/su的变化幅值受结构性影响减弱。

图3为不同超固结比和结构性参数下孔周有效径向、切向和竖向应力分布规律。如图3所示,无论是否考虑结构性影响,扩孔稳定后孔周土体分为弹性、塑性和破坏3种状态,靠近孔壁区域的土体达到了临界状态而破坏,应力不再变化,塑性区土体应力变化急剧,弹性区土体随着距孔壁越远应力状态越接近初始应力状态。由图3a可知,当忽略土体结构性影响时,本方法计算结果与Chen

21的计算结果一致,表明采用修正剑桥模型分析柱孔扩张是本方法的特例(C = 0)。由图3b、3c可知,当考虑土体结构性影响时,相比于Chen21的计算结果,本方法计算的塑性区半径与破坏区半径均减小,并且有效径向、切向和竖向应力也均减小,但有效径向应力相比于有效切向应力减小的幅度要小,即(σr' - σθ')增大,表明扩孔过程中考虑结构性参数影响提高了土体抵抗剪切变形的能力,还可知修正剑桥模型的计算结果低估了扩孔过程结构性土体的强度。

图4为不同超固结比和结构性参数下孔周土体有效平均应力p'、偏应力q与超静孔压Δu的分布规律。

图4可以看出,无论是否考虑结构性影响,在破坏区,有效平均应力与偏应力均不变,但随着距孔壁越近超静孔压呈线性急剧增大。在塑性区,随着距孔壁越近,有效平均应力变化较小,对于轻度超固结土(R = 1.2),p'逐渐减小且q逐渐增加到临界值,超静孔压急剧增大,但对于中度超固结土(R = 3.0)和重度超固结土(R = 10.0),p'逐渐增大且q先增大到最大值再逐渐减小至临界值,并且超固结比越大,此种变化趋势越明显。在重度超固结土中超静孔压先减小至最小值再增大,表明在重度超固结土体中出现了应变软化特性伴随着剪胀效应,但是考虑结构性影响能够减小土体的剪胀效应,如图4c所示。在弹性区超静孔压与有效平均应力保持不变,随着距孔壁越近,偏应力逐渐增加。由图4a可知,当忽略土体结构性影响时,本方法计算结果与Chen

21计算结果一致,表明采用修正剑桥模型分析柱孔扩张是本研究的特例(C = 0),这与上述分析一致。由图4b、c可知,当考虑土体结构性影响时,相比于Chen21的计算结果,本方法计算的超静孔压与偏应力均增大,有效平均应力减小,但是超固结比增大,三者的变化幅值均减小,表明在重度超固结下的结构性土体中,结构性参数对土体的应力与变形影响减弱,这与图3中分析结果一致。

图5为不同超固结比与结构性参数影响下孔壁处黄土微单元在p'‒q平面上的有效应力路径(ESP),图中EFG分别表示初始加载点、初始屈服点和临界状态点。

图5 孔壁处不同超固结比结构性黄土在p'q平面上的应力路径(a/a0=2)

Fig.5 Stress path in p'q plane around cavity wall for different over-consolidated structured loesses (a/a0=2)

图5a可知,对于轻度超固结土(R = 1.2),采用的黄土结构性修正剑桥模型的初始屈服面与Chen

21采用的修正剑桥模型的初始屈服面相同,并且土体应力路径也一致,进一步说明当不考虑结构性影响时本方法解析解退化为Chen21的解。对于中度超固结土(R = 3.0)和重度超固结土(R = 10.0),从图5b、c可以看出,考虑土体结构性参数时,两者初始屈服面均向左平移,并且屈服面形状发生变化,临界状态线也向左平移相同距离且斜率不变。无论是否考虑土体结构性影响,在3种超固结土体扩孔时孔壁土体均先经历弹性阶段(p'不变,q增加),在应力点运动到初始屈服面并产生塑性体应变后则进入塑性阶段(p'增大或减小,q增加),直到应力点达到临界状态线后进入破坏阶段(p'不变,q不变)。

4 结论

(1)考虑了土体结构性,推导出的半解析解能够更加真实地反映扩孔过程中土体应力状态。与采用修正剑桥模型的计算结果相比,本方法中结构性黄土的初始屈服面向左平移并且椭圆尺寸发生改变,平移量与尺寸改变量取决于结构性参数的取值。当忽略结构性影响时,本方法计算结果与Chen

21提出方法的计算结果一致。

(2)孔周土体分为3种状态:外侧土体处于弹性状态,并且距孔壁越远土体应力状态越接近扩孔前的初始状态;中间土体处于塑性状态,有效应力分量与超静孔压急剧变化;内侧土体处于临界状态,土体有效应力保持不变,但超静孔压仍变化较大。

(3)土体结构性与超固结比对孔周应力与变形影响显著,弹塑性半径与临界区半径均随着超固结比和综合结构性参数增大而减小。对于轻度与中度超固结土,当产生塑性体应变后土体出现应变硬化现象。当超固结比较大时出现负超静孔压,表明出现塑性软化现象,但考虑结构性能够减小土体软化程度。

作者贡献声明

李镜培:给出论文研究方法,稿件的审核和监督。

周 攀:给出论文研究方法,软件实施,研究结果验证,稿件初稿写作和返修。

李 亮:稿件审阅和修改。

谢 峰:稿件审阅和修改。

崔纪飞:稿件审阅和修改。

参考文献

1

李宁程国栋谢定义. 西部大开发中的岩土力学问题[J]. 岩土工程学报2001233): 268. [百度学术

LI NingCHENG GuodongXIE Dingyi. Geomechanics development in civil construction in Western China[J]. Chinese Journal of Geotechnical Engineering2001233): 268. [百度学术

2

张伟杜太生张彬. 钢渣桩在湿陷性黄土地基中的应用探讨[J]. 哈尔滨工业大学学报2003356): 698. [百度学术

ZHANG WeiDU TaishengZHANG Binet al. Discussion on application of steel slag pile in collapsible loess foundation[J]. Journal of Harbin Institute of Technology2003356): 698. [百度学术

3

李林李镜培龚卫兵. K0固结天然饱和黏土中柱孔扩张弹塑性解[J]. 哈尔滨工业大学学报2017496): 90. [百度学术

LI LinLI JingpeiGONG Weibinget al. Elasto-plastic solution to expansion of a cylindrical cavity in K0-consolidated natural saturated clay[J]. Journal of Harbin Institute of Technology2017496): 90. [百度学术

4

CHEN S LABOUSLEIMAN Y N. Drained and undrained analyses of cylindrical cavity contractions by bounding surface plasticity[J]. Canadian Geotechnical Journal2016539): 1398. [百度学术

5

LIU KCHEN S LGU X Q. Analytical and numerical analyses of tunnel excavation problem using an extended Drucker-Prager model[J]. Rock Mechanics and Rock Engineering2020531): 1777. [百度学术

6

张亚国李镜培胡志平. 基于CPTU确定软黏土超固结比的理论改进方法[J]. 岩石力学与工程学报2017366): 1488. [百度学术

ZHANG YaguoLI JingpeiHU Zhiping. Modified theory approach to determination of OCR from piezocone penetration tests in soft clays[J]. Chinese Journal of Rock Mechanics and Engineering2017366): 1488. [百度学术

7

彭文祥张旭曹佳文. 充气锚杆极限承载力计算方法[J]. 岩土力学2013346): 1696. [百度学术

PENG WenxiangZHANG XuCAO Jiawen. Calculation method for ultimate bearing capacity of inflatable anchor[J]. Rock and Soil Mechanics2013346): 1696. [百度学术

8

VESIC A S. Expansion of cavities in infinite soil mass[J]. Journal of Soil Mechanics and Foundation Engineering1972983): 265. [百度学术

9

CARTER J PBOOKER J RYEUNG S K. Cavity expansion in cohesive frictional soils[J]. Géotechnique1986363): 345. [百度学术

10

邹金锋罗强李亮. 考虑体变与大变形时高应力砂性土扩孔问题能量分析方法比较[J]. 岩石力学与工程学报2008272): 426. [百度学术

ZOU JinfengLUO QiangLI Lianget al. Comparison of energy dissipation analytical methods for cavity expansion considering large deformation and volumetric change in sand soil under high stresses[J]. Chinese Journal of Rock Mechanics and Engineering2008272): 426. [百度学术

11

饶平平王道远. 基于广义Lade‒Duncan准则的圆柱孔扩张问题分析[J]. 地下空间与工程学报201285): 946. [百度学术

RAO PingpingWANG Daoyuan. Analysis of cylindrical cavity expansion based on extended Lade-Duncan criterion[J]. Chinese Journal of Underground Space and Engineering201285): 946. [百度学术

12

曹黎娟赵均海魏雪英. 基于统一强度理论的灰土挤密桩应力分析[J]. 岩土力学20062710): 1786. [百度学术

CAO LijuanZHAO JunhaiWEI Xueying. Stress analysis of lime-soil compacted pile based on the unified strength theory[J]. Rock and Soil Mechanics20062710): 1786. [百度学术

13

綦春明莫斌聂春龙. 饱和土体柱形扩孔时大变形不排水统一解析解[J]. 岩石力学与工程学报2009284): 827. [百度学术

QI ChunmingMO BinNIE Chunlonget al. Unified analytical solutions for cylindrical cavity expansion in saturated soil under large deformation and undrained conditions[J]. Chinese Journal of Rock Mechanics and Engineering2009284): 827. [百度学术

14

CAO L FTEH C ICHANG M F. Undrained cavity expansion in modified Cam clay I: theoretical analysis[J]. Géotechnique2001514): 323. [百度学术

15

肖昭然张昭杜明芳. 饱和土体小孔扩张问题的弹塑性解析解[J]. 岩土力学2004259): 1373. [百度学术

XIAO ZhaoranZHANG ZhaoDU Mingfang. An elastoplastic closed-form approach of cavity expansion in saturated soil based on modified Cam clay model[J]. Rock and Soil Mechianics2004259): 1373. [百度学术

16

胡伟刘明振. 非饱和土中球形孔扩张的弹塑性分析[J]. 岩土工程学报20062810): 1292. [百度学术

HU WeiLIU Mingzhen. Elastic-plastic solution of expansion of sphere cavity in unsaturated soils[J]. Chinese Journal of Geotechnical Engineering20062810): 1292. [百度学术

17

周茗如卢国文王腾. 结构性黄土劈裂注浆力学机理分析[J]. 工程力学2019363): 169. [百度学术

ZHOU MingruLU GuowenWANG Tenget al. Mechanism analysis of structured loess fracturing grouting[J]. Engineering Mechanics2019363): 169. [百度学术

18

宋勇军胡伟王德胜. 基于修正剑桥模型的挤密桩挤土效应分析[J]. 岩土力学2011323): 811. [百度学术

SONG YongjunHU WeiWANG Deshenget al. Analysis of squeezing effect of compaction piles based on modified Cam-clay model[J]. Rock and Soil Mechianics2011323): 811. [百度学术

19

李镜培唐剑华李林. 饱和黏土中柱孔三维弹塑性扩张机制研究[J]. 岩石力学与工程学报2016352): 378. [百度学术

LI JingpeiTANG JianhuaLI Linet al. Mechanism of three dimensional elastic-plastic expansion of cylindrical cavity in saturated clay[J]. Chinese Journal of Rock Mechanics and Engineering2016352): 378. [百度学术

20

李林李镜培孙德安. 剪胀性砂土中球孔扩张弹塑性解[J]. 岩土工程学报2017398): 1453. [百度学术

LI LinLI JingpeiSUN Deanet al. Elasto-plastic solution to expansion of a spherical cavity in dilatant sand[J]. Chinese Journal of Geotechnical Engineering2017398): 1453. [百度学术

21

CHEN S LABOUSLEIMAN Y N. Exact undrained elasto-plastic solution for cylindrical cavity expansion in modified Cam clay soil[J]. Géotechnique2012625): 447. [百度学术

22

CHEN S LABOUSLEIMAN Y N. Exact drained solution for cylindrical cavity expansion in modified Cam clay soil[J]. Géotechnique2013636): 510. [百度学术

23

郑金辉齐昌广王新泉.考虑砂土颗粒破碎的柱孔扩张问题弹塑性分析[J]. 岩土工程学报20194111): 2156. [百度学术

ZHENG JinhuiQI ChangguangWANG Xinquanet al. Elasto-plastic analysis of cylindrical cavity expansion considering particle breakage of sand[J]. Chinese Journal of Geotechnical Engineering20194111): 2156. [百度学术

24

邓国华邵生俊佘芳涛. 结构性黄土的修正剑桥模型[J]. 岩土工程学报2012345): 834. [百度学术

DENG GuohuaSHAO ShengjunSHE Fangtao. Modified Cam-clay model of structured loess[J]. Chinese Journal of Geotechnical Engineering2012345): 834. [百度学术

25

邓国华. 真三轴条件下黄土的结构性参数及结构性本构关系研究[D]. 西安西安理工大学2009. [百度学术

DENG Guohua. Research on structure parameter of loess and structure constitutive relations under true tri-axial condition[D]. Xi’anXi’an University of Technology2009. [百度学术

26

邓国华邵生俊. 基于真三轴试验的黄土结构性变化规律研究[J]. 岩土力学2013343): 679. [百度学术

DENG GuohuaSHAO Shengjun. Research on change structural characteristics of loess based on true triaxial tests[J]. Rock and Soil Mechanics2013343): 679. [百度学术

27

邵生俊周飞飞龙吉勇. 原状黄土结构性及其定量化参数研究[J]. 岩土工程学报2004264): 531. [百度学术

SHAO ShengjunZHOU FeifeiLONG Jiyong. Structural properties of loess and its quantitative parameter[J].Chinese Journal of Geotechnical Engineering2004264): 531. [百度学术