近河谷地铁车站振动台模型试验研究

禹海涛; 王治坤; 张中杰; 宋毅

doi:10.11779/CJGE20230379

近河谷地铁车站振动台模型试验研究 English Version

1.
同济大学土木工程防灾国家重点实验室, 上海 200092
2.
同济大学土木工程学院, 上海 200092
3.
上海市城市建设设计研究总院(集团)有限公司, 上海 200125

基金项目:

国家自然科学基金项目 42177134

中央高校基本科研业务费专项资金项目

详细信息

作者简介:
禹海涛(1983—)，男，工学博士，教授，博士生导师，主要研究方向为地下工程抗震防灾减灾。E-mail: yuhaitao@tongji.edu.cn

中图分类号: TU924
计量
- 文章访问数: 440
- HTML全文浏览量: 62
- PDF下载量: 118
出版历程
- 收稿日期: 2023-05-03
- 网络出版日期: 2024-05-05
- 刊出日期: 2024-08-31

Shaking table tests on a near-valley subway station

1.
State Key Laboratory for Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China
2.
College of Civil Engineering, Tongji University, Shanghai 200092, China
3.
Shanghai Urban Construction Design & Research Institute, Shanghai 200125, China

摘要

摘要: 近年来沿海沿江城市建设涌现出大量紧邻河谷地形修建的地铁车站，由于河谷地形具有地震动放大效应，会对近河谷结构抗震安全性构成威胁，但目前关于河谷地形效应对近邻结构的影响机制尚不清晰。提出了近河谷地铁车站1g振动台模型试验方法，分别设计了3组深宽比为1/8，1/4和1/3的圆弧形河谷地形，开展了不同地形、不同地震动输入条件下的近河谷地铁车站振动台模型试验研究，并以平整场地条件（无地形效应）下的地铁车站振动台试验工况为参照，对比分析了河谷地形效应对车站结构地震响应的影响规律。试验结果表明：与平整场地车站结构响应相比，在水平横向地震动作用下近河谷车站结构的动力响应较大，尤其是近河谷车站顶板处的加速度放大效应高达1.27倍；受河谷地形效应影响，车站中柱弯矩显著增加，而且邻近河谷侧的车站侧墙弯矩放大更为明显，即与平整场地相比放大效应可达2.3倍，同时该侧墙顶部的动土压力峰值也显著增强；随着河谷地形深宽比的增加，地形引起的局部地震动放大效应逐渐增加，导致近河谷车站结构的地震响应更为显著。研究成果可为近河谷地形地铁车站结构的抗震设计和安全性评估提供科学依据。
- 地铁车站 /
- 地震响应 /
- 振动台试验 /
- 河谷地形效应
Abstract: A large number of subway stations close to valley topography have been built in coastal and riverine cities in recent years, and the seismic safety of near-valley structures is threatened due to the seismic amplification effects of valley topography during earthquakes. However, the influence mechanism of valley topography on the adjacent structures has not been well understood. In this study, a 1g shaking table test model method for a near-valley subway station is proposed. Three sets of circular valley topography with depth-to-width ratios of 1/8, 1/4 and 1/3 are designed, the tests on a near-valley subway station under different topography and ground motion input conditions are performed, and a subway station test without valley site is used as a benchmark to investigate the seismic response of the valley topography on the adjacent structures. The results show that the dynamic response of the near-valley structures is greater under horizontal transverse motion compared to the response of the station structures without valley, particularly the acceleration amplification effects are up to 1.27 times at the structural top slab. The bending moment of the structural column increases due to the valley topography, and the amplification of the bending moment in the sidewall adjacent to the valley is more pronounced, i.e., up to 2.3 times compared to the tests without valley site. Additionally, the peak dynamic earth pressure at the top of the side wall is higher. The local seismic amplification effects of valley are more pronounced as the depth-to-width ratio of the valley topography increases, which leads to a greater seismic response of the near-valley station structures. The research results can provide a scientific basis for the seismic design and safety assessment of station structures around valley topography.
- subway station /
- seismic response /
- shaking table test /
- valley topographic effect

HTML全文

0. 引言

自1960年Henri Vidal^[1]首次提出了“加筋土”的概念和设计理论以来,土工合成材料广泛应用于岩土工程领域^[1]。在边坡中铺设土工合成材料能够优化土体应力分布,提高边坡的稳定性和承载力^[2-4]。有别于二维土工合成材料如土工格栅、土工布。张孟喜^[5]首先提出H-V立体筋的概念。作为一种新型加筋方式,H-V立体筋在筋–土相互作用上展现出更好的优势,除了具有水平筋作用外,竖筋能约束土体的变形并在竖筋间形成“加强土楔体”或“加强土柱”,增强土体的抗拉强度和抗剪强度^[6],具有较好的应用前景。

常用的加筋边坡稳定性分析方法主要有极限平衡法、有限元法等。改进最小势能法作为边坡稳定性分析方法之一,研究成果表明,相对极限平衡法,该方法计算方便快捷,具有不需要条块划分、不需要反复迭代求解等优点^[7-9]。目前,改进最小势能法主要应用于有加固措施如锚杆、锚索、抗滑桩的边坡稳定性分析中。邹治来^[10]考虑土工格栅与土体接触面之间摩擦势能和土工格栅弹性拉伸势能,建立土工格栅加筋边坡的最小势能稳定性分析方法。但是分析过程中,由于土工格栅与土体之间摩擦阻力计算过程复杂,参数众多,需要进行拉拔试验才能确定参数值,实现有很大难度。采用最小势能法分析了普通土工格栅加筋边坡稳定性,李铀等^[11]、孙加平等^[12]提出了考虑滑床剪切势能影响的近似计算模型,改进最小势能对锚杆（索）加固边坡稳定分析方法。温树杰等^[13]分析桩体变形特征,考虑抗滑桩与滑体相互接触处的弹性压缩势能、挠曲势能提出基于最小势能法加固边坡稳定性分析方法。未检索到最小势能法应用于H-V立体筋加筋边坡稳定性分析。

H-V立体筋边坡稳定性分析主要采用条分法,张孟喜等^[14]采用瑞典条分法对立体加筋边坡的稳定性进行探讨,在计算过程中未考虑条间力的作用,容易导致计算结果偏小。马学宁等^[15]建立水平条间力的假设条件,推导出水平条分瑞典法和简布法分析H-V立体加筋边坡的稳定性,该方法需要在水平筋材处于土条的中线上的条件下进行,且水平条分简布法需要进行反复迭代,计算量较大。其它对于H-V立体筋边坡稳定性分析文献尚不多见。

因此,本文基于改进最小势能法,对H-V立体筋作用的基本原理和受力特性进行分析,构建出水平筋和竖筋的弹性势能方程,并对加筋边坡单元体虚位移方向的剪切势能进行推导,通过计算边坡最危险滑裂面系统势能最小时的虚位移,求解出边坡安全系数,并通过算例对该方法的适用性和可行性进行验证。

