竖向受荷横截面异形桩弹塑性变分解

李籼橙; 周航; 亓戈平

doi:10.11779/CJGE20211472

竖向受荷横截面异形桩弹塑性变分解 English Version

李籼橙^{1, 2,},
周航^{1, 2, ,},
亓戈平^{1, 2}

1.
重庆大学土木工程学院, 重庆 400045
2.
重庆大学山地城镇建设与新技术教育部重点实验室, 重庆 400045

基金项目:

国家自然科学基金面上项目 51978105

详细信息

作者简介:
李籼橙(1998—)，男，硕士研究生，主要从事桩-土相互作用方面的理论研究。E-mail: cqulxc1113@163.com

通讯作者:
周航, E-mail: zh4412517@163.com

中图分类号: TU473
计量
- 文章访问数: 199
- HTML全文浏览量: 58
- PDF下载量: 66
出版历程
- 收稿日期: 2021-12-10
- 网络出版日期: 2023-02-03
- 发布日期: 2021-12-10
- 刊出日期: 2022-12-31

Elasto-plastic variational solution for vertically loaded noncylindrical piles

LI Xiancheng^{1, 2,},
ZHOU Hang^{1, 2, ,},
QI Geping^{1, 2}

1.
School of Civil Engineering, Chongqing University, Chongqing 400045, China
2.
Key Laboratory of New Technology for Construction of Cities in Mountain Area, Chongqing University, Chongqing 400045, China

摘要

摘要: 现有计算理论无法严格考虑弹塑性土中竖向受荷横截面异形桩的异形效应。针对该问题，基于虚功原理推导出桩-弹塑性土模型的控制方程，采用保角变换技术克服了复杂边界条件控制方程求解问题，结合双曲形D-P本构模型的本构积分算法获得弹塑性土中竖向受荷横截面异形桩的半解析算法，建立了能够考虑横截面异形效应的竖向受荷横截面异形桩荷载传递理论模型。将理论模型的预测结果与有限元计算结果对比，验证了理论模型的可靠性和计算的高效性。最后，通过参数分析探讨了横截面异形效应对荷载沉降曲线的影响，结果表明：在正常工作荷载下，横截面异形效应对桩顶沉降影响不大；横截面异形效应主要对桩的极限承载力大小产生影响。
- 弹塑性土 /
- 横截面异形桩 /
- 虚功原理 /
- 保角变换 /
- D-P本构模型 /
- 荷载传递
Abstract: The existing analytical or semi-analytical methods cannot rigorously capture the geometrical effects of noncylindrical piles in elasto-plastic soils. To solve this issue, based on the principle of virtual work, the governing equations for the pile-soil system are derived. Solving the governing equation with a complex boundary is overcome by the conformal mapping technique. A semi-analytical algorithm for the vertically loaded noncylindrical piles in elasto-plastic soils is developed by combing the use of integral algorithm for the hyperbolic D-P constitutive model. A general theoretical model for the load transfer of the vertically loaded noncylindrical piles considering the geometrical effects is proposed. The reliability and efficiency of the proposed semi-analytical method is validated by comparing the predicted results with those of FEM. Finally, detailed parametric studies are conducted to investigate the geometrical effects on the influences of load-settlement curve. It is found that the non-circular cross-section has insignificant influences on the settlement of the pile head under working loads, which principally affects the magnitude of the ultimate bearing capacity of the pile.
- elasto-plastic soil /
- noncylindrical pile /
- principle of virtual work /
- conformal mapping /
- D-P constitutive model /
- load transfer

HTML全文

0. 引言

近年来，海洋开发和利用不断加深，特别是中国海上风电、跨海大桥、海上石油平台、海底隧道、海底管道等重大海洋工程投入使用。与此同时，中国海域及邻区所处地震地质环境复杂，地震活动性较强，海域地震频发^[1]，海洋工程的抗震设防能力及重要性凸显，海域场地工程特性得到充分关注^[2]。针对中国海域及邻近区域的场地地震反应研究工作仍不足^[1]，对海域土层条件，场地特性及分类等研究还处在起步阶段^[3]，所以研究海域场地效应和海域场地分类可以为海洋工程抗震设防，海洋结构安全性评价提供帮助。

