CN105825027A

CN105825027A - 顶管顶进过程的多体系统动力学数值模拟方法

Info

Publication number: CN105825027A
Application number: CN201610192677.9A
Authority: CN
Inventors: 吉小明; 龙兵; 许江海
Original assignee: Guangdong University of Technology
Current assignee: Guangdong University of Technology
Priority date: 2016-03-30
Filing date: 2016-03-30
Publication date: 2016-08-03

Abstract

本发明公开了一种顶管顶进过程的多体系统动力学数值模拟方法，隧道工程管幕预加固的研究领域。本发明方法的操作步骤为：第一步确定直线、曲线顶管顶进过程中顶管上的荷载分布和计算；第二步建立直、曲线钢顶管多体系统动力学模型；第三步确定模型计算物理力学参数；第四步计算模型和连接器单元运动约束属性的确定；第五步进行直线顶管多体模型分析；第六步进行直线顶管转为曲线顶管多体模型分析；第七步曲线顶管多体系统动学模型顶管顶进力预计计算。

Description

顶管顶进过程的多体系统动力学数值模拟方法

技术领域

本发明涉及隧道工程管幕预加固的研究领域，特别涉及一种顶管顶进过程的多体系统动力学数值模拟方法。

背景技术

顶管法是一种非开挖技术，主要用于地下管道的铺设，顶进管道时借助主顶油缸的顶推力,在长度较大时也可借助中继间来传递顶推力，把顶管掘机从始发工作井顶到接收井内并吊起，在此过程中顶管掘进后的管道直接铺设在始发井与接受井之间的土层中。1896年美国北太平洋铁路铺设工程中第一次使用顶管施工技术，并取得满意的效果。顶管施工技术在我国最早开始使用是在1953年的北京，经过多年的发展，顶管技术在我国已得到大量地实际工程应用，且保持着高速的增长势头，无论在技术上、顶管设备还是施工工艺上取得了很大的进步，在某些方面甚至已达到了世界领先水平。

发明内容

本发明的主要目的在于克服现有技术的缺点与不足，提供一种顶管顶进过程的多体系统动力学数值模拟方法，通过该方法探讨顶管顶进时管节之间的力学行为，提出了基于多体模型预测顶力的新方法，对顶进力与摩擦力和迎面阻力的关系进行讨论。

为了达到上述目的，本发明采用以下技术方案：

本发明提供了一种顶管顶进过程的多体系统动力学数值模拟方法，包括下述步骤：

S1、确定直线、曲线顶管顶进过程中顶管上的荷载分布并计算荷载；

S2、建立直、曲线钢顶管多体系统动力学模型，介绍多体动力学基本理论，了解动力学普遍方程、广义速度公式、动力学控制方程和多体系统的动力学控制方程；

S3、确定计算模型和连接器单元运动约束属性，操作Abaqus使用两节点连接单元在系统各部分之间建立联系，并通过定义连接属性来描述各部分之间的相对运动约束关系；

S4、确定模型计算物理力学参数；假定地层荷载、摩擦阻力、管端正面阻力都是均匀分布的，假定钢材为理想弹塑性钢材，管节模型取两节，通过连续单元连接两管节；

S5、进行直线顶管多体模型分析，钢管采用壳单元，选用S4R类型进行网格划分；

S6、进行直线顶管转变为曲线顶管多体模型分析，模拟中采用多节管道，曲线管道钢管通过接头连接而成；

S7、曲线顶管多体系统动学模型顶管顶进力预计计算，钢管采用壳单元，选用S4R类型单元进行网格划分。

作为优选的技术方案，步骤S1中，确定直线、曲线顶管顶进过程中顶管上的荷载分布并计算荷载具体为：

确定顶管上的土压力，分别采用土柱理论、普氏理论和太沙基理论进行分析，与实验对比确定采用太沙基理论；

