力学与实践, 2021, 43(3): 371-379 DOI: 10.6052/1000-0879-20-487

应用研究

考虑邻近建筑物影响的基坑支护结构受力响应分析1)

赵世永,2)

中铁十五局集团有限公司,上海 200070

中铁十五局集团城市轨道交通工程有限公司,河南洛阳 471000

MECHANICAL BEHAVIORS OF SUPPORTING STRUCTURE OF FOUNDATION PIT WITH CONSIDERATION OF THE INFLUENCE OF ADJACENT BUILDING1)

ZHAO Shiyong,2)

China Railway 15th Bureau Group Limited Company, Shanghai 200070, China

Urban Rail Transit Engineering Limited Company of China Railway 15th Bureau, Luoyang 471000, Henan, China

通讯作者: 2)赵世永,高级工程师,研究方向为地基基础工程。E-mail:150302137@qq.com

责任编辑: 王永会

收稿日期: 2020-11-19   网络出版日期: 2021-06-08

基金资助: 1)河南省住房城乡建设科技计划资助项目(K-1817)
河南省住房城乡建设科技计划资助项目(K-1818)
河南省住房城乡建设科技计划资助项目(K1816)
河南省住房城乡建设科技计划资助项目(K-1940)

Received: 2020-11-19   Online: 2021-06-08

作者简介 About authors

摘要

为了深入了解临近建筑物对基坑支护结构受力响应的影响,通过学习前人的研究成果和基于自身经验采用三参数地基弹性梁模型,来考虑临近建筑物与支护结构间有限宽土压力、基底应力等的影响,根据支护结构受力特点,得到一种考虑临近建筑影响的支护结构计算分析模型;采用桩身离散的方法创建了支护结构的内力及其变形计算的矩阵表达式,得到一种半解析解答的方法;并根据工程实例对该方法进行计算验证,对建筑基底应力、距支护结构距离的影响规律进行了分析,验证了本方法的可行性,可为城镇建筑密集区域的基坑工程类似支护设计提供参考。

关键词: 基坑工程; 邻近建筑; 受力响应; 理论计算; 半解析解

Abstract

In order to better understand the mechanical behaviors of the supporting structure of the deep foundation pit, a tri-parameter foundation model is used, according to different stress patterns of the supporting structure, to build a computational analysis model with consideration of the influence of the adjacent building. In this computational analysis model, the additional stress caused by the adjacent building and the pressure of a certain width soil between the supporting structure and the adjacent building are considered. The matrix expressions for the internal forces and the deformations of the supporting structure are derived by using the pile discrete method and the matrix transfer method, and the semi-analytical solution for the supporting structure is obtained. For an excavation project, the influences of the additional stress caused by the adjacent building and the distance between the supporting structure and the adjacent building are analyzed, this method is verified, and it can be used for the design of this type of supporting projects in dense urban areas.

Keywords: foundation pit engineering; adjacent building; mechanical response; theoretical calculation; semi-analytical solution

PDF (724KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

赵世永. 考虑邻近建筑物影响的基坑支护结构受力响应分析1). 力学与实践, 2021, 43(3): 371-379 DOI:10.6052/1000-0879-20-487

ZHAO Shiyong. MECHANICAL BEHAVIORS OF SUPPORTING STRUCTURE OF FOUNDATION PIT WITH CONSIDERATION OF THE INFLUENCE OF ADJACENT BUILDING1). MECHANICS IN ENGINEERING, 2021, 43(3): 371-379 DOI:10.6052/1000-0879-20-487

目前,我国的基础设施建设正处于快速发展的时期,有许多的城市开始进入到地下空间的高速建设阶段,其"规模化、系统化、综合化"成为了这一阶段较为显著的特征。而且随着国家的节能省地型城市建设政策的不断落实,建筑密度不断加大,城市中心地带建筑密集区高层、超高层建筑、地下商场、地下通道等工程大量涌现,与之对应的基坑开挖工程与周围建构筑物的相互影响问题不可避免。目前,针对基坑开挖对临近建筑基础的影响研究逐渐深入。Son等[1]通过数值模拟分析了地层开挖对邻近建筑物变形的影响;李志伟等[2]采用有限元分析的手段,给出了基坑开挖影响下周边建筑物相对挠曲变形的评估方法;张治国等[3-4]和王卫东等[5]依据上海地区土质等实际工程中存在的特点,研究了深基坑开挖与周围浅基础沉降之间的规律,并给出了基坑开挖对邻近浅基础影响的简化评价方法。然而,这些研究多集中于基坑开挖对临近建筑基础的影响,而在基坑开挖中如何合理考虑临近既有建筑物基础的影响来进行支护结构的设计,成为工程设计人员必须合理处理并进行深化研究的课题。

考虑到周边建筑基础对基坑支护结构的影响,金亚兵等[6]通过理论假定给出了基坑坡顶局部超载作用下土压力的简化计算图式;张莲花等[7]基于该方法,考虑超载大小、距基坑边缘距离及其作用宽度和相对深度等因素的影响,依据实际工程对采用不同基坑支护结构的受力响应进行了对比分析研究;张浩等[8]针对不平衡堆载工况,通过研究给出了邻近结构工程桩的受力计算分析方法,为基坑支护结构的受力响应理论分析提供了可供借鉴的经验;李峰等[9]通过研究临近建筑物基础埋深较深的工况(如地下室、地铁车站等),提供了一种有限宽土体作用在基坑支护结构上土压力的计算方法,进一步深化了基坑开挖对附近建筑基础影响的相关研究。然而,这些成果大多是和悬臂式基坑支护方面相关的受力研究,没有考虑受附近建筑基础影响下的多支点基坑支护方式(如:墙锚支护等)的相关受力分析,但在我国应用此类支护结构的基坑工程很常见。张浩等[10]依据弹性地基梁模型,给出了考虑到基坑内预留土影响下的多支点支护结构受力变形计算模型,但未考虑基坑附近建筑物基础局部超载对支护结构产生的影响。因此,考虑临近基础影响的基坑支护结构的设计计算理论仍有待进一步深入的研究。