Nakamura^[4]提出基于地脉动的水平和竖向谱比（horizontal-to-vertical spectral ratio, HVSR）作为场地特征评估方法，国内外涌现很多基于HVSR的场地地震效应和场地分类的研究。其中Boore等^[5]认为水平向地震动与竖向地震动之比（H/V）可以消除震源、断层类型、传播路径等因素，用以分析场地效应；同时Boore等^[5]认为H/V谱比法可同时应用于陆域和海域场地的研究。Field等^[6]利用HVSR法研究了1994年北岭地震场地非线性放大效应。Bonilla等^[7]利用南加州地区的井下阵列钻孔记录对比了HVSR法和传统的场地特征研究方法，发现两种方法得到的场地卓越周期基本一致。Wen等^[8]利用HVSR法分析了集集地震的场地非线性反应。任叶飞等^[9]将广义线性反演法和HVSR法进行对比，得出HVSR方法可以应用于评估场地卓越周期。荣棉水等^[10]利用美国Garner Valley Downhole Array（GVDA）竖向观测台阵的记录，对比了传递函数法和HVSR法的差异，并指出在可忽略场地竖向放大的情况下，可以利用HVSR研究场地特性。Ren等^[11]分析了2008年汶川地震的场地非线性效应，并提出了5个关键参数来表征场地非线性。Dhakal等^[12]利用K-NET六座海底台站记录到的315次地震事件分析了海底台站的场地效应，当PGA > 50gal时，海底部分台站所在场地出现非线性。Kubo等^[13]基于日本DONET1观测网记录到的2016年日本Off-Mie地震，利用HVSR法分析了台站所在场地的场地效应，认为部分台站场地的非线性反应抑制了地震动的放大效应。Ji等^[14]将HVSR计算出的非线性响应程度（DNL）和非线性响应百分比（PNL）联合PGA和V_S30（或V_S20）作为变量，采用k-means方法识别非线性场地响应。Yaghmaei-Sabegh等^[15]利用2017年伊朗Ezgeleh地震所记录到的地震数据，提出了基于广义回归神经网络（generalized regression neural networks，GRNN）的HVSR分类方法，该方法仅利用地震记录就可进行场地类别的划分，并将其与Vs30场地分类方法进行对比，展现了GRNN方法的可靠性。周旭彤等^[16]利用HVSR方法对日本DONET1台网20座海底台站场地效应进行分析，研究结果表明，部分台站具有相似的场地特征，海底台站的布设方式和海底地形会对谱比结果产生影响。

本文基于日本K-NET台网的海底台站收集的地震数据，利用HVSR联合动态聚类的方法对海底台站的场地特征进行识别，给出了台站所在场地的分类结果，为研究海域场地效应提供参考。

1. 海域强地震动观测数据与谱比方法

1.1 海域强地震动观测现状及数据处理

海域强地震动观测与陆域相比，存在设备安置难度大、回收效率低及观测稳定性较差等技术难题，发展远落后于现有陆域观测水平及海洋工程建设需求。美国和日本较早开展海域地震动观测，中国台湾地区也布设了海底地震观测设备^[17]，美国在20世纪70年代（1977年）针对海底地震记录观测需求，提出了SEMS计划（seafloor earthquake monitor system，以下简称SEMS）^[3]，日本于90年代（1996年）在东京圈近海开展ETMC计划（earthquake and tsunami monitoring cable，以下简称ETMC）^[18]。美国SEMS计划开发了4代海底地震仪，安装在南加州海域，水深区间为数10~300 m。日本ETMC系统布置在相模湾，观测系统由海底强震仪及海啸压力传感器组构，水深区间在900~2000 m不等。本文收集、整理了日本K-net台网6个海底台站2000—2018年的海域地震动数据。海底台站信息见表 1，分布见图 1。

表 1 海底台站信息

Table 1. Information of seabed stations

台站编号	东经/(°)	北纬/(°)	海水深度/m
KNG201	139.9183	34.5956	2197
KNG202	139.8393	34.7396	2339
KNG203	139.6435	34.7983	902
KNG204	139.5711	34.8931	933
KNG205	139.4213	34.9413	1486
KNG206	139.3778	34.0966	1130

下载: 导出CSV

| 显示表格

图 1 海底台站与水深分布图

Figure 1. Distribution map of submarine stations and water depth

下载: 全尺寸图片幻灯片

为研究6个台站所在海域的场地效应，考虑噪声对地震记录的影响，本文对上述数据进行筛选，筛选原则如下：①震级大于4级且PGA < 20 cm/s²；②震源位置均在海域。据此选取338组地震记录，每组记录均包含三分量，震级震中距分布见图 2。由图 2可知，震级主要分布在4.0~5.5级，震中距主要分布在0~300 km。

图 2 震级震中距分布图

Figure 2. Distribution map of magnitude and distance

下载: 全尺寸图片幻灯片