确定顶进过程中的单位面积摩擦阻力，摩擦阻力受土的性质、管道的材料、润滑浆液的配合比、管-润滑浆液-土三者的接触状态的影响；

确定顶进时的顶力，顶力由管周摩擦阻力和管端正面阻力的大小决定。

作为优选的技术方案，步骤S2中，动力学普遍方程为：

\frac{d}{d t} (\frac{\partial T}{\partial \overset{\cdot}{q}}) - \frac{\partial T}{\partial q} = Q;

其中，q为广义坐标向量；T为系统动能；Q为广义力；

系统动能T由广义速度表示为：

T = \frac{1}{2} {\overset{\cdot}{q}}^{T} M \overset{\cdot}{q};

其中，q为广义速度向量；M为质量矩阵；

物体的动力学控制方程为：

M^{i} {\overset{. .}{q}}^{i} + K^{i} q^{i} = Q_{F}^{i} + Q_{V}^{i};

其中，为广义加速度向量；M为质量矩阵；K为刚度矩阵；Q_F主动力对应的广义力；Q_V为速度的二次项有关的广义力；

多体系统的动力学控制方程为：

M \overset{\cdot\cdot}{q} + K q + C_{q}^{T} = Q_{F} + Q_{V}

其相应的约束方程为：

C(q,t)＝0

式中，λ为拉格朗日乘子列阵；C为约束矩阵；为约束的雅克比转置矩阵。

作为优选的技术方案，步骤S3具体为：

S31、在PART、ASSEMBLY或INTERACTION功能模块中，定义连接单元和约束所要用到的参考点和基准坐标系；

S32、在INTERACTION模块中，设置连接单元、连接属性和约束；

S33、在STEP模块中，设置单元的历史变量输出；如果模型中出现较大的位移或转动，应将几何非线性参数NLGEOM设置为ON；

S34、在LOAD模块中，定义边界条件和载荷，以及连接单元的边界条件和载荷；

S35、在VISUALIZATION模块中，查看连接单元的历史变量输出、控制连接单元的显示方式。

作为优选的技术方案，步骤S5中，直线顶管采用Slide-Plane和旋转连接属性Cardan；计算出每增加一个管节推力值增加为一个定值，主要原因是摩擦力的影响；在直线顶管顶进过程中管节与管节之间会出现一定的相对转动，且相对转角处于动态变化的过程；管节之间连接时允许一定的转动自由度，可以减少管节之间的转动弯矩。

作为优选的技术方案，步骤S6中，还包括确定连接属性的步骤，两管节之间采用四个连接单元对称分布在管周，曲线内侧连接单元为平移连接属性Join和旋转连接属性Cardan,其余连接单元为平移连接属性Slot和旋转连接属性Cardan；顶管直线转变为曲线的过程中，随着顶进距离的增加，管节之间的夹角在增大，且呈现出线性关系。

作为优选的技术方案，步骤S7中，连接属性采用平移连接属性Join和旋转连接属性Revolute；由于摩擦力的影响，每增加一节管节，推力值增加为一个定值，且与直线顶管增加值相近；管节之间的相对转角在顶进时表现出动态变化，但相对转动角都非常小，比直线顶管小很多；管道在曲线顶进过程中，为保持管节轨迹为曲线，管节之间会产生水平力，且呈现出水平力两头大中间小的规律；在长距离顶管中，单位面积摩擦阻力是影响顶力大小的主要因素。

本发明与现有技术相比，具有如下优点和有益效果：

1、本发明可计算出顶管顶进的总顶力，为顶进设备的选购提供依据。

2、本发明可计算出各管节之间的顶进力和相对角度，为管节的顶进提供可靠度，避免管节顶进过程中出现偏差而影响施工。

附图说明

图1为本发明顶管施工过程中受力模型图；

图2(a)、图2(b)分别为本发明直线顶管进受力、管节断面受力图；

图3(a)、图3(b)分别为本发明曲线顶管进受力、管节断面受力图；

