软弱土层下预制拼装综合管廊受力性能分析
0 引言
综合管廊是指将两种以上的城市管线集中设置于同一地下空间中, 所形成的一种现代化、集约化的城市基础设施
综合管廊本质上属于长线形隧道结构。因纵向地基土层一般无法保证在综合管廊整个纵向长度范围内均为匀质土层, 而纵向地基土层不均匀将导致综合管廊出现纵向附加内力, 从而威胁结构安全。因此, 对综合管廊在软弱土层中受力性能开展系统的有限元参数分析, 对于进一步了解综合管廊结构的受力特点、明确其结构安全性具有重要意义。
目前, 国内外有关综合管廊在土压力作用下的受力性能有限元分析成果较少, 而对预制拼装综合管廊的相关研究则几乎为空白。但针对混凝土地下箱涵结构, 有学者开展了一些探索性研究。段旻罡
鉴于此, 本文采用大型非线性有限元软件ANSYS, 考虑预制拼缝界面特性以及土-结构共同工作, 对预制拼装混凝土综合管廊在软弱土层下的受力性能进行系统分析, 从而为预制拼装综合管廊的推广应用及相关技术标准的编制提供参考。
1 有限元建模
本文建立的非线性有限元分析模型以上海市某综合管廊工程为背景。采用如图1所示整舱预制拼装方案, 纵向预制节段之间通过螺栓进行连接。
1.1 单元类型及材料
1.1.1 单元选择
由于综合管廊结构外形对称, 为提高计算效率, 基于对称约束, 本文仅建立综合管廊的半结构模型进行分析。综合管廊结构的混凝土底板、侧壁及土体均采用八节点四面体实体单元solid65, 钢筋采用线单元link8, 连接螺栓采用link10单元。本文建立的基于ANSYS软件的有限元模型如图2所示。
1.1.2 材料的基本力学特性
预制拼装综合管廊分析模型中, 混凝土强度等级为C40, 钢筋强度等级为HRB335。为研究软弱土层弹性模量对预制拼装综合管廊结构受力性能的影响, 本文选择了5种典型软弱土层的弹性模量, 非软弱土层地基弹性模量取30MPa。

图2 预制拼装综合管廊及周围土层有限元分析模型Fig.2 Finite element analysis model of precast assembled utility tunnel structure and subsoil
有限元模型中, 混凝土的受压应力-应变本构采用曲线加直线形式, 当混凝土应变εc≤ε0 (最大应力时的混凝土压应变) 时, 应力-应变曲线关系为:

当应变ε0<εc<εcu时, σc=f0c。εc为0.002με, εcu为0.003 3με, 因此取极限压应变值0.003 3με。
土的本构模型采用Drucker-Prager准则, 考虑了静水压力对屈服的影响。根据上海市土质情况, 本文取内摩擦角为27.9°, 黏聚力为24.4k Pa。
1.2 关键问题处理
1.2.1 边界条件
根据周围土层特点及地下结构有限元分析的常用处理方式, 本文将计算土层底面设置为竖向固定约束, 两侧设置为水平固定约束, 对称面则设置为ANSYS系统给定的对称边界条件。本文通过整体受力并协调变形的结构模型来模拟管廊和土体的共同作用, 即不考虑土与结构之间的相互滑动和错开。该方法是目前地下结构分析中普遍采用的方法, 能够保证足够的计算精度
1.2.2 预制拼缝界面
各预制分块之间的接触界面采用ANSYS中的面-面接触单元来模拟, 其中一个面为“目标面”, 另一个面为“接触面”, 分别用targe170和conta173来做接触分析。
2 试验验证
为验证本文建模方法的合理性, 基于文献
由表1, 2分析可见, 虽然本文模型的计算结果与文献中的实测结果误差较大, 但二者变化趋势接近, 且本文模型的计算结果与文献中的计算结果较为接近。这主要是由于有限元模型一般难以完整地模拟实际情况。总体而言, 本文建模方法建立的有限元分析模型能够较好地反映地下结构考虑土-结构共同工作影响的受力性能变化规律。
表2 混凝土应变计算值与实测值对比Table 2 Comparison between calculated values and measured values of concrete strain

3 参数分析
以上海市某综合管廊工程为原型 (见图1) , 建立了预制拼装综合管廊非线性有限元分析模型, 如图4所示。为充分反映地基土层的影响, 有限元模型中土层模型的长、宽、高分别为50, 17.3, 20.8m。有限元模型中预制拼装综合管廊结构的上覆土深度为3m, 上部路面的车辆荷载产生大小为40k Pa的面荷载, 均与实际工程一致。
本文拟重点分析软弱土层对预制拼装综合管廊受力性能的影响规律, 因此在有限元模型中部设置了局部软弱土层。有限元分析参数为软弱土层弹性模量。不同软弱土层弹性模量下预制拼装综合管廊受力性能的有限元分析参数如表3所示。
表3 软弱土层弹性模量影响有限元分析参数Table 3 Parameters of finite element analysis on influence of modulus of soft soil

