振动作用下碎石桩复合地基渗流理论分析与物理模拟

刘凯; 周燕国; 孙政波; 陈云敏

doi:10.11779/CJGE202009019

振动作用下碎石桩复合地基渗流理论分析与物理模拟 English Version

刘凯^{1, 2,},
周燕国^{1, 2, ,},
孙政波³,
陈云敏^{1, 2}

1.
浙江大学软弱土与环境土工教育部重点实验室,浙江　杭州　310058
2.
浙江大学岩土工程研究所,浙江　杭州　310058
3.
浙江省建筑设计研究院,浙江　杭州　310006

基金项目:

国家自然科学基金项目 51988101

国家自然科学基金项目 51578501

国家自然科学基金项目 51778573

国家自然科学基金项目 51978613

高等学校学科创新引智计划项目 B18047

详细信息

作者简介:
刘凯(1993—),男,博士研究生,主要从事土动力学与地震工程方面的研究工作。E-mail：11612027@zju.edu.cn

通讯作者:
周燕国, E-mail：qzking@zju.edu.cn

中图分类号: TU435;TU411
计量
- 文章访问数: 284
- HTML全文浏览量: 40
- PDF下载量: 179
出版历程
- 收稿日期: 2019-11-17
- 网络出版日期: 2022-12-07
- 刊出日期: 2020-08-31

Theoretical analysis and physical modeling of stone columns-improved ground under shaking loading

LIU Kai^{1, 2,},
ZHOU Yan-guo^{1, 2, ,},
SUN Zheng-bo³,
CHEN Yun-min^{1, 2}

1.
Key Laboratory of Soft Soils and Geoenvironmental Engineering of Ministry of Education, Hangzhou 310058, China
2.
Institute of Geotechnical Engineering, Zhejiang University, Hangzhou 310058, China
3.
Zhejiang Province Institute of Architectural Design and Research Institute, Hangzhou 310006, China

摘要

摘要: 碎石桩技术是常用的抗液化处理措施之一。在地震荷载作用时,由于桩体材料比桩间土渗透性更高,能够加快地基土体超静孔压的消散从而部分或完全消除地基的液化可能性。首先将单根碎石桩及其影响范围的地基土当做一个表征体元,将碎石桩影响边界处的孔压作为已知条件,基于已有解析解推导了碎石桩复合地基竖向和水平向渗流的计算表达式。并开展了一组碎石桩处理地基超重力振动台模型试验,用提出的碎石桩复合地基径竖向渗流解得到总渗流量,并基于地基土体的体积守恒关系得到计算沉降值,与试验结果吻合,从而验证了提出的碎石桩复合地基渗流计算公式的正确性。进一步基于所推导的渗流计算表达式,分析了本次超重力振动台模型试验碎石桩复合地基水平向和竖向渗流演化规律。
- 碎石桩复合地基 /
- 振动荷载 /
- 解析解 /
- 渗流排水 /
- 超重力模型试验
Abstract: The stone column technique is one of the liquefaction mitigation methods used in liquefiable deposite. Owing to the higher permeability of stone column materials than that of the surrounding soil, the stone columns may expedite the drainage process of the excess pore pressure under seismic loading so that to partly or completely eliminate the triggering of liquefaction. One stone column with the influenced surrounding soil is treated as a representative element, and based on the previous analytical solution to this problem, the analytical solutions for vertical and radical drainages are obtained by considering the excess pore water pressure on the boundary of single stone column-influenced area as the known condition. The centrifugal model tests on stone column-improved ground are carried out. The total drainage water in the centrifuge model improved by stone columns is calculated using the above analytical solution, and is further converted into the settlement based on the equivalent relationship between the seepage discharge and the volumetric change of the surrounding soil. The calculated settlement is found to be in good agreement with the measured one in experiment, which verifies the proposed solution. Furthermore, the evolution of horizontal and vertical drainages of stone column composite deposite in the model tests is analyzed by combining the derived analytical solutions for water drainage with the test results.
- stone column-improved ground /
- shaking loading /
- analytical solution /
- water drainage /
- hypergravity physical modeling

HTML全文

0. 引言

地震液化是国内外历次大地震中导致土工建构筑物失效破坏和造成严重经济损失的主要因素之一^[1-3]。地震液化导致的地面过大沉降和水平侧移是造成房屋、桥梁和地下生命管线等基础设施破坏的主要原因。当前国内外采用的可液化地基处理技术主要包括换填法、加密处理、固化技术（水泥、微生物固化等）、导排法（加速超静孔压消散）等,碎石桩法被认为是最经济和抗液化内涵最全面的技术手段,其抗液化效果在以往国内外历次大地震中得到检验^[4-5]。在碎石桩复合地基中,桩体材料渗透性比地基土更大,在地震过程中和震后可以作为排水通道,缩短渗流路径,改善排水条件。因此,可加速地震时场地超静孔隙水压力消散,使振动时孔压增长和消散同时发生,降低孔压峰值,提高处理地基的抗液化能力。