为此,考虑到临近建筑基础超载产生的影响,本文依据三参数弹性地基梁模型的方法,结合多支点支护结构的受力特点,建立支护结构受力变形计算分析模型,然后支护结构受力响应的半解析解答可通过桩身离散与矩阵传递法推导得出,从而为合理计算临近建筑基础超载对多支点支护结构产生的作用提供参考。

1 计算模型与基本假定

1.1 支护结构受力控制方程

考虑临近建筑(如地下室、地铁车站等)影响的典型支护结构受力特征见图1。在图1中,$h_{\rm c}$表示支护结构与临近建筑间有限宽土体影响作用区域;$ h_{\rm u}$为支护结构的下部临空段;$q$为临近建筑基底超载应力;$\Delta q$为基底局部超载对支护结构的附加土压力;$ q_{\rm d}$为嵌固段坑底土对支护结构产生的水平向抗力。实际工程中,基坑开挖势必引起周围土体发生向坑内的位移变形,而对支护结构产生土压力作用($p_{\rm a}$);支护结构受土压力作用亦会发生向坑内的侧向位移,而受到坑内被动区土的水平抗力作用($q_{\rm d}$)。由此可见,支护结构与土的相互作用属于典型的被动桩桩土相互作用模式[8,11]

图1

图1   支护结构内力计算模型


据此,考虑临近建筑影响,根据开挖支护结构受力特征,可将支护桩划分为两段,即开挖深度范围的临空段和坑底以下的嵌固段。其中,考虑受力不同,临空段又可分为:①上部受有限宽土体作用的临空段,受坑外有限宽土体土压力作用和支点(锚杆或内支撑等)集中力作用;②下部临空段,不仅受坑外主动土压力的作用,还受到因建筑基底超载引起的附加载荷($\Delta q$)作用;嵌固段则同时受到坑外主动土压力和坑内土抗力作用。进而,以临空段和嵌固段顶点为原点分别建立独立坐标系$y_{\rm c}-z_{\rm c}$和$y_{\rm d}-z_{\rm d}$,根据静力平衡,忽略轴力影响,可建立各桩段的控制微分方程

(1)临空段

$ \begin{eqnarray} EI\frac{{\rm d}^{4}y_{\rm c} }{{\rm d}z_{\rm c}^{4}}=p_{\rm a}\left( {z_{\rm c} } \right) \end{eqnarray} $

(2)嵌固段

$ \begin{eqnarray} EI\frac{{\rm d}^{4}y_{{\rm d}} }{{\rm d}z_{{\rm d}}^{4}}+q_{{\rm d}}\left( {y_{{\rm d}} ,z_{{\rm d}} } \right)=p_{\rm a} \left({z_{{\rm d}} } \right) \end{eqnarray} $

式中,$EI$为支护桩刚度,$q_{\rm d}(y_{\rm d}$, $z_{\rm d})$为坑底被动区土的抗力,$p_{\rm a}(z_{\rm c})$和$p_{\rm a}(z_{\rm d})$为支护结构外侧土压力。

1.2 计算参数的确定

图1所示,由于临近建构筑物(地下室、地铁车站等)的存在,支护结构受其影响承受的是有限土体的土压力,若仍根据朗肯土压力理论计算,易导致土压力计算偏大,造成浪费。因此,需对该区域土压力做进一步修正。假定临近建构筑物埋深为$D$,受基坑开挖影响,有限土体极限滑裂面倾角与水平面呈$\theta$开展,基于滑楔体平衡理论,可求得该状态下有限宽土体作用于基坑支护桩上的主动土压力[9]

$ \begin{eqnarray} && E_{\rm a} = {\left( {\frac{1}{2}\gamma S-k} \right)\left( {2h_{\rm c} -S{\rm ctg}\theta } \right)}\frac{\cos \left( {\theta +\varphi } \right)}{\sin \left( {\theta +\delta +\varphi } \right)}-\\ && \qquad \frac{cS}{\cos \delta }\left[ {1+\frac{\cos \left( {\theta +\delta } \right)\cos \left( {\theta +\varphi } \right)}{\sin \theta \sin \left( {\theta +\delta +\varphi } \right)}} \right] \end{eqnarray} $

式中,$S$为有限土体宽度,$\gamma $为土体容重,$c$为黏聚力,$\varphi$为内摩擦角,$\delta $为土体的外摩擦角,$k$为黏着力。

为便于计算,根据土压力强度的线性关系,由式(3)可求得有限宽土体作用段分布土压力

$ \begin{eqnarray} p_{\rm a} \left( {z_{\rm c} } \right)={{2E_{\rm a} z_{\rm c}}/{h_{\rm c}^{2}}} \end{eqnarray} $

考虑建构筑物基底超载的影响,基底应力$q$沿$\theta $扩散,在区间[$h_{\rm c}$,$h_{\rm c}+B+2(h_{\rm c}-D)\tan\theta$]内对支护结构产生附加土压力

$ \begin{eqnarray} \Delta q=\frac{qB}{B+2\left( {h_{\rm c} -D} \right)\tan \theta } \end{eqnarray} $

