Ansys转子动力学

更新时间：2024-05-31 15:08:01 阅读量：综合文库文档下载

说明：文章内容仅供预览，部分内容可能不全。下载后的文档，内容与下面显示的完全一致。下载之前请确认下面内容是否您想要的，是否完整无缺。

基于ANSYS的转子动力学分析

1、题目描述

如图1-1所示,利用有限原原理计算转子临界转速以及不平衡响应。

图 1-1 转子示意图及尺寸

2、题目分析

采用商业软件ANSYS进行分析，转子建模时用beam188三维梁单元,该单元基于Timoshenko梁理论，考虑转动惯量与剪切变形的影响。每个节点有6个(三个平动,三个转动)或7各自由度（第七个自由度为翘曲，可选）。

轴承用combine214单元模拟。该单元可以模拟交叉刚度和阻尼。只能模拟拉压刚度，不能模拟弯曲或扭转刚度。该单元如图2-1所示，其有两个节点组成，一个节点在转子上，另一个节点在基础上。

图 2-1 combine214单元

对于质量圆盘，可以用mass21单元模拟，该单元有6个自由度，可以模拟X,Y,Z三个方向的平动质量以及转动惯性。

3、计算与结果分析 3.1 转子有限元模型

建模时，采用钢的参数，密度取7800kg/m3，弹性模量取2.11?1011pa，泊松比取0.3。轴承刚度与阻尼如表1所示，不考虑交叉刚度与阻尼，且为各项同性。

表 3-1 轴承刚度与阻尼参数 Kxx Kyy Cxx Cyy 4e7N/m 4e7N/m 4e5N.s/m 4e5N.s/m 将转子划分为93个节点共92个单元。有限元模型如图3-1所示。

图 3-1 转子有限元模型

施加约束时，由于不考虑纵向振动与扭转振动，故约束每一节点的纵向与扭转自由度，同时约束轴承的基础节点。施加约束后的模型如3-2所示。

图 3-2 施加约束后的有限元模型

3.1 转子临界转速计算

在ANSYS中可以很方便的考虑陀螺力矩的影响。考虑陀螺力矩时，由于陀螺矩阵是反对称矩阵，所以求取特征值时要用特殊的方法。本文考虑陀螺力矩的影响，分析了在陀螺力矩的影响下，转子涡动频率随工作转速的变化趋势，其Campell图如图3-3所示。同时给出了转子的前四阶正进动涡动频率与反进动涡动频率以及固有频率。如表3-2所示。

?(rpm) 0 54.738 54.738 174.12 174.12 301.97 301.97 484.00 484.00 表 3-2 转子涡动频率随转速的变化 10000 20000 30000 54.833 54.131 174.85 173.31 303.56 300.35 488.60 479.44 55.027 53.938 175.61 172.55 305.18 298.76 493.24 474.92 55.248 53.718 176.38 171.78 306.82 297.19 497.93 470.45 40000 55.478 53.489 177.14 171.02 308.46 295.63 502.65 466.02 ?F1(Hz) ?B1(Hz) ?F2(Hz) ?B2(Hz) ?F3(Hz) ?B3(Hz) ?F4(Hz) ?B4(Hz)

图 3-3 转子Campell图

从表3-2与图3-3可以看出,陀螺力矩提高了转子的正向涡动频率,降低了转子的反向涡动频率。同时，陀螺力矩对高阶涡动频率影响更显著。

表3-3对比了不考虑陀螺力矩与考虑陀螺力矩时转子前三阶正临界转速。

表 3-3 临界转速对比不考虑陀螺力矩 ?CF1(Hz) 54.49 ?CB1(Hz) 54.49 ?CF2(Hz) 174.21 ?CB2(Hz) 174.21 ?CF3(Hz) 302.06 ?CB3(Hz) 302.06 考虑陀螺力矩 54.77 54.20 174.89 173.28 304.91 299.09 从表3-3可以看出,陀螺力矩导致了转子系统出现正向临界转速与反向临界转速。同时，陀螺力矩增加了转子正向临界转速，降低了转子反向临界转速。从这个意义上说，陀螺力矩似乎增加了转子正向涡动的刚度。

图3-4至图3-6分别给出了转子在10000rpm时前三阶正向涡动振型图。

图 3-4 转子第一阶正向涡动轨迹

图 3-5 转子第二阶正向涡动轨迹

图 3-6 转子第三阶正向涡动轨迹

3.2 不平衡响应计算

有转子几何尺寸得盘1与盘3质量为7.12Kg，盘2质量为14.24Kg。假设三个圆盘均有质量偏心，其偏心按表3-4分布。所以不平衡量按me?ej?计算。表 3-4 各圆盘偏心分布偏心e(mm) 相位φ(°) 盘1 30 45 盘2 15 60 盘3 20 135 计算2000rpm时转子的不平衡响应。采用ANSYS的瞬态动力学分析模块。图3-7与图3-8分别代表了圆盘1的在2000rpm时的不平衡响应与轴心轨迹图，由图可以看出不平衡响应达2.45mm。

图 3-7 2000rpm时圆盘1不平衡响应

图 3-8 2000rpm时圆盘1轴心轨迹

图3-9与图3-10分别代表了圆盘2的在2000rpm时的不平衡响应与轴心轨迹图，由图可以看出不平衡响应达3.2mm。

图 3-9 2000rpm时圆盘2不平衡响应

图 3-10 2000rpm时圆盘2轴心轨迹

图3-11与图3-12分别代表了圆盘3的在2000rpm时的不平衡响应与轴心轨迹图，由图可以看出不平衡响应达2.4mm。

图 3-11 2000rpm时圆盘3不平衡响应

图 3-12 2000rpm时圆盘3轴心轨迹

同时，计算转子随转速变化的不平衡响应图。采用ANSYS的谐响应分析模块。由于转子各向同性,只给出三个圆盘转子Y向的振动幅值随转速的变化关系。

图3-13到图3-15分别给出了三个圆盘不平衡响应随转速的变化关系。其中盘2的不平衡响应在第二阶固有频率处为0，这是因为转子系统完全对称，盘2在第二阶固有频率处刚好为节点。

图 3-13 圆盘1 Y向不平衡响应随转速变化关系

图 3-14 圆盘2 Y向不平衡响应随转速变化关系