兆瓦级风力发电复合材料叶片疲劳寿命的预测方法技术

本发明专利技术公开了一种兆瓦级风力发电复合材料叶片疲劳寿命的预测方法，其包括：建立复合材料温度场计算模型；开发叶片‑风速有限元耦合动力学模型；开发风速‑裂纹‑缺口的相互作用模型；对上述模型中的应力时程曲线相加，得到组合载荷下叶片各部的应力时程曲线；用雨流计数法对应力时程曲线进行计算，得到组合载荷下叶片各部二维疲劳应力谱；建立等幅疲劳中值S

Prediction method of fatigue life of composite blades for megawatt wind power generation

【技术实现步骤摘要】
兆瓦级风力发电复合材料叶片疲劳寿命的预测方法
本专利技术涉及风力发电机叶片疲劳寿命预测领域。更具体地说，本专利技术涉及一种兆瓦级风力发电复合材料叶片疲劳寿命的预测方法。
技术介绍
现有的风力机叶片疲劳寿命的预测，主要通过测量叶片表面风速，计算叶片所受载荷，或者在风力机上加装传感器来获得疲劳载荷，确定应力循环次数，最后得到疲劳寿命。然而，该预测方法未对疲劳裂纹和外物损伤缺口所带来的破坏性进行考虑，致使最终预测的疲劳寿命存在较大误差，不符合实际生产的需要，降低了叶片疲劳寿命预测在工程应用中的价值。
技术实现思路
本专利技术的一个目的是解决至少上述问题和/或缺陷，并提供至少后面将说明的优点。本专利技术还有一个目的是提供一种兆瓦级风力发电复合材料叶片疲劳寿命的预测方法，其通过综合考虑风速中的风速、叶片的自重、叶片的疲劳裂纹以及外物损伤缺口的相互作用，减小复合材料叶片疲劳寿命的预测误差，提高其预测结果在实际生产中应用价值。为了实现本专利技术的这些目的和其它优点，提供了一种兆瓦级风力发电复合材料叶片疲劳寿命的预测方法，其包括：步骤一、根据气象局数据，结合复合材料温度场计算模型进行叶片温度场的计算，得到时间-叶片温度及温度梯度荷载时程曲线；步骤二、开发叶片-风速-风向有限元耦合动力学模型，采用预处理共轭梯度法求解耦合系统振动方程，在微机上进行叶片-风速-风向耦合系统动力仿真计算，结合步骤一中时间-叶片温度及温度梯度荷载时程曲线，得到时间梯度温度荷载下叶片各部的应力时程曲线；>步骤三、开发考虑叶片疲劳裂纹、外物损伤缺口与风速相互作用的风速-裂纹-缺口的相互作用模型，用所述风速-裂纹-缺口的相互作用模型进行计算风速作用下叶片应力时程曲线计算，结合步骤一中时间-叶片温度及温度梯度荷载时程曲线的计算结果，得到风速-裂纹-缺口作用下叶片上各部的应力时程曲线；步骤四、用叶片-风速-风向有限元耦合动力学模型进行叶片自重荷载下叶片各部应力时程曲线的计算，得到叶片自重载荷下叶片各部应力时程曲线；步骤五、步骤二、步骤三、步骤四中的应力时程曲线相加，得到组合载荷下叶片各部的应力时程曲线；步骤六、用雨流计数法对步骤五中的应力时程曲线进行计算，得到组合载荷下叶片各部二维疲劳应力谱；建立等幅疲劳中值Sn-Sm-N曲面，根据Miner线性疲劳累积损伤准则，基于二维疲劳应力谱对叶片的疲劳寿命进行估算。优选的是，从气象局网站下载包括年、月、日、平均风速、平均气温、日最低气温、日最高气温、日最大风速；借鉴工程结构物温度场研究成果，建立复合材料温度场计算模型。优选的是，所述疲劳裂纹的预测方法为：获取复合材料叶片的三维CAD模型和材料参数，确定疲劳裂纹的位置、裂纹的深度和裂纹的宽度，对健康叶片的三维CAD模型进行切割和剖分，对裂纹叶片进行网格划分，构建出含真实裂纹扩展路径的叶片有限元模型，对裂纹叶片进行加速转动过程中的振动相应分析，获取裂纹裂片在叶尖处尾缘的位移振动响应以及裂纹面上的接触压力响应方程；建立健康叶片的有限元模型，对健康叶片进行加速转动过程中的振动相应分析，获取裂纹裂片在叶尖处尾缘的位移振动响应以及裂纹面上的接触压力响应方程，将裂纹叶片和健康叶片在加速过程中的位移振动响应进行比较，确定叶片疲劳裂纹的位置、宽度、深度。优选的是，所述外物损伤缺口的预测方法为：对运行多年的复合材料叶片进行扫描分析，获得外物损伤缺口的形貌特征以及几何尺寸，对外物损伤缺口的形貌特征以及几何尺寸进行时间均值化，构建时间-风速-外物损伤缺口模型，根据时间-风速-外物损伤缺口模型采用叠加缺口进行计算，获得健康叶片在风速中运行产生的外物损伤缺口。优选的是，所述疲劳裂纹包括叶背面疲劳裂纹和叶片前缘疲劳裂纹。优选的是，所述的剖分具体为分隔成四边形并利用四边形网格进行划分。优选的是，步骤二、步骤三采用APDL二次开发语言，ANSYS通用有限元软件环境条件下进行开发。本专利技术至少包括以下有益效果：本专利技术综合考虑了风速、风向、叶片自重、疲劳产生的裂纹以及外物损伤缺口的综合作用，考虑的载荷较多、较完善，大大提高了复合材料叶片疲劳应力谱模拟精度，提高了复合材料叶片的疲劳寿命预测精度，使其具有更好的工程应用价值；基于二维疲劳应力谱的疲劳寿命预测方法，充分考虑到应力谱的均值与幅值的共同影响，理论上更为合理、完善，使预测结果更加具备实际应用价值。本专利技术的其它优点、目标和特征将部分通过下面的说明体现，部分还将通过对本专利技术的研究和实践而为本领域的技术人员所理解。具体实施方式下面对本专利技术做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。应当理解，本文所使用的诸如“具有”、“包含”以及“包括”术语并不配出一个或多个其它元件或其组合的存在或添加。实施例一种兆瓦级风力发电复合材料叶片疲劳寿命的预测方法，其包括：步骤一、根据气象局数据，结合复合材料温度场计算模型进行叶片温度场的计算，得到时间-叶片温度及温度梯度荷载时程曲线；步骤二、开发叶片-风速-风向有限元耦合动力学模型，采用预处理共轭梯度法求解耦合系统振动方程，在微机上进行叶片-风速-风向耦合系统动力仿真计算，结合步骤一中时间-叶片温度及温度梯度荷载时程曲线，得到时间梯度温度荷载下叶片各部的应力时程曲线；步骤三、开发考虑叶片疲劳裂纹、外物损伤缺口与风速相互作用的风速-裂纹-缺口的相互作用模型，用所述风速-裂纹-缺口的相互作用模型进行计算风速作用下叶片应力时程曲线计算，结合步骤一中时间-叶片温度及温度梯度荷载时程曲线的计算结果，得到风速-裂纹-缺口作用下叶片上各部的应力时程曲线；步骤四、用叶片-风速-风向有限元耦合动力学模型进行叶片自重荷载下叶片各部应力时程曲线的计算，得到叶片自重载荷下叶片各部应力时程曲线；步骤五、步骤二、步骤三、步骤四中的应力时程曲线相加，得到组合载荷下叶片各部的应力时程曲线；步骤六、用雨流计数法对步骤五中的应力时程曲线进行计算，得到组合载荷下叶片各部二维疲劳应力谱；建立等幅疲劳中值Sn-Sm-N曲面，根据Miner线性疲劳累积损伤准则，基于二维疲劳应力谱对叶片的疲劳寿命进行估算。在上述实施例的基础上，在一个实施例中，从气象局网站下载包括年、月、日、平均风速、平均气温、日最低气温、日最高气温、日最大风速；借鉴工程结构物温度场研究成果，建立复合材料温度场计算模型。温度的变化会使复合材料叶片内部产生应力，使得复合材料叶片的力学性能受到的影响，在叶片疲劳寿命预测时考虑温度变化，有利于提高叶片疲劳寿命预测的精确度。在上述实施例的基础上，在一个实施例中，所述疲劳裂纹的预测方法为：获取复合材料叶片的三维CAD模型和材料参数，确定疲劳裂纹的位置、裂纹的深度和裂纹的宽度，对健康叶片的三维CAD模型进行切割和剖分，对裂纹叶片进行网格划分，构建出含真实裂纹扩展路径的叶片有限元模型，对裂纹叶片进行加速转动过程中的振动相应分析，获取裂纹裂片在叶尖处尾缘的本文档来自技高网...

