上海交通大学学报, 2022, 56(9): 1247-1255 doi: 10.16183/j.cnki.jsjtu.2021.144

船舶海洋与建筑工程

含复合材料结构的非黏结柔性立管弯曲特性

刘庆升1,2, 薛鸿祥,1,2, 袁昱超1,2, 唐文勇1,2

1.上海交通大学 海洋工程国家重点实验室

2.高新船舶与深海开发装备协同创新中心,上海 200240

Bending Properties of Unbonded Flexible Risers with Composite Materials

LIU Qingsheng1,2, XUE Hongxiang,1,2, YUAN Yuchao1,2, TANG Wenyong1,2

1. State Key Laboratory of Ocean Engineering

2. Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration, Shanghai Jiao Tong University, Shanghai 200240, China

通讯作者: 薛鸿祥,男,教授,博士生导师,电话(Tel.):021-34204470;E-mail:hongxiangxue@sjtu.edu.cn.

责任编辑: 陈晓燕

收稿日期: 2021-04-30  

基金资助: 国家自然科学基金(51979166)
国家自然科学基金(51909163)

Received: 2021-04-30  

作者简介 About authors

刘庆升(1990-),男,黑龙江省哈尔滨市人,博士生,从事非黏结柔性立管研究.

摘要

非黏结柔性立管常用于将海底油气资源输送到平台,有时因设计需要会包含复合材料圆柱壳层结构.以一类8层的非黏结柔性立管为研究对象,分析了含有复合材料圆柱壳层结构的非黏结柔性立管的弯曲特性.提出了一种作用于非黏结柔性立管的轴对称载荷和弯曲载荷联合作用的理论模型,同时建立了计及详细几何特性的非黏结柔性立管数值模型进行了验证.分析了两种典型复合材料及其纤维体积分数对非黏结柔性立管弯曲特性的影响.理论与数值结果吻合较好,计算结果表明:非黏结柔性立管的弯曲特性受复合材料圆柱壳层轴向弹性模量影响较大,轴向弹性模量较大的复合材料会较大地增强完全滑移阶段的非黏结柔性立管的弯曲刚度.非黏结柔性立管的轴向抗拉刚度和完全滑移阶段的弯曲刚度都随着复合材料圆柱壳层的纤维体积分数增大,接近线性增强.

关键词: 非黏结柔性立管; 复合材料; 圆柱壳层; 弯曲刚度; 纤维体积分数

Abstract

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: unbonded flexible risers; composite materials; cylindrical layers; bending stiffness; fibre volume fraction