3.1 沉降变形
预制拼装综合管廊纵向沉降变形随软弱土层弹性模量变化的关系曲线如图5所示, 预制拼装综合管廊纵向不同位置处沉降变形随软弱土层弹性模量比 (软弱土层弹性模量与周围非软弱土层弹性模量的比值) 的变化关系曲线如图6所示。由图可见:

图5 纵向沉降变形随软弱土层弹性模量变化的关系曲线Fig.5 The relationship curves between settlement and elasticity modulus of soft subsoil
1) 随着软弱土层弹性模量的减小, 预制拼装综合管廊的不均匀沉降明显加剧。
2) 由图6分析可知, 沿预制拼装综合管廊纵向不同位置处的沉降变形与软弱土层弹性模量的关系函数可近似简化为对数函数y=alnx+b, 其中y为沉降变形量, x为软弱土层弹性模量比, a, b均为计算常数, 与荷载情况及综合管廊结构形式有关。
3.2 纵横向内力
当存在软弱土层时, 预制拼装综合管廊结构内部将产生纵向弯矩以及横向附加弯矩。不同软弱土层弹性模量下管廊纵向不同位置处横截面方向侧壁底部截面附加弯矩随软弱土层弹性模量变化的关系曲线如图7所示, 管廊不同位置处纵向弯矩随软弱土层弹性模量变化的关系曲线如图8所示。

图6 纵向典型位置处的沉降变形随软弱土层弹性模量比变化的关系曲线Fig.6 The relationship curves between settlement of typical section and elasticity modulus of soft subsoil

图7 各横向截面侧壁底部截面附加弯矩随软弱土层弹性模量变化关系曲线Fig.7 The relationship curves between increased section moment at different positions in longitudinal direction and elasticity modulus of soft subsoil

图8 不同位置处纵向弯矩随软弱土层弹性模量关系曲线Fig.8 The relationship curves between increased longitudinal moment at different positions and elasticity modulus of soft subsoil
由图7, 8所示曲线分析可知:
1) 软弱土层对综合管廊横截面方向的弯矩有较大影响。其中, 纵向中间区域所受的影响最大, 使其侧壁底部截面的弯矩最大减小55.1k N·m, 约占无软弱土层下该截面弯矩的31.9%。此外, 在综合管廊纵向20m处, 软弱土层对横截面方向的弯矩也具有较大影响, 使其侧壁底部截面的弯矩最大增加33.8k N·m, 约占无软弱土层下该截面弯矩的19.7%。因此, 在预制拼装综合管廊结构设计时应关注软弱土层对横截面方向结构内力的影响。
2) 软弱土层对预制拼装综合管廊纵向内力有一定影响。影响最大的位置在综合管廊纵向中部, 即软弱土层中部。该位置处的截面纵向弯矩随着软弱土层弹性模量的降低而不断增大。当软弱土层的弹性模量降低为5MPa (约为周围非软弱土层弹性模量的17%) 时, 纵向弯矩最大达3 500k N·m, 约为该截面开裂弯矩的25%。
4 结语
本文基于ANSYS软件建立了考虑预制拼缝界面和土-结构共同工作的预制拼装综合管廊非线性有限元分析模型, 并通过已有文献的实测值与计算值验证该模型的合理性。在此基础上, 开展软弱土层下预制拼装综合管廊受力性能的有限元参数分析, 重点分析了软弱土层弹性模量的影响规律。
1) 随着软弱土层弹性模量的减小, 预制拼装综合管廊的不均匀沉降明显加剧。预制拼装综合管廊纵向不同位置处的沉降变形与软弱土层弹性模量的关系函数可近似简化为对数函数。
2) 软弱土层对综合管廊横截面方向的弯矩有较大影响。对于纵向中间区域, 侧壁底部截面的弯矩最大可减小31.9%;对于纵向20m处, 侧壁底部截面的弯矩最大可增加19.7%。因此, 在预制拼装综合管廊结构设计时应关注软弱土层对横截面方向结构内力的影响。
3) 软弱土层对预制拼装综合管廊纵向内力有一定影响。对于管廊纵向中部, 当软弱土层的弹性模量降为周围非软弱土层弹性模量的17%时, 纵向弯矩最大可达到该截面开裂弯矩的25%。
参考文献
[1]胡翔, 薛伟辰.预制拼装综合管廊受力性能试验研究[J].土木工程学报, 2010, 43 (5) :29-37.
[2]薛伟辰, 胡翔, 王恒栋.上海世博园区预制拼装综合管廊力学性能试验研究[J].特种结构, 2009, 26 (1) :105-108.
[3]段旻罡.钢筋混凝土箱涵非线性有限元分析技术研究[D].天津:天津大学, 2004.
[4]刘蒙娜.地下箱涵有限元结构分析[D].郑州:郑州大学, 2014.
[5]鞠杨, 徐广泉, 谢和平.钢筋混凝土盾构衬砌结构的三维数值分析[J].计算力学学报, 2005, 22 (4) :437-442.
[6]李元松, 李新平, 代翼飞, 等.隧道围岩与衬砌受力特性测试与数值分析[J].岩土力学, 2007, 28 (7) :1348-1352.