自升式钻井平台桩腿安装工况风载荷随机振动分析

1、自升式平台桩腿坞外接长过程风载荷频域分析陈峰 , 马骏(大连理工大学）摘 要：本文主要研究三桩腿桁架式自升式钻井平台在船厂进行桩腿连接过程脉动风载荷作用下随机振动分析。应用ANSYS软件功率谱分析（PSD）功能分析随着脉动风速的增大，平台本体升降不同高度时，平台振动响应的变化情况。分析表明，随着平台本体爬升和风速的增大，平台受脉动风载荷的影响增大，影响平台总装过程的安全性。关键词：自升式 脉动风载荷 加速度均方根 功率谱 ANSYS 自升式钻井平台在桩腿连接过程中受动风载荷影响较大，需要进行随机振动分析。结构的风振响应分析分为确定性和非确定性两类。确定性分析是指结构在确定的风荷载作用下的响应分

2、析；非确定性分析是考虑风荷载的随机性，研究自升式钻井平台在脉动风荷载作用下的随机振动分析。研究风作用下高耸钢结构动力响应的理论方法有两种，第一种是频率域法（frequency domain)，它按随机振动理论，建立了输入风荷载频谱的特性与输出结构响应之间的直接关系。第二种是时间域法(Time domain)，它是基于将随机的风荷载模拟成时间函数，然后直接求解运动微分方程。本文主要是通过第一种方法研究自升式钻井平台的振动响应。1 结构上的随机风荷裁的描述由于风和建筑物的相互作用引起多种气流，使风成为一极复杂的现象。但在风工程中作了一定的简化，使得可籍着识别下列特性得到风特性的有意义的预测：风速随

3、高度而变化；风的紊流性质；概率假设；风结构相互作用的动力性质。1.1风的脉动性质对结构工程而言，可以认为风速由两部分组成：随高度增大的平均速度部分和脉动风速部分。也可形象地认为平均风速随高度增大并作用有随机方向的涡流(在空间转动的小气流)。不同的涡流尺寸和转速有相应的不同波长和频率的紊流分量。由以上风的脉动性质，可得任意瞬时的速度 （1）即为平均风速和脉动风速，结构振动速度之和。1.2脉动风速功率谱脉动风的风速谱密度模型主要有Davenport谱、Kaimal谱、Vickery谱等，本文选择Davenport谱，其表达形式为： （2）式中：K为地面粗糙度系数，L为湍流积分尺度，为高度10m处的

4、平均风速。2 风振随机响应分析方法由连续体系微元的动力平衡关系，可以推导出连续体系动力学方程： （3）式中为某一高度处的水平位移，为沿高度方向的线质量，为阻尼，E为弹性模量，为截面惯性矩。由随机振动的模态分解理论，可得位移响应的谱密度： (4)式中：第i模态的频响函数（9），为i、j两模态广义力的互谱密度，*表示复共轭根据(5)(6)可得到位移响应的均方根1： (5)加速度的均方根为 （6）根据加速度响应的均方根值可以求出相应的风振力2 （7）其中为第i节点的集中质量，为峰值因子，一般取2-3，故结构在脉动风荷载激励下的风振系数为: （8）式中，分别为风振力和静风力。根据统计理论，结构在平均风

5、和脉动风作用下的总位移2为: （9）式中为平均风引起的位移，为位移均方根，为峰值保证因子，本章取2.5。3 自升式钻井平台抗风分析3.1 风振响应计算参数确定我国建筑结构荷载规范3GBJ987是从原GBJ974把地面粗糙度分为海洋和陆地二类修订为海洋(A)、空旷平坦地区(B)和大城市中心(C)三类，其中c类是根据当时中国大多数城市为依据而制定的，一般偏于安全。由于近年来，中国城市得到了飞速发展，因此GBJ9 87的修订稿中，将地面粗糙度增加为四类，除A、B二类不动外，调整后的C类是适应中等城市，D类适应大城市中心网。对于这四类地貌，其相应的地面粗糙度系数和梯度风高度,如表1所示： 表 地面粗糙

6、度对照表Tab.1roughness of terrene地貌分类 A B C D 0.12 0.16 0.2 0.3 300 350 400 450 地面粗糙度类别取类，Davenport脉动风速谱的公式中的地貌系数K为0.03，L为湍流积分尺度取为常数1200m。3.2 自升式平台平台本体高度对平台风振响应的影响本文以JU2000自升式钻井平台为例，研究自升式钻井平台平台本体在顶部和底部时 ，如图1所示，桩腿不同高度处的节点，顺风向位移、加速度均方根、风振系数以及横风向的加速度均方根值的变化趋势，从而全面把握结构的抗风性能。计算时考虑十阶振型参与计算，选用Davenport谱，不考虑气动阻

7、尼，风攻角为0度，即10m高处平均风速为20m/s。根据式（3），得到风载加速度功率谱表，在平台高度151米和48米两种工况下建立有限元模型，如图1, 应用ANSYS功率谱分析功能，进行随机振动（PSD）分析。 图1 自升式平台风载工况模型 Fig1 the jack-up modal with wind loading1）平台升至151米平台本体上的节点（点5）顺风向风振响应功率谱，如图2，3 图2 顺风方向位移响应功率谱 图3 顺风方向加速度响应功率谱Fig2 The displacement PSD spectrum of tail wind Fig3 The acceleration

