真实复杂地层大直径钢顶管三维数值模拟研究.docx

PAGE

PAGE

真实复杂地层大直径钢顶管三维数值模拟研究

前言

随着城市建设的发展，地面建（构）筑物越来越密集，地上交通越来越繁忙，在城市建成区进行市政管线建设时采用明挖法施工的可能性逐渐减少，取而代之的是采用对地面交通影响较小的非开挖法施工，其中顶管法是一种典型的非开挖市政管线施工方法，应用逐步增多。因此，顶管法施工对周围土体的扰动及地表沉降的影响，愈来愈受学术界和工程界关注。

目前国内外一些学者通过现场测试[1~3]、理论分析[4~8]、数值模拟[9~14]等方法对顶管施工造成的地表沉降进行了研究。由于顶管施工中的管土相互作用是一个复杂的动态三维过程，因此仅开展单纯的理论研究难以解决具体的工程问题，而数值模拟并结合现场的实测结果进行分析，更有助于理解顶管施工对周围土体扰动的机理，也可以为工程施工中遇到的问题提供合理的解决方案。实际的顶管工程中往往穿越的是错综复杂的地层环境，而目前的顶管工程数值模拟往往在简化后的均匀地层模型中进行，与实际的工程地质环境不相符。参照盾构隧道模拟研究成果，采用符合实际的复杂地层模型开展数值模拟，更加符合工程实际[15~17]。

所谓复杂地层，不仅仅指顶管工程横断面上的地层分布的复杂性（成层分布、不均匀分布），更重要的是指顶管工程纵断面上地层的不均匀分布。现有的数值模拟对于复杂地层的处理一般均是横断面层状均匀分布，而纵断面则分段层状均匀分布，导致各分段边界呈台阶状，与工程实际的差距仍然较大。为此，结合某市北线引水工程隐患整改工程，本文提出真实复杂地层模型建模方法，用于模拟顶管施工对周围地层扰动及地面沉降的影响，作为对比，同时建立常规简化均匀地层模型并进行数值模拟，通过数值模拟结果与现场实测结果的比较，说明采用本文提出的真实复杂地层模型的合理性和必要性。相关研究结果可供工程实际及数值模拟时参考。

数值建模及顶管施工模拟方法

工程概况

某市北线引水工程隐患整改工程E段采用顶管施工，采用NPD3400泥水平衡顶管机顶进，顶管为DN3400×32mm钢管。顶管外径为3.4m，壁厚为0.032m，顶管机头外径为3.46m。选择的模拟区段位于工作井和接收井中间，且距工作井和接收井距离较远，可忽略工作井和接收井施工造成的土体扰动，因此本文中的地层变形及地表沉降量可认为均由顶管施工引起。

真实复杂地层三维数值建模方法

本文基于有限差分软件FLAC3D进行数值模拟。虽然FLAC3D6.0具备强大的运算能力和较好的后处理功能，但其前处理能力较弱，在建立复杂地层的三维数值模型时存在一定困难。因此本文提出采用有限元软件ABAQUS来建立复杂地层的三维数值模型。首先根据实际工程地质勘测确定的地层纵断面图，对数值模型进行空间划分，然后在ABAQUS中将模型剖分为四面体四节点网格单

元。将地层信息和剖分网格数据导入FLAC3D6.0，然后进行边界条件设置、物理力学参数赋值及地应力平衡。最终得到真实物理力学参数及地层分布的三维数值模型。

为了进行对比分析真实复杂地层的准确性，在模型纵向采用各土层体积等量代换的方法对土层进行简化处理，生成纵向地层均匀成层分布的简化地层模型，并导入FLAC3D6.0进行顶管顶进施工模拟。

由于顶管工程的钢管每节长度为6m，因此数值模型的纵向模拟范围取6m的整数倍，本文取

60m，数值建模时将开挖区和注浆层单元均分为10段，每段长6m。相应地，数值模型在横向模拟范围也取60m，顶管两侧模拟范围约为8倍的顶管直径，可以认为模型两侧边界条件对模拟结果没有影响，故在两侧边界上均匀施加法向（水平）约束。垂直方向，模型上表面取为地面，顶管轴线埋深约12m，取模型下部模拟范围为20m，超过顶管直径的5倍，底部边界的影响可忽略不计，故模型上表面为自由边界，下表面施加法向（垂直）和切向（水平两向）约束。此外，在模型轴向，也在相对两表面施加法向（水平）约束。

根据实际工程地质勘测资料，数值模拟段地层由上至下分别为素填土、粉质黏土、粉砂岩残积土、全风化风砂岩，各地层沿顶管纵向分布不均，尤其是素填土层和粉质黏土层层厚差距大。在模型纵向模拟范围内，素填土层厚2.2m~11.4m，粉质黏土层厚3.2m~16.8m，粉砂岩残积土层厚

m~10.1m，全风化风砂岩层厚7.4m~10.0m（部分）；采用各土层体积等量代换的方法对土

层进行简化处理后，在模型纵向模拟范围内各地层等效厚度分别为素填土层厚7.0m，粉质黏土层

厚8.0m，粉砂岩残积土层厚9.1m，全风化风砂岩层厚7.9m（部分）。

建立的真实复杂地层数值模型如图1(a)所示，相应简化均匀地层数值模型如图1