本文以长沙中建梅溪一号(3B地块)基坑工程为例,通过前期的地质调查、多方案论证,确定了初步的土钉墙支护结构和放坡挂网素喷的基坑支护方案,应用FLAC3D对支护方案模拟,模拟分析了采用土钉墙支护条件下,基坑开挖过程中基坑土体的应力场分布规律、支护结构受力状况,模拟结果显示:1) 从土钉受力情况来看,土钉轴力基本为拉力。同时,随着开挖工程的开展,土钉受力情况是随之而发生变化的,数值模拟结果显示,开挖过程中,土钉所承受的最大轴力为5.2 t。2) 土钉在长度方向上轴力分布遵循“中性点”原则,即中间数值大、两边数值小。最大轴力位置在潜在滑移面处。靠近坡面土钉段的轴力大于土钉深部段轴力这说明,土钉除了承受土体变形后施加的动力外,还受到喷射混凝土面层施加在土钉上的抗拔力作用。3) 随着开挖的进行,土钉受力逐渐增大。这主要是由于开挖深度增加,基坑边坡土体变形增加,将土体的位移力和土体膨胀力传递给土钉所致。 In this paper, taking the foundation pit project of No. 1 Meixi (3b block) of Changsha Zhongjian as an example, through the preliminary geological investigation and multi scheme demonstration, the preliminary soil nailing wall supporting structure and the foundation pit supporting scheme of sloping hanging net and plain spraying are determined. The FLAC3D is used to simulate the supporting scheme, and the stress field distribution rule and supporting structure of the foundation pit soil during the excavation are simulated and analyzed. The results show that: 1) The axial force of soil nail is tensile force. At the same time, with the development of excavation engineering, the stress of soil nails changes accordingly. The numerical simulation results show that the maximum axial force of soil nails is 5.2 t. 2) The axial force distribution of soil nail in the length direction follows the principle of “neutral point”, that is, the middle value is large and the two sides value is small. The maximum axial force is located at the potential slip surface. The axial force of the soil nail near the slope is greater than that of the deep part of the soil nail, which indicates that the soil nail is not only subjected to the force exerted by the deformation of the soil, but also to the pull-out force exerted by the shotcrete surface on the soil nail. 3) With the progress of excavation, the stress of soil nail increases gradually. This is mainly due to the increase of excavation depth and soil deformation of foundation pit slope, which transfers the displacement force and expansion force of soil to soil nails.
——以长沙中建梅溪湖项目基坑支护为例
李庆年
中国煤炭地质总局广西煤炭地质局,广西 柳州
收稿日期:2020年6月25日;录用日期:2020年7月21日;发布日期:2020年7月28日
本文以长沙中建梅溪一号(3B地块)基坑工程为例,通过前期的地质调查、多方案论证,确定了初步的土钉墙支护结构和放坡挂网素喷的基坑支护方案,应用FLAC3D对支护方案模拟,模拟分析了采用土钉墙支护条件下,基坑开挖过程中基坑土体的应力场分布规律、支护结构受力状况,模拟结果显示:1) 从土钉受力情况来看,土钉轴力基本为拉力。同时,随着开挖工程的开展,土钉受力情况是随之而发生变化的,数值模拟结果显示,开挖过程中,土钉所承受的最大轴力为5.2 t。2) 土钉在长度方向上轴力分布遵循“中性点”原则,即中间数值大、两边数值小。最大轴力位置在潜在滑移面处。靠近坡面土钉段的轴力大于土钉深部段轴力这说明,土钉除了承受土体变形后施加的动力外,还受到喷射混凝土面层施加在土钉上的抗拔力作用。3) 随着开挖的进行,土钉受力逐渐增大。这主要是由于开挖深度增加,基坑边坡土体变形增加,将土体的位移力和土体膨胀力传递给土钉所致。
关键词 :基坑支护,土钉墙支护,应力状态
Copyright © 2020 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
土钉是一种基于新奥隧道法原理发展起的一种基坑支护技术 [
本文长沙中建梅溪一号(3B地块)基坑支护方案为例,利用FLAC3D程序对其位移场及稳定性进行数值模拟分析,并与现场实测数据结果进行比较,以期得到一些有参考意义的结论,为以后本地区及相似工程地质条件的基坑支护研究工作提供一定的参考。
场地内共涉及11层土,分别为第1层:杂填土(Q4-3m1);第2层:新近沉积粉土(Q4-3a1);第3层:新近沉积粉质粘土夹粉土(Q4-3a1);第4层:新近沉积粉质粘土(Q4-3a1);第5层(Q4-21):粉土;第6层:粉质粘土(Q4-21);第7层:粉土(Q4-21);第8层:粉质粘土(Q4-21);第8夹层:粉土(Q4-21);第9层:粉土(Q4-21);10-1层:粉砂(Q4-1a1);第10层:细砂(Q4-1a1);第11层:粉质粘土(Q3a1)。
前期工程勘查资料表明,勘探深度范围内场地的地下水可分为两层,上层为潜水,下部为承压水,两者被相对隔水层第8层粉质粘土隔开。场地内自然水位约为地表下1.5 m。根据规范要求以及周边工程施工经验,建议上部潜水含水层的综合渗透系数取K = 0.5 m/d,下部承压水含水层的综合渗透系数取K = 8.0 m/d。
场地内各主要土层的物理力学指标见表1。
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | |
---|---|---|---|---|---|---|---|---|---|
土层名称 | 杂填土 | 沉积粉土 | 粉质粘土 | 沉积粉质粘土 | 粉土 | 粉质粘土 | 粉土 | 粉质粘土 | 粉土 |
厚度(m) | 0.3~2.35 | 2.43 | 1.84 | 2.46 | 3.39 | 1.69 | 1.83 | 3.12 | 1.44 |
体积模量(MPa) | 4.3 | 4.04 | 5.59 | 6.44 | 8.05 | 8.33 | 6.35 | 9.32 | 5.59 |
剪切模量(MPa) | 1 | 1.1 | 2.29 | 2.63 | 3.29 | 3.20 | 2.44 | 4.3 | 2.29 |
质量密度 (×103 Kg/m3) | 1.7 | 1.9 | 1.9 | 2.0 | 2.0 | 2.0 | 1.9 | 2.0 | 1.9 |
粘聚力(KPa) | 15 | 20 | 30 | 27 | 15 | 32 | 12 | 31 | 30 |
内摩擦角(˚) | 15 | 20 | 20 | 21 | 21 | 20 | 20 | 20 | 20 |
剪胀角 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
表1. 主要土层物理力学指标
本工程的基坑支护对三个1-1’、2-2’、3-3’剖面(图1,图2,图3)采用土钉支护的方案,对4-4’剖面采用放坡挂网素喷的支护方式,土钉支护方案的支护参数见表2,表3,表4;素喷坡面的设计参数见图4。
当各个三角形内的应变张量确定以后,然后初始设定的本构规律计算相应的应力。在此之前,需要对原有的应力进行两个方面的修正:1) 转动应力作修正;2) 扣除孔隙压力。应变与应力间的本构关系给成增量形式:
Δ σ i j = f ( Δ l i j , σ i j , ⋯ ) (1)
式中: Δ σ i j ——应力增量张量, i , j = 1 , 2 ;
f ( ) ——包含有增量应变,现存总应力状态,材料常数等的本构方程。
图1. 剖面1-1’示意图
图2. 剖面2-2’示意图
图3. 剖面3-3’示意图
排号 | 埋置深度(m) | 土钉长度(m) | 钢筋配置 | 水平间距(m) | 倾角(˚) | 成孔直径(mm) |
---|---|---|---|---|---|---|
1 | 1.3 | 8.0 | 1φ18 | 1.3 | 5~15 | 100 |
2 | 2.8 | 10.0 | 1φ18 | 1.3 | 5~15 | 100 |
3 | 4.3 | 9.0 | 1φ18 | 1.3 | 5~15 | 100 |
4 | 5.8 | 9.0 | 1φ18 | 1.3 | 5~15 | 100 |
5 | 7.3 | 6.0 | 1φ18 | 1.3 | 5~15 | 100 |
表2. 剖面1-1’土钉支护参数
排号 | 埋置深度(m) | 土钉长度(m) | 钢筋配置 | 水平间距(m) | 倾角(˚) | 成孔直径(mm) |
---|---|---|---|---|---|---|
1 | 1.3 | 8.0 | 1φ18 | 1.3 | 5~15 | 100 |
2 | 2.7 | 9.0 | 1φ18 | 1.3 | 5~15 | 100 |
3 | 4.1 | 8.0 | 1φ18 | 1.3 | 5~15 | 100 |
4 | 5.5 | 6.0 | 1φ18 | 1.3 | 5~15 | 100 |
表3. 剖面2-2’土钉支护参数
排号 | 埋置深度(m) | 土钉长度(m) | 钢筋配置 | 水平间距(m) | 倾角(˚) | 成孔直径(mm) |
---|---|---|---|---|---|---|
1 | 1.3 | 6.0 | 1φ18 | 1.3 | 5~15 | 100 |
2 | 2.8 | 9.0 | 1φ18 | 1.3 | 5~15 | 100 |
3 | 4.3 | 8.0 | 1φ18 | 1.3 | 5~15 | 100 |
4 | 5.8 | 6.0 | 1φ18 | 1.3 | 5~15 | 100 |
5 | 7.3 | 4.5 | 1φ18 | 1.3 | 5~15 | 100 |
表4. 剖面3-3’土钉支护参数
一旦四边形域内的应力张量已经求得,便可以用它来计算各网点处的不平衡力。重写运动方程式 [
p [ ∂ u ˙ i ∂ τ ] = ∂ σ i j ∂ x j + p g i (2)
其结果为
p [ ∂ u ˙ i ∂ τ ] = 1 A ∑ σ i j ε j k Δ x k + p g i (3)
或
∂ u ˙ i ∂ τ = F i m + g i
式中:
F i = ∑ σ i j ε j k Δ x k , m = p A
因此网点力变为:
F i = σ i j a ε j k ( x k ( 1 ) − x k ( 4 ) ) + σ i j b ε j k ( x k ( 2 ) − x k ( 1 ) ) + σ i j c ε j k ( x k ( 3 ) − x k ( 2 ) ) + σ i j d ε j k ( x k ( 4 ) − x k ( 3 ) ) (4)
这个力在所有应力和坐标,在时刻t,t + Δt等时都是已知的。在下一个半时间步长时的网点速度为:
u ˙ i ( t + Δ t 2 ) = u ˙ i ( t − Δ t 2 ) + F i t m Δ t + g i (5)
应该指出,如果域中出现任何一个错误,那么出现错误的域中的各项将由式(4)中直接省略掉。
如果模拟过程中,选择尺度较大的应变模型,则把修正以后的网格坐标用于式(5)的,求得各新速度:
x i ( t + Δ t ) = x i ( t ) + u ˙ i ( t + Δ t 2 ) Δ t (6)
当利用代码中的大应变逻辑时,由于物体的转动,必须修正应力。应力的正规转动方程式为:
σ ¯ i j = σ α β J i α J j β (7)
J i j = [ cos θ sin θ − sin θ cos θ ] (二维)
式中: σ ¯ i j ——大应变转动时修正了的应力张量。
当角度 Δ θ 非常小的条件下,因为 cos θ ≈ 1 , sin θ ≈ Δ θ ,则
J i j = [ 1 Δ θ − Δ θ 1 ] (二维)
因此,
σ ¯ 11 = σ 11 J 11 J 11 + σ 21 J 12 J 11 + σ 12 J 11 J 12 + σ 22 J 12 J 12 = σ 11 + σ 21 Δ θ + σ 12 Δ θ + σ 22 Δ θ 2 Δ σ 11 = 2 σ 12 Δ θ ( 因 为 σ 21 = σ 12 , Δ θ 很 小 时 , Δ θ 2 = 0 ) (8)
特殊条件下, σ 11 的修正项为
Δ σ 11 = − 2 σ 12 Δ θ (9)
同理,有
σ ¯ 22 = σ 11 J 21 J 21 + σ 21 J 22 J 21 + σ 12 J 21 J 22 + σ 22 J 22 J 22 = σ 11 Δ θ 2 − σ 21 Δ θ − σ 12 Δ θ + σ 22 2 = σ 22 − 2 σ 12 Δ θ + σ 11 Δ θ 2 (10)
反号并略去 Δ θ 2 项后,有
Δ σ 22 = − 2 σ 12 Δ θ (11)
还有
σ ¯ 12 = σ 11 J 11 J 21 + σ 21 J 12 J 21 + σ 12 J 11 J 22 + σ 22 J 12 J 22 = − σ 11 Δ θ − σ 21 Δ θ 2 + σ 12 + σ 22 Δ θ (12)
Δ σ 12 = ( σ 11 − σ 22 ) Δ θ (13)
当采用大应变逻辑时,把这些修正项加到各应力中。
本基坑位于长沙市岳麓区梅溪湖片区境内,南邻梅溪湖路,东邻E2路,场地交通便利,地理位置优越。实际基坑为长方形,基坑外轴线面积约102 × 371.9 m2,南北侧开挖深度为地表下6 m,中部为8 m。根据场地工程地质条件,基坑支护方案拟采用以下两种形式:1) 土钉墙支护;2) 放坡挂网素喷支护方案。据工程经验与现场实际情况,本次数值模拟模型尺寸设置为30 m × 20 m × 20 m。
模型中岩土物理力学参数按基坑支护方案设计中参数取值。
本次共建立三个力学模型,分别是工况一基坑开挖6 m (设计方案剖面2,位于基坑北侧)、工况二基坑开挖深8 m (设计方案剖面1,位于基坑东侧)、工况三基坑开挖深8 m (设计方案剖面3,位于基坑西侧)。各剖面数值模拟模型图分别见图4~6。
图4. 剖面二模型图
图5. 剖面一模型图
图6. 剖面三模型图
土钉采用FLAC3D中的锚索结构(cable),喷射混凝土用壳结构(shell-type)。土钉及混凝土层面的参数如表5和表6所示。
弹性模量/GP | 截面积/m2 | 屈服强度/KPa | 水泥黏结刚度/N/M | 水泥黏结强度/KPa |
---|---|---|---|---|
20 | 2.54 × 10−4 | 800 | 4 × 109 | 89 |
表5. 土钉材料参数
弹性模量/GP | 泊松比 |
---|---|
2.1 | 0.28 |
表6. 混凝土参数
在长期自重应力作用下,基坑在开挖施工前,施工区的地层土体已经处于稳定状态。所以在模型过程中,需要对模型进行应力初始化。开挖计算前,位移场和速度场需要清零,以避免出现误差。见图7所示。
从模拟结果来看,研究区基坑深度范围内土体垂直应力呈现由上至下呈递增分布的规律,在地表附近约为−0.05 MPa (负号表示为压力),直至基坑深处时,其数值大约为−0.38 MPa。模拟结果与现场土层的应力分布状态基本吻合。
1) 剖面1土钉轴力(图8~13)
图7. 应力初始化云图
图8. 开挖1.5 m土钉轴力图
图9. 开挖3 m土钉轴力图
图10. 开挖4.5 m土钉轴力图
图11. 开挖6 m土钉轴力图
图12. 开挖8 m土钉轴力图
图13. 开挖稳定后土钉轴力曲线
2) 剖面2土钉轴力(图14~18)
图14. 开挖1.5 m土钉轴力图
图15. 开挖3 m土钉轴力图
图16. 开挖4.5 m土钉轴力图
图17. 开挖6 m土钉轴力图
图18. 开挖稳定后土钉轴力曲线
3) 剖面3土钉轴力(图19~24)
图19. 开挖1.5 m土钉轴力图
图20. 开挖3 m土钉轴力图
图21. 开挖4.5 m土钉轴力图
图22. 开挖6 m土钉轴力图
图23. 开挖8 m土钉轴力图
图24. 开挖稳定后土钉轴力曲线
对基坑开挖过程中土钉受力进行了实时监测如图8~24所示,得出如下结论:
1) 从土钉受力情况来看,土钉轴力基本为拉力。同时,随着开挖工程的开展,土钉受力情况是随之而发生变化的,数值模拟结果显示,开挖过程中,土钉所承受的最大轴力为5.2 t。
2) 土钉在长度方向上轴力分布遵循“中性点”原则,即中间数值大、两边数值小。最大轴力位置在潜在滑移面处。靠近坡面土钉段的轴力大于土钉深部段轴力这说明,土钉除了承受土体变形后施加的动力外,还受到喷射混凝土面层施加在土钉上的抗拔力作用。
3) 随着开挖的进行,土钉受力逐渐增大。这主要是由于开挖深度增加,基坑边坡土体变形增加,将土体的位移力和土体膨胀力传递给土钉所致。
1) 从土钉受力情况来看,土钉轴力基本为拉力。数值模拟结果显示,开挖过程中土钉所承受的最大轴力为5.2 t。
2) 土钉在长度方向上轴力分布遵循“中性点”原则,即中间数值大、两边数值小。
3) 随着开挖的进行,土钉受力逐渐增大。
李庆年. 土钉墙支护及其应力特征数值模拟分析——以长沙中建梅溪湖项目基坑支护为例Numerical Simulation Analysis of Soil Nailing Wall Support and Its Stress Characteristics—Taking the Foundation Pit Support of Meixi Lake Project of Changsha Zhongjian as an Example[J]. 土木工程, 2020, 09(07): 715-730. https://doi.org/10.12677/HJCE.2020.97077