Patran与Abaqus技术在海底管道屈曲分析中的应用研究

针对海底管道屈曲稳定性问题,应用Patran二次开发技术,提出海底管道屈曲分析的参数化自动建模方法,其中核心问题为在Patran中直接建立楔形体单元组成止屈器几何形状以及管道表面skin单元的创建等,生成Abaqus可识别的数据文件。通过算例分析,将计算结果与全尺寸试验和DNV-OS-F101规范进行对比,验证了方法的可靠性。     

双拱初始缺陷海底管线水平向整体屈曲数值模拟分析

为了研究具有双拱反对称初始缺陷海底管线的整体屈曲特性,采用模态分析法将最可能出现的缺陷形态引入数值分析模型中。针对管线在高温高压作用下发生整体屈曲的动态变形特征,运用显式动力数值模拟方法研究了管线整体屈曲过程中水平向变形与轴向变形随温度和内压的变化规律,建立了在整体屈曲过程中屈曲管段与滑动管段轴力的变化过程与初始缺陷形态的关系。将数值模拟结果同经典解析解和室内模型实验结果进行对比,验证了本方法的可靠性。工程算例的分析结果表明,管线整体屈曲的发生是一个由低阶向高阶发展的过程,具有双拱缺陷的管线首先发生二阶模态的整体屈曲,而后过渡到四阶模态;管线整体屈曲的变形包括屈曲段的水平向变形和滑动管段的轴向缩进,其中水平变形释放了管壁内的轴力,轴力的释放量随初始缺陷尖锐程度的降低而增大;轴向缩进变形由于受到地基土的摩阻力使滑动管段内的轴力发生累积,轴力的累积量随初始缺陷的尖锐程度的降低而增加。以上研究成果对指导实际工程具有现实意义。     

椭圆度-凹坑双缺陷海底管道局部屈曲特性研究

海底管道在制造、埋设以及使用过程中极易产生椭圆度-凹坑双缺陷,双缺陷影响管道局部屈曲,对含椭圆度-凹坑双缺陷海底管道的局部屈曲特性研究十分必要。现行规范中采用等效椭圆度对含椭圆度-凹坑海底管道进行评估,该方法无法准确评估不同缺陷形式的屈曲特性。采用形状系数对含椭圆度-凹坑双缺陷的海底管道进行评估,运用有限元软件ABAQUS进行数值模拟,并进行试验验证。在此基础上对含有不同凹坑深度、不同椭圆度的海底管道进行局部屈曲的数值模拟,计算不同形状椭圆度、含有不同凹坑深度海底管道的形状系数,对其进行敏感性分析。计算结果表明:形状系数对海底管道椭圆度、凹坑深度、径厚比敏感性较强;对凹坑宽度敏感性较弱。     

Numerical study on lateral buckling of pipelines with imperfection and sleeper

Lateral buckling is an important issue in unburied high-temperature and high-pressure (HT/HP) subsea pipelines systems. The imperfection–sleeper method is one of the most well-known methods used to control lateral buckling of HT/HP pipelines. Pipelines–sleeper–seabed numerical models are established and verified to analyze the buckling behavior of pipelines using the imperfection–sleeper method. The influence of six main factors on lateral buckling behavior is investigated in details based on the numerical results. Equations of buckling displacement (buckling displacement is defined by the final displacement of the middle point of the pipelines), critical buckling force, and buckling stress (Mises stress) are proposed using the gene expression programming technique. These equations show good accuracy and can be used to assist in the design of sleepers and assess the compressive and stress levels of pipelines.     

管道热屈曲动力过程数值模拟

运输高温高压油气的海底管道会发生整体热屈曲现象。管道热屈曲过程中可能会产生平衡状态的跃迁(snapthrough),且这样的跃迁过程必然会伴随着动力响应。管道热屈曲动力过程中侧向弹出的速度以及轴向缩进的速度对管土相互作用参数的取值有很大影响,然而关于管道热屈曲动力过程的研究却很少。本文给出了数值模拟过程中管道系统阻尼值和升温速率的确定方法,研究了管道初始几何缺陷以及海床参数对管道热屈曲动力过程的影响。     

