ANSYS中的转子动力学计算

时间:2026-01-23

ANSYS中的转子动力学计算

安世亚太 许明财

转子动力学是固体力学的一个重要分支,它主要研究旋转机械的“转子-支承”系统在旋转状态下的振动、平衡和稳定性问题,其主要研究内容有几个方面:临界转速、动力响应、稳定性、动平衡技术和支承设计。在旋转机械研究设计中,转子动力学的性能分析是极其重要的一个方面。

旋转机械广泛应用于以下领域: y 涡轮机械 y 能量站 y 机械工具 y 汽车

y 家用器械 y 航空领域 y 船舶推进系统 y 医疗器械 y 计算机设备

传统的转子动力学分析采用传递矩阵方法进行,由于将大量的结构信息简化为极为简单的集中质量—梁模型,不能确保模型的完整性和分析的准确度;而有限元在处理转子动力学问题时,可以很好地兼顾模型的完整性和计算的效率,但多年来转子的“陀螺效应”一直是

ANSYS 很好地解决了动力特性分析中“陀螺效制约转子动力学有限元分析的“瓶颈”问题。

应”影响的问题,而且陀螺效应的考虑不受计算模型上的限制,使得转子动力学有限元分析变得简单高效。

本文对ANSYS 的转子动力学计算功能进行简要介绍。

1 ANSYS转子动力学的理论基础

ANSYS 转子动力学分析中,两种参考坐标系可供选择:静止坐标系和旋转坐标系。空间点P在静止坐标系(其原点在O′)下的位置矢量为r′,在旋转坐标系(其原点在O)下的位置矢量为r。

&}+([C]+[Cgyr]){u&}+[K]{u}={F} [M]{&u

在静止坐标系下转子的动力方程为:

式中:[Cgyr]为陀螺效应矩阵

在旋转坐标系下转子的动力方程为:

&}+([C]+[Ccor]){u&}+([K] [K[M]{&uspin]){ur}={F} rr

式中:[Ccor]为哥氏效应矩阵, [Kspin]为旋转软化效应刚度矩阵

2 ANSYS转子动力学的计算功能和新技术

ANSYS 转子动力学计算包含如下功能:

y 无阻尼临界转速分析 y 不平衡响应分析 y 阻尼特征值分析 y 涡动和稳定性预测 典型的应用包括: y 轴的弯曲变形 y 扭转振动 y 转子轴未对准 y 旋转部分的平衡 y 流制振动

为了分析时计入哥氏效应、陀螺效应和支承的影响,ANSYS 发展了下列新技术单元:

SHELL181 PLANE182 PLANE183 SOLID185 SOLID186 SOLID187 BEAM188 BEAM189 SOLSH190 COMBIN214

4节点有限应变壳单元 二维4节点结构实体单元 二维8节点结构实体单元 三维8节点结构实体单元 三维20节点结构实体单元 三维10节点四面体结构实体单元

三维一次有限应变梁单元 三维二次有限应变梁单元 三维8节点层合实体-壳单元 二维轴承单元(可变刚度和阻尼)

ANSYS考虑陀螺效应时没有计算模型上的限制,故可选择一维(梁、管)、二维(轴对称)和三维复杂计算模型进行分析。同时,ANSYS 还提供了一系列功能以完善转子动力学的计算,包括:

y CMOMEGA 可以通过组(component CM_NAME)对多个转子指定不同的转速 y CORIOLIS 可以考虑哥氏效应在不同参考坐标系下的影响 y PLCAMP 可绘制坎贝尔图,为临界转速的确定提供了方便 y PRCAMP 可打印固有频率和临界转速 y CAMPB 可绘制预应力结构的坎贝尔图。 y PRORB 可打印出转子涡动幅值

y PLORB 可绘制转子不同截面的涡动轨迹

y SYNCHRO可以指定激励频率是否与结构速度同步 y 多轴转子

y 正向/反向旋转和不稳定性

3 实例

实例1:梁单元与实体单元比较

coriolis,on

omega,2*62.832,0,0

不同单元结果比较

First Bending Second Bending

Torsion Third Bending

185单元

196.42 236.28 658.52 698.06 782.58 1340.9 1380.0

188单元 195.61 235.34 666.36 705.42 782.79 1385.3 1423.5

实例2:圆柱壳模型采用壳单元和实体单元比较

coriolis,on

omega,0,0,314.16,1

不同单元在转速为314.16时的频率值比较

实例3:Campbell图表

验证模型图 分析结果:Campbell图(FW为正转,BW为反转)

ANSYS与文献[1]计算阻尼自然频率比较

Damped Natural Frequencies(Hz)

F(Hz) 1 2 3 4 5 6

BW FW BW FW BW FW

Whirl ANSYS

[1] BW FW BW FW BW FW

ANSYS 271.2 271.2 808.8 808.8 1272.0 1272.0

0rpm

[1] 271.1 271.1 806.4 806.4 1273.0 1273.0

70,000rpm ANSYS 214.5 329.8 762.4 844.9 1068.7 1516.2

[1] 213.6 330.6 760.0 842.6 1066.5 1522.0

ANSYS与文献[1]计算临界速度比较

Critical speeds(rpm)

ANSYS 15.494 17.146 46.729 50.114 64.924 95.747

实例4:多轴转子

y 自振频率

y 不同旋转速度下的临界速度

[1] 15.470 17.159 46.612 49.983 64.752 96.457

[1] The dynamics of rotor-bearing systems using finite elements, Joumal of Eng. For industry-May 1976

模型图 一阶模态

内轴Campbell

图表及不同旋转速度下的临界速度

外轴Compbell图表及不同旋转速度下的临界速度

实例5:鼓风机轴,转速:10,000rpm

模态分析

模型图 频率及相应模态轨迹

频率结果 稳定性结果

Campbell

图表

第一阶FW极限速度

圆盘的不平衡响应

y 圆盘偏心距为.002”

y 圆盘质量为.0276 lbf-s2/in

y 扫频范围为

0-10000 rpm

圆盘位移幅值 极限速度轨迹

轴承反作用力

4 结束语

对于旋转机械,转子动力学 …… 此处隐藏:1058字,全部文档内容请下载后查看。喜欢就下载吧 ……

ANSYS中的转子动力学计算.doc 将本文的Word文档下载到电脑

    精彩图片

    热门精选

    大家正在看

    × 游客快捷下载通道(下载后可以自由复制和排版)

    限时特价:4.9 元/份 原价:20元

    支付方式:

    开通VIP包月会员 特价:19元/月

    注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
    微信:fanwen365 QQ:370150219