基坑开挖位移及土体水平抗力比例系数的反分析

李佳

摘要：基坑设计土体参数的正确选取直接影响支护结构的位移和内力计算结果，特别是土体水平抗力比例系数m值对支护结构位移的影响最为显著。结合基坑开挖工程实例，通过材料性状与基坑开挖施工过程的平面有限差分法模拟，模拟计算基坑分布开挖过程结构的水平位移；通过分析诸多土体参数对基坑支护结构水平位移的影响程度，确定杨氏模量作为位移反分析所需获取的参数，构建实测位移与模拟计算位移之间的目标函数，反分析给出各层土的杨氏模量；基于弹性地基梁法推导的m值与杨氏模量的关系，给出反映实际侧向变形特性各层土的m值。研究结果为基坑工程设计参数的合理确定提供了一定理论参考依据，对有效控制基坑位移与变形具有重要的实践意义。

关键词：基坑；平面有限差分法；水平抗力比例系数m值；杨氏模量；弹性地基梁；位移反分析

中图分类号：TU470 文献标志码：A 文章编号：1672-1683（2017）02-0185-07

在基坑工程中，一般都是已知边界条件、力学参数来求解结构的内力与位移。然而，许多基坑工程中材料参数往往是未知的，通过试验测定费时费力。为此，工程界提出了反分析方法。反分析是通过现场获取的少量易测数据（如位移），以这些数据为已知量，通过理论方法反求不易量测的材料参数，获得所需合适的参数用于结构设计或变形控制计算。国家基坑规范和许多地方规范已将弹性抗力法即m法作为基坑支护结构设计的推荐方法土体水平抗力比例系数m值成为关键参数。m值受诸多因素的影响，如基坑开挖深度、基坑内土体加固、开挖卸载、施工扰动等。规范对m值仅给出了较大的经验取值范围，选取具有较大的随机性和不确定性。位移反分析为m值的准确确定提供了一种重要的途径。熊巨华川用试验结果和实测数据证明采用位移反分析方法是确定m值的有效方法。徐中华提出了一种结合U code反分析软件和Abaqus有限元分析软件，根据基坑围护墙实测变形反分析m值的方法。张军平通过理论计算和实测数据的对比，认为对于浅基坑，采用基于弹性支点法的m值的位移反分析方法可达到工程需要参数的反演精度。陆勇根据基坑工程施工各工况实测位移，采用单纯形法位移反分析，对m值进行分析，反演得到优化位移与实测位移基本接近，建议m值取下限来控制基坑的最大变形。王旭东根据基坑施工监测资料提出利用非线性单纯形方法反演地基土的m值。戴加东、宋建学等根據基坑工程变形实测值采用相关软件或理论方法对m值进行了反演分析。

本文依托深港西部通道深圳侧接线工程道路部分城区地段基坑工程采用平面应变有限差分法模拟基坑分布开挖施工过程，通过构建实测位移和计算模拟位移之间的目标函数，采用位移反分析方法分析基坑土层变形参数，利用基于弹性地基梁法推导的杨氏模量与m值的关系确定m值。为合理确定基坑支护结构设计参数m值及基坑工程变形控制提供了一定的理论依据。

1工程简介

1.1工程概况

深港西部通道深圳侧接线工程是西部通道过境车辆专用通道，全长约4.48 km，道路部分区段采用下穿地道形式。城区地段道路两侧建（构）筑物密集，沿线埋设有市政、电力和通讯管线，周边环境对基坑结构的变形、位移和地下水的控制有很高的要求。该段线路基坑开挖深度范围内有一层5.0～8.0 m厚的砾砂，具有强透水性，且地下水极为丰富。基坑支护工程的重点是地下水治理、基坑支护结构变形控制以及周边建（构）筑物的保护问题。

城区地段周边复杂路段对应里程K2+100-K3+460，线路长为1.36 km。地道工程基坑开挖深度10.0～12.0 m，采用钻孔咬合排桩加多排钢支撑的支护方案，安全等级为一级。钻孔咬合排桩采用全长跟套管钻进法施工，具有较好的挡土、截水功能，避免因基坑开挖而造成地下水大量流失，发生流砂等事故；钢支撑具有较大的支撑刚度，通过施加预应力，可以将基坑的位移和变形控制在一定范围内。

1.2工程地质与水文地质条件

该工程城区地段所经场地原始地貌为海冲积平原，根据地质钻探揭露，地层自上而下为4层。（1）人工填土层。道路填土层，经碾压后呈中密-密实状态，填筑时间超过10年。（2）粗砾砂层。厚度5.0～8.0 m，为较纯净的粗砾层，呈中密状态，透水性好。（3）冲洪积层。主要为黏土层，局部有淤泥质黏土分布。（4）残积土层。混合花岗岩风化残积而成，主要为砂质黏性土，压缩性小、强度高、透水性小，为良好的持力层。城区地段地下水埋深较浅，埋深为1.0～2.0 m，粗砾砂层透水性大，地下水由海水补给，地下水丰富。

