盾构隧道接头构造的模拟方法

1.一种盾构隧道接头构造的模拟方法，其特征是按如下步骤实施：
(1)设置接头单元：将接头接触面离散为多个节点，节点之间通过模拟弹簧的接头单元连接，环向上接头单元有多层，分别模拟实际接头中的传力衬垫、螺栓和密封垫，同层接头单元类型相同；
(2)设置初始弹簧刚度、边界条件、接头初始变位δ1；
(3)施加荷载，计算固定时步；
(4)查看第一个接头单元；
(5)根据该接头单元类型设置弹簧法向刚度；
(6)判断是否设置完所有接头单元，如果是，继续下个步骤(7)；如果不是，查看下一个接头单元，回到步骤(5)；
(7)计算固定时步
(8)计算接头实际变位δ2；
(9)比较δ1和δ2，如果两者之差在允许范围to1内，完成计算，算出隧道的接头变形；
如果两者相差超过允许范围，根据δ2计算隧道的附加荷载，并回到步骤(3)重新施加荷载。
2.根据权利要求1所述的盾构隧道接头构造的模拟方法，其特征是：环向上接头单元有四层，从内向外第一层弹簧模拟内侧传力衬垫，第二层弹簧模拟螺栓，第三层弹簧模拟外侧传力衬垫，第四层弹簧模拟弹性密封垫。
3.根据权利要求2所述的盾构隧道接头构造的模拟方法，其特征是步骤(5)中根据接头单元类型设置弹簧法向刚度的步骤如下：
(1)查看接头单元的法向向量，根据法向向量的纵向分量判断接头单元属于环缝还是纵缝；
(2)根据坐标判断接头单元属于哪层弹簧类型；
(3)查看弹簧应力，根据弹簧应力设置弹簧法向刚度。
4.根据权利要求3所述的盾构隧道接头构造的模拟方法，其特征是设置弹簧法向刚度的方法如下：对于传力衬垫层弹簧，若法向弹簧受压，则取实际压缩刚度，否则取零；对于螺栓层弹簧，若法向弹簧受压，则取实际压缩刚度，否则取实际拉伸刚度；对于密封垫层弹簧，若法向弹簧受压，则取实际压缩刚度，否则取零。

