基于围岩变形时空效应的顶管顶推力计算方法专利检索-时空数学与统计专利检索查询-专利查询网

积极推动地理标志专门立法

2022-03-10 地理标志，立法，知识产权
保护知识产权是对创新最大的激励

2022-03-10 保护知识产权，创新，激励
谢商华：加快制定知识产权基本法

2022-03-10 知识产权基本法
擦亮“双奥之城”品牌

2022-03-10 双奥，知识产权
让冰雪运动“热”力全开

2022-03-10 冰雪运动，知识产权
携手共奋进　走好强国路

2022-03-10 强国，知识产权
坚持创新引领　方能稳中求进

2022-03-10 创新，稳中求进，知识产权
答好“两张卷” 奋进新征程

2022-03-10 知识产权
专家解读政府工作报告中的创新和知识产权相关部署

2022-03-10 政府工作报告，创新，知识产权
今年政府工作报告指出：加强知识产权保护和运用

2022-03-10 政府工作报告，知识产权保护

基于围岩变形时空效应的顶管顶推力计算方法

阅读：778发布：2021-02-28

IPRDB可以提供基于围岩变形时空效应的顶管顶推力计算方法专利检索，专利查询，专利分析的服务。并且本发明涉及一种基于围岩变形时空效应的顶管顶推力计算方法，其是将侧向摩阻力再叠加上迎面阻力得到总的顶推力，作为基于围岩变形时空效应的顶管顶推力。可以考虑顶管施工过程中围岩变形的时空效应，对顶管不同位置的荷载进行分析，最后积分得到总的顶推力，所得结果更加符合实际情况，可得到更为精细的顶推力计算结果，有利于对顶管结构、始发井、中继间和后靠背进行更为精细的设计。，下面是基于围岩变形时空效应的顶管顶推力计算方法专利的具体信息内容。

权利要求

1.一种基于围岩变形时空效应的顶管顶推力计算方法，其特征在于：将侧向摩阻力再叠加上迎面阻力得到总的顶推力，作为基于围岩变形时空效应的顶管顶推力，侧向摩阻力包括顶管不同位置径向荷载产生的侧向摩阻力，顶管不同位置径向荷载产生的侧向摩阻力是根据施工过程中所采取的减阻措施，对管土之间摩擦系数进行确定后得到的，顶管不同位置的径向荷载是根据顶进速度得到顶管线路方向上不同位置的接触情况计算得到的，顶管不同位置的径向荷载的计算还根据管土相互作用关系进行；

其包括以下步骤：

S1、获取包括顶管材料、埋深、断面尺寸、地层物理力学参数在内的基础参数，基础参数包括：地层压缩模量Es、泊松比μ、内摩擦角φ、顶管材料与地层摩擦系数fs、顶管埋深h、顶管机外径D0、管节外径D1、单位长度管节自重ω；

S2、在以上参数基础上对围岩变形特性进行预测，根据顶进速度得到顶管线路方向上不同位置的接触情况，并根据管土相互作用关系对顶管不同位置的径向荷载进行计算；

S3、根据施工过程中所采取的减阻措施，对管土之间摩擦系数进行确定，并由此得到顶管不同位置径向荷载产生的侧向摩阻力；

S4、对顶管摩阻力沿顶进方向进行积分，得到总的侧向摩阻力，再叠加上迎面阻力就得到总的顶推力。

2.如权利要求1所述的基于围岩变形时空效应的顶管顶推力计算方法，其特征在于，步骤S2中，顶管不同位置的径向荷载的表征包括单位长度内总的径向荷载Pr，其计算公式为：当管节只在底部与土体接触，管节与土体之间的接触荷载为结构自重时，pr＝ω，式中，ω为单位长度管节自重；

当围岩变形发展到已经与管节接触时，pr＝f1(ES,μ,h,D0,D1,t)+ω，式中，表示围岩的形变荷载部分，是与围岩变形时空效应相关的函数，自变量ES,μ,h,D0,D1是步骤S1所获取的基础参数，t是围岩开挖后的时长。

3.如权利要求2所述的基于围岩变形时空效应的顶管顶推力计算方法，其特征在于，表征围岩变形程度的围岩变形量为：其中，P(t)为释放荷载与时间的关系，E为ES。

