不同地震下单桩式风电机地震响应分析

（文末免费领取风电塔ABAQUS模型）

1.背景介绍

为了提高风能利用效率，增加单位时间内风力发电数量，海上风电塔向着叶轮更大、高度更高的趋势发展。如图1所示，2010年我国风电塔总高度在90m左右，到2020年，风电塔总高度上升至220m，预计到2025年，海上风电塔装机容量可达15~20MW，高度可达240m以上。

图1海上风电塔发展趋势

单桩式海上风电塔属于较为敏感的高耸薄壁结构，因其高度较高，支撑结构径厚比较大，因而具有较大的柔性。因此，在地震作用下极易发生风机停机或倒塌破坏。

图2 单桩式风电塔

    单桩式风电塔支撑结构为锥形塔筒，塔筒顶部直径3.07m，底部直径4.5m，顶部塔筒厚度20mm，底部塔筒厚度50mm。塔筒分上段塔筒、中段塔筒、下段塔筒三部分，材料为Q345B钢材，塔筒具体结构参数如上图所示。

2.风机有限元建模

采用ABAQUS建立风机整体有限元模型，结构与土相互作用采用假想嵌固点法。

图3嵌固点示意图

风电塔筒为锥形，分上、中、下三部分，在ABAQUS中利用旋转命令建立锥形塔筒模型，打开草图绘制界面，依次输入图示坐标点，点击完成，进行塔筒部件的创建。桩基础按照假想嵌固点法，截取12.346m钢管桩并将底部固定，材料为Q345钢材，网格采用C3D8R实体单元。

图4 塔筒部件建模

上部结构包括机舱、轮毂及叶片，机舱内主要为齿轮箱、发电机等，在建模时将机舱及轮毂简化为整体结构，且不考虑其塑性变形，弹性模量按照钢材输入206GPa，材料密度按照机舱总重进行等效机舱总质量为131.427t，等效密度为361kg/m³。叶片由玻璃复合纤维材料制成，计算时不考虑其塑性破坏，其弹性模量为28GPa，泊松比0.3，密度为66.18 kg/m³。叶片、轮毂及机舱之间采用“Tie”约束。网格划分时，叶片及轮毂采用C3D10单元，机舱采用C3D8R单元，风电塔整体有限元模型如图5所示。

图5风电塔整体结构有限元模型

3.地震加载

按照日本《道桥示方书》《建筑抗震设计规范》GB50011-2010的相关规定，选取El Centro、Nihonkai-Chubu、Hyogoken-Nanbu三种不同类型地震动作为数值模拟输入地震动。将所选地震记录调整为不同的PGA值(0.07g、0.2 g、0.40 g)，对于El Centro和Nihonkai-Chubu地震，仅考虑水平地震动输入，对于Hyogoken-Nanbu地震，考虑水平地震动和竖向地震动的双向输入，按照规范要求，竖向地震加速度峰值取水平峰值的0.65倍。

图6地震动加速度时程曲线

在分析步设置中创建重力分析步Step-1和动力隐式分析步Step-2，重力加载时间为1s，地震加载时间按照所选取的地震波时长进行更改。

图7 分析步设置

在ABAQUS中，取消塔筒底部X方向平动自由度约束，并在X方向添加加速度激励，加速度幅值曲线选取相应地震波。

图8 地震加载设置

4.结果对比

在El Centro地震作用下，当PGA=0.07g、0.2g和0.4g时，风电塔顶部位移最大值MDT（The Maximum Displacement of the Top）分别为0.051m， 0.122m和0.245m，地震峰值加速度增加分别导致结构MDT值增加136%和378%，且峰值出现时间不同；在Nihonkai-Chubu地震作用下，三种峰值加速度下结构的顶部位移分别为0.211m， 0.493m和0.989m；在Hydrogen地震作用下，结构MDT值分别为0.031m， 0.085m和0.170m；比较以上结果可以看出，El Centro地震作用下结构出现最大位移的时间存在变化，但是其他两类地震作用下结构最大位移值未发生变化。

图9风电塔结构顶部位移时程曲线

从上述结果中可以看出，近场地震和远场地震对风电塔这种高耸薄壁结构的影响差别较大。由于远场地震卓越频率与风电塔自振频率接近，因此在远场地震作用下的加速度及位移响应值更大，但近场地震PGA的增加对结构响应的变化幅度影响更大。

三种地震作用下PGA=2g时的结构破坏形态，如图10所示。由图可知，在El Centro及Hydrogen地震作用下，结构均表现为塔筒顶部弯曲过大导致上部结构倒塌，在Nihonkai地震作用下不仅塔筒顶部出现弯曲破坏，在塔筒与桩基础相交处也出现了屈曲和塑性变形。由此可知，远场地震作用对于风电塔结构的破坏性更大。

图10 不同地震作用下结构的倒塌形式

相关课程推荐

手把手教你ABAQUS不同方法的钢框架地震响应分析

ABAQUS螺栓连接T 型钢节点滞回模拟

钢结构S1 ABAQUS经典金属弹塑性本构及模拟应用

钢结构S2 ABAQUS常见金属循环本构对比及应用