固定式海洋平台的时域振动响应研究

谢　卓，张火明，方贵盛，孔令滨（. 中国计量学院 浙江流量计量技术重点实验室，浙江 杭州 3008；. 浙江水利水电学院 机械与汽车工程学院，浙江 杭州 3008）

固定式海洋平台的时域振动响应研究

谢卓1，张火明1，方贵盛2，孔令滨1
（1. 中国计量学院 浙江流量计量技术重点实验室，浙江 杭州 310018；2. 浙江水利水电学院 机械与汽车工程学院，浙江 杭州 310018）

近海结构的研究是海洋工程领域的重要专题。本文以 Airy 线性波理论和莫里森方程为基础推导平台结构的运动方程。将波浪力简化为二节点集中载荷，从动力学的角度分析结构的响应特征。以一简化的固定平台为算例，以不同的入射波为激励条件分析比较结构的位移响应、速度响应以及加速度响应。研究结果表明，随着入射波强度的增大，结构的振动变得更为剧烈，且结构的振动响应不仅与入射波的周期和波高有关，与波浪的形式也有一定的关系。

固定式海洋平台；时域；振动响应

0　引　言

作为海洋石油、天然气资源开发的重要基地，海洋平台的开发、研究和利用越来越受到人们的广泛关注［1］；从渤海 2 号被推倒，我国相继开展了对于固定式平台的研究并取得一系列研究成果。近海结构的受力研究是海洋工程领域的重要专题，田宏升［2］研究计算了小构件的波浪载荷。张金平［3］综述了近年来海洋平台波浪理论及波浪载荷计算方法的研究与发展概况。何懋华［4］分析研究了冰载荷作用下的固定式平台的极限承载能力以及波浪荷载作用下的固定式平台的极限承载能力。刘为民等［5］从入射波角度、波高以及周期方面探讨了规则波下张力腿平台的动力特性。Sang-Yun BAE［6］采用高阶边界元法对非线性波浪流共同作用下的 FPSO 进行了数值模拟，并对衍射和辐射波进行了时域计算。Jimmy Ng K.T 等［7］于 2010 年对只有在波、流以及波流联合作用下的多立柱进行试验并进行频谱、受力分析，研究了上游立柱尾流对下游立柱的影响及下游立柱的动力响应。海洋平台结构的可靠性尤为重要。一旦结构失效，不仅会造成巨大的经济损失，而且还会有严重的人员伤亡和环境污染［8 - 10］。

本文从动力学的角度分析时域内结构的响应特征，以不同的入射波作为激励条件，分析比较结构的位移、速度以及加速度响应，并得出一些有意义的结论，从而可以为固定式海洋平台的相关设计提供理论计算依据。

1　理论分析

Airy 线性波理论是波浪理论中最基本最常用的波浪理论。Airy 波浪理论将非线性的波浪自由面条件，近似以线性的边界条件代替，这一线性边界适应于波高很小的情况。Airy 线性波的速度势为：

波面表达式为：

式中：H 为波高；k 为波数，k = 2π/L；L 为波长；d 为水深。

图 1　波浪示意图Fig. 1　Schematic diagram of wave

设 θ = kx - ωt，由线性波速度势可以得到波浪水质点速度和加速度的表达式为：

莫里森方程作为计算海洋结构物上波浪力的常用方法，其计算的波浪力包括与水平速度的平方成正比的阻力项以及与水平加速度成正比的惯性力项。阻力项是由于水质点的速度引起的，而惯性力项则是由于水质点的加速度引起的。其表达式为：

式中：AI= 0.25 ρπD2；AD= 0.5 ρD；CD和 CM为拖曳力系数和惯性力系数；D 为桩柱直径；u 为海水质点水平速度；为海水质点的水平加速度。

将式（4）和式（5）带入式（6），整理得到线性波理论波浪力的表达式为：

式中：

式（7）同样适用于不规则的波浪场，可以将不规则的波浪场的波面表达函数 η（x，t）看作是多个不同频率和不同幅值的规则波浪场相叠加的结果。对波浪力 f（x，t）进行积分，可求得整根柱体上所受的波浪力。

固定式平台的振动方程为

采用模态叠加法求解振动位移响应，通过各个模态振动响应的叠加可求得几何坐标表示的位移即 y1= y1（t）和 y2= y2（t）解。其中 λdi和 λi分别为有阻尼固有频率和无阻尼固有频率，ξi为结构阻尼比。

进一步可求得平台振动的速度和加速度为：

2　算例分析

以一简化的具有上下 2 层甲板的固定平台为算例，文中不牵涉流载荷分布问题，故简化为 2 层甲板完全可取。其结构示意图如图 2 所示，工作水深为d，集中质量到上下 2 层甲板，其中 2 层甲板质量为m1和 m2，平台具有 2 个质量独立的响应位移分别为 y1和 y2。