式中,$B$为临近建构筑物基础宽度。

下部临空段外侧的主动土压力

$ \begin{eqnarray} p_{\rm a} \left( {z_{\rm c} } \right)=\gamma \left( {z_{\rm c} -D} \right)K_{\rm a} -2c\sqrt {K_{\rm a} } \end{eqnarray} $

式中,$K_{\rm a}$为主动土压力系数,$K_{\rm a} =\tan^{2}(45^{\circ}-\varphi /2)$。受基底局部载荷影响区域的土压力为

$ \begin{eqnarray} p_{\rm a} \left( {z_{\rm c} } \right)=\gamma \left({z_{\rm c} -D} \right)K_{\rm a} -2c\sqrt {K_{\rm a} }+\Delta q \end{eqnarray} $

据此,综合式(4) $\sim\!$式(7),基坑范围内支护结构临空段外侧土压力可统一表示为

$ \begin{eqnarray} p_{\rm a} \left( {z_{\rm c} } \right)=a'z_{\rm c} +b' \end{eqnarray} $

式中,$a'$和$b'$为常数。

参考文献[10],坑底以下支护结构外侧土压力为一定值,其大小与坑底处土压力相等,则嵌固段外侧土压力为

$ \begin{eqnarray} p_{\rm a} \left( {z_{{\rm d}} } \right)=c' \end{eqnarray} $

式中,$c'$为常数。

基坑底部的被动区土体会对基坑支护桩的受力产生影响,考虑到这一点,根据弹性地基梁法,将坑内土体和支护桩间存在的相互作用通过建立被动受压弹簧来模拟,只需选择符合工程实际的土层弹簧刚度值,就可以通过支护桩位移模式及其大小,来求得基坑内分布在支护结构上的土体抗力

$ \begin{eqnarray} q_{{\rm d}} \left( {y_{{\rm d}} ,z_{{\rm d}} } \right)=k_{{\rm d}} \left( {z_{{\rm d}} } \right)b_{1} y_{{\rm d}} \end{eqnarray} $

式中,$b_{1}$是基坑支护桩的计算宽度;$ k_{\rm d}(z_{\rm d})$为坑底嵌固段的地基抗力系数,假定坑内土抗力与支护桩侧向位移成正比,采用三参数($m$,$z_0$,$n$)地基抗力模型,有:$k_{{\rm d}} \left( {z_{{\rm d}} } \right)=m\left( {z_{\rm s} +z_{{\rm d}} } \right)^{n}$,$m$为地基比例系数,$n$为土体抗力的深度指数,$z_{\rm s}$为坑底处的当量深度,此处考虑基坑开挖后坑底土体变成超固结土,虽然发生应力释放,但$z_{\rm s}$仍具有一定的刚度。

2 支护结构内力响应计算

2.1 控制方程求解

为了简化计算和便于编程,对支护结构进行离散,采用矩阵传递系数法进行统一求解。

2.1.1 临空段控制方程求解

考虑支护结构支点集中力与外侧土压力分界的影响,将临空段离散为$N_{\rm c}$份,每段长度$h_{1}=h_{\rm c}/N_{\rm c}$,并保证集中力作用于离散节点处;取任意段$i$建立独立坐标系进行分析,如图2所示。为简化分析,将微段外侧主动土压力取为

$ \begin{eqnarray} p_{{\rm a}i} ={{a\cdot \left( {2i-1} \right)h_{1} }/2}+b \end{eqnarray} $

图2

图2   临空段离散示意图


由式(1)将临空段第$i$桩段控制方程转化为

$ \begin{eqnarray} \frac{{\rm d}^{4}y_{{\rm c}i} }{{\rm d}z_{{\rm c}i}^{4}}=\overline{p_{{\rm a}i} } \end{eqnarray} $

式中,$\overline {p_{{\rm a}i} } ={{p_{{\rm a}i} }/({EI})}$。对微分方程(11)求解可得

$ \begin{eqnarray} y_{{\rm c}{i}} =C_{{\rm o}1} +C_{{\rm o}2} z_{{\rm c}i} +C_{{\rm o}3} z_{{\rm c}i}^{2}+C_{{\rm o}4} z_{{\rm c}i}^{3}+\frac{\overline {p_{{\rm a}i} } }{24}z_{{\rm c}i}^{4} \end{eqnarray} $

式中,$y_{{\rm c}{i}}$为第$i$段任意位置$z_{{\rm c}{i}}$处的桩身挠曲;$C_{\rm o1}$, $C_{\rm o2}$,$C_{\rm o3}$,$C_{\rm o4}$为常系数。

设$\varphi_{{\rm c}{i}}$,$M_{{\rm c}{i}}$和$V_{{\rm c}{i}}$分别为$z_{{\rm c}{i}}$处桩身转角、弯矩和剪力,则由材料力学基本理论可得

$ \begin{eqnarray} \varphi_{{\rm c}i} =\frac{{\rm d}y_{{\rm c}i} }{{\rm d}z_{{\rm c}i} },\ \ M_{{\rm c}i} =EI\frac{{\rm d}^{2}y_{{\rm c}i} }{{\rm d}z_{{\rm c}i} ^{2}},\ \ V_{{\rm c}i} =EI\frac{{\rm d}M_{{\rm c}i} }{{\rm d}z_{{\rm c}i} } \end{eqnarray} $

设第$i$段顶端($z_{{\rm c}{i}}=0$)响应参量为$y_{{\rm c}{i}0}$, $\varphi_{{\rm c}{i}0}$, $M_{{\rm c}{i}0}$和$V_{{\rm c}i0}$,联合式(12)和式(13)可求出$C_{{\rm o}1}$, $C_{{\rm o}2}$, $C_{{\rm o}3}$, $C_{{\rm o}4}$;进而将其回代入式(12)和式(13)可得