图4为本发明直线钢顶管多体动力学模型；

图5为本发明曲线钢顶管多体动力学模型；

图6为本发明连接器单元位置示意图。

具体实施方式

下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。

实施例

本实例选取港珠澳大桥珠海段拱北暗挖隧道，全长255m，全线穿越珠海进入澳门的通关口岸拱北口岸和茂盛围军事管理区等重要区域，沿线地面建筑多，地下管线及邻近桩基础密集以及地质情况复杂。拱北暗挖隧道管幕工程由36根外径1620mm钢管组成，每根管道有64节管节，每节钢管长4m，管节之间采用F型承插口连接。管道的壁厚以整个管幕的中板为界，中板上部管道的壁厚为20mm，共17根，下部管道的壁厚为24mm，共19根。本隧道管道采用顶管法，顶管顶进过程的多体系统动力学数值模拟方法的操作步骤如下：

第一步，以港珠澳大桥拱北暗挖隧道29#管为研究对象，确定直线，曲线顶管顶进过程中顶管上的荷载分布和计算；确定管道上的土压力、顶进过程中的单位面积摩擦阻力和顶进时的顶力。可查看附图1的顶管施工过程中受力模型图、附图2(a)、图2(b)的直线顶管荷载示意图和附图3(a)、图3(b)的曲线顶管荷载示意图。

表1土体材料物理力学参数表

1)确定管端正面阻力：

P_f＝P_w+P_s+P

式中P_f为管端正面阻力，Pa；

P_w为地下水压力，Pa；

P_s为土体压力，Pa；

P为常数，一般取20～30kPa。

2)确定管道上的垂直土压力：

其中，N_z为垂直土压力，kPa；

k₀为朗肯主动土压力系数，

c为土体粘聚力，kPa；

b₁为土体破坏线宽度的一半，m；

γ为管道上覆土重度，kN/m³；

为土体内摩擦角，°；

H为拱高,m；

B为拱跨的一半,m。

3)确定管周法向土压力

管道上部法向土压力取管道顶部的垂直土压力N_z，管道底部的法向土压力是在管道顶部法向土压力的基础上加上管道自重为：

N_z2＝N_z1+G，其中，G为管道自重，kg；

管道左侧法向土压力为：N_x1＝k₀N_z1，其中，k₀为水平侧压力系数；

管道右侧法向土压力为：N_x2＝k₀N_z1，其中，k₀为水平侧压力系数。

在多体模型分析中为简化计算模型，可将管道顶部、底部、两侧法向土压力的平均值作为管周法向土压力加载到多体模型中，管周法向土压力可取为：

N＝(N_z1+N_z2+N_x1+N_x2)/4

4)确定顶进过程中单位面积摩擦阻力

F＝f·N·A＝M·A；

M = \frac{F}{A};

其中，F为单位长度摩擦阻力，kN；

f为土与管壁摩擦系数；

N为单位面积管道上的法向土压力，kN；

A为管道外表面的面积，m²；

M为单位面积摩擦阻力，kN/m²，M＝f·N。

5)确定顶进时的顶力

P＝F·L+P_f

其中，P为顶力，kN；

F为单位面积摩擦阻力，kN/m；

L为顶进距离，m；

P_f为管端正面阻力，KN。

以上公式，对于没有施工现场监测的情况下比较适用，有现场监测数据时，可以根据现场监测数据确定各种荷载。

第二步，建立直、曲线钢顶管多体系统动力学模型(如图4、图5所示)，介绍多体动力学基本理论，了解动力学普遍方程、广义速度公式、动力学控制方程和多体系统的动力学控制方程。

动力学普遍方程为：

\frac{d}{d t} (\frac{\partial T}{\partial \overset{\cdot}{q}}) - \frac{\partial T}{\partial q} = Q

其中，q为广义坐标向量；

T为系统动能；

Q为广义力。

系统动能T由广义速度表示为：

