一种风电叶片主梁褶皱损伤失效分析方法

本发明属于风电叶片复合材料力学性能测试与分析，具体涉及一种基于机理与静动态实验数据融合的风电叶片主梁褶皱损伤失效分析方法。

背景技术：

1、随着全球可再生能源产业的蓬勃发展，风力发电装机容量持续攀升。为降低发电成本，风电机组单机容量不断增大，叶片尺寸已突破百米级。叶片主梁作为主要承受弯曲与拉压交变载荷的核心骨架，其结构安全性至关重要。目前，玻璃纤维增强复合材料凭借优异的轻质高强、耐腐蚀性及成本优势，仍是风电叶片主梁最广泛使用的材料。然而，为了满足大型叶片的刚度需求，玻纤主梁的铺层厚度往往达到数十毫米甚至更厚。在大厚度层合板的真空灌注或预浸料铺放成型工艺中，极易受树脂流动冲击、纤维张力不均及固化收缩应力等因素影响，极易产生褶皱缺陷。

2、褶皱是玻纤主梁制造中最为常见的结构缺陷。它破坏了纤维的平直度，导致局部坐标系发生偏转。当主梁承受巨大的轴向压缩载荷时，褶皱区域会产生显著的层间剪切应力和横向拉伸应力。由于玻璃纤维复合材料的层间性能远低于其纤维方向性能，这些离面应力极易诱发基体开裂、局部脱粘，进而导致层间分层与纤维屈曲的耦合失效。这种失效模式往往具有突发性，严重削弱了主梁的极限承载能力和疲劳耐久性。

3、目前，针对含褶皱玻纤主梁的失效分析方法主要面临以下三大挑战：

4、第一，现有的工程折减算法过于保守，导致结构重量过大。工业界通常依据褶皱角度查表选取固定的强度折减系数。然而，这种方法忽略了褶皱波形几何特征（如波长/波高比）对应力集中的非线性影响。为了覆盖工艺的不确定性，设计人员往往被迫采用极大的安全系数。这对于本就厚重的玻纤主梁而言，这种保守设计无疑会导致叶片重量显著增加，进而增大机组的惯性载荷与系统的制造成本，违背了轻量化设计的初衷。

5、第二，单纯的实验测试难以全面覆盖复杂的缺陷形态。玻璃纤维复合材料的性能离散性本身较大，且实际生产中产生的褶皱形态随机多变。若要通过实验穷举所有不同角度、波长及应力比下的失效行为，将耗费难以估量的时间与经费。此外，常规的宏观力学测试只能获得最终的断裂载荷，难以精准捕捉材料内部损伤从基体开裂到分层扩展的动态演化过程，无法揭示失效背后的物理机理。

6、第三，缺乏数据校准的数值模型预测精度不足。现有的有限元模型在模拟褶皱失效时，往往直接采用通用的复合材料本构关系和标准材料手册参数。然而，玻纤主梁的实际层间剪切强度及损伤演化速率受制造工艺影响极大（如树脂富集区的影响）。未经具体实验数据融合修正的"纯理论模型"，难以准确模拟褶皱区域在静载下的非线性剪切行为以及在动载下的刚度退化规律，导致分析结果与真实失效模式存在较大偏差。

7、综上所述，迫切需要一种能够将褶皱诱发层间失效的物理机理与小样本静动态实验特征数据深度融合的分析方法。通过实验数据校准机理模型，消除参数不确定性，从而精准解析玻璃纤维主梁褶皱缺陷的损伤机理、演化规律及失效阈值，这对提升风电叶片的结构安全性与设计经济性具有重要意义。

技术实现思路

1、针对现有技术中，玻璃纤维风电叶片主梁褶皱缺陷的失效分析存在工程折减算法过于保守、纯实验测试覆盖度不足以及传统数值模型缺乏实验数据校准导致预测精度低的问题，本发明提出了一种风电叶片主梁褶皱损伤失效分析方法。

2、为实现上述技术目的，本发明采取的技术方案为：

3、一种风电叶片主梁褶皱损伤失效分析方法，包括以下步骤：

4、步骤一：制备具有可控波高h和波长l的褶皱缺陷的玻璃纤维层合板试样；

5、步骤二：对试样进行静态力学实验，获取静态失效机理基准数据，静态失效机理基准数据至少包括试样的极限承载能力和整体刚度特征；

6、步骤三：对试样进行动态力学实验，获取动态损伤演化特性数据，动态损伤演化特性数据至少包括试样刚度随循环次数的退化曲线；

7、步骤四：建立基于物理机理的参数化有限元模型，模型至少包含表征褶皱几何特征的纤维波纹形态、树脂富集区以及层内与层间损伤起始及演化准则；

8、步骤五：采用三级参数融合校准策略，利用静态失效机理基准数据和动态损伤演化特性数据对参数化有限元模型进行修正；

9、步骤六：将待测风电叶片主梁的实际褶皱几何参数输入修正后的参数化有限元模型，进行失效行为分析并输出预测结果。

10、为优化上述技术方案，采取的具体措施还包括：

