DOI:10.19476/j.ysxb.1004.0609.2002.01.032
饱和粘土中平均塑性体积应变对压密注浆的影响
温世游 陈云敏 李夕兵 陈仁朋
浙江大学岩土工程研究所
中南大学资源环境与建筑工程学院
浙江大学岩土工程研究所 杭州310027
长沙410083
摘 要:
针对工程上通常不考虑平均塑性体积应变Δ的影响或者简单取Δ为 0 .0 1的情况 , 分析了饱和粘土中平均塑性体积应变对压密注浆的影响范围、极限孔压、超孔隙水压力的影响并给出考虑其影响的计算公式。研究表明 :在工程中平均塑性体积应变的影响通常是应该考虑的 , 它是反映土体变形规律的重要参数 , 并提出了平均塑性体积应变的简化确定方法 , 操作简单 , 判断快捷 , 适合工程应用。当G/Cu≤ 4 0时 , 可以不考虑平均塑性体积应变的影响。
关键词:
饱和粘土 ;平均塑性体积应变 ;压密注浆 ;剪切刚度比 ;
中图分类号: TU472
收稿日期: 2001-04-28
基金: 浙江省自然科学基金青年人才专项基金资助项目 (RC960 9);
Effects of average plastic volumetric strain on compaction grouting in saturated clay
Abstract:
Aiming at the situation of neglecting the effect of average plastic volumetric strain or giving it the value 0.01 in engineering, the author analyzed its effects on the effect scope, displacement in plastic zone, ultimate pressure, excess pore water pressure in compaction grouting in saturated clay and presented the corresponding formulas. It is shown that the effect of average plastic volumetric strain should be considered generally in engineering. It's an important parameter reflecting the deformation law of soil. The average plastic volumetric strain can be obtained by the inversion with the heave height on the surface. In general, when G/C u is less than 40, its effect can be neglected.
Keyword:
saturated clay; average plastic volumetric strain; compaction grouting; shear rigidity ratio;
Received: 2001-04-28
在地质工程、 岩土工程中, 压密注浆技术日益被大量用来处理饱和软土地基。 通常对平均塑性体积应变Δ的影响不予考虑或者简单取Δ=0.01作近似处理。 当扩张的塑性区半径与孔 (桩) 径之比大于10时, 这种近似处理有一定的合理性和工程精度。 但实际工程中, 无法事先知道塑性区的半径, 而且Δ的变化范围在10-3 ~10-1 之间。 工程上认为当Δ大于0.01时土体的结构性受到破坏。 不考虑平均塑性体积应变的影响主要是因为难以确定Δ。 Vesic提出用迭代法数值求解
[1 ]
或者实验室测定土体参数来确定
[1 ,2 ]
。 Skempton
[3 ]
曾尝试采用数值法来评价其重要性。 文献
[
4 ,
5 ]
认为在饱和粘土压密注浆的过程中, 其影响范围只与被加固土体的刚度比和注浆孔径有关。 利用实验室测定方法的致命弱点在于可能扰动原状土, 而一旦Δ大于0.01时, 实测值与真实值的吻合性值得怀疑。 根据塑性流动法由Carter
[6 ]
可以推导出Δ的计算公式, 但对于饱和粘土来说, Δ=0.5, 无法计算式中的D E , D -1 E , 也就是说理论上可行, 实际上无法操作。 研究中我们发现Δ是反映土体变形的重要参数, 在大多数的工程中是不可忽略的, 必须加以考虑, 这与文献
[
1 ,
7 ]
的研究结果是一致的。 本文作者分析了饱和粘土中平均塑性体积应变对压密注浆的影响范围、 塑性区位移、 极限孔压、 超孔隙水压力的影响以及确定平均塑性体积应变的简单实用方法。
1饱和粘土中平均塑性体积应变对压密注浆的影响
在分析中作如下假定:土体是均质各向同性的理想弹塑性材料, 服从摩尔-库仑屈服准则, 不计体积力的影响。 为研究方便, 假设注浆体在土中呈圆柱形, 压密注浆过程就相当于在无限土体中扩张一个半径为r 的小孔, 在孔穴四周形成一个应力影响区。 这样压密注浆过程就模拟为半无限土体中的圆柱孔的扩张, 当压力增加时, 围绕着圆柱的区域将由弹性状态进入塑性状态。 随着内压力的逐渐增大, 塑性区不断扩大, 直至内压力增大到极限孔压力p u , 此时圆柱孔半径为R u 。 设弹、 塑性交界处的半径为R , 在半径R 以内为塑性区, R 以外土体仍处于弹性状态, 如图1所示。
图1 圆孔的扩张
Fig.1 Cylindrical cavity expansion
1.1 压密注浆影响范围的确定
圆柱孔扩张按平面轴对称问题来研究, 其平衡微分方程为
d
σ
r
d
r
+
σ
r
-
σ
t
r
=
0
?
?
?
(
1
)
式中 σ r , σ t 分别为径向应力、 切向应力。
在弹塑性边界上必须满足摩尔-库仑破坏准则, 对饱和粘土 (φ =0) 在不排水条件下, 有
σ r -σ t =2C u (2)
根据式 (1) , (2) 可求得塑性区内任意一点在3个方向的应力分量为
σ
r
=
[
2
ln
(
R
/
r
)
+
1
]
C
u
σ
t
=
[
2
ln
(
R
/
r
)
-
1
]
C
u
σ
z
=
2
ln
(
R
/
r
)
C
u
}
?
?
?
(
3
)
式中 C u 为饱和粘土的不排水抗剪强度, R 为弹塑性交界处半径, σ z 为垂直向应力。
同理可得出塑性区交界处位移为
u
p
=
1
+
ν
E
R
C
u
?
?
?
(
4
)
考虑到周围土体变形后的体积变化等于弹性区的体积变化和塑性区的体积变化之和, 可得
πR
u
2
=π R2 -π (R-up ) 2 +π (R2 -R
u
2
) Δ (5)
式中 Δ 为平均塑性体积应变。 联立式 (4) 和 (5) , 略去高阶微量, 得塑性区半径为
R
=
1
+
Δ
C
u
G
+
Δ
R
u
?
?
?
(
6
)
式中 G为剪切模量, G=E/2 (1+ν) , ν为泊松比, G/Cu 称为剪切刚度比。 当不考虑平均塑性体积应变的影响时, 式 (6) 退化为文献
[
4 ]
中的式 (7) 。 可见在饱和粘土压密注浆的过程中, 其影响半径不仅与被加固土体的剪切刚度比和注浆孔径有关, 而且与平均塑性体积应变有关。
又根据文献
[
4 ]
可得塑性区的位移满足:
u
r
=
ε
R
[
A
(
R
r
)
1
+
α
+
B
(
R
r
)
1
-
β
+
C
]
?
?
?
(
7
)
式中 α, β, A, B和C的含义见文献
[
6 ]
。 对于饱和粘土, A=1, B=0, C=0, α=1, β=1。 由文献
[
5 ]
得:
ε
R
=
C
u
2
G
?
?
?
(
8
)
综合式 (7) 和 (8) 可得:
u
=
C
u
2
G
R
2
r
=
R
u
2
2
r
?
?
?
(
9
)
可见塑性区位移随着距离的增大而减小。 考虑平均塑性体积应变的影响时,
u
=
C
u
2
G
1
+
Δ
C
u
G
+
Δ
R
u
2
r
=
1
+
Δ
1
+
Δ
G
/
C
u
R
u
2
2
r
?
(
R
u
≤
r
≤
R
)
?
?
?
(
1
0
)
可见实际的塑性区位移比不考虑平均塑性体积应变的影响时要小。
1.2 平均塑性体积应变对塑性区半径的影响
从式 (6) 发现塑性区半径R 与土的剪切刚度比、 注浆孔径和平均塑性体积应变有关。 下面我们分别对不同数量级的剪切刚度比来研究平均塑性体积应变对压密注浆塑性区半径的影响。 表1所示为G /C u =1 000, 100, 50, 40时, 考虑平均塑性体积应变与不考虑平均塑性体积应变对压密注浆塑性区半径的影响。 图2所示为不同刚度比下R /R u -Δ的关系。
表1 平均塑性体积应变对压密注浆塑性区半径的影响
Table 1 Effects of average plastic volumetric strain on radius R at different shear rigidity ratio in compaction grouting
Δ
R /R u
G /C u =1 000
G /C u =100
G /C u =50
G /C u =40
0 0.001 0.002 0.003 0.004 0.005 0.007 0.01 0.02
31.62 22.37 18.28 15.84 14.17 12.94 11.22 9.58 6.97
10 9.58 9.14 8.78 8.45 8.19 7.7 7.11 5.83
7.1 6.9 6.75 6.6 6.47 6.34 6.11 5.8 5.0
6.3 6.2 6.1 6.0 5.9 5.8 5.6 5.4 4.8
图2 不同刚度比下R/Ru-Δ关系图
Fig.2 Relationship between R /R u and Δat different shear rigidity ratio
从表1和图2中可见:随着剪切刚度比G /C u 的增大, 平均塑性体积应变对压密注浆塑性区半径的影响越显著。 在同一剪切刚度比下, 随着平均塑性体积应变的增大, 其对压密注浆塑性区半径的影响越显著。 当剪切刚度比G /C u 的数量级在103 时, 平均塑性体积应变的影响必须加以考虑。 当G /C u ≤40时, 一般可以不考虑平均塑性体积应变的影响, 如文献
[
8 ]
中的桩D 1 , D 2 , D 3 。 以上海地区软土为例, 淤泥质粘土层的剪切刚度比为50~130, 所以应该考虑平均塑性体积应变的影响。
1.3 平均塑性体积应变的简化确定方法
饱和粘土中, 当灌入土中的注浆深度为H 时, 排开的土体体积为πR 2 u H 。 由于周围土体受到挤压作用产生弹塑性变形, 设塑性体积应变为ε V 1 , 弹性体积应变为ε V 2 , 总的体积应变为ε V =ε V 1 +ε V 2 。 由于塑性范围相对R u 来说, 径向挤土效应扩散到弹性区已经变得很小, 所产生的弹性应变相对于塑性应变来说, 可以忽略不计, ε V 2 ?ε V 1 , 所以近似地认为总的体积应变为ε V =ε V 1 =Δ。
注浆时土体受到径向挤压后, 使得土体向地表面挤出, 土体的侧向挤出主要是由于土体塑性流动产生的, 因此可近似地认为注浆区周围土体由于受到挤压而产生的隆起发生在塑性区的范围内, 得:
πR 2 u H (1-Δ) =πR 2 w (11)
定义m =w /H 为隆起比, 则由式 (6) 可得:
w
=
1
-
Δ
1
+
Δ
(
C
u
G
+
Δ
)
Η
?
?
?
(
1
2
)
可见对一定的土体来说, w 与H , Δ, G /C u 有关。 上式变化为
Δ2 + (m -1+C u /G ) Δ+ (m -C u /G ) =0 (13)
这样可以通过测定地表土体的隆起高度来反演确定平均塑性体积应变。 这种方法简单快捷, 不必考虑土体的结构性是否扰动, 因为测定的地表土体隆起高度能够间接反映土体的变形规律, 土体变形是一个有机统一的自平衡系统, 通过这种途径反演平均塑性体积应变能够满足工程需要。 也可利用文献
[
4 ]
式 (23) 来反算平均塑性体积应变Δ:
w
max
=
2
C
u
ln
(
R
/
R
u
)
Η
E
-
q
Η
E
+
2
ν
(
C
u
E
+
Κ
0
q
E
)
Η
?
?
?
(
1
4
)
式中 K 0 地基土体的静止侧压系数, q 为地面荷载, 当q =0时, 对饱和粘土来说, 可以简化为
w
max
=
1
3
[
ln
(
1
+
Δ
C
u
G
+
Δ
)
+
1
]
C
u
G
Η
?
?
?
(
1
5
)
同样说明地面隆起高度w 是注浆深度H 、 平均塑性体积应变Δ、 剪切刚度比G /C u 的函数。
1.4平均塑性体积应变对极限孔压、 超孔隙水压力影响
在饱和粘土中进行压密注浆挤土效应十分明显, 引起的超孔隙水压力由于土体的渗透系数小, 因而消散慢, 对工程建设和周围环境会带来较严重的不利影响。 因此研究平均塑性体积应变对极限孔压、 超孔隙水压力影响尤为重要。 考虑初始应力p 0 , 饱和土体中, 极限孔压p u 的计算公式为:
p
u
=
p
0
+
[
1
+
2
ln
(
R
R
u
)
]
C
u
=
p
0
+
[
1
+
2
ln
1
+
Δ
C
u
G
+
Δ
)
]
C
u
?
?
?
(
1
6
)
由文献
[
1 ]
得超孔隙水压力的计算公式为
Δ
p
C
u
2
ln
R
r
+
1
.
7
3
A
f
-
0
.
5
8
?
?
?
(
1
7
)
式中 Δp 为孔隙水压力增量, A f 为Skempton孔隙压力系数
[9 ]
。 可见离注浆孔中心距离越近, 则超孔隙水压力上升越高, 最大超孔隙水压力发生在注浆孔边界, 其值为2ln (R /R u ) +1.73A f -0.58。 随着距离的增大, 孔压是按对数曲线变化, 直至减小到0。 考虑平均塑性体积应变的影响时,
Δ
p
Δ
C
u
=
2
ln
R
Δ
r
+
1
.
7
3
A
f
-
0
.
5
8
?
?
?
(
1
8
)
式中 Δp Δ , R Δ 分别为考虑平均塑性体积应变的影响时超孔隙水压力和塑性区半径。 两者压力差为
Δ
p
-
Δ
p
Δ
=
C
u
ln
(
1
+
G
C
u
Δ
1
+
Δ
)
?
?
?
(
1
9
)
2 算例
2.1 算例1
取文献
[
4 ]
中的工程为例。 注浆深度为5.0 m, 注浆量为0.78 m3 /孔, C u =150 kPa, E =5 MPa, K 0 =0.75, ν =0.5, 现场地面无超载。 根据式 (14) 算出w max =8.5 cm, 现场实测隆起为7.6 cm, 误差为11.8%。 根据m =w /H =0.015 2, 代入式 (13) 算得平均塑性体积应变Δ=0.007, 应该考虑平均塑性体积应变的影响。 当考虑Δ时, 由 (15) 计算出地面隆起高度为7.7 cm, 误差为1.2%。
2.2 算例2
有一正常固结粘土, 初始应力p 0 =220 kPa, R u =250 cm, C u =60 kPa, E =12 MPa, ν =0.5, 实测到极限孔压为p u =0.52 MPa。 考虑平均塑性体积应变的影响, 取Δ=0.005时, 计算出极限孔压为0.515 MPa, Δp Δ =255 kPa; 不考虑平均塑性体积应变的影响, 计算出极限孔压为0.53 MPa, Δp =290 kPa。 可见在工程中考虑其影响是很有必要的。
参考文献
[1] VesicAS .Expansionofcavitiesininfinitesoilmass[J].JSMFEASCE , 1972, 98 (SM 3) :265-290.
[2] PalmerAC .Undrainedplane strainexpansionofacylin dricalcavityinclay:asimpleinterpretationofthepres suremetertest[J].Geotechnique, 1972, 22 (3) :451-457.
[3] SkemptonAW , YassinAA , GibsonRE .Th啨oriedelaforceportantedespieux.Journ啨esdeM啨caniqiedesSols[J].AnnalesduB timentetdesTravauxPublics, 1953, 6 (63-64) :285-290.
[4] 王广国, 杜明芳, 苗兴城, 等.压密注浆机理研究及效果检验[J].岩石力学与工程学报, 2000, 5:670-673. WANGuang guo, DUMing fang, MIAOXing cheng, etal.Mechanismofcompactiongroutingandeffectex amination[J].ChineseJournalofRockMechanicsandEngineering, 2000, 5:670-673.
[5] 胡中雄.土力学与环境土工学[M].上海:同济大学出版社, 1997. HUZhong xiong.SoilMechanicsandEnvironmentalGeotechnics[M ].Shanghai:TongjiUniversityPress, 1997.
[6] CarterJP , BookerJR , KeungSK .Cavityexpansionincohesivefrictionalsoils[J].Geotechnique, 1972, 36 (2) :349-358.
[7] 王启铜, 龚晓南, 曾国熙.考虑土体拉、压模量不同时静压桩的沉桩过程[J].浙江大学学报自然科学版, 1992, 6:678-687. WANGQi tong, GONGXiao nan, ZENGGuo xi.Pilejackinginasoilmasswithdifferentelasticmoduliinten sionandcompression[J].JournalofZhejiangUniversity (NaturalScience) , 1992, 6:678-687.
[8] 樊良本, 朱国元.桩周土应力状态的圆柱孔扩张理论试验研究[J].浙江大学学报自然科学版, 1998, 2:228-235. FANLiang ben, ZHUGuo yuan.Experimentalstudyontheoryofexpansionofacylindricalcavitytoanalysingstressconditioninsaturatedclayaroundasinglepile[J].JournalofZhejiangUniversity (NaturalScience) , 1998, 2:228-235.
[9] 钱家欢, 殷宗泽.土工原理与计算[M].北京:中国水利水电出版社, 1996. QIANJia huan, YINZong ze.GeotechnicalTheoryandComputation[M ].Beijing:WaterResourcesandHydrowaterPress, 1996.