针对原始地震记录可能出现的基线漂移、仪器噪声及环境噪声等信号偏差，对筛选后的海域地震动记录采用基线校正和滤波的方法进行处理，其中采用Boore提出的分段式基线对地震动记录进行校正^[19]，构建4阶Butterworth带通滤波器，截止频率为0.1~20 Hz，处理前后的地震记录对比见图 3，由图可知，未处理地震动将出现基线漂移，中、长周期信号失真。处理后的地震记录更符合实际物理过程。

图 3 地震记录处理

Figure 3. Processing of seismic records

下载: 全尺寸图片幻灯片

1.2 海域强地震动记录HVSR计算

在利用强震记录研究场地效应时，需尽量选取S波时段^[20]。姚鑫鑫等^[21]关于S波的截取对HVSR谱比结果影响的研究表明，截取S波可以更好地体现高频段场地放大效应。因此本文对筛选后的海域地震动记录截取S波到时，将Arias强度能量^[22]达到90%所对应的持时作为S波的结束时间。

对截取后的S波数据进行快速傅里叶变换（FFT），选取0.2~20 Hz频带范围，并利用Konno平滑方法^[23]对傅里叶振幅谱进行平滑，平滑系数b取20。图 4（a）表明，利用该平滑方法既不低估场地放大效应的幅值，又不影响卓越频率的识别，可以较好地体现各频段的信息。将平滑后的地震记录S波的EW、NS向利用几何平均合成水平向，并除以相应的UD向平滑后的傅里叶振幅谱，得到的HVSR结果见图 4（b）。

图 4 HVSR计算过程

Figure 4. Calculation process of HVSR

下载: 全尺寸图片幻灯片

2. 动态聚类学习原理与评价

2.1 动态k-means聚类学习原理

Zhang等^[24]完善了动态 $k$ -means聚类方法，针对反应谱数据开展了聚类分析，本文通过添加谱比序列标签，采用动态聚类算法实现HVSR谱比分类。动态聚类学习原理如下：针对 $m$ 个数据序列 $\boldsymbol{Z} =$ $\left\{\boldsymbol{z}_1, \boldsymbol{z}_2, \cdots, \boldsymbol{z}_m\right\}^{\mathrm{T}} \in \mathbb{R}^{m \times n}$ ，各序列由 $n$ 个样本点组成行向量，即 ${{\bf{z}}_i} = \{ {z_{i1}}, {z_{i2}}, \cdot \cdot \cdot , {z_{in}}\}$ ，其中 $i = 1, 2, \cdot \cdot \cdot , m$ 。通过聚类、迭代方法，将数据序列实现分类，即将计算序列 ${\bf{Z}}$ 和各类数据中心 $\boldsymbol{C}=\left\{\boldsymbol{c}_1, \boldsymbol{c}_2, \cdots, \boldsymbol{c}_k\right\}^{\mathrm{T}} \in \mathbb{R}^{k \times n}$ 的欧式距离聚类：

$P(\boldsymbol{U}, \boldsymbol{W}, \boldsymbol{C})=\sum\limits_{p=1}^k \sum\limits_{i=1}^m \sum\limits_{j=1}^n \lambda_{i p} w_{p j}\left(z_{i j}-c_{p j}\right)^2 \text{，}$

(1)

式中： $\boldsymbol{c}_p=\left\{c_{p 1}, c_{p 2}, \cdots, c_{p n}\right\}$ 是 $p$ 个类中心，由n个样本点组成的行向量。 $\boldsymbol{U} \in {\mathbb{Z}^{z \times k}}$ 为二元矩阵，当 ${\lambda _{ip}} = 1$ 表示第 $i$ 个序列分配给 $p$ 类，并且每个序列仅可分配给同类，即 $\sum\nolimits_{p = 1}^k {{\lambda _{ip}} = 1}$ 。通过设置权重矩阵 $\boldsymbol{W}$ ，可更好识别局部数据跨度内时间序列的相关性，提高聚类性能，每个序列的权重之和为1，即 $\sum\nolimits_{j = 1}^n {{w_{pj}} = 1}$ 。为了平滑在相邻序列间的权重，防止陷入局部最小，设置约束正则项，定义如下：

$Q_w(\boldsymbol{W})=\frac{\alpha}{2} \sum\limits_{p=1}^k \sum\limits_{j=1}^{n-1}\left(w_{p j}-w_{p(j+1)}\right)^2 \text{，}$

(2)

式中， $\alpha$ 为平滑参数， $\alpha = \sum\nolimits_{i = 1}^m {\sum\nolimits_{j = 1}^n {({z_{ij}} - \sum\nolimits_{i = 1}^m {{z_{ij}}/} } }$ $m{)^2}$ 。