图 2　平台结构示意图Fig. 2　Schematic diagram of platform structure

2.1参数设置

算例采用固定式平台的参数如表 1 所示，其中 Kij为平台的刚度系数，ξ 为阻尼比，Nl= 4 表明平台为四垂直立柱架构，Nc= 2 为平台的 2 个平行支撑。采用的入射波参数如表 2 所示。由表 1 和 表 2 可知：此平台垂直立柱外径 Dl= 5.5 m，入射波波长 Lw＞ 224 m，Dl/Lw＜ 0.024 ＜ 0.2，满足 Morison 方程的使用条件。结构的直径尺寸相对于波浪要小得多，假设结构的运动相对于波浪的运动小得多，这样 Morison 方程中的惯性系数 CM和流体阻力中的 CD相比，前者处于主导地位，即惯性体系下，流体流动占主导地位。

2.2计算结果

以规则波为入射波，根据某海域气象资料，取波浪周期 12.39 s，波高 12.5 m，计算得到波长 L = 224.243 m，如表 2 中 η1所示，进行计算并得到平台上甲板和下甲板的位移 y、速度 v 及加速度 a，对比得到结果如图 4 ～图 6 所示。

表 1　固定式平台参数Tab. 1　stationary platform parameters

表 2　波浪参数Tab. 2　Wave parameters

图 3　一年一遇海况规则波面高程Fig. 3　Once a year the sea surface elevation rule

图 4　平台位移对比图Fig. 4　Platform displacement contrast diagram

图 5　平台速度对比图Fig. 5　Platform speed comparison chart

图 6　平台加速度对比图Fig. 6　Platform acceleration contrast diagram

由图 4-图 6 可以观察到，平台架构的位移响应、速度响应及加速度响应在经过 2 个波浪周期时间的振荡后逐渐趋于稳定，且响应均在 t = 4 s 时达到峰值。同一时刻平台上甲板的位移响应、速度响应和加速度响应均要大于下甲板的。由此得出，平台在某一特定入射波激励下上甲板的振动相比下甲板更为剧烈。

采用 3 种不同条件的规则波为入射波进行计算，3种波浪即表 2 中的 η1，η2及 η3，周期及波高等参数如表 2 所列。计算得到平台结构的位移响应，速度响应以及加速度响应，所得结果如图 7 ～ 图 12 所示。

以 3 种不同条件的规则波为入射波，由上、下甲板的位移响应对比图 7 和 图 8 可知，结构上甲板和下甲板的位移响应均在 t = 4 s 处达到峰值，且波浪 η1激励下位移响应峰值最大，波浪 η3激励下最小；由图 9 和 图 10 可以看出结构上甲板和下甲板的速度响应均在 t = 3 s 处达到响应峰值，且 3 种波浪条件中 η3激励下的平台速度响应峰值大于其它波浪条件激励下的响应，η1时响应的峰值最小；由图 11 和图 12 得到结构上、下甲板的加速度响应均在 t = 4 s 处达到响应峰值，且同位移响应和速度响应一样，波浪 η3激励下的加速度响应峰值最大，η1时的响应峰值最小。综上可得结构的位移响应、速度响应及加速度响应随着入射波强度的增大而趋于增大，故随着入射波强度的增大，结构的振动也变得更为剧烈。

图 7　三种海况下结构上甲板位移响应图Fig. 7　Comparison of displacement response of structure on deck

图 8　三种海况下结构下甲板位移响应图Fig. 8　Comparison of displacement response of structure under deck

图 9　三种海况下结构上甲板速度响应图Fig. 9　Comparison of speed response of structure on deck

图 10　三种海况下结构下甲板速度响应图Fig. 10　Comparison of speed response of structure under deck

图 11　三种海况下结构上甲板加速度响应图Fig. 11　Comparison of acceleration response of structure on deck

图 12　三种海况下结构下甲板加速度响应图Fig. 12　Comparison of acceleration response of structure under deck

Pierson-Moskowitz 波高谱于 1966 年被国际船模水池会议定为标准海浪谱，目前广泛应用于海浪研究及有关工程问题。采用以 P-M 谱为随机波浪谱作为入射波，取波浪周期 12.39 s，有义波高 12.5 m，其波面高程如图 13 所示。计算结构位移、速度和加速度响应并与同条件下以规则波为入射波所得振动响应结果进行比较，限于篇幅这里仅给出位移响应的对比结果，如图 14 和图 15 所示。

图 13　P-M 谱的波面高程Fig. 13　Wave height of P-M spectrum

图 14　结构上甲板位移响应对比图Fig. 14　Comparison of displacement response of structure on deck

图 15　结构上甲板位移响应对比图Fig. 15　Comparison of displacement response of structure under deck

