深厚软土区桩柱式桥墩临界荷载简化计算方法

Simplified Calculation Method for Critical Load of Pile Piers in Deep Soft Soil Regions

摘要

现有桩柱式桥墩临界荷载简化计算方法在大多数情况下计算误差较大，适用性有限.在利用软土区现场桩基屈曲实验结果标定的基础上，采用基于特征值屈曲分析的有限元方法对有典型桩柱式桥墩的临界荷载进行了计算并与现有简化计算方法结果进行了对比，指出了现有简化计算方法的适用范围，发现现有简化方法不适用于软土区桩柱式桥墩计算.在大量桩柱式桥墩有限元特征值屈曲分析数据的基础上，采用随机森林算法对计算结果进行分析，得到可以计算桩柱式桥墩计算长度的随机森林模型以及各影响参数重要性指数分析结果.最后以现常用简化公式为基本形式，纳入随机森林算法给出的重要参数，通过回归分析，提出了深厚软土区桩柱式桥墩临界荷载的简化计算方法.

Abstract

Existing simplified calculation methods for critical load of pile piers in bridge structures exhibit significant errors in most cases， limiting their applicability. Based on calibration using field load tests on pile foundations in soft soil regions， a finite element method based on eigenvalue buckling analysis was employed to calculate the critical load of typical pile piers and compared with the results obtained from existing simplified calculation methods. The study identified the applicability range of the existing simplified calculation methods and found that they are not suitable for calculating pile piers in soft soil regions. Using a large dataset of finite element eigenvalue buckling analysis results for pile piers， the random forest algorithm was applied to analyze the calculation results， yielding a random forest model capable of calculating the length of pile piers and determining the importance index of various influencing parameters. Finally， based on the form of commonly used simplified calculation formulas and incorporating the important parameters identified by the random forest algorithm， a simplified calculation method for determining the critical load of pile piers in deeply soft soil regions was proposed through regression analysis.

Graphical abstract

关键词

桩柱式桥墩 / 临界荷载 / 有限元 / 随机森林

桥梁桩基屈曲问题是桥梁设计中不可避免的关键问题之一，对于深厚软土地区尤其如此.20世纪20年代以来，学者们对桩基稳定性问题进行了大量的理论、试验和数值分析研究，得到了丰富的成果.赵明华等^［1-3］基于能量法推导出桩柱式桥墩在不同边界条件下的计算长度公式，开展了室内模型箱试验验证其准确性，并提出了有限元分析桩基屈曲问题的方法.彭锡鼎^［4］基于伽辽金法考虑了桩侧土抗力，推导出两端固支桩基的计算长度公式.邹新军等^［5］对桩基屈曲问题引入新型无单元Galerkin数值分析法，建立了无单元网格计算程序.雷勇等^［6］依据能量法建立了穿越单层溶洞基桩总势能方程，并引入尖点突变理论，导出了分岔集方程，提出了穿越单层溶洞基桩屈曲临界荷载计算方法.但是，要将这些科研成果应用于工程，需要开发复杂的程序或软件，因此设计中往往倾向于简化方法.对于桩柱式桥墩，基于欧拉公式导出的简化计算公式^［7］简单明了，是目前设计中最常用的计算方法［因其被写入《建筑桩基技术规范》（JGJ 94—2008）^［8］，下文简称该方法为规范法］，但是分析表明，它在理论上存在缺陷，计算结果并不准确.

经过多年发展，基于特征值屈曲分析的有限元方法已经在结构稳定性分析中得到广泛的应用^［9］，其结果也被证明与实际相符.本文拟在深厚软土区桩柱式桥墩临界荷载分析中引入特征值屈曲分析的有限元方法对桩柱式桥墩进行分析，在对比说明规范法的适用范围后，基于机器学习方法整理出简化计算公式供设计同行参考.

1 有限元法的特征值屈曲分析

1.1 计算原理

特征值屈曲分析的有限元方法是将桩身划分为若干单元，考虑每个单元所受的节点荷载、桩侧土抗力、自重、摩阻力，建立相应于该单元的桩土体系总势能方程，利用势能驻值原理，获得各单元的刚度矩阵［ k ］^e，将单元刚度矩阵组成与桩顶荷载P有关的总刚度方程后得到：

