穿黄隧洞纵向变形三维有限元分析-计算机应用-长江科学院.PDFVIP

　 第２７卷第７期 长　江　科　学　院　院　报 Ｖｏｌ．２７　Ｎｏ．７　　 　　２０１０年 ７月 ＪｏｕｒｎａｌｏｆＹａｎｇｔｚｅＲｉｖｅｒＳｃｉｅｎｔｉｆｉｃＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ Ｊｕｌ．２０ １０　 　 　　文章编号：１００１－５４８５（２０１０）０７－００６０－０５ 穿黄隧洞纵向变形三维有限元分析 １，２ １，２ １，３ 谢小玲 ，苏海东 ，林绍忠 （１．水利部 水工程安全与病害防治工程技术研究中心，武汉　４３００１０； ２．长江科学院材料与结构研究所，武汉　４３００１０；３．长江科学院非连续变形分析实验室，武汉　４３００１０） 摘要：穿黄隧洞深埋于河床摆幅大、冲淤变化剧烈的黄河河道，地处７度地震区，是一个由数百个环拼装联结的非 连续变形柔性结构。河床冲淤变化以及地震荷载对隧洞纵向变形以及纵向变形缝的影响问题，是隧洞方案是否可 行的关键技术问题之一。为此，采用自编的考虑接触的三维非线性有限元分析程序，计算隧洞结构纵向变形、变形 缝张开度、错动位移等。针对隧洞模型计算范围大，接触面多等特点，通过采取多项措施控制计算规模，保证计算 精度，结合程序的方程组快速求解方法，解决了穿黄隧洞纵向变形这一技术难题。计算结果表明，在最不利荷载组 合作用下，隧洞变形缝的张开度、错动位移和相应跨缝螺栓应力以及混凝土应力均在设计控制的范围内，满足设计 要求。 关　键　词：穿黄隧洞；纵向变形；接触非线性；三维有限元 中图分类号：ＴＵ３１１．４；ＴＶ３１４　　　文献标识码：Ａ 面接触状态以及螺栓应力等，为隧洞设计提供依据。 程序中的接触算法和改进的对称逐步超松弛预处理 １　概　况 共轭梯度迭代（ＳＳＯＲ－ＰＣＧ）的方程组解法，解决了 南水北调中线工程穿黄隧洞穿越的河段属典型 计算模型规模大、接触面多等造成的计算量庞大的 游荡性河道，河势演变复杂，河床冲淤变化剧烈，主 关键技术问题。 槽摆幅较大。河床主槽摆动产生的冲淤荷载以及地 震荷载等因素，可能使隧洞纵向发生不均匀的沉降 ２　隧洞介绍 或回弹变形。而过大的不均匀变形会对衬砌结构、 纵向变形缝止水等产生不利影响，因此隧洞设计必 穿黄隧洞设计为双线布置，单线长数千米，轴线 须将隧洞纵向的不均匀变形控制在允许的范围内。 间距３０多米，埋深在３０ｍ以下，穿越的河床地基由 对于隧洞这种非连续变形柔性结构的受力特性 ２ 亚粘土、粉细砂、中砂、粘土岩组成，其地质系为Ｑ， ４ 分析，国内外学者也进行了很多研究，提出了多种计 １ Ｑ，Ｑ，Ｎ，其中亚粘土、粉细砂的力学性质较差。 ４ ２ 算模型［１，２］，如梁－弹簧模型、梁－接头模型、纵向等 隧洞衬砌结构设计为双层，分外衬和内衬，外衬 效连续化模型以及三维骨架模型等，但这些模型应用 由管片拼装成整环，内衬为现浇的整环结构。盾构 于实际工程还存在很多不足，或不能准确模拟管片接 施工过程中，首先由盾壳撑住已开挖的周围土体，管 头的非线性特性，或模