考虑路基弹塑性特性的条形超载下水平附加应力修正方法研究

《Buildings》：Study on the Modification Method of Horizontal Additional Stress Under Strip Surcharge Considering Elastoplastic Characteristics of the Subgrade

　　

针对沿海软土地基中条形超载引起侧向挤压并危及邻近既有桥桩安全的问题，传统布辛尼斯克（Boussinesq）弹性理论无法反映土体的真实弹塑性特性，易低估实际水平附加应力。本文建立了一个二维平面应变有限元模型，并在纯弹性理论解标定的基础上，于弹塑性地基条件下开展了广泛的对比分析。通过皮尔逊（Pearson）相关性和随机森林（Random Forest）敏感性分析，明确了内摩擦角、荷载比和归一化距离比是影响水平附加应力重分布的核心控制变量，从而证明了弹性模量和黏聚力影响的局限性。研究揭示了地基深层塑性区贯入所引起的水平应力传递的非线性放大机制，以及应力修正系数呈现先高斯（Gaussian）峰值增强、随后随空间距离呈指数衰减的物理演化规律。基于这些机制，提出了水平附加应力修正系数的组合预测公式，其决定系数R²达到0.903，并构建了修正系数的安全评价图以量化高风险区域。结果表明，当归一化距离比大于或等于4时，弹塑性挤压效应基本消散。所提出的修正方法有效界定了弹性解的适用边界，为软土地区被动受荷桩基的承载力计算与安全控制提供了理论依据。

### **沿海软土地区条形超载下水平附加应力修正方法研究解读**

#### **一、 研究背景与问题**

在沿海地区，桥梁作为关键的交通生命线，其安全至关重要。然而，随着沿海地区大规模填海造陆与城市建设的不断推进，一系列复杂的岩土工程灾害随之而来。其中，软土地基中邻近既有桥桩的回填或堆载超载问题尤为突出。地表附加超载破坏了土体的初始应力平衡，迫使高压缩性的软土产生显著的侧向挤压变形。这种侧向位移的土体直接推挤邻近的桥桩基础，使桩基处于被动受荷状态，进而可能导致桩身产生过大的弯曲变形甚至剪切破坏风险，严重威胁桥梁的安全运营。国内外因不当堆载引发的桩基失效与工程灾害屡见报道，例如英国的瑟蒂斯桥（Surtees Bridge）以及上海某工业仓库倒塌和“莲花河畔”楼房倒塌等事故，其核心原因均与附加超载触发土体弹塑性破坏及其对邻近桩基产生的巨大侧向推力密切相关。

针对此类超载下被动桩基的受力特性研究，主要存在两大难点：一是桩身所受被动推力的确定，二是超载下水平附加应力的计算。目前，国内外对于水平附加应力的计算主要依赖于布辛尼斯克（Boussinesq）弹性理论体系，并已发展出成熟的分析方法，为构建更精确的被动桩理论计算模型提供了可靠的应力边界条件。然而，上述研究均基于弹性理论体系，当超载荷载较小或地基土质较好时近似可行，但对于软土地基则存在显著的应用风险。其根本原因在于软土地基的弹塑性特性与弹性体系差异巨大，这也是软土地区超载导致桥桩基础失效工程事故频发的原因。尽管先前研究为被动桩受荷分析奠定了坚实基础，但它们往往侧重于特定案例下的宏观结构响应，且严重依赖定制化的有限元分析或复杂的经验法则。在桩基设计之前，对于水平应力场重分布的直接、通用的基线定义仍存在明确的知识空白。本研究通过解耦条形分布超载下的自由场水平附加应力修正系数（ξ），填补了这一空白，使设计工程师能够直接从传统的布辛尼斯克解中即时评估弹塑性应力放大的边界，而无需进行大规模定制化的数值运算。

#### **二、 主要技术方法**

