删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

K0固结天然饱和黏土中柱孔扩张弹塑性解

本站小编 哈尔滨工业大学/2019-10-24

K0固结天然饱和黏土中柱孔扩张弹塑性解

李林1,2,李镜培1,2,龚卫兵1,2,张凌翔1,2

(1. 岩土及地下工程教育部重点实验室(同济大学),上海 200092; 2. 同济大学 地下建筑与工程系,上海 200092)



摘要:

为研究天然饱和黏土中圆孔不排水扩张引起的土体力学响应,考虑天然饱和黏土的K0固结特性、应力历史及柱孔扩张过程中径向、环向和竖向应力之间的关系,基于修正剑桥模型和大变形理论,结合弹塑性边界条件推导了K0固结饱和天然黏土中柱孔不排水扩张问题的弹塑性解析解,并与既有数值解答进行对比验证.在此基础上,通过与以往基于初始应力各向同性假设下的柱孔扩张解答进行对比分析,研究了土体K0固结特性和应力历史对柱孔扩张过程中土体应力和塑性区半径的影响.研究结果表明:本文提出的柱孔扩张弹塑性解答与既有基于严格本构关系的数值解基本一致,且可以退化为基于初始应力各向同性假设条件下的解答;土体K0固结特性和应力历史对柱孔扩张过程中周围土体应力场和塑性区半径存在一定程度的影响,特别对极限扩张压力的影响较为明显,且其影响程度均随初始各向异性的增加而增大.研究结果为天然饱和黏土地基静压沉桩、旁压试验提供了一定的理论依据.

关键词:  初始应力各向异性  应力历史  修正剑桥模型  柱孔扩张  弹塑性解析解

DOI:10.11918/j.issn.0367-6234.201602038

分类号:TU473

文献标识码:A

基金项目:国家自然科学基金(41272288)



Elasto-plastic solution to expansion of a cylindrical cavity in K0-consolidated natural saturated clay

LI Lin1,2,LI Jingpei1,2,GONG Weibing1,2,ZHANG Lingxiang1,2

(1.Key Laboratory of Geotechnical and Underground Engineering (Tongji University), Ministry of Education, Shanghai 200092, China; 2.Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China)

Abstract:

To study the mechanical responses induced by undrained expansion of a cavity in natural clay soils, considering the K0-consolidated property and the stress history of the natural saturated clay, as well as the relationship among the radial, tangential, and vertical stresses, an elasto-plastic analytical solution to the undrained cylindrical cavity expansion in K0-consolidated saturated natural clay, which was under the initial stress anisotropic condition, was derived based on the modified Cam-clay (MCC) model, large-deformation theory and the elastic-plastic boundary conditions. The solution presented in this paper was verified by the well-documented numerical solution and compared with the previous initial stress isotropy assumption based solution. The effect of the K0-consolidated property and stress history on the stress fields and the plastic radius during the cavity expansion process was discussed in detail. The results show that the present closed-form elasto-plastic solution matches the numerical solution fairly well, and the present solution can be reduced to the initial stress isotropy assumption based solution. The K0-consolidated property has a significant effect on the stress fields around the cavity during the cavity expansion process, especially on the limit cavity expansion pressure. The effect of the initial stress anisotropy on cavity expansion increases with the degree of the initial stress anisotropy. This paper provides a theoretical benchmark for pile installation and pressuremeter test in natural clay soils.

Key words:  initial stress anisotropy  stress history  MCC model  cylindrical cavity expansion  elasto-plastic analytical solution


李林, 李镜培, 龚卫兵, 张凌翔. K0固结天然饱和黏土中柱孔扩张弹塑性解[J]. 哈尔滨工业大学学报, 2017, 49(6): 90-95. DOI: 10.11918/j.issn.0367-6234.201602038.
LI Lin, LI Jingpei, GONG Weibing, ZHANG Lingxiang. Elasto-plastic solution to expansion of a cylindrical cavity in K0-consolidated natural saturated clay[J]. Journal of Harbin Institute of Technology, 2017, 49(6): 90-95. DOI: 10.11918/j.issn.0367-6234.201602038.
基金项目 国家自然科学基金 (41272288) 作者简介 李林 (1986—),男,博士研究生;
李镜培 (1963—),男,教授,博士生导师 通信作者 李镜培,lijp2773@tongji.edu.cn 文章历史 收稿日期: 2016-02-26



Contents            -->Abstract            Full text            Figures/Tables            PDF


K0固结天然饱和黏土中柱孔扩张弹塑性解
李林1,2, 李镜培1,2, 龚卫兵1,2, 张凌翔1,2    
1. 岩土及地下工程教育部重点实验室 (同济大学),上海 200092;
2. 同济大学 地下建筑与工程系,上海 200092

收稿日期: 2016-02-26
基金项目: 国家自然科学基金 (41272288)
作者简介: 李林 (1986—),男,博士研究生;
李镜培 (1963—),男,教授,博士生导师
通信作者: 李镜培,lijp2773@tongji.edu.cn