$ \begin{eqnarray} {\mathbf U}_{{\rm c}i} \left( {h_{1} } \right)={\mathbf S}_{{\rm c}i} \left( {h_{1} } \right)\cdot {\mathbf U}_{{\rm c}i0} \end{eqnarray} $

式中,${\mathbf U}_{{\rm c}i} \left( {h_{1} } \right)=\left[ y_{{\rm ch}},\varphi_{{\rm ch}} ,M_{{\rm ch}} ,V_{{\rm ch}} ,1\right]^{\rm T}$, 其中$y_{{\rm c}{h}} $, $\varphi_{{\rm c}{h}}$, $M_{{\rm c}{h}} $, $V_{{\rm c}{h}}$分别为第$i$段底部($z_{{\rm c}i}=h_{1}$)的水平位移、转角、弯矩和剪力;${\mathbf U}_{{\rm c}i0} =[y_{{\rm c}i0},\varphi_{{\rm c}i0}$, $M_{{\rm c}i0} ,V_{{\rm c}i0},1]^{\rm T}$, $y_{{\rm c}i0} $, $\varphi_{{\rm c}i0}$, $M_{{\rm c}i0} $, $V_{{\rm c}i0}$不仅为第$i$段顶部的水平位移、转角、弯矩和剪力,还是第$i-1$段底部的水平位移、转角、弯矩和剪力,即${\mathbf U}_{{\rm c}\left( {i-1} \right)} \left( {h_{1} } \right)$;${\mathbf S}_{{\rm c}i} \left( {h_{1} } \right)$为第$i$微段的系数矩阵

$$\begin{eqnarray*} {\mathbf S}_{{\rm c}i} \left( {h_{1} } \right)=\left[ {{\begin{array}{*{20}c} {A_{{\rm c}1} } & {B_{{\rm c}1} } & {C_{{\rm c}1} } & {D_{{\rm c}1} } & {E_{{\rm c}1} } \\ {A_{{\rm c}2} } & {B_{{\rm c}2} } & {C_{{\rm c}2} } & {D_{{\rm c}2} } & {E_{{\rm c}2} } \\ {A_{\rm c3}} & {B_{\rm c3} } & {C_{\rm c3} } & {D_{\rm c3} } & {E_{\rm c3} } \\ {A_{\rm c4} } & {B_{\rm c4} } & {C_{\rm c4} } & {D_{\rm c4} } & {E_{\rm c4} } \\ 0 & 0 & 0 & 0 & 1 \\ \end{array} }} \right] \end{eqnarray*}$$

$A_{{\rm c}j}$, $B_{{\rm c}j}$,$C_{{\rm c}j}$,$D_{{\rm c}j}$,$E_{{\rm c}j}$ ($j=1\sim4$)为矩阵方程的20个系数。

考虑变形受力连续性,基于式(14)可得

$ \begin{eqnarray} {\mathbf U}_{{\rm c}i} \left( {h_{1} } \right)={\mathbf S}_{{\rm c}i} \left( {h_{1} } \right)\cdot {\mathbf U}_{{\rm c}\left( {i-1} \right)} \left( {h_{1} } \right) \end{eqnarray} $

考虑到锚杆、内支撑等支点集中力对支护结构的作用,假设$F_{j}$为支点$x_{j}$处的集中力,则$x_{j}$点处上下截面的变形内力关系为

$ \begin{eqnarray} {\mathbf U}_{x}^{{\rm d}} ={\mathbf S}_{{\rm c}x} {\mathbf U}_{x}^{\rm u} \end{eqnarray} $

式中,${\mathbf U}_{x}^{{\rm d}} $和${\mathbf U}_{x}^{\rm u}$分别为临空段$x_{j}$节点处下截面和上截面的内力变形参量矩阵,${\mathbf S}_{{\rm c}x} $为集中力引起的突变矩阵

$$\begin{eqnarray*} {\mathbf S}_{{\rm c}x} =\left[ {\begin{array}{c@{\ \ }c@{\ \ }c@{\ \ }c@{\ \ }c} 1&0&0&0&0 \\ 0&1&0&0&0 \\ 0&0&1&0&0 \\ 0&0&0&1&F_{j} \\ 0&0&0&0&1 \\ \end{array}} \right] \end{eqnarray*}$$

进而,可得临空段受力响应的矩阵传递方程

$ \begin{eqnarray} &&{\mathbf U}_{{\rm c}N_{\rm c} } \left( {h_{1} } \right)={\mathbf S}_{{\rm c}N_{\rm c} } \left( {h_{1} } \right)\cdot \cdot \cdot {\mathbf S}_{{\rm c}x} \cdot \cdot \cdot {\mathbf S}_{{\rm c}1} \left( {h_{1} } \right)\cdot {\mathbf U}_{{\rm c}0} =\\&&\qquad {\mathbf S}_{{\rm c}} \cdot {\mathbf U}_{{\rm c}0} \end{eqnarray} $

式中,${\mathbf S}_{{\rm c}} $为支护桩临空段总的系数矩阵;${\mathbf U}_{{\rm c}0}=\left[ y_{{\rm c}0} ,\varphi_{{\rm c}0} ,M_{{\rm c}0}, V_{{\rm c}0} ,1\right]^{\rm T}$,为桩顶处的内力变形参量矩阵;${\mathbf U}_{{\rm c}N_{\rm c} } \left( {h_{1} } \right)=\left[ y_{{\rm c}N_{\rm c} } ,\varphi_{{\rm c}N_{\rm c} } ,M_{{\rm c}N_{\rm c} } ,V_{{\rm c}N_{\rm c} } ,1\right]^{\rm T}$,为临空段和嵌固段交界面处的参量矩阵。