由图 14 和图 15 上、下甲板位移响应对比可看出在周期相同，并取有义波高与规则波波高相同的 2 种不同入射波激励时，结构位移响应会在同一时刻达到峰值，且以规则波为入射波时的位移响应峰值要大于P-M 谱激励下的位移响应峰值。由此可知，结构在规则波的激励下的振动响应要比相同强度以 P-M 谱为随机波浪谱的振动响应更为剧烈。综上所述，结构的振动响应不仅与入射波的周期和波高有关，与波浪的形式也有一定的关系。

3　结　语

本文分析了固定式平台的振动响应问题，将立柱上分布的波浪力处理为集中作用在上甲板和下甲板 2个结点上的集中载荷。对比分析了上、下甲板的振动响应；研究了以不同海况（一年一遇，十年一遇和百年一遇）的规则波为入射条件，固定式平台结构的响应以及相同海况不同形式入射波（规则波和以 P-M 谱为随机波浪谱）下结构的运动响应。主要结论如下：平台在入射波激励下上甲板比下甲板的振动更为剧烈；随着入射波强度的增大，结构的振动也变得更为剧烈；结构的振动响应不仅与入射波的周期和波高有关，与波浪的形式也有一定的关系。

［1］刘与嘉，张宝琳. 含控制时滞的海洋平台前馈反馈最优跟踪控制［J］. 中国计量学院学报，2014，25（1）：99-106. LIU Yu-jia，ZHANG Bao-lin. Feedforward and feedback optimal tracking control for offshore steel jacket platforms with control delay［J］. Journal of China Jiliang University，2014，25（1）：99-106.

［2］田宏升. 固定式海洋结构物载荷分析与研究［D］. 大连：大连理工大学，2005.

［3］张金平，段艳丽，刘学虎. 海洋平台波浪载荷计算方法的分析和建议［J］. 石油矿场机械，2006，35（3）：10-14.

［4］何懋华. 南海固定式平台极限承载能力分析研究［D］. 广州：华南理工大学，2013.

［5］刘为民，谷家扬，卢燕祥. 规则波下张力腿平台时域响应分析［J］. 舰船科学技术，2014，36（8）：27-32. LIU Wei-min，GU Jia-yang，LU Yan-xiang. Research on timedomain dynamic response of tension leg platform in regular wave［J］. Ship Science and Technology，2014，36（8）：27-32.

［6］BAE S Y. Numerical simulation of inviscid wave -current interaction with an FPSO［J］. Dissertations & Theses - Gradworks，2008.

［7］NG J K T，HALKYARD J E，SOON C E. Statistical characteristics of flow in the wake region of a vertical bluff cylinder in waves，currents and combined wave-current flows［C］//Proceedings of the ASME 2010 29th International Conference on Ocean，Offshore and Arctic Engineering，Shanghai，China：American Society of Mechanical Engineers，2010：989-998.

［8］修宗祥. 深水导管架海洋平台安全可靠性分析及优化设计［D］. 青岛：中国石油大学（华东），2010.

［9］张淑华，徐磊，钱进. 海洋导管架平台疲劳问题分析［J］. 石油矿场机械，2012，41（11）：16-19.

［10］吕秀艳. 固定式海洋平台结构风险评估及应用［D］. 青岛：中国海洋大学，2005.

Time domain vibration response of the fixed offshore platform

XIE Zhuo1，ZHANG Huo-ming1，FANG Gui-sheng2，KONG Ling-bin1
（1. Zhejiang Provincial Key Laboratory of Flow Measurement Technology，China Jiliang University，Hangzhou 310018，China 2. Zhejiang Water Conservancy and Hydropower College，Mechanical and Automotive Engineering，Hangzhou 310018，China）

The study of offshore structures is an important subject in the field of ocean engineering. Based on Airy linear wave theory and Morrison equation，the equation of motion of the platform was derived. The wave force was simplified to two nodes，and the response characteristics of the structure were analyzed from the point of view of dynamics. A simplified fixed platform was used as an example. The displacement response，velocity response and acceleration response of the structure were analyzed and compared with different incident waves. Research results show that the intensity of incident wave increases the structure vibration becomes more intense and structural vibration response is not only related to the incident wave period and wave height，and the wave form of also have certain relationship.

fixed offshore platform；time domain；vibration response

P751

A

1672 - 7619（2016）08 - 0052 - 05

10.3404/j.issn.1672 - 7619.2016.08.011

2016 - 02 - 02；

2016 - 03 - 15

浙江省自然科学基金资助项目（Y14E090034，Y13F020140）；浙江省青年科学家培养计划资助项目（2013R60G7160040）

谢卓（1990 - ），男，硕士研究生，研究方向为流体检测与仿真技术。