本文利用改进最小势能法从能量的角度对H-V立体筋边坡的稳定性展开分析,在系统中涉及的能量种类主要有滑体储存的弹性势能,H-V筋因水平筋线弹性变形产生弹性势能、筋土之间因摩擦储存摩擦势能、竖筋产生侧向变形储存弹性势能和滑面剪切势能,概念清晰,可直接计算加筋边坡稳定性分析结果,相对于计算复杂的有限元法和需反复迭代、需对超静定问题进行多次假定的极限平衡法更有实用价值,为H-V立体筋边坡稳定性分析提供新的思路。

1. 滑体储存的弹性势能

如图1所示的H-V立体筋均质边坡,设潜在滑裂面为AC,其曲线方程表示为 $y = f (x)$ ,土体黏聚力为 $c$ ,内摩擦角为 $φ$ ,重度为 $γ$ ,土体的重力为 $W$ ,第 $i$ 排筋条与水平方向的夹角为 $α_{i}$ ,水平方向布置则 $α_{i} = 0$ ；第 $i$ 排筋条与竖直方向的夹角为 $β_{i}$ ,每排筋条的单位方向向量记为 $\vec{T} = (\cos α_{i}, \cos β_{i})$ 。为合理建立力学模型,做如下说明：①假设滑动体整体为刚性体,在滑裂面上存在弹性压缩变形和剪切变形,储存的弹性压缩势能为 $V_{e}$ ,剪切势能为 $V_{τ}$ ；②在合外力 $\vec{R}$ 的作用下加筋边坡与滑体产生相同的虚位移 $\vec{d} =$ $(d_{1}, d_{2})$ ；③微分段的长度为 $d l$ ,弹性压缩变形可用刚度系数为 $K_{i}$ 的弹簧模拟,其大小与微段 $d l$ 呈正比,即 $K_{i} = m_{i} d l$ , $m_{i}$ 为地基系数,可通过试验测定或查阅相关资料获得。

图 1 H-V加筋均质边坡计算图示

Figure 1. Calculation of H-V ribs-reinforced homogeneous slope

下载: 全尺寸图片幻灯片