Model test based soil spring model and application in pipeline thermal buckling analysis

The buckling of submarine pipelines may occur due to the action of axial soil frictional force caused by relative movement of soil and pipeline, which is induced by the thermal and internal pressure. The likelihood of occurrence of this buckling phenomenon is largely determined by soil resistance. A series of large-scale model tests were carried out to facilitate the establishment of substantial data base for a variety of burial pipeline relationships. Based on the test data, nonlinear soil spring can be adopted to simulate the soil behavior during the pipeline movement. For uplift resistance, an ideal elasticity plasticity model is recommended in the case of H/D (depth-to-diameter ratio)>5 and an elasticity softened model is recommended in the case of H/D≤5. The soil resistance along the pipeline axial direction can be simulated by an ideal elasticity plasticity model. The numerical analyzing results show that the capacity of pipeline against thermal buckling decreases with its initial imperfection enlargement and increases with the burial depth enhancement.     

海底管线水平向整体屈曲及失效判断方法研究

海底管线是海洋石油的重要输运手段。为满足输送工艺的需要,正常工作条件下管线往往被施加较高的温度和压强,高温高压使管线内产生附加应力,当附加应力大于土体对管线的约束力时,管线就会发生整体屈曲。过度的水平向整体屈曲会导致截面产生较大的弯曲应力和压缩应变,对管线系统的安全运行造成威胁,因此需要对发生水平向整体屈曲后的管线进行验算。采用解析解法、规范法和有限元法对管线的整体屈曲进行分析,提出了应用临界屈曲荷载值域空间和值域下限来判断不同缺陷大小下管线是否发生水平向整体屈曲的方法。结合工程实例,分别采用内力控制标准和位移控制标准对管线水平向整体屈曲后是否失效进行了验算。研究指出,相较于位移标准,内力控制标准更为严格。     

海底埋设高温管道隆起屈曲数值模拟研究

高温是引发海底管道整体屈曲失效的主要因素,而海床上存在的局部隆起使得高温埋设管道更加容易发生隆起屈曲。这里重点研究海底埋设高温管道发生隆起屈曲的临界温度载荷及其影响因素,提出了一种简化的数值模拟分析模型,同已有的相关实验结果比较表明,本方法可以较好地近似计算高温管道的隆起屈曲。基于本方法开展的参数分析,得到了管道覆土高度、混凝土配重层厚度、海床不平整对海底高温管道发生隆起屈曲的影响趋势。     

海底管线整体屈曲过程中塑性变形研究

海底管线是海洋油气工程中主要的输送手段。在工作状态下,受高温高压的影响,深海管线可能会发生水平向整体屈曲。随海洋油气作业水深的增大,施加于管线的温度和压强也逐渐增加,导致管线产生较大的屈曲位移和截面应力,使得截面产生塑性应变。本文采用数值模拟方法,对海底管线整体屈曲过程中塑性区的分布及其与整体屈曲影响因素的关联性、塑性变形对水平向变形的影响和塑性变形造成的截面椭圆度的变化规律进行分析,研究塑性变形对整体屈曲过程影响的规律。     

后挖沟深度对深海海底管道屈曲影响数值分析

鉴于海底管道的服役水深越来越深,主要采用犁式挖沟机对预铺设于海床之上的海底管道采取后挖沟的方式将海底管道埋设于海床之下,以保护其免受不必要的损伤。针对后挖沟深度H是海底挖沟机的重要设计参数,也是影响管道悬跨的重要因素的问题,对SMD(UK)犁式挖沟机展开参数优化,确保作业过程中悬跨段管道在外部静水压力作用下,海底管道不会发生屈曲破坏。采用ABAQUS软件,分别建立了作业前和作业中两种工况下的悬跨模型,分析机械手对接触部分管道的损伤,结果显示,作业中的机械手对悬跨管道的损伤更大;同时,建立了作业中不同管径下,后挖沟深度对管道损伤的安全裕量关系曲线。进一步,结合作业中不同挖沟深度下的管跨段屈曲数值模型,对处于外部静水压力作用下的悬跨管的屈曲失效展开分析,结果显示,随着后挖沟深度的加大,不同管径下的悬跨段管道局部出现塑性压溃的临界压力值不断降低;管道外径的增大,降低了同一后挖沟深度下发生屈曲失效的压力值。最后,在后挖沟深度与外部静水压力组成的区域内,建立屈曲失效临界关系曲线,并划分出工作区和压溃区,为深海管道后挖沟埋管的施工提供工程参考。     

