土压平衡式盾构施工三维有限元模拟方法研究.docxVIP

土压平衡式盾构施工三维有限元模拟方法研究 1 土体抗侧力的数值模拟 随着结构设计技术的完善，结构法广泛应用于地下铁道的施工中。由于地质条件和施工工艺的限制，盾构推进对周围土体的扰动影响仍是不可避免的。如何预测盾构穿越所引起的地层位移，确保已有建(构)筑物正常使用和盾构的顺利掘进，是盾构隧道设计与施工中非常关键的问题。 目前工程中用于预测隧道施工引起的土体变形的主要方法是建立在现场实测资料基础上的经验公式法，提供的信息大都限于地表面处的位移，很难满足实际需要，因而在应用上具有很大的局限性。有限元法能更充分地反映土体性质的影响，可以考虑地层损失及土与衬砌的相互作用。有限元法不仅能够计算地面沉降，还可得到盾构推进过程中的沉降变化情况，不同深度处的沉降、侧向位移，地层应力变化以及对临近建筑物的影响。有限元法在盾构隧道施工分析中得到越来越广泛的应用。 盾构隧道掘进尤其是靠近盾构所在位置及其前方土体的位移是一个三维问题，平面有限元分析无法模拟盾构推进的实际情况，无法模拟开挖面周围土体的三维移动，也无法反映在到达平面应变状态之前盾构周围土体所经历的应力路径，计算结果与实际应力、位移有较大差异。尽管人们都认识到这一点，但是盾构隧道施工对周围土体的影响因素很多，已有的盾构施工三维有限元模拟方法大都作了较多的简化：文开发了一种三维弹塑性有限元技术对浅埋隧道无衬砌和完全衬砌2种极限情况下隧道周围土体不排水情况下的应力场和位移场进行了计算，并根据三维有限元计算结果提出用间隙系数G来反映盾构施工引起的地层损失；文在模拟盾构推进时在开挖面前方设置开挖单元，并考虑了盾构自重、千斤顶推力、土压仓压力等因素，但没有考虑盾构刀盘超挖和盾尾脱空引起建筑空隙的闭合以及盾尾注浆的影响；文在对土压平衡式盾构推进的三维模拟时考虑了开挖面土压力、盾尾空隙闭合、盾尾注浆等因素，但没有考虑盾构自重、千斤顶推力的影响；文用刚度迁移法模拟盾构推进，采用权刚度修正Goodman单元处理存在两种材料的混合接触刚度，提出了分3步模拟盾构前行一个单元距离的模拟步骤。 本文在前人工作的基础上，提出一种盾构施工三维有限元模拟方法，综合考虑盾构推进、开挖面前方及隧道壁面土体的剪切扰动、盾构刀盘超挖空隙和盾尾空隙的闭合、盾尾注浆等因素。 2 土体受力状态变化的有限元分析 隧道盾构施工中盾构机渐进向前，周围土体则是相对静止的。有限元法难以做到完全模拟盾构的连续推进过程，必须作一定的简化。通常将盾构推进作为一个非连续的过程来研究，假设盾构一步一步跳跃式向前推进，每次向前推进的长度(纵向)恰好为一个衬砌单元宽度，用改变单元材料类型的方法(刚度迁移法)来反映盾构的向前，而一次向前的过程中盾构周围土体受力状态也发生变化，将这种变化转化为相应的结点荷载作用于结点，进行有限元计算。这与在模拟基坑开挖时每一级开挖掉一层网格以及模拟土石坝、路堤填筑时每一级荷载增加一层网格的处理方法是相似的。 本文所考虑的盾构施工过程中引起土层移动的主要因素有以下6个部分(如图1(a)所示)： (1) 单元应力改变的开挖面结论 土压平衡式盾构开挖面支护压力难以和土体原位土压力达到完全平衡，当支护压力小于原位土压力时，开挖面前方土体向土压舱位移；反之，前方土体远离土压舱位移(图1(b))。 由于盾构施工中开挖面存在支护力，因而不能像模拟基坑开挖那样用上一步得到的单元应力直接计算新开挖面处的开挖释放荷载。常用的方法是每向前开挖一步将开挖面处的土体原位土压力与土压舱支护压力的差值作用在开挖面结点上，实际上这种处理方法并不合理，因为它使得开挖面结点荷载沿纵向(沿推进方向)的总和与每一步推进的距离有关，即在模拟盾构前进相同距离的情况下，每一步推进距离越小推进步数越大，开挖面结点荷载沿纵向总和也随之增加，这显然与实际情况不符。 本文建议根据开挖面处单元应力改变量{?σ}计算相应的结点荷载，即 这样处理后，不论每一步推进距离多少，盾构前进相同距离时开挖面结点荷载沿隧道纵向的总和保持不变。由于开挖面处土体受到盾构刀盘的搅动，开挖面上土体状态可以近似认为由原来的静止土压力状态变为各向等压状态(等于土压舱压力)，由此来计算单元应力改变量。 (2) 硬脆材料的抗拉受力 为了减小盾构壳上的摩擦阻力以使盾构能够顺利前行，通常盾构的刀盘外径要大于盾构壳的外径，从而在盾构壳外围产生一定厚度的间隙δ，盾构机由于自重发生下沉到底部，在横断面上形成超挖间隙2δ，随后周围土体由于填充超挖间隙而产生指向盾构内部的径向位移(图1(c))。 (3) 盾构推进方向的作用力 盾构向前掘进时盾构壳与周围土体之间发生错位，周围土体对盾构产生摩擦力，同时盾构也对周围土体施加指向盾构推进方向的摩擦力，周围土体受到盾构壳的剪切和挤压作用向前方和远离盾构的方向移动(图1(d))。盾构推进时