\{\delta \}={[\mathit{K}]}^{-\mathrm{1}}\{P\}（1）

式中：［ K ］为总刚度矩阵，由单元刚度矩阵［ k ］^e组装而成；｛δ｝为结构在荷载作用下的位移；｛P｝为结构上的荷载.

式（1）中，总刚度矩阵［ K ］的逆矩阵等于其伴随矩阵除以系数的行列式，当总刚度系数行列式等于零时，即：

\left|\mathit{K}\right|=\mathrm{0}（2）

此时，位移｛δ｝将趋于无穷大，说明结构丧失了稳定性，将发生失稳破坏，故式（2）通常称为桩身屈曲稳定的有限元特征方程.求解该n阶方程可得到n个特征值，对应n阶的屈曲模态，在实际工程中，当结构的荷载达到特征值中的最小值时即发生破坏，故其中最小特征值即桩身的临界荷载P_cr.

1.2 现场试验简介

采用作者早年完成的位于洞庭湖深厚软土区的茅草街大桥超长嵌岩灌注桩屈曲试验对有限元模型进行标定.现场地质条件如表1所示^［10］，桩周土体水平抗力系数m≈5 500 kN·m^-4，桩径1 m，计算宽度b₁=1.8 m，混凝土等级C30，弹性模量E=31 GPa，计算可得桩土变形系数α=0.381 9 m^-1.限于篇幅，试验过程不再赘述，可参阅相应文献.

试验加载至17 280 kN并维持加载约5 min时，桩身沉降突然增加，地下同时发出一声很闷的巨响，之后试桩偏位，地面出现裂缝，桩身沉降无法稳定，可判定桩身折断，其屈曲的临界荷载可取16 720 kN.

1.3 有限元计算

使用有限元软件ABAQUS对该现场试验进行建模分析.模型采用三维可变形的C3D8R单元建立，桩基本构模型采用弹性模型，弹性模量E设为3.1 GPa；土层的本构模型采用摩尔库仑模型，具体参数取值如表1所示.桩土之间设置接触面，法向行为设置为硬接触、切向行为设置为罚函数，摩擦系数设为0.5.土体四周及底面设置位移边界条件以约束土体位移，桩底边界条件为嵌固，桩端自由.同时为提高计算精度，将桩与桩周土体进行进一步的加密划分，模型建立如图1所示.

模型的材料属性、网格、接触属性、边界条件建立后模拟现场试验的各级加载，创建静力学功能下的geo分析步对地应力进行平衡，之后创建静力通用分析步，将荷载分级加于桩顶，得到各级荷载下桩身轴力，绘制深度-轴力曲线，并将现场试验实测数据以散点形式在图上对应位置标出，如图2所示.图2曲线显示，现场试验实测数据点基本落于深度-轴力曲线上，说明了模型参数的正确性以及模型计算的可行性与准确性.

然后分析该地层条件下桩的临界荷载，创建线性摄动功能下的buckle分析步计算桩的临界荷载，计算得其临界荷载P_cr=17 006 kN，桩发生屈曲破坏时，桩、土模态如图3所示.

计算表明，采用有限元方法计算结果与现场试验实测值接近，说明有限元方法计算桩基屈曲临界荷载是可靠的.

1.4 按规范简化公式计算

常用的简化计算方法^［7］（下文称：常用简化公式）过程如下：首先按式（3）确定计算长度l_p的取值.

l_{\mathrm{p}}=\mathrm{0.7}(l_{\mathrm{0}}+\frac{\mathrm{4}}{\alpha })（3）

式中：l_{\mathrm{0}}为柱高， m； α为桩土变形系数，m^-1.

然后将l_p代入式（4），即可得到规范法下的临界荷载值.

P_{\mathrm{c}\mathrm{r}}=\frac{{\mathrm{\pi }}^{\mathrm{2}}EI}{{l_{\mathrm{p}}}^{\mathrm{2}}}（4）