Numerical analysis of cross-section ovalization in the deep-sea pipeline lateral buckling process

Abstract

The deep-water pipeline is the main means of transportation in offshore oil and gas development engineering. The deep-water pipeline may incur lateral global buckling due to the high temperature and pressure that are applied on the pipeline to ensure the contents’ liquidity. With the increasing operating water depth, a higher temperature and pressure are applied to the pipeline, causing large lateral deformation and a large bending moment. Due to the inhomogeneous distribution of the bending moment on the cross-section, different points on the cross-section will deform differently. This kind of deformation causes the cross-section to turn into an oval ring. The cross-section ovalization caused by global buckling was rarely analyzed in former engineering practice since the load is relatively low. With the increase in operation water depth and operation load, the ovality caused by global buckling is noticeable. This article analyzed cross-section ovalization caused by pipeline lateral global buckling with a numerical simulation method. The pipelines with different initial cross-section shapes were simulated, and the influence of several impact factors, including load, pipeline and soil factors on the ovality of the cross-section, were analyzed. The results show that the initial cross-section shape type has little effect on the pipeline ovalization pattern. The initial ovality of the pipeline with an oval ring cross-section shape has little influence on the residual ovality. Among all the factors analyzed in this paper, the pressure difference is the primary factor that should be considered in a pipeline ovalization check.     

含初始缺陷海底管道侧向屈曲临界力的参数解

未埋设的海底管道在高温高压运行条件下可能发生侧向屈曲,情况严重时影响海底管道的结构安全。侧向屈曲临界力作为判定海底管道发生侧向屈曲的重要依据,主要影响因素有初始缺陷、管土相互作用等。现有关于侧向屈曲临界力的公式并未考虑管土相互作用、缺陷不平直度和管道自身材料特性对侧向屈曲临界力的综合影响。建立含有通用几何初始缺陷海底管道的数值模型,使用Riks算法进行参数分析以研究极限侧向土壤阻力、管道缺陷不平直度和截面几何尺寸对海底管道侧向屈曲的具体影响。基于量纲分析法和多元线性回归,推导出海底管道侧向屈曲临界力关于上述3个影响参数的一般公式,并对该公式进行了检验,结果表明文中推导的公式在参数涵盖的研究范围内有效。     

海底管道整体屈曲的枕木-浮力耦合法研究

不埋海底管道在高温高压作用下,易发生水平向整体屈曲。实际工程中,常通过在管道路由上设置整体屈曲触发装置,实现对水平向整体屈曲的有效控制,其中以枕木法的成功应用最为多见。本文分析了枕木法的主要影响因素并验证了采用枕木法会出现管道屈曲段应力集中的现象,对比了枕木法、分布浮力法和枕木-浮力耦合法对管道整体屈曲变形规律的影响,采用数值模拟方法系统研究了枕木及浮力参数对管道水平向屈曲和后屈曲的影响规律。研究表明,在枕木两侧设置浮力段的人工触发装置可有效触发管道整体屈曲,同时促使管道虚拟锚固点间轴力的释放,降低了管道中屈曲段的应力,相较枕木法,枕木-浮力耦合法可将管道中的最大应力降低23%。     

蛇形铺设管道触发水平向整体屈曲效果研究

