基于围岩变形时空效应的顶管顶推力计算方法
【技术领域】
[0001] 本发明涉及一种隧道与地下工程建设中的顶管顶推力计算方法，尤其是基于围岩 变形时空效应的顶管顶推力计算方法。
【背景技术】
[0002] 在隧道和地下管线建设过程中，顶管施工作为一种非开挖方法，可以在不开挖地 表土的情况下完成工程施工，具有无可比拟的优点，得到了越来越广泛的运用。顶管施工过 程中，顶推力计算一直以来都倍受关注，其直接关系到顶管结构、始发井、中继间和后靠背 的设计。
[0003] 当前针对顶推力的计算主要有两种模型：开挖面稳定模型和管土全接触模型。开 挖面稳定模型认为:顶管机外径与管道外径差值所产生的缝隙一直存在，管道顶进过程中 围岩上部是稳定的，管道只在底部一定宽度的表面上滑动，在其它部分由于挖掘面保持稳 定，管土之间没有接触，因此，顶力主要由管道自重产生的摩擦阻力组成；管土全接触模型 认为:管道顶进过程中管道周围均与土体接触，顶力主要用来克服由作用在管周的土压力 引起的摩擦阻力。
[0004] 现有顶管顶推力计算公式很多，但几乎都是在以上两种模型的基础上发展而来 的，总的顶推力为迎面阻力与管节摩擦力之和，主要区别在于摩擦力的计算。当前顶推力计 算方法都视顶管沿线路方向的径向荷载为与施工过程和时间无关的不变量，即管土之间的 接触状态从始至终是不变的。
[0005] 这样简化虽然使得计算简单，易于理解和实际操作，但随着技术的发展，岩土行业 精细化设计要求越来越高的今天，亟需一种更加切合实际的顶管顶推力计算方法。本发明 涉及内容正是针对以上技术难点，可得到较为精细的顶推力计算结果。
【发明内容】
[0006] 为了克服当前顶管施工过程中顶推力计算没有充分考虑围岩变形时空效应的不 足，本发明提供了一种基于围岩变形时空效应的顶管顶推力计算方法。该方法能够考虑顶 管施工过程中沿线路方向围岩变形的不同特征，从而使计算结果更加切合实际，能有效指 导顶管结构、工作井、中继间以及后靠背等的设计。
[0007] 本发明解决其技术问题所采用的主要技术方案包括：
[0008] 一种基于围岩变形时空效应的顶管顶推力计算方法，将侧向摩阻力再叠加上迎面 阻力得到总的顶推力，作为基于围岩变形时空效应的顶管顶推力，该方法能够考虑顶管施 工过程中沿线路方向围岩变形的不同特征，从而使计算结果更加切合实际。
[0009] 优选为侧向摩阻力包括顶管不同位置径向荷载产生的侧向摩阻力。
[0010] 更优选为顶管不同位置径向荷载产生的侧向摩阻力是根据施工过程中所采取的 减阻措施，对管土之间摩擦系数进行确定后得到的。
[0011]更优选为顶管不同位置的径向荷载是根据顶进速度得到顶管线路方向上不同位 置的接触情况计算得到的更优选为顶管不同位置的径向荷载的计算还根据管土相互作用 关系进行。
[0012] 本发明一个实施例中，顶管不同位置的径向荷载的计算包括首先在基础参数的基 础上对围岩变形特性进行预测。
[0013] 优选为基础参数包括顶管材料、埋深、断面尺寸、地层物理力学参数。
[0014] 本发明一个实施例中，侧向摩阻力指的是，对顶管摩阻力沿顶进方向进行积分，得 到的总的侧向摩阻力。
[0015] 本发明一个实施例中，其包括以下步骤：
[0016] S1、根据工程基本情况获取包括顶管材料、埋深、断面尺寸、地层物理力学参数在 内的基础参数；
[0017] S2、在以上参数基础上对围岩变形特性进行预测，根据顶进速度得到顶管线路方 向上不同位置的接触情况，并根据管土相互作用关系对顶管不同位置的径向荷载进行计 算；
[0018] S3、根据施工过程中所采取的减阻措施，对管土之间摩擦系数进行确定，并由此得 到顶管不同位置径向荷载产生的侧向摩阻力；
[0019] S4、对顶管摩阻力沿顶进方向进行积分，得到总的侧向摩阻力，再叠加上迎面阻力 和管土粘滞力就得到总的顶推力。
[0020] 步骤S1中，基础参数包括:地层压缩模量Es、泊松比μ、内摩擦角Φ、顶管材料与地 层摩擦系数fs、顶管埋深h、顶管机外径Do、管节外径0:、单位长度管节自重ω
[0021] 本发明一个实施例中，步骤S2中，顶管不同位置的径向荷载的表征包括单位长度 内总的径向荷载Pr，其计算公式为：
[0022] 当管节只在底部与土体接触，管节与土体之间的接触荷载主要为结构自重时，Pr =ω，式中，ω为单位长度管节自重；
[0023]当围岩变形发展到已经与管节接触时，口1^ =负^，以，11，0()，〇1，1:)+〇，式中， 系我表示围岩的形变荷载部分，是与围岩变形时空效应相关的函数，自变 量Ε，μ，h，Do，Di是步骤S1所获取的基础参数，t是围岩开挖后的时长。
[0024] 本发明一个实施例中，表征围岩变形程度的围岩变形量为：
[0025]
，其中，P(t)为释放荷载与时间的关系。
[0026] 本发明一个实施例中，管土作用荷载的计算公式为
[0027]
，其中，P(t)为释放荷载与时间的关系，to为管土刚接 触的时刻。
[0028] 本发明一个实施例中，释放荷载与时间的关系P(t)为
艾 中
V为工作面平均推进速度。
[0029] 本发明一个实施例中，当％ <@4时，判定土体与顶管没有完全接触；当 Σ
，判断土体与顶管已经接触。
[0030] 优选为，步骤S2中，先计算得到 斤对应的时刻to,再根据顶管顶进速度 V(即工作面平均推进速度），计算得到:管工按删方开长度lo = Vto。
[0031] 本发明一个实施例中，步骤S3中，顶管上某截面单位长度的侧摩阻力为:Ps(l)=3i DlPrfso
[0032] 优选为总的侧摩阻力为
[0033] 本发明一个实施例中，步骤S4中，总的顶推力P的计算公式）
，式
中:Pf为顶管迎面阻丈 总的侧摩阻力。
[0034]本发明的有义术疋，η」以考虑顶管施工过程中围岩变形的时空效应，对顶管不 同位置的荷载进行分析，最后积分得到总的顶推力，所得结果更加符合实际情况，可得到更 为精细的顶推力计算结果，有利于对顶管结构、始发井、中继间和后靠背进行更为精细的设 计。
【附图说明】
[0035]图1是本发明的分析流程图；
[0036]图2是本发明第一实施例的纵剖面构造图；
[0037]图3是图2中的截面示意图（其中，Α为Α-Α截面，Β为Β-Β截面，C为C-C截面）；
[0038] 图4是依照现有技术的方法得到的本发明第一实施例的不同顶进长度的顶推力分 布示意图；
[0039] 图5是本发明第一实施例的方法得到的不同顶进长度的顶推力分布示意图。
[0040] 【附图标记说明】
[0041 ] 1、地表;2、工作井；3、后背墙;4、顶推系统；5、顶推块;6、管节;7、中继间；8、顶管 机。
【具体实施方式】
[0042]为了更好地解释本发明，以便于理解，下面结合附图通过【具体实施方式】对发明作 进一步描述。
[0043] 参见图1，本发明的分析流程图，本发明的基于围岩变形时空效应的顶管顶推力计 算方法，包括如下步骤：
[0044] S1、对于某一需要分析的特定工程，根据工程基本情况获取顶管材料、埋深、断面 尺寸、地层物理力学参数等基础信息；
[0045] S2、在以上参数基础上对围岩变形特性进行预测，根据顶进速度得到顶管线路方 向上不同位置的接触情况，并根据管土相互作用关系（管土相互作用关系即围岩与顶管之 间的相互作用关系)对顶管不同位置的径向荷载进行计算；
[0046] S3、根据施工过程中所采取的减阻措施，对管土之间摩擦系数进行确定，并由此得 到顶管不同位置径向荷载产生的侧向摩阻力；
[0047] S4、对顶管摩阻力沿顶进方向进行积分，得到总的侧向摩阻力，再叠加上迎面阻力 就得到总的顶推力。
[0048] 实际实施过程中，根据本发明计算流程图的步骤和顺序，具体如下：
[0049] 步骤S1:初始参数获取。
[0050] 根据实际工程获取如下初始参数：
[0051 ]地层压缩模量Es、泊松比μ、内摩擦角Φ、顶管材料与地层摩擦系数fs、顶管埋深h、 顶管机外径Do、管节外径以、单位长度管节自重ω。
[0052]步骤S2:径向荷载计算。
[0053] 通常情况下顶管机外径Do会大于管节外gDKDADi)，从而导致，在地层变形之前， 管节（即顶管)上部与地层无接触，随着地层变形的逐渐发展，土体慢慢与管节接触，该过程 管节所受径向荷载逐渐增大。
[0054