式中：E为桩弹性模量，Pa； I为桩横截面极惯性矩，m⁴.

对于茅草街大桥试验基桩，根据式（3）可得l_p=7.33 m，再根据式（4）可得临界荷载值P_cr=223 610 kN，这与实测值（P_cr=16 720 kN）相差很大，表明常用简化公式计算结果并不准确，这是因为该公式假定桩基在4/α深度处已经嵌固，故将计算长度定为柱高与嵌固深度之和再乘边界条件系数，桩嵌固深度以下的部分不受影响.事实上，不同种类土的力学性质差异很大，这一假定显然是不严谨的，因此本文将对该方法的适用范围进行研究，然后提出新的简化计算方法.

2 常用简化公式适用范围分析

2.1 数值实验方案

一般来说，深厚软土区的地层分布特点为：上层的软土层（包括淤泥、淤泥质土、流塑性黏土、软塑黏性土等），中层的硬土层（包括砂层、粉砂层、硬塑黏土层、卵石层等），下层的岩石层.图4所示地层可代表典型深厚软土区的地层.本文基于该图所示地层研究典型软土区桥梁桩基临界荷载变化规律.

考虑对临界荷载有影响的土层水平抗力系数m值、桩径d、桩长l、柱高l₀ 这4种因素进行数值实验.其中，桩长l由软土层中桩长l₁与硬土层中桩长l₂两部分组成（l=l₁+l₂）；m由软土层中m₁与硬土层中m₂两部分组成.改变各参数的取值水平，分别采用常用简化公式法和基于特征值分析的有限元法计算各参数组合下桩柱式桥墩临界荷载值，模型参数示意图如图4所示.模型中各参数取值范围如表2所示，计算时将各个因素取不同水平组合进行计算.

由于有限元软件中没有直接定义土的水平抗力系数的参数，故需要将土层水平抗力系数m值与弹性模量E值进行换算，其方法为将有限元软件各弹性模量下桩的水平载荷试验产生的桩顶水平位移与规范法各m值下计算的桩顶水平位移进行对比，将计算结果非常接近的弹性模量E与m值对应起来，形成若干组数据，运用最小二乘法进行回归分析，以得到弹性模量E与m值的对应关系^［11］，本文模型将桩周土层简化为单层土，这样能将土层参数与某一确定m值对应起来，本文中m值与弹性模量E的对应关系如表3所示.

2.2 计算结果分析

将相同参数组合下，两种不同方法计算所得临界荷载绘制成散点图（规范法P_crs‒有限元法P_cra），如图5所示（长细比=l₁/d+l₂/d）.

计算结果显示，长细比为7.5~20时，两种方法的计算结果接近；长细比为20~45时，两种方法的计算结果已经较为接近；长细比为60~45时，两种方法的计算结果偏差很大.这说明随着桩长细比的增大，规范法计算结果与有限元计算的结果偏差越来越大，故需要对规范法的计算公式进行修正.

3 临界荷载简化计算公式

3.1 随机森林算法简介

近年来，人工智能发展迅速，机器学习逐渐应用到越来越多的专业领域，其中随机森林算法（Random Forest）对于数据的适应能力强，实现方法较为简便，被广泛地应用于数据处理当中.随机森林是一种可用于回归分析的非参数集成学习算法^［12］，起源于对于决策树的研究^［13-14］，通过随机选择训练样本与特征变量的方法生成多个CART决策树^［15］，之后通过投票计分的方式将这些决策树的训练结果结合起来，以获得最终的回归模型，训练好的模型将可以通过输入自变量参数，较为准确地预测出因变量参数.

评价随机森林模型的回归效果，可采用平均绝对值误差（MAE）、平均均方根误差（RMSE）和拟合优度（R²）这3个评价指标^［13］，本文采取拟合优度（R²）指标进行评价.评估方法采用“留出法”，将用有限元方法计算出来的数据划分为两个数据集合，把其中一个集合作为训练集S，另一个作为测试集T，在训练集S上训练好模型之后，再在测试集T上评估其预测精度，若两个集合的预测精度均在0.99以上，则说明模型训练效果较好，模型预测精度较高.