综上，动态聚类可以被定义为具有正则化目标函数的优化问题：

$\begin{aligned} & \{\hat{\boldsymbol{U}}, \hat{\boldsymbol{W}}, \hat{\boldsymbol{C}}\}=\underset{\{\boldsymbol{U}, \boldsymbol{W}, \boldsymbol{C}\}}{\arg \min } P(\boldsymbol{U}, \boldsymbol{W}, \boldsymbol{C})+Q_{\mathrm{w}}(\boldsymbol{W}) \\ & \mathrm{s}.\text{t}.\left\{\begin{array}{l}{\displaystyle \sum _{i=1}^{m}{\lambda }_{ip}\ge {m}_{0}}且{\displaystyle \sum _{p=1}^{k}{\lambda }_{ip}=1};{\lambda }_{ip}\in \{0, 1\}\\ {\displaystyle \sum _{j=1}^{n}{w}_{pj}=1;}\qquad \qquad\qquad 0\le {w}_{pj}\le 1\end{array}\right\} \text{，} \end{aligned}$

(3)

式中， ${m_0}$ 为每类所存在的最小序列数。 $\sum\nolimits_{i = 1}^m {{\lambda _{ip}} \geqslant {m_0}}$ 代表 $p$ 类中的每一类的序列数至少是m₀。式（3）的第一项使得所有类内的加权总距离最小，第二项使相邻类间隔的权重更平滑，即类间的区别更明显。通过设置计算序列到各中心的欧式距离，当此距离越近表示序列越相似。

2.2 动态聚类学习流程

该方法旨在最小化目标函数式（3），式（3）的变量由二元矩阵 $\boldsymbol{U}$ ，权重矩阵 $\boldsymbol{W}$ 和聚类中心矩阵 $\boldsymbol{C}$ 构成。采用优化方法，即将其中两个变量作为已知量，采用信息更新方法优化第三变量，当迭代次数达到40次，目标函数（3）收敛条件为10^-6时停止更新。整个优化过程可以分为如下3步。

（1）给定假设的权重矩阵 $\boldsymbol{W}$ 和聚类中心 $\boldsymbol{C}$ ，当式（1）满足式（4），序列被分配到假设聚类中心的类中，各类中心与序列之间的权重距离被最小化。

${\displaystyle \sum _{i=1}^{m}{\lambda }_{ip}\ge {m}_{0}}且{\displaystyle \sum _{p=1}^{k}{\lambda }_{ip}=1};{\lambda }_{ip}\in \{0, 1\} 。$

(4)

（2）当确定权重矩阵 $\boldsymbol{W}$ 和二元矩阵 $\boldsymbol{U}$ 后，通过令式（1）的导数为0求解信息更新后的新聚类中心，具体过程如下：

$\frac{\partial P(\boldsymbol{U}, \boldsymbol{W}, \boldsymbol{C})}{\partial c_{p j}}=2 \sum\limits_{i=1}^n \lambda_{i p} w_{p j}\left(z_{i j}-c_{p j}\right)=0 \text{，}$

(5)

变换得

${c_{pj}} = {{\sum\limits_{i = 1}^m {{\lambda _{ip}}{c_{ij}}} } \mathord{\left/ {\vphantom {{\sum\limits_{i = 1}^m {{\lambda _{ip}}{c_{ij}}} } {\sum\limits_{i = 1}^m {{\lambda _{ip}}} }}} \right. } {\sum\limits_{i = 1}^m {{\lambda _{ip}}} }} \text{，}$

(6)

通过式（6）可得到更新后的聚类中心。

（3）通过第一、二步确定二元矩阵 $\boldsymbol{U}$ 和聚类中心 $\boldsymbol{C}$ ，此时当式（3）满足

$\sum\limits_{j = 1}^n {{w_{pj}} = 1} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} (0 \leqslant {w_{pj}} \leqslant 1) 。$

(7)

式（3）转化为计算权重矩阵 $\boldsymbol{W}$ 的二次规划问题，具体求解方法可以参考文献[25]。

2.3 动态k-means聚类方法验证

为验证该方法在HVSR分类中的性能，选取5%阻尼比的反应谱来进行方法验证。依据《建筑抗震设计规范：GB50011—2010》^[26]，在场地分类为III类条件下，生成两组设计反应谱，两类反应谱见图 5，图 5中第一类反应谱的特征周期为0.9 s，地震影响系数最大值为0.9；第二类反应谱特征周期为0.3 s，地震影响系数最大值为0.79。