滑裂面方程为 $y = f (x)$ ,则求得曲线上任意切线上法线的斜率为 $k = - 1 / f^{'} (x)$ 。滑裂面外法线单位向量和滑面切向单位向量可表示为

$\vec{n} (x) = [n_{1} (x), n_{2} (x)] = {\frac{f^{'} (x)}{\sqrt{1 + {[f^{'} (x)]}^{2}}}, \frac{- 1}{\sqrt{1 + {[f^{'} (x)]}^{2}}}}$ ,

(1)

$\vec{t} (x) = [t_{1} (x), t_{2} (x)] = {\frac{1}{\sqrt{1 + {[f^{'} (x)]}^{2}}}, \frac{f^{'} (x)}{\sqrt{1 + {[f^{'} (x)]}^{2}}}}$ 。

(2)

假设在各种力的作用下,滑体发生虚位移,使系统势能最小,滑体此时的势能即是系统达到极限平衡状态时的最小势能,则任意微段储存的弹性势能为

$V_{e1} = \sum_{i = 1}^{N} \frac{1}{2} K_{i} {(\vec{d} \cdot \vec{n})}^{2} = \frac{1}{2} \int_{l} m {(\vec{d} \cdot \vec{n})}^{2} d l$ 。

(3)

2. H-V立体筋势能构建

H-V立体筋结构在传统水平筋的基础上设置竖筋,竖筋宽度可与水平筋相同,其截面形状可根据工程需要采用不同厚度变化的矩形、多边形等,形成一种空间形式的三维组合加筋结构,如图2所示。

图 2 H-V立体筋三维结构示意图

Figure 2. Three-dimensional structural diagram of H-V ribs

下载: 全尺寸图片幻灯片

H-V立体筋结构的基本原理：将一定量具有抗拉强度的H-V立体筋铺设在土体中,水平筋的拉力作用会逐渐延伸到路堤边坡自由面的筋条中,水平筋土之间摩擦力提供锚固力；竖筋的作用则表现为增加筋材拔出过程中竖筋的抗阻力,使土体结合更紧密,将整个H-V筋锚固在了土体中,提高边坡承载力。相对普通均质边坡,H-V筋的这种拉力状态也会使上部土体传递来的荷载较为集中,相对更均匀地分布到筋材下部的土体,使下部土体受力更合理。

2.1 筋条特征坐标的确定

构建H-V立体筋势能函数时需要确定水平筋与滑面交点坐标和竖筋有效数目,因此首先需要确定筋条的特征坐标。加筋边坡在发生失稳破坏时,根据坡体的滑动状态可分为滑动区和锚固区,如图3所示,设坡面倾角为 $θ$ ,筋材各层的垂直间距为 $h_{i}$ ,第一层水平筋距离坡脚的竖直距离为 $h_{1}$ 。则第 $i$ 层水平筋与滑动面的交点纵坐标 $y_{i}$ 可表示为

图 3 筋条的特征坐标计算简图

Figure 3. Simplified calculation of characteristic coordinates of geogrid

下载: 全尺寸图片幻灯片

$y_{i} = y_{1} + h_{i} (i - 1)$ 。

(4)

将式（4）与滑动面方程 $y = f (x)$ 联立即可确定各排筋条与滑面的交点坐标,若方程组有解,则筋条产生变形储存了势能,确认筋条储存势能且取 $x > 0$ 时,第 $i$ 层水平筋与滑动面的交点横坐标 $x_{i}$ 可表示为

$x_{i} = f^{'} (y_{i}) (x > 0)$ 。

(5)

土工格栅长度为 $l$ ,则第 $i$ 层水平筋右末端横坐标 ${x^{'}}_{i}$ 表示为

${x^{'}}_{i} = x_{1} + l + h_{i} (i - 1) / \tan θ$ 。

(6)

由式（5）,（6）可得第 $i$ 层土工格栅的锚固区长度为

$l_{i} = {x^{'}}_{i} - x_{i}$ 。

(7)

设竖筋间距为D,则锚固区内第 $i$ 层内竖筋的个数 $n_{i}$ 可表示为

$n_{i} = l_{i} / D + 1$ 。

(8)

则可计算出锚固区内竖筋的总个数 $n = \sum_{i = 1}^{m} n_{i}$ 。

2.2 水平筋弹性势能

在土体滑动的过程中,水平筋与土体之间摩擦力作功会产生摩擦势能,水平筋自身也会由于拉伸变形产生线弹性应变能。

对水平筋进行作用力分析时,由于筋土之间接触复杂,水平筋和填土分开考虑计算筋材的摩擦力参数较多,实现起来有很大的难度。加筋路堤稳定性计算时,土工合成材料应用手册将水平筋与路堤填土分开考虑,所以在计算稳定性系数时水平筋的作用力一般直接取其抗拉强度值。本文将水平筋与其附近土体考虑为一种复合体,所以筋条的抗拉力应该分散到其作用范围以内的土体中^[14]。假设筋条的抗拉力平均分配到每层筋条的土体中,若路堤中每层筋条水平间距为 $S_{x}$ ,竖向间距为 $S_{y}$ （图4）。

图 4 水平筋复合抗拉力计算示意图

Figure 4. Schematic diagram of calculation of horizontal ribs

下载: 全尺寸图片幻灯片

则第 $i$ 个土条中单位厚度复合土体中抗拉力为

$T_{st} = \frac{T_{h}}{S_{x} \cdot S_{y}} l_{i} \sin θ_{i}$ ,

(9)

式中, $T_{h}$ 为单根水平筋条抗拉力, $S_{x}$ 为筋条水平间距, $S_{y}$ 为筋条竖向间距, $l_{i}$ 为第 $i$ 条底弧长度, $θ_{i}$ 为第 $i$ 条底弧仰角。

如图5所示,对单位长度为 $l_{i}$ 的土体单元进行分析,水平筋的摩擦阻力 $T_{k i}$ 与水平筋的作用力 $T_{st}$ 平衡,所以满足 $T_{k i} = T_{st}$ ,由于土工格栅所受摩擦阻力与位移呈钝角关系,所以摩擦所储存的势能为力所做功的负值。

图 5 土体单元水平筋受力分析图

Figure 5. Analysis diagram of horizontal reinforcement of soil element

下载: 全尺寸图片幻灯片

$V_{e2} = - \sum_{i = 1}^{n} T_{k i} (\vec{d} \cdot \vec{T})$ ,

(10)

式中, $T_{k i}$ 为第 $i$ 排水平筋摩擦阻力, $\vec{d}$ 为滑体虚位移, $\vec{T}$ 为水平筋水平轴向单位向量。

根据以往文献资料结果显示,基于筋-土界面加筋摩擦机理的圆弧滑动分析方法可以计算加筋后的安全系数,在抗滑项中加入界面摩擦阻力,但结果可能偏小^[16]。本文将水平筋与其附近土体考虑为一种复合体,因为筋土间的界面被忽略,在计算水平筋所受摩擦阻力时缺少界面摩擦阻力,抗滑力中缺少界面摩擦项,且土工格栅摩擦阻力是在抗拉强度条件下计算的,因此可能导致安全系数计算结果偏大,具体影响还需要进一步研究。

分析水平筋的线应变势能时,滑体在向最小势能方向移动时会产生虚位移,因筋材均必须要伸入到较深坡体内部,因此可假设筋材与右端与滑裂面的交点g处是固定的^[11]。设布置了 $m$ 排筋材,筋条的刚度为 $K^{'}$ ,如图6所示,在合外力作用下滑体发生虚位移 $\vec{d}$ ,水平筋的坡面端头由 $e$ 移至 $e^{'}$ ,伸长量为 $Δ_{k^{'}}$ ,假设 $e$ 坐标为（ $x_{e}$ , $y_{e}$ ）, $g$ 的坐标为（ $x_{g}$ , $y_{g}$ ）。可得

图 6 立体筋材结构受力示意图

Figure 6. Schematic stress diagram of three-dimensional rib structure

下载: 全尺寸图片幻灯片

$\begin{matrix} Δ_{k^{'}} = | g e^{'} - g e | \\ = \sqrt{{(d_{1} + v_{k^{'} x})}^{2} + {(d_{2} + v_{k^{'} y})}^{2}} - \sqrt{v_{k^{'} x}^{2} + v_{k^{'} y}^{2}} \end{matrix}$ ,

(11)

式中, $v_{k^{'} x} = x_{g} - x_{e}, v_{k^{'} y} = y_{g} - y_{e}$ 。

水平筋的线应变能为

$\begin{matrix} V_{e3} = \sum_{i = 1}^{m} [\frac{1}{2} K^{'} Δ_{k^{'}}^{2}] \\ = \frac{K^{'}}{2} \sum_{i = 1}^{m} [{(d_{1} + v_{k^{'} x})}^{2} + {(d_{2} + v_{k^{'} y})}^{2} + (v_{k^{'} x}^{2} + v_{k^{'} y}^{2})] - \\ K^{'} \sum_{i = 1}^{m} [\sqrt{{(d_{1} + v_{k^{'} x})}^{2} + {(d_{2} + v_{k^{'} y})}^{2}} \cdot \sqrt{v_{k^{'} x}^{2} + v_{k^{'} y}^{2}}] \end{matrix}$ 。

(12)

2.3 竖筋弹性势能

当H-V立体筋在荷载作用下竖筋首先发挥作用,在其自身刚度下对土体产生侧阻力,阻挡土体向两侧运动,这相当于给土体增加多个侧向约束,竖筋对土体的势能主要是侧阻力对竖筋作用产生的变形势能。

对竖筋的作用力进行分析,竖筋的侧向抵抗作用在竖筋间将形成“土体加强区”,以有效改变加筋土的受力状况^[6]。如图7所示,在竖筋两侧分别作用有主动土压力 $σ_{a}$ 和被动土压力 $σ_{p}$ ,沿竖筋表面作用有剪切应力 $τ_{v}$ , $τ_{h}$ ,在竖筋顶部作用有竖向压力 $σ_{b}$ 。由于竖筋表面积与厚度及其顶部面积相对路堤尺寸可忽略,因此竖筋的主要作用力为竖向筋条两侧的主动土压力及被动土压力,在稳定性分析时,忽略竖筋表面的剪应力作用,竖筋的侧阻力可由被动土压力与主动土压力之差计算得到^[14]。

图 7 竖筋单元受力图

Figure 7. Force diagram of three-dimensional tendons

下载: 全尺寸图片幻灯片

则主动、被动土压力计算公式为

$T_{a} = (γ H_{i} K_{a} - 2 c \sqrt{K_{a}}) A$ ,

(13)

$T_{P} = (γ H_{i} K_{P} + 2 c \sqrt{K_{P}}) A$ 。

(14)

竖筋作用力为 $T_{sh}$ ,

$\begin{matrix} T_{s h} = T_{p} - T_{a} = [T_{p} - T_{a}] \\ = [γ H_{i} (K_{p} - K_{a}) + 2 c (\sqrt{K_{p}} + \sqrt{K_{a}})] A \end{matrix}$ ,

(15)

式中, $c$ 为土体黏聚力, $γ$ 为土体重度, $H_{i}$ 为竖筋上覆土高度, $φ$ 为内摩擦角, $K_{a}$ 为主动土压力系数, $K_{p}$ 为被动土压力系数,且 $K_{a} = \tan^{2} (45^{°} - φ / 2)$ , $K_{a} =$ $\tan^{2} (45^{°} + φ / 2)$ , $A$ 为竖筋侧面积, $A = B h$ ,其中B为竖筋宽度, $h$ 为竖筋高度。

由于竖筋的侧向变形对H-V筋中竖筋作用发挥起重要作用,所以尽管竖筋高度相对边坡尺寸较小,也需要进行势能分析。在进行H-V加筋边坡稳定性分析时,考虑到竖筋的排列会对微分圆弧作用面的受力产生影响,所以在计算抗滑力时只计算锚固区（即破裂面外部）的有效竖筋阻力,锚固区内竖筋总个数为 $n$ ,则竖筋在侧阻力作用下产生的变形势能为

$V_{e4} = \sum_{i = 1}^{n} T_{sh} Δ_{x} = \sum_{i = 1}^{n} T_{sh} d_{1}$ ,

(16)

式中, $T_{sh}$ 为竖筋作用力, $Δ_{x}$ 为竖筋水平方向的变形量,大小为虚位移 $\vec{d}$ 在水平方向的分量。

通过式（15）可看出竖筋作用力 $T_{sh}$ 的作用点位于竖筋高度的中部,但由于竖筋高度 $h$ 较小,在计算立体筋合力时,竖筋作用力的方向可近似认为作用在水平筋上。则立体筋的合力为竖筋作用力 $T_{sh}$ 和水平筋作用力 $T_{st}$ 的和,即

$T_{i} = T_{s t} + n T_{sh}$ 。

(17)

3. 滑面剪切势能的构建

边坡的失稳往往是沿着土体抵抗力最小的方向滑动,即为使坡体系统势能最小的虚位移方向。由文献[8]可知采用弹性势能以及合外力做功的方式表示边坡的剪切势能,可大大简化计算量,同时便于程序的实现。为便于势能公式的表达,令 $Δ = \sqrt{R_{1}^{2} + R_{2}^{2}}$ , $Δ^{'} = \sqrt{d_{1}^{2} + d_{2}^{2}}$ ,则合外力 $\vec{R}$ 和虚位移 $\vec{d}$ 的方向向量可表示为

$\begin{array}{l} {\vec{R}}_{r} = \frac{1}{Δ} (R_{1}, R_{2}), \\ {\vec{d}}_{d} = \frac{1}{Δ^{'}} (d_{1}, d_{2}) 。 \end{array}}$

(18)

从图1中任取土条作为研究对象,土条的受力分析如图8所示。

图 8 边坡土条受力图

Figure 8. Force diagram of slope strips

下载: 全尺寸图片幻灯片

作用在土条上的力有法向力 $N_{i}$ ,剪应力 $τ_{i}$ ,合外力为 $R_{i}$ ,H-V立体筋合力为 $T_{i}$ ,土条左右面上的条间力为 $E_{i}$ , $E_{i + 1}$ 。由于边坡失稳沿着岩土体抵抗力最小的方向,即虚位移 $\vec{d}$ ,因此各个土条左右面上的条间力作用方向与虚位移方向一致,在计算过程中,将第 $i$ 个土条条间力合力记为 $E_{i}$ , $E_{i + 1}$ ,将第 $i + 1$ 个土条条间力合力记为 ${E^{'}}_{i + 1}$ , $E_{i + 2}$ , $Δ E = E_{i + 1} - E_{i}$ ,则各个土条沿虚位移方向满足平衡方程：

$N_{i} \cos α + τ_{i} l_{i} \cos β + T_{i} \cos ω + Δ E_{i} = R_{i} \cos ϕ$ 。

(19)

式中 $α$ 为任意微面上法向力与虚位移方向的夹角,且 $\cos α = \vec{n} \cdot (- {\vec{d}}_{d})$ ； $l_{i}$ 为滑裂面微段长度； $β$ 为剪应力与虚位移的夹角,且 $\cos β = \vec{t} \cdot (- {\vec{d}}_{d})$ ； $ω$ 为水平筋与虚位移的夹角, $\cos ω = \vec{T} \cdot (- {\vec{d}}_{d})$ ； $ϕ$ 为合外力与虚位移的夹角, $\cos ϕ = {\vec{R}}_{r} \cdot {\vec{d}}_{d}$ 。

由式（19）可知滑面上剪应力的解析解为

$τ_{i} = \frac{\frac{R_{i}}{Δ} {\vec{R}}_{i} \cdot {\vec{d}}_{d} + T_{i} \vec{T} \cdot {\vec{d}}_{d} - N_{i} {\vec{n}}_{i} \cdot {\vec{d}}_{d} - Δ E_{i}}{- {\vec{d}}_{d} \cdot \vec{t} l_{i}}$ 。

(20)

为得到线弹性范围内剪切应变势能产生的函数,需要先构建剪应变和剪应力的力学求解模型,土体由于剪切变形发生的剪应变为δ_i,在剪切力 $τ_{i}$ 的作用下滑体沿着滑裂面切向会产生相对位移,即由滑裂面上的点 $f$ 移动至点 $f^{'}$ ,受剪切变形影响的土体深度为 $h_{i}$ ,剪切位移记为 $d^{'}$ ,则

$δ_{i} = \frac{d^{'}}{h_{i}}$ 。

(21)

滑面上的剪切势能为

$\begin{array}{l} V_{τ} = \frac{1}{2} \sum_{i = 1}^{n} τ_{i} δ_{i} h_{i} \\ = - \frac{1}{2} \sum_{i = 1}^{n} \frac{\frac{R_{i}}{△} {\vec{R}}_{i} \cdot {\vec{d}}_{d} + T_{i} \vec{T} \cdot {\vec{d}}_{d} - N_{i} {\vec{η}}_{i} \cdot {\vec{d}}_{d} - Δ E_{i}}{- {\vec{d}}_{d} \cdot \vec{t} l_{i}} \frac{\vec{d} \cdot \vec{t}}{h_{i}} h_{i} l_{i} \\ = - \frac{1}{2} \sum_{i = 1}^{n} (\frac{{\vec{R}}_{i}}{△} \vec{R} \cdot \vec{d} + T_{i} \vec{T} \cdot \vec{d} - N_{i} {\vec{η}}_{i} \cdot \vec{d} - Δ E_{i}) \\ = V_{e 1} - \frac{1}{2} \vec{R} \cdot \vec{d} - \frac{1}{2} \sum_{i = 1}^{n} T_{i} \vec{T} \cdot \vec{d} + \frac{1}{2} \sum_{i = 1}^{n} Δ E_{i} \end{array}$ ,

(22)

式中, $\sum_{i = 1}^{n} Δ E_{i}$ 的物理意义为各个条块条间力的矢量和,由于各个条块之间的作用力属于作用力与反作用力,因此 $\sum_{i = 1}^{n} Δ E_{i} = 0$ ,则上式可化简为

$V_{τ} = V_{e1} - \frac{1}{2} \vec{R} \cdot \vec{d} - \frac{1}{2} \sum_{i = 1}^{n} T_{i} \vec{T} \cdot \vec{d}$ 。

(23)

综合式（3）、（10）、（12）、（17）、（23）可得H-V立体筋边坡系统的总势能为

$\begin{matrix} V = V_{e1} + V_{e2} + V_{e3} + V_{e4} + V_{τ} - \vec{R} \cdot \vec{d} \\ = 2 V_{e1} + V_{e 2} + V_{e 3} + V_{e4} - \frac{3}{2} \vec{R} \cdot \vec{d} - \frac{1}{2} \sum_{i = 1}^{n} T_{i} \vec{T} \cdot \vec{d} \end{matrix}$ 。

(24)

4. 安全系数的计算

安全系数 $F_{s}$ 为沿着虚位移 $\vec{d}$ 方向投影得到的抗滑力 $F_{m}$ 与下滑力 $F_{n}$ 的比值^[8],即

$F_{s} = \frac{F_{m}}{F_{n}},$

(25)

式中,抗滑力主要由微分圆弧 $d l$ 的极限抗滑力和H-V立体筋复合抗拉力产生的抗滑力提供,下滑力由合外力（包括滑体自重）提供。本文的抗滑力 $F_{m}$ 及下滑力 $F_{n}$ 均与虚位移相关,因此需先求解虚位移,再进行安全系数求解。

4.1 虚位移 $\vec{d}$ 的求解

根据最小势能原理可知,H-V立体筋边坡系统势能最小时,总势能函数在虚位移的一阶偏导处需满足

$\begin{array}{l} \frac{\partial V}{\partial d_{1}} = 0, \\ \frac{\partial V}{\partial d_{2}} = 0 。 \end{array}}$

(26)

将式（24）代入式（26）,即可求出系统最小势能时坡体的虚位移 $\vec{d}$ 。其中令

$a_{1} = \sum_{i = 1}^{n} k_{i} η_{x}^{2} = m \int_{x 1}^{x 2} \frac{{[f^{'} (x)]}^{2}}{\sqrt{1 + {[f^{'} (x)]}^{2}}} d x$ ,

(27)

$a_{2} = a_{3} = \sum_{i = 1}^{n} k_{i} η_{x}^{2} = m \int_{x 1}^{x 2} \frac{- f^{'} (x)}{\sqrt{1 + {[f^{'} (x)]}^{2}}} d x$ ,

(28)

$a_{4} = \sum_{i = 1}^{n} k_{i} η_{y}^{2} = m \int_{x 1}^{x 2} \frac{1}{\sqrt{1 + {[f^{'} (x)]}^{2}}} d x$ ,

(29)

$A = K^{'} [d_{1} + v_{k^{'} x} - \frac{2 (v_{k^{'} x}^{2} + v_{k^{'} y}^{2}) (d_{1} + v_{k^{'} x})}{\sqrt{[{(d_{1} + v_{k^{'} x})}^{2} + {(d_{2} + v_{k^{'} y})}^{2}] \cdot (v_{k^{'} x}^{2} + v_{k^{'} y}^{2})}}],$

(30)

$B = K^{'} [d_{2} + v_{k^{'} y} - \frac{2 (v_{k^{'} x}^{2} + v_{k^{'} y}^{2}) (d_{2} + v_{k^{'} x})}{\sqrt{[{(d_{1} + v_{k^{'} x})}^{2} + {(d_{2} + v_{k^{'} y})}^{2}] \cdot (v_{k^{'} x}^{2} + v_{k^{'} y}^{2})}}],$

(31)

$C = (\frac{1}{2} \sum_{i = 1}^{n} T_{i} + \sum_{i = 1}^{n} T_{k i} - \sum_{i = 1}^{n} T_{sh}) \cos α_{i} + \frac{3}{2} R_{x}$ ,

(32)

$D = (\frac{1}{2} \sum_{i = 1}^{n} T_{i} + \sum_{i = 1}^{n} T_{k i}) \cos β_{i} + \frac{3}{2} R_{y}$ 。

(33)

则式（26）可表示为

$\begin{array}{l} 2 (a_{1} d_{1} + a_{2} d_{2}) + A = C, \\ 2 (a_{4} d_{2} + a_{3} d_{1}) + B = D 。 \end{array}}$

(34)

根据方程组计算得出 $d_{1}$ , $d_{2}$ 。

4.2 抗滑力 $F_{m}$ 的求解

抗滑力 $F_{m}$ 由滑裂面上各微段的法向力N_i、极限剪切力τ_il_i以及水平加筋等效应力或H-V立体加筋应力沿着虚位移方向提供的抗滑力投影之和组成。本文研究的土工格栅对抗滑力的贡献主要体现在增加了滑体的剪切力。整个滑裂面上任意微段的法向力 $N_{i}$ 提供的下滑力为

$\begin{matrix} F_{m1} = \sum_{i = 1}^{n} N_{i} \cos α = \sum_{i = 1}^{n} m (\vec{d} \cdot \vec{n}) l_{i} \cos α \\ = - \sum_{i = 1}^{n} m [\frac{d_{1} f^{'} (x) - d_{2}}{\nabla}] {\vec{d}}_{d} \cdot \vec{n} l_{i} \end{matrix}$ ,

(35)

式中, $T_{s h}$ 为任意微面上法向力与虚位移方向的夹角,且 $\cos α = \vec{n} \cdot (- {\vec{d}}_{d})$ 。

由莫尔-库仑破坏准则,微分圆弧 $d l$ 滑动面上极限抗剪力在虚位移方向上所能提供的抗滑力 $F_{m2}$ 为

$\begin{matrix} F_{m2} = \int_{x 1}^{x 2} d F_{m2} \\ = \int_{x 1}^{x 2} [c + m (\vec{d} \cdot \vec{n}) \tan φ] d l \cos β \end{matrix}$ ,

(36)

式中, $β$ 为剪应力与虚位移的夹角, $\cos β = \vec{t} \cdot (- {\vec{d}}_{d})$ 。

任意微段上的H-V立体筋复合抗拉力提供的抗滑力 $F_{m3}$ 为

$F_{m3} = \sum_{i = 1}^{n} T_{i} \cdot \cos ω = \sum_{i = 1}^{n} (T_{st} + T_{sh}) \cdot \cos ω$ ,

(37)

式中, $ω$ 为水平筋与虚位移的夹角, $\cos ω = \vec{T} \cdot (- {\vec{d}}_{d})$ 。

在虚位移方向上总的抗滑力 $F_{m}$ 为

$F_{m} = F_{m1} + F_{m2} + F_{m3}$ 。

(38)

4.3 下滑力 $F_{n}$ 的求解

作用在滑体上的所有力中,只有合外力 $\vec{R}$ 沿着虚位移 $c = 10$ 方向上提供下滑力 $F_{n}$ ,则总下滑力为

$F_{n} = \sum_{i = 1}^{n} {\overset{\leftarrow}{R}}_{i} \cos ϕ = | \vec{R} | \frac{R_{1} d_{1} + R_{2} d_{2}}{Δ Δ^{'}}$ 。

(39)

5. 算例

如图9所示的H-V立体筋均质边坡^[15],边坡坡度为1∶1,边坡高度6 m,土体黏聚力 $c = 10$ kPa,重度 $γ = 18$ kN/m³,内摩擦角 $φ = 20^{\circ}$ 。该边坡采用的H-V立体筋材料为单向拉伸高密度聚乙烯拉筋,筋材的的抗拉刚度为10000 kN/m,土体的地基系数为100 kN/m³,共铺设5层,每层筋材长为6 m,第一层距坡顶0.6 m,最后一层距坡底0.6 m,竖筋间距D=20 cm,筋材之间的间距为1.2 m,竖筋高度 $H$ =10 mm。

图 9 H-V立体筋边坡示意图

Figure 9. Schematic diagram of reinforced slope

下载: 全尺寸图片幻灯片

建立如图9所示直角坐标系,通过4.5H法搜索滑动面圆心线,通过改变边坡坡度,得到不同的滑动面的方程,利用MATLAB编程试算求解出安全系数的数值最小的滑裂面即为最危险滑动面。则试算得出的边坡最危险滑动面圆心 $o$ 的坐标为（-1.02,10.72）,半径R=10.77 m,滑动面方程AC为

${(x + 1.02)}^{2} + {(y - 10.72)}^{2} = {10.77}^{2}$

。

5.1 算例计算结果分析

本文采用改进最小势能法计算未加筋、加水平筋及加H-V立体筋边坡的安全系数,并与文献[14,15]的水平条分瑞典法的计算结果进行比较,计算结果如表1所示。

表 1 与文献计算结果对比

Table 1. Comparison between calculated results and those in published papers

项目	改进最小势能法	文献[14]的方法	文献[15]的水平条分瑞典法
未加筋	1.178	1.204	1.215
加水平筋	1.563	1.576	1.749
加立体筋	1.728	1.719	2.121

下载: 导出CSV

| 显示表格

由表1可见,对未加筋边坡,本文采用的改进最小势能法计算的安全系数与文献[14,15]的计算结果相差最大值为0.037,相差最小值0.026,相对误差小于3.05%,证明本文的计算方法是合理可行的。对于加水平筋,加立体筋边坡稳定性提高了14%,说明立体筋加筋效果更明显。对比文献[15]的水平条分瑞典法计算结果,本文的计算结果明显偏小,原因在于对加筋边坡稳定性分析时考虑了单个条块水平条间合力为0和整体力矩平衡的条件,考虑了竖向条间力作用,使计算结果偏大。

5.2 不同土体参数算例计算结果对比

目前,文献[15]采用的水平条分法还是一种相对新的计算方法,为了更好地验证本文计算方法的可靠性,针对算例在表2中设置不同的参数与极限平衡法（Bishop法和Janbu法）在无筋（不加筋）、加筋（加水平筋和H-V立体筋）状态下的计算结果进行对比。

表 2 算例参数

Table 2. Parameters of examples

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

下载: 导出CSV

| 显示表格

根据以往文献资料表明,当边坡为单一地层且不考虑加固措施时,通过最小势能法计算边坡稳定性系数过程中可将地基系数m约去^[8],即单一均质边坡稳定性不受地基系数的影响。因为极限平衡法分析边坡稳定性参数中没有地基系数m,为避免地基系数测试结果取值误差产生干扰,以及为后文验证H-V立体筋对提高边坡稳定性的作用,设置表3进行改进最小势能法和极限平衡法无筋边坡计算结果对比。

表 3 无筋边坡算例计算结果对比

Table 3. Comparison of calculated results of examples of unreinforced slopes

对比分组	表2编号	改进最小势能法	极限平衡法		备注
对比分组	表2编号	改进最小势能法	Bishop	Janbu	备注
A	1	1.178	1.117	1.098	$φ$ 渐增
	2	1.180	1.144	1.137
	3	1.182	1.179	1.160
	4	1.184	1.207	1.192
B	5	1.117	1.242	1.231	c渐增
	6	1.238	1.337	1.329
	7	1.344	1.415	1.411
	8	1.594	1.588	1.584

下载: 导出CSV

| 显示表格

从表3计算结果可知,其他参数保持不变,设置两组对比方式,A组内摩擦角 $φ$ 逐渐增大,B组黏聚力c逐渐增大,从计算结果可以看出无筋边坡利用改进最小势能法与极限平衡法计算结果相对误差保持在10%范围以内,证明本文分析计算边坡稳定性的方法是合理可行的,且改进最小势能法计算结果偏大于Janbu法,偏小于Bishop法。

由表4可看出：①改进最小势能法计算结果跟极限平衡法对比,施加水平筋相对误差小于4.7%,施加立体筋相对误差小于2.2%,证明改进最小势能法与极限平衡方法的计算结果具有很好的一致性。②加筋边坡的安全系数均随土体强度指标的增加呈现递增的趋势,其他参数保持不变,对照组C内摩擦角 $φ$ 逐渐增大,对照组D黏聚力c逐渐增大,分析数值变化规律可知时D组加筋边坡安全系数值变化更大,可知黏聚力c对边坡稳定性影响较大。③综合表3,4计算结果还可分析得出,对比无筋边坡,施加水平筋边坡的安全系数提高33.4%,对比施加水平筋,施加立体筋边坡的安全系数提高16.5%,说明施加立体筋在提高边坡稳定性方面的效果更显着。在不同的条件下,采用改进最小势能法的计算结果对比极限平衡法（Bishop法和Janbu法）还存在微小差异,主要因为本文的计算方法与极限平衡法采用的计算模型不一样,且对安全系数的定义和计算仍存在差异性。

表 4 加筋边坡算例计算结果对比

Table 4. Comparison of calculated results of examples of reinforced slope

对比分组	表2编号	改进最小势能法		极限平衡法				备注
		改进最小势能法		Bishop		Janbu
		水平筋	H-V筋	水平筋	H-V筋	水平筋	H-V筋
C	1	1.563	1.749	1.585	1.735	1.559	1.723	φ渐增
	2	1.565	1.761	1.613	1.762	1.558	1.758
	3	1.567	1.784	1.645	1.791	1.561	1.782
	4	1.570	1.792	1.678	1.830	1. 576	1.804
D	5	1.630	1.919	1.646	1.927	1.627	1.930	$c$ 渐增
	6	1.824	2.078	1.844	2.105	1.815	2.111
	7	1.879	2.125	1.902	2.173	1.880	2.123
	8	1.979	2.221	2.020	2.250	1.972	2.218

下载: 导出CSV

| 显示表格

6. 结论

本文在改进最小势能边坡稳定性分析基础上,提出H-V立体筋边坡稳定性分析新方法,并通过算例验证该方法的有效性和合理性,得出以下两点结论。

（1）本文基于改进的最小势能算法,考虑水平筋的摩擦势能、线弹性应变能以及竖筋变形势能的影响,进一步分析了加筋边坡的剪切势能,并将其表示为弹性势能的函数,简化了计算过程。通过计算边坡最危险滑裂面沿虚位移方向静力平衡方程,得到适用于施加水平筋和H-V立体筋的边坡稳定性分析方法。

（2）利用改进最小势能法求解文献算例,将其计算结果与极限平衡法（Bishop法和Janbu法）的计算的结果相比较,相对误差在10%以内,表明本文分析方法是可行且合理的。相比施加水平筋,施加立体筋时边坡安全系数提高了14%,说明H-V立体筋在提高边坡稳定性方面效果显着。

图 1 层状剪切模型箱

Figure 1. Layered shear model box

下载: 全尺寸图片幻灯片

图 2 近河谷地铁车站分布形式

Figure 2. Distribution of near-valley subway station

下载: 全尺寸图片幻灯片

图 3 模型土动力剪切特性曲线

Figure 3. Dynamic shear characteristics of model soil

下载: 全尺寸图片幻灯片

图 4 河谷地形制作过程

Figure 4. Production of valley topography

下载: 全尺寸图片幻灯片

图 5 材料性能试验

Figure 5. Material performance tests

下载: 全尺寸图片幻灯片

图 6 地铁站模型设计

Figure 6. Structural design of model station

下载: 全尺寸图片幻灯片

图 7 结构模型的施工过程

Figure 7. Construction process of model structures

下载: 全尺寸图片幻灯片

图 8 不同频率的正弦波记录及Chuetsu地震记录的时程和频谱曲线

Figure 8. Time histories of acceleration of Sin waves with different frequencies and time histories and frequency spectra of Chuetsu seismic record

下载: 全尺寸图片幻灯片

图 9 传感器布置示意图

Figure 9. Layout of sensors

下载: 全尺寸图片幻灯片

图 10 V0场地和V2场地下结构的峰值加速度放大因子分布

Figure 10. Distribution of peak acceleration amplification factor of structures for V0-site and V2-site

下载: 全尺寸图片幻灯片

图 11 V0场地和V2场地下结构右侧壁的峰值动土压力

Figure 11. Peak dynamic earth pressures of structural right sidewalls for V0-site and V2-site

下载: 全尺寸图片幻灯片

图 12 结构顶底板与相邻土的加速度时程记录

Figure 12. Time histories of acceleration of top base slab of structures and adjacent soils

下载: 全尺寸图片幻灯片

图 13 M-O理论的主动和被动土压力示意图

Figure 13. Schematic diagram of active and passive earth pressures of Mononobe-Okabe theory

下载: 全尺寸图片幻灯片

图 14 V0场地和V2场地下结构的动应变分布

Figure 14. Distribution of dynamic strain of structures for V0-site and V2-site

下载: 全尺寸图片幻灯片

图 15 V0场地和V2场地结构中心柱和侧墙的最大弯矩值

Figure 15. Maximum bending moments of center columns and sidewall for V0-site and V2-site

下载: 全尺寸图片幻灯片

图 16 不同深宽比河谷近邻结构的峰值加速度放大因子比值

Figure 16. Distribution of peak acceleration amplification factor ratio of the structure near the valley with different depth-to-width ratios

下载: 全尺寸图片幻灯片

图 17 不同深宽比河谷近邻结构右侧壁的峰值动土压力比

Figure 17. Ratios of peak dynamic earth pressure along right sidewall of station near valley with different depth-to-width ratios

下载: 全尺寸图片幻灯片

图 18 不同深宽比河谷近邻车站结构C2中心柱的剪力

Figure 18. Shear forces of C2 column of station near valley with different depth-to-width ratios

下载: 全尺寸图片幻灯片

表 1 振动台系统主要技术参数

Table 1 Main technical parameters of shaking table system

项目	技术参数
台面尺寸	4 m×6 m
额定承载能力	70 t
控制自由度	水平纵向、水平横向、水平转动
额定行程	水平横、纵向±500 mm
额定速度	水平横、纵向±1000 mm/s
额定加速度	水平横、纵向±1.5g
工作频率	0.1~50 Hz
额定倾覆力矩	400 t·m

下载: 导出CSV

表 2 相似比设计

Table 2 Similarity ratios

物理量	相似关系	相似比
应变 $\varepsilon$	${S_\varepsilon }$	1
应力 $\sigma$	${S_\sigma } = {S_E}{S_\varepsilon }$	1/70
加速度 $a$	${S_a} = {S_E}S_l^{ - 1}S_\rho ^{ - 1}$	1
时间 $t$	${S_t} = {S_l}S_\rho ^{0.5}S_E^{ - 0.5}$	0.17
速度 $v$	${S_v} = S_\rho ^{ - 0.5}S_E^{0.5}$	0.17
频率 $w$	${S_w} = S_{_l}^{ - 1}S_\rho ^{ - 0.5}S_E^{0.5}$	5.91
压力 $P$	${S_P} = {S_\rho }{S_a}{S_l}$	1/70
质量 $m$	${S_m} = {S_\rho }S_l^3$	1.17×10^-5

下载: 导出CSV

表 3 试验河谷地形对比

Table 3 Comparison of test valley topography

河谷编号	河谷尺寸(深×宽)	宽深比	对应圆弧直径/m
V1	0.22 m×1.78 m	1∶8	3.77
V2	0.44 m×1.78 m	1∶4	2.22
V3	0.59 m×1.78 m	1∶3	1.93

下载: 导出CSV

表 4 试验加载工况

Table 4 Test conditions

工况	地震动类型	加速度峰值/g
1	白噪声(WN1)	0.05
2	Sin-2Hz	0.2
3	Sin-4Hz
4	Sin-8Hz
5	Sin-10Hz
6	Chuetsu-Kashiwazaki-2007
7	白噪声(WN2)	0.05

下载: 导出CSV

参考文献(15)

[1]	ANDERSON J G, BODIN P, BRUNE J N, et al. Strong ground motion from the Michoacan, Mexico, earthquake[J]. Science, 1986, 233(4768): 1043-1049. doi: 10.1126/science.233.4768.1043
[2]	IIDA H, HIROTO T, YOSHIDA N, et al. Damage to Daikai subway station[J]. Soils and Foundations, 1996, 36: 283-300. doi: 10.3208/sandf.36.Special_283
[3]	李平, 刘红帅, 薄景山, 等. 汶川Ms8.0地震河谷地形对汉源县城高烈度异常的影响[J]. 地球物理学报, 2016, 59(1): 174-184. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201601015.htm LI Ping, LIU Hongshuai, BO Jingshan, et al. Effects of river valley topography on anomalously high intensity in the Hanyuan town during the Wenchuan Ms8.0 earthquake[J]. Chinese Journal of Geophysics, 2016, 59(1): 174-184. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201601015.htm
[4]	CAO H, LEE V W. Scattering and diffraction of plane P waves by circular cylindrical canyons with variable depth-to-width ratio[J]. Soil Dynamics and Earthquake Engineering, 1990, 9(3): 141-150. doi: 10.1016/S0267-7261(09)90013-6
[5]	高玉峰, 代登辉, 张宁. 河谷地形地震放大效应研究进展与展望[J]. 防灾减灾工程学报, 2021, 41(4): 734-752. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXK202104007.htm GAO Yufeng, DAI Denghui, ZHANG Ning. Progress and prospect of topographic amplification effects of seismic wave in canyon sites[J]. Journal of Disaster Prevention and Mitigation Engineering, 2021, 41(4): 734-752. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DZXK202104007.htm
[6]	庄海洋, 龙慧, 陈国兴, 等. 可液化地基中地铁车站周围场地地震反应分析[J]. 岩土工程学报, 2012, 34(1): 81-88. http://cge.nhri.cn/article/id/14492 ZHUANG Haiyang, LONG Hui, CHEN Guoxing, et al. Seismic responses of surrounding site of subway station in liquefiable foundation[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(1): 81-88. (in Chinese) http://cge.nhri.cn/article/id/14492
[7]	WU W F, GE S, YUAN Y, et al. Seismic response of a cross interchange metro station in soft soil: physical and numerical modeling[J]. Earthquake Engineering and Structural Dynamics, 2021, 50(9): 2294-2313. doi: 10.1002/eqe.3446
[8]	CHEN G X, CHEN S, ZUO X, et al. Shaking-table tests and numerical simulations on a subway structure in soft soil[J]. Soil Dynamics and Earthquake Engineering, 2015, 76: 13-28. doi: 10.1016/j.soildyn.2014.12.012
[9]	ZHANG S, YUAN Y, YANG Y, et al. Dynamic performance of soil-tunnel system under transverse sinusoidal excitations[J]. Applied Sciences, 2021, 11(11): 5079. doi: 10.3390/app11115079
[10]	WOOD D M, CREWE A, TAYLOR C. Shaking table testing of geotechnical models[J]. International Journal of Physical Modelling in Geotechnics, 2002, 2: 1-13. doi: 10.1680/ijpmg.2002.020401
[11]	张中杰, 禹海涛, 宋毅, 等. 圆形河谷地形效应振动台模型试验研究[J]. 岩土工程学报, 2024, 46(5): 1102-1111. doi: 10.11779/CJGE20221586 ZHANG Zhongjie, YU Haitao, SONG Yi, et al. Shaking table test on the topographic effect of circular arc valley[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(5): 1102-1111. (in Chinese) doi: 10.11779/CJGE20221586
[12]	建筑砂浆基本性能试验方法标准: JGJ/T 70—2009[S]. 北京: 中国建筑工业出版社, 2009. Standard for Test Method of Basic Properties of Construction Mortar: JGJ/T 70—2009[S]. Beijing: China Architecture & Building Press, 2009. (in Chinese)
[13]	WANG J N. Seismic design of tunnels: a simple state-of-the-art design approach[J]. New York: Parsons Brinckerhoff, 1993.
[14]	MONONOBE H. Considerations into earthquake vibrations and vibration theories[J]. Journal of the Japan Society of Civil Engineers, 1924, 10(5): 1063-1094.
[15]	牌立芳. 高烈度地震区隧道下穿滑坡体系灾变演化机制与控制技术研究[D]. 北京: 中国铁道科学研究院, 2022. PAI Lifang. Research on Catastrophic Evolution Mechanism and Control Technology of Tunnel Underpass Landslide System in High Intensity Earthquake Area[D]. Beijing: China Academy of Railway Sciences, 2022. (in Chinese)

施引文献(7)

期刊类型引用(4)

1.	邢亮. 基于判别分析法和机器学习对边坡稳定性的预测研究. 广东土木与建筑. 2024(12): 25-29+50 . 百度学术
2.	郝广杰. 考虑土层变形特征的边坡稳定性分析方法研究. 粉煤灰综合利用. 2023(02): 62-67 . 百度学术
3.	温树杰，赖光甜，孙自立. 改进的边坡最小势能法及潜在滑动方向研究. 铁道科学与工程学报. 2022(08): 2249-2258 . 百度学术
4.	何军. 深层岩土开采中边坡水土流失控制方法研究. 环境科学与管理. 2022(09): 148-153 . 百度学术

其他类型引用(3)

资源附件(1)

其他相关附件
- PDF格式
  02-202409-G23-0379⼀⽂审稿意⻅与作者答复点击下载(514KB)

图(18) / 表(4)

计量

文章访问数: 440
HTML全文浏览量: 62
PDF下载量: 118
被引次数: 7

0. 引言
1. 滑体储存的弹性势能
2. H-V立体筋势能构建
2.1 筋条特征坐标的确定
2.2 水平筋弹性势能
2.3 竖筋弹性势能
3. 滑面剪切势能的构建
4. 安全系数的计算
4.1 虚位移d→的求解
4.2 抗滑力Fm的求解
4.3 下滑力Fn的求解
5. 算例
5.1 算例计算结果分析
5.2 不同土体参数算例计算结果对比
6. 结论

0. 引言
1. 滑体储存的弹性势能
2. H-V立体筋势能构建
2.1 筋条特征坐标的确定
2.2 水平筋弹性势能
2.3 竖筋弹性势能
3. 滑面剪切势能的构建
4. 安全系数的计算
4.1 虚位移 $\vec{d}$ 的求解
4.2 抗滑力 $F_{m}$ 的求解
4.3 下滑力 $F_{n}$ 的求解
5. 算例
5.1 算例计算结果分析
5.2 不同土体参数算例计算结果对比
6. 结论

参考文献(15)

施引文献

资源附件(1)

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

近河谷地铁车站振动台模型试验研究 English Version

作者简介: 禹海涛(1983—)，男，工学博士，教授，博士生导师，主要研究方向为地下工程抗震防灾减灾。E-mail: yuhaitao@tongji.edu.cn

计量

出版历程

Shaking table tests on a near-valley subway station

0. 引言

1. 滑体储存的弹性势能

2. H-V立体筋势能构建

2.1 筋条特征坐标的确定

2.2 水平筋弹性势能

2.3 竖筋弹性势能

3. 滑面剪切势能的构建

4. 安全系数的计算

4.1 虚位移d→<math xmlns="http://www.w3.org/1998/Math/MathML" id="M163"><mover accent="true"><mi>d</mi><mo>→</mo></mover></math>的求解

4.2 抗滑力Fm<math xmlns="http://www.w3.org/1998/Math/MathML" id="M176"><mrow><msub><mi>F</mi><mtext>m</mtext></msub></mrow></math>的求解

4.3 下滑力Fn<math xmlns="http://www.w3.org/1998/Math/MathML" id="M193"><mrow><msub><mi>F</mi><mi>n</mi></msub></mrow></math>的求解

5. 算例

5.1 算例计算结果分析

5.2 不同土体参数算例计算结果对比

6. 结论

期刊类型引用(4)

其他类型引用(3)

其他相关附件

PDF格式

计量

出版历程

目录

0. 引言

1. 滑体储存的弹性势能

2. H-V立体筋势能构建

2.1 筋条特征坐标的确定

2.2 水平筋弹性势能

2.3 竖筋弹性势能

3. 滑面剪切势能的构建

4. 安全系数的计算

4.1 虚位移d→<math xmlns="http://www.w3.org/1998/Math/MathML" id="M163"><mover accent="true"><mi>d</mi><mo>→</mo></mover></math>的求解

4.2 抗滑力Fm<math xmlns="http://www.w3.org/1998/Math/MathML" id="M176"><mrow><msub><mi>F</mi><mtext>m</mtext></msub></mrow></math>的求解

4.3 下滑力Fn<math xmlns="http://www.w3.org/1998/Math/MathML" id="M193"><mrow><msub><mi>F</mi><mi>n</mi></msub></mrow></math>的求解

5. 算例

5.1 算例计算结果分析

5.2 不同土体参数算例计算结果对比

6. 结论

作者简介:
禹海涛(1983—)，男，工学博士，教授，博士生导师，主要研究方向为地下工程抗震防灾减灾。E-mail: yuhaitao@tongji.edu.cn

4.1 虚位移 $\vec{d}$ 的求解

4.2 抗滑力 $F_{m}$ 的求解

4.3 下滑力 $F_{n}$ 的求解

4.1 虚位移 $\vec{d}$ 的求解

4.2 抗滑力 $F_{m}$ 的求解

4.3 下滑力 $F_{n}$ 的求解

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200

编号	黏聚力 $c$ /kPa	重度 $γ$ /(kN·m^-3)	内摩擦角 $φ$ /(°)	地基系数 $m$ /(kN·m^-3)
1	10	18	20	100
2	10	18	22	100
3	10	18	24	100
4	10	18	26	100
5	5	20	30	200
6	10	20	30	200
7	15	20	30	200
8	20	20	30	200