11、步骤一中的玻璃纤维层合板试样采用单向玻璃纤维布作为基体铺层，通过在预定铺层深度嵌入异向预浸料作为致褶皱单元，利用真空灌注成型工艺诱导产生平面外弯曲变形，从而形成具有可控波高h和波长l的褶皱缺陷。

12、步骤二中的静态力学实验为单轴静态压缩实验，采集加载全过程的载荷-位移曲线与局部应变数据，将极限承载能力作为校准机理模型的宏观基准。

13、步骤三中的动态力学实验为压缩-压缩疲劳加载实验，在不同应力水平下记录试样弹性模量随循环次数增加而下降的变化曲线，将刚度退化曲线作为校准机理模型中损伤演化参数的物理基准。

14、步骤四中建立基于物理机理的参数化有限元模型的具体实现方式包括：

15、基于abaqus/explicit平台，建立包含褶皱几何特征的代表性体积单元模型，还原波浪状纤维铺层并在波纹空隙处建立树脂富集区；在层间界面插入零厚度内聚力单元；

16、定义正交各向异性弹性本构关系：在vumat用户子程序中，读取材料的工程常数，构建未损伤状态下的柔度矩阵并求逆得到刚度矩阵c；引入损伤变量纤维损伤和基体损伤，构建损伤后的本构方程，其中，损伤刚度矩阵  中的主对角元素及剪切模量根据当前的损伤状态进行折减；

17、植入3d hashin失效准则：在每个增量步中，提取积分点的应力张量，基于物理机理将失效模式划分为纤维拉伸失效、纤维压缩失效、基体拉伸失效和基体压缩失效四种类型，分别建立判据进行并行判定；

18、通过编写vumat用户子程序，建立基于3d hashin准则的正交各向异性损伤本构模型，用于判定层内纤维和基体的损伤起始；采用牵引-分离法则描述层间界面的损伤起始与扩展，采用bk能量准则判定分层扩展，实现层内损伤与层间分层的混合失效模拟。

19、vumat用户子程序中嵌入的3dhashin失效准则将失效模式划分为四种类型，分别建立如下判据：

20、纤维拉伸失效，即时：

21、

22、纤维压缩失效，即：

23、

24、基体拉伸失效，即：

25、

26、基体压缩失效，即：

27、

28、其中，为纤维方向拉/压强度，为横向拉/压强度，为各方向剪切强度，以上强度参数由步骤二中的静态力学实验校准得到。

29、vumat用户子程序中设定状态变量用于记录每个积分点的损伤状态；当任一失效判据满足时，对应的损伤状态变量置为1；基于更新后的损伤变量重新计算损伤刚度矩阵并更新应力状态；引入数值稳定性控制算法，当单元畸变超过阈值时将其删除以模拟物理断裂。

30、步骤五中的三级参数融合校准策略包括：第一级几何与弹性校准：以静态实验测得的局部应变数据为目标，调整模型中的树脂模量和几何特征参数，使仿真局部应变响应逼近实验值；第二级强度校准：以静态实验测得的极限破坏载荷为目标，反向修正模型中的剪切强度参数，使仿真预测的失效载荷与实验值一致；第三级动态演化校准：以动态实验测得的刚度退化曲线为目标，优化损伤演化方程中的退化系数，使仿真复现的刚度衰减趋势与实验值吻合。

31、步骤六中输出的预测结果包括褶皱区域层间分层与层内基体损伤耦合的失效模式云图，通过云图特征区分失效模式是由层间分层主导还是由纤维微屈曲主导，并评估结构安全性。

32、玻璃纤维层合板为单向玻璃纤维增强环氧树脂基复合材料层合板，试样铺层总数、致褶皱单元的预埋位置和厚度根据实际风电叶片主梁的铺层工艺确定。

33、与现有技术相比，本发明具有以下有益效果：

34、1.本发明克服了现有常规有限元模型直接套用理想化标准手册参数所带来的误差，通过将"层间剪切物理机理"与"真实小样本实验数据"相结合，既保留了模型的可解释性，又消除了参数的不确定性，显著提高了对玻纤主梁褶皱失效阈值及刚度退化规律的分析精度。

35、2.相比于仅能提供极限强度的传统静态分析，本发明利用动态刚度退化数据校准模型，能够精准复现褶皱区域在长期载荷下的“变软”过程（刚度衰减），从而更准确地识别失效的前兆特征。

36、3.通过本发明建立的"缺陷特征提取-模型构建-数据校准"闭环分析方法，能够将实际生产中任意实测尺寸的褶皱缺陷高保真地映射到有限元模型中。该方法仅需借助小样本静动态实验即可完成核心力学参数的标定，克服了纯实验方法无法穷举所有复杂缺陷形态的局限性，有效避免了对各类含缺陷结构件进行大量且昂贵的破坏性测试，大幅降低了研发成本。

37、4.通过精准量化褶皱对强度的具体折减程度，替代了传统粗糙的折减因子法，为风电叶片主梁的结构优化设计提供了精确的数据支撑，有助于在保证安全的前提下减轻叶片重量。