摘要: 为研究天然饱和黏土中圆孔不排水扩张引起的土体力学响应,考虑天然饱和黏土的K0固结特性、应力历史及柱孔扩张过程中径向、环向和竖向应力之间的关系,基于修正剑桥模型和大变形理论,结合弹塑性边界条件推导了K0固结饱和天然黏土中柱孔不排水扩张问题的弹塑性解析解,并与既有数值解答进行对比验证.在此基础上,通过与以往基于初始应力各向同性假设下的柱孔扩张解答进行对比分析,研究了土体K0固结特性和应力历史对柱孔扩张过程中土体应力和塑性区半径的影响.研究结果表明:本文提出的柱孔扩张弹塑性解答与既有基于严格本构关系的数值解基本一致,且可以退化为基于初始应力各向同性假设条件下的解答;土体K0固结特性和应力历史对柱孔扩张过程中周围土体应力场和塑性区半径存在一定程度的影响,特别对极限扩张压力的影响较为明显,且其影响程度均随初始各向异性的增加而增大.研究结果为天然饱和黏土地基静压沉桩、旁压试验提供了一定的理论依据.
关键词: 初始应力各向异性    应力历史    修正剑桥模型    柱孔扩张    弹塑性解析解    
Elasto-plastic solution to expansion of a cylindrical cavity in K0-consolidated natural saturated clay
LI Lin1,2, LI Jingpei1,2, GONG Weibing1,2, ZHANG Lingxiang1,2    
1. Key Laboratory of Geotechnical and Underground Engineering (Tongji University), Ministry of Education, Shanghai 200092, China;
2. Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China


Abstract: To study the mechanical responses induced by undrained expansion of a cavity in natural clay soils, considering the K0-consolidated property and the stress history of the natural saturated clay, as well as the relationship among the radial, tangential, and vertical stresses, an elasto-plastic analytical solution to the undrained cylindrical cavity expansion in K0-consolidated saturated natural clay, which was under the initial stress anisotropic condition, was derived based on the modified Cam-clay (MCC) model, large-deformation theory and the elastic-plastic boundary conditions. The solution presented in this paper was verified by the well-documented numerical solution and compared with the previous initial stress isotropy assumption based solution. The effect of the K0-consolidated property and stress history on the stress fields and the plastic radius during the cavity expansion process was discussed in detail. The results show that the present closed-form elasto-plastic solution matches the numerical solution fairly well, and the present solution can be reduced to the initial stress isotropy assumption based solution. The K0-consolidated property has a significant effect on the stress fields around the cavity during the cavity expansion process, especially on the limit cavity expansion pressure. The effect of the initial stress anisotropy on cavity expansion increases with the degree of the initial stress anisotropy. This paper provides a theoretical benchmark for pile installation and pressuremeter test in natural clay soils.
Key words: initial stress anisotropy    stress history    MCC model    cylindrical cavity expansion    elasto-plastic analytical solution    
柱孔扩张理论由于形式简单,力学原理明确,因此在岩土工程领域得到广泛应用[1-4].然而,受理论及本构模型的限制,早期柱孔扩张解答多数假设土体为各向同性的理想弹塑性材料,且在解答过程中忽略竖向应力作用,仅考虑径向应力和环向应力来求解孔周土体的应力和位移[5-9].剑桥模型由于不仅可以综合考虑土体各应力之间的相互关系,同时可以体现土体剪切变形和剪切破坏的一致性,因此被广泛应用于柱孔弹塑性扩张研究中.如:Collins等[10-11]基于相似分析方法,采用原始剑桥模型和修正剑桥模型推导了初始应力各向同性土体中圆孔排水和不排水扩张问题的弹塑性解答;Yu[12]系统总结了初始应力各向同性土体中圆孔扩张问题的各类解答,并详细介绍了圆孔扩张理论在岩土工程中的应用;Cao等[13]、李镜培等[14]在初始应力各向同性假设的基础上,基于修正剑桥模型推导了圆孔扩张问题的半解析半数值解答;Chen等[15]基于修正剑桥模型,从严格的弹塑性本构关系出发,推导了K0固结土体中柱孔不排水扩张问题的数值解答,但其解答过程较为繁琐,不易推广和应用.可以看出,尽管目前基于剑桥模型的柱孔扩张解答考虑了柱孔扩张过程中各应力之间的相互关系,但其解答均没有考虑到K0固结特性对柱孔扩张过程中土体力学行为的影响,进而无法体现天然饱和黏土中柱孔弹塑性扩张的本质.

由于受沉积环境和固结环境的影响,天然状态下黏性土多处于K0固结状态.由于初始应力各向异性对土体的力学特性存在较大影响,因此实际工程中对应的圆孔扩张问题需考虑土体的K0固结特性.鉴于此,本文采用修正剑桥模型,在充分考虑土体K0固结特性和应力历史的基础上,推导了天然饱和黏土中柱孔不排水扩张问题的弹塑性解析解答,并在与Chen等[15]基于精确本构关系的数值解答对比验证的基础上,详细研究了土体初始应力各向异性对柱孔扩张过程的影响.

1 柱孔扩张力学模型如图 1所示,在内压力σa的作用下,初始半径为a0的柱孔在水平向应力和竖向应力分别为σh0σv0K0固结初始应力各向异性饱和黏土体中不排水扩张,随着内压力σa的不断增加,孔径从a0扩张至a,孔壁土体将首先发生屈服,由弹性状态进入塑性状态,进而在柱孔周围形成一定范围的塑性区.以au表示柱孔最终半径,以ry表示柱孔扩张过程中对应的塑性区半径,塑性区边界的径向位移以Ury表示.在整个求解过程中假定应力以压为正,且弹性区土体服从胡克定律和小变形假定,土体进入塑性后采用大变形理论和MCC模型来描述土体的力学特性.