2基坑开挖的平面有限差分模拟方法

基坑开挖施工是一个动态的过程，为了能够近似模拟基坑的开挖响应，假设基坑施工每一个步骤都近似达到平衡，即节点的不平衡力为零或者小于给定的某个阈值，将基坑开挖看作轴对称二维平面应变问题，采用平面有限差分法进行模拟。

2.1材料性状的模拟

目前，土工计算中被广泛采用各向同性的本构模型，以Duncan&Chang非线性弹性模型与MohrCoulomb弹塑性模型为典型代表。Mohr Coulomb弹塑性模型能够反映土体不可恢复的塑性应变，计算参数较易确定，能够反映土体的剪切破坏，在实际应用中比较普遍。Vermeer&deBorst的试验报告证明，砂土和混凝土的应力应变关系与MohrCoulomb模型较为一致。模拟中土体采用MohrCoulomb弹塑性本构模型。

钻孔灌注桩采用梁单元来模拟，单元参数根据混凝土的弹性模量和泊松比μ来选取。由于模拟分析采用平面应变问题，对于平面应力的结构单元应将弹性模量除以（1-μ²）。水平钢支撑也选用梁单元，单元参数也根据混凝土或钢材的弹性模量和泊松比来选取。

基坑开挖过程中，每步开挖都要卸掉一部分土层，处理的方法是直接从模型中将其去掉，设置为空单元。基坑开挖前或者开挖到一定深度后，需要对坑底软土进行注浆加固，直接将该加固土层设置为独立对象，加固土层材料参数主要依靠现场数据和经验选取。

土体与墙体的交接界面相对来说是光滑的，界面的强度低于邻近土体的强度，界面摩擦角取12.5°，黏聚力取为2.5 kPa。对于界面的法向刚度和切向刚度的选取，根据有关文献选取相邻土体最大刚度的10倍左右是合理的。这样选取的目的是让界面的变形对整个模拟的变形产生最小的影响，然而又不能取过大的值，这样会产生相邻节点刚度差过大，产生畸变。相邻土体最大刚度可用下式得到：

max[（K+（4/3）G）/△z_min] （1）式中：K为体积模量；G为剪切模量；△z_min为与界面单元邻近土体的最小法向距离。

2.2基坑开挖的模拟

自重应力的模拟。通常情况下，地基土是一个半无限体，任何垂直面都是对称面，在自重应力的作用下，不可能发生侧向剪切变形，地面以下任何深度z处的垂直应力σ_y0等于上覆土层的白重应力。对于不考虑超载作用的地基其初始应力场只考虑土自重应力的作用。水平应力q₀为侧压力系数K₀乘以垂直应力。K₀取0.5。

围护桩施工过程的模拟。灌注桩的置入过程，主要是通过在安置桩的位置将土体单元设置为空单元，然后置入梁单元。桩的自重会使土体的应力场和位移场变化，主要是通过界面单元与土体发生相互作用产生。

基坑降水的模拟。基坑降水有坑内降水和坑外降水，本次模拟的是坑内降水。坑内降水主要有两种模拟路径，一种是开挖一定深度前降水到该深度以下一定范围，然后开挖，开挖后又继续降水到下一个指定的开挖面以下位置，这样循环直到基坑开挖至基坑底部标高；另一种是在开挖前一次性降水到基坑底以下一定位置，然后开挖，在开挖过程中保持基坑内水位不变。根据实际情况，选用第二种降水路径模拟基坑降水，一次性降水至坑底1 m。

基坑开挖土体和支撑结构施工的模拟。在基坑降水后，开挖土体，先开挖到第一道水平支撑中心下0.5 m处，计算系统稳定后，架设第一道水平支撑，预加一定轴力，计算系统稳定后，继续开挖至第二道水平支撑下0.5 m处，计算系统稳定后，架设第二道水平支撑，预加一定轴力，计算系统稳定后，开挖到基坑底设计标高，计算系统稳定。模拟中的每一个步骤能够计算应力和位移的重分布，稳定后再模拟下一个步骤，能够模拟多支撑支护结构的内力和变形随开挖过程的变化情况。

2.3位移反分析模拟过程

土体材料参数泊松比μ对位移的影响不敏感变化通常很小，均在0.25～0.35之间变化。内摩擦角φ、和剪胀角φ也不敏感。对于松散的砂，一般剪胀角φ为刻度1°-2°密实的砂砾、粗（砾）砂等通常不考虑剪胀角的影响，φ为零。内摩擦角φ、泊松比μ和杨氏模量E对结构水平位移影响分析结果见图1-图3。

从图1-图3中水平位移影响因素分析可见，杨氏模量对水平位移影响最为显著，因此主要通过反分析得到杨氏模量E，而將内摩擦角φ、泊松比μ等参数作为已知值取白土工试验或现场试验。以基坑开挖监测得到的水平位移与有限差分数值模拟得到的位移值的误差函数为目标函数，即：