海底管道在深海石油开发工程中有着广泛的应用,管道在工作时受到高温高压会触发水平向整体屈曲变形,蛇形铺管法是控制管道水平向整体屈曲变形的有效手段。采用数值模拟方法,对蛇形铺设管道的关键参数进行研究,分析关键参数对临界屈曲荷载和屈曲后截面应力应变状态的影响。对每一组蛇形铺设管道都设置了直线型铺设管道的对照组,对照组采用引入初始挠曲的方法激发水平向整体屈曲,初始挠曲程度与蛇形铺设管道相同。通过与直线型铺设管道的对比,显示了蛇形铺设管道在激发和控制屈曲方面的优势,并对不同土体阻力情况下蛇形铺设管道的实际效果进行了评估。结果表明,蛇形铺设管道的临界屈曲力和屈曲后的弯矩皆远小于直线型管道。蛇形铺设管道的临界屈曲力随圆心角θ增大而减小,随曲率半径R的增大而增大。增大跨度L、曲率半径R和圆心角θ都能有效减小蛇形铺设管道的截面弯矩。水平向土体阻力对蛇形铺设管道影响较大,水平向土体阻力较小时,蛇形铺设管道控制屈曲的效果更为明显。     

Numerical evaluation of a subsea equipment installation method designed to avoid resonant responses

Several methods have been employed to install the different types of subsea equipment that are required for offshore oil and gas production systems. More recently, the increasing development of oil fields at remote, ultra-deepwater scenarios has been requiring heavier and more complex subsea equipment. In such scenarios the usual installation methods have been facing increasing challenges, including resonance effects associated to the dynamic behavior of the cable-equipment system. In this context, after presenting a review of the state-of-the-art of installation methods for subsea equipment, this work presents a methodology for the design and evaluation of an installation method based on the combination of steel wire and polyester ropes. The goal is to obtain a feasible and safe installation method that can avoid resonance effects along all phases of the lowering procedure, without the need of specialized devices such as heave compensators, and using only simpler installation vessels such as standard tugboats instead of more expensive and specialized vessels (such as drilling rigs). The method is evaluated by numerical simulations for case studies considering representative metocean data for Brazilian ultra-deepwater scenarios, following a methodology that considers the random characteristic of actual sea states by employing an irregular spectral representation. The results indicate that the method may indeed be feasible for the deepwater installation of heavier equipment, potentially being more cost-effective especially for remote locations where other more complex methods might become inefficient.     

Non-linear wave-induced transient response of soil around a trenched pipeline

Submarine pipelines are always trenched within a seabed for reducing wave loads and thereby enhancing their stability. Based on Biot's poroelastic theory, a two-dimensional finite element model is developed to investigate non-linear wave-induced responses of soil around a trenched pipeline, which is verified with the flume test results by Sudhan et al. [Sudhan, C.M., Sundar, V., Rao, S.N., 2002. Wave induced forces around buried pipeline. Ocean Engineering, 29, 533–544] and Turcotte et al. [Turcotte, B.R., Liu, P.L.F., Kulhawy, F.H., 1984. Laboratory evaluation of wave tank parameters for wave-sediment interaction. Joseph H. Defree Hydraulic Laboratory Report 84-1, School of Civil and Environmental Engineering, Cornell University]. Non-linear wave-induced transient pore pressure around pipeline at various phases of wave loading is examined firstly. Unlike most previous investigations, in which only a single sediment layer and linear wave loading were concerned, in this study, the influences of the non-linearity of wave loading, the physical properties of backfill materials and the geometry profile of trenches on the excess pore pressures within the soil around pipeline, respectively, were explored, taking into account the in situ conditions of buried pipeline in the shallow ocean zones. Based on the parametric study, it is concluded that the shear modulus and permeability of backfill soils significantly affect the wave-induced excess pore pressures around trenched pipeline, and that the effect of wave non-linearity becomes more pronounced and comparable with that of trench depth, especially at high wave steepness in shallow water.     

Global buckling of pipelines in the vertical plane with a soft seabed