Figure 1
图 1 柱孔扩张力学模型 Figure 1 Mechanical model for cylindrical cavity expansion


基于以上条件可知,对于初始应力各向异性天然饱和黏土,柱孔周围初始径向应力σr0、环向应力σθ0和竖向应力σz0之间存在如下关系:

${{\sigma '}_{{\rm{r0}}}} + {{\sigma '}_{{\rm{\theta 0}}}} = {K_0}{{\sigma '}_{{\rm{z}}0}},$ (1)

式中K0为土体静止侧压力系数,其大小反映了土体初始应力各向异性程度.

在柱孔不排水扩张的任何阶段,孔周土体中任一单元体都满足平衡方程:

$\frac{{\partial {\sigma _{\rm{r}}}}}{{\partial r}} + \frac{{{\sigma _{\rm{r}}} - {\sigma _{\rm{\theta }}}}}{r} = 0.$ (2)

2 柱孔扩张弹塑性解答 2.1 弹性区解答由弹性应力、应变和位移之间的相互关系,并考虑土体初始水平向应力σh0和初始竖向应力σv0之间的关系,可得弹性区位移解答[13]

${U_{\rm{r}}} = \frac{{\left( {{\sigma _{\rm{a}}} - {\sigma _{{\rm{h0}}}}} \right)}}{{2G}}{\left( {\frac{a}{r}} \right)^2}r,$ (3)

式中:G=3(1-2ν′)υp′/[2(1+ν′)κ]为土体的剪切模量, κ为卸荷曲线在υ-ln p′平面上的斜率, υ=1+e为土体比体积, e为土体孔隙比, p′=σ/3为平均有效主应力, ν'为土体有效泊松比.

进而,由胡克定律和小变形关系可得孔周弹性区应力分量σrσθσz的表达式[13]

${\sigma _{\rm{r}}} = {\sigma _{{\rm{h0}}}} + \left( {{\sigma _{\rm{a}}} - {\sigma _{{\rm{h0}}}}} \right){\left( {\frac{a}{r}} \right)^2},$ (4)

${\sigma _{\rm{\theta }}} = {\sigma _{{\rm{h0}}}} - \left( {{\sigma _{\rm{a}}} - {\sigma _{{\rm{h0}}}}} \right){\left( {\frac{a}{r}} \right)^2},$ (5)

${\sigma _{\rm{z}}} = {\sigma _{{\rm{v0}}}}.$ (6)

在柱孔不排水扩张条件下,土体体积不发生变化,因此,弹性区土体体应变可表示为

${\rm{d}}v = - \kappa \frac{{{\rm{d}}p'}}{{p'}} = 0.$ (7)

由式 (7) 可知,在扩孔过程中弹性区平均有效主应力p′保持不变,且由式 (4)~(6) 可知弹性区平均总主应力p亦保持不变,则根据有效应力原理可得弹性区超孔隙水压力为零,即:

$u = {u_0},$ (8)

式中u0为初始孔隙水压力.

2.2 弹塑性边界随着扩张内压力的不断增加,孔壁处土体首先屈服进入塑性状态.文献[15]根据屈服面之间的几何关系,给出了可以考虑土体K0固结特性的修正剑桥模型,其屈服准则可表示为[15]

$q = M\sqrt {p'\left[ {{{p'}_{\rm{c}}}\left( {1 + {{\left( {{\eta _0}/M} \right)}^2}} \right) - p'} \right]} .$ (9)

式中:η0=|3(1-K0)|/(2K0+1) 为初始应力比,pc为土体前期平均有效固结应力, Mp′-q平面中CSL线的斜率,其值等于6sin φ′/(3-sin φ′),φ′为土体有效内摩擦角, q为广义剪应力,可表示为:

$q = \sqrt {\frac{3}{2}{s_{{\rm{ij}}}}{s_{{\rm{ij}}}}} ,$ (10)

${s_{{\rm{ij}}}} = {\sigma _{{\rm{ij}}}} - p{\delta _{{\rm{ij}}}}.$ (11)

式中sij为偏应力张量,δij为克朗内克符号.

由于柱孔扩张过程中其周围弹性区平均有效主应力p′保持不变,故弹塑性边界处土体平均有效主应力py仍保持不变,即:

${{p'}_{\rm{y}}} = {{p'}_0}.$ (12)

由式 (12) 可知,柱孔弹性扩张阶段孔周土体单元在p′-q平面上有效应力路径从初始应力点 (p0, q0) 开始,竖直上升与初始屈服面相交于 (p0, qy) 点,由式 (9) 和式 (12) 可得弹塑性边界处的广义偏应力qy

${q_{\rm{y}}} = M{{p'}_0}\sqrt {{\rm{OCR}}\left( {1 + \frac{{\eta _0^2}}{{{M^2}}}} \right) - 1} .$ (13)

式中OCR=pc/p0,为土体超固结比.

在弹塑性边界处,弹性解答仍然成立,故由式 (4)~(6) 和式 (8) 可得弹塑性边界处应力分量σryσθyσzy之间存在如下关系:

${{\sigma }_{\text{ry}}}+{{\sigma }_{\text{ }\text{ }\!\!\theta\!\!\text{ }\text{ y}}}={{\sigma }_{\text{r0}}}+{{\sigma }_{\text{ }\theta \text{ 0}}}=2{{\sigma }_{\text{h0}}},$ (14)

${\sigma _{{\rm{zy}}}} = {\sigma _{{\rm{z0}}}}.$ (15)

由式 (10) 可知,弹塑性边界处的广义偏应力qy同样可表示为

${q_{\rm{y}}} = \sqrt {\frac{3}{2}\left( {{{\sigma '}_{{\rm{ijy}}}} - p'{\delta _{{\rm{ij}}}}} \right)\left( {{{\sigma '}_{{\rm{ijy}}}} - p'{\delta _{{\rm{ij}}}}} \right)} ,$ (16)

式中σijy为弹塑性边界处的应力张量.

联立式 (14)~(16) 可解得:

${{\sigma '}_{{\rm{ry}}}} = {{\sigma '}_{{\rm{r0}}}} + \frac{1}{{\sqrt 3 }}\sqrt {q_{\rm{y}}^2 - {{\left( {{\eta _0}{{p'}_0}} \right)}^2}} ,$ (17)

${{{\sigma }'}_{\text{ }\!\!\theta\!\!\text{ y}}}={{{\sigma }'}_{\text{r0}}}-\frac{1}{\sqrt{3}}\sqrt{q_{\text{y}}^{2}-{{\left( {{\eta }_{0}}{{{{p}'}}_{0}} \right)}^{2}}},$ (18)

${{\sigma '}_{{\rm{zy}}}} = {{\sigma '}_{{\rm{z0}}}}.$ (19)

将式 (17) 代入式 (3) 可得弹塑性边界处的位移:

${U_{{\rm{ry}}}} = \frac{{\sqrt {q_{\rm{y}}^2 - {{\left( {{\eta _0}{{p'}_0}} \right)}^2}} }}{{2\sqrt 3 G}}{r_{\rm{y}}}.$ (20)

由于半无限体中柱孔扩张为平面应变问题,故柱孔扩张过程中不产生竖向应变εz,即:

${\varepsilon _{\rm{z}}} = 0.$ (21)

在不排水条件下,体应变εv=0,故轴向应变εr和环向应变εθ之间存在如下关系:

${\varepsilon _{\rm{r}}} = - {\varepsilon _{\rm{\theta }}}.$ (22)

为考虑桩周土体进入塑性状态后的大变形效应,采用对数应变将柱孔周围塑性区内土体应变分量表示为:

${\varepsilon _{\rm{r}}} = - \ln \left( {\frac{{{\rm{d}}r}}{{{\rm{d}}{r_0}}}} \right),$ (23)

${{\varepsilon }_{\text{ }\!\!\theta\!\!\text{ }}}=-\ln \left( \frac{r}{{{r}_{0}}} \right).$ (24)

将式 (21) 和式 (22) 代入式 (24),在塑性区内积分,并考虑孔壁的初始条件和弹塑性边界条件可得:

$r_{\rm{y}}^2 - r_{{\rm{y0}}}^2 = {a^2} - a_0^2,$ (25)

式中ry0ry对应的初始半径位置.

考虑到在不排水条件下土体体积守恒,故弹塑性边界处位移Ury=ryry0,则式 (25) 可进一步表示为

$r_{\rm{y}}^2 - {\left( {{r_{\rm{y}}} - {U_{{\rm{ry}}}}} \right)^2} = {a^2} - a_0^2.$ (26)

将式 (20) 代入式 (26),并忽略式中$\sqrt{q_{\text{y}}^{\text{2}}-{{\left( {{\eta }_{0}}p_{0}^{'} \right)}^{2}}} $/G的高阶项,可得塑性区半径ry与当前孔径a及初始孔径a0之间的关系:

${\left( {\frac{{{r_{\rm{y}}}}}{a}} \right)^2} = \frac{{\sqrt 3 G}}{{\sqrt {q_{\rm{y}}^2 - {{\left( {{\eta _0}{{p'}_0}} \right)}^2}} }}\left[ {1 - {{\left( {\frac{{{a_0}}}{a}} \right)}^2}} \right].$ (27)

2.3 塑性区解答对于MCC模型,土体屈服后体应变增量dεv由弹性体应变增量dεve和塑性体应变dεvp两部分组成,且由于不排水条件下体积守恒,上述关系可表示为

$\kappa \frac{{{\rm{d}}p'}}{{\nu p'}} + \left( {\lambda - \kappa } \right)\frac{{{\rm{d}}{{p'}_c}}}{{\nu {{p'}_c}}} = 0,$ (28)

式中λ为各向同性加载曲线在υ-ln p′平面上的斜率.

对式 (28) 在塑性区内积分,并考虑到弹塑性边界处pc=pc0p′=p0,则可得:

${{p'}_{\rm{c}}} = p'{\rm{OCR}}{\left( {\frac{{p'}}{{{{p'}_0}}}} \right)^{\frac{{ - \lambda }}{{\lambda - \kappa }}}}.$ (29)

将式 (9) 代入式 (29) 可得:

$q = p'\sqrt {\left( {{M^2} + \eta _0^2} \right){\rm{OCR}}{{\left( {\frac{{p'}}{{{{p'}_0}}}} \right)}^{ - \frac{1}{\mathit{\Lambda }}}} - {M^2}} ,$ (30)

式中Λ=1-κ/λ为塑性体应变比.

由MCC模型对应的临界状态条件η=M,根据式 (30) 可求得孔周围土体单元进入临界状态时对应的平均有效主应力pf和偏应力qf分别为:

${{p'}_{\rm{f}}} = {{p'}_0}{\left[ {\frac{{\left( {{M^2} + \eta _0^2} \right){\rm{OCR}}}}{{2{M^2}}}} \right]^\mathit{\Lambda }},$ (31)

${q_{\rm{f}}} = M{{p'}_0}{\left[ {\frac{{\left( {{M^2} + \eta _0^2} \right){\rm{OCR}}}}{{2{M^2}}}} \right]^\mathit{\Lambda }}.$ (32)

由于平衡微分方程式 (2) 在塑性区仍成立,故柱孔扩张后孔周塑性区内任意位置处总应力σrx可通过积分式 (2) 得到:

${\sigma _{{\rm{rx}}}} = {\sigma _{{\rm{ry}}}} - \int_{{r_{\rm{y}}}}^{{r_{\rm{x}}}} {\frac{{\left( {{\sigma _{\rm{r}}} - {\sigma _{\rm{\theta }}}} \right)}}{r}} {\rm{d}}r,$ (33)

式中 (σrσθ) 与计算点的位置rx有关,因此,式 (33) 无法直接积分得出解析表达式.然而,在孔周塑性区内广义偏应力q基本为一定值,且孔周临界状态区内土体单元径向应力σrf,环向应力σθf与竖向应力σzf之间存在如下关系[15]

${{\sigma '}_{{\rm{zf}}}} = \frac{{{{\sigma '}_{{\rm{rf}}}} + {{\sigma '}_{{\rm{\theta f}}}}}}{2}.$ (34)

基于上述关系,如图 2所示,假设式 (34) 在整个塑性区内均成立,使 (σrσθ) 独立于计算点的位置rx,则塑性区半径ry可重新表示为

${\left( {\frac{{{r_{\rm{y}}}}}{a}} \right)^2} = \frac{{\sqrt 3 G}}{{\sqrt {q_{\rm{y}}^2 - {{\left( {{\eta _0}{{p'}_0}} \right)}^2}} }}\left[ {1 - {{\left( {\frac{{{a_0}}}{a}} \right)}^2}} \right].$ (35)

Figure 2
图 2 孔周假设偏应力分布与实际情况对比 Figure 2 Comparisons of derivative stress distributions around cavity between assumed and actual cases


同时,基于上述假设,式 (34) 可进一步表示为

${\sigma _{{\rm{rx}}}} = {\sigma _{{\rm{ry}}}} - \frac{{2{q_{\rm{f}}}}}{{\sqrt 3 }}\int_{{r_{\rm{y}}}}^{{r_{\rm{x}}}} {\frac{1}{r}{\rm{d}}r} .$ (36)

将式 (35) 代入式 (36),并在塑性区内积分可得孔周塑性区内径向总应力为

${\sigma _{{\rm{rx}}}} = {\sigma _{{\rm{h0}}}} + \frac{{\bar q}}{{\sqrt 3 }} + \frac{{{q_{\rm{f}}}}}{{\sqrt 3 }}\ln \frac{{\sqrt 3 G}}{{\bar q}}\left( {\frac{{{a^2} - a_0^2}}{{r_{\rm{x}}^2}}} \right),$ (37)

式中q=$\sqrt{q_{\text{f}}^{\text{2}}-{{\left( {{\eta }_{0}}p_{0}^{'} \right)}^{2}}}$.

将式 (34) 和式 (37) 代人式 (32),可得塑性区环向总应力和竖向总应力的解析式:

${{\sigma }_{\text{ }\!\!\theta\!\!\text{ x}}}={{\sigma }_{\text{h0}}}+\frac{{\bar{q}}}{\sqrt{3}}+\frac{{{q}_{\text{f}}}}{\sqrt{3}}\left[ \ln \frac{\sqrt{3}G}{{\bar{q}}}\left( \frac{{{a}^{2}}-a_{0}^{2}}{r_{\text{x}}^{2}} \right)-2 \right],$ (38)

${\sigma _{{\rm{zx}}}} = {\sigma _{{\rm{h0}}}} + \frac{{\bar q}}{{\sqrt 3 }} + \frac{{{q_{\rm{f}}}}}{{\sqrt 3 }}\left[ {\ln \frac{{\sqrt 3 G}}{{\bar q}}\left( {\frac{{{a^2} - a_0^2}}{{r_{\rm{x}}^2}}} \right) - 1} \right].$ (39)

根据有效应力原理,结合式 (37)~(39) 和式 (31),可得塑性区内任意位置处的超孔隙水压力Δurx

$\begin{array}{l}\Delta {u_{{\rm{rx}}}} = \frac{{\left( {{K_0} - 1} \right){{p'}_0}}}{{1 + 2{K_0}}} + \frac{{{q_{\rm{f}}}}}{{\sqrt 3 }}\ln \frac{{\sqrt 3 G}}{{\bar q}}\left( {\frac{{{a^2} - a_0^2}}{{r_{\rm{x}}^2}}} \right) + \\\;\;\;\;\;\;\;\;\frac{{\bar q - {q_{\rm{f}}}}}{{\sqrt 3 }} + {{p'}_0}\left[ {1 - {{\left[ {\frac{{\left( {{M^2} + \eta _0^2} \right){\rm{OCR}}}}{{2{M^2}}}} \right]}^\mathit{\Lambda }}} \right].\end{array}$ (40)

令式 (37) 和式 (40) 中rx=a,可得孔壁处的扩孔压力σa和孔壁处的超孔隙水压力Δua

${\sigma _{\rm{a}}} = {\sigma _{{\rm{h0}}}} + \frac{{\bar q}}{{\sqrt 3 }} + \frac{{{q_{\rm{f}}}}}{{\sqrt 3 }}\ln \frac{{\sqrt 3 G}}{{\bar q}}\left( {1 - \frac{{a_0^2}}{{{a^2}}}} \right),$ (41)

$\begin{array}{l}\Delta {u_{\rm{a}}} = \frac{{\left( {{K_0} - 1} \right){{p'}_0}}}{{1 + 2{K_0}}} + \frac{{{q_{\rm{f}}}}}{{\sqrt 3 }}\ln \frac{{\sqrt 3 G}}{{\bar q}}\left( {1 - \frac{{a_0^2}}{{{a^2}}}} \right) + \\\;\;\;\;\;\;\;\;\frac{{\bar q - {q_{\rm{f}}}}}{{\sqrt 3 }} + {{p'}_0}\left[ {1 - {{\left[ {\frac{{\left( {{M^2} + \eta _0^2} \right){\rm{OCR}}}}{{2{M^2}}}} \right]}^\mathit{\Lambda }}} \right].\end{array}$ (42)

a/a0→∞时,从式 (41) 和式 (42) 可得孔壁极限扩孔压力σrult和孔壁处的极限超孔隙水压力Δuault

$\begin{array}{l}\Delta {u_{{\rm{ault}}}}:\\\;\;\;\;\;\;\;\;\;{\sigma _{{\rm{rult}}}} = {\sigma _{{\rm{h0}}}} + \frac{{\bar q}}{{\sqrt 3 }} + \frac{{{q_{\rm{f}}}}}{{\sqrt 3 }}\left( {\ln \frac{{\sqrt 3 G}}{{\bar q}} + 1} \right),\end{array}$ (43)

$\begin{array}{l}\Delta {u_{{\rm{ault}}}} = \frac{{\left( {{K_0} - 1} \right){{p'}_0}}}{{1 + 2{K_0}}} + \frac{{{q_{\rm{f}}}}}{{\sqrt 3 }}\ln \frac{{\sqrt 3 G}}{{\bar q}} + \frac{{\bar q - {q_{\rm{f}}}}}{{\sqrt 3 }} + \\\;\;\;\;\;\;\;\;\;\;{{p'}_0}\left[ {1 - {{\left[ {\frac{{\left( {{M^2} + \eta _0^2} \right){\rm{OCR}}}}{{2{M^2}}}} \right]}^\mathit{\Lambda }}} \right].\end{array}$ (44)

值得注意的是:若令式 (37)~(44) 中K0=1,则土体为各向同性体,本文解答可退化为各向同性体中柱孔不排水扩张解答,与文献[13]所得出的各向同性土体中柱孔扩张解答完全一致,说明文献[13]的解答是本文解答的一个特例.

3 验证与分析为验证本文所得解答的正确性和合理性,同时为分析K0固结特性和应力历史对柱孔扩张过程的影响,本文选取不同K0固结条件下的3种饱和黏土进行分析,其土体参数取值均引自文献[15],具体数值见表 1.

表 1
表 1 土体参数取值[15] Table 1 Properties of soils[15] OCRK0σv0/kPaG0/kPaυ0u0/kPa

1.20.6251604 3022.06100

3.01.0001204 1131.97100

10.02.000723 7561.80100



表 1 土体参数取值[15] Table 1 Properties of soils[15]


此外,根据文献[15],4种土体临界状态线CSL斜率M=1.2,初始K0固结曲线斜率λ=0.15,回弹曲线的斜率κ=0.03,泊松比ν=0.278.

3.1 对比验证图 3(a)~(c)为本文解析解与文献[15]所得不同初始应力各向异性土体中柱孔不排水扩张后孔周土体应力分布的对比,图中横坐标轴均为对数坐标.

Figure 3
图 3 孔壁周围应力分布解析解与数值解对比 Figure 3 Comparisons of stress distributions around cavity between numerical solution and analytical solution


从图 3可见,本文解析解答与基于严格本构关系的数值解在靠近孔壁的塑性区内吻合良好,仅在远离孔壁的弹塑性边界处两者存在一些差异.然而,对于静压沉桩和静力触探等岩土工程问题,多关注靠近孔壁处的应力分布,因此本文解答可以满足一般工程问题的需要.

同时,从图 3可见,初始应力各向异性饱和黏土中柱孔不排水扩张引起孔壁周围土体应力状态发生较大变化,且随距孔壁距离的增加,土体应力变化逐渐减小,在弹性区内趋于初始应力状态.

3.2 柱孔扩张过程分析图 4~6分别给出了初始应力各向异性黏土中柱孔不排水扩张过程中孔壁压力σa、孔壁超孔压Δua及塑性区半径ry随当前孔径a的变化规律,图中变量均采用无量纲化形式表示.

Figure 4
图 4 柱孔扩张过程中孔壁应力变化 Figure 4 Variations of cavity stresses during cavity expansion


Figure 5
图 5 柱孔扩张过程中孔壁超孔隙水压力变化 Figure 5 Variations of excess pore water pressures at cavity wall during cavity expansion


Figure 6
图 6 柱孔扩张过程中塑性区半径变化 Figure 6 Variations of plastic radius during cavity expansion


从图 4可见,孔壁压力σa在柱孔扩张初期 (a/a0 < 2) 迅速增大,之后 (a/a0≥2) 孔壁压力σa增长缓慢且逐步趋于临界状态,当a/a0≥5时,孔壁压力趋于稳定,土体进入临界状态.

从图 5可知,与孔壁压力σa相对应,孔壁处的超孔压在柱孔扩张初期随孔径的增大而迅速增大,此后逐步趋向于恒定值.值得注意的是,对于超固结土 (OCR>1), 在柱孔扩张初期孔壁处首先产生负超孔压,之后随孔径的增大,超孔压由负值转变为正值,且逐步趋向于临界值,这一现象与文献[13, 15]分析所得规律一致.

从图 6可看出,塑性区半径ry在柱孔扩张过程中随孔径a的变化规律与孔壁压力σa和孔壁超孔隙水压力Δua随孔径的变化规律一致,均是在扩孔初期变化较大,之后逐步趋于定值.

除以上规律外,图 4~6亦给出了基于各向同性假设的解答[13]和本文考虑初始应力各向异性解答之间的差异.可以看出,初始应力各向异性和应力历史对柱孔扩张过程有明显的影响,特别是对柱孔扩张过程中极限扩张压力的影响最为明显.

4 结论本文针对K0固结初始应力各向异性土体中柱孔不排水扩张问题,采用修正剑桥模型在考虑土体K0固结特性和应力历史的基础上,推导了柱孔不排水扩张的弹塑性解析解.在与既有基于严格本构关系的数值解对比验证的基础上,分析了K0固结特性及应力历史对柱孔扩张的影响,得出如下结论:

1)K0固结黏土中柱孔不排水扩张引起孔壁周围土体应力状态发生较大变化,且孔周应力随距径向距离的增加呈对数规律衰减,在弹性区逐渐趋向于初始应力状态.

2) 孔壁压力、超孔隙水压力和塑性区半径在柱孔不排水扩张的初始阶段变化较大,之后随孔径的增大而逐步趋向于临界状态.

3) 土体初始应力各向异性对孔壁压力、超孔隙水压力及塑性区半径存在一定程度的影响,特别对极限扩张压力存在较大影响,且其影响程度随土体各向异性的增加而增大.


参考文献
[1]VESIC A S. Expansion of cavity in infinite soil mass[J].Journal of Soil Mechanics Foundation Division, American Society of Civil Engineering, 1972, 98(3): 265-289.

[2]RANDOLPH M F, CARTER J P, WROTH C P. Driven piles in clay-the effects of installation and subsequent consolidation[J].Géotechnique, 1979, 29(4): 361-393.DOI: 10.1680/geot.1979.29.4.361

[3]SAGASET A C. Analysis of undrained soil deformation due to ground loss[J].Géotechnique, 1987, 37(3): 301-320.DOI: 10.1680/geot.1987.37.3.301

[4]SAGASET A C. Prediction of ground movements due to pile-driving in clay[J].Journal of Geotechnical and Geoenvironmental Engineering, 2001, 127(1): 55-66.DOI: 10.1061/(ASCE)1090-0241(2001)127:1(55)

[5]RANDOLPH M F. Science and empiricism in pile foundation design[J].Géotechnique, 2003, 53(10): 847-875.DOI: 10.1680/geot.2003.53.10.847

[6]HILL R. The mathematical theory of plasticity[M]. London: Oxford University Press, 1950.

[7]CHADWICK P. The quasi-static expansion of spherical cavity in metals and ideal soils[J].Quarterly Journal of Mechanics and Applied Mathematics, 1959, 12: 52-71.DOI: 10.1093/qjmam/12.1.52

[8]CATER J P, RANDOLPH M F, WROTH C P. Stress and pore pressure changes in clay during and after the expansion of a cylindrical cavity[J].International Journal for Numerical and Analytical Methods in Geomechanics, 1979, 3(4): 305-322.DOI: 10.1002/(ISSN)1096-9853

[9]YU H S, HOULSBY G T. Finite cavity expansion in dilatants soils: loading analysis[J].Géotechnique, 1991, 41(2): 173-183.DOI: 10.1680/geot.1991.41.2.173

[10]COLLINS I F, STIMPSON J R. Similarity solutions for drained and undrained cavity expansions in soils[J].Géotechnique, 1994, 44(1): 21-34.DOI: 10.1680/geot.1994.44.1.21

[11]COLLINS I F, YU H S. Undrained cavity expansions in critical state soils[J].International Journal for Numerical and Analytical Methods in Geomechanics, 1996, 20(7): 489-516.DOI: 10.1002/(ISSN)1096-9853

[12]YU H S. Cavity expansion methods in geomechanics[M]. Kluwer Academic Publishers: 2000.

[13]CAO L F, TEH C I, CHANG M F. Undrained cavity expansion in modified Cam Clay Ⅰ: theoretical analysis[J].Géotechnique, 2001, 51(4): 323-334.DOI: 10.1680/geot.2001.51.4.323

[14]李镜培, 唐剑华, 张亚国, 等. 饱和粘土中球孔扩张问题弹塑性解析[J].哈尔滨工业大学学报, 2014, 46(12): 71-77.
LI Jingpei, TANG Jianhua, ZHANG Yaguo, et al. elastic-plastic solution of sphere cavity expansion in saturated clay[J].Journal of Harbin institute of technology, 2014, 46(12): 71-77.DOI: 10.11918/j.issn.0367-6234.2014.12.012

[15]CHEN S L, ABOUSLEIMAN Y N. Exact undrained elasto-plastic solution for cylindrical cavity expansion in modified Cam Clay soil[J].Géotechnique, 2012, 62(5): 447-456.DOI: 10.1680/geot.11.P.027


相关话题/塑性 历史 同济大学 上海 过程

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 热轧过程中高纯钴微观组织及织构演变
    热轧过程中高纯钴微观组织及织构演变李震1,宋克睿1,韩彦鹏1,肖柱1,贺昕2,陈志永1(1.中南大学材料科学与工程学院,长沙410083;2.北京有色金属研究总院有研亿金新材料股份有限公司,北京102200)摘要:金属钴具有同素异构转变特性。为探究热轧工艺对高纯钴的微观组织及织构演变规律的影响,对纯 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • AZ31B镁合金带材热轧过程组织均匀性及性能研究
    AZ31B镁合金带材热轧过程组织均匀性及性能研究曹东东1,2,梅瑞斌1,2,包立1,侯铮1,黄芸1(1.东北大学秦皇岛分校资源与材料学院,河北秦皇岛066004;2.东北大学材料科学与工程学院,沈阳110819)摘要:本文开展了变形温度为300、350、400℃和总压下率分别为15%、30%、45% ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 考虑高温水泄漏影响的凝水泵启动过程仿真分析
    考虑高温水泄漏影响的凝水泵启动过程仿真分析覃海波,金家善,倪何(海军工程大学动力工程学院,武汉430033)摘要:针对某船用给水机组凝水泵在备用状态和备用转换启动过程中,由于高温除氧水进入凝水泵及其吸入管路而在机组启动后发生汽蚀,导致凝水泵出口压力长时间达不到规定要求的问题,在基于差异演化算法辨识得 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 污水生物处理过程中溶解性有机氮分布和转化特征
    污水生物处理过程中溶解性有机氮分布和转化特征王小东1,陈明飞1,王子文1,王燕1,王硕1,2,3,李激1,2,3(1.江南大学环境与土木工程学院,江苏无锡214122;2.江苏省厌氧生物技术重点实验室(江南大学),江苏无锡214122;3.江苏省高校水处理技术与材料协同创新中心,江苏苏州215009 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 卷积神经网络多变量过程特征学习与故障诊断
    卷积神经网络多变量过程特征学习与故障诊断陈淑梅,余建波(同济大学机械与能源工程学院,上海201804)摘要:为提取复杂多变量过程的有效特征,提高故障诊断性能,提出一种基于卷积神经网络(convolutionalneuralnetwork,CNN)特征学习的多变量过程故障诊断模型.将高维过程信号归一化 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 具有荷载作用历史的CFRP-钢界面粘结性能试验
    具有荷载作用历史的CFRP-钢界面粘结性能试验徐佰顺1,姚亚东1,钱永久2,马明3,宋帅4(1.内蒙古大学交通学院,呼和浩特010070;2.西南交通大学土木工程学院,成都610031;3.中交第一公路勘察设计研究院有限公司,西安710075;4.太原理工大学建筑与土木工程学院,太原030024)摘 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 钠质膨润土渗透性与膨胀性及可塑性的相关性
    钠质膨润土渗透性与膨胀性及可塑性的相关性于泽溪1,2,李育超1,2,陈冠年1,2(1.软弱土与环境土工教育部重点实验室(浙江大学),杭州310058;2.浙江大学岩土工程研究所,杭州310058)摘要:盐溶液与钠质膨润土的化学相容性研究对评估防污屏障的服役性能起到重要作用,为此,采用渗透试验、自由膨 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • IN718单晶微柱体锯齿状变形的晶体塑性分析
    IN718单晶微柱体锯齿状变形的晶体塑性分析张旭1,3,祝明昊1,张万瑞2,尚福林3(1.郑州大学力学与工程科学学院,郑州450001;2.西北大学机械工程系,埃文斯顿60201,美国;3.西安交通大学航天航空学院,西安710049)摘要:为较为准确可靠地预测IN718合金服役于各种复杂工况下的力学 ...
    本站小编 哈尔滨工业大学 2019-10-24
  • 典型台风登陆过程平均风时距转换系数分析
    典型台风登陆过程平均风时距转换系数分析潘晶晶1,赵林1,梁旭东2,朱乐东1,葛耀君1(1.土木工程防灾国家重点实验室(同济大学),上海200092;2.中国气象局北京城市气象研究所,北京100089)摘要:为掌握中国东南沿海台风登陆过程近地面风速变化规律及影响因素,根据2005年浙江省东海塘观测塔和 ...
    本站小编 哈尔滨工业大学 2019-10-24
  • 臭氧氧化过程中羟基自由基产率测定与分析
    臭氧氧化过程中羟基自由基产率测定与分析刘永泽,江进,马军,罗从伟,皇甫小留,郭忠凯(城市水资源与水环境国家重点实验室(哈尔滨工业大学),150090哈尔滨)摘要:为研究臭氧氧化过程中OH·的产生,利用“叔丁醇方法”测定了芳香类、胺类、烯烃类等活性物质与O3反应过程中OH·的产率,并进一步测定了不同水 ...
    本站小编 哈尔滨工业大学 2019-10-24