轴固有频率计算.docVIP

转子固有频率计算方法比照

本文通过理论计算与ansys模拟两种方法计算转子的固有频率，分别对单盘与多盘情况下作了计算，本文中转子与轴的材料参数如下：

单盘时计算与比照

1、理论计算

中点C处挠度

推出轴的刚度,其中为轴总长度，E为弹性模量，为惯性矩，F为外力

，d为轴的轴径

得：

代入数据有：

质量

2、ansys模态计算固有频率

约束方式：A端铰支，即约束X、Y、Z平动自由度，不约束转动自由度，B端只约束Y、Z自由度

用mass21单元：

combin14

直接约束

Maxtrix27

10.0000

20.0000

30.0000

40.1E-03

584.643

6596.12

71810.9

82798.3

93446.5

105266.5

10.4E-04

286.050

386.050

41526.6

52180.2

62180.2

74131.4

84131.4

97288.2

1011965.

10.0000

20.2E-04

381.439

481.439

51552.1

61552.1

72413.3

82413.3

96868.7

106868.7

3、结论：

1).不加集中质量结果偏差较大

2).直接约束与用combin14和matrix27单元模拟与理论计算结果差不多

二、多盘时计算与比照

模型结构图

考虑多个盘时比照拟复杂，先画出本文结构如下列图：

轴系统固有频率计算

轴系统固有频率计算

ANSYS中模态分析

直接得出固有频率

通过柔度计算刚度，求固有频率

根据轴挠度公式计算得柔度，得固有频率

ANSYS中静力分析求出柔度，推出固有频率

理论推导示意图

1、理论推导

其中：C、D两点为转盘所在位置，AC=CD=DB=，为轴长，A处铰支，B处限制y、z方向自由度。

挠曲轴方程：

将该系统视为两质量弹簧-阻尼系统，通过挠度公式推导柔度矩阵，继而推出刚度矩阵。

柔度计算：〔参考机械振动，张义民，第135页〕

在C点施加单位力，那么

，

其中，E为弹性模量，I为轴惯性矩。

，

在D点施加单位力，那么

所以，柔度矩阵为

根据柔度矩阵与刚度矩阵的关系，有

根据机械振动课本第80页得固有频率计算公式：

〔1〕

其中：

得

2、ansys静力分析计算柔度

约束方式：A端铰支，即约束X、Y、Z平动自由度，不约束转动自由度，B端只约束Y、Z自由度

在ansys中使用BEAM188单元，将轴分6段，分别在节点3和节点5施加沿y轴负方向的单位力，分别读取节点3和节点5沿y方向的位移，即柔度，形成柔度矩阵：

继而，得

用公式〔1〕计算得到固有频率：

结论：柔度的理论计算与ansys中柔度模拟结果根本一致，所得固有频率根本一致。

3、Ansys中模态分析直接计算系统固有频率

在ansys中用mass21单元，运行结果如下：

combin14

直接约束

Maxtrix27

10.0000

20.0000

30.3E-04

40.4E-04

5252.47

61160.9

72539.0

84069.2

95019.3

1013981.

10.0000

2252.84

3252.84

41165.0

51165.0

62189.7

75783.8

813981.

931263.

1036051.

10.0000

20.35E-04

3234.53

4234.53

5869.06

6869.06

74898.3

813163.

913163.

1013644.

4、结论

1).采用施加集中质量方法时，ANSYS模态分析得出的固有频率与通过柔度计算得出的固有频率根本一致；

2).不施加集中质量时，结果差异较大。