r40轴流风机动叶片的动力学分析.docx

r40轴流风机动叶片的动力学分析 1 叶片的动态特性 作为一种旋转机械，轴流轴的运动特性不同于其他传统机械。叶片式流体机械具有复杂的结构和力学特性, 具有大量的过流曲面零部件, 特别是航空翼型应用于叶片之中, 叶片的剖面较为复杂, 同时整个叶片存在一定的扭角, 使得叶片的空间结构更加复杂。在流体机械的运动过程中, 叶片结构受力比较复杂, 有叶轮旋转过程中产生的离心力, 同时也存在着“流-固”耦合作用, 例如气流对叶片的冲击、叶片脱流激振等。 叶片在工作时的受力较为复杂, 一方面由于叶轮在旋转过程中会产生离心力;另一方面, 气流在通过叶片时会对叶片产生气流力的作用。在叶片受到的几种力当中, 最主要的是气流激振力, 非稳态的气流力是气动激振力的主要来源。同时, 气流在经过叶片时会与气流之间产生“流-固”耦合作用, 叶片在气流的作用下会发生弯曲和扭转变形, 这些变形将会导致叶片气动性能的改变。叶片的振动将产生较大的动态应力, 导致叶片的疲劳, 缩短叶片的使用寿命, 严重时甚至发生断裂。因此, 对风机叶片结构谐响应分析具有十分重要的意义。 2 有限模型与网格划分 2.1 叶片叶型的表征 本文研究对象为R40风机叶片, 叶片选用CLARK-Y翼型, 相对厚度为12%, 工作时的转速为1440rpm。叶片各控制截面的参数如表1所示。叶片材料为玻璃纤维增强塑料, 弹性模量为MPa, 密度2600kg/m3, 泊松为0.3。 对于风机的叶型是通过多个截面的参数来定义的, 通过曲面连接得到扭曲的风机叶型形状。所以在建模过程中, 采用了由点到线再到面, 最后到体的建模方法。 2.2 u3000叶片有限元方程 本文选用了Solid186单元对结构进行了离散, Solid186是一个高阶3维20节点固体结构单元, 具有二次位移模式可以更好的模拟不规则的网格。单元通过20个节点来定义, 每个节点有3个沿着xyz方向平移的自由度。Solid186可以具有任意的空间各向异性, 单元支持塑性、超弹性、蠕变、应力钢化以及大变形和大应变能力。 叶片在离心力场下被离散后的平衡方程为 式中, [B]为单元的应变矩阵, {σ}为单元内点的应力向量, [N]为单元形函数矩阵, {f}为单元单位体积内的惯性力向量, [P]为单元单位面积上的压力, {R}e为单元的节点集中力载荷向量。叶片的有限元方程为: 从图中可以看出, 叶片在气流激振力下的响应位移最大值为0.919mm, 最大应力值为14.8MPa, 远低于玻璃钢材料强度许用应力(255MPa) , 叶片在工作转速下满足设计要求。 图3给出了叶片上A点随激振力频率变化的响应位移大小, 从图中可以看出, 当激振力的频率从110Hz增加到120Hz时, A处的径向响应位移急剧增加, 激振力的频率再增加到140Hz时, A处的径向响应位移急剧下降, 可以进一步确定, 共振发生在120Hz附进, 此时A处的径向响应位移为5.72mm, 是激振力的频率为96Hz时对应响应位移 (0.155mm) 的36.9倍。图4为叶片上B点随激振力频率变化的响应应力大小, 最大响应应力为25.7226Mpa, 是频率为96Hz时对应响应应力 (0.90423MPa) 的27.65倍。 4 叶片固有频率及气流激振力分析 本文应用Ansys软件对风机叶片进行了动力学分析, 得到了叶片在静止和工作转速下的前10阶固有频率, 以及在气流激振力下的谐波响应, 可得出如下结论: 式中, [M][C][K][MC]分别为叶片的总质量矩阵、阻尼矩阵、刚度矩阵以及离心力质量矩阵。 3 气流激振力作用下叶片的响应 应用计算流体力学Fluent软件对风机流场进行了数值模拟, 得到作用在叶片上的气流激振力大小为140.9N, 将该力加载到叶片上, 计算了叶片在工作转速下的响应位移和响应应力, 结果如图1和图2所示。 4.1 r40轴流轴通过指定的旋转次数工作，运行稳定，不受共振。叶片最大高差小于材料的允许重量，满足设计要求 4.2 抗振剂最大限度地响应在120hz附近。在这种情况下，a点的径向响应位移和b点的响应能力分别为96hz，对应值的36.9倍和27.65倍