含复合材料结构的非黏结柔性立管弯曲特性
Bending Properties of Unbonded Flexible Risers with Composite Materials
通讯作者: 薛鸿祥,男,教授,博士生导师,电话(Tel.):021-34204470;E-mail:hongxiangxue@sjtu.edu.cn.
责任编辑: 陈晓燕
收稿日期: 2021-04-30
基金资助: |
|
Received: 2021-04-30
作者简介 About authors
刘庆升(1990-),男,黑龙江省哈尔滨市人,博士生,从事非黏结柔性立管研究.
非黏结柔性立管常用于将海底油气资源输送到平台,有时因设计需要会包含复合材料圆柱壳层结构.以一类8层的非黏结柔性立管为研究对象,分析了含有复合材料圆柱壳层结构的非黏结柔性立管的弯曲特性.提出了一种作用于非黏结柔性立管的轴对称载荷和弯曲载荷联合作用的理论模型,同时建立了计及详细几何特性的非黏结柔性立管数值模型进行了验证.分析了两种典型复合材料及其纤维体积分数对非黏结柔性立管弯曲特性的影响.理论与数值结果吻合较好,计算结果表明:非黏结柔性立管的弯曲特性受复合材料圆柱壳层轴向弹性模量影响较大,轴向弹性模量较大的复合材料会较大地增强完全滑移阶段的非黏结柔性立管的弯曲刚度.非黏结柔性立管的轴向抗拉刚度和完全滑移阶段的弯曲刚度都随着复合材料圆柱壳层的纤维体积分数增大,接近线性增强.
关键词:
Unbonded flexible risers are widely applied to transport oil and gas resources from seabed to platform, and composite cylindrical layers are sometimes contained due to design requirements. Based on an 8-layer unbonded flexible riser model, the bending properties of the unbonded flexible risers with composite cylindrical layers are studied. A theoretical method of combining axisymmetric and bending loads acting on unbonded flexible risers is proposed, and a numerical method with detailed geometric characteristics taken into account is established for verification. The influence of two typical composite materials as well as the fibre volume fraction on the bending properties of the unbonded flexible risers is analyzed. The theoretical and the numerical results are in good agreement, which shows that the bending properties of the unbonded flexible risers are greatly affected by the axial Young’s modulus of the composite cylindrical layers, and composite material with larger axial Young’s modulus would greatly enhance the bending stiffness of the unbonded flexible risers, expecially in the full-slipping stage. In addition, the axial tensile stiffness and bending stiffness of full-slipping stage are proportional to the fibre volume fraction of the composite cylindrical layers.
Keywords:
本文引用格式
刘庆升, 薛鸿祥, 袁昱超, 唐文勇.
LIU Qingsheng, XUE Hongxiang, YUAN Yuchao, TANG Wenyong.
非黏结柔性立管被广泛应用于海洋油气业,常用于从海底到平台的油气资源传输,有时会根据设计的需要增加耐腐蚀、隔热层或是防止发生径向失效的复合材料圆柱壳层结构[1],例如,Epsztein等[2]给出了含有耐H2S气体腐蚀的复合材料层的非黏结柔性立管结构,Malta等[3]给出含有以凯夫拉49复合材料构成的强结构带层,用于防止非黏结柔性立管发生轴向压力作用下的径向失效,但一般学者都直接将复合材料简化为各向同性材料或是等效为忽略复合材料圆柱壳层厚度变化的正交各向异性结构进行研究.非黏结柔性立管的中的螺旋层结构在弯矩载荷作用下可以克服层间的摩擦力并在层间自由滑动,导致其弯曲刚度具有高度非线性[4⇓⇓-7],一直以来是研究的热点和难点所在.基于一些严格的假设,一些学者给出了非黏结柔性立管的非线性弯曲模型.Féret等[8]较早地给出螺旋条带弯矩-曲率迟滞曲线,但未考虑部分滑移阶段的弯矩-曲率关系.文献[5-6,9]基于功能原理给出了螺旋条带未滑移、部分滑移和完全滑移3个阶段的弯矩-曲率表达式.此后,Dong等[7]在此基础上进一步给出了部分滑移阶段的弯曲刚度的显式表达式,并合理考虑了螺旋条带局部弯曲和扭转的影响.Ramos等[10-11]也给出了螺旋条带未滑移和完全滑移阶段的弯曲刚度计算方法,可用于快速计算非黏结柔性立管未滑移和完全滑移阶段的弯曲刚度.与理论方法相比,数值和试验方法的相关研究相对较少.由于非黏结柔性立管造价较高,一般只有立管厂商会开展试验研究[12-13],但通常不会给出详细的立管参数.Witz[14]在1996年开展了关于非黏结柔性立管截面力学性能的试验,不仅给出了详细的立管参数和试验结果,而且给出了试验的开展过程和边界条件,此后的理论和数值方法研究大都基于Witz的试验进行验证.数值方法经历了从等效简化数值模型到考虑非黏结柔性立管详细几何特性的数值模型发展[15],文献[4,16]都建立了考虑骨架层和抗压铠装层实际几何形状的数值模型,研究了非黏结柔性立管的弯曲特性.
以往研究大多只关注螺旋层结构的弯曲非线性模型,并将作用于非黏结柔性立管的轴对称载荷作为计算非黏结柔性立管层间接触压力的基础,而忽略了轴对称载荷对弯曲特性的影响,在立管的弯曲刚度较小时,作用于立管的轴对称载荷可能会对立管的弯曲特性产生较大影响.本文以一类8层的非黏结柔性立管为研究对象,考虑作用于非黏结柔性立管的轴对称载荷和弯曲载荷的联合作用,提出了一种轴对称载荷对非黏结柔性立管弯曲特性影响的理论模型.考虑了复合材料圆柱壳层结构对非黏结柔性立管弯曲特性的影响,通过将原防摩擦层结构替换为复合材料圆柱壳层结构,采用功能原理推导了轴对称载荷作用下计及复合材料层结构厚度变化的理论模型,在此基础上求解了非黏结柔性立管的层间相互接触压力,并推导出弯曲载荷作用下的理论模型.建立了含有详细几何特性的非黏结柔性立管数值模型,并进行了验证.基于本文提出的理论和数值方法,计算了含有复合材料圆柱壳层结构的弯曲特性,并分析了复合材料纤维体积分数对非黏结柔性立管弯曲特性的影响.
1 理论方法
1.1 复合材料层弹性常数计算
1.2 复合材料圆柱壳层在轴对称载荷作用下的理论模型
非黏结柔性立管轴对称载荷作用下的理论模型是弯曲载荷作用下理论模型的基础,需要基于轴对称载荷作用下的理论模型求解层间接触力,并基于该结果求解螺旋层结构的临界曲率,判断螺旋层是否发生滑移[19].本文对复合材料圆柱壳层模型进行了一定简化,忽略了复合材料的分层以及制造过程中的损伤,将宏观的复合材料圆柱壳层结构视为一个整体,在此基础上基于受力和变形关系推导相应的理论模型.
假设作用于圆柱壳层结构的外力载荷为轴向力F、绕中心轴转矩T、内外压力
图1
式中:
1.3 复合材料圆柱壳层在弯曲载荷作用下的理论模型
先假设作用于圆柱壳层的轴对称载荷和弯曲载荷可以进行解耦,只考虑在单一方向的弯矩作用,并忽略圆柱壳层其他方向的曲率和挠率变化,相应受到的弯矩作用Mc和曲率κ的关系可以简化为
式中:I为惯性矩.
1.4 轴对称载荷与弯曲载荷联合作用下的理论模型
图2
式中: r为曲率半径;φ为绕圆柱壳层单元周向的角度;β为圆柱壳层曲率中心到该位置的夹角.
利用式(4)对角度φ和β进行积分,可以得到微元上的压力在图2坐标轴各方向的分量:
式中:θ为微元中心到曲率中心的角度;圆柱壳层结构外径横截面积Ao=π
忽略轴对称载荷引起的圆柱壳层结构的轴向伸长,假设圆柱壳层结构在轴对称载荷和弯曲载荷联合作用后各处曲率保持一致,如图3所示.图中,RP1和RP2分别代表圆柱壳层两端的几何中心端点; O为曲率中心.并考虑弯曲后的圆弧角度δ与曲率的关系为δ=Lκ,由轴对称载荷引起的对弯曲特性的影响
图3
对于单位长度的圆柱壳层结构,定义变量a为
求解曲率0到1范围内的变量a并进行线性拟合,得到如图4所示的曲线,忽略线性拟合后的曲线截距得到简化后的变量a的表达式为
图4
相应的外压载荷引起的弯矩可以表示为
同理可以推导内压和轴向力载荷的影响.
2 数值方法
通过Abaqus软件建立了含有复合材料圆柱壳层结构的8层非黏结柔性立管数值模型,各层的剖面如图5(a)所示,完整的有限元模型示意如图5(b)所示.根据立管厂商的建议,为了减少边界效应的影响,非黏结柔性立管模型需要有足够的长度,通常为两倍的抗拉铠装层螺距的长度[22],本文数值模型长度为1 m,其中,骨架层和抗压铠装层结构计及了实际的截面几何形状(见图5(b)),非黏结柔性立管有限元模型整体包含超过30万个单元.在有限元模型两端建立两个几何参考点RP1和RP2,并将所有层截面的端点与之耦合,层间摩擦因数取为0.1[4].考虑到数值模型中存在大量的非线性接触的情况,以及尽可能真实地模拟非黏结柔性立管在复杂载荷作用下的结构响应,所有层结构都采用非协调体单元(C3D8I)进行模拟,并采用显式算法进行求解,采用通用接触模拟立管内部结构的复杂非线性接触.采用放大计算时间步的方法提高数值方法的计算效率,本文的最小时间步长取为5×10-7s,经验证,该时间步长可以在确保计算精度和稳定性的同时提高计算效率.同时,在复合材料圆柱壳层中定义局部柱坐标系赋予材料属性.采用准静态的方式对载荷进行加载,并释放如图5(b)中绕x轴旋转的自由度,以及顶端RP1参考点的轴向(z方向)自由度.非黏结柔性立管的几何和材料特性如表1所示,其中几何特性来自文献[14].
图5
图5
非黏结柔性立管剖面和有限元模型示意图
Fig.5
Schematic of longitudinal sectional and finite element model of unbonded flexible risers
Tab.1 Geometric and material properties of unbonded flexible risers
层数 | 名称 | 条带数 | 内径/mm | 铺设角/(°) | 材料 |
---|---|---|---|---|---|
1 | 骨架层 | 1 | 31.60 | -87.5 | AISI304 |
2 | 内部护套层 | - | 35.10 | - | 尼龙12 |
3 | 抗压铠装层 | 2 | 40.05 | -85.5 | FI-15 |
4 | 加强层 | - | 46.25 | - | 复合材料 |
5 | 抗拉铠装层 | 40 | 47.75 | -35.0 | FI-41 |
6 | 加强层 | - | 50.75 | - | 复合材料 |
7 | 抗拉铠装层 | 44 | 52.25 | 35.0 | FI-41 |
8 | 外部护套层 | - | 55.25 | - | - |
3 计算结果及分析
3.1 数值和理论方法验证
非黏结柔性立管在弯曲载荷作用下的理论模型以轴对称载荷作用下的理论模型为基础进行研究.首先通过轴对称载荷作用下的理论模型求解非黏结柔性立管不同层之间的接触压力,在此基础上求解螺旋层结构的弯曲临界曲率,之后加载弯矩载荷的作用,同时结合本文提出的轴对称载荷和弯曲载荷联合作用下的理论模型,求解非黏结柔性立管在未滑移、部分滑移和完全滑移阶段的非黏结柔性立管的弯曲特性.对于数值方法,载荷的加载方式与理论方法保持一致,同样地首先施加轴对称载荷,之后再在两端几何参考点RP1和RP2施加弯矩为M的作用.采用显式算法进行求解时需要考虑惯性效应对结果的影响[23],因此本文采用准静态的方法对载荷进行加载,并模拟足够的加载时长以避免惯性效应可能对结果造成的影响[4],并始终控制动能与内能的比率在5%以内.
通过本文提出的理论和数值方法计算Witz[14]试验中的工况,相应的计算结果如图6、7所示,得到相应的截面力学性能如表2所示,表中Ea为抗拉刚度;Ei为完全滑移阶段的弯曲刚度.从图7和表2中可以看出:轴对称载荷作用下的数值和理论方法吻合较好,且轴向延伸率随着轴向张力的增大基本呈线性变化,而试验结果存在一定的迟滞效应,将试验结果进行拟合得到的轴向抗拉刚度与本文的结果吻合较好,证明了理论和数值方法的有效性.在弯曲载荷作用下,分别通过理论方法计算了未考虑外压载荷和考虑外压载荷的情况,其中,考虑了外压载荷的完全滑移阶段的弯曲刚度与数值结果吻合较好,而未考虑外压载荷的弯曲刚度与试验结果吻合较好,侧面验证了本文数值方法的有效性.
图6
图6
非黏结柔性立管轴向延伸率-张力曲线
Fig.6
Axial elongation-tension curves of unbonded flexible risers
图7
图7
非黏结柔性立管弯矩-曲率曲线
Fig.7
Bending moment-curvature curves of unbonded flexible risers
表2 不同方法计算结果
Tab.2
方法 | 力学性能 | |
---|---|---|
Ea/MN | Ei/(kN·m2) | |
理论 | 105.88 | 1.19, 1.42 |
数值 | 99.19 | 1.41 |
试验 | 91.19 | 1.19, 1.04 |
需要指出的是,相比于数值方法,本文所提理论方法具有极高的计算效率,且可以得到与数值方法近似的计算结果.本文的数值方法在某具有48核的工作站上开展,弯曲载荷工况的载荷加载时间为 1.8 s,单一工况计算需要约72 h.本文的理论计算方法通过MATLAB软件进行编程,计算时间约为 15 s.
3.2 含复合材料圆柱壳层立管弯曲特性
表3 复合材料的弹性常数
Tab.3
材料 | 弹性模量/GPa | μ12 | μ13 | μ23 | G12/GPa | ||
---|---|---|---|---|---|---|---|
E1 | E2 | E3 | |||||
1 | 9.40 | 4.02 | 4.01 | 0.36 | 0.36 | 0.49 | 4.08 |
2 | 16.0 | 1.07 | 3.86 | 0.44 | 0.44 | 0.45 | 0.56 |
同样对于含有复合材料圆柱壳层的非黏结柔性立管模型首先施加0.2 MPa的外压载荷作用,之后在模型的两端几何参考点施加弯矩,求解弯曲特性.以含有复合材料1的非黏结柔性立管模型为例,分别给出模型在仅加载外压载荷以及加载外压和弯矩耦合载荷作用时外部抗拉铠装层结构在立管轴向的应力(σ)云图,如图8所示.图8(a)为立管在仅施加外压载荷时的应力云图,可以看出应力在立管轴向方向并未完全均匀分布,但数值计算结果比较稳定,外压载荷和立管的轴向位移基本呈线性变化,在外压加载结束后,非黏结柔性立管的轴向延伸率仅为 -0.0105%,对计算曲率影响较小,可以在计算时忽略.图8(b)给出了在外压载荷和弯矩载荷联合作用下的应力云图,可以看出螺旋条带有着明显的弯曲变形,边界处有明显的应力集中现象,此外,与理论方法的假设一致,单一螺旋条带的应力沿着其轴线方向呈现周期性变化.分别通过本文提出的理论和数值方法求解含有复合材料圆柱壳层的非黏结柔性立管的弯矩-曲率曲线,相应的计算结果如图9所示,可以看出数值方法计算得到的未滑移和部分滑移阶段不明显,但在完全滑移阶段理论结果和数值结果吻合较好,相应的完全滑移阶段的弯曲刚度计算结果如表4所示.对比原非黏结柔性立管模型,强度较高的复合材料会对弯曲刚度增强较大,相应的弯曲变形减小.同时,给出各层结构对非黏结柔性立管整体的弯曲刚度贡献,如表5所示.可以看出,复合材料圆柱壳层结构贡献了大部分的非黏结柔性立管模型整体的弯曲刚度,在承受弯矩载荷作用时,含有轴向强度较高的复合材料圆柱壳层结构的非黏结柔性立管模型将不再容易发生大挠度弯曲变形.
图8
图8
不同受力情况下外部抗拉铠装层沿立管轴向应力云图
Fig.8
Stress contours of outer tensile armor layers in the axial direction of risers
图9
图9
含复合材料圆柱壳层非黏结柔性立管弯矩-曲率曲线
Fig.9
Bending moment-curvature curves of unbonded flexible risers with composite cylindrical layers
Tab.4 Calculation results of bending stiffness of unbonded flexible risers with composite cylindrical layers
复合材料 | Ei/(kN·m2) | ||
---|---|---|---|
理论方法 不考虑外压 | 理论方法 考虑外压 | 数值方法 | |
玻纤/环氧树脂 | 11.51 | 11.74 | 11.52 |
钢丝/高密度聚乙烯 | 18.92 | 19.15 | 18.57 |
表5 完全滑移阶段各层结构对非黏结柔性立管模型的弯曲刚度贡献
Tab.5
层结构及轴对称载荷影响 | Ei/(kN·m2) | ||
---|---|---|---|
含防摩擦层尼龙-11模型 | 含复合材料玻纤/环氧树脂模型 | 含复合材料钢丝/高密度聚乙烯模型 | |
第1, 2, 3, 5, 7, 8层结构的综合作用 | 0.85 | — | — |
第4层 | 0.15 | 4.61 | 7.81 |
第6层 | 0.19 | 6.05 | 10.26 |
轴对称载荷作用 | 0.23 (0.2 MPa外压载荷作用) | — | — |
总计 | 1.42 | 11.74 | 19.15 |
3.3 纤维体积分数的影响分析
表6 材料特性
Tab.6
材料类型 | 弹性模量/GPa | 泊松比 |
---|---|---|
普通碳纤维 | 221 | 0.30 |
钢丝 | 210 | 0.26 |
凯夫拉49 | 131 | 0.36 |
玻璃纤维 | 88 | 0.36 |
环氧树脂 | 3.5 | 0.33 |
考虑复合材料圆柱壳层的纤维体积分数从0%至70%变化,通过理论方法计算不同纤维体积分数的复合材料对非黏结柔性立管完全滑移阶段弯曲刚度的影响,结果如图10所示.从图中可以看出,含有复合材料圆柱壳层的非黏结柔性立管完全滑移阶段的弯曲刚度随着纤维体积分数的增大接近线性增长,但与原非黏结柔性立管模型相比,完全滑移阶段的弯曲刚度已经发生了巨大的变化,在复合材料材料强度最大的工况(普通碳纤维材料,纤维体积分数为70%),相应的完全滑移阶段的弯曲刚度为 176.78 kN·m2,远大于原含有防摩擦层非黏结柔性立管模型完全滑移阶段的弯曲刚度.
图10
图10
含有不同纤维体积分数的复合材料立管模型的完全滑移阶段的弯曲刚度
Fig.10
Full-slipping bending stiffness of riser model with different fibre volume fraction composite materials
同时,以碳纤维复合材料为例,考虑含有不同纤维体积分数复合材料圆柱壳层结构的非黏结柔性立管模型的轴向抗拉刚度的影响,给出轴向抗拉刚度与完全滑移阶段的弯曲刚度的关系,如图11所示.从图中可以看出,不同含有纤维体积分数复合材料圆柱壳层非黏结柔性立管的轴向抗拉刚度和完全滑移阶段的弯曲刚度正相关,但完全滑移阶段的弯曲弯曲刚度比轴向抗拉刚度对纤维体积分数更敏感,以碳纤维复合材料为例,当纤维体积分数从0%至70%变化时,轴向抗拉刚度增大了116.29%,而完全滑移阶段的弯曲刚度增大了 3582.92%.
图11
图11
含有碳纤维复合材料的立管模型完全滑移阶段弯曲刚度-轴向抗拉刚度曲线
Fig.11
Full-slipping bending stiffness-axial tensile stiffness of riser model with carbon fibre composites
4 结论
采用复合材料圆柱壳层结构替换了原非黏结柔性立管中的防摩擦层结构,考虑复合材料的材料特性,并基于功能原理推导了复合材料圆柱壳层结构在轴对称载荷作用下的平衡方程,在此基础上推导了含有复合材料圆柱壳层结构在弯曲载荷作用下的理论方程,提出了施加于非黏结柔性立管的轴对称载荷和弯曲载荷联合作用的理论模型,同时建立了相应的含有详细几何特性的数值模型,并通过与试验结果进行对比验证了所提方法的准确性.本文提出的理论模型可以有效地解释轴对称载荷和弯曲载荷的联合作用,并快速地计算得到复杂载荷作用下的非黏结柔性立管的结构响应,准确地计算得到复杂载荷作用下的非黏结柔性立管的弯曲特性.主要结果有:
(1) 外压载荷会增大非黏结柔性立管的弯曲刚度.
(2) 采用强度较高的复合材料替换原防摩擦层结构,非黏结柔性立管的完全滑移阶段的弯曲刚度会有较大提升,并且不再容易发生大挠度弯曲变形.
(3) 复合材料的纤维体积分数对含有复合材料圆柱壳层结构的非黏结柔性立管的轴向抗拉刚度和完全滑移阶段的弯曲刚度都有一定提升,并且随着纤维体积分数增大,两种截面力学性能都基本呈线性变化.
(4) 完全滑移阶段的弯曲刚度对复合材料纤维体积分数的敏感性较高,轴向抗拉刚度对复合材料纤维体积分数的敏感性较低.
参考文献
海洋非粘结型柔性软管的性能结构及其研究热点
[J]. ,
Structure and research focus of marine unbonded flexible pipes
[J]. ,
New anti H2S layer for flexible pipes
[C]// .
Finite element analysis of flexible pipes under compression
[C]// ,
非黏结柔性立管骨架层失效数值计算方法
[J]. ,
Numerical method to predict failure of carcass layer of unbonded flexible risers
[J]. ,
Slip initiation and progression in helical armouring layers of unbonded flexible pipes and its effect on pipe bending behaviour
[J]. ,DOI:10.1243/0309324011514458 URL [本文引用: 3]
An analytical model to predict the bending behavior of unbonded flexible pipes
[J]. ,DOI:10.5957/jsr.2013.57.3.171 URL [本文引用: 2]
Calculation of stresses and slip in structural layers of unbonded flexible pipes
[J]. ,DOI:10.1115/1.3257019 URL [本文引用: 1]
Non-linear bending behaviour of offshore flexible pipes
[C]// .
A consistent analytical model to predict the structural behavior of flexible risers subjected to combined loads
[J]. ,DOI:10.1115/1.1710869 URL [本文引用: 1]
A new analytical expression to estimate the bending stiffness of flexible risers
[C]// .
Experimental estimation of physical properties of a flexible riser
[C]// .
An experimental investigation on the bending behaviour of flexible pipes
[C]// .
A case study in the cross-section analysis of flexible risers
[J]. ,DOI:10.1016/0951-8339(95)00035-6 URL [本文引用: 5]
轴向承载力作用下非粘结柔性立管的侧向失效分析
[J]. ,
Lateral failure analysis of unbonded flexible risers under axial compression force
[J]. ,
Theoretical and numerical analysis of bending behavior of unbonded flexible risers
[J]. ,DOI:10.1016/j.marstruc.2015.10.001 URL [本文引用: 1]
Theoretical and numerical methods to predict the behaviour of unbonded flexible riser with composite armour layers subjected to axial tension
[J]. ,
Review of riser analysis techniques
[J]. ,DOI:10.1016/S0141-1187(82)80002-3 URL [本文引用: 1]
Mathematics of flexible risers including pressure and internal flow effects
[J]. ,DOI:10.1016/0951-8339(92)90025-K URL [本文引用: 1]
/
〈 | 〉 |