2.1.2 嵌固段控制方程求解

图3所示,将嵌固段分为$N_{\rm d}$份,每段长$h_{2}= h_{\rm d}/N_{\rm d}$,取任一段$i$建立独立坐标系进行分析,微段外侧土压力取为

$ \begin{eqnarray} p_{{\rm a}i} =c' \end{eqnarray} $

图3

图3   嵌固段离散示意图


考虑坑底地基土的抗力作用,基于式(10)采用中值定理,将第$i$微段的土体抗力系数简化为

$ \begin{eqnarray} &&k_{{\rm d}i} ={{\int_{\left( {i-1} \right)h_{2} }^{ih_{2} } {m\left( {z_{\rm s} +z_{{\rm d}} } \right)^{n}} {\rm d}z_{{\rm d}} }/{h_{2} }}= \\&&\qquad \frac{m\left( {z_{\rm s} +ih_{2} } \right)^{n+1}-m\left[ {z_{\rm s} +\left( {i-1} \right)h_{2} } \right]^{n+1}}{\left( {n+1} \right)h_{2} }\qquad \end{eqnarray} $

由式(2)可将第$i$微段控制方程转化为

$ \begin{eqnarray} \frac{{\rm d}^{4}y_{{\rm d}i} }{{\rm d}z_{{\rm d}i}^{4}}+\alpha_{{\rm d}i}^{4}y_{{\rm d}i} =t \end{eqnarray} $

式中,$\alpha_{{\rm d}i}^{4}={{k_{{\rm d}i} b_{1} }/({EI})}$,$t={{p_{{\rm a}i}}/({EI})}$,$b_{1}$为计算宽度。进而,对该微分方程求解可得

$ \begin{eqnarray} &&y_{{\rm d}i} ={\rm e}^{gz_{{\rm d}i} }\left( {C_{{\rm d}1} \cos gz_{{\rm d}i} +C_{{\rm d}2} \sin gz_{{\rm d}i} } \right)+ \\&&\qquad {\rm e}^{-gz_{{\rm u}i} }\left( {C_{{\rm d}3} \cos gz_{{\rm d}i} +C_{{\rm d}4} \sin gz_{{\rm d}i} } \right)+\\&&\qquad {t/{\alpha_{{\rm d}i}^{4}}} \end{eqnarray} $

式中,$y_{{\rm d}i}$为$z_{{\rm d}i}$处桩身挠曲变形;$g=\alpha_{{\rm d}i}/\sqrt{2}$; $C_{{\rm d}1}$,$C_{{\rm d}2}$,$C_{{\rm d}3}$,$C_{{\rm d}4}$为常系数。

与临空段推导方法相同,可得$i$微段矩阵方程

$ \begin{eqnarray} {\mathbf U}_{{\rm d}i} \left( {h_{2} } \right)={\mathbf S}_{{\rm d}i} \left( {h_{2} } \right)\cdot {\mathbf U}_{{\rm d}i0} \end{eqnarray} $

式中, ${\mathbf U}_{{\rm d}i} \left( {h_{2}} \right)=\left[{y_{{\rm d}h} ,\varphi_{{\rm d}h} ,M_{{\rm d}h} ,V_{{\rm d}h} ,1}\right]^{\rm T}$, $y_{{\rm d}h} $, $\varphi_{{\rm d}h}$, $M_{{\rm d}h} $, $V_{{\rm d}h}$为第$i$段底部($z_{{\rm d}i}=h_{2}$)的水平位移、转角、弯矩和剪力; ${\mathbf U}_{{\rm d}i0} =\big[ y_{{\rm d}i0}$, $\varphi_{{\rm d}i0}$, $M_{{\rm d}i0}$, $V_{{\rm d}i0}, 1 \big]^{\rm T}$, $y_{{\rm d}i0} $, $\varphi_{{\rm d}i0} $, $M_{{\rm d}i0} $, $V_{{\rm d}i0}$不仅为第$i$段顶部的水平位移、转角、弯矩和剪力, 还是第$i-1$段底部的水平位移、转角、弯矩和剪力, 即${\mathbf U}_{{\rm d}\left( {i-1} \right)} \left( {h_{2} } \right)$;${\mathbf S}_{{\rm d}i} \left( {h_{2} } \right)$为第$i$微段的系数矩阵

$$\begin{eqnarray*} {\mathbf S}_{{\rm d}i} \left( {h_{2} } \right)=\left[ {{\begin{array}{*{20}c} {A_{{\rm d}1} } & {B_{{\rm d}1} } & {C_{{\rm d}1} } & {D_{{\rm d}1} } & {E_{{\rm d}1} } \\ {A_{\rm d2} } & {B_{\rm d2} } & {C_{\rm d2} } & {D_{\rm d2} } & {E_{\rm d2} } \\ {A_{\rm d3} } & {B_{\rm d3} } & {C_{\rm d3} } & {D_{\rm d3} } & {E_{\rm d3} } \\ {A_{\rm d4} } & {B_{\rm d4} } & {C_{\rm d4} } & {D_{\rm d4} } & {E_{\rm d4} } \\ 0 & 0 & 0 & 0 & 1 \\ \end{array} }} \right] \end{eqnarray*}$$