图 5 设计反应谱

Figure 5. Design response spectra

下载: 全尺寸图片幻灯片

依据上述不同特性反应谱，从PEER强震动数据库中^[27]，分别针对各分类选取50条地震记录。图 6给出了数据的选取情况。

图 6 不同记录反应谱特性

Figure 6. Two clusters of response spectra

下载: 全尺寸图片幻灯片

将两组数据利用动态聚类学习方法进行分类测试，结果见图 7。由图 7可知，采用动态聚类学习方法划分的中心和反应谱均值无偏差，展现出良好的聚类性能及鲁棒性。

图 7 反应谱均值与中心匹配

Figure 7. Mean values of response spectra and cluster centroid

下载: 全尺寸图片幻灯片

3. 海域地震动谱比法分类

3.1 海域地震动谱比法分类依据

将HVSR结果利用上述动态聚类方法进行聚类。利用误差平方和准则评价聚类性能，表示聚类的畸变程度。给定数据序列 $\boldsymbol{Z}$ ，假设 $\boldsymbol{Z}$ 中的序列可分为 $k$ 类，即 ${\boldsymbol{Z}_1}, {\boldsymbol{Z}_2}, \cdot \cdot \cdot , {\boldsymbol{Z}_k}$ ；各类序列数分别为 ${m_1}, {m_2}, \cdot \cdot \cdot , {m_k}$ ；各类中心分别为 ${c_1}, {c_2}, \cdot \cdot \cdot , {c_k}$ ，则误差平方和准则函数为

$D = \sum\nolimits_{i = 1}^k {\sum\limits_{z \in {Z_i}} {{{(d(z, {c_i}))}^2}} } \text{，}$

(8)

式中， $D$ 为数据序列分为 $k$ 类时的每类数据序列与其中心序列的距离和。 $D$ 越小表示每类数据序列到对应中心的距离越小，越紧凑，分类性能越好。分别令 $k = 1, 2, \cdot \cdot \cdot , 8$ 进行聚类计算并计算 $D$ ，计算结果见图 8。从中可以看出，随着分类数增大，误差平方和趋于稳定，结合肘部法则^[28]确定最优分类数为 $k = 3$ 。此外为保证每类结果具有统计意义，设置每类最少存在6组HVSR数据序列，即 ${m_0} = 6$ 。

图 8 误差平方和与分类数关系

Figure 8. Identification of optimal number of cluster at the elbow point

下载: 全尺寸图片幻灯片

3.2 海域场地特性分析

HVSR动态聚类结果见图 9，对应台站统计如图 10所示。由图 9，10可知：海域场地HVSR曲线无明显峰值或多峰值现象，结合肘部算法，相模湾海区最优化场地可分为三类，结合相模湾海区地貌^[29-30]，综合分为海底峡谷区（KNG203，KNG204）及海洋盆地区（KNG201，KNG202，KNG205，KNG206），两类地貌分类的HVSR谱比呈显著差异，海底峡谷区高频衰减慢且优势频段主要集中在中高频，台站KNG204代表的第一类主要集中在12 Hz，体现了海底峡谷中少量沉积的薄覆盖层环境；台站KNG203代表的第三类优势频带跨度较大（5~10 Hz），台站KNG201，KNG202，KNG205，KNG206则表现为中长周期发育，高频衰减快，呈现出典型的深厚覆盖层的海洋盆地特性。

图 9 HVSR聚类结果

Figure 9. Clustering results of HVSR

下载: 全尺寸图片幻灯片

图 10 各海域台站HVSR与类别对应

Figure 10. HVSR of each station corresponding to category

下载: 全尺寸图片幻灯片

其次，局部场地地形效应也显著影响海域地震动传播，以台站作为中心点，以类似地貌作为依据，通过采样点与标准点之间的距离（500 m）与高程计算局部场地地形坡度：

$i = h/l \times 100\% \text{，}$

(9)

式中，h为采样点与标准点的平均高程差，l为采样点与标准点距离。

各台站所在局部场地地形坡度见表 2。由表 2可知：海底峡谷区台站KNG203较KNG204场地局部地形平缓度高，存在局部沉积形成的浅覆盖层；海洋盆地区台站KNG201，KNG205，KNG206所在场地局部地形坡度基本一致，KNG202所在场地处于盆地边缘，存在基岩面拉伸导致的局部场地不均匀分布特性。

表 2 局部场地地形坡度

Table 2. Topographic slopes of local site

