09基于三维有限元复合材料桨叶动力学特性分析.docVIP

第二十八届（2012）全国直升机年会论文 基于三维有限元的复合材料桨叶动力学特性分析 郝启东 李建波 韩 东 （南京航空航天大学直升机旋翼动力学重点实验室，南京，210016） 摘 要：本文基于三维有限元模型，分析符合材料桨叶在不同转速下的固有频率及模态动力学仿真。仿真结果与实验数据相对比表明，误差较小，低阶挥舞和摆振频率都在5%以内，扭转误差较大，大于10%。通过本文研究表明利用三维有限元对直升机旋翼桨叶进行动力学建模是完全可行的，可用于分析计算旋翼桨叶固有特性。 关键词： 复合材料桨叶，动力学特性，三维有限元模型，固有频率 1 引言 分析直升机旋翼动力学特性主要是对其进行固有频率和振型分析。这些特性对直升机的动力学问题起决定性作用，通过调整桨叶的固有频率不与工作转速的激振力发生共振，这对减小直升机的振动水平十分有利[1]。对桨叶进行合理的频率配置和调整，是降低旋翼和机体振动载荷、保证桨叶和机体结构疲劳寿命的重要工程手段。 旋翼桨叶固有频率的计算方法有很多，一般分析计算梁的固有振动特性的方法都能适用。要想对旋翼桨叶进行动力学特性分析，首先需要建立桨叶模型，然后对其进行动力学仿真。而模型的建立有很多方法，由于桨叶的各组成部分的材料不同，也就是其材料的参数不同，因此可以利用材料力学的方法，将不同材料部分当量化成一种材料，即求出这种材料的当量模量、当量泊松比和当量密度，建立模型,然后再对其进行分析计算[3]。 典型桨叶动力分析是把桨叶简化为一维梁，即通过将无穷多个自由度的振动问题简化为有限自由度问题。通过选择有限个数位移函数逼近桨叶弹性振动的振型或把桨叶离散成有限的单元段[5], 具体有雷利法、伽辽金法、三弯矩法、传递矩阵法和动力有限元法等[3]。 对于固有特性的计算，目前流行的做法就是计算桨叶剖面参数，根据计算得到的剖面数据，利用一维梁理论计算桨叶固有频率[6,7]。旋翼桨叶是带有负扭转和尖削的三维弹性结构，因此在有限元软件中进行三维实体建模比较困难。本文主要用有限元软件并借助其他软件对复合材料桨叶的三维有限元模型进行各组成部分分别建模，并对其进行有限元分析。 2 复合材料桨叶三维模型的建立 直升机旋翼桨叶的构成比较复杂，那么准确的建模对以后的动态分析是很有利的。先进直升机的桨叶一般是由复合材料构成，主要有：大梁、泡沫、后缘条、前缘配重以及蒙皮组成，如图1所示。由于这些部件的材料和纤维方向不同，因此桨叶是各向异性结构。 三维有限元软件的分析过程主要有三个主要的步骤：1）创建有限元模型（预处理阶段）；2）施加载荷并求解（求解阶段）；3）查看结果（后处理阶段）。其中创建模型时，我们可以直接输入其它三维软件画好的实体模型，而本文是在有限元软件中进行三维有限元模型的建立。其中旋翼桨叶的主要数据如表1所示。三维模型如图2 所示，剖面网格见图3。 图1 旋翼桨叶剖面表1 旋翼主要参数 翼型 弦长 旋翼直径 工作转速 桨尖速度 单片桨叶长 单片桨叶质量 负扭转 OA212 0.178m 5m 76rad/s 190m/s 2.25m 4.0Kg 0° 图2 桨叶的三维模型图 图3 模型剖面网格本模型的蒙皮是用复合材料进行正负45°铺层而得到的，其模型采用8节点的复合材料体单元，此单元每节点有3个自由度，每个单元最多允许250层的等厚度复合材料，同样允许125层厚度在单元面内成双线性变化的不等厚度材料层。该单元的另一个特点是可以用几个单元叠加的方式对多于250层的复合材料建模并允许沿厚度方向的横向变形斜率可以不连续，而且用户可以输入自定义的本构矩阵。 模型其他部分均采用8节点的实体单元建立模型，此类单元的每个节点有三个自由度，该单元具有塑性、超弹、应力刚化、大变形以及大应变等功能。也可利用混合公式模拟几乎或完全不可压的超弹性材料的变形。剖面建好后，用三角形单元对剖面进行网格划分，然后再将网格拖拽成体，最终得到137060个单元和94066个节点。 3 模型的分析 3.1 复合材料属性的确定 本模型中的C型梁、后缘条以及蒙皮均是复合材料的，因此需要确定他们的材料弹性常数。用有限元软件对复合材料进行分析时要对材料进行属性设置，这就要用到几个主要的参量：E1、E2、ν21和G12。根据复合材料的细观力学分析方法 （1） （2） （3） （4） （5） 其中 ，