4.如权利要求2所述的基于围岩变形时空效应的顶管顶推力计算方法，其特征在于：管土作用荷载的计算公式为其中，P(t)为释放荷载与时间的关系，t0为管土刚接触的时刻。

5.如权利要求2所述的基于围岩变形时空效应的顶管顶推力计算方法，其特征在于，步骤S3中，顶管上某截面单位长度的侧摩阻力为：Ps(l)＝πD1Prfs。

说明书全文

基于围岩变形时空效应的顶管顶推力计算方法

技术领域

[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、顶管机外径D0、管节外径D1、单位长度管节自重ω

[0021] 本发明一个实施例中，步骤S2中，顶管不同位置的径向荷载的表征包括单位长度内总的径向荷载Pr，其计算公式为：

[0022] 当管节只在底部与土体接触，管节与土体之间的接触荷载主要为结构自重时，pr＝ω，式中，ω为单位长度管节自重；

[0023] 当围岩变形发展到已经与管节接触时，pr＝f1(E,μ,h,D0,D1,t)+ω，式中，表示围岩的形变荷载部分，是与围岩变形时空效应相关的函数，自变量E,μ,h,D0,D1是步骤S1所获取的基础参数，t是围岩开挖后的时长。

[0024] 本发明一个实施例中，表征围岩变形程度的围岩变形量为：

[0025] 其中，P(t)为释放荷载与时间的关系。

[0026] 本发明一个实施例中，管土作用荷载的计算公式为

[0027] 其中，P(t)为释放荷载与时间的关系，t0为管土刚接触的时刻。

[0028] 本发明一个实施例中，释放荷载与时间的关系P(t)为：式中， V为工作面平均推进速度。

[0029] 本发明一个实施例中，当时，判定土体与顶管没有完全接触；当时，判断土体与顶管已经接触。

[0030] 优选为，步骤S2中，先计算得到所对应的时刻t0，再根据顶管顶进速度V(即工作面平均推进速度)，计算得到：管土接触分界长度l0＝Vt0。

[0031] 本发明一个实施例中，步骤S3中，顶管上某截面单位长度的侧摩阻力为：Ps(l)＝πD1Prfs。

[0032] 优选为总的侧摩阻力为：

[0033] 本发明一个实施例中，步骤S4中，总的顶推力P的计算公式为式中：PF为顶管迎面阻力，为总的侧摩阻力。

[0034] 本发明的有益效果是，可以考虑顶管施工过程中围岩变形的时空效应，对顶管不同位置的荷载进行分析，最后积分得到总的顶推力，所得结果更加符合实际情况，可得到更为精细的顶推力计算结果，有利于对顶管结构、始发井、中继间和后靠背进行更为精细的设计。

附图说明

[0035] 图1是本发明的分析流程图；

[0036] 图2是本发明第一实施例的纵剖面构造图；

[0037] 图3是图2中的截面示意图(其中，A为A-A截面，B为B-B截面，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、顶管机外径D0、管节外径D1、单位长度管节自重ω。

[0052] 步骤S2：径向荷载计算。

[0053] 通常情况下顶管机外径D0会大于管节外径D1(D0>D1)，从而导致，在地层变形之前，管节(即顶管)上部与地层无接触，随着地层变形的逐渐发展，土体慢慢与管节接触，该过程管节所受径向荷载逐渐增大。

[0054] 如图2、图3所示，对于A-A截面位置：管节只在底部与土体接触，管节与土体之间的接触荷载主要为结构自重，此时，单位长度内总的径向荷载为：pr＝ω。

[0055] 如图2、图3所示，B-B截面为围岩变形发展到一定程度，已经与管节接触，此时单位长度内总的径向荷载为：pr＝f1(E,μ,h,D0,D1,t)+ω。

[0056] 式中：表示围岩的形变荷载部分，是与围岩变形时空效应相关的函数。其中自变量E,μ,h,D0,D1就是步骤1所获取的参数，t是围岩开挖后的时长。

[0057] 随着围岩变形持续发展，会在某一时刻与管节接触，在接触前后顶管受到的荷载不同，因此首先必须计算管土接触的分界点，确定分界长度l0。

[0058] 围岩变形量为：

[0059] 上式中，释放荷载与时间的关系P(t)为：

[0060] 式中： V是工作面平均推进速度。

[0061] 当：

[0062] 时，土体与顶管没有完全接触，该情况为截面A-A情形；

[0063] 时，土体与顶管已经接触。

[0064] 计算得到所对应的时刻t0。

[0065] 根据顶管顶进速度V(即工作面平均推进速度)，进一步得到：管土接触分界长度l0＝Vt0。

[0066] 在得到分界长度l0和对应的时刻t0后，就可以按下式对管土作用荷载进行计算：

[0067]

[0068] 步骤S3：侧摩阻力计算。

[0069] 顶管上某截面单位长度的侧摩阻力为：Ps(l)＝πD1Prfs。

[0070] 总的侧摩阻力为

[0071] 步骤S4：总的顶推力计算。

[0072] 总的顶推力

[0073] 式中：PF为顶管迎面阻力；为总的侧摩阻力。

[0074] 参见图2和图3，本发明的顶管顶推力计算的一个实施例。

[0075] 假设某顶管工程初始参数如下：顶管机外径D0＝505cm，混凝土管的外径D1＝500cm，每米管的重力ω＝29.4kN/m，混凝土容重γ＝26kN/m3，土的天然容重γ＝20kN/m3，内摩擦角 Es＝10MPa，μ＝0.3，管与土的粘着力c＝10kPa，标准贯入数N＝4，覆土深度H＝10m，无地下水，顶程L＝100m，日推进距离3m。

[0076] 1)按照现有技术，顶管规范计算过程与结果如下：

[0077] 按《给水排水工程顶管技术规程》和《顶管工程施工规程》中采用的公式P＝πD1Lfs+PF

[0078] 根据规范和实施例的上述参数可得到，上式中：

[0079] fs＝5kN/m2

[0080]

[0081] Pl＝πD1lfs+PF＝78.5l+3926(kN)

[0082] 不同顶进长度的顶推力分布如图4所示。

[0083] 2)本发明计算过程与结果