为克服传统弹性理论体系无法考虑地基真实弹塑性变形特性带来的不确定性影响，研究人员首先建立了弹性有限元模型以与弹性理论解进行标定验证。随后，综合考虑地基土弹性模量E、黏聚力c、内摩擦角φ、距离比s/B及地基承载力荷载比R等因素的影响，在弹性与弹塑性有限元地基模型之间，对条形分布荷载下的水平附加应力进行了广泛的对比分析。此过程旨在提取一个考虑地基弹塑性特性影响的水平附加应力修正系数ξ的计算方法。作为一个基础性方法论研究，本研究在受控的连续介质框架内隔离了纯力学机制，以建立一个适用于广泛岩土材料和地基土在条形超载下的通用基线，从而侧重于针对经典解析解进行严格的理论验证，而非针对特定场地的经验性现场报告。

具体技术方法包括：1. **有限元模型建立与验证**：采用ABAQUS软件建立二维平面应变模型，模型尺寸、边界条件（两侧约束水平位移，底部约束全部位移）及网格划分（CPS4R单元）经过精心设计。通过将纯弹性有限元结果与经典弹性理论解（Boussinesq解）进行对比，验证了有限元建模方案的正确性，确认弹性有限元结果可完全替代理论解析解，为后续对比分析提供了可靠基准。2. **弹塑性本构与参数设定**：弹塑性模型统一采用莫尔-库仑（Mohr-Coulomb）本构模型，并考虑非关联流动法则以确保计算收敛。模型参数（如E, c, φ）根据《工程地基手册》中软土至硬土的典型分布范围选取宽泛变量，以保障结论的普适性。3. **正交试验设计与批量计算**：采用正交试验设计方法组合材料参数，结合多个荷载比R水平，生成了总计100组的有限元模拟工况进行批量计算。利用维西克（Vesic）公式计算地基极限承载力p_u，进而确定实际加载值p。4. **数据提取与敏感性分析**：从每个模型中提取不同s/B值下的水平附加应力数据。利用皮尔逊（Pearson）相关系数和随机森林（Random Forest）算法对多元数据进行参数敏感性分析，量化各参数对修正系数ξ的影响权重和方向。5. **修正系数模型构建**：基于敏感性分析结果，聚焦于内摩擦角φ、荷载比R和距离比s/B三个核心变量。针对ξ > 1的高风险区间，采用高斯函数与指数函数叠加的组合公式结构进行高精度拟合，建立了修正系数ξ的预测模型，并评估其拟合精度（R², RMSE）。

#### **三、 研究结果**

**1. 有限元方法对弹性理论解的验证**
研究首先验证了纯弹性有限元模型的准确性。将有限元计算的竖向和水平向附加应力结果与经典弹性理论解进行对比，结果显示两者基本一致。这证实了本研究所采用的有限元建模、边界约束和网格划分方案的正确性，并且弹性有限元模型的结果可以完全替代理论解析解，为后续考虑弹塑性特性的附加应力修正分析提供了可靠的分析基准。

**2. 条形超载下水平附加应力修正方法研究**
**（1）参数敏感性分析**
通过对批量有限元计算得到的多元数据进行相关性分析和随机森林重要性评估，明确了影响水平附加应力修正系数ξ的核心因素：
* 内摩擦角φ是影响修正系数的核心控制因素，其皮尔逊相关系数高达-0.83，随机森林特征重要性为0.7027，远超其他变量。φ越大，土体的非线性应力重分布效应越强，导致ξ显著减小。
* 荷载比R具有中等强度的影响，与ξ呈负相关。R越大（荷载越接近地基极限承载力），塑性区范围越大，实际应力与弹性解的偏差也越大。
* 距离比s/B表现出空间衰减效应。当s/B超过3后，ξ基本趋于稳定并回归至1附近，表明修正系数具有显著的空间局限性。
* 弹性模量E和黏聚力c的影响可忽略。E的绝对值在应力修正比ξ的除法形式中基本被抵消；c的影响极其有限，对工程应用的实际修正意义不大。
因此，在构建修正系数ξ的预测模型时，只需关注内摩擦角φ、荷载比R和距离比s/B这三个核心自变量。