盾构隧道接头构造的模拟方法 \n技术领域\n[0001] 本发明涉及一种属隧道结构模型计算领域，具体涉及一种盾构隧道接头构造的模拟方法。 \n背景技术\n[0002] 盾构隧道由管片拼装而成，在隧道的长期运营过程中，由于外界各种因素的影响，隧道沿着纵向会发生不均匀沉降，隧道不均匀沉降引起衬砌结构接头变形，进而引起隧道的渗水漏泥，目前隧道的渗漏水是运营隧道的主要病害之一。因此，如何合理模拟隧道纵向不均匀沉降引起的接头变形对研究运营隧道的性能十分重要，其中接头的模拟合理与否决定了接头变形计算的合理性和准确性。目前关于隧道纵向变形的模拟主要有以下几种方法： \n[0003] (1)纵向梁-弹簧模型 \n[0004] 采用梁单元模拟衬砌环，以弹簧的轴向、剪切和转动效应模拟环间接头和螺栓。纵向梁-弹簧模型将隧道作为一根梁来考虑。因此，不能从空间上描述接头的变形特征以及对隧道防水的影响。 \n[0005] (2)纵向等效连续化模型 \n[0006] 纵向等效连续化模型将隧道在横向上视为一均质圆环，在纵向以刚度等效的方法把结构转化为由接头和管片组成的梁-接头不连续模型与管片环间剪切模型，即采用梁单元来模拟管片，而借助一维接头单元来模拟管片间的接头效应。研究隧道的纵向沉降有一定的适用性，但是不能描述隧道接头的变形。 \n[0007] (3)梁-接头不连续模型 \n[0008] 梁-接头不连续模型考虑到接头变形的不连续性，将管片离散为梁单元，管片间的接头考虑成接头单元，借助一维接头单元来模拟管片间的接头效应。 \n[0009] (4)实体模型 \n[0010] 随着计算机技术的发展，利用数值计算方法将隧道衬砌按照实体单元进行计算，接头部位采用接触单元来进行，接头螺栓采用梁单元或者杆单元进行模拟，目前的大型有限元分析软件均可以进行混合单元的计算和分析。 \n[0011] 现有的计算方法主要存在以下两个方面的问题：1、上述计算方法主要针对隧道的纵向受力和纵向沉降展开，因此不能合理地描述接头的变形特征。2、接头模拟参数不容易确定。目前接头的计算模型一般需要输入接头的转动刚度和剪切刚度参数，由于衬砌整环试验或管片 接头试验成本较大，参数的确定一般凭经验或工程类比确定，存在较大的任意性。为了确定这些参数，一些学者根据接头构造建立了管片接头的力学模型，计算得到接头在荷载作用下的转动刚度，这些力学模型对于理解接头在已知荷载作用下的受力特征和变形模式具有重要意义。在实际衬砌结构中，不同位置的接头受到的荷载是未知并且不断变化的。不同的荷载组合条件下，各接头虽然构造相同，接头的刚度参数却存在较大的差别。\n因此，直接应用接头分析得到的参数计算衬砌结构的内力和接头变形会产生较大的误差。 发明内容\n[0012] 本发明要解决的技术问题是提供一种盾构隧道接头构造的模拟方法，能更准确地反映管片变形和受力状况。 \n[0013] 本发明采用以下技术方案： \n[0014] 一种盾构隧道接头构造的模拟方法，其特征是按如下步骤实施： [0015] (1)设置接头单元：将接头接触面离散为多个节点，节点之间通过模拟弹簧的接头单元连接，环向上接头单位有多层，分别模拟实际接头中的传力衬垫、螺栓和密封垫，同层接头单元类型相同； \n[0016] (2)设置初始弹簧刚度、边界条件、接头初始变位δ1； \n[0017] (3)施加荷载； \n[0018] (4)查看第一个接头单元； \n[0019] (5)根据该接头单元类型设置弹簧法向刚度； \n[0020] (6)判断是否设置所有接头单元，如果是，继续下个步骤(7)；如果不是，查看下一个接头单元，回到步骤(5)； \n[0021] (7)计算接头实际变位δ2； \n[0022] (8)比较δ1和δ2，如果两者之差在允许范围to1内，完成计算，算出隧道的接头变形；如果两者相差超过允许范围，根据δ2计算隧道的附加荷载，并回到步骤(3)重新施加荷载。 \n[0023] 进一步地，环向上接头单位有四层，从内向外第一层弹簧模拟内侧传力衬垫，第二层弹簧模拟螺栓，第三层弹簧模拟外侧传力衬垫，第四层弹簧模拟弹性密封垫。 [0024] 进一步地，步骤(5)中根据接头单元类型设置弹簧法向刚度的步骤如下： [0025] (1)查看接头单元的法向向量，根据法向向量的纵向分量判断接头单元属于环缝还是纵缝； \n[0026] (2)根据坐标判断接头单元属于哪层弹簧类型； \n[0027] (3)查看弹簧应力，根据弹簧应力设置弹簧法向刚度。 \n[0028] 进一步地，设置弹簧法向刚度的方法如下：对于传力衬垫层弹簧，若法向弹簧受压，则取实际压缩刚度，否则取零；对于螺栓层弹簧，若法向弹簧受压，则取实际压缩刚度，否则取实际拉伸刚度；对于密封垫层弹簧，若法向弹簧受压，则取实际压缩刚度，否则取零。 [0029] 本发明为了建立更真实反映接头的力学模型，按照接头构造排列顺序由多组弹簧来模拟，每组弹簧代表一类接头构造材料，根据得到的节点弹簧应力或法向位移在计算过程中动态调整弹簧参数，将接头模型应用到整环衬砌结构计算中。本发明采用考虑隧道接头构造特点的模型来模拟接头在隧道纵向沉降过程中的变形特点，和其他接头模型相比较，具有如下特点： \n[0030] (1)和纵向梁-弹簧模型和梁-接头不连续模型相比，采用考虑接头构造的隧道接头模拟方法，接头的刚度参数比较容易确定，可以根据弹簧的受力状态调整弹簧的刚度，使计算结果更加合理； \n[0031] (2)和实体模型相比较，采用考虑接头构造的隧道接头模拟方法，能够有效提高计算效率，而且可以避免接触单元参数确定的随意性； \n[0032] (3)采用考虑接头构造的隧道接头模拟方法，不仅能合理的计算隧道接头的张开量，并且能够合理反映隧道接头的错台量，以及管片之间的相互位移和约束； [0033] (4)采用考虑接头构造的隧道接头模拟方法，能够考虑接头刚度变化对隧道管片受力的影响，能够比较合理反映隧道的受力； \n[0034] (5)采用考虑接头构造的隧道接头模拟方法，解决了纵向梁-弹簧模型和纵向等效连续化模型中的纵横向刚度不匹配的问题。 \n附图说明\n[0035] 图1为本发明的流程图。 \n[0036] 图2为单个组合弹簧的示意图。 \n[0037] 图3为接头力学模型。 \n[0038] 图4、图5为管片接头受力前后的构造示意图。 \n[0039] 图6为节点对应的接触面面积。 \n[0040] 图7为弹簧法向刚度参数的设置和调整原则。 \n[0041] 图8为内侧传力衬垫受拉后的状态示意。 \n[0042] 图9为实例环向接头的模拟。 \n具体实施方式\n[0043] 参见图1，实施步骤如下： \n[0044] (1)开始； \n[0045] (2)设置接头单元，采用四组弹簧模拟管片四层构造，由内向外第一层弹簧模拟内侧传力衬垫D，第二层弹簧模拟螺栓C，第三层弹簧模拟外侧传力衬垫B，第四层弹簧模拟弹性密封垫A，每组弹簧包括多个子弹簧，环向上子弹簧的数量可根据管片的尺寸和计算精度确定，轴向上的数量同样由计算精度确定； \n[0046] (3)设置初始弹簧刚度，包括法向刚度和剪切刚度，根据实际条件设置边界条件，拟定接头初始变位δ1； \n[0047] (4)施加荷载，计算固定时步； \n[0048] (5)查看第一个接头单元法向向量； \n[0049] (6)根据法向向量的纵向分量判断接头单元属于环缝还是纵缝； \n[0050] (7)根据坐标判断接头单元属于哪层弹簧类型； \n[0051] (8)查看弹簧应力，根据弹簧应力设置弹簧法向刚度； \n[0052] (9)判断是否设置完所有接头单元参数，如果是，继续下个步骤(10)；如果不是，查看下一个接头单元法向向量，回到步骤(6)； \n[0053] (10)计算固定时步； \n[0054] (11)计算接头实际变位δ2； \n[0055] (12)比较δ1和δ2，如果两者之差在允许范围to1内，完成计算，算出隧道的接头变形；如果两者相差超过允许范围，根据δ2计算隧道的附加荷载，并回到步骤(4)重新施加荷载。 \n[0056] 本发明的接头单元模拟弹簧，按照实际接头构造分为四组弹簧，见图2～5。每组弹簧代表一类接头构造材料，每组弹簧再由一系列子弹簧构成，子弹簧的数量根据管片的尺寸和计算精度确定，每个子(组合)弹簧由剪切弹簧ks和法向弹簧kn组成。将接头接触面离散为一系列对应的节点p，节点p之间通过组合弹簧产生相互作用，每个节点对应的接触面的面积称为特征面积，见6图中阴影部分。每个节点对应一个组合弹簧，沿环向共有四层弹簧。 \n[0057] 接头单元的类型由接头单元法向向量的纵向分量来确定。无论是环缝接头单元还是纵缝 接头单元，在判断弹簧法向刚度时，都需要确定该单元所处四层弹簧中的哪一层以及查看法向应力。具体判断原则如下，参见图7：对于传力衬垫层弹簧，若法向弹簧受压，则取实际压缩刚度，不然取零(见图8)；对于螺栓层弹簧，若法向弹簧受压，则取实际压缩刚度，不然取实际拉伸刚度；对于密封垫层弹簧，若法向弹簧受压，则取实际压缩刚度，不然取零。 \n[0058] 调整弹簧参数的方法和设置一样，对于传力衬垫和弹性密封垫，如果出现弹簧拉伸现象，弹簧退出工作状态。 \n[0059] 本发明基于衬砌结构接头的构造特征，建立了衬砌结构接头的力学模型，根据得到的节点弹簧应力或法向位移在计算过程中动态调整弹簧参数，将接头模型应用到整环衬砌结构计算中。可以同时考虑衬砌的环向、纵向结构性能，计算衬砌管片和接头的内力和变形，更加准确地反映管片的变形和受力状况。 \n[0060] 下面通过实例进一步描述本发明。 \n[0061] 隧道衬砌圆环内径5.5m，管片厚度0.35m，环宽1m。一环内衬砌分6块，其中1块封顶块，2块邻接块，2块标准块，1块封底块，衬砌环间以16根M30的纵向螺栓相连。管片混凝土强度等级C50，弹性模量3.45E10pa，泊松比0.2。传力衬垫材料的压缩刚度2.48E11Pa/m；弹性密封垫材料的压缩刚度为9E7Pa/m。 \n[0062] 以环向接头为例介绍接头的模拟方法： \n[0063] 管片单元将衬砌结构在环向96等分，每个单元对应的圆心角θ为3.75°，沿径向分为3个单元，对应4层弹簧，因此环向接触面采用了384个组合弹簧(每个组合弹簧包括压缩弹簧和剪切弹簧)进行模拟，如图9所示。由内向外第1层弹簧模拟内侧传力衬垫，第\n2层弹簧模拟螺栓，第3层模拟外侧传力衬垫，第4层弹簧模拟弹性密封垫。 [0064] 1、弹簧压缩刚度的确定 \n[0065] 以环缝内侧传力衬垫为例，衬垫的实际面积为1.59m2，对应衬垫弹簧的节点特征\n2\n面积为1.02m，节点的特征面积对弹簧对应的单元面积，则计算模型中衬垫的压缩刚度参\n2 2\n数为2.48E11Pa/m×1.59m÷1.02m ＝3.87E11Pa/m。 \n[0066] 2、弹簧剪切刚度的确定 \n[0067] 考虑到接头发生错台时的变形特点，接头的剪切刚度可以取为同一数值，针对实例中的隧道剪切刚度可取为9.78E7Pa/m。