支撑结构频率对海上风机动力响应的影响机理研究.pdf

第 41 卷第 4 期 2023 年 7 月Vol.41 No.4Jul.2023海 洋 工 程THE OCEAN ENGINEERING支撑结构频率对海上风机动力响应的影响机理研究李彦娥，李涛，彭驰，李强，张杰（中国船舶重工集团海装风电股份有限公司，重庆 401120）摘要：海上支撑结构的优化设计是海上风机技术发展的必然趋势，降低支撑结构的载荷是保证风机安全运行的有效途径。海上支撑结构受到风浪复杂环境荷载作用，风、浪载荷决定着塔底承受较大的剪力和倾覆力矩，同时风浪的随机性和周期性会影响塔架的疲劳载荷。基于海上风机支撑结构频率对载荷影响的研究，分析海上支撑结构频率对叶片根部挥舞和摆振响应的影响，探究频率对风机响应的影响机理。研究表明：波浪频率诱导是基础疲劳载荷响应的主要原因；开展单桩基础设计，当整机频率确定时，基础变径段可设置于浪溅区下部区域范围；叶片摆振响应受1P频率影响较大，在风机设计时可增加叶片摆振方向的阻尼；当整机频率邻近3P频率时，塔筒刚度变化对基础载荷响应的影响大于基础刚度变化；海上支撑结构设计时可优先考虑塔筒刚度。关键词：海上支撑结构；整机频率；疲劳载荷；极端载荷；波浪频率中图分类号：TK81 文献标志码：A DOI:10.16483/j.issn.1005-9865.2023.04.004Research on the influence mechanism of support structure frequency on the dynamic response of offshore wind turbineLI Yane，LI Tao，PENG Chi，LI Qiang，ZHANG Jie（China State Shipbuilding Corporation Haizhuang Wind Power Co.，Ltd.，Chongqing 401120，China）Abstract:The optimal design of the offshore support structure is an inevitable trend in the development of offshore wind turbine technology,and reducing the load of the support structure is an effective way to ensure the safe operation of the wind turbine.The offshore support structure is subjected to the complex environmental loads of wind and wave loads,which determines that the tower bottom bears a large shear force and overturning moment.Meanwhile,the randomness and periodicity of wind and wave will affect the fatigue load of the structure.The influence of offshore wind turbine support structure frequency on blade root flapwise and degewise response is analyzed based on the research on the influence mechanism of the frequency on the load of the support structure,and the influence mechanism of frequency on wind turbine response is explored.The research shows that the wave frequency induction is the main cause of foundation fatigue load response.When carrying out the monopile foundation design,as the machine natural frequency is determined,the foundation variant section can be set in the lower area of the splash zone.The blade edgewise response is greatly affected by the 1P frequency,and the damping of the blade fore-aft direction can be increased in the design of the wind turbine.As the natural frequency is close to the 3P frequency,the influence of the tower stiffness change on the load response is greater than that of the foundation.Tower stiffness can be given priority in the design of offshore support structures.Keywords:offshore support structure；machine natural frequency；fatigue load；extreme load；wave frequency海上风机的主要组成部分包括：叶片与轮毂、动力传动系统（齿轮箱）、发电机、电控系统以及支撑上部组件的塔筒。随着技术的不断更新和发展，质量更轻、成本更低是海上风机发展的基调。塔筒的质量占据文章编号：1005-9865（2023）04-0038-11收稿日期：2022-09-26作者简介：李彦娥（1993），女，工程师，主要从事土木及海洋工程结构方面的研究。E-mail：通信作者：李 涛。E-mail：第 4 期李彦娥，等：支撑结构频率对海上风机动力响应的影响机理研究了上部风机总质量的60%左右1，是设计的关键点。海上风机的成本比陆上风机高20%，塔架和基础的成本比陆上基础高350%2，基础结构成本通常占一个风电场总成本的20%30%3。海上风机结构受到风浪流耦合荷载作用，其支撑结构（塔筒+基础）受到较大倾覆力矩进而影响风机正常运行。塔筒和基础作为风机的支撑结构，其整体结构的优化设计是海上风机技术发展的必然趋势4-6，降低支撑结构的载荷是保证风机安全的有效途径。海上风机属于高度非线性和弹性系统，风机载荷计算采用考虑系统弹性响应的时域分析，非线性主要影响因素包括：空气动力、控制系统、波浪荷载以及桩土相互作用7-8。非线性响应进一步受到海上支撑结构高动态响应的影响。海上风机支撑结构的响应受到风谱、波浪谱以及空气动力阻尼的影响，空气动力阻尼又受到风机的空气动力和控制系统的影响9。海上风机支撑结构载荷计算通常通过模拟湍流风场和波谱的时域方法，采用叶素动量理论和莫里森方程计算风浪环境荷载作用下每个时间步长的结构响应10-11。吴加文等12基于Bladed软件研究海上风机单桩基础载荷影响因素，发现支撑结构的频率对基础载荷影响较大。孔德森等13研究了水平循环荷载作用下不同因素对桩身水平位移、剪力和弯矩的影响规律，表明不同荷载频率时桩身位移在零点以上变化较大，桩身弯矩随着频率的增加逐渐增大。海上风机支撑结构需要保证风电机组的正常运行外，还需确保整机频率满足风机允许运行频率范围，避免整机出现共振现象。风电系统风剪切塔影产生的转矩也将造成风机构件疲劳损伤14。Dolan和Lehn15对风剪切及塔影效应进行了分析，认为由此引起3P及其谐波扰动对系统影响最大。Zhang等16针对电功率的3P扰动，分别采用独立变桨控制及带通滤波器的转矩控制减小了功率波动。海上单桩基础相对于其他基础型式整机刚度较小，可能和波浪频率发生共振，导致结构动力效应放大。黄扬等17分析叶片挥舞与扭转变形对风力机气动载荷影响，李亮等18和楼文娟等19分析风机叶片挥舞和摆振影响机理。Klose等9开展了风浪流联合作用下海上风机基础极端荷载和疲劳荷载的时域分析，分析各节点处荷载响应，进而优化支撑结构。针对海上风机单桩基础支撑结构，分析支撑结构极端荷载和等效疲劳荷载的主要影响因素，研究支撑结构对基础载荷的影响机理，进一步分析海上支撑结构频率对叶片根部挥舞和摆振响应的影响，探究频率对风机响应的影响机理。1海上支撑结构模态分析法海上风机由叶轮、机舱、塔筒和基础组成，属于复杂的多自由度系统。在考虑阻尼时，海上风机的单自由度运动方程为20-21：x+2iixi+2ix=FiMi（1）式中：Mi=rotorm(r)2i(r)dr为模型质量矩，m(r)为模型质量；Fi=rotorf(r)2i(r)dr为模型的动力荷载；i(r)为模态振型；f(r)为风轮和塔架的分布荷载；x，x，x分别为实践变量的模型位移、速度和加速度；i为模型频率；i为模型阻尼系数。多自由度模型的耦合运动方程为：Mx+Cx+Kx=F（2）式中：M，C，K分别为模型的质量矩阵、阻尼矩阵和刚度矩阵。海上风机结构作用力的平衡方程为：Mcrr McrMcTr McDcr Dc DcTr DcT 0vc cc+0Cc0 c+0Kc0c=fca+fc0ca+c0ac2-fcici0（3）式中：vc为节点的速度矢量；c为应变矢量；c为拉格朗日乘子；fci和ci分别为非线性惯性力和应变矢量；fca为施加的节点力矢量；ca为与应变率对偶的应力矢量；fc0为节点反力；c0为约束应力；ac2为加速度变换的对流项；Mcrr，Mcr，Mc为质量矩阵分项；Cc为阻尼矩阵；Kc为刚度矩阵；Dcr，Dc为约束力矩阵。基于式（3）表达的平衡方程，结构荷载通过节点（如叶根、动力传动系统和塔基）处的反力计算。根据虚功原理，惯性力由每一处的质量特性和加速度矢量积分计算；加速度矢量包括模态分量、离心分量、科里奥39第 41 卷 海 洋 工 程利分量和重力分量。同时，考虑了由于结构变形引起的轴向力的二阶效应，即作用在形变结构上施加轴向力对弯矩的贡献。对于海上风机单桩基础结构，需要保证支撑结构在外部极端荷载作用下的极限承载能力，同时要保证风机在25年服役期间，在风浪循环荷载作用下风电结构不发生疲劳破坏。因此，研究主要针对极端荷载和等效疲劳荷载。2海上支撑结构载荷影响因素2.1环境荷载海上支撑结构受到风浪复杂环境荷载作用，风、浪载荷决定着塔底承受较大的剪力和倾覆力矩，同时风浪的随机性和周期性会影响塔架的疲劳载荷。风载作用于桨叶传递到塔架，波浪载荷作用于基础传递到塔架。基于IEC61400-3规范22和EN 1993-1-9规范23，采用载荷计算软件Bladed分析海上支撑结构在不同状态下极端载荷和疲劳载荷的影响因素。Bladed计算结构荷载包括风机正常运行状态、运行+故障状态、启停机状态、空转状态以及停机维护状态24。选取的 DLC6.2 工况为极端风浪耦合工况，风机处于空转状态。该工况下风机载荷较大，计算风速37.5 m/s，采用Kaimal湍流风；波浪采用Jonswap谱，谱峰升高因子取3.3，有效波高5.63 m，谱峰周期9.62 s；水深32 m。图1为DLC6.2工况下的风速和波浪的功率谱密度曲线。可知风的能量集中在0.1 Hz以下；波浪能量集中在波浪频率0.1 Hz 附近。表1为海上风机不同结构的频率，叶片转动范围考虑10%的裕度。图2为海上风机结构示意。图1风速谱和波浪谱Fig.1The wind speed spectrum and wave spectrum表1风机频率Tab.1The wind turbine frequency结构叶片转动1P范围叶片转动3P范围频率/Hz0.0810.1850.2430.556结构叶片一阶挥舞频率叶片一阶摆振频率频率/Hz0.5540.692图2海上风机结构示意Fig.2The structural diagram of offshore wind turbine40第 4 期李彦娥，等：支撑结构频率对海上风机动力响应的影响机理研究2.2基础结构对载荷的影响海上风电场在同一场区一般采用相同的塔筒结构，进而根据不同地质条件选择不同基础结构。研究支撑结构的塔筒部分保持不变，塔筒底部直径7 m，基础底部直径分别为7.0、7.5、8.0以及8.5 m时海上支撑结构不同位置处载荷变化，文中载荷均为弯矩载荷。图3为整机模型示意，为对比明显，将单桩基础变径段设置于浪溅区中上部区域。塔筒底部固定状态初始一阶频率为0.32 Hz。风电机组支撑结构设计时需要考虑强度及稳定性影响，支撑结构在环境荷载作用下不同位置处的极端荷载是其设计的必要条件。图4和图5为DLC6.2工况下不同基础底部直径时泥面处和底法兰处的极端荷载时程曲线。图6为不同位置处极端荷载最大值统计。随着基础底部直径增大，基础泥面处和底法兰处荷载最大值依次增大；基础底部直径7.0 m和8.5 m时泥面处极端荷载最大值分别为268 MN m和318 MN m，增大18.6%；底法兰处极端荷载最大值分别为138 MN m和146 MN m，增大5.8%。基础底部直径增大对泥面处荷载影响较大。图3整机模型Fig.3The holistic model图4泥面处极端荷载时程曲线Fig.4Time history curve of extreme load at mud surface图5法兰处极端载荷时程曲线Fig.5Time history curve of extreme load at flange41第 41 卷 海 洋 工 程图7和图8分别为DLC6.2工况下不同基础底部直径下泥面处和底法兰处极端荷载频谱。海上风电结构在运行状态会受到风诱导的低频响应、波浪频率、机组1P、3P、叶片挥舞及摆振频率以及固有频率共振的影响。可以看出，风频率对支撑结构泥面处和底法兰处极端荷载的影响相当，泥面处受到波浪荷载影响远大于底法兰处，且随着基础底部直径增大，基础受到波浪频率的影响增大，泥面处荷载增大。泥面处极端荷载受波浪频率诱导作用大于风频率。基础底部直径分别为7.0、7.5、8.0和8.5 m时整机结构的固有频率分别为0.22、0.225、0.23、0.234 Hz。根据IEC 61400-3规定，疲劳荷载包括正常运行状态、运行+故障状态、启停机状态以及维修状态。图8为不同基础底部直径时泥面处和底法兰处等效疲劳荷载。图9为等效疲劳荷载统计值。随着基础直径增大，基础受到波浪荷载作用较大，且对于维护工况，此时风机处于停机状态，波浪荷载效应远大于风荷载，因此基础泥面处和底法兰处的等效疲劳荷载均增大。单桩基础设计时增加直径会增加结构刚度和强度，从而可能降低结构应力，同时增加直径会增加波浪载荷，从而可能增加结构应力，在塔筒结构确定的条件下，波浪荷载增大引起的结构应力增加对疲劳荷载影响更大，桩基疲劳寿命减小，因此设计时在保证基础刚度的前提下可采用较小直径的基础。图10为DLC6.2工况下不同基础底部直径时风电机组叶片根部弯矩的统计特征值。由于DLC6.2工况为机组空转工况，叶片受到的非对称气动载荷小，且波浪荷载作用对风机机组荷载影响很小，因此不同基础底部直径时风机叶片根部弯矩最大值和标准差相差很小。图11为DLC6.2工况下不同基础底部直径时风电机组叶片根部弯矩频谱。叶片根部弯矩主要受到低频风荷载和叶片挥舞及摆振频率影响，波浪荷载对其影响甚微。不同基础底部直径下风荷载和叶片挥舞及摆振频率响应幅度基本相同，因此，不同基础底部直径下叶片根部弯矩变化较小。图6不同基础直径下极端荷载Fig.6Extreme load under different foundation diameters图7泥面处极端荷载频谱Fig.7Spectrum of extreme load at mud surface图8法兰处极端荷载频谱Fig.8Spectrum of extreme load at flange图9不同基础直径下等效疲劳荷载Fig.9Equivalent fatigue load under different foundation diameters42第 4 期李彦娥，等：支撑结构频率对海上风机动力响应的影响机理研究3支撑结构频率对载荷的影响3.1基础频率对载荷的影响风机频率是影响风机支撑结构荷载的重要因素。塔筒和基础是海上风机支撑结构的组成部分，其荷载受到塔筒和基础频率的影响。为避免机组共振响应，风电机组运行时风机整机频率应避开其1P和3P频率范围，海上支撑结构初始设计时可根据风资源、海洋水文参数以及1P和3P频率初步确定支撑结构的目标频率范围25。图12为海上风机支撑结构频率谱。本节首先针对同一风电场区塔筒结构保持不变，仅改变基础频率分析频率变化对基础荷载的影响。选用风电机组频率范围为0.1850.25 Hz之间，针对DLC6.2工况风机运行状态研究基础频率对基础极端荷载的影响。3.1.1极端荷载图13为DLC6.2工况不同基础频率下泥面处和底法兰处极端荷载统计。可以看出，当风机结构整机频率为0.188 Hz时，由于邻近1P频率，基础泥面处和底法兰处极端荷载最大；频率为0.205 Hz时，基础荷载达到最小值；继而随着频率增大，基础载荷出现增长趋势，当频率增大至0.220 Hz后，基础泥面处和底法兰处荷载出现微幅减小趋势。由于DLC6.2工况下风速频率主要在0.1 Hz以下，波浪荷载频率在0.1 Hz左右范围，因此，随着整机频率增大，风荷载和波浪荷载诱导作用对基础荷载影响减小，基础荷载主要受到风电机组共振频率影响。整机频率远离1P和3P范围边界值时，基础荷载较小。海上单桩基础荷载计算工况包括风机运行不同状态，极端荷载和疲劳荷载是海上支撑结构设计的关键因素。表2为不同计算工况下塔筒和基础频率。图14为全工况不同频率下单桩基础泥面处和底法兰处的极端荷载统计。由于风浪的随机性和周期性，不同频率下基础极端荷载最大值工况并不一定出现在风机的同一运行状态，因此，整体工况基础极端荷载和DLC6.2工况基础极端荷载随频率变化趋势不同。极端荷载最大值工况风机受到风浪荷载夹角均为30；基础频率在 0.188 Hz 时，基础响应受 1P 共振影响，泥面处和底法兰处极端荷载最大；基础频率增大至 0.220 Hz时，载荷随频率变化趋于平稳；工况五频率下基础泥面处频率比工况一减小12.7%。图10叶片根部弯矩统计值Fig.10Statistical value of bending moment at blade root图11叶片根部弯矩频谱Fig.11Spectrum of bending moment at blade root图12海上风机支撑结构频率谱Fig.12Frequency spectrum of offshore wind turbine support structure43第 41 卷 海 洋 工 程图15（a）为不同频率下风机动态响应最大运行状态叶片根部摆振位移统计特征值。工况一时叶片根部摆振位移远大于其他工况，工况一时叶片根部摆振位移最大幅度为2.7 m，其他工况最大幅度为0.55 m。图16为叶片根部摆振位移频谱。叶片根部摆振响应主要受固有频率和1P频率的影响，工况一下风机响应远大于其他工况。图15（b）为叶片根部挥舞位移统计特征值。随着风机频率增大，叶片根部挥舞位移标准差值逐渐减小，工况五下标准差值比工况一减小52%。图17为不同频率下风机动态响应最大运行状态时叶片根部摆振和挥舞加速度统计特征值。图18为其对应的频谱。工况一下叶片根部摆振加速度远大于其他工况，叶片根部挥舞加速度幅度随频率变化较小，表2塔筒和基础频率Tab.2The tower and foundation frequency单位：Hz工况塔筒底部固定频率整机频率工况一0.3240.188工况二0.3240.205工况三0.3240.220工况四0.3240.223工况五0.3240.238图13DLC6.2工况不同频率下极端荷载Fig.13Extreme load at different frequencies in DLC6.2 case图14全工况下不同频率下极端荷载Fig.14Extreme load at different frequencies under all cases图15叶片根部位移统计特征值Fig.15Statistical eigenvalue of blade root displacement图16叶片根部摆振位移频谱Fig.16Blade root edgewise displacement spectrum44第 4 期李彦娥，等：支撑结构频率对海上风机动力响应的影响机理研究标准差值随频率增大逐渐减小。叶片根部加速度响应受到1P、叶片挥舞和摆振频率及结构固有频率的诱导作用。叶片荷载在整个风力机系统中占很大的比率，其对于塔架底部弯矩的影响显著。风力机叶片是气动敏感结构，在进行风机降载设计时可增加叶片的结构阻尼，特别是摆振方向的阻尼。3.1.2疲劳荷载图19为不同频率下泥面处和法兰处的等效疲劳荷载统计。泥面处和法兰处的等效疲劳荷载受风机1P和3P影响较小，因此随着风机频率增大，基础等效疲劳荷载出现增长趋势；当整机频率增大至0.238 Hz时，基础等效疲劳荷载呈现微幅减小趋势。海上风机支撑结构等效疲劳荷载受风浪频率诱导影响较大，疲劳荷载工况主要为DLC1.2、DLC6.4和DLC8.3工况，DLC8.3工况为安装工况正常湍流风浪耦合工况，计算风速由326.25 m/s，波浪采用JONSWAP谱，最大有效波高5.92 m，谱峰周期7.71 s，此工况为未接入电网工况，波浪荷载为主要影响因素。图20为支撑结构泥面处和底法兰处的弯矩荷载范围谱。累计循环次数较小时，支撑结构弯矩荷载范围边界随累计循环次数增大呈线性减小趋势，继而缓慢降低。由于波浪疲劳荷载以惯性诱导为主，可以假设图17叶片根部加速度统计特征值Fig.17Statistical eigenvalue of blade root acceleration图18叶片根部加速度频谱Fig.18Blade root acceleration spectrum图19不同频率下等效疲劳荷载Fig.19Equivalent fatigue load at different frequencies45第 41 卷 海 洋 工 程相应的波浪荷载与响应成正比，支撑结构波浪荷载是疲劳损伤的主要因素。3.1.3基础结构优化由于塔筒结构保持不变，因此基础顶部直径确定，当基础底部直径大于顶部直径时，基础存在变直径段，当变直径段在浪溅区下部区域及水下区时，由于基础受到波浪荷载较小，两种结构形式下基础泥面处和底法兰处的荷载相当。图21为整机频率分别为0.220、0.225、0.230以及0.234 Hz下不同桩径基础极端荷载对比。图中实线代表不同频率下单桩基础底部直径7.5 m，变径段在浪溅区中下部区域极端荷载；虚线代表不同频率下单桩基础底部直径7.08.5 m，变径段在浪溅区中上部区域极端荷载。可以看出，在该整机频率范围内，同一基础结构泥面处和底法兰处的极端荷载随着频率增大而减小；相同频率下，当直径较大时，由于波浪荷载作用较大，因此泥面处和底法兰处的荷载较大，且直径增大对基础泥面处荷载影响程度大于底法兰处。因此，在开展基础设计，当整机频率确定时，基础变径段可设置于浪溅区下部区域范围。3.2塔筒频率对载荷的影响保证基础结构不变，改变塔筒刚度进而分析塔筒结构频率对整机荷载的影响。表3为不同计算工况下塔筒和基础频率。图22和图23分别为不同塔筒频率下基础泥面和法兰处等效疲劳荷载对比，当整机频率大于0.205 Hz后，泥面处极端荷载和等效疲劳荷载均随着频率增大呈现增长趋势。而图14和图19表明仅通过改变基础刚度工况下整机频率达到0.223 Hz后，基础泥面处和法兰处极端荷载和等效疲劳荷载均出现微幅减小趋势。图24和图25分别为工况五和工况十下支撑结构泥面处和法兰处极端荷载频谱。两者的整机频率一致，工况十下风频率和3P频率诱导作用略大于工况五，进而基础固有频率响应幅度大于工况五。因此，当整机频率接近3P频率，同一频率下仅改变塔筒刚度较仅改变基础刚度对基础荷载的响应幅度明显。图20支撑结构弯矩荷载范围谱Fig.20Bending moment load spectra for support structure图21不同频率下不同桩径基础极端荷载Fig.21Extreme loads of foundations with different pile diameters at different frequencyt46第 4 期李彦娥，等：支撑结构频率对海上风机动力响应的影响机理研究4结 语文中研究的主要结论如下：1）基于同一塔筒结构，随着基础直径增大，基础受到波浪荷载作用较大，且对于维护工况，此时风机处于停机状态，波浪荷载效应远大于风荷载，因此基础泥面处和底法兰处的等效疲劳荷载均增大。2）叶片根部响应受到1P、叶片挥舞和摆振频率及结构固有频率的诱导作用，叶片摆振响应受1P频率影响较大，且风力机叶片是气动敏感结构，在进行设计时可增加叶片的结构阻尼，特别是摆振方向的阻尼。3）在开展基础设计时，当整机频率确定时，基础变径段可设置于浪溅区下部区域范围。4）当整机频率邻近3P频率时，塔筒刚度变化对基础载荷响应的影响大于基础刚度变化，同一频率范围基础和塔筒刚度改变对载荷的影响程度不同，海上支撑结构设计时可优先考虑塔筒刚度。参考文献：1 ANCONA D，MCVEIGH J.Wind turbine-materials and manufacturing fact sheetJ.Princeton Energy Resources International，LLC，2001：1-8.2 张浦阳，李彦娥，丁红岩，等.海上风电基础下放通过浪溅区响应特性研究 J.海洋工程，2022，40（4）：77-87.（ZHANG P Y，LI Y E，DING H Y，et al.Response analysis of lowering through the splash zone operation for offshore wind turbine foundation J.The Ocean Engineering，2022，40（4）：77-87.（in Chinese）3 GASCH R.TWELE J.Wind power plants：fundamentals，design，construction and operationM.Berlin，Heidelberg：Springer-Verlag，2011.表3塔筒和基础频率Tab.3The tower and foundation frequency单位：Hz工况塔筒底部固定频率整机频率工况六0.2360.188工况七0.2730.205工况八0.3120.220工况九0.3240.223工况十0.3710.238图22极端荷载统计Fig.22Extreme load statistics图23等效疲劳荷载统计Fig.23Equivalent fatigue load statistics图24泥面处极端荷载频谱Fig.24Spectrum of extreme load at mud surface图25法兰处极端荷载频谱Fig.25Spectrum of extreme load at flange47第 41 卷 海 洋 工 程4 翟恩地，张新刚，李荣富.海上风电机组塔架基础一体化设计 J.南方能源建设，2018，5（2）：1-7.（ZHAI E D，ZHANG X G，LI R F.Integrated design of offshore wind tower and foundation J.Southern Energy Construction，2018，5（2）：1-7.（in Chinese）5 唐浩渊.海上风电机组单桩与导管架基础结构一体化分析及选型方法研究 D.重庆：重庆大学，2021.（TANG H Y.Research on integrated analysis and selection method of monopile and jacket substructure for offshore wind turbineD.Chongqing：Chongqing University，2021.（in Chinese）6 张博.海上风机一体化载荷仿真方法研究 D.上海：上海交通大学，2016.（ZHANG B.Study on integrated loads simulation of offshore wind turbine D.Shanghai：Shanghai Jiao Tong University，2016.（in Chinese）7 Petroleum and natural gas industries-Fixed steel offshore structures：ISO 19902：2020S.International Organization for Standardization，2020.8 孟祥然，刘润，练继建，等.海上大直径钢管桩打桩过程中桩周土强度弱化研究 J.海洋工程，2022，40（2）：88-96.（MENG X R，LIU R，LIAN J J，et al.Study on the weakening of soil resistance around large diameter steel pipe pile during pile driving J.The Ocean Engineering，2022，40（2）：88-96.（in Chinese）9 KLOSE M，DALHOFF P，ARGYRIADIS K.Integrated load and strength analysis for offshore wind turbines with jacket structures C/Proceedings of the 17th ISOPE.2007.10DNV.Environmental conditions and environmental loads：DNV-RP-C205 S.Norway：Det Norske Veritas，2010.11韩彦青.海上潜式张力腿平台浮式风机耦合动力特性研究 D.天津：天津大学，2016.（HAN Y Q.Coupled dynamic analysis of a submerged tension leg platform wind turbine D.Tianjin：Tianjin University，2016.（in Chinese）12吴加文，冀卫东，王爱国.海上风电单桩基础载荷影响因素研究 J.水电与新能源，2020，34（3）：17-22.（WU J W，JI W D，WANG A G.Influencing factors on monopile foundation load in offshore wind farms J.Hydropower and New Energy，2020，34（3）：17-22.（in Chinese）13孔德森，刘一，邓美旭，等.海上风电单桩基础土相互作用特性影响因素分析 J.海洋工程，2021，39（1）：100-111.（KONG D S，LIU Y，DENG M X，et al.Analysis of influencing factors of monopile foundation-soil interaction characteristics for offshore wind power J.The Ocean Engineering，2021，39（1）：100-111.（in Chinese）14张兴，滕飞，谢震，等.风力机的风剪塔影模拟及功率脉动抑制 J.中国电机工程学报，2014，34（36）：6506-6514.（ZHANG X，TENG F，XIE Z，et al.Simulation and power pulsation suppression of wind turbines with wind shear and tower shadow effects J.Proceedings of the CSEE，2014，34（36）：6506-6514.（in Chinese）15DOLAN D，LEHN P.Simulation model of wind turbine 3p torque oscillations due to wind shear and tower shadow J.IEEE Transaction on Energy Conversion，2006，21（3）：717-724.16ZHANG Y，HU W，CHEN Z.Flicker mitigation strategy for a doubly fed induction generator by torque controlJ.IET Renewable Power Generation，2014，8（2）：91-99.17黄扬，赵伟文，万德成.浮式风机气动水动气弹性耦合响应数值模拟 J.海洋工程，2022，40（4）：88-101.（HUANG Y，ZHAO W W，WAN D C.Numerical simulation of coupled aero-hydro-elastic performance of floating offshore wind turbine J.The Ocean Engineering，2022，40（4）：88-101.（in Chinese）18李亮，李映辉，杨鄂川.风力机叶片挥舞摆振气弹稳定性分析 J.噪声与振动控制，2015，35（5）：30-34.（LI L，LI Y H，YANG E C.Flap/lead-lag aeroelastic stability analysis for wind turbine blades J.Noise and Vibration Control，2015，35（5）：30-34.（in Chinese）19楼文娟，余江，潘小涛.风力机叶片挥舞摆振气弹失稳分析 J.工程力学，2015，32（11）：236-242.（LOU W J，YU J，PAN X T.Calculating for aerodynamic stability response of wind turbine blade in flapwise and edgewise directionJ.Engineering Mechanics，2015，32（11）：236-242.（in Chinese）20BOSSANYI A.Bladed theory manual R.Garrad Hassan，2006.21吴智深，李素贞，ADEWUYI A P.基于应变分布响应的模态分析理论与应用 J.科技导报，2010，28（8）：94-103.（WU Z S，LI S Z，ADEWUYI A P.Modal analysis based on distributed strain me