类别	台站	地貌	坡度	地形	场地特征
一	KNG204	海洋峡谷	4.62	陡峭	薄覆盖层
二	KNG201	海洋盆地	1.39	平缓	厚覆盖层
	KNG205		1.85
	KNG206		1.02
	KNG202		3.81	陡峭	薄覆盖层
三	KNG203	海洋峡谷	0.38	平缓	浅覆盖层

下载: 导出CSV

| 显示表格

4. 结论

本文采用动态聚类学习方法结合海域场地HVSR谱比，分析了动态聚类方法在地震动分类的适用性；给出了相模湾海区场地工程特性，结合海洋工程地质给出了海域地形坡度与场地分类关系。希望通过对有限的海域场地地震记录研究，为近海以及海洋工程结构和设施的设计提供合理、科学的工程场地参考。得到以下2点结论。

（1）动态聚类学习方法针对时序数据的聚类性能及鲁棒性强，可满足局部及全局分辨率，应用于场地工程特性分析中，可结合HVSR谱比，最优化、自动化的给出场地分类特性；

（2）HVSR动态聚类的海洋场地特性表明：台站KNG203和KNG204位于海底峡谷区，台站KNG201、KNG202、KNG205、KNG206位于海洋盆地区，两类地貌分类对海域地震动传播影响较大，海底峡谷区高频衰减慢且优势频段主要集中在中高频，通过局部场地地形坡度分析给出了各台站所在场地的工程特性。

图 1 子午面上原D-P与双曲形D-P屈服面

Figure 1. Original D-P and hyperbolic D-P yield surfaces in meridional plane

下载: 全尺寸图片幻灯片

图 2 理想弹塑性土中竖向受荷横截面异形桩力学模型

Figure 2. Mechanical model for axially loaded noncylindrical piles embedded in elasto-plastic soils

下载: 全尺寸图片幻灯片

图 3 水平面内区域与边界的定义

Figure 3. Definition of region and boundary in xy plane

下载: 全尺寸图片幻灯片

图 4 相平面的定义

Figure 4. Definition of phase plane

下载: 全尺寸图片幻灯片

图 5 竖向受荷桩-土体系弹塑性分析求解方案

Figure 5. Solution scheme for elasto-plastic analysis of axially loaded pile-soil system

下载: 全尺寸图片幻灯片

图 6 对比本文解与有限元法中矩形桩和土的响应

Figure 6. Comparison of responses of rectangular piles and soils from present analysis and FE method

下载: 全尺寸图片幻灯片

图 7 矩形桩与等横截面面积圆桩的荷载沉降曲线

Figure 7. Load-settlement curves for rectangular and circular piles with same cross-sectional area

下载: 全尺寸图片幻灯片

图 8 对比本文解与有限元法中XCC桩和土的响应

Figure 8. Comparison of responses of XCC piles and soils from present analysis and FE method

下载: 全尺寸图片幻灯片

图 9 XCC桩与等横截面面积圆桩的荷载沉降曲线

Figure 9. Load-settlement curves for XCC and circular piles with same cross-sectional area

下载: 全尺寸图片幻灯片

图 10 不同参数对归一化的荷载-沉降曲线的影响

Figure 10. Influence of different parameters on normalized load-settlement curves

下载: 全尺寸图片幻灯片

图 11 不同参数对横截面异形桩承载力增加量的影响

Figure 11. Influences of different parameters on increment of bearing capacity of noncylindrical piles

下载: 全尺寸图片幻灯片

参考文献(26)