8、PSD spectrum of tail wind 图4 横风方向位移响应功率谱 图5 横风方向加速度响应功率谱Fig2 The displacement PSD spectrum of beam wind Fig3 The acceleration PSD spectrum of beam wind2）平台降至48米平台本体上的节点（点2）顺风向风振响应功率谱 图6 顺风方向位移响应功率谱 图7 顺风方向加速度响应功率谱Fig4 The displacement PSD spectrum of tail wind Fig5 The acceleration PSD spectrum of t

9、ail wind 图8 横风方向位移响应功率谱 图9 横风方向加速度响应功率谱Fig6 The displacement PSD spectrum of beam wind Fig7 The acceleration PSD spectrum of tail wind利用ANSYS程序得到图1中桩腿不同高度处顺风、横风的加速度谱功率及位移功率谱，应用式（13）计算桩腿不同高度出顺风加速度均方根、横风加速度均方根、顺风位移均方根如表2，3，平台不同高度处的上述指标见图8图10。 表2平台升至装腿上部风振响应指标值Tab.2 wind vibration response when the mai

10、n platform raising up风振响应指标12 3 4 5 6 顺风加速度均方根0.119 0.322 0.569 0.757 0.820 0.824 横风加速度均方根0.0113 0.0306 0.0545 0.0727 0.0786 0.0787 顺风位移均方根0.951.131.421.621.892.03加速度均方根单位0.01g， 位移均方根单位cm表3平台降至装腿底部风振响应指标Tab.3 wind vibration response when the main platform down to the bottom风振响应指标1 2 3 4 5 6 顺风加速度均方根

11、0.075 0.142 0.206 0.318 0.401 0.505横加速度均方根3.57 3.33 2.86 2.50 2.00 1.67 顺风位移均方根0.007 0.013 0.022 0.037 0.046 0.058 加速度均方根单位0.01g， 位移均方根单位cm 图10顺风加速度均方根 图11横风加速度均方根 图12 顺风位移均方根Fig8 acceleration mean square Fig9 acceleration mean square Fig10 displacement mean squareroot of tail wind root of beam wind

12、 root of tail wind分析上面的图表可知，桩腿上节点，高度越高，风振力越大；随着平台本体高度的爬升，顺风向的风振力逐渐增加，顺风向位移逐渐加大；当平台本体高度较低时，横风向的风振力很小，工程上可以不考虑，随着其爬升高度增加，横风向加速度显著增加，在进行平台抗风计算时要考虑其对平台安全性的影响。3.3 风速对自升式平台风振响应的影响当平台升至安装过程的最大高度151米，脉动风平均分速分别取20m/s、25m/s、25m/s时，装腿不同高度处的节点，顺风向位移、加速度均方根、横风风向位移、加速度均方根见表4、5，以及各风振指标的变化趋势见图11-14。表4风速取20m/s时风振响应指

13、标值Tab4 wind vibration response at the speed of 20m/s风振响应指标1 2 3 4 5 6 顺风加速度均方根0.119 0.322 0.569 0.757 0.820 0.824 顺风位移均方根0.093 0.252 0.449 0.600 0.650 0.653 横风加速度均方根0.011 0.031 0.055 0.073 0.079 0.079 横风位移均方根0.009 0.024 0.043 0.057 0.061 0.061 表5风速取25m/s时风振响应指标值Tab5 wind vibration response at the sp

14、eed of 25m/s风振响应指标123456顺风加速度均方根0.180 0.484 0.854 1.135 1.230 1.235 顺风位移均方根0.138 0.375 0.670 0.895 0.971 0.975 横风加速度均方根0.017 0.046 0.082 0.109 0.118 0.118 横风位移均方根0.013 0.036 0.064 0.085 0.092 0.092 表6风速取30m/s时风振响应指标Tab6 wind vibration response at the speed of 30m/s风振响应指标1 2 3 4 5 6 顺风加速度均方根0.216 0.5

15、81 1.025 1.361 1.475 1.481 顺风位移均方根0.166 0.450 0.803 1.074 1.164 1.169 横风加速度均方根0.020 0.055 0.098 0.131 0.141 0.142 横风位移均方根0.016 0.043 0.077 0.102 0.110 0.111 图13 顺风加速度均方根 图14 顺风加位移均方根Fig11 acceleration mean square root of tail wind Fig12 displacement mean square root of tail wind 图15 横风加速度均方根 图16 横风加

16、位移均方根Fig11 acceleration mean square root of beam wind Fig12 displacement mean square root of beam wind通过比较，发现结构顺、横风向加速度均方根及顺风向位移均方根随风速的变化而显著变化，风速大则风振响应大，但顺风向风振系数却变化不大，这是因为动力风响应加大的同时静力风也随风速加大而加大，因此其比值风振系数变化不大，三种情况的平台本体上风振系数如下表所示：表7 不同风速下平台本体顺风风振系数Tab.7 vibration coeficient of plateform in different wind velocity风速20m/s25m/s30m/s风振系数（）1.721.691.714 结论(1) 自升式钻井平台桩腿上节点，高度越高，风振力越大。随着平台本体高度的爬升，顺