半潜式平台撑杆结构波浪载荷计算方法研究

1、半潜式平台撑杆结构波浪载荷计算方法研究李辉，任慧龙，吴东伟，张健(哈尔滨工程大学船舶工程学院 哈尔滨市 150001)摘 要：考虑平台与波浪的相对运动，对适用于空间倾斜构件的Morison公式进行了线性化处理。将依据该公式算得的Morison力计入基于三维势流理论建立的平台整体运动方程中，采用迭代法进行求解。对一艘半潜式平台的计算结果表明，平台撑杆所遭受的Morison力对平台整体运动的影响较小。关键词：半潜式平台；撑杆；Morison公式；线性化随着全球经济的发展，人类对油气资源的需求与日俱增。陆地上资源的日益枯竭，人们逐渐把目光投向了海洋，并正在向深海区域发展。半潜式平台因其自身明显的优势

2、逐渐成为海洋油气资源开发最主要的工具之一。由于半潜式平台的工作环境一般比较恶劣，发生事故的可能性较大，因而在设计阶段就要对其结构进行合理的强度评估，以保证其结构的安全性。而准确的预报平台在使用期间可能遭受的波浪载荷是强度评估正确进行的关键。对于海洋平台波浪载荷的计算，一般可采用两种方法1。一种是采用适用于大尺度构件的基于势流理论的波浪载荷计算方法，另一种是采用基于经验的适用于小尺度构件的Morison公式。半潜式平台在水面以下既有尺度较大的构件如旁通和立柱，又有尺度较小的构件如撑杆，因而在确定其波浪载荷时，需要将两种方法结合起来使用2。由于Morison公式中含有非线性项（拖曳力项），要与目前

3、较为通用的基于三维线性势流理论的波浪载荷计算方法联合使用，需要对其进行线性化处理。为了合理预报半潜式平台的运动与波浪载荷，将Morison公式与三维线性势流理论联合使用。采用考虑平台与波浪相对运动的适用于倾斜构件的Morison公式，应用Borgman方法3进行线性化。将采用Morison公式算得的撑杆结构的波浪载荷计入平台的整体运动微分方程。结合实际算例对半潜式平台撑杆结构受力对于平台整体运动的贡献进行了分析。1空间倾斜构件的Morison公式及线性化方法1.1 相对速度下倾斜构件的Morison公式半潜式平台的撑杆结构是一组空间倾斜构件，在与波浪具有相对速度的情况下，作用于其上的波浪载荷可

4、采用如下的Morison公式进行计算4。 （1）其中：平台撑杆的直径；为撑杆的截面积；=水质点的速度；=水质点的加速度。=平台的运动速度；=平台运动的加速度；相对速度， ；构件的法向相对速度：+为转换矩阵，将速度和加速度变换到与构件轴线垂直的方向。在转换矩阵中，为杆轴线方向单位矢量，夹角与定义如图1所示。图1. 空间杆件示意图1.2 相对速度下倾斜构件线性化的Morison公式为了与三维线性势流理论联合使用，准确计算半潜式平台的波浪载荷，需要对Morison公式中的非线性项（拖曳力项）进行线性化处理。本文采用Borgman线性化方法对（1）式进行线性化处理，得到相对速度下倾斜构件线性化的Mor

5、ison公式式（2）。 （2）式中为杆件法向相对速度均方差。均方差的求解方法是线性化的关键，本文考虑的是在随机波浪下的线性化问题。利用谱分析方法速度方差表示为： （3）其中，为杆的法向相对速度谱。当求得后，也就解决了线性化的问题。2平台在规则波中运动微分方程基于三维势流理论的波浪载荷计算方法目前已发展比较成熟，在船舶与海洋工程领域有着非常广泛的应用。在此方法中，浮体运动微风方程可以表示为式（4）的形式5。 （4）其中，波浪频率；浮体的质量矩阵；浮体的附加质量矩阵；浮体的阻尼矩阵；浮体的静恢复力矩阵；流体载荷列向量；当计入小尺度构件的受力时，需要在方程的右端增加线性化的Morison力项，即：

6、（5）对于式（5）的求解，可采用迭代法，其基本求解步骤如图1所示。图2迭代程序示意图3算例分析本文以LiuHua半潜式平台为例，对平台撑杆的波浪载荷进行了计算。并与平台主体的波浪载荷叠加，重新对平台的整体运动进行了求解，通过对计算结果的分析，考察了撑杆结构的受力对平台整体运动的影响。3.1 平台基本参数LiuHua半潜式平台的主要尺度参数见表1。平台结构模型见图3，平台撑杆结构的模型见图4。表1 平台各项尺度参数尺度参数米（m）英尺（ft）主甲板长90.2273.46296宽73.46241高39.62130纵向浮筒长89.93295宽15.2450高6.4021横向浮筒宽6.4021高3.6

7、612角立柱直径9.1430中间立柱直径5.4918工作吃水22.8675图3 平台整体结构图4 平台撑杆的布置形式3.2 计算结果选取平台的工作海况（HS=9m，TZ=12.5s）对Morison公式进行线性化处理，进而对LiuHua半潜式平台的运动响应进行了迭代计算，迭代计算过程及结果见表2表4。由于半潜式平台的摇摆运动幅值较小，因而表中只列出了平动位移幅值。表2 0度浪向运动幅值变化OmeSurgeSwayHeave原始运动0.5980.208500.3332第一次迭代0.5980.200900.3247第二次迭代0.5980.200300.3245第三次迭代0.5980.200300.

8、3245运动偏差3.93%0.00%2.61%表3 90度浪向运动幅值变化OmeSurgeSwayHeave原始运动0.96700.310.0687第一次迭代0.96700.30290.0682第二次迭代0.96700.30310.0682第三次迭代0.96700.30310.0682运动偏差0.00%2.20%0.73%表4 150度浪向运动幅值变化OmeSurgeSwayHeave原始运动0.5030.3850.21580.4167第一次迭代0.5030.36590.21570.4012第二次迭代0.5030.36530.21570.4007第三次迭代0.5030.36530.21570.

9、4007运动偏差5.31%0.46%3.83%通过计算结果可以看出，Morison力对于平台整体运动的影响是很小的。由于计算此力的方法繁琐，在平台整体载荷计算与运动响应分析中，可以忽略不计。4.结论本文将经线性化处理的Morison公式与三维线性势流理论联合使用，建立了平台在规则波中的运动微分方程，采用迭代方法进行求解。采用该方法对一艘半潜式平台的运动响应进行了分析。计算结果表明半潜式平台的撑杆结构遭受的波浪力对平台整体运动的影响很小，在平台整体载荷计算和运动响应分析中，可以忽略不计。参考文献：1 李润培,王志农. 海洋平台强度分析M. 上海:上海交通大学出版社,1992.2 刘海霞,肖熙. 半潜式平台结构强度分析中的波浪载荷计算,中国海洋平台J, 2003,18(2):1-4.3 徐继祖,李维扬,汪克让,张振衡,聂武. 海洋工程结构动力分析M. 天津:天津大学出版社，19