黄茂松等结合实际工程对振冲碎石桩加固饱和粉砂地基各施工过程的孔压规律进行了全面研究^[6],并建议了可应用于工程分析的的碎石桩排水效应简化分析方法^[7]。此外,国内外许多学者从模型试验和理论研究层面对碎石桩复合地基排水效应或固结过程展开研究。李立军等^[8]通过常重力振动台模型试验对比了碎石桩和水泥土桩加固液化地基的效果,发现碎石桩模型中地基土的超静孔压明显更小,由此得出碎石桩复合地基的排水效应对地基土抗液化贡献更为显著；Sasaki等^[9]将饱和砂箱固定在常重力振动台上研究砾石桩的排水性能和地基土超静孔压增长和消散特性,发现砾石排水桩可以迅速降低土体中的超静孔压水平,防止地基液化；Huang等^[10]通过常重力振动台试验发现,经碎石桩处理后的场地在振动过程中超静孔压累积速率明显降低,地基土体的整体刚度增强,地表沉降显著减小；Dashti等^[11]开展的超重力振动台试验表明,碎石桩加速了地震过程中地基土体三维排水过程,显著减小液化持时并进一步减小了震后残余沉降量；Badanagki等^[12]通过离心模型试验发现碎石桩的排水效应显著减小了倾斜场地的水平和竖向位移。

Seed等^[13]最早建立了考虑碎石桩水平向排水的微分方程,通过引入室内三轴不排水试验得到的孔压增长模型,得到考虑地震荷载的碎石桩复合地基的孔压控制方程,通过数值求解得到了不同碎石桩设计参数下超静孔压增长和消散规律并用于指导工程实践；徐志英^[14]建立了考虑地震过程中砾石排水桩径竖向排水的处理场地超静孔压演化偏微分方程,给出了超孔隙水压力的一般解析式。其后,国内学者在理论层面对碎石桩复合地基的堆载固结问题研究较多：Tang等^[15]建立了变荷载作用下考虑井阻效应和涂抹效应的径竖向组合渗流解析解；郭彪等^[16]推导了考虑桩体和土体径竖向渗流、上部荷载逐级施加、扰动区渗透系数线性变化的较全面的散体材料桩复合地基解析解；卢萌盟等^[17-18]给出了考虑桩体固结与土体渗透性抛物线变化和地基中附加应力沿深度线性分布的荷载效应的碎石桩复合地基固结解析解。

由此可见,模型试验主要集中在碎石桩排水效应和抗液化效果的宏观现象观测,尚未和理论研究相结合；而碎石桩理论研究多集中于对碎石桩复合场地堆载预压固结的孔压解析表达,而考虑地震荷载的研究成果多把建立的超静孔压解析式与和按照Seed简化方法得到的等效振次相联系,并未对碎石桩复合地基的排水性能进行深入研究。振动荷载作用下碎石桩复合地基的渗流与排水性能值得进一步研究。

本文在已有解析解的基础上,进一步推导了任意荷载作用下考虑径竖向渗流的碎石桩复合地基水平向和竖向渗流计算表达式,并开展了一组碎石桩处理地基超重力振动台模型试验,结合本文推导的径竖向渗流量表达式,发现计算得到的沉降量与模型试验实测沉降量吻合良好,验证了本文提出的碎石桩复合地基径竖向渗流计算表达式的正确性。

1. 微分方程的建立

地震引起的振动荷载作用于散体材料桩复合地基时,可认为各桩体和桩周土的受力性状相同,任选一受力变形复合体作为分析计算模型,如图1所示。本文做以下假设：①桩体和桩周土遵循等应变假设,即在同一深度处地基土体和桩体的竖向变形相等；②复合地基中水的渗流符合达西定律；③忽略桩体内的水平向渗流；④在任一深度z处,从土体流入桩体的水量等于桩体中向上水流的增量；⑤土体饱和、土颗粒和水不可压缩。

图 1 碎石桩处理地基计算模型示意图

Figure 1. Schematic diagram of numerical model for stone column-improved ground

下载: 全尺寸图片幻灯片

根据有效应力原理,任意时刻变荷载作用下土体的应力应变关系满足下式：