3.2 随机森林模型预测结果

为得到深厚软土区桩柱式桥墩计算长度的简化公式，以第2节中的数值模型的计算结果为基础，将结果换算为计算长度，与输入参数构成数据集.随机选取数据集66.7%的数据作为训练集S，33.3%的数据作为测试集T^［16］，以长细比、柱高、桩土变形系数α值、土层水平抗力系数m值作为输入特征变量，采用随机森林回归算法对计算长度进行预测.将随机森林模型得到的预测值与有限元法的计算值绘制散点图，如图6所示.

结果显示，训练集与测试集的计算结果均处于 1∶1线附近，测试集的拟合优度也达到了0.99以上，深度学习的效果非常好，证明输入参数与输出参数之间具有很强的相关性，可使用公式的形式对其进行拟合.

3.3 基于随机森林模型的计算长度简化公式

对于随机森林模型，根据Gini指数计算特征变量，可分析各个属性的重要性^［15］，以各因素的重要性指数为基础可以确定拟合公式中参数的形式从而导出公式.计算本文随机森林模型中各特征的重要性，并绘制柱状图如图7所示.

计算结果显示，在深厚软土区的桩柱式桥墩的稳定性问题中，桩的长细比对于计算长度的影响占主导地位，其中软土层中桩长细比l₁/d占44.79%，硬土层中桩长细比l₂/d占24.76%，合计69.55%.m值的影响次之，其中软土层m₁占15.15%，硬土层m₂占5.19%，合计20.34%.柱高l₀对计算长度也有一定的影响，占10.11%；因此以式（3）为基本形式，纳入l/d参数可得到更准确的简化计算公式.本文采用不同形式的公式进行分析，得到桩柱式桥墩计算长度简化公式如（5）所示，将拟合结果绘制在图8中，可见计算较好.

l_{\mathrm{p}}=\mathrm{0.7}(l_{\mathrm{0}}+\frac{\mathrm{4}}{\alpha }{\mathrm{e}}^{\mathrm{0.046}·\frac{l_{\mathrm{1}}}{d}+\mathrm{0.037}·\frac{l_{\mathrm{2}}}{d}-\mathrm{0.68}})（5）

根据式（5）计算各参数水平下计算长度，然后采用式（4）计算桩柱式桥墩的临界荷载，再与有限元法计算值比较，绘制散点图（2 520个数据点）如图9所示.

计算结果表明，本文提出的简化计算方法效果良好，具有实用价值.

4 结 论

本文分别采用基于特征值屈曲分析的有限元数值软件法、规范法，对桩柱式桥墩的屈曲临界荷载进行计算和比较，并使用随机森林法对有限元计算结果进行分析，得到结论如下：

1）采用有限元数值软件方法计算桩柱式桥墩的临界荷载，其计算结果与现场试验结果较为吻合，说明基于特征值屈曲分析的有限元法对桩柱式桥墩的临界荷载计算是适用的.

2）现有常用简化计算方法在桩的长细比较大、时计算结果误差较大，适用性差.

3）基于有限元计算以及随机森林回归分析，拟合得到了深厚软土区桩柱式桥墩的临界荷载计算简化方法，可供设计人员使用与参考.

参考文献

[1] 赵明华 ．桥梁桩基的屈曲分析及试验［J］．中国公路学报，1990， 3（4）： 47-56．

[2] ZHAO M H ． Buckling analysis and tests of bridge piles［J］．China Journal of Highway and Transport， 1990， 3（4）： 47-56． （in Chinese）

[3] 赵明华， 王季柏 ．摩擦桩的屈曲分析［J］．湖南大学学报（自然科学版）， 1995， 22（6）： 109-114．

[4] ZHAO M H， WANG J B ．Buckling analysis of friction piles ［J］．Journal of Hunan University （Natural Sciences），1995，22（6）：109-114．（in Chinese）

[5] 赵明华， 汪优， 黄靓 ．基桩屈曲的几何非线性有限元分析［J］．岩土力学， 2005，2 6（）： 184-188．

