大口径X80管道沉管应力响应规律与临界深度分析.pdf

“十四五”期间，国内大口径高钢级长输管道建设将进入高峰期。沉管下沟是大口径管道铺设的主要施工方式之一，当管道完全沉管施工时，所需管沟开挖长度、特别是管道应力状态难以估计，可能造成管沟开挖长度不足或局部应力过大等情况，难以保证施工进度与质量。应用非线性有限元方法，以D1219mm大口径X80直管段管道为研究对象，分析不同土壤、管道壁厚、支墩位置、有无支墩及下沟深度5 项影响因素对直管段管道沉管敷设完全沉管时悬空长度与管顶应力状态的影响规律。结果表明：不同土壤类型对管道完全沉管敷设时悬空长度与管顶应力状态有一定影响，但影响较小；当管沟深度在5 m以内时，壁厚与支墩位置变化时，悬空长度和管顶应力的变化趋势相同，且完全沉管时悬空长度与最大拉应力基本一致；相同管沟深度完全沉管时，管道悬空长度与管顶最大应力在有支墩情况下大于无支墩情况；沉管敷设管沟越深，悬空长度与管顶应力越大。对于D1219mmX80管道，无支墩情况下，双侧单向与双侧背向沉管临界深度分别为6.7 4m和6.8 0 m；有支墩情况下，双侧单向与双侧背向沉管临界深度分别为6.2 5 m和6.31m。分析结果可为管道铺设施工作业提供参考。关键词：大口径管道；沉管下沟；有限元分析；应力状态；临界下沟深度中图法分类号：TE973D0I:10.19459/ki.61-1500/te.2023.05.007Analysis of Stress Response Trends and Critical Depth ofLarge Diameter X80 Pipe Lowering-inSHI Tong,WANG Yanbing,JIN Guangyi,LIU Xiaoben,ZHANG Hong,HUANG Qiyu(1.National Engineering Laboratory for Pipeline Safety/MOE Key Laboratory of Petroleum Engineering/Beijing Key Laboratory ofUrban Oil and Gas Distribution Technology,China University of Petroleum-Bejing,Bejing 102249,China;Abstract:During the 14th Five-Year Plan period,the construction of large-diameter high grade long distance pipelines will enter the peakperiod.Pipe lowering-in is one of the main ways to lay large-diameter pipelines.When the pipeline is completely lowered-in,the requiredtrench excavation length,especially the pipeline stress state is dfficult to estimate,which may lead to insufficient trench excavation lengthor excessive local stress,and it is difficult to ensure the construction schedule and quality.With the D1 219 mm large diameter X80 straightpipe as the research object,the nonlinear finite element method is used to analyze the influence of the five influencing factors,including dif-ferent soil,pipe wall thickness,the position of pipe sleeper,the presence or absence of pipe sleeper,and the depth of the trench,on thehanging length and the stress state of the pipe top when the straight pipe is completely lowered in.The results show that different soil types初投稿收稿日期：2 0 2 3-0 2-2；修改稿收稿日期：2 0 2 3-0 3-0 1*基金项目：国家重点研发计划重点专项项目“中俄管道重大风险防控与安全保障关键技术”课题一“高压力高钢级管道失效机理与全生命周期可靠性评价技术研究”（编号：2 0 2 2 YFC3070101）；国家自然科学基金项目“逆断层作用下X80管道屈曲演化与韧性破损机理研究”（编号：5 2 0 0 4314）；北京市科协“青年人才托举工程”项目（编号：BYESS2023145）；国家管网科学研究与技术开发项目“高钢级管道环焊缝失效机理研究”（编号：WZXGL202105）；国家管网科学研究与技术开发项目“高钢级管道环焊缝缺陷检测评价技术研究”（编号：WZXGL202104）第一作者简介：石彤，男，19 9 7 年生，中国石油大学（北京）在读博士研究生，主要从事油气装备失效分析与完整性管理方向的研究工作。E-mail;shitong0409 文献标识码：A2.Jilin Enthalpy Wing Technology Co.,Ltd.,Changchun,Jilin 130033,China)文章编号:2 0 9 6-0 0 7 7(2 0 2 3)0 5-0 0 32-0 92023年第9 卷第5 期have certain influence on the suspended length and stress state of pipe top when the pipe is completely lowered in,but the influence issmall.When the trench depth is less than 5 m,the variation trend of the suspended length and the stress at the top of the pipe is the samewhen the wall thickness and the position of the supporting sleeper change,and the suspended length and the maximum tensile stress are ba-sically the same when the pipe is completely lowered in.When the pipe is completely lowered-in at the same trench depth,the pipe suspen-sion length and the maximum stress at the top of the pipe with pipe sleeper are greater than those without pipe sleeper.The deeper the pipetrench is laid,the greater the suspended length and the pipe top stress.For D1 219 mm X80 pipeline without pipe sleeper,the criticaldepth of bilateral unidirectional and bilateral back-to-back pipe lowering-in is 6.74 m and 6.80 m,respectively.In the case of pipe sleep-er,the critical depths of bilateral unidirectional and bilateral back-to-back pipe lowering-in are 6.25 m and 6.31 m,respectively.The a-nalysis results can provide reference for pipeline lowering-in construction.Key words:large diameter pipeline;pipe lowering-in;finite element analysis;stress response;critical lowering-in depth0引言“十四五”期间，随着对油气需求量的不断加大，我国大口径高钢级长输管道建设将进人高峰期1-2 。而作为管道建设重要的一环，下沟铺设的安全性就显得尤为重要。沉管法是利用管道自身重力将管道自然降落到管沟内的施工方法，因为沉管施工具有征地协调量小、操作简便、施工作业效率高等优点，在我国管道建设中得到了越来越多的应用3-6 。但是，由于沉管施工过程中存在较大的管沟落差，管道发生弯曲变形,容易在局部产生较大的拉压应变，严重时将导致拉裂或屈曲失效。针对沉管下沟管道铺设方法，国内外学者开展了一系列研究。在沉管施工方法与解析计算研究方面，康万平、刘长江、李艳华等人7 研究了不同施工环境下管道铺设方法的选择以及沉管下沟在油气长输管道施工中的应用。张英奎等人8 介绍了西气东输管道工程在水网地区采用大口径管道施工的特点及方法。对施工便道修筑和作业带加固、河流和水塘及道路的穿越、管材运输等进行了论述，为今后类似的管道施工建设积累了一定的经验。季蓓蕾、刘啸奔等人9 研究了不同沉降形式对管道轴向应力状态的影响，掌握大口径管道的轴向应力应变状态，保障了管道的安全运行。王付会等人10 将模型简化成超静定梁，最终推导出不同深度管道最大应力公式。夏梦莹、张宏等 基于弹性梁理论,给出了挖沟过程中管道挠曲线微分方程，通过边界条件与协调方程,给出了考虑几何非线性的挖沟过程中管道应力与变形解析计算方法。刘啸奔等人12 1 提出一种改进的位移变形作用下海底管道应变解析分析方法，根据线性强化模型考虑了管材的非线性本构关系，最终基于平衡方程和迭代计算，可以精确计算管道应力应变。在数值计算方面，程浩力等人【13 对高地下水位段采用单侧沉管下沟的施工工艺进行了有限元模型验证，得到了小口径油石彤等：大口径X80管道沉管应力响应规律与临界深度分析气管线在高地下水位段施工的相关结论。孟佳等14 简化出管道和土体相互作用的物理模型，利用ANSYS软件进行沉管下沟的开挖长度分析。谷青悦等人15 使用有限元模型对管道沉管下沟所受应力进行模拟计算和比较分析，论证该长输管道项目采用管段端部沉管下沟或管段中间沉管下沟的施工方法。刘玉卿等人16-17 利用有限元方法分析了中俄东线采用沉管下沟施工作业方式的可行性。综上所述,虽然国内外对沉管施工方法进行了一定的探究，综合采用解析法和有限元方法进行了管道沉管下沟过程中管道过程的仿真分析，但未得到X80管道沉管下沟的临界悬空长度与沉管深度，现场施工缺乏直接指导。针对此，本文基于非线性有限元法，建立我国D1219mmX80直管段管道沉管下沟应变数值计算模型，通过分析X80管道在沉管下沟时的变形过程，得到管道完全沉管时悬空长度与管沟深度的关联关系，进而得到完全沉管施工临界的悬空长度与管沟深度。研究成果可以直接指导大口径X80管道管沟开挖过程的控制，保障沉管施工安全、高效进行。1沉管下沟简介1.1沉管过程简介沉管下沟指的是沿管道开挖管沟,利用管道自身重力作用缓慢将管道自然下沉到管沟内的施工方法。Q/SYGDJ03872014油气输送管道沉管下沟施工规范中将沉管下沟分为单侧沉管下沟和双侧沉管下沟，其中双侧下沟的方式在工程实际中应用更广泛。双侧沉管下沟可分为双侧单向沉管下沟及双侧背向沉管下沟。双侧单向沉管下沟是从管段的一端开始，逐步向管段另一端进行管沟开挖和沉管下沟作业，如图1(a)所示;双侧背向沉管下沟则是从管段的中间一点开始，同时向两端进行管沟开挖和沉管下沟作业，如图1(b)所示。当管道完全接触管沟沟底时，视为完全3334沉管。挖洲进前进力情道合地面下沟管沟深度(a)双侧单向沉管下沟挖掘机前进方向管道地面管沟(b)双侧背向沉管下沟图1双侧沉管下沟示意图1.2沉管施工中管道应力控制要求沉管过程中因管道会产生较大范围的垂向弯曲，进而形成较大的轴向弯曲应力，根据Q/SYGDJ0387一2014油气输送管道沉管下沟施工规范的规定，我国管道沉管敷设时管道的许用应力应满足：0。=h-0 a =0.8g,式中：。为管道中的等效应力,MPa;Qh为管道环向应力,MPa;oa为管道轴向应力,MPa；为管道设计的许用应力,MPa;。为管材的最低屈服强度,MPa。2管道沉管下沟有限元模型应用非线性有限元分析软件，建立D1219mmX80直管段管道沉管下沟的数值计算模型，分析X80直管段管道在沉管下沟施工时的变形过程，经计算得到管道完全沉管时悬空长度与应力响应规律，并计算得到临界深度，可直接指导沉管施工管沟开挖并校核油气长输管道沉管下沟的安全性。下沟过程中，管道长度达数百米，针对管道沉管下沟数值计算模型，将管道本体简化为三维薄壁管，使用有限元软件ABAQUS中的PIPE31空间线性管单元进行模拟。管道材料采用线弹性,其弹性模量为2 10 GPa、泊松比0.3,最小屈服强度5 5 5 MPa。考虑到下沟管段足够长，管道两端自由约束。在下沟敷设过程中，管道受到重力和土体支撑力的作用，重力载荷即为管道自重，以惯性力的方式施加。在下沟过程中，在管沟边缘支点处会对管道进行挤压，产生较大应力集中现象，管道最大拉应力出现在管沟边缘管顶处。计算时,基于建立的双侧单向沉管下沟与双侧背向沉管下沟有限元模型，通过将不同位置土块下移模拟挖掘机开挖产生的沟底深度，模拟完全沉管时悬空长度与管顶最大应力。建立的双侧单向沉管下沟与双侧背向沉管下沟有限元模型如图2 所示。五油管材与仪器S,MisesAngle=180.0000,(1-fraction=-1.000000,2-fraction=0.000000):t3:751e+083:4380+08H3:126e+08沉管方向H:813e+08+2:501e+08E2:188e+08H+:8750+081363e+08+9:377e+026:252e+07H+3:126e+07=+2.403e-05挖掘机前进方向沉管方向下沟立深度2023年10 月(a)双侧单向沉管下沟S,MisesAngle=180.0000,(1-fraction=-1.000000,2-fraction=0.000000)(Avg:75%)+3.582e+08三3:32 4e+08+2:866e+08+2:507e+08+2.149e+08+1:791e+08+1.433e+08E+1:075e+07三+2:16 4e+07+5.790e+00(1)图2 双侧沉管下沟有限元模型3管道应力和悬空长度影响分析3.1不同土壤对管道应力的影响在进行沉管敷设施工时，不同土壤类型会对沉管敷设时管道应力有一定影响。为探究大口径高钢级长输管道不同土壤情况下双侧单向与背向完全沉管施工时管道悬空长度与管顶最大应力状态分布规律，选取D1219mmX80管道，以软黏土、中黏土、硬黏土、松砂、中砂和密砂六种常见的土壤类型及3m沟深施工工况为例，进行数值仿真计算。依据BSEN1997-12004欧洲法规7：土木设计总则中给出的土壤类型,进行土壤参数取值,6 种常见土壤类型弹性模量见表1。表1不同土壤类型弹性模量软黏土中黏土硬黏土松砂1.02.0D1219mm管道3m沉降时管道纵向位移分布曲线如图3所示，双侧单向与双侧背向完全沉管时，管道悬空长度随土壤弹性模量增大而增加，而后趋于稳定，土壤类型为弹性模量较大的密砂时,单向下沟在9 2.7 m左右,背向下沟在9 3.5 m左右。计算结果说明在管径不变条件下，土壤参数对沉管施工中管道悬空长度有一定影(b)双侧背向沉管下沟7.015.0MPa中砂密砂45.075.02023年第9 卷第5 期响，但影响较小。不同土壤参数、同一下沟深度工况管道悬空长度详细计算结果见表2。0.5沟顶0.0-0.5-1.0-2.0-2.5-3.0-3.500.5沟顶0.0-0.5管沟边缘-1.0-1.5-2.0-2.5-3.0-3.5100图3不同土壤管道3m沉降时双侧单向和双侧背向管道纵向位移曲线表2 D1219mm管道不同土壤参数双侧单向、双侧背向沉管下沟悬空长度计算表下沟方式软黏土中黏土硬黏土松砂中砂密砂双侧单向88.7双侧背向89.5D1219mm管道3m沉降时管顶轴向应力分布曲线如图4所示。500400300200松砂100密砂中砂0100-20050100200300400管道坐标/m(a)双侧单向下沟轴向应力分布曲线5004003002001000-100-200100图4不同土壤管道3m沉降时双侧单向和双侧背向管顶轴向应力分布曲线石彤等：大口径X80管道沉管应力响应规律与临界深度分析位于下沟施工边界处，管顶最大拉应力随着土壤弹性模量增大而趋于稳定。土壤类型为密砂时,双侧单向沉管软黏土硬黏土中黏土-2.80管沟边缘松砂中砂-2.90密砂6-2.953-3.00沟底3.05-3.10-3.154260270280290300100200300400管道坐标/m(a)双侧单向下沟纵向位移曲线沟顶软黏土中黏土-2.90硬黏土-2.92松砂3-2.94-2.96中砂密砂位2.9 8包-3.0 0管沟边缘-3.02沟底-3.04-3.06管道坐标/m(b)双侧背向下沟纵向位移曲线89.691.190.792.0-0.80.软黏土管沟边缘中黏土硬黏土管沟迈缘誓沟迈缘0.8 0.一软黏土中黏土3150硬黏土310305300295290松砂285中砂2755280密砂490.49550050551051500400200管道坐标/m管道坐标/m(b)双侧背向下沟轴向应力分布曲线35双侧单向与双侧背向完全沉管时，管顶最大拉应力管道管顶最大应力为30 1.3MPa左右，双侧背向沉管管-2.85道管顶最大应力为30 2.2 MPa左右。由计算结果可知，同一直径、壁厚管道，管道双侧单向、双侧背向完全沉管时管顶最大拉应力随土壤弹性模量增大而增加并趋于稳定。不同土壤弹性模量、同一下沟深度工况管顶轴向管道坐标/m应力详细计算结果见表3。表3D1219mm管道不同土壤参数双侧单向、双侧背向沉管下沟管顶最大拉应力计算表MPa下沟方式软黏土中黏土硬黏土松砂中砂密砂双侧单向295.1双侧背向296.33.2管道壁厚对管道应力的影响410420430440管道坐标/mm91.792.492.792.693.393.53102852755i90195.200,205管道坐标/m297.3298.6为探究大口径高钢级不同壁厚长输管道在双侧单向与背向完全沉管施工时管道悬空长度与管顶最大应力状态分布规律，选取壁厚18.4、2 2、2 7.5 mm的3种D1219mmX80管道，以密砂、下沟深度3、4、5 m施工工况为例，进行数值仿真计算。D1219mm管道3m沉降时管道纵向位移分布曲线如图5 所示。1沟顶0-1-2-302101沟顶0-1-2管沟边缘-310100200300 400500600700管道坐标/m(b)双侧背向下沟纵向位移曲线图5不同壁厚管道3m沉降时双侧单向和双侧背向管道纵向位移曲线双侧单向与双侧背向完全沉管时，不同壁厚管道悬299.9300.918.4mm壁厚22mm壁厚27.5mm壁厚管沟边缘沟底100200管道坐标/m(a)双侧单向下沟纵向位移曲线沟底，300.5301.6300400沟顶-18.4mm壁厚22mm壁厚27.5mm壁厚管沟边缘300.9 301.3302.1302.236空长度基本相当，在9 2 m左右。说明在管径不变条件下,管道壁厚在沉管施工中对管道悬空长度影响较小，同一管径条件下管道完全沉管时，管沟开挖长度基本一致便满足施工要求。不同壁厚、同一下沟深度工况管道悬空长度详细计算结果见表4。表4D1219mm管道不同深度双侧单向、双侧背向沉管下沟悬空长度计算表D1 219 mm D1 219 mm 沟深/m下沟方式18.4 mm管道22mm管道27.5 mm管道悬空长度/m悬空长度/m悬空长度/m双侧单向92.73.0双侧背向双侧单向4.0双侧背向5.0双侧单向双侧背向D1219mm管道3m沉降时的管顶轴向应力分布曲线如图6 所示。双侧单向与双侧背向完全沉管时，管顶最大拉应力位于下沟施工边界处。双侧单向沉管时，不同壁厚管道管顶最大应力为30 1.3MPa左右。双侧背向沉管时，不同壁厚管道管顶最大应力为30 2.2 MPa左右。由此可见，同一直径、不同壁厚管道双侧单向、双侧背向完全沉管时管顶最大拉应力基本相当，说明在管径不变条件下，管道壁厚对完全沉管时管顶最大拉应力影响较小。不同壁厚、同一下沟深度工况，管顶轴向应力详细计算结果见表5。5004003002001000-100-2000(a)双侧单向下沟轴向应力分布曲线500400管沟边缘3002001000-100-2000100200300400500600700管道坐标/m(b)双侧背向下沟轴向应力分布曲线图6不同壁厚管道3m沉降时双侧单向和双侧背向管顶轴向应力分布曲线五油管材与仪器表5 D219mm管道不同深度双侧单向、双侧背向沉管下沟管顶最大拉应力计算表D1 219 mm D1 219 mm 18.4mm管道22mm管道27.5mm管道沟深/m下沟方式管顶最大应力/MPa双侧单向301.33.0双侧背向D1 219 mm 双侧单向4.0双侧背向91.090.093.592.099.098.0100.599.5104.5103.5106.5105.5-0.80,管沟边缘18.4mm壁厚22mm壁厚-27.5mm壁厚100200管道坐标/m管沟边缘-27.5mm壁厚2023年10 月D1 219 mm 管顶最大管顶最大应力/MPa应力/MPa300.8300.7302.2302.0342.0341.8342.5342.2双侧单向382.492.05.097.098.5103.0104.5300400-0.8018.4mm壁厚22mm壁厚301.9341.7342.2382.3382.3双侧背向382.23.3支墩位置对管道应力的影响管道下沟施工时,支墩一般间隔12 m设置,随着管沟开挖，支墩会被人工清除，支墩与管沟的距离在0 12m之间连续变化。为探究大口径高钢级长输管道在不同支墩位置时，双侧单向与背向完全沉管施工管道悬空长度与管顶最大应力状态分布规律，选取D1219mm18.4mmX80管道，以密砂、下沟深度3m、支墩位置距管沟0、4、8、12 m工况为例，进行数值仿真计算。D1219mm18.4mm管道不同支墩位置3m沟深沉降时，管道纵向位移分布曲线如图7 所示，双侧单向与双侧背向完全沉管时，管道悬空长度保持不变。双侧单向与双侧背向沉管时，不同支墩位置管道悬空长度相差小于0.5 m。由此可见，不同支墩位置施工时，管道的悬空长度变化很小。不同下沟深度工况管道悬空长度详细计算结果见表6。0-12-3-50050100150200250300350400450管道坐标/m(a)双侧单向管道纵向位移曲线0.8030.7-10.6距管沟OmS0.4位0.5-2鹿管沟4m0.3F距管沟8 m距管沟12 m0.2-34755005255505756000100200300400500600700管道坐标/m(b)双侧背向管道纵向位移曲线图7 D1219mm18.4mm管道沟深3m不同支墩位置双侧单向和双侧背向管道纵向位移曲线382.10.90.7一距管沟0 m0.6距管沟4m0.5距管沟8 m0.4距管沟12 m0.30.2125150,175,200225管道坐标/m管道坐标/m381.92023年第9 卷第5 期表6 D121918.4mm管道不同支墩位置双侧单向和双侧背向管道悬空长度m支墩距管支墩距管下沟方式沟0 m双侧单向95.5双侧背向96D1219mm18.4mm管道3m沉降时，管顶轴向应力分布曲线如图8 所示。5004003002001000-100-20050（050100150200250300350400450管道坐标/m(a)双侧单向管顶应力分布曲线5000.8054003303003202003101003000290-100280-2004804905005105205305400100200300400500600700管道坐标/m(b)双侧背向管顶应力分布曲线图8 D1219mm18.4mm管道沟深3m不同支墩位置双侧单向和双侧背向沉管管顶应力曲线双侧单向与双侧背向完全沉管时，管顶与管底最大应力位置会随着支墩位置的改变而移动，最大应力在320.5MPa左右。由计算数据可以判断，支墩距管沟位置增加时，管道最大应力值几乎一致，当距离增加到一定值时，管道与管沟边缘接触，土壤给管道一支撑作用，管顶最大应力减小，在距管沟12 m工况下，管顶最大应力为30 8.4MPa,因此可以认为支墩位置对管道完全沉管时管顶最大拉应力影响不大。不同支墩位置管顶轴向应力详细计算结果见表7。表7 D1219mm18.4mm管道不同支墩位置双侧单向、双侧背向管顶最大拉应力MPa支墩距管沟支墩距管沟支墩距管沟支墩距管沟下沟方式0m双侧单向320.6双侧背向320.53.4有无支墩对管道应力的影响为探究大口径高钢级长输管道在有无支墩时双侧单向与背向完全沉管施工管顶最大应力状态与悬空长度分布规律，以不同支墩位置计算结果为基础，选取D1219mm18.4mmX80管道,以密砂、下沟深度3m、有无支墩施工工况为例，进行数值仿真计算。石彤等：大口径X80管道沉管应力响应规律与临界深度分析道纵向位移曲线如图9 所示。双侧单向与双侧背向完全支墩距管支墩距管沟4m沟8 m沟12 m95.595.59696一距管沟0 m0.86距管沟4m350距管沟8 m325距管沟12 m3002752502251754m8 m320.5320.5320.5320.437D1219mm18.4mm管道有无支墩时3m沉降管沉管时，有支墩的管道悬空长度要大于无支墩的。由计算结果可知，相同管道尺寸条件下，管道完全沉管时，有96.5支墩情况下需要加长管沟开挖长度。有无支墩时管道96.5悬空长度计算结果见表8。沟顶0-1-2200225管道坐标/m距管沟0 m距管沟4m距管沟8 m距管沟12 m管道坐标/m12 m308.4308.3无支墩有支墩管沟边缘沟底-30(a)双侧单向管道纵向位移曲线1沟顶0-1-2管沟边缘-301002003004005500600700管道坐标/m(b)双侧背向管道纵向位移曲线图9Dl219mm18.4mm管道不同深度有无支墩双侧单向和双侧背向沉管下沟悬空长度计算曲线表8D1219mm18.4mm管道不同深度有无支墩双侧单向、双侧背向沉管下沟悬空长度D1 219 mm 18.4 mm D1 219 mm 18.4 mm沟深/m下沟方式双侧单向3.0双侧背向双侧单向4.0双侧背向双侧单向5.0双侧背向D1219mm18.4mm管道有无支墩时3m沉降管顶轴向应力分布曲线如图10 所示。有支墩工况下管道应力呈波动分布，双侧单向与双侧背向完全沉管时，有支墩的管顶最大拉应力明显要大于无支墩的。由计算结果可知，相同管道尺寸条件下，管道完全沉管时，有支墩情况下管顶最大应力更大。有无支墩时管顶应力详细计算结果见表9。100管道坐标/m无支墩一有支墩管沟边缘沟底管道（无支墩）悬空长度/m92.793.599.0100.5104.5106.5200300沟顶管道（有支墩）悬空长度/m95.596.0102.0102.0107.0107.540038五油管材与仪器2023年10 月500-0.80,4003002020010000-5-10010-20002550 75100125150管道坐标/m20A5-1002000255075100125150管道坐标/m图10 D1219mm18.4mm管道不同深度有无支墩双侧单向和双侧背向沉管下沟管顶轴向应力分布曲线表9D121918.4mm管道不同深度有无支墩双侧单向和双侧背向沉管下沟管顶最大拉应力D1 219 mm 18.4 mm D1 219 mm 18.4 mm沟深/m下沟方式管顶最大应力/MPa双侧单向301.33.0双侧背向双侧单向4.0双侧背向双侧单向5.0双侧背向3.5管沟深度对管道应力的影响基于上述不同土壤、不同壁厚、支墩位置、有无支墩影响因素分析，探究大口径高钢级长输管道在双侧单向与背向不同管沟深度完全沉管施工时管道悬空长度与管顶最大应力状态分布规律，得到管道完全沉管施工临界沟深。选取D1219mm18.4mmX80管道，并以密砂、不同下沟深度施工工况为例，进行数值仿真计算。D1219mm18.4mm管道不同沟深沉降时管道纵向位移分布曲线如图11所示。双侧单向与双侧背向完全沉管时，管道悬空长度随着沟深增加而增大。对于D1219mm18.4mm管道，下沟深度每增加0.5 m，管道的悬空长度增加值在4m以下,且随着沟深增加,增幅越来越小。不同下沟深度工况管道悬空长度详细计算结果见表10。无支墩有支墩0100管道坐标/m(a)双侧单向管顶轴向应力分布曲线500400管沟迈缘3002001000100200,300400500600700管道坐标/m(b)双侧背向管顶轴向应力分布曲线10管顶-1管道(无支墩)管道(有支墩)管顶最大应力/MPa320.6302.2320.5342.0363.7342.5362.5382.4401.2382.2400.0330320310300290280270200300管沟迈缘无支墩有支墩-2管沟边缘-4-50100200300 400管道坐标/m(a)双侧单向管道纵向位移曲线沉降3.0 m管顶管顶沉降3.5 m0沉降4.0 m沉降4.5 m-1沉降5.0 m0.4-20.3管底0.2-30.10.0-4管沟边缘50100200300400500600700管道坐标/m(b)双侧背向管道纵向位移曲线图11D1219mm18.4mm管道不同深度双侧单向和双侧背向沉管下沟悬空长度计算曲线由D1219mm18.4mm管道不同沉降时管顶轴向应力分布曲线如图12 所示，同一直径、不同沉管深度管道双侧单向、双侧背向完全沉管时管顶轴向应力随着管沟深度增加而增大。可以看出，下沟深度是影响管道安全的最重要因素之一，针对不同尺寸管道的最大下沟深度需经计算得到。对于D1219mm18.4mmX80管道，由不同工况管顶应力详细计算结果见表11。400-0.80,管沟边缘180330320310300290470480490500510520530管道坐标/m沉降3.0 m沉降3.5 m沉降4.0 m沉降4.5 m一沉降5.0 m管底0.00-0.25-0.50150175200225管道坐标/m-0.1-0.2475500525550575管道坐标/m190管道坐标/m2002102202023年第9 卷第5 期表10 D1219mm18.4mm管道不同深度双侧单向和双侧背向沉管下沟悬空长度D1 219 mm 18.4 mm D1 219 mm18.4 mm沟深/m下沟方式双侧单向3.0双侧背向双侧单向3.5双侧背向双侧单向4.0双侧背向双侧单向4.5双侧背向双侧单向5.0双侧背向对模型进一步计算临界沟深与悬空长度可知,无支墩情况下，双侧单向沉管沟深达到6.7 4m、悬空长度达113.0m，管道最大应力达到临界许用应力；双侧背向沉管沟深达到6.8 0 m、悬空长度达115.5 m，管道最大应力达到临界许用应力。有支墩情况下，双侧单向沉管沟深达到6.2 5 m、悬空长度达113.0 m，管道所受最大应力达到临界许用应力；双侧背向沉管沟深达到6.31m时,悬空长度达到113.5 m，管道所受最大应力达到临界许用应力。50004003002001000100管沟边缘-20050100200300400管道坐标/m(a)双侧单向管顶轴向应力分布曲线500-0.80,400375300沉降3.0 m沉降3.5 m350200沉降4.0 m100沉降4.5 m盒32 5八0-100管沟边缘-2000100200300400500600700管道坐标/m(b)双侧背向管顶轴向应力分布曲线图12D121918.4mm管道不同深度双侧单向和双侧背向沉管下沟管顶应力曲线石彤等：大口径X80管道沉管应力响应规律与临界深度分析沟深/m下沟方式管道（无支墩）管道(有支墩)悬空长度/m悬空长度/m92.795.593.596.095.599.097.099.099.0102.0100.5102.0102.0104.5103.5105.0104.5107.0106.5107.50.80.沉降3.0 m沉降3.5 m沉降4.0 m沉降4.5 m沉降5.0 m管沟边缘39表11D121918.4mm管道不同深度双侧单向、双侧背向沉管下沟管顶最大拉应力D1 219 mm 18.4 mm D1 219 mm 18.4 mm管道（无支墩）管顶最大应力/MPa双侧单向301.33.0双侧背向双侧单向3.5双侧背向双侧单向4.0双侧背向双侧单向4.5双侧背向双侧单向5.0双侧背向4 结 论针对大口径X80直管段管道沉管下沟敷设问题，本文以D1219mmX80管道为例，建立了直管段管道沉管下沟敷设参数化数值计算模型,计算了工程中4类影响因素作用下的管道管顶最大拉应力变化及管道纵向位移分布，可以为D1219mmX80管道沉管施工提供指导。1)不同土壤类型条件下，双侧单向与双侧背向完全沉管时，管顶应力和管道悬空长度随土壤弹性模量增大而有小幅增加,并趋于稳定。说明在其他条件不变的情况下，土壤类型对沉管施工中管道悬空长度有一定影响，但这种影响较小。2不同管道壁厚条件下，管道应力和悬空长度的变化趋势相同，同一直径、不同壁厚管道双侧单向、双侧背375向完全沉管时悬空长度与管顶最大拉应力基本相当，同350一管径、不同壁厚条件下，管道完全沉管时管沟开挖长325度一致便满足施工要求。3003)不同支墩位置条件下，当支墩位置发生改变时，190200管道坐标/m480490.500.510520530管道坐标/m管道（有支墩）管顶最大应力/MPa320.6302.2320.5319.8342.2320.7342.1342.0363.7342.5362.5362.8382.3362.8381.8382.4401.2382.2400.0210220管顶最大应力位置随支墩位置移动，管道悬空长度与管顶应力基本一致，认为支墩位置对管道施工影响较小。4)无支墩和有支墩情况下，悬空长度和管顶应力的变化趋势相同，无支墩情况下呈连续性变化，有支墩情况下呈阶梯状或波动状变化；相同管沟深度沉管，管道中最大应力在有支墩