[1]	LI X C, ZHOU H, LIU H L, et al. Three‐dimensional analytical continuum model for axially loaded noncircular piles in multilayered elastic soil[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2021, 45(18): 2654-2681. doi: 10.1002/nag.3281
[2]	SEED H B, REESE L C. The action of soft clay along friction piles[J]. Transactions of the American Society of Civil Engineers, 1957, 122(1): 731-754. doi: 10.1061/TACEAT.0007501
[3]	MINDLIN R D. Force at a point in the interior of a semi-infinite solid[J]. Physics, 1936, 7(5): 195-202. doi: 10.1063/1.1745385
[4]	POULOS H G, DAVIS E H. The settlement behaviour of single axially loaded incompressible piles and piers[J]. Géotechnique, 1968, 18(3): 351-371. doi: 10.1680/geot.1968.18.3.351
[5]	HIRAI H. Settlement analysis of rectangular piles in nonhomogeneous soil using a Winkler model approach[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2014, 38(12): 1298-1320. doi: 10.1002/nag.2270
[6]	雷国辉, 洪鑫, 施建勇. 壁板桩承载特性的近似三维分析[J]. 岩土力学, 2004, 25(4): 590-594, 600. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200404020.htm LEI Guohui, HONG Xin, SHI Jianyong. Approximate three-dimensional analysis of load-carrying behaviour of barrettes[J]. Rock and Soil Mechanics, 2004, 25(4): 590-594, 600. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200404020.htm
[7]	VALLABHAN C V G, MUSTAFA G. A new model for the analysis of settlement of drilled piers[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1996, 20(2): 143-152. doi: 10.1002/(SICI)1096-9853(199602)20:2<143::AID-NAG812>3.0.CO;2-U
[8]	SHEN W Y, CHOW Y K, YONG K Y. A variational approach for vertical deformation analysis of pile group[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1997, 21(11): 741-752. doi: 10.1002/(SICI)1096-9853(199711)21:11<741::AID-NAG898>3.0.CO;2-D
[9]	SHEN W Y, CHOW Y K, YONG K Y. Variational solution for vertically loaded pile groups in an elastic half-space[J]. Géotechnique, 1999, 49(2): 199-213. doi: 10.1680/geot.1999.49.2.199
[10]	LEE K M, XIAO Z R. A new analytical model for settlement analysis of a single pile in multi-layered soil[J]. Soils and Foundations, 1999, 39(5): 131-143. doi: 10.3208/sandf.39.5_131
[11]	SEO H, PREZZI M. Analytical solutions for a vertically loaded pile in multilayered soil[J]. Geomechanics and Geoengineering, 2007, 2(1): 51-60. doi: 10.1080/17486020601099380
[12]	BASU D, PREZZI M, SALGADO R, et al. Settlement analysis of piles with rectangular cross sections in multi-layered soils[J]. Computers and Geotechnics, 2008, 35(4): 563-575. doi: 10.1016/j.compgeo.2007.09.001
[13]	TEHRANI F S, SALGADO R, PREZZI M. Analysis of axial loading of pile groups in multilayered elastic soil[J]. International Journal of Geomechanics, 2016, 16(2): 04015063. doi: 10.1061/(ASCE)GM.1943-5622.0000540
[14]	雷国辉, 詹金林, 洪鑫. 壁板桩群桩竖直荷载沉降关系的变分法分析[J]. 岩土力学, 2007, 28(10): 2071-2076. doi: 10.3969/j.issn.1000-7598.2007.10.012 LEI Guohui, ZHAN Jinlin, HONG Xin. Variational analysis of load-settlement relationship of vertically loaded barrette groups[J]. Rock and Soil Mechanics, 2007, 28(10): 2071-2076. (in Chinese) doi: 10.3969/j.issn.1000-7598.2007.10.012
[15]	周航, 李籼橙, 刘汉龙, 等. 均质土中竖向受荷X形混凝土桩的三维弹性变分解[J]. 岩土力学, 2021, 42(4): 1012-1024. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202104014.htm ZHOU Hang, LI Xiancheng, LIU Hanlong, et al. Three-dimensional variational elastic solution for axially loaded X-section cast-in-place concrete pile in homogeneous soil[J]. Rock and Soil Mechanics, 2021, 42(4): 1012-1024. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202104014.htm
[16]	JEONG S, LEE J, LEE C J. Slip effect at the pile-soil interface on dragload[J]. Computers and Geotechnics, 2004, 31(2): 115-126. doi: 10.1016/j.compgeo.2004.01.009
[17]	吕亚茹, 丁选明, 刘汉龙, 等. 刚性荷载下现浇X形桩复合地基桩侧摩阻力数值分析[J]. 岩土工程学报, 2012, 34(11): 2134-2140. http://cge.nhri.cn/cn/article/id/14886 LÜ Yaru, DING Xuanming, LIU Hanlong, et al. Numerical analysis of side resistance of composite foundation with X-section cast-in-place concrete piles under rigid load[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(11): 2134-2140. (in Chinese) http://cge.nhri.cn/cn/article/id/14886
[18]	ZHANG Q Q, FENG R F, YU Y L, et al. Simplified approach for prediction of nonlinear response of bored pile embedded in sand[J]. Soils and Foundations, 2019, 59(5): 1562-1578. doi: 10.1016/j.sandf.2019.07.011
[19]	HAN F, SALGADO R, PREZZI M. Nonlinear analyses of laterally loaded piles-A semi-analytical approach[J]. Computers and Geotechnics, 2015, 70: 116-129. doi: 10.1016/j.compgeo.2015.07.009
[20]	GUPTA B K, BASU D. Nonlinear solutions for laterally loaded piles[J]. Canadian Geotechnical Journal, 2020, 57(10): 1566-1580. doi: 10.1139/cgj-2019-0341
[21]	周永强, 盛谦, 刘芳欣, 等. 一种修正的Drucker-Prager屈服准则[J]. 岩土力学, 2016, 37(6): 1657-1664. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201606016.htm ZHOU Yongqiang, SHENG Qian, LIU Fangxin, et al. A study of modified Drucker-Prager yield criterion[J]. Rock and Soil Mechanics, 2016, 37(6): 1657-1664. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201606016.htm
[22]	ABBO A J, SLOAN S W. A smooth hyperbolic approximation to the Mohr-Coulomb yield criterion[J]. Computers & Structures, 1995, 54(3): 427-441.
[23]	CHEN H H, LI L, LI J P, et al. A rigorous elastoplastic load-transfer model for axially loaded pile installed in saturated modified Cam-clay soils[J]. Acta Geotechnica, 2022, 17(2): 635-651. doi: 10.1007/s11440-021-01248-z
[24]	ZHOU H, LIU H L, YUAN J R. A novel analytical approach for predicting the noncylindrical pile penetration-induced soil displacement in undrained soil by combining use of cavity expansion and strain path methods[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2018, 42(11): 1270-1305.
[25]	詹云刚, 袁凡凡, 栾茂田. 纯摩擦型岩土介质本构积分算法及其在ABAQUS中开发应用[J]. 岩土力学, 2007, 28(12): 2619-2623, 2628. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200712029.htm ZHAN Yungang, YUAN Fanfan, LUAN Maotian. Integration algorithm of constitutive equation for cohesionless-frictional geomaterial and its application to ABAQUS[J]. Rock and Soil Mechanics, 2007, 28(12): 2619-2623, 2628. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200712029.htm
[26]	RANDOLPH M F. The response of flexible piles to lateral loading[J]. Géotechnique, 1981, 31(2): 247-259. http://www.onacademic.com/detail/journal_1000034855799910_85ab.html

施引文献(15)

期刊类型引用(13)

1.	侯鹏，刘尧，黄毅. 地铁区间隧道防沉“微扰动”注浆加固技术研究. 城市建设理论研究(电子版). 2025(02): 173-175 . 百度学术
2.	梁发云，魏圣明，陈可，袁强，闫静雅，顾晓强. 基于UH模型的盾构隧道非线性地基反力确定方法. 岩土工程学报. 2025(02): 284-295 . 本站查看
3.	张玄，肖颖，梁寿，王禹玺，王睿，胡黎明. 地热源作用下盾构地铁隧道管片应力响应规律. 公路. 2025(02): 443-449 . 百度学术
4.	梁荣柱，张志伟，李忠超，符宇坤，曹世安，康成，肖铭钊，吴文兵. 考虑环间接头影响的盾构隧道纵向变形简化解. 岩土工程学报. 2024(02): 254-263 . 本站查看
5.	王朋路，王微笑，冯畅. 城市砂土隧道注浆效果评价研究. 武汉理工大学学报(交通科学与工程版). 2024(01): 180-184 . 百度学术
6.	申子龙. 地铁盾构隧道掘进过程注浆技术研究. 工程机械与维修. 2024(02): 162-164 . 百度学术
7.	王观璐，徐德鹏，张辉. 市政配套工程施工对既有地铁车站的影响分析. 现代城市轨道交通. 2024(04): 97-105 . 百度学术
8.	黄大维，罗仲睿，罗文俊，朱碧堂，刘家璇，赵梽錡. 地层注浆附加土压力形成及影响因素单元体试验研究. 岩石力学与工程学报. 2024(S1): 3520-3529 . 百度学术
9.	袁钊. 城市复杂环境条件下盾构隧道变形治理监控关键措施. 建设监理. 2024(07): 131-134 . 百度学术
10.	刘禹，刘性锋，李世堂，周静. 考虑多因素的深埋管线施工对临近地铁影响分析——以龙泉驿区玉竹路为例. 科技和产业. 2024(18): 327-332 . 百度学术
11.	王虎，叶胜，罗天霖，宋敦安，李栋. 深厚软土区地表附加荷载对顶管隧道变形的影响分析. 现代隧道技术. 2024(S1): 1011-1017 . 百度学术
12.	袁强，赵旭，陈诗艾，李建平，袁炳祥. 盾构隧道纵缝接头抗弯性能评价指标及极限承载分析. 现代隧道技术. 2024(S1): 432-444 . 百度学术
13.	涂玖林，王春森，何永辉. 市政道路斜交上跨既有地铁隧道影响研究. 交通节能与环保. 2023(05): 275-279+284 . 百度学术