叶片模态和刚度的一种计算方法

由天下分享时间：2024/12/19 16:47:58 加入收藏我要投稿点赞

复合材料风电叶片有限元刚度分析

靳交通1，梁鹏程1，曾竟成2，杨跃华1，江翼1

（1. 株洲时代新材料科技股份有限公司，湖南株洲，412007；2. 国防科学技术大学航天与材料工程学院，

湖南长沙，410073）

摘要：运用大型通用软件PROE建立了某兆瓦级风电叶片的三维模型，将所建立的三维模型通过输出输入接口导入到通用有限元软件ANSYS环境下进行有限元分析。分析采用壳单元shell99模拟风电叶片并进行了模态分析和静力分析。风电叶片flapwise和edgewise两个方向上一阶固有频率以及两种载荷工况下的位移分析结果与实测结果基本吻合，从而证明了该方法的可行性。该方法缩短了建模时间，提高了工作效率，对工程上叶片结构校核及新产品开发具有一定的指导意义。关键词：复合材料叶片，有限元，模态，刚度中图分类号：TK83 文献识别码：A

引言

复合材料叶片是风机设备中将风能转化为机械能的关键部件[1]。目前，叶片尺寸正在朝着大型化的方向发展，而其结构性能试验的成本也随之增加，因此，找到一种有效的结构计算分析方法对于节约成本以及结构校核和开发新型叶片就显得尤为重要。

随着计算机技术的发展，有限元法在结构分析中得到了广泛的应用。有限元强大的建模和结构分析功能适用于复合材料叶片的应力、变形、频率、屈曲、疲劳及叶根强度分析。ANSYS是一款著名的商业化大型通用有限元软件，广泛应用于航空航天、机械制造等领域。ANSYS多物理场仿真及耦合的独特功能，以及200多种单元类型，可以对复合材料结构设计、材料研究及制造工艺提供完整的解决方案。总之，对于复合材料结构计算分析，完全可以通过ANSYS程序来实现[2]。

但是，因复合材料叶片结构的特殊性，例如：①形状不规则（每个截面都不同）；②铺层复杂，过渡层很多；③大量夹层结构（“三明治”结构）；④大量粘结区域。所以，其有限元模型的建立是叶片结构有限元分析中的一大难题，而单元类型的选择又决定着建立有限元模型的难易。目前，复合材料风电叶片有限元模型在单元类型的选择上主要采用三种单元类型：shell99壳单元、shell91壳单元、Solid46实体单元。选择实体单元，虽然能提高有限元的计算精度，但是建立叶片的有限元模型会花费大量的工作时间，且很难定义单元坐标，这非常不利于工程上叶片的结构校核及分析；选择壳单元，可以方便地设置和修改铺层厚度，单元坐标的设置容易实现，建模和计算时间比采用实体单元少，这极大地提高了工作效率，而且其计算精度完全可以满足工程需要。

基金项目:国家863项目(编号:2007AA03Z563),湖南省重大科技专项(编号:No.2006GK1002). 第一作者:靳交通(1980.9-),工学硕士,主要从事复合材料风电叶片结构设计工作. Email:jinjiaotong@yahoo.com.cn

因此，本文使用shell99壳单元，通过三维建模，建立了叶片的有限元模型，并以悬臂梁的方式，对叶片的模态和静力变形进行了计算分析，通过计算，得出了叶片的重量、振型及最大变形，并与试验数据进行了对比。 1 有限元模型的建立

通常，在整个有限元求解过程中最重要的环节是有限元前处理模型的建立。一般包括几何建模、定义材料属性和实常数（要根据单元的几何特性来设置，有些单元没有实常数）、定义单元类型，网格划分、添加约束与载荷等。

由于叶片形状复杂，而一般有限元软件所提供的几何建模工具功能相当有限，所以在ANSYS中难以快速方便地对其建模。因此，针对较复杂的结构，可以先在三维CAD软件（如在PROE中）建立几何模型，然后在有限元分析软件ANSYS中通过输入接口读入实体模型，最后，在ANSYS环境下，通过几何修补和简化、板壳中面抽取、节点偏置、网格自动划分等技术对叶片模型进行处理，并形成高效准确的有限元模型，使之适用于CAE分析。 1.1 单元设置与材料属性

针对叶片中的梁、壳等复合材料层合结构，ANSYS提供了一系列的特殊单元——结构多层复合材料单元，以模拟各种复合材料。铺层单元中可以考虑复合材料特有的铺层特性和各向异性特性。

本计算采用的是相对简单的线性铺层单元Shell99。该单元是一种八节点3D壳单元，每个节点有六个自由度，主要适用于薄到中等厚度的板和壳结构，一般要求宽厚比应大于10。Shell99可实现多达250层的等厚材料层，或者125层厚度在单元面内呈现双线性变化的不等厚材料层。如果材料层大于250，用户可通过输入自己的材料矩阵形式来建立模型，还可以通过一个选项将单元节点偏置到结构的表层或底层。

单元铺层主要是确定纤维方向和纤维量，是复合材料风电叶片结构设计的一个重要环节。铺层设计的优劣在很大程度上决定着结构设计的成败[3]。

本计算的铺层完全按照工艺铺层进行设计。在ANSYS环境下，针对Shell99单元，通常有两种方法来定义材料层的配置：①通过定义各层材料的性质；②通过定义表示宏观力、力矩与宏观应变、曲率之间相互关系的本构矩阵。第一种方法是由下到上一层一层定义材料层的配置，底层为第一层，后续的层沿单元坐标系的Z轴正方向自底向上叠加，对于每一层材料，由单元实常数表来定义材料性质、铺层方向角、厚度，如图1所示为叶片某部分的单元铺层；第二种方法是定义各层材料性质的另一种方式，矩阵表示了单元的力-力矩与应变-曲率的关系，必须在ANSYS外进行计算。

图1.单元铺层图

Fig.1 element laminated diagram

叶片的材料体系为玻纤/环氧，叶片制作采用真空灌注工艺，所用复合材料有：三轴向玻璃布、双轴向玻璃布、单轴向玻璃布、PVC泡沫、Balsa木、毡等。玻璃钢复合材料与泡沫材料的主要力学性能见表1、表2所示，其中，Ex为材料的纤维方向，玻璃钢复合材料的密度取ρ=1888kg/m3，Balsa木密度取ρ=150kg/m3，PVC密度取ρ=80kg/m3。

表1.玻璃钢复合材料力学性能

名称玻纤/环氧

Table 1. Mechanical properties of FRP materials

UD Biaxial 符号单位

Ex Ey Ez

Mpa Mpa Mpa

39000 8920 8920

11400 11400 8920

Triaxial 28500 13500 8920

表2.泡沫材料力学性能

Table 2. Mechanical properties of foam materials 名称 Balsa木 PVC

符号 Ex