（2）式中：n为监测点总数；ui为第i测点的实测位移值：ui（X）为第i测点的计算位移值；X为所需要确定的参数即设计变量；D⁰为设计变量的区间，根据实际情况而定。

分步开挖模拟反分析按照如下过程进行：根据实测位移不断试算，按照每一步开挖平衡后计算位移值来确定目标函数。即从开始就反分析设计变量，把目标函数控制在允许范围内，不断丢掉最坏点，代之以既能使目标函数有所改进，又能满足约束条件的新点，不断试算，逐步调向“最优点”。当第一步开挖上层土体得到的参数所带来的误差结果在给定范围内时，认为上层土体的实际参数就在此计算参数附近，缩小搜索范围。继续第二步开挖，第三步开挖，一直到基坑开挖到底部。位移反分析模拟过程见图4。

3基坑开挖过程模拟反分析

3.1计算模型与计算参数

选取分析断面ZXK3+315，平面有限差分计算模型网格划分见图5，土层参数见表1。表1中杨氏模量根据标贯击数及两者的经验关系得到，注浆固结体的参数选择参考《工程地质手册》，由标贯击数得到，其余参数由土工试验得到。围护结构灌注桩与钢支撑计算参数分别见表2和表3。

3.2分析步骤

有限差分网格模拟基坑开挖的程序可分为如下步骤：（1）生成模型网格，设置材料参数与边界条件；（21调用程序确定初始应力，确定结构桩安装后的最初应力；（3）降低地下水位至基坑底1 m；（4）分布开挖第一步，开挖到第一道钢支撑中心线下0.5 m；（5）安装第一道钢支撑并预加轴力240 kN；（6）分布开挖第二步，开挖到第二道钢支撑中心线下0 5 m；（7）安装第二道钢支撑并预加轴力280 kN；（8）分布开挖第三步，开挖到基坑底面标高处。模型建立后对上述步骤分别计算并修正参数，得到每一个施工阶段支护结构的内力与水平位移等计算结果。

3.3模拟反分析结果分析

模拟过程从第一步开始，通过不断调整初始参数得到结构的水平位移，计算目标函数式直到达到允许的范围。开挖阶段目标函数中的位移实测值的选择见表4，反分析模拟开挖最终得到符合目标函数的土工参数杨氏模量见表5。

3.3.1水平基床系数

在水平荷载作用下，桩水平位移u越大，侧压力（即土的弹性抗力）σ越大，侧压力大小还取决于土的性质、桩身刚度及截面形状和入土深度，可用如下公式表示：

σ=Cu （3）式中：C为土的水平基床系数（kN/m³），可通过试验得到，如可以对桩在不同类别土质及不同深度实测水平位移u及侧压力σ后反算得到。大量试验表明，基床系数C不仅与土的类别及其性质有关而且也随深度而变化。关于基床系数分布规律种类繁多，其中被广泛采用的是假设基床系数C随深度成正比例增长即

C=mZ （4）式中：m为土体水平抗力比例系数。

3.3.2水平基床系数的反分析

如图6，在区域D内作用均布荷载p，某点o处的地面沉降量可用集中荷载作用下的地面弹性位移量计算公式：

（5）式中：r为计算点至集中荷载作用点的距离。以微面积dA上作用的荷载作为集中荷载Q，然后在区域D内积允即可求o点处的沉降量：

（6）

（7）式中：B为荷载面宽度；Iρ是各种形状的面积积分，称为位移影响值。

对于宽为B，长为L的长方形荷载，将式（6）在x方向从0到B，y方向从0到L积分，得到长方形荷载角点下的沉降量

（8）

（9）

（10）

（11）

（12）

（13）

（14）

（15）

（16）

不同土层的m值可由杨氏模量E来反映，硬土层E较大，m值也较大，故由此确定的m值可以考虑软土和硬土层的差异，符合实际工程。

3.3.3m值反分析结果

利用反分析得到的杨氏模量E和式（15），得到比例系数m值，计算结果见表6。

4结语

（1）结合基坑工程实例，假设基坑施工每一步骤近似达到平衡状态，将基坑分布开挖过程看作平面应变问题处理，建立了二维有限差分计算模型，利用平面有限差分法模拟了基坑分布开挖过程。

（21分析内摩擦角φ、泊松比μ和杨氏模量E等土体参数对基坑支护结构水平位移的影响程度确定杨氏模量E作为位移反分析获取的参数。

（3）构建了以基坑开挖实测位移与有限差分数值模拟计算位移之间的目标函数，结合基坑工程典型开挖断面，反分析模拟基坑分布开挖施工过程，获取符合目标函数的各土层的杨氏模量。

（4）基于弹性地基梁弹性理论推导了土体杨氏模量与土体水平抗力比例系数m值之间的关系利用反分析获取各土层的杨氏模量获得了各土层的m值的建议取值。

上述分析结果为基坑工程设计参数的合理确定提供了一种重要途径，对基坑结构位移与变形控制提供了一定的理论依据和实践意义。