基于非线性本构关系的复合材料风机叶片有限元极限分析与设计.docVIP

基于非线性本构关系的复合材料风机叶片有限元极限分析与设计 同济大学 张春丽? 黄争鸣? 董国华?? 为了研究具有三维复杂构形的复合材料风机叶片的逐次破坏过程和极限承载能力，将复合材料细观力学非线性本构理论桥联模型与有限元软件ABAQUS通过用户子程序UGENS结合起来对风力发电机叶片结构进行极限强度分析。只需提供纤维和基体的材料性能参数、纤维体积含量以及蒙皮和增强筋的铺层数据包括铺设角、层厚和铺层数，就可预报出复合材料复杂叶片结构的整体承载能力以及叶片破坏所处的位置，为正确评估和合理设计风机叶片结构提供了一种简便有效的分析方法。本文以一种20kW风机叶片为例，用此方法实现了新型复后材料叶片结构的极限分析 风能是一种清洁的可再生能源，取之不尽、用之不竭。叶片是风力发电机中最主要的部件。目前的大、中型风机叶片基本上采用蒙皮与主梁的构造形式，须通过多步成型工艺制备，即先分别制作叶片的上、下外壳和龙骨梁（腹板）后，再粘成一体。由于粘接处的强度远低于壳体本身的强度，使叶壳性能得不到充分发挥，类似开口薄壁梁远不及闭口薄壁梁的承载能力。单腹板支撑的叶壳易发生失稳破坏，多个梁或腹板则须增添更多模具。这无疑会增加成本，降低叶片的利用率。因此，本文作者提出采用整体一次成型技术制备中空叶片，可有效减轻重量，降低成本，提高叶片的整体力学性能。这就需要对这种新型叶片结构进行极限分析并在此基础上实现合理设计，因为传统叶片以龙骨梁（腹板）为主承力件。 ??? 在目前的设计中，叶片的结构参数（蒙皮厚度，腹板厚度、宽度、位置等）一般是通过有限元法(FEM)分析后决定，直到设计的叶片满足规范。Oh等人运用复合材料梁理论结合有限元法预测了叶片的静态响应和无阻尼动态响应。Saravanos进一步采用梁单元预测了复合材料叶片的有阻尼动态响应。Maheri运用壳单元对叶片进行结构划分，并通过考察不同区域的单元密度来考查收敛性并避免应力集中现象。Kong等运用有限元法对叶片进行静强度分析，将叶片视作蒙皮-腹板-泡沫芯结构。用壳单元对蒙皮进行离散，用12节点三维“夹心”单元模拟腹板，通过对蒙皮、腹板的厚度设计使整个结构满足叶片设计规范对刚度和强度的要求。这些分析都是借助有限元软件基于材料线弹性本构关系进行的，并没考虑复合材料的非线性，使叶片极限承载能力的计算与实际情况存在差异。实际上，复合材料层合板在逐层破坏过程中基体材料表现出非线性特性,导致层合板的刚度矩阵呈非线性变化。桥联模型充分考虑了基体材料的非线性特性对复合材料本构方程的影响,将其应用于有限元分析能较好解决上述问题。 1 桥联模型 ??? 经典层板理论中，复合材料结构第k层的平面应力增量与平面应变增量之间的关系是 ??? 其中：G表示总体坐标，是该复合材料层在总体坐标系下的当前刚度系数。???? ??? 以往分析都是采用初始线弹性刚度矩阵并在整个加载过程中保持不变，因此会产生计算误差。而复合材料非线性本构理论—桥联模型为该问题的解决提供了有效途径。? ??? 在单向复合材料受到外力增量作用时，纤维和基体内产生的应力增量和间存在一个非奇异矩阵相联系，即 ??? 其中： [A]称为桥联矩阵，桥联模型也由此而得。 ??? 应用关系式(2)，可导出单向复合材料的当前柔度矩阵为 ??? 其中：和分别是纤维和基体的体积含量，和分别是纤维和基体的当前柔度矩阵， [I]是单位矩阵。纤维可以看作直到破坏都是线弹性的。从而，复合材料的当前柔度矩阵依赖于基体的当前柔度矩阵，后者取决于基体的当前应力。根据式(3) ,可得?????????????????? ??? 其中：[T]c 是坐标变换矩阵，L表示局部坐标，上标T代表转置。 ??? 假设为施加到单向复合材料在局部坐标系下的外应力增量，则可得基体以及纤维中的内应力增量为： ??? 基体、纤维和单层板中的总应力按下式更新： ??? 其中，初始。基于当前的可以确定基体的当前材料参数，用于更新基体的柔度矩阵。由式(6)得到的可代入经典的第一强度理论，根据纤维或基体是否破坏来判断单层板是否达到破坏。更详细内容参见文献。 2 桥联模型的程序化 ??? 桥联模型本构理论须通过用户子程序与ABAQUS实现接口，才可用于风机叶片的极限分析。本文中使用平面壳单元模拟叶片结构,在有限元分析中运用增量理论,对壳单元的各个积分点进行计算。假定单位长度上的内力和内力矩增量分别为dNxx 、dNyy 、dNxy 、dMxx 、dMyy 、dMxy，它们须与截面上的应力合力平衡。据此得到： ??? 其中：是层合板的厚度，和分别是第k层的上顶面和下底面的z坐标，这些都是已知数据，在有限元模拟