T = \frac{1}{2} {\overset{\cdot}{q}}^{T} M \overset{\cdot}{q}

其中，q为广义速度向量；

M为质量矩阵。

物体的动力学控制方程为：

M^{i} {\overset{\cdot\cdot}{q}}^{i} + K^{i} q^{i} = Q_{F}^{i} + Q_{V}^{i}

其中，为广义加速度向量；

M为质量矩阵；

K为刚度矩阵；

Q_F主动力对应的广义力；

Q_V为速度的二次项有关的广义力。

多体系统的动力学控制方程为：

M \overset{\cdot\cdot}{q} + K q + C_{q}^{T} = Q_{F} + Q_{V}

其相应的约束方程为：

C(q,t)＝0

式中，λ为拉格朗日乘子列阵；

C为约束矩阵；

为约束的雅克比转置矩阵。

第三步，确定计算模型和连接器单元运动约束属性，操作Abaqus使用两节点连接单元在系统各部分之间建立联系，并通过定义连接属性来描述各部分之间的相对运动约束关系。其基本分析步骤为：

(1)在PART、ASSEMBLY或INTERACTION功能模块中，定义连接单元和约束所要用到的参考点和基准坐标系；

(2)在INTERACTION模块中，设置连接单元、连接属性和约束；

(3)在STEP模块中，设置单元的历史变量输出；如果模型中出现较大的位移或转动，应将几何非线性参数NLGEOM设置为ON；

(4)在LOAD模块中，定义边界条件和载荷，以及连接单元的边界条件和载荷；

(5)在VISUALIZATION模块中，查看连接单元的历史变量输出、控制连接单元的显示方式。

第四步，确定模型计算物理力学参数；首先假定地层荷载、摩擦阻力、管端正面阻力都是均匀分布的；假定钢材为理想弹塑性钢材，管节模型取两节，通过连续单元连接两管节。模型主要取两节管节，29#钢管采用Q235,管节长度为4m，管节直径为1.62m，管壁厚度为20mm，通过连接单元连接两管节，钢管的弹性模量E＝210GPa，管道的密度为ρ＝7850kg/m³，泊松比ν＝0.3，屈服强度为3.45×10⁸Pa。

第五步，进行直线顶管多体模型分析，钢管模拟中采用壳单元，选用S4R类型进行网格划分，六节管节共划分12792个单元。直线顶管施工顶进中考虑限制两管节之间管轴方向的位移，因此本次直线顶管多体模拟中连接属性采用平移连接属性Slide-Plane和旋转连接属性Cardan，该组合连接属性中仅限制管轴方向的位移，可查阅附图4直线顶管连接属性图；计算出每增加一个管节推力值增加为34586N，主要原因是摩擦力的影响；在直线顶管顶进过程中管节与管节之间会出现一定的相对转动，且相对转角处于动态变化的过程，最大转动角出现在管节3与管节4处，最大值为4.27×10^-7度；管节之间连接时允许一定的转动自由度，可以减少管节之间的转动弯矩。

第六步，进行直线顶管转变为曲线顶管多体模型分析，模拟中采用6节管道，曲线管道钢管通过F接头连接而成；确定连接属性，两管节之间采用四个连接单元对称分布在管周，曲线内侧连接单元为平移连接属性Join和旋转连接属性Cardan,其余三个连接单元为平移连接属性Slot和旋转连接属性Cardan，可查阅附图5直线顶管转变为曲线顶管连接属性图；顶管直线转变为曲线的过程中，随着顶进距离的增加，管节之间的夹角在增大，且呈现出线性关系。