Offshore pipelines are usually buried to avoid damage from fishing activities and to provide thermal insulation. Provided that the buried pipelines are sufficiently confined in the lateral direction by the passive resistance of the trench walls, they may be subject to vertical buckling caused by a rise in temperature. Vertical buckling is usually called upheaval buckling because the heated pipeline is assumed to move upwards conventionally. However, the seabed may be very soft, especially where a pockmark or abyssal ooze appears. Consequently, under thermal compressive force, the pipeline may buckle downward and penetrate into the seabed because the downward soil resistance is small. In this study, we extended an analytical solution for vertical pipeline buckling on a rigid seabed to a soft seabed, and the effects of soil resistance on pipeline stability, buckling mode and amplitude are illustrated and analyzed.     

Feasible numerical method for analysis of offshore pipeline in installation

Feasible numerical method for a structural analysis of a pipeline configuration during the installation process is presented. The method considers the whole pipeline, which is partially suspended and partially laid-on a seabed, as a single continuous segment, and is valid for a complete range of laying angles between 0°–90°, i.e., valid for both S-lay and J-lay configurations. The method accounts for a pipeline–seabed interaction and the pipeline is modeled by means of nonlinear large deformation beam theory. The numerical solution is carried out in an incremental-iterative manner by following the actual pipeline installation process, and thus allowing efficient treatment of pipeline-seabed interaction circumventing the further complexities with contact detection. At each increment, the length of the pipeline is increased and new sequential equilibrium configuration is assessed by direct minimization of a total potential energy approximated as a Riemann sum, which yields algebraic system of nonlinear finite difference equations that is further solved by iterations with Newton-Raphson technique. The simplicity, flexibility and robustness of the proposed method allow to enhance the efficiency of engineering calculations and design. Accounting for a bending stiffness in a suspended part allows analyzing variations in laying angle and lay tension independently. The method convergence is validated and compared with Abaqus. The results are in an excellent agreement. Moreover, the comparison with Abaqus shows that for the selected parameters the assumption that the pipeline is inextensible and unshearable is very reasonable. Representative parametric study is conducted to demonstrate the feasibility of the method. Parametric study considers the effects of laying angle (0°–90°), lay tension, laying water depth (up to 3000 m) and seabed stiffness.     

深水S型铺管中残余变形对屈曲压力的影响

随着S型海管铺设逐步走向深水,管道在铺设过程中承受的荷载增加,发生一定程度的塑性变形。本文讨论S型铺设引起的残余塑性变形对管道屈曲承载能力的影响。首先基于壳单元建立"管道-托辊"相互作用的局部有限元模型,分析了管道在铺设过程中的受力状态,获得了管道的残余塑性变形。然后以该残余变形作为管道非线性屈曲分析的初始缺陷,基于改进的RIKS方法计算了管道的临界屈曲压力。研究结果表明,铺设残余塑性变形在一定程度上削弱管道的承载能力,在深水铺设中应予以考虑。     

Modeling for the optimization of layout scenarios of cluster manifolds with pipeline end manifolds

The pipeline end manifold (PLEM) is an important subsea facility, which can greatly reduce the cost and risk of the development scenarios of deepwater oil and gas fields by declining the number of export pipelines and risers. However, the employment of PLEMs is a multidiscipline task involving substantial financial and technical factors. Due to various uncertainties of influencing factors, the evaluation process may take several months or years by the engineers with rich project experience. Thus, how to develop quantified reference tools using mathematical models to assist engineers in efficiently making their crucial decisions is essential. In this paper, the optimization of the layout scenarios of cluster manifolds with PLEMs is discussed, where a proposed mathematical model and its dedicated algorithm are illustrated. The optimal solution at the lowest cost can be obtained through in-house routine in MATLAB, including the optimal layout scenario, the number and locations of PLEMs, and the connection relations. Besides, the numerical simulations are performed to demonstrate the validity of the proposed mathematical model and its algorithm. The results show that this optimization layout problem in engineering can be described accurately by the presented mathematical model and the convergence rate of the given algorithm is efficient.