其中$A_{{\rm d}j}$,$B_{{\rm d}j}$,$C_{{\rm d}j}$,$D_{{\rm d}j}$,$E_{{\rm d}j}$ ($j=1,2,3,4$)为矩阵方程的20个系数,是微段长度$h_{2}$的函数。

考虑各微段变形受力连续性,由式(22)可得

$ \begin{eqnarray} {\mathbf U}_{{\rm d}i} \left( {h_{2} } \right)={\mathbf S}_{{\rm d}i} \left( {h_{2} } \right)\cdot {\mathbf U}_{{\rm d}\left( {i-1} \right)} \left( {h_{2} } \right) \end{eqnarray} $

进而,可得整个桩段受力响应的矩阵传递方程

$ \begin{eqnarray} && {\mathbf U}_{{\rm d}N_{\rm d} } \left( {h_{2} } \right)={\mathbf S}_{{\rm d}N_{{\rm d}} } \left( {h_{2} } \right)\cdot \cdot \cdot {\mathbf S}_{{\rm d}2} \left( {h_{2} } \right){\mathbf S}_{{\rm d}1} \left( {h_{2} } \right)\cdot {\mathbf U}_{{\rm d}0} =\\&& \qquad {\mathbf S}_{{\rm d}} \cdot {\mathbf U}_{{\rm d}0} \end{eqnarray} $

式中, ${\mathbf U}_{{\rm d}0} =\left[ {y_{{\rm c}N_{\rm c}} ,\varphi_{{\rm c}N_{\rm c} } ,M_{{\rm c}N_{\rm c} } ,V_{{\rm c}N_{\rm c}} ,1} \right]^{\rm T}$,为临空段和嵌固段交界面处的内力变形参量矩阵;${\mathbf U}_{{\rm d}N_{\rm d} } \left( {h_{2} } \right)={\mathbf U}_{{\rm d}L}=\left[ {y_{{\rm d}L} ,\varphi_{{\rm d}L} ,M_{{\rm d}L},V_{{\rm d}L} ,1}\right]^{\rm T}$,为支护桩底端处的内力变形参量矩阵;${\mathbf S}_{{\rm d}i}\left( {h_{2} } \right)$为第$i$段的系数矩阵;${\mathbf S}_{{\rm d}}$为嵌固段总的系数矩阵。

2.2 连续条件与求解方法

由临空段与嵌固段的连续条件可知

$ \begin{eqnarray} {\mathbf U}_{{\rm d}0} ={\mathbf U}_{{\rm c}N_{\rm c} } \left( {h_{1} } \right) \end{eqnarray} $

据此,联合式(18)和式(25)可得整个支护桩的内力变形矩阵方程

$ \begin{eqnarray} {\mathbf U}_{{\rm d}L} ={\mathbf S}_{{\rm d}} {\mathbf S}_{{\rm c}} \cdot {\mathbf U}_{{\rm c}0} ={\mathbf S}\cdot {\mathbf U}_{{\rm c}0} \end{eqnarray} $

式中,${\mathbf S}$即为整个支护桩的总系数矩阵。

矩阵方程(27)中涉及到支护桩顶边界参量$y_{{\rm c}0}$, $\varphi_{{\rm c}0}$, $M_{{\rm c}0}$, $V_{{\rm c}0}$和桩端边界参量$y_{{\rm d}L}$, $\varphi_{{\rm d}L}$, $M_{{\rm d}L}$, $V_{{\rm d}L}$。考虑到不同的桩顶和桩端约束条件, 那么桩顶自由时

$ \begin{eqnarray} M_{{\rm c}0} =M_{0},\ \ V_{{\rm c}0} =H_{0} \end{eqnarray} $

桩顶固定时

$ \begin{eqnarray} y_{{\rm c}0} =0,\ \ \varphi_{{\rm c}0} =0 \end{eqnarray} $

桩端自由时

$ \begin{eqnarray} M_{{\rm d}{L}} =0,\ \ V_{{\rm d}{L}} =0 \end{eqnarray} $

桩端嵌固时

$ \begin{eqnarray} y_{{\rm d}{L}} =0,\ \ \varphi_{{\rm d}{L}}=0 \end{eqnarray} $

据此,可通过以下步骤计算支护结构的变形内力。

(1)将支护桩顶已知条件式(28)或式(29)、桩端已知条件式(30)或式(31)代入到矩阵方程(26)中,可得到有关4个未知量的4个方程组成的方程组,通过求解该方程组可得未知的桩顶和桩端边界参量;

(2)根据桩顶边界${\mathbf U}_{{\rm c}0} =\left[ {y_{{\rm c}0} ,\varphi_{{\rm c}0} ,M_{{\rm c}0} ,V_{{\rm c}0} ,1}\right]^{\rm T}$和递推式(16)可求得临空段中任一微段$i$ ($1\leqslant i \leqslant N_{\rm c}$)的下截面内力变形

$ \begin{eqnarray} {\mathbf U}_{{\rm c}i} \left( {h_{1} } \right)={\mathbf S}_{{\rm c}i} \left( {h_{1} } \right)\cdot \cdot \cdot {\mathbf S}_{{\rm c}x} \cdot \cdot \cdot {\mathbf S}_{{\rm c}1} \left( {h_{1} } \right)\cdot {\mathbf U}_{{\rm c}0} \end{eqnarray} $

当$i=N_{\rm c}$时,式(32)可求得临空段底端的内力变形${\mathbf U}_{{\rm c}N_{\rm c} } \left( {h_{1} } \right)=\left[ {y_{{\rm c}N_{\rm c} },\varphi_{{\rm c}N_{\rm c} } ,M_{{\rm c}N_{\rm c} },V_{{\rm c}N_{\rm c} } ,1} \right]^{\rm T}$;

