双盘转子动力学仿真c语言程序,[转载]基于ANSYS经典界面的双转子电机的转子动力学仿真...

【问题描述】

一个双转子电机如图所示

该电机含有两个转子：内转子和外转子。

内转子是一根实心轴，较长；它的两端通过轴承与机架相连；在两端距离轴承不远的地方装有两个圆盘(图中没有绘制，在有限元分析中圆盘会用质量单元表示)，而且右边的圆盘上存在不平衡质量，该不平衡质量产生了不平衡的力。

外转子是一根空心轴，它套在内转子外面。外转子的左端与机架通过轴承相连，右端面通过轴承与内转子连接(图中没有表示出来)。在外转子上也有两个圆盘，这两个圆盘不存在偏心质量的问题。

内转子的转速是14000转每分，而外转子的转速是21000转每分。

所有的相关几何尺寸，轴承的参数，以及圆盘的质量和惯性量，在下面建模的时候给出。

现在要对该双转子电机进行转子动力学仿真，具体是做谐响应分析，目的是考察：

(1)7号节点(内转子上)和12号节点(外转子上)的幅值与频率的关系图。也就是要绘制这两个点的幅频关系曲线。

(2)在某一个给定频率处的转轴轨迹图。

(3)在某一个给定频率处转轴的涡动动画。

《注》该算例来自于ANSYS APDL转子动力学部分的帮助实例。

【问题分析】

1.

对于内转子用梁单元BEAM188建模,对于外转子也用BEAM188建模。由于这里涉及到圆盘的位置，集中质量的位置，准备用直接建模法。这就是说，先创建节点，然后由节点创建单元。

2.

对于4个轴承，使用COMBI214建模，该单元是二维的弹簧/阻尼单元，支持在两个方向上定义刚度和阻尼特性。

3.对于4个圆盘，使用MASS21建模。质点单元创建在相应的转轴上，设置其质量和转动惯量。

4.由于内外转子的转速不同，需要分别定义两个组件，并对每个组件给以不同的转速。这两个组件只包含轴以及轴上的集中质量。

5.位移边界条件的处理。

轴承是用COMBI214建模的，该单元有两个节点，其中一个节点在转子轴上，另外一个在转子轴外，对于转子轴外的这个节点完全固定，以此来模拟机架。

此外，所有的节点都限制围绕轴线方向的转动。以限制扭转自由度。

5.对于不平衡力的处理。基于对内转子右边圆盘的偏心质量的计算，得到不平衡力，把该不平衡力使用力的滑移定理移动到轴上。因此直接在轴上施加集中力就可以表示该不平衡力。

6.分析设置。使用谐响应分析方法，开启科里奥利效应，并对两个转子组件施加不同的转速，使用固定坐标系进行计算。

7.后处理。使用时间历程后处理得到两个点的频率响应曲线，使用通用后处理得到转轴轨迹图和涡动动画。

8.关于GUI和命令方式的选择。由于本算例过程相对复杂，如果用GUI方式，篇幅很长。本文准备对重复的部分用命令方式，对简单的部分用GUI方式。对于命令方式，笔者会阐述该命令的含义，读者只需要把相关的命令拷贝到ANSYS的命令窗口中执行就好。

【求解过程】

一 建模

本部分的操作均用命令流的方式进行。

进入到前处理器

1. 定义材料属性

在命令窗口输入下列命令

mp,EX ,1,2.1e+11

mp,DENS,1,7800

mp,PRXY,1,0.3

上面三个命令分别定义了弹性模量，密度和泊松比。

即指定只有一种材料，该材料是钢材。

2. 创建梁单元并定义截面

在命令窗口输入下列命令

et,1,188,,,2

sectype,1,beam,csolid secdata,0.01524,32

sectype,2,beam,ctube secdata,0.0254,0.03048,32

第一行定义了单元BEAM188；

第二行定义了一种实心圆截面梁，代表内转子；

第三行定义了一种空心管截面梁，代表外转子。

3. 创建质量单元并定义其实常数

在命令窗口输入下列命令

et,2,21

r,3,10.51,10.51,10.51,8.59e-2,4.295e-2,4.295e-2

r,4,7.01 ,7.01 ,7.01 ,4.29e-2,2.145e-2,2.145e-2

r,5,3.5 ,3.5 ,3.5 ,2.71e-2,1.355e-2,1.355e-2

r,6,7.01 ,7.01 ,7.01 ,6.78e-2,3.390e-2,3.390e-2

其中第一行定义了MASS21单元，用它来模拟圆盘。

接连下来的四行定义了四个圆盘的质量和转动惯量。

4. 创建轴承单元并定义实常数

在命令窗口输入下列命令

et,3,214,,1

r,7 ,2.63e+7 ,2.63e+7

r,8 ,1.75e+7 ,1.75e+7

r,9 ,0.875e+7,0.875e+7

r,10,1.75e+7 ,1.75e