【技术保护点】
1.一种兆瓦级风力发电复合材料叶片疲劳寿命的预测方法，其特征在于，包括：/n步骤一、根据气象局数据，结合复合材料温度场计算模型进行叶片温度场的计算，得到时间-叶片温度及温度梯度荷载时程曲线；/n步骤二、开发叶片-风速-风向有限元耦合动力学模型，采用预处理共轭梯度法求解耦合系统振动方程，在微机上进行叶片-风速-风向耦合系统动力仿真计算，结合步骤一中时间-叶片温度及温度梯度荷载时程曲线，得到时间梯度温度荷载下叶片各部的应力时程曲线；/n步骤三、开发考虑叶片疲劳裂纹、外物损伤缺口与风速相互作用的风速-裂纹-缺口的相互作用模型，用所述风速-裂纹-缺口的相互作用模型进行风速作用下叶片各部的应力时程曲线计算，结合步骤一中时间-叶片温度及温度梯度荷载时程曲线的计算结果，得到风速-裂纹-缺口作用下叶片上各部的应力时程曲线；/n步骤四、用叶片-风速-风向有限元耦合动力学模型进行叶片自重荷载下叶片各部应力时程曲线的计算，得到叶片自重载荷下叶片各部应力时程曲线；/n步骤五、步骤二，步骤三，步骤四中的应力时程曲线相加，得到组合载荷下叶片各部的应力时程曲线；/n步骤六、用雨流计数法对步骤五中的应力时程曲线进行计算，得到组合载荷下叶片各部二维疲劳应力谱；建立等幅疲劳中值S...

