涡轮泵涡轮叶片振动特性有限元分析程序实现研究.pdfVIP

涡轮泵涡轮叶片振动特性有限元分析程序实现 高伟田爱梅陈海燕 (北京航空航天大学宇航学院，北京100083) 摘要：引八空间八节点等参数单元，对某涡轮叶片进行了参数化三维造型和有限元网格划 分．利用自行编制的有限元自由振动分析程序对该型叶片进行了模态分析与仿真，给出了前十阶固有 频率以及振型的仿真结果．实现了对满轮叶片进行三维实体造型，网格自动生成和振动特性分析的一 体化流程．最后对计算结果进行了分析，经过对比，本文程序的计算结果与ANSYS软件的计算结果吻 合良好，从而验证了本文的计算程序具有良好的工程计算精度，具有工程实用价值． 主题词：涡轮叶片，有限元，振动，程序 l引言 涡轮泵是液体火箭发动机的关键部件，其涡轮在高转速、高温的恶劣环境下工作．振动是导致叶 片产生裂纹、以致发生断裂的重要原因之一。因此，振动特性是涡轮叶片设计的重要内容。对叶片进 行模态分析，是动力分析的基础。 有限元法被广泛应用于叶片强度振动特性分析中，并取得了良好的效果。在叶片模型建立过程中， 粱单元、壳单元都曾被采用，并取得了一定进展，具有计算量小、程序简单、方法成熟等优点，但其 模型略显粗糙【IJ。三维实体单元则适应性很强，可以较为准确的拟合具有复杂几何外形的叶片的形状， 但箕有限元阿格的生成一直是涡轮叶片实体建模的建点。本文采用空间八节点等参元模型对叶片进行 了参数化三维实体造型与有限元网格自动划分，然后利用自行编制的振动特性有限元分析程序对生成 的叶片有限元模型进行了模态分析，实现了涡轮叶片的参数化三维造型、分网和振动特性分析的一体 化。 2叶片参数化造型及网格自动划分 涡轮叶片是一种包含复杂曲面的空间结构，其复杂的成型规律决定了其实体建模的困难“J。本文 对比工程上常用的解析造型法，绘制弯曲叶型的抛物线中弧线，将厚度叠加到抛物线中弧线之上，运 用四次多项式对叶片流线上的点进行插值拟台逼近使畸型曲线重构，解决了组合曲线曲率倒数不连续 问题，使叶型具有良好的流体动力性能，采用空间八节点六面体等参元模型实现了对涡轮叶片叶型的 三维造型和有限元网格划分，为下一步的振动特性分析做好了准备工作．涡轮叶片三维有限元网格自 动生成的过程如图1所示．对某涡轮叶片实倒进行三维实体造型和网格自动划分。 某液体火箭发动机氢涡轮的基本参数为；叶片数Ⅳ=80，叶高Ⅳ=17ram，几何入口角 届々。66。，几何出口角屈I 成结果如图2所示。采用空间八节点等参元，如图2所示，将涡轮叶片划分为832个单元，共1190个节 21 圉l涡轮时片隔椿生成流程图 其造型及网格生成结果如图3所示。 圈2空间八节点等参元模型 图3涡轮叶片造型及有限元网格划分模型 3振动特’陛分析模型 3．1数学模型 根据动力学基本方程【3】； ， 【M静(f)}+【c水o)】+【足】{a(f))={Q(f)) (1) 不计外载荷和阻尼因素，则上式演变为反映叶片固有振动特性的自由振动方程 、 阻】{fi(哪+【足】{a(f))--{0) (2) 【^，】和[足】分别为系统的质量矩阵和刚度矩阵；{占(r)}和{口(f)}分别是系统的节点加速度向量 和节点位移向量．微分方程(2)的基本解的形式为： {口0)}={钟sin耐 (3) 22 将式(3)代入式(2)得： 旺量卜国2【材D{≯}={o} (4) 记五=m2，则式(2)转化为广义特征值问题的求解： I五渺}--,i1M1{≯) (5) 对叶片旅加约束条件后，求解式(5)即可得到疗个广义特征值暑0=1，2，…