第七步，曲线顶管多体系统动学模型顶管顶进力预计计算；钢管采用壳单元，选用S4R类型单元进行网格划分，六节管节共划分12792个单元；连接属性采用平移连接属性Join和旋转连接属性Revolute，该连接属性仅在管轴方向有旋转自由度。为了能更好地反映钢管顶进过程中直线、曲线顶管在顶进过程中的运动特性和受力特性，在建立钢管顶进多体力学模型中选取的不同的连接器单元及运动约束关系来反映顶管接头的受力和变形，运动约束。直线、曲线顶管接头连接器单元选择如表2所示，其中向接头单元的位置如图6所示。由于摩擦力的影响，每增加一节管节，推力值增加为34579N，且与直线顶管增加值相近；管节之间的相对转角在顶进时表现出动态变化，但相对转动角都非常小，最大转角为1.16×10^-5度，比直线顶管小很多；管道在曲线顶进过程中，为保持管节轨迹为曲线，管节之间会产生水平力，且呈现出水平力两头大中间小的规律；在长距离顶管中，单位面积摩擦阻力是影响顶力大小的主要因素。

表2钢管接头连接器单元的运动约束属性

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims

1.顶管顶进过程的多体系统动力学数值模拟方法，其特征在于，包括下述步骤：

2.根据权利要求1所述的顶管顶进过程的多体系统动力学数值模拟方法，其特征在于，步骤S1中，确定直线、曲线顶管顶进过程中顶管上的荷载分布并计算荷载具体为：

3.根据权利要求1所述的顶管顶进过程的多体系统动力学数值模拟方法，其特征在于，步骤S2中，动力学普遍方程为：

\frac{d}{d t} (\frac{\partial T}{\partial \overset{\cdot}{q}}) - \frac{\partial T}{\partial q} = Q;

其中，q为广义坐标向量；T为系统动能；Q为广义力；

系统动能T由广义速度表示为：

T = \frac{1}{2} {\overset{\cdot}{q}}^{T} M \overset{\cdot}{q};

其中，q为广义速度向量；M为质量矩阵；

物体的动力学控制方程为：

M^{i} {\overset{\cdot\cdot}{q}}^{i} + K^{i} q^{i} = Q_{F}^{i} + Q_{V}^{i};

多体系统的动力学控制方程为：

M \overset{\cdot\cdot}{q} + K q + C_{q}^{T} = Q_{F} + Q_{V}

其相应的约束方程为：

C(q,t)＝0

4.根据权利要求1所述的顶管顶进过程的多体系统动力学数值模拟方法，其特征在于，步骤S3具体为：

S32、在INTERACTION模块中，设置连接单元、连接属性和约束；

5.根据权利要求1所述的顶管顶进过程的多体系统动力学数值模拟方法，其特征在于，步骤S5中，直线顶管采用Slide-Plane和旋转连接属性Cardan；计算出每增加一个管节推力值增加为一个定值，主要原因是摩擦力的影响；在直线顶管顶进过程中管节与管节之间会出现一定的相对转动，且相对转角处于动态变化的过程；管节之间连接时允许一定的转动自由度，可以减少管节之间的转动弯矩。

6.根据权利要求1所述的顶管顶进过程的多体系统动力学数值模拟方法，其特征在于，步骤S6中，还包括确定连接属性的步骤，两管节之间采用四个连接单元对称分布在管周，曲线内侧连接单元为平移连接属性Join和旋转连接属性Cardan,其余连接单元为平移连接属性Slot和旋转连接属性Cardan；顶管直线转变为曲线的过程中，随着顶进距离的增加，管节之间的夹角在增大，且呈现出线性关系。

7.根据权利要求1所述的顶管顶进过程的多体系统动力学数值模拟方法，其特征在于，步骤S7中，连接属性采用平移连接属性Join和旋转连接属性Revolute；由于摩擦力的影响，每增加一节管节，推力值增加为一个定值，且与直线顶管增加值相近；管节之间的相对转角在顶进时表现出动态变化，但相对转动角都非常小，比直线顶管小很多；管道在曲线顶进过程中，为保持管节轨迹为曲线，管节之间会产生水平力，且呈现出水平力两头大中间小的规律；在长距离顶管中，单位面积摩擦阻力是影响顶力大小的主要因素。