[6] ZHAO M H， WANG Y， HUANG L ．Geometrically nonlinear finite element analysis of pile buckling［J］．Rock and Soil Mechanics，2005，26（Sup.2）：184-188．（in Chinese）

[7] 彭锡鼎 ．考虑桩侧土壤弹性抗力时桩的临界荷载计算［J］．土木工程学报，1996，29（5）：43-48．

[8] PENG X D ． Calculation of critical load of pile considering resistance of soil on the side of pile［J］．China Civil Engineering Journal，1996，29（5）：43-48．（in Chinese）

[9] 邹新军，赵明华，刘光栋 ．基于无单元Galerkin法的软土区高承台超长嵌岩基桩屈曲分析［J］．工程力学，2009，26（1）：125-130．

[10] ZOU X J， ZHAO M H， LIU G D ． Buckling analysis of super-long rock-socketed filling piles in soft soil area by element free Galerkin method［J］．Engineering Mechanics，2009，26（1）：125-130．（in Chinese）

[11] 雷勇，李鹏甲，刘泽宇，等 ．岩溶区穿越溶洞基桩屈曲临界荷载计算方法与试验研究［J］．岩土力学，2022，43（12）：3347-3356．

[12] LEI Y， LI P J， LIU Z Y，et al ．Method for calculation of buckling critical load of pile foundation crossing Karst cave in Karst area［J］．Rock and Soil Mechanics，2022，43（12）：3347-3356．（in Chinese）

[13] 铁道部第三勘测设计院. 铁路工程设计技术手册-桥梁设计通用资料［M］. 北京： 人民铁道出版社， 1983.

[14] The Third Survey and Design. Technical manual for railway engineering design-general information for bridge design ［M］. Beijing： People’s Railway Publishing House， 1983. （in Chinese）

[15] 建筑桩基技术规范： JGJ 94—2008 ［S］．北京： 中国建筑工业出版社， 2008．

[16] Technical code for building pile foundations：JGJ 94—2008 ［S］．Beijing：China Architecture & Building Press，2008．（in Chinese）

[17] 王凤，王文达，史艳莉 ．钢管混凝土框架柱计算长度研究［J］．工程力学， 2015， 32（1）： 168-175．

[18] WANG F， WANG W D， SHI Y L ． Research on the effective length of concrete-filed steel tube frame columns［J］．Engineering Mechanics，2015， 32（1）： 168-175．（in Chinese）

[19] 刘齐建， 赵明华， 李俞， 等 ．茅草街大桥桩基试验研究［J］．湖南大学学报（自然科学版）， 2004， 31（4）： 51-54．

[20] LIU Q J， ZHAO M H， LI Y，et al ．In-site load test of the piles in the foundation of Maocaojie Bridge［J］．Journal of Hunan University （Natural Sciences），2004，31（4）：51-54．（in Chinese）

[21] 周昌林， 朱爱军. 地基土水平抗力系数的比例系数m与其弹性模量E关系的数值试验研究［J］. 建筑科学， 2015， 31（）：127-131.

[22] ZHOU C L， ZHU A J. Numerical test research of elastic modulus with proportional coefficient of horizontal resistance coefficient of soil［J］. 2015， Building Science， 2015，31（Sup.1）：127-131. （in Chinese）

[23] BREIMAN L. Random forests［J］. Machine Learning，2001，451：5-32.

[24] QUINLAN J R ．Induction of decision trees［J］．Machine Learning， 1986， 1（1）： 81-106．

[25] QUINLAN J R ．Probabilistic decision trees［M］//Machine Learning．Amsterdam：Elsevier，1990：140-152．

[26] ABELLÁN J， MANTAS C J， CASTELLANO J G ．A random forest approach using imprecise probabilities［J］．Knowledge-Based Systems，2017，134：72-84．

[27] RUDŽIANSKAITĖ-KVARACIEJIENĖ R， APANAVIČIENĖ R， GELŽINIS A ．Modelling the effectiveness of PPP road infrastructure projects by applying random forests［J］．Journal of Civil Engineering and Management，2015， 21（3）： 290-299．

基金资助

国家自然科学基金资助项目(52178331)