$ε_{v} = \frac{σ^{'}}{E_{s}} = \frac{q - {\bar{u}}_{s}}{E_{s}}$ ,

(1)

式中, $ε_{v}$ 为土体体变, $E_{s}$ 为土体压缩模量,q为任意荷载作用下产生的总应力,一般用经验孔压模型替代^[13], ${\bar{u}}_{s}$ 为土体任一深度的平均超静孔压。

式（1）两边对时间求偏导得

$\frac{\partial ε_{v}}{\partial t} = - \frac{1}{E_{s}} [\frac{\partial {\bar{u}}_{s}}{\partial t} - \frac{\partial q}{\partial t}] 。$

(2)

根据文献[15,19],对于径竖向渗流情况,由饱和土体变与固结排水量等量关系以及达西定律,可得

$- \frac{k_{sh}}{γ_{w}} \frac{1}{r} \frac{\partial}{\partial r} (r \frac{\partial u_{s}}{\partial r}) - \frac{k_{sv}}{γ_{w}} \frac{\partial^{2} {\bar{u}}_{s}}{\partial z^{2}} = \frac{\partial ε_{v}}{\partial t}$ ,

(3)

式中,k_sh和k_sv分别为地基土的水平向和竖向渗透系数, $γ_{w}$ 为水的重度。

根据假设④,且根据碎石桩与土体交界面处连续条件可得

$\frac{\partial^{2} u_{p}}{\partial z^{2}} = - \frac{2 k_{sh}}{r_{p} k_{p}} \frac{\partial u_{s}}{\partial r} |_{r = r_{p}},$

(4)

式中,u_p为任一深度处桩体内的超静孔压,r_p为桩体半径,k_p为桩体渗透系数。

任一深度地基土的平均超静孔压值定义参考已有研究^[15-17]：

${\bar{u}}_{s} = \int_{r_{p}}^{r_{e}} \frac{1}{π (r_{e}^{2} - r_{p}^{2})} \times 2 π r u_{s} d r$ ,

(5)

式中,r_e为碎石桩单桩影响半径。

式（2）～（5）为碎石桩复合地基微分控制方程。方程求解的初始条件和边界条件为

$\frac{\partial u_{s}}{\partial r} |_{r = r_{e}} = 0,$

(6)

$u_{p} |_{z = 0} = 0, {\bar{u}}_{s} |_{z = 0} = 0,$

(7)

$\frac{\partial u_{p}}{\partial z} |_{z = H} = 0, \frac{\partial {\bar{u}}_{s}}{\partial z} |_{z = H} = 0,$

(8)

${\bar{u}}_{s} |_{t = 0} = u_{0} (z) = 0$ ,

(9)

$u_{s} |_{r = r_{p}} = u_{p}$ 。

(10)

2. 微分方程的求解

参照Tang等^[15,19]对该类问题的求解过程和解的形式,可得如下 ${\bar{u}}_{s}$ 的解析表达式：

${\bar{u}}_{s} (z, t) = \int_{0}^{t} \frac{d q}{d ζ} \sum_{m = 1}^{\infty} \frac{2}{M} \sin \frac{M z}{H} e^{- A_{rv} (t - ζ)} d ζ,$

(11)

式中, $ζ$ 为任意荷载施加的时刻, $A_{rv} = A_{r} + A_{v}$ , $A_{v} = \frac{k_{sv} E_{s}}{γ_{w}} \cdot \frac{M^{2}}{H^{2}}, A_{r} = \frac{β}{1 + ψ^{2} {(\frac{H}{M})}^{2} - A ψ^{2}}, M = \frac{2 m + 1}{2} π$ （m=0,1,2,…）, $β =$ $2 E_{s} k_{h} / (γ_{w} r_{e}^{2} C_{n})$ , $ψ^{2} = 2 (N^{2} - 1) \cdot$ $k_{h} /$ $(k_{w} r_{e}^{2} C_{n})$ , $N = \frac{r_{e}}{r_{p}}$ 为井径比。

式（3）两边对r积分两次,并利用边界条件（6）和（10）可得

$u_{s} = \frac{γ_{w}}{2 k_{sh}} (r_{e}^{2} \ln \frac{r}{r_{p}} - \frac{r^{2} - r_{p}^{2}}{2}) (\frac{\partial ε_{v}}{\partial t} + \frac{k_{sv}}{γ_{w}} \frac{\partial^{2} {\bar{u}}_{s}}{\partial z^{2}}) + u_{p}$ 。

(12)

将式（11）和式（2）代入式（12）可得

$u_{s} = R (r) \cdot G (q, t, z) + u_{p}$ ,

(13)

式中,

$\begin{matrix} G (q, t, z) = \frac{γ_{w}}{2 k_{sh}} \cdot {\frac{1}{E_{s}} {\frac{\partial q}{\partial ζ} - \sum_{m = 1}^{\infty} \frac{2}{M} \sin \frac{M z}{H} \cdot \\ \int_{0}^{t} [\frac{d}{d t} (\frac{d q}{d ζ}) e^{- A_{rv} (t - ζ)} - A_{rv} \frac{d q}{d ζ} e^{- A_{rv} (t - ζ)}] \cdot d ζ} - \\ \frac{k_{sv}}{γ_{w}} \int_{0}^{t} \frac{d q}{d ζ} \sum_{m = 1}^{\infty} \frac{2 M}{H^{2}} \sin \frac{M z}{H} e^{- A_{rv} (t - ζ)} d ζ}, \\ R (r) = (r_{e}^{2} \ln \frac{r}{r_{p}} - \frac{r^{2} - r_{p}^{2}}{2}) \end{matrix}$

。

由式（13）可见,在碎石桩桩径影响范围内（影响半径为r_e）,土体任一深度的超静孔压u_s可以写成R与G的乘积与该深度桩体内孔压值u_p的和。其中,函数R仅仅是关于半径r的函数,而函数G与z,t和q相关。以上说明土体中某一深度距离碎石桩中心距为r (r_p ≤r≤r_e)的超静孔压u_s分布形式仅与r有关,不同外荷载仅是通过函数G和u_p影响超静孔压的大小。

在dt时间内,深度z处桩周土径向（水平向）通过碎石桩排水的流量为

$Δ {(q_{sc})}_{z} = k_{sh} i A = \frac{k_{sh}}{γ_{w}} \frac{\partial u_{s}}{\partial r} |_{r = r_{p}} 2 π r_{p} d z d t$ 。

(14)

进一步考虑式（13）,式（14）可化为

$Δ {(q_{sc})}_{z} = \frac{k_{sh}}{γ_{w}} G (z, t, q) (\frac{r_{e}^{2}}{r_{p}} - r_{p}) 2 π r_{p} d z d t$ 。

(15)

式（15）对时间积分即可得到超静孔压累积和消散过程中桩周土深度z处通过碎石桩的水平向排水量。

假设碎石桩处理地基中单桩影响范围边界处的超静孔压为 $u_{e} |_{r = r_{e} = u_{e} (z, t)}$ ,且假设在地基土体中超静孔压从产生到消散过程中满足u_e(z,t)>>u_p,则式（13）可化为

$u_{s} \approx G (q, z, t) (r_{e}^{2} \ln \frac{r}{r_{p}} - \frac{r^{2} - r_{p}^{2}}{2}),$

(16)

且,

$u_{e} (z, t) \approx G (q, z, t) (r_{e}^{2} \ln \frac{r_{e}}{r_{p}} - \frac{r_{e}^{2} - r_{p}^{2}}{2})$ 。

(17)

结合式（15）,（16）,（17）可得

$Δ {(q_{sc})}_{z} = \frac{k_{sh}}{γ_{w}} \frac{u_{e} (z, t)}{(r_{e}^{2} \ln \frac{r_{e}}{r_{p}} - \frac{r_{e}^{2} - r_{p}^{2}}{2})} (\frac{r_{e}^{2}}{r_{p}} - r_{p}) 2 π r_{p} d z d t 。$

(18)

式（18）即为在dt时间内深度z处地基土通过桩体水平向排水的流量计算式。

将式（16）代入式（5）,并结合式（17）,可得某一深度土体的平均孔压为

${\bar{u}}_{s} \approx \frac{\frac{r_{e}^{4}}{r_{e}^{2} - r_{p}^{2}} \ln \frac{r_{e}}{r_{p}} - \frac{3}{4} r_{e}^{2} + \frac{1}{4} r_{p}^{2}}{(r_{e}^{2} \ln \frac{r_{e}}{r_{p}} - \frac{r_{e}^{2} - r_{p}^{2}}{2})} u_{e} (z, t) 。$

(19)

由此得到在dt时间内,在深度z_i和z_i+1处地基土土竖向渗流流量公式为

$\begin{array}{l} Δ {(q_{s})}_{z = z_{i} \sim z_{i + 1}} = \frac{π k_{sv} r_{e}^{2} (1 - A_{r})}{γ_{w} (z_{i + 1} - z_{i}) (r_{e}^{2} \ln \frac{r_{e}}{r_{p}} - \frac{r_{e}^{2} - r_{p}^{2}}{2})} \cdot \\ [\frac{r_{e}^{4}}{r_{e}^{2} - r_{p}^{2}} \ln \frac{r_{e}}{r_{p}} - \frac{3}{4} r_{e}^{2} + \frac{1}{4} r_{p}^{2}] [u_{e} (z_{i}, t) - u_{e} (z_{i - 1}, t)] d t 。 \end{array}$

(20)

3. 超重力模型试验验证

3.1 试验设备和试验材料

利用浙江大学ZJU-400超重力离心机振动台开展了一组碎石桩处理地基的模型试验,以验证本文所提出计算公式的正确性。超重力离心机通过高速旋转使得模型土体恢复原型自重应力,然后通过机载振动台在模型底部实现地震动输入,从而模拟原型场地地震响应和液化灾变过程^[20]。本次试验在50g离心加速度下进行,选用粘滞系数为水的50倍的甲基硅油保证模型动力时间与渗流时间相似比尺一致。

采用层状剪切模型箱模拟自由场地边界条件,内部尺寸（长×宽×高）为730 mm×330 mm×425 mm。为了合理模拟现场地基土的渗透特性,试验地基土体选用福建细砂掺10%的粉土（钱塘江粉土）,碎石桩采用粗粒径的福建砂,两种材料的基本物性参数和颗分曲线分别见表1和图2。由表1可见,碎石桩材料渗透系数约为地基土体的800倍。Seed等^[13]建议对于碎石桩处理地基,碎石桩材料的渗透系数一般为处理地基土体的两个数量级以上。因此,模型试验设计地基土体和碎石桩材料的渗透特性符合工程实际情况。

表 1 试验材料基本物理参数

Table 1. Physical properties of test materials

指标	G_s	ρ_max/(g·cm^-3)	ρ_min/(g·cm^-3)	ρ_min/(g·cm^-3)	k/(m/s, 27℃)
福建砂掺10%粉土	2.647	1.765	1.386	1.695	1.864×10^-5
福建粗砂	2.644	1.713	1.489	1.605	1.49×10^-2

下载: 导出CSV

| 显示表格

图 2 试验材料级配曲线

Figure 2. Grain-size distribution curves of test materials

下载: 全尺寸图片幻灯片

3.2 模型布置

模型试验传感器布置如图3所示。模型土层厚度为400 mm。在振动台台面布置一个三向加速度计记录模型底层振动输入,在中轴线左侧50 mm沿深度布置6个水平加速度计,在中轴线右侧50 mm沿深度布置6个孔压计,分别在距离中轴线左右两侧各150 mm对称布置土压力计和弯曲元,在模型地表布置1个激光位移传感器监测模型的表面沉降。

图 3 离心模型试验布置图

Figure 3. Schematic configurations of centrifuge model

下载: 全尺寸图片幻灯片

在模型地基中设计了3×7共21根碎石桩,桩体贯穿整个液化土层。碎石桩的桩径d采用和Adalier等^[21]的离心模型试验原型尺寸一致,为1.5 m,对应模型尺寸为3 cm,桩间设计为10 cm,计算得到碎石桩处理地基的置换率为7%。

考虑到前述公式推导假设地基土的渗透系数为常数,本文选用频率1 Hz、峰值加速度小于0.05g的等幅正弦波作为台面输入（如图4所示）,以避免过大的振动引起过高的超静孔压甚至触发液化,导致地基土体渗透系数发生显著变化^[22]。试验共进行两次振动,振动1结束后,所有通道稳定15 min以上等待振动1产生的超静孔压完全消散,再开始振动2。

图 4 两次台面振动输入

Figure 4. Input motions at base of two shaking events

下载: 全尺寸图片幻灯片

3.3 典型试验结果

图5,6分别给出了两次振动的模型地表沉降和模型内部超静孔压时程曲线。由图5可见,振动初期地表沉降出现“不降反升”的异常现象。有两个可能的原因：首先,在振动初期的近似不排水循环剪切过程中,由于模型土体需要满足变形协调导致模型地表出现上下起伏的动态变形,但平均（残余）位移是向下发展的；此外,振动造成超静孔压上升,模型土体骨架有瞬态卸载、有限膨胀的响应,也会造成模型地表短时的微小隆起。

图 5 模型地表沉降时程曲线

Figure 5. Time histories and prediction of surface settlement

下载: 全尺寸图片幻灯片

图 6 超静孔压时程曲线

Figure 6. Time histories of excess pore water pressure

下载: 全尺寸图片幻灯片

由图6可见,随着振动的发生,模型不同深度土层的超静孔压随着振动迅速累积,并在振动结束时达到峰值。振动1中,不同深度超静孔压峰值在6～13 kPa；振动2中,不同深度的超静孔压峰值在16～22 kPa。振动停止后,模型深层的超静孔压消散比浅层更快,说明振后模型内孔隙水向上渗流。

对照图5,6,可以发现地表沉降发展与超静孔压消散存在对应关系。地表沉降与超静孔压累积主要发生在振动过程,振后随着超静孔压消散沉降继续发生,但振后固结沉降比例相对较小,这与以往离心机振动台试验结果一致^[23]。主要原因在于模型地基地表是自由排水边界,振动过程中就发生了显著的孔隙水向上消散。由图6可见,振动1情况模型全部深度的超静孔压比（r_u=u/ ${σ^{'}}_{v0}$ ,u为任意时刻的超静孔压, ${σ^{'}}_{v0}$ 为初始上覆有效应力）都小于0.4,振动对模型地基造成的损伤较小,对应模型地表沉降较小；振动2情况模型浅层r_u接近0.8,其余深度r_u小于0.3,模型地表沉降大于振动1引起的地表沉降。

已有研究^[24]表明,砂土液化时土体渗透系数比未液化时增大数量级的倍数。Shahir等^[25]提出的随超静孔压比变化的砂土变渗透系数表达式如下：

$\frac{k}{k_{i}} = {\begin{cases} 1 + (α - 1) r_{u}^{β_{1}}, r_{u} < 1 (孔压累积阶段) \\ α, r_{u} = 1 \\ 1 + (α - 1) r_{u}^{β_{2}}, r_{u} < 1 (孔压消散阶段) \end{cases}$ 。

(21)

式中, $k_{i}$ 为未液化前渗透系数初始值,α为液化时渗透系数增大的倍数,β₁和β₂分别为孔压累积阶段和消散阶段的常数,与土体材料有关。

结合图6,对式（21）中参数取值α=10,β₁=2和β₂=10,可得当r_u分别等于0.3,0.4和0.8时,对应的渗透系数增大倍数分别为1.8,2.4和6.7。由图6可见,r_u处于峰值的时间占整个超静孔压从产生到消散的时间比例不足5%。因此,两次振动引起的超静孔压比不足以引起整个模型深度土体渗透系数的显著变化,本文采用常渗透系数假设基本合理。

3.4 渗流公式正确性验证

由图3可见,孔压计布置在4根碎石桩围成的正方形的中点,可近似认为孔压计测值为碎石桩桩径影响范围边界r_e上的值。为验证式（18）和（20）的正确性,利用地基地表沉降量对振动过程中场地超静孔压从产生到消散时间T内的总排水量Q进行校核。假设土颗粒和饱和流体不可压缩,则总排水量等于地基沉降引起的体变值。注意到本文计算模型是基于单桩影响范围的场地沉降,排水量和地表沉降的计算都是基于单桩影响范围,即得到的总排水量Q需除以单桩影响范围的地基土面积才得到对应的沉降量。碎石桩水平排水总量Q_L为地基各深度水平向排水量之和,而竖向排水总量Q_V只考虑最浅层竖向排水量,即

$Q = Q_{L} + Q_{V},$

(22)

$Q_{L} = \int_{0}^{H} \int_{0}^{T} Δ {(q_{s c})}_{z} d t d z,$

(23)

$Q_{V} = \int_{0}^{T} Δ {(q_{s})}_{top_layer} d t$ 。

(24)

图7为根据两次振动超静孔压时程,通过式（22）计算沉降值和激光位移计实测沉降值。由图7可见,振动1情况下本文计算值与实测值接近,而振动2情况下的计算值稍小于实测值。由前述可知,振动2作用下模型浅层土体的渗透系数增大了6.7倍,本文“振动过程中地基土体渗透系数不变”的假设略微低估了模型浅层桩体的渗流量。

图 7 模型地表沉降实测值和计算值对比

Figure 7. Comparison between measured and predicted surface settlements

下载: 全尺寸图片幻灯片

图8给出了根据式（22）反算得到的沉降时程预测曲线。由图8可见,振动阶段实测沉降比计算值显著偏大；对于振后渗流固结过程,实测值与计算值接近。由于本文提出的渗流量计算公式仅考虑渗流固结效应,对模型在振动过程中可能出现的诸如渗流速率效应和液化时渗透系数变化等^[26-27]无法合理描述,可以认为本文公式适用于预测土体振后沉降规律。

图 8 模型地表沉降时程实测值和计算值对比

Figure 8. Time histories of measured and predicted surface settlement

下载: 全尺寸图片幻灯片

3.5 碎石桩复合地基振动排水特性分析

结合本次离心模型试验土体参数、碎石桩布置和不同深度超静孔压时程,根据式（18）和式（20）得到如图9所示的两次振动下碎石桩复合地基竖向和水平渗流流量随时间变化关系。

图 9 碎石桩复合地基渗流时程

Figure 9. Time histories of discharge water in stone column-improved ground

下载: 全尺寸图片幻灯片

由图9可见,水平向和竖向渗流排水在振动一开始就有,表明一旦产生超静孔压,在水力梯度作用下地基就开始进行竖向和水平向排水。在振动阶段地基土体不同深度通过碎石排水的渗流速率逐渐加快,振动停止后渗流速率基本保持稳定,随着超静孔压消散渗流速率逐渐降低。以上径竖向渗流随时间发展规律与图5所示的沉降时程一致。另外,两次振动下地基土体通过碎石桩的水平向渗流量随深度的增加而逐渐增大,不同深度土体通过碎石桩水平向排水量与相应的超静孔压水平正相关。两次振动下竖向渗流量都在浅层最大,且比较接近浅层水平向渗流量,但远小于深层水平向排水量。

4. 结论

本文针对振动荷载作用下碎石桩复合地基的渗流与排水性能问题,基于已有的解析解,进一步推导了任意荷载作用下碎石桩复合地基水平向和竖向渗流计算表达式,并开展了一组碎石桩处理地基超重力振动台模型试验进行验证,得到如下3点结论。

（1）将单根碎石桩及其影响范围的地基土作为一个表征体元来研究,以碎石桩单桩影响边界处的孔压值为已知条件,得到了任意荷载作用下碎石桩复合地基竖向和水平向渗流量计算公式。

（2）假设土颗粒和水不可压缩,将计算得到的模型地基排水量换算为地表沉降值,发现最终沉降的计算值与模型试验实测值一致,验证了本文推导的碎石桩复合地基渗流量计算公式的正确性。

（3）将模型试验测得的超静孔压与本文推导的渗流量计算公式相结合,给出了本次模型试验中不同深度水平向和竖向渗流时程,与实测沉降时程基本吻合。碎石桩复合地基的水平向渗流起主要作用,竖向渗流贡献较小；受下大上小的孔压梯度场的驱动,深层的水平向渗流量明显大于浅层。

在实际工程中,如果能基于地震动强度指标等获得易液化场地的超静孔压水平,则可以利用本文提出的计算公式估算处理地基的震后排水量和地表沉降值,为工程设计和抗震性能调控提供依据。

图 1 碎石桩处理地基计算模型示意图

Figure 1. Schematic diagram of numerical model for stone column-improved ground

下载: 全尺寸图片幻灯片

图 2 试验材料级配曲线

Figure 2. Grain-size distribution curves of test materials

下载: 全尺寸图片幻灯片

图 3 离心模型试验布置图

Figure 3. Schematic configurations of centrifuge model

下载: 全尺寸图片幻灯片

图 4 两次台面振动输入

Figure 4. Input motions at base of two shaking events

下载: 全尺寸图片幻灯片

图 5 模型地表沉降时程曲线

Figure 5. Time histories and prediction of surface settlement

下载: 全尺寸图片幻灯片

图 6 超静孔压时程曲线

Figure 6. Time histories of excess pore water pressure

下载: 全尺寸图片幻灯片

图 7 模型地表沉降实测值和计算值对比

Figure 7. Comparison between measured and predicted surface settlements

下载: 全尺寸图片幻灯片

图 8 模型地表沉降时程实测值和计算值对比

Figure 8. Time histories of measured and predicted surface settlement

下载: 全尺寸图片幻灯片

图 9 碎石桩复合地基渗流时程

Figure 9. Time histories of discharge water in stone column-improved ground

下载: 全尺寸图片幻灯片

表 1 试验材料基本物理参数

Table 1 Physical properties of test materials

指标 G_s ρ_max/(g·cm^-3) ρ_min/(g·cm^-3) ρ_min/(g·cm^-3) k/(m/s, 27℃)

福建砂掺10%粉土 2.647 1.765 1.386 1.695 1.864×10^-5
福建粗砂 2.644 1.713 1.489 1.605 1.49×10^-2

下载: 导出CSV

参考文献(27)

[1]	黄雨, 于淼, BHATTACHARYA S. 2011年日本东北地区太平洋近海地震地基液化灾害综述[J]. 岩工工程学报, 2013, 35(5): 834-840. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201305005.htm HUANG Yu, YU Miao, BHATTACHARYA S. Review on liquefaction induced damages of soils and foundations during the 2011 off the Pacific Coast of Tohoku Earthquake(Japan)[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(5): 834-840. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201305005.htm
[2]	ZHOU Yan-guo, CHEN Yun-min, LING Dao-sheng. Shear wave velocity-based liquefaction evaluation in the great Wenchuan Earthquake: a preliminary case study[J]. Earthquake Engineering and Engineering Vibration, 2009, 8(2): 231-239. doi: 10.1007/s11803-009-9049-9
[3]	周燕国, 李永刚, 丁海军, 等. 砂土液化后再固结体变规律表征与离心模型试验验证[J]. 岩土工程学报, 2014, 36(10): 1838-1845. doi: 10.11779/CJGE201410011 ZHOU Yan-guo, LI Yong-gang, DING Hai-jun, et al. Characterization of reconsolidation volumetric strain of liquefied sand and validation by centrifuge model tests[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(10): 1838-1845. (in Chinese) doi: 10.11779/CJGE201410011
[4]	BHATTACHARYA S, HYODO M, GODA K, et al. Liquefaction of soil in the Tokyo Bay area from the 2011 Tohoku (Japan) earthquake[J]. Soil Dynamics Earthquake Engineering, 2011, 31(11): 1618-1628. doi: 10.1016/j.soildyn.2011.06.006
[5]	陈国兴, 顾小锋, 常向东, 等. 1989～2011期间8次强地震中抗液化地基处理成功案例的回顾与启示[J]. 岩土力学, 2015, 36(4): 1102-1118. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201504029.htm CHEN Guo-xing, GU Xiao-feng, CHANG Xiang-dong, et al. Review and implication of successful ground improvement cases about mitigating soil liquefaction induced by 8 strong earthquakes from 1989 to 2011[J]. Rock and Soil Mechanics, 2015, 36(4): 1102-1118. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201504029.htm
[6]	黄茂松, 陈云敏, 吴世明. 振冲碎石桩加固饱和粉砂粉土地基试验研究[J]. 岩土工程学报, 1992, 14(6): 69-73. doi: 10.3321/j.issn:1000-4548.1992.06.008 HUANG Mao-song, CHEN Yun-min, WU Shi-ming. Strengthening of saturated silty soils by vibro replacement stone columns[J]. Chinese Journal of Geotechnical Engineering, 1992, 14(6): 69-73. (in Chinese) doi: 10.3321/j.issn:1000-4548.1992.06.008
[7]	黄茂松, 吴世明, 赵竹占. 振动挤密砂桩与振冲碎石桩抗液化分析[J]. 浙江大学学报(工学版), 1992, 26(2): 165-171. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC199202004.htm HUANG Mao-song, WU Shi-ming, ZHAO Zhu-zhan. Liquefaction analysis of subsoil strengthed by vibroflotation stone pile method and vibro-compaction sand pile method[J]. Journal of Zhejiang University (Engineering Science), 1992, 26(2): 165-171. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC199202004.htm
[8]	李立军, 牛琪瑛, 梁仁旺, 等. 碎石排水桩抗地震液化的试验研究与理论分析[J]. 工程力学. 2010(增刊2): 226-230. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX2010S2028.htm LI Li-jun, NIU Qi-yinng, LIANG Ren-wang, et al. Experimental study and theoretical analysis on drainage gravel pile resistance of earthquake liquefaction[J]. Engineering Mechanics, 2010(S2): 226-230. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX2010S2028.htm
[9]	SASAKI Y, TANIGUCHI E. Large scale shaking table tests on the effectiveness of gravel drains for liquefiable sand deposits[C]//Proceedings of the Conference on Soil Dynamics and Earthquake Engineering, 1982, Southampton: 843-857.
[10]	HUANG Chun-xia, SUI Zhi-long, WANG Lei, et al. Mitigation of soil liquefaction using stone columns: an experimental investigation[J]. Marine Georesource and Geotechnology, 2016, 34: 244-251.
[11]	DASHTI S, BRAY J D, PESTANA J M, et al. Centrifuge testing to evaluate and mitigate liquefaction-induced building settlement mechanisms[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2010, 136(7): 918-929. doi: 10.1061/(ASCE)GT.1943-5606.0000306
[12]	BADANAGKI M, DASHTI S, KIRKWOOD P. Influence of dense granular columns on the performance of level and gently sloping liquefiable sites[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2018, 144: 0401806.
[13]	SEED H B, MARTIN P P, LYSMER J. Generation and Dissipation of Pore Pressure Water Pressures During Liquefaction[R]. California: Earthquake Engineering Research Center, University of California, 1975.
[14]	徐志英. 用砾石排水桩抗地震液化的砂基孔压计算[J]. 地震工程与工程振动, 1992, 12(4): 88-92. https://www.cnki.com.cn/Article/CJFDTOTAL-DGGC199204011.htm XU Zhi-ying. Calculation of pore pressure of liquefiable sand deposits using gravel drains[J]. Earthquake Engineering and Engineering Vibration, 1992, 12(4): 88-92. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DGGC199204011.htm
[15]	TANG Xiao-wu, KATSUTADA O K. Consolidation by verticaldrains under time-dependent loading[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2000, 24(9): 739-751.
[16]	郭彪, 龚晓南, 李亚军. 考虑桩体径竖向渗流的碎石桩复合地基固结解析解[J]. 岩土工程学报, 2017, 39(8): 1485-1492. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201708020.htm GUO Biao, GONG Xiao-nan, LI Ya-jun. Analytical solution for consolidation of stone column-reinforced foundations[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(8): 1485-1492. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201708020.htm
[17]	卢萌盟, 谢康和, 刘干斌, 等. 考虑桩体固结和土体渗透性抛物线分布固结解[J]. 岩土工程学报, 2009, 31(7): 1069-1074. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200907017.htm LU Meng-meng, XIE Kang-he, LIU Gan-bin, et al. Solution for consolidation considering column consolidation and parabolic distribution of soil permeability[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(7): 1069-1074. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200907017.htm
[18]	卢萌盟, 谢康和, 张玉国, 等. 考虑施工扰动和荷载效应的复合地基固结解[J]. 岩土工程学报, 2008, 30(4): 549-554. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200804016.htm LU Meng-meng, XIE Kang-he, ZHANG Yu-guo, et al. Theoretical solution for consolidation of composite foundations considering construction disturbance and loading effect[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(4): 549-554. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200804016.htm
[19]	唐晓武, 柳江南, 杨晓秋, 等. 开孔管桩动孔压消散特性的理论研究[J]. 岩土力学, 2019, 40(9): 3335-3343. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201909006.htm TANG Xiao-wu, LIU Jiang-nan, YANG Xiao-qiu, et al. Theoretical study on dynamic pore water pressure dissipation characteristics of the open-hole pipe pile[J]. Rock and Soil Mechanics, 2019, 40(9): 3335-3343. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201909006.htm
[20]	陈云敏, 韩超, 凌道盛, 等. ZJU400离心机研制及其振动台性能评价[J]. 岩土工程学报, 2011, 33(12): 1887-1894. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201112013.htm CHEN Yun-min, HAN Chao, LING Dao-sheng, et al. Development of geotechnical centrifuge ZJU400 and performance assessment of its shaking table system[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(12): 1887-1894. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201112013.htm
[21]	ADALIER K, ELGAMAL A, MENESES J, et al. Stone columns as liquefaction counter-measure in non-plastic silty soils[J]. Soil Dynamics and Earthquake Engineering, 2003, 23(7): 571-584. http://www.sciencedirect.com/science/article/pii/S0267726103000708
[22]	SU D, LI X S, XING F. Estimation of the apparent permeability in the dynamic centrifuge tests[J]. Geotechnical Testing Journal, 2009, 32(1): 22-30.
[23]	周燕国, 谭晓明, 梁甜, 等. 利用地震动强度指标评价场地液化的离心模型试验研究[J]. 岩土力学, 2017, 38(7): 1869-1877. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201707004.htm ZHOU Yan-guo, TAN Xiao-ming, LIANG Tian, et al. Evaluation of soil liquefaction by ground motion intensity index by centrifuge model test[J]. Rock and Soil Mechanics, 2017, 38(7): 1869-1877. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201707004.htm
[24]	SHAHIR H, MOHAMMADI-HAJI B, GHASSEMI A. Employing a variable permeability model in numerical simulation of saturated sand behavior under earthquake loading[J]. Computers and Geotechnics, 2014, 55(1): 211-223.
[25]	SHAHIR H, PAK A, TAIEBAT M, JEREMIC′ B. Evaluation of variation of permeability in liquefiable soil under earthquake loading[J]. Computers and Geotechnics, 2012, 40(3): 74-88.
[26]	ZHANG J M, WANG G. Large post-liquefaction deformation of sand, part I: physical mechanism, constitutive description and numerical algorithm[J]. Acta Geotechnica, 2012, 7(2): 69-113.
[27]	ELGAMAL A, YANG Z, PARRA E, et al. Modeling of cyclic mobility in saturated cohesionless soils[J]. International Journal of Plasticity, 2003, 19(6): 883-905.

施引文献(20)

期刊类型引用(3)

1.	邵帅，赵梓君，邵生俊，田珮琦，刘小康，张彬. 复杂应力作用的黄土动剪切破坏强度研究. 岩土工程学报. 2025(03): 559-568 . 本站查看
2.	高常辉，刘松玉，杜广印，庄仲旬，杨泳，何欢. 气动振杆密实法加固湿陷性黄土的模型试验研究. 岩土工程学报. 2024(02): 325-334 . 本站查看
3.	邱毅清，郭志军，金鑫，司胜洁，郭福贵. 土壤振动切削减阻技术研究现状分析. 农机化研究. 2024(08): 8-18 . 百度学术