[0084]

[0085] α根据不同地层取不同值：

[0086]

[0087]

[0088] 当t＝t0：

[0089] 得：t0＝2.5天

[0090] 进一步得到l0＝7.5m

[0091] 因此在最前面7.5m，顶管竖向荷载由顶管结构自重组成，在7.5m之后，竖向荷载为：

[0092]

[0093] 最后得到总的顶推力随顶推长度关系为：

[0094]

[0095] 根据前面分析结果，代入各计算参数，利用数值积分方法，可得最终顶推力如图5所示。

[0096] 结合图4、图5对比可知，本发明计算结果与现有规范计算结果有一定差距，且本发明计算所得结果与实际中顶推力的非线性变化趋势较为相符，而且，本发明计算结果所需的顶推力小于现有技术计算得到的数值，即采用本发明的方法可以得到较小的合适的顶推力。

[0097] 综上所述，本发明的顶管顶推力计算方法，与现有方法相比最大的优势在于：能够考虑实际工程中围岩变形的时空效应，从而获得更为精确的顶推力计算结果。

标题	发布/更新时间	阅读量
时空棋-专利编号CN1076132A	2020-05-11	107
时空表盘-专利编号CN104597740A	2020-05-11	493
一种基于时空场景的信息推送系统和方法-专利编号CN104166929A	2020-05-13	516
一种时空仪-专利编号CN1834824A	2020-05-13	319
时空调节子-专利编号CN110088285A	2020-05-13	830
一种新型时空色彩学习平台系统-专利编号CN109938486A	2020-05-11	128
时空钟表-专利编号CN2884270Y	2020-05-12	753
时空计-专利编号CN2256110Y	2020-05-12	1049
时空钟-专利编号CN2273430Y	2020-05-11	158
时空仪-专利编号CN202443688U	2020-05-12	421

基于围岩变形时空效应的顶管顶推力计算方法

基于围岩变形时空效应的顶管顶推力计算方法

技术领域

背景技术

发明内容

附图说明

具体实施方式

IPRDB

热门服务

关于我们

友情链接

联系方式