考虑涡激振动变形影响的海底管道静态悬跨分析.doc

考虑涡激振动变形影响的海底管道静态悬跨分析 崔少敏 （海洋石油工程股份有限公司，天津 300451） 摘 要 根据DNV在2006年提出的海底管道悬跨长度评估指导性规范DNV-RP-F105，本文给出了一种更加切合实际的海底管道临界静态悬跨长度的计算方法。该方法考虑了由于涡激振动（VIV）））可用顺流向响应模型表示： (1) 式中：为流速比的修正系数；为与应力幅值相乘的安全因子，一般取1.3； 为单位应力振幅（单位直径顺流向模态变形引起的应力）； （2） 式中：为钢管外径；为壁厚；为管道外径（包含所有涂层）；为边界因子；为混凝土增强因子；为有效跨长；为钢材的杨氏模量。 为最大顺流向VIV响应振幅，它可根据衰减速度从下图插值得到。 图1 顺流向VIV响应振幅 上图中四个顶点坐标值可根据以下公式进行计算： (3) 式中：为稳定性系数；和为衰减因子。 横向VIV引起的应力幅值可用横向响应模型表示： 式中：为横向单位应力振幅， 为波流联合作用下的横向VIV振幅，它可根据从下图的得到。 图2 横向VIV响应振幅 上图中五个顶点坐标值可根据以下公式进行计算： (4) 式中；为两相邻横向模态频率比值。 1.2 临界静态悬跨计算 根据以上公式可计算得出给定波流条件下的顺流向应力幅值和横向应力幅值。将结果代入下面公式中计算最大环境应力： （5） （6）为由直接波浪力产生的最大环境应力。 环境弯矩由下式计算得出： （7） 此处弯矩值为最终环境弯矩。另计算临界跨长前还需计算出功能性荷载（如重力、路由弯曲）产生的弯矩值和管道的有效轴力值。将这些荷载值全部计算得出后，就可根据DNV-OS-F101中荷载组合校核公式进行临界跨长的计算了。相关公式如下： ，时 （8） ，时 （9） 式中：为考虑了环境弯矩和功能弯矩的设计弯矩；为设计有效轴力；为设计过压压差；为塑性弯矩抗力；为特征塑性轴向抗力；为破裂压力；为考虑应变硬化的流动应力参数。 通过以上分析已经获得管道的临界静态悬跨，设计中还需考虑管道的临界涡激振动悬跨，最终结果取二者小值，工程中则需要把最大实际悬跨控制在临界悬跨之下来防止管道不发生静态屈曲和涡激振动。 2 数值算例 下面以某海底管道为例，按照上述方法计算其临界静态悬跨并与其涡激振动悬跨进行比较。 海底管道参数：外径323.9mm，壁厚12.7mm，内径219.1mm，壁厚14.3mm，外管涂层2.8mm，密度940kg/m3。 设计波高：一年一遇最大波高7.3m，最大波周期10.7s，有效波高12.5m，有效波周期11.2s；十年一遇最大波高9.6m，最大波周期12.1s，有效波高16.6m，有效波周期12.6s；百年一遇最大波高12.1m，最大波周期13.8s，有效波高20.8m，有效波周期14.4s。 设计流速：一年一遇海底流速0.371m/s；十年一遇海底流速0.532m/s；百年一遇海底流速0.742m/s。 设计水深为121.8m。 根据计算分析，临界悬跨长度为： 工况 静态悬跨(m) 涡激振动悬跨(m) 临界悬跨 安装工况 126.747 27.382 27.382 水压试验工况 36.692 23.253 23.253 操作工况 19.507 15.258 15.258 根据上述结果分析得出以下结论： （1）本算例中，安装工况、水压试验工况和操作工况下的临界悬跨长度依次递减。 （2）海底管道三种工况下的静态临界悬跨长度小于临界涡激振动悬跨长度，即最终结果都是由涡激振动临界悬跨长度控制的。 （3） 同时也可看出在操作工况下，临界静态悬跨长度与临界涡激振动悬跨长度已比较接近，有可能成为控制准则。实际上在其他项目中也有临界静态悬跨成为控制准则的情况。 （4） 临界静态悬跨长度随着工况的不同差别较大，这主要是由于内部介质和安全系数的不同及腐蚀情况的考虑与否造成的。 3 结语 本文根据DNV-RP-F105规范介绍了管道临界静态悬跨的计算方法，该方法考虑了波流联合作用下管道振动对临界静态悬跨的影响，更加符合规范的要求，可供工程设计人员参考。 工程实践表明，即使相应规范及计算方法的不断优化改进，仍很难对海底管道临界悬跨长度进行准确的评估。实际情况中，管土相互作用的复杂情况、实际涡激状况及管道内在缺陷等都对管道临界悬跨有着较大的影响，此类影响点都有待于进一步研究。 参考文献 1 董丽丽，徐慧，郭振邦：海底管