【技术特征摘要】
1.一种兆瓦级风力发电复合材料叶片疲劳寿命的预测方法，其特征在于，包括：
步骤一、根据气象局数据，结合复合材料温度场计算模型进行叶片温度场的计算，得到时间-叶片温度及温度梯度荷载时程曲线；
步骤二、开发叶片-风速-风向有限元耦合动力学模型，采用预处理共轭梯度法求解耦合系统振动方程，在微机上进行叶片-风速-风向耦合系统动力仿真计算，结合步骤一中时间-叶片温度及温度梯度荷载时程曲线，得到时间梯度温度荷载下叶片各部的应力时程曲线；
步骤三、开发考虑叶片疲劳裂纹、外物损伤缺口与风速相互作用的风速-裂纹-缺口的相互作用模型，用所述风速-裂纹-缺口的相互作用模型进行风速作用下叶片各部的应力时程曲线计算，结合步骤一中时间-叶片温度及温度梯度荷载时程曲线的计算结果，得到风速-裂纹-缺口作用下叶片上各部的应力时程曲线；
步骤四、用叶片-风速-风向有限元耦合动力学模型进行叶片自重荷载下叶片各部应力时程曲线的计算，得到叶片自重载荷下叶片各部应力时程曲线；
步骤五、步骤二，步骤三，步骤四中的应力时程曲线相加，得到组合载荷下叶片各部的应力时程曲线；
步骤六、用雨流计数法对步骤五中的应力时程曲线进行计算，得到组合载荷下叶片各部二维疲劳应力谱；建立等幅疲劳中值Sn-Sm-N曲面，根据Miner线性疲劳累积损伤准则，基于二维疲劳应力谱对叶片的疲劳寿命进行估算。


2.根据权利要求1所述的兆瓦级风力发电复合材料叶片疲劳寿命的预测方法，其特征在于，从气象局网站下载包括上一年度的日平均风速、日平均气温、日最低气温、日最高气温、日最大风速；借鉴工程结构物温度场研究成果，建立复合材料温度场计算模型。


3.根据权利要求1所述的兆瓦级风力发电复合材料...

【专利技术属性】
技术研发人员：杨林，张国勇，
申请(专利权)人：通标标准技术服务天津有限公司，
类型：发明
国别省市：天津;12