**（2）水平附加应力修正系数**
研究提出了一个针对高风险区间（ξ > 1）的水平附加应力修正系数预测公式。该公式采用高斯函数与指数函数叠加的组合结构，以模拟ξ随距离先快速达到峰值后衰减收敛的规律。通过全局多元回归分析得到的最终预测公式决定系数R² = 0.903，均方根误差RMSE = 0.0407，在工程高度重视的ξ > 1风险区域实现了高精度预测，可用于可靠评估传统弹性解在特定条件下的不安全性。

**3. 变量影响规律与修正系数图表**
**（1）变量影响规律研究**
通过对比验证，表明预测公式在不同工况下总体预测效果良好，能可靠捕捉修正系数ξ > 1的危险区域。研究发现：
* 修正系数ξ随s/B的变化规律与预测公式形态一致：先经历高斯增强阶段，随后进入指数衰减阶段。当s/B ≥ 4时，修正系数ξ值基本趋于1，即超过此范围后无需修正，可直接使用经典布辛尼斯克弹性解。这一规律与Zhou等人的离心机试验观测到的非线性挤压与消散现象相印证。
* 荷载比R增加会导致曲线整体上抬和拓宽，表明荷载强度越高，弹塑性条件下的水平附加应力值越大，且大于弹性解的区域也逐渐扩大。
* 内摩擦角φ对ξ的影响复杂：在s/B = 0处，ξ随φ增大而减小；在s/B = 2~8范围内，ξ随φ增大而增大，但增幅渐缓。

**（2）修正系数图表及其意义**
基于预测公式，研究构建了条形分布荷载下水平附加应力修正系数ξ图表。图中ξ = 1的等值线为分界线，其上方的区域（ξ > 1）为需要修正的危险区，下方（ξ < 1）为安全区。该图表的核心功能在于能够快速判断当前工程参数组合是否需要进行水平附加应力系数修正。例如，当图表显示工程处于安全区时，则无需考虑弹塑性影响下的修正；反之，则需按预测公式计算修正系数值，对弹性解进行修正。图表还揭示了不同s/B值下安全区与危险区分布形态的差异，为工程加固处理提供了直观的定量指导。

#### **四、 讨论与结论**

本研究针对沿海软土地基中条形超载引起的水平附加应力计算问题，通过建立二维平面应变有限元模型，对比分析了弹性和弹塑性地基的力学特性，提出了一种考虑地基弹塑性特性的水平附加应力修正方法。主要结论如下：

1. 传统弹性理论解难以反映软土的真实弹塑性特性，在修正系数ξ > 1的高风险工程设计条件下，会显著低估实际水平附加应力。参数敏感性分析表明，内摩擦角φ、荷载比R和归一化距离比s/B是影响水平附加应力重分布的核心控制变量，而弹性模量E和黏聚力c对修正系数的影响极弱，在计算中予以排除。
2. 地基深层塑性区的贯入对水平应力传递具有非线性放大机制。水平附加应力修正系数ξ沿空间距离呈现非单调演化规律：先呈现高斯峰值增强，随后呈指数衰减。当空间距离达到s/B ≥ 4时，弹塑性挤压效应基本消散，修正系数ξ趋近于1.0，表明超出此范围后可直接采用经典布辛尼斯克弹性解。
3. 基于上述变量规律，本文提出了一个高斯函数与指数函数叠加的组合预测公式。该公式在ξ > 1的关键区间具有较高的预测精度，R² = 0.903，RMSE = 0.0407。基于此构建的修正系数图表可直观界定并量化需要应力修正的危险区域边界，为软土地区被动受荷桩基的承载力计算与安全处治提供了直接的理论与定量依据。