【2017年整理】浅论断层运动下管图接触面费线性反应分析.docVIP

　摘要：对埋地管道结构分析，考虑管道与土的相互作用问题是非常必要的。其相互作用问题包括管土界面滑移、分离现象及界面处接触应力的确定。采用非线性接触的处理方法建立了管土相互作用分析模型，管道及其周围土体分别以空间薄壳单元和实体单元进行离散，采用有限元手段分析了断层错动下管土界面的非线性反应情况。通过算例，分析了位错量、埋深、管道径厚比及场地条件对管土界面非线性反应的影响；此外对管土界面接触应力进行了分析，发现近断层处管土发生强烈的非线性反应，在距离断层较远处，管土切向摩擦应力的是以一种“抛物线”式的方式递减，而非以线性递减的方式减少。　　关键词：埋地管道，非线性接触，断层运动，接触应力 ? ??? 1.引言??? 埋地管线是重要的生命线工程，影响人们日常的生产生活。西气东输和川气东送同三峡工程、青藏铁路、南水北调并称为中国五大工程，建成后将形成覆盖全国的天然气管道网络。这些重大的天然气管道工程遍及全国各地的自来水、供燃气、排水等管道工程，这些管道工程构成了维系现代社会正常运转的基石。管线工程多埋于地下，形成的网络系统庞大复杂，需要覆盖广阔的地域，对于一些不利甚至危险的地域也必须经过。例如，西气东输工程跨越了多个地震活动地区，需要穿越100多个断层。埋地管线担负着输送生产和生活物质的功能，它需要覆盖广阔的区域，遇到不利甚至危险的地质条件(断层、液化、震陷、滑坡等)也无法避让而只能穿过。这些地质灾害将导致地面产生永久变形，严重威胁着埋地管线工程的安全性，有可能造成管线系统服务功能的下降甚至丧失，对社会生产、生活产生影响。更为严重的是，油、气、水等泄露不仅经济损失巨大，而且可能导致环境污染，甚至引发火灾等次生灾害。因此，研究在地面永久变形下埋地管道的工作性状，是不得不面临和需要解决的重要问题之一。　　埋地管道在上覆土层和外荷载等作用下，管道将因受力而变形，管周土体既是作用在管道上的荷载，同时又是增强管道强度和刚度的一种介质。因此在研究地埋管道的工作机理时，必须把管道周围一定范围内的土体作为结构的一部分加以考虑，即考虑管土相互作用问题。地面永久变形下管土相互作用问题主要包括管土接触界面的相对滑移、分离及闭合等非线性反应现象以及界面处分布接触应力的确定[1-3]。　　地面大变形下埋地管道的反应分析研究从最初的索及梁模型发展到壳模型，是一个巨大的改进；但是目前处理管土相互作用时，多将土体对管道的作用简化为一系列弹簧参数[4-6]来替代，而土体的强非线性也导致了弹簧参数的难确定。本文采用非线性接触[7-12]的处理方法建立了管土相互作用模型，通过有限元计算分析了位错量、埋深、管道径厚比及场地条件对管土界面非线性反应的影响，并对管土界面上接触应力及反应变化规律进行了简单的分析研究。 　　2.有限元模型及参数选择??? 大量的动力和静力实验及研究分析表明管线一般情况下受震和受到断层作用运动形式和原因是不同的，跨越断层埋地管线运动主要是由于地表的永久性位移产生，其动力放大作用微小，可以忽略管道惯性力的影响，因此对管线的地震反应可以采用拟静力方法进行研究。　　建立了地面永久变形下管土相互作用的三维有限元分析模型，如图1所示。管道及周围土体从半无限地球介质中取出，采用非线性接触单元来模拟管土之间的相对滑移、分离及闭合等非线性现象。由模型可知，计算区域包括管道结构及周围土体介质，由6面边界条件组成，除了地表为自由面外，其他5个面为人工边界条件。边界条件为：地表为自由面，左右端面、模型前后两面为水平约束，左侧底面采用竖向约束，右侧底面为位移输入边界；管道两侧为自由边界。　　圆柱形中空薄壁管道采用四节点薄壳单元，土体采用八节点三维实体单元，考虑初始应力场的影响，采用线性位移加载方式来表征断层位错作用，对管土相互作用体系进行非线性分析，研究管土交界面上非线性接触的反应过程。　　取管径为0.4m的钢质管道作为案例分析对象，上、下盘土体纵向长度均取为160m，即为管径的400倍长度；土体宽度取为8m（管径的20倍长度）；管道计算长度取为324m。在断层两侧各20m区域内网格进行细分。　　模型中钢质管道采用三折线本构模型，其弹性阶段弹性模量MPa，当应力达到418MPa管道开始进入弹塑性阶段，当应力与塑性应变分别达到516MPa、0.016时完全进入塑性状态，泊松比52.110E=×0.3υ=。土体采用Mohr-Coulomb模型。 ??? 3.管土界面反应现象分析??? 为了了解认识地震导致的地面永久变形下管土交界面上非线性接触反应现象，主要包括管土界面的滑移及管土相对分离现象，对正断层作用下管土相互作用反应进行了参数影响分析。以接触面下侧的竖向位移为指标来衡量管土间的相对分离反应，以接触面上侧的轴向位移为指标来表征管土间的相对滑移行为。依