PDF (3186KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

刘庆升, 薛鸿祥, 袁昱超, 唐文勇. 含复合材料结构的非黏结柔性立管弯曲特性[J]. 上海交通大学学报, 2022, 56(9): 1247-1255 doi:10.16183/j.cnki.jsjtu.2021.144

LIU Qingsheng, XUE Hongxiang, YUAN Yuchao, TANG Wenyong. Bending Properties of Unbonded Flexible Risers with Composite Materials[J]. Journal of Shanghai Jiaotong University, 2022, 56(9): 1247-1255 doi:10.16183/j.cnki.jsjtu.2021.144

非黏结柔性立管被广泛应用于海洋油气业,常用于从海底到平台的油气资源传输,有时会根据设计的需要增加耐腐蚀、隔热层或是防止发生径向失效的复合材料圆柱壳层结构[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 复合材料层弹性常数计算

复合材料通常可以简化为正交各向异性的材料进行研究,基于复合材料截面的等应力假设可以求解复合材料的弹性常数,相应的推导过程可以参考Liu等[17-18]的研究成果.

1.2 复合材料圆柱壳层在轴对称载荷作用下的理论模型

非黏结柔性立管轴对称载荷作用下的理论模型是弯曲载荷作用下理论模型的基础,需要基于轴对称载荷作用下的理论模型求解层间接触力,并基于该结果求解螺旋层结构的临界曲率,判断螺旋层是否发生滑移[19].本文对复合材料圆柱壳层模型进行了一定简化,忽略了复合材料的分层以及制造过程中的损伤,将宏观的复合材料圆柱壳层结构视为一个整体,在此基础上基于受力和变形关系推导相应的理论模型.

假设作用于圆柱壳层结构的外力载荷为轴向力F、绕中心轴转矩T、内外压力PiPo,如图1所示.图中,x、y、z方向分别为立管的周向、径向和轴向;x'y'z'方向分别为立管受力变化后的坐标轴方向;定义该层结构的初始长度L和轴向变化量ΔL;平均半径Rm和径向变化量ΔR;厚度t和厚度变化量Δt;绕轴扭转角度Δϕ.通过胡克定律和功能原理可以求解轴对称载荷作用下的平衡方程:

图1

图1   圆柱壳层受力及变形示意图

Fig.1   Forced and deformation schematic of cylindrical layer


k11k12k13k14k21k22k23k24k31k32k33k34k41k42k43k44ΔL/LΔϕ/LΔR/RmΔt/t=  F+πPiRi2-πPoRo2T2πRm(PiRi-PoRo)-πt(PiRi+PoRo)

式中:

k11=(1-μ23μ32)/(E2E3λ)k13=k31=(μ21+μ31μ23)/(E2E3λ)k14=k41=(μ31+μ21μ32)/(E2E3λ)k12=k21=k23=k32=k24=k42=0k22=G12Jk33=(1-μ13μ31)/(E1E3λ)k44=(1-μ12μ21)/(E1E2λ)k34=k43=(μ32+μ12μ31)/(E1E3λ)J=π(Ro4-Ri4)/2λ=π(Ro2-Ri2)E1E2E31-μ21-μ31-μ121-μ32-μ13-μ231

E1E2E3为材料x、y、z向弹性模量;G12为切变模量;相关参数μ为泊松比;RiRo分别为内、外半径;J为扭转惯性矩;λ为自定义的参数.

其他层结构的平衡方程推导可以参考任少飞等[4-5]的研究成果,此外,考虑层间的几何关系可以得到轴对称载荷作用下的平衡方程,求解层间接触压力大小[4,19].

1.3 复合材料圆柱壳层在弯曲载荷作用下的理论模型

先假设作用于圆柱壳层的轴对称载荷和弯曲载荷可以进行解耦,只考虑在单一方向的弯矩作用,并忽略圆柱壳层其他方向的曲率和挠率变化,相应受到的弯矩作用Mc和曲率κ的关系可以简化为

Mc=E1Iκ
I=π(Ro4-Ri4)/4

式中:I为惯性矩.

结合其他圆柱壳层的弯矩-曲率关系,并考虑螺旋层结构的未滑移、部分滑移和完全滑移阶段,建立非黏结柔性立管在弯曲载荷作用下的理论模型[4-5,19].

1.4 轴对称载荷与弯曲载荷联合作用下的理论模型

假设圆柱壳层结构在发生弯曲后,内外压力仍沿管长方向均匀分布,如图2所示,假设仅作用外压载荷,微元面积上受到的压力作用可以表示为[20-21]

图2

图2   弯曲圆柱壳层单元受力示意图

Fig.2   Forced schematic of curved cylindrical layer


dFo=PoRo(r+Rocosφ)dφdβ

式中: r为曲率半径;φ为绕圆柱壳层单元周向的角度;β为圆柱壳层曲率中心到该位置的夹角.

利用式(4)对角度φβ进行积分,可以得到微元上的压力在图2坐标轴各方向的分量:

 Fxo=β=θ-dβ2β=θ+dβ2φ=0φ=2πdFxo=2AoPosindθ2cosθFyo=β=θ-dβ2β=θ+dβ2φ=0φ=2πdFyo=-2AoPosindθ2sinθFzo=0

式中:θ为微元中心到曲率中心的角度;圆柱壳层结构外径横截面积AoRo2.

忽略轴对称载荷引起的圆柱壳层结构的轴向伸长,假设圆柱壳层结构在轴对称载荷和弯曲载荷联合作用后各处曲率保持一致,如图3所示.图中,RP1和RP2分别代表圆柱壳层两端的几何中心端点; O为曲率中心.并考虑弯曲后的圆弧角度δ与曲率的关系为δ=Lκ,由轴对称载荷引起的对弯曲特性的影响Mp可以表示为

图3

图3   圆柱壳层弯曲后构型示意图

Fig.3   Schematic of curved cylindrical layer after bending


Mp=-L/20dMp+0L/2dMp=-Lκ/20r-FxorcosLκ2-rcosθ+Fyorsinθ+0Lκ/2r[-FxorcosLκ2-rcosθ+Fyorsinθ+sinLκ2]=AoPoκ22sinLκ2-sinLκ

对于单位长度的圆柱壳层结构,定义变量a

a=2sinκ2-sinκ/κ2

求解曲率0到1范围内的变量a并进行线性拟合,得到如图4所示的曲线,忽略线性拟合后的曲线截距得到简化后的变量a的表达式为

a0.118κ

图4

图4   a-κ关系及简化曲线

Fig.4   a-κ curves and simplified curves


相应的外压载荷引起的弯矩可以表示为

Mp=0.118κPoAo

同理可以推导内压和轴向力载荷的影响.

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   非黏结柔性立管剖面和有限元模型示意图

表1 非黏结柔性立管几何特性及材料特性

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骨架层131.60-87.5AISI304
2内部护套层-35.10-尼龙12
3抗压铠装层240.05-85.5FI-15
4加强层-46.25-复合材料
5抗拉铠装层4047.75-35.0FI-41
6加强层-50.75-复合材料
7抗拉铠装层4452.2535.0FI-41
8外部护套层-55.25--

新窗口打开| 下载CSV


3 计算结果及分析

3.1 数值和理论方法验证

非黏结柔性立管在弯曲载荷作用下的理论模型以轴对称载荷作用下的理论模型为基础进行研究.首先通过轴对称载荷作用下的理论模型求解非黏结柔性立管不同层之间的接触压力,在此基础上求解螺旋层结构的弯曲临界曲率,之后加载弯矩载荷的作用,同时结合本文提出的轴对称载荷和弯曲载荷联合作用下的理论模型,求解非黏结柔性立管在未滑移、部分滑移和完全滑移阶段的非黏结柔性立管的弯曲特性.对于数值方法,载荷的加载方式与理论方法保持一致,同样地首先施加轴对称载荷,之后再在两端几何参考点RP1和RP2施加弯矩为M的作用.采用显式算法进行求解时需要考虑惯性效应对结果的影响[23],因此本文采用准静态的方法对载荷进行加载,并模拟足够的加载时长以避免惯性效应可能对结果造成的影响[4],并始终控制动能与内能的比率在5%以内.

根据Witz[14]试验对轴向张力作用下和弯曲载荷作用下的非黏结柔性立管特性进行验证.需要指出的是,Witz[14]的弯曲试验中未施加轴对称载荷的作用,但是非黏结柔性立管由于制造过程中存在的初始缺陷等因素,层间始终存在初始的接触压应力作用[6],本文通过对最外部的护套层施加初始外压载荷,给出层间的接触压力作用,初始的外压载荷取为0.2 MPa[4],并先于弯矩载荷进行加载.

通过本文提出的理论和数值方法计算Witz[14]试验中的工况,相应的计算结果如图67所示,得到相应的截面力学性能如表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  Calculation results of different methods

方法力学性能
Ea/MNEi/(kN·m2)
理论105.881.19, 1.42
数值99.191.41
试验91.191.19, 1.04

新窗口打开| 下载CSV


需要指出的是,相比于数值方法,本文所提理论方法具有极高的计算效率,且可以得到与数值方法近似的计算结果.本文的数值方法在某具有48核的工作站上开展,弯曲载荷工况的载荷加载时间为 1.8 s,单一工况计算需要约72 h.本文的理论计算方法通过MATLAB软件进行编程,计算时间约为 15 s.

3.2 含复合材料圆柱壳层立管弯曲特性

选取两种海洋油气业常用复合材料替换原防摩擦层结构,相应的纤维和基体材料特性来自文献[18,24],整体材料特性如表3所示.两种复合材料的纤维体积分数(φvf)都在7%左右,其中,第1种材料以环氧树脂为基体,玻璃纤维为增强材料;第2种材料以高密度聚乙烯为基体,钢丝为增强材料,两种复合材料都表现出了明显的正交各向异性.

表3   复合材料的弹性常数

Tab.3  Elastic constants of composite materials

材料弹性模量/GPaμ12μ13μ23G12/GPa
E1E2E3
19.404.024.010.360.360.494.08
216.01.073.860.440.440.450.56

新窗口打开| 下载CSV


同样对于含有复合材料圆柱壳层的非黏结柔性立管模型首先施加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   含复合材料圆柱壳层非黏结柔性立管弯矩-曲率曲线

表4 含复合材料圆柱壳层弯曲刚度计算结果

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.5111.7411.52
钢丝/高密度聚乙烯18.9219.1518.57

新窗口打开| 下载CSV


表5   完全滑移阶段各层结构对非黏结柔性立管模型的弯曲刚度贡献

Tab.5  Contribution of each layer to bending stiffness of unbonded flexible riser model in full-slipping stage

层结构及轴对称载荷影响Ei/(kN·m2)
含防摩擦层尼龙-11模型含复合材料玻纤/环氧树脂模型含复合材料钢丝/高密度聚乙烯模型
第1, 2, 3, 5, 7, 8层结构的综合作用0.85
第4层0.154.617.81
第6层0.196.0510.26
轴对称载荷作用0.23 (0.2 MPa外压载荷作用)
总计1.4211.7419.15

新窗口打开| 下载CSV


3.3 纤维体积分数的影响分析

考虑纤维体积分数对含有复合材料圆柱壳层的非黏结柔性立管弯曲特性的影响,选取几种海洋油气业常用的增强纤维材料进行研究,包括常用于组成防止非黏结柔性立管发生径向失效的强结构带的凯夫拉49增强纤维[4],所有的基体材料都采用环氧树脂,相应的材料参数如表6所示.

表6   材料特性

Tab.6  Material properties

材料类型弹性模量/GPa泊松比
普通碳纤维2210.30
钢丝2100.26
凯夫拉491310.36
玻璃纤维880.36
环氧树脂3.50.33

新窗口打开| 下载CSV


考虑复合材料圆柱壳层的纤维体积分数从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]. 油气储运, 2016, 35(12): 1255-1260.

[本文引用: 1]

YU Ronghua, YUAN Pengbin.

Structure and research focus of marine unbonded flexible pipes

[J]. Oil & Gas Storage and Transportation, 2016, 35(12): 1255-1260.

[本文引用: 1]

EPSZTEIN T, DEMANZE F, LEFEBVRE X, et al.

New anti H2S layer for flexible pipes

[C]// Offshore Technology Conference. Houston, USA: OTC, 2011.

[本文引用: 1]

MALTA E R, MARTINS C D A.

Finite element analysis of flexible pipes under compression

[C]// International Conference on Ocean, Offshore and Arctic Engineering. San Francisco, USA: ASME, 2014.

[本文引用: 1]

任少飞, 唐文勇, 薛鸿祥.

非黏结柔性立管骨架层失效数值计算方法

[J]. 上海交通大学学报, 2016, 50(3): 466-471.

[本文引用: 9]

REN Shaofei, TANG Wenyong, XUE Hongxiang.

Numerical method to predict failure of carcass layer of unbonded flexible risers

[J]. Journal of Shanghai Jiao Tong University, 2016, 50(3): 466-471.

[本文引用: 9]

ELIZBAR K. Theoretical modelling of unbonded flexible pipe cross-sections[D]. London: South Bank University, 2000.

[本文引用: 4]

KRAINCANIC I, KEBADZE E.

Slip initiation and progression in helical armouring layers of unbonded flexible pipes and its effect on pipe bending behaviour

[J]. The Journal of Strain Analysis for Engineering Design, 2001, 36(3): 265-275.

DOI:10.1243/0309324011514458      URL     [本文引用: 3]

DONG L L, HUANG Y, ZHANG Q, et al.

An analytical model to predict the bending behavior of unbonded flexible pipes

[J]. Journal of Ship Research, 2013, 57(3): 171-177.

DOI:10.5957/jsr.2013.57.3.171      URL     [本文引用: 2]

FÉRET J J, BOURNAZEL C L.

Calculation of stresses and slip in structural layers of unbonded flexible pipes

[J]. Journal of Offshore Mechanics and Arctic Engineering, 1987, 109(3): 263-269.

DOI:10.1115/1.3257019      URL     [本文引用: 1]

KEBADZE E, KRAINCANIC I.

Non-linear bending behaviour of offshore flexible pipes

[C]// International Offshore and Polar Engineering Conference. Brest, France: ISOPE, 1999: 226-233.

[本文引用: 1]

RAMOS R J, PESCE C P.

A consistent analytical model to predict the structural behavior of flexible risers subjected to combined loads

[J]. Journal of Offshore Mechanics and Arctic Engineering, 2004, 126(2): 141-146.

DOI:10.1115/1.1710869      URL     [本文引用: 1]

RAMOS R, PESCE C P, MARTINS C A.

A new analytical expression to estimate the bending stiffness of flexible risers

[C]// International Conference on Offshore Mechanics and Arctic Engineering. Honolulu, USA: ASME, 2009: 513-522.

[本文引用: 1]

MAGLUTA C, ROITMAN N, VIERO P F, et al.

Experimental estimation of physical properties of a flexible riser

[C]// International Conference on Offshore Mechanics and Arctic Engineering. Rio de Janeiro, Brazil: ASME, 2001: 255-264.

[本文引用: 1]

TROINA L M B, ROSA L F L, VIERO P F, et al.

An experimental investigation on the bending behaviour of flexible pipes

[C]// International Conference on Offshore Mechanics and Arctic Engineering. Honolulu, USA: ASME, 2009: 637-645.

[本文引用: 1]

WITZ J A.

A case study in the cross-section analysis of flexible risers

[J]. Marine Structures, 1996, 9(9): 885-904.

DOI:10.1016/0951-8339(95)00035-6      URL     [本文引用: 5]

刘晓媛, 薛鸿祥, 唐文勇.

轴向承载力作用下非粘结柔性立管的侧向失效分析

[J]. 上海交通大学学报, 2018, 52(9): 1017-1022.

[本文引用: 1]

LIU Xiaoyuan, XUE Hongxiang, TANG Wenyong.

Lateral failure analysis of unbonded flexible risers under axial compression force

[J]. Journal of Shanghai Jiao Tong University, 2018, 52(9): 1017-1022.

[本文引用: 1]

ZHANG M M, CHEN X Q, FU S X, et al.

Theoretical and numerical analysis of bending behavior of unbonded flexible risers

[J]. Marine Structures, 2015, 44: 311-325.

DOI:10.1016/j.marstruc.2015.10.001      URL     [本文引用: 1]

LIU Q S, XUE H X, TANG W Y, et al.

Theoretical and numerical methods to predict the behaviour of unbonded flexible riser with composite armour layers subjected to axial tension

[J]. Ocean Engineering, 2020, 199: 107038.

[本文引用: 1]

朱彦聪. 钢丝缠绕增强塑料复合管外压失稳研究[D]. 杭州: 浙江大学, 2007.

[本文引用: 2]

ZHU Yancong. Buckling analysis of plastic pipe reinforced by winding steel wires under external pressure[D]. Hangzhou: Zhejiang University, 2007.

[本文引用: 2]

董磊磊. 非粘合柔性立管截面特性的理论计算及BSR区域的疲劳分析[D]. 大连: 大连理工大学, 2013.

[本文引用: 3]

DONG Leilei. Theoretical prediction of cross-sectional properties and fatigue analysis in the BSR area for unbonded flexible risers[D]. Dalian: Dalian University of Technology, 2013.

[本文引用: 3]

CHAKRABARTI S K, FRAMPTON R E.

Review of riser analysis techniques

[J]. Applied Ocean Research, 1982, 4(2): 73-90.

DOI:10.1016/S0141-1187(82)80002-3      URL     [本文引用: 1]

SEYED F B, PATEL M H.

Mathematics of flexible risers including pressure and internal flow effects

[J]. Marine Structures, 1992, 5(2/3): 121-150.

DOI:10.1016/0951-8339(92)90025-K      URL     [本文引用: 1]

ALI B. Development of a constitutive model to simulate unbonded flexible riser pipe elements[D]. London: Brunel University, 2008

[本文引用: 1]

庄茁, 张帆, 岑松, . ABAQUS非线性有限元分析与实例[M]. 北京: 科学出版社, 2005.

[本文引用: 1]

ZHUANG Zhuo. Nonlinear finite element analysis and examples of ABAQUS[M]. Beijing: Science Press, 2005.

[本文引用: 1]

王人杰.

纤维增强复合材料横向弹性常数

[J]. 复合材料学报, 1996, 13(2): 98-104.

[本文引用: 1]

WANG Renjie.

Transverse elastic constants for fiber reinforced composite materials

[J]. Acta Materiae Compositae Sinica, 1996, 13(2): 98-104.

[本文引用: 1]

/