(3)由临空段与嵌固段交界面处的连续条件(${\mathbf U}_{{\rm d}0} ={\mathbf U}_{{\rm c}N_{\rm c} } \left( {h_{1} }\right)$)可得嵌固段顶端处的内力变形参量${\mathbf U}_{{\rm d}0} =\left[{y_{{\rm c}N_{\rm c} } ,\varphi_{{\rm c}N_{\rm c} } ,M_{{\rm c}N_{\rm c} },V_{{\rm c}N_{\rm c} } ,1}\right]^{\rm T}$,进而代入递推公式(24)可求得嵌固段中任一微段$i$ $(1\leqslant i\leqslant N_{\rm d}$)的下截面内力变形

$ \begin{eqnarray} &&{\mathbf U}_{{\rm d}i} \left( {h_{2} } \right)=\prod\limits_{j=1}^i {{\mathbf S}_{{\rm d}j} \left( {h_{2} } \right)\cdot } {\mathbf U}_{{\rm d}0} =\\&&\qquad \prod\limits_{j=1}^i {{\mathbf S}_{{\rm d}j} \left( {h_{2} } \right)\cdot } {\mathbf S}_{{\rm c}} {\mathbf U}_{{\rm c}0} \end{eqnarray} $

综上所述,通过对支护桩的合理离散,由式(32)和式(33)可求得各节点处支护桩的变形内力值。

3 工程实例分析

郑州市黄河路西延隧道工程位于郑州市金水西路西站北街以东,黄河路下穿北编组站隧道以西,呈西南至东北走向,与郑州轨道交通5号线西站街站至沙口路站区间走向基本重合,采用明挖顺筑法施工。区间场地平坦开阔,场地地貌属于山前冲洪积缓倾平原,地层主要为第四系上更新统地层,主要有砂质粉土、黏质粉土、粉质黏土,夹有粉细砂,主要地层物理力学指标如表1所示。其中,区段DK0$+$943.450$\sim$DK0$+$987.395临近有一多层建筑,地上6层,地下1层为地下室,基础埋深4.85 m,与围护桩净距3.8 m。临近该建筑的围护结构采用钻孔灌注桩$+$混凝土支撑/预应力锚索支护体系,采用1道混凝土支撑和混凝土支撑下设10道竖向等距锚索进行支护,如图4所示。基坑开挖深度为23 m,钻孔灌注桩的直径为1.2 m,其长度为34 m;锚索竖向间距为2.0 m,采用抗拉设计强度值为1320 MPa的钢绞线。

表1   土层物理力学参数

新窗口打开| 下载CSV


图4

图4   断面示意图


利用本文提出的方法对该区段围护结构的受力变形进行计算分析。在计算坑外分布在支护结构上的土压力时,考虑到支护结构深度范围内土层较多,可以对土层指标进行加权平均来合理简化计算。经过对土层参数加权平均计算后可得:$\gamma=19.6$ kN/m$^{3}$,$c=17.2$ kPa,$\varphi=24.6^\circ$;嵌固段土层主要为黏质粉土层和粉质黏土层,考虑坑内卸荷影响,取其综合抗力系数为150 kN/m$^{4}$。考虑到内支撑受力和锚索集中力在基坑开挖过程中是动态变化的,为简化计算,取内支撑设计轴力和锚索的锁定值。根据设计资料,内支撑设计轴力为1441 kN,锚索锁定值从上向下依次取值为300 kN, 350 kN, 350 kN, 400 kN, 400 kN, 450 kN, 500 kN, 500 kN, 500 kN和450 kN;临近建筑每层超载取20 kPa,基底总应力为140 kPa;围护桩桩身模量为38 GPa,临空段长度为23 m,嵌固段长度为11 m,围护桩的桩顶和桩端都采取自由边界条件来合理简化实际工程。

该区段现场设置了两个水平位移监测点LD-1和LD-2,现场实测与采用本文方法计算所得的围护结构位移分布曲线如图5所示。可以看出,围护结构实测最大位移为7.12 cm和7.81 cm,采用本文方法计算得到的围护结构最大位移为6.58 cm,两者较为接近,且随深度的变化趋势基本一致,验证了本文提出的半解析解答的合理性。

图5

图5   实测结果与计算结果比较图


为进一步考虑周边临近建筑对基坑支护结构受力变形的影响,在该工程案例计算模型的基础上,计算地下室基底总应力分别为140 kPa,280 kPa和420 kPa情况下的支护桩的水平位移和弯矩,如图6所示。为进一步将基坑内预留土横断面的尺寸对基坑支护结构受力产生的作用考虑在内,在本工程算例的前提下,保持预留土高度和坡度不变,分别计算了预留土上宽为4 m,6 m,8 m,10 m和12 m时支护桩的水平位移和弯矩,见图6(a)和图6(b)。可以看出,围护桩水平位移和弯矩均随着临近建筑地基基底应力的增加而增加,当临近建筑基底总应力由140 kPa增加至280 kPa时,桩身最大水平位移增加了约14.9%,其弯矩增加了约4.6%,说明临近建筑基底应力对围护结构受力变形的影响不容忽视。

图6

图6   临近建筑基底应力的影响


