主页 › 范文 › 地质学 › 本文

地层结构对能源桩热力学响应的影响

摘要

文章通过解析解与实验数据对比，研究地层结构对能源桩热力学行为的影响。基于弹性土-桩界面假设，建立均质与四层土壤的解析模型，预测桩头位移、轴向应变和应力。离心机试验（均质土壤）与瑞士现场试验（分层土壤）验证表明：均质模型能准确捕捉整体响应趋势，而分层模型可识别层间界面处的应力/应变突变。研究发现，相邻土层刚度差异（如B/C层界面）显著影响桩身热应力分布，刚度突变导致层边界处出现极值点。

1. 引言

能源桩作为低碳建筑关键技术，通过桩基结构实现地热交换，但其热变形受周围土层约束产生附加应力。文章聚焦地层结构复杂性对土-桩相互作用的影响，基于Cossel（2019）与Perić（2020）的解析解，对比均质与分层土壤中能源桩的热力学响应。研究通过离心机试验（均质Bonny粉土）和瑞士洛桑现场试验（四层粘土）验证模型精度。

2. 解析模型构建

2.1 基本控制方程

考虑热-机械荷载的能源桩满足以下控制方程：

运动学关系：\(\varepsilon = \frac{du}{dx}\)
热弹性本构：\(\sigma = E(\varepsilon - \alpha \Delta T)\)
土-桩界面剪应力：\(|\tau| = k_s |u|\)

推导得控制微分方程：

\[ \frac{d^2 u}{dx^2} - \psi^2 u = 0, \quad \psi^2 = \frac{p k_s}{A E} \]

其中\(\psi\)为特征参数，\(k_s\)为土层剪切刚度。

2.2 边界条件与解

半浮式能源桩边界条件：

\[ \sigma(0) = k_b u(0), \quad \sigma(L) = \frac{F}{A} \]

解的形式为双曲函数组合，零位移点（null point）位置由下式确定：

\[ x_0 = \frac{1}{\psi} \tanh^{-1} \left[ \frac{\cosh(\psi L) - 1}{\sinh(\psi L) + \frac{k_b}{E\psi}} \right] \]

端承桩特例（\(k_b \to \infty\)）下，零位移点位于桩端（\(x_0=0\)），解简化为：

\[ \sigma(x) = E\alpha\Delta T \left[ \frac{\cosh(\psi x)}{\cosh(\psi L)} - 1 \right] + \frac{F \cosh(\psi x)}{A \cosh(\psi L)} \]

3. 地层结构影响分析

3.1 均质土壤验证

Stewart离心机试验（端承桩，\(\Delta T=13-40^\circ C\)）显示：

解析解精确预测轴向应变/应力分布（图2），最大误差<5%。
桩端位移近似为零的假设在硬质基岩条件下成立。

3.2 分层土壤响应

瑞士现场四层粘土（A1/A2/B/C）试验表明：

分层模型预测桩头位移与实测数据吻合最佳（图4），均质A1/B土壤预测次之。
层间界面处出现响应突变：B/C层边界因\(k_s\)差异最大（18.2→121.4 MPa/m），产生应变极值点（图5）。
均质模型无法捕捉层间不连续性，导致桩端应力预测偏差达25%（图6）。

3.3 热-机械耦合效应

机械荷载（F=-1000 kN）与热荷载（\(\Delta T=14^\circ C\)）共同作用下：

分层模型在C层预测应变递增/应力递减趋势（图7-8），与实测定性一致。
桩端弹簧刚度\(k_b=6675\) MPa/m对纯热荷载（T1）优化，但机械荷载（T7）下需调整。

4. 结论

地层结构通过剪切刚度\(k_s\)的空间变化显著影响能源桩热力学响应：

均质模型适用于单一土层或刚度渐变地层，可高精度预测整体趋势。
分层模型能捕捉层间界面处的应力/应变突变，B/C类高刚度差异界面需重点建模。
桩端约束条件对响应分布形态敏感，硬质基岩假设在深部地层需进一步验证。