同时,在基底应力保持不变的情况下,考虑临近建筑地下室距围护结构距离的影响,计算距离$S$分别为3.8 m,7.6 m和11.4 m情况下围护桩的内力变形,如图7所示。可以看出,当临近建筑物与围护结构之间的距离不断增加时,围护桩的水平位移和弯矩值均呈逐渐减小的变化趋势,但距离对围护桩水平位移的影响更为显著。例如,当距离由3.8 m增加到7.6 m和11.4 m时,桩身最大水平位移减小了约10.5%和20%,且最大水平位移位置呈下移变化趋势;桩身最大弯矩减低了4.6%和11.6%。这主要是因为随着临近建筑距围护桩距离的增加,基底应力的影响范围逐渐下移,且坑内内支撑的轴力相对较大。

图7

图7   临近建筑距离的影响


4 结论

针对工程实践中经常遇到的临近建筑影响下基坑支护结构上受力分布的问题,进行基坑支护结构受力变形相关的理论计算研究:

(1)考虑临近建筑与支护结构间有限宽土体土压力、建筑基底压力等的影响,根据三参数弹性地基梁模型,得到支护结构合理的内力和变形的计算模型,并给出了较为合理的求解方法;

(2)依托典型工程案例,采取预测值和现场实测数据进行对比研究的方法,检验了本文半解析解答方法的合理性,可为类似建筑密集区基坑支护结构的设计计算提供借鉴;

(3)案例与算例分析表明,临近建筑对基坑支护结构的受力响应明显,不容忽视,在设计阶段宜充分考虑。

参考文献

Son M, Cording EJ.

Evaluation of building stiffness for building response analysis to excavation-induced ground movements

Journal of Geotechnical and Geoenvironmental Engineering, 2007,133(8):995-1002

DOI      URL     [本文引用: 1]

李志伟, 郑刚.

基坑开挖对邻近不同刚度建筑物影响的三维有限元分析

岩土力学, 2013,34(6):1807-1814

[本文引用: 1]

Li Zhiwei, Zheng Gang.

Finite element analysis of response of building with different stiffnesses adjacent to excavation

Rock and Soil Mechanics, 2013,34(6):1807-1814 (in Chinese)

[本文引用: 1]

张治国, 赵其华, 鲁明浩.

邻近深基坑开挖的历史保护建筑物沉降实测分析

土木工程学报, 2015,48(z2):137-142

[本文引用: 1]

Zhang Zhiguo, Zhao Qihua, Lu Minghao.

Analysis on settlement monitoring of historical protective buildings adjacent to deep foundation pit excavation

China Civil Engineering Journal, 2015,48(z2):137-142 (in Chinese)

[本文引用: 1]

张治国, 杨轩, 赵其华 .

浅基础框架建筑受基坑开挖影响简化分析

岩土工程学报, 2017,39(z2):224-227

[本文引用: 1]

Zhang Zhiguo, Yang Xuan, Zhao Qihua, et al.

Simplified analysis of frame buildings with shallow foundation induced by excavation of adjacent foundation pit

Chinese Journal of Geotechnical Engineering, 2017,39(z2):224-227 (in Chinese)

[本文引用: 1]

王卫东, 徐中华.

预估深基坑开挖对周边建筑物影响的简化分析方法

岩土工程学报, 2010,32(z1):32-38

[本文引用: 1]

Wang Weidong, Xu Zhonghua.

Simplified analysis method for evaluating excavation-induced damage of adjacent buildings

Chinese Journal of Geotechnical Engineering, 2010,32(z1):32-38 (in Chinese)

[本文引用: 1]

金亚兵, 彭虹洋.

深基坑边坡超载计算方法探讨及其工程应用

水文地质工程地质, 2000,2:4-7

[本文引用: 1]

Jin Yabing, Peng Hongyang.

Calculation method of overload for deep foundation pit slop and its application

Hydrogeology & Engineering Geology, 2000,2:4-7 (in Chinese)

[本文引用: 1]

张莲花, 赵毅.

超载作用对深基坑支护设计影响的研究

地质灾害与环境保护, 2009,20(1):66-69

[本文引用: 1]

Zhang Lianhua, Zhao Yi.

The analysis of the influence caused by overloading action on the design of deep foundation pit retaining structure

Journal of Geological Hazards and Environment Preservation, 2009,20(1):66-69 (in Chinese)

[本文引用: 1]

张浩, 石名磊, 郭院成 .

不平衡堆载作用下邻近结构桩的侧向受力机制

岩土工程学报, 2016,38(12):2226-2236

[本文引用: 2]

Zhang Hao, Shi Minglei, Guo Yuancheng, et al.

Lateral mechanical behaviors of structural piles adjacent to imbalanced surcharge loads

Chinese Journal of Geotechnical Engineering, 2016,38(12):2226-2236 (in Chinese)

[本文引用: 2]

李峰, 郭院成.

基坑工程有限土体主动土压力计算分析研究

建筑科学, 2008,24(1):15-18

[本文引用: 2]

Li Feng, Guo Yuancheng.

Analytical study on active soil pressure from finite soil body in construction pit

Building Science, 2008,24(1):15-18 (in Chinese)

[本文引用: 2]

张浩, 郭院成, 石名磊 .

坑内预留土作用下多支点支护结构的变形内力计算

岩土工程学报, 2018,40(1):162-168

[本文引用: 2]

Zhang Hao, Guo Yuancheng, Shi Minglei, et al.

Calculation of deformation and internal force of multi-pivot retaining structure considering influence of earth berm

Chinese Journal of Geotechnical Engineering, 2018,40(1):162-168 (in Chinese)

[本文引用: 2]

Viggiani C.

Ultimate lateral load on piles used to stabilize landslides

The 10th International Conference on Soil Mechanics and Foundation Engineering, Stockholm, Sweden, 1981

[本文引用: 1]

/