高速旋转轮盘模态分析.docx
《高速旋转轮盘模态分析.docx》由会员分享,可在线阅读,更多相关《高速旋转轮盘模态分析.docx(18页珍藏版)》请在冰点文库上搜索。
高速旋转轮盘模态分析
密度 DENS=7.8E-9Tn/mm 。
高速旋转轮盘模态分析
在进行高速旋转机械的转子系统动力设计时,需要对转动部件进行模态分析,求解出其
固有频率和相应的模态振型。
通过合理的设计使其工作转速尽量远离转子系统的固有频率。
而对于高速部件,工作时由于受到离心力的影响,其固有频率跟静止时相比会有一定的变化。
为此,在进行模态分析时需要考虑离心力的影响。
通过该实验掌握如何用ANSYS进行有预应
力的结构的模态分析。
一.问题描述
本实验是对某高速旋转轮盘进行考虑离心载荷引起的预应力的模态分析,求解出该轮盘
的前 5 阶固有频率及其对应的模态振型。
轮盘截面形状如图所示,该轮盘安装在某转轴上以
12000 转/分的速度高速旋转。
相关参数为:
弹性模量 EX=2.1E5Mpa,泊松比 PRXY=0.3,
3
1-5 关键点坐标:
1(-10, 150, 0)
2(-10, 140, 0)
3(-3, 140, 0)
4(-4, 55, 0)
5(-15, 40, 0)
L=10+(学号×0.1)
RS=5
二.分析具体步骤
1. 定义工作名、工作标题、过滤参数
①定义工作名:
Utility menu > File > Jobname
②工作标题:
Utility menu > File > Change Title(个人学号)
2. 选择单元类型
本实验将选用六面体结构实体单元来分析,但在建模过程中需要使用四边形平面单元,
所有需要定义两种单元类型:
PLANE42 和 SOLID45,具体操作如下:
1
Main Menu >Preprocessor > Element Type > Add/Edit/Delete
①“ Structural Solid”→“ Quad 4node 42” →Apply( 添加 PLANE42 为 1 号单元)
②“ Structural Solid”→“brick 8node 45” →ok( 添加六面体单元 SOLID45 为 2 号单元)
在 Element Types (单元类型定义 )对话框的列表框中将会列出刚定义的两种单元类型:
PLANE42、 SOLID45,关闭 Element Types (单元类型定义)对话框,完成单元类型的定义。
3. 设置材料属性
由于要进行的是考虑离心力引起的预应力作用下的轮盘的模态分析,材料的弹性模量 EX
和密度 DENS 必须定义。
① 定义材料的弹性模量EX
Main Menu >Preprocessor > Material Props > Material Models> Structural > Linear >
Elastic >Isotropic
弹性模量 EX=2.1E5
泊松比 PRXY=0.3
② 定义材料的密度DENS
Main Menu >Preprocessor > Material Props > Material Models>density
DENS =7.8E-9
4. 实体建模
对于本实例的有限元模型,首先需要建立轮盘的截面几何模型,然后对其进行网格划分,
最后通过截面的有限元网格扫描出整个轮盘的有限元模型。
具体的操作过程如下。
① 创建关键点操作:
Main Menu > Preprocessor > Modeling > Create > Keypoints > In Active CS
列出各点坐标值 Utility menu >List > Keypoints >Coordinate only
2
1
POINTS
TYPE NUM
MAR 28 2013
1
2
12
3 10 11
20:
33:
54
4
9
5
6
Y
Z X
8
7
1006024223chenli
② 由关键点生成线的操作:
Main Menu > Preprocessor > Modeling > Create > Lines > Lines > In
Active Coord
3
1
LINES
TYPE NUM
APR 2 2013
12:
58:
41
112
23 10 11
4
9
5
6
Y
Z X
8
7
1006024223chenli
③ 建立圆角:
Preprocessor > Modeling > Create > Lines > Lines > Lines Fillet
4
1
LINES
TYPE NUM
APR 2 2013
13:
01:
02
112
23 10 11
4
9
5
13
14
Y
Z
8
15
16
X
1006024223chenli
④ 生成面:
Main Menu > Preprocessor > Modeling > Create >Areas >Arbitrary >By Lines(逆时
针选线)
5
1
AREAS
1 12
TYPE NUM
MAR 28 2013
23 10 11
20:
41:
44
49
58
1316
1415
Y
ZX
1006024223chenli
5.划分网络
Main Menu >Preprocessor >Meshing >MeshTool
①对全局进行设置。
单击 Size Controls (尺寸控制区)全局设置项(Global)的 Set 按钮,将弹出
Global Element Sizes 单元尺寸全局设置对话框在对话框中输入 Element edge Length (单元边长
度)为 6。
②单击 Size Controls (尺寸控制区) Lines (线设置项)的按钮,将弹出 Element Size on Picked
Lines (在所选线上定义单元尺寸)的拾取对话框。
用鼠标左键在图形输出窗口中拾取圆角对应
的线。
单击 ok 按钮,将弹出 Element Size on Picked Lines (在所选线上定义单元尺寸)对话框,
在对话框中输入 No. of element divisions ( 每条线将要分成的单元数)为 1,设定圆角处对应的
线只分一个单元。
(由于是模态分析,只要能反应出需要知道的前几阶模态就行,而不需要知
道具体的应力值,所以不需要对此处进行单元细化。
)
③对分网进行控制。
在分网控制区的 Mesh 下拉框中选定分网类型为 Area (面),Shape (网格
形状)设置为 Quad (四边形),分网方式设置为 Free (自由分网)。
④对面进行分网。
在 MeshTool (分网工具)对话框中单击 Mesh 按钮,将弹出 Mesh Aeras (对面
6
划分网格)拾取对话框。
从图形输出窗口中拾取创建的面,单击 ok 按钮。
完成网格划分。
1
ELEMENTS
APR2 2013
13:
25:
50
Y
ZX
1006024223chenli
6. 出整个轮盘的有限元模型
通过将面绕轴旋转成有限元实体模型的功能,将前面建立的轮盘截面有限元网格,围绕
定义的旋转轴扫掠成整个轮盘的实体有限元模型。
具体的操作过程如下。
① 定义旋转轴。
可以通过定义旋转轴所在轴线上的两个关键点来,指定旋转轴的位置。
Main Menu >Preprocessor >Modeling >Create >Key points >In Active CS
生成两个关键点 20、21。
关键点 20:
X,Y,Z 位置分别为-10,0,0
关键点 21:
X,Y,Z 位置分别为 10,0,0
7
② 设置单元生成选项
Main Menu >Preprocessor >Modeling >Operate >Extrude >Elem Ext Opts,
弹出 Element Extrusion Options (单元挤出选项)对话框,在对话框中的 Element type number (单
元类型序号)下拉框中选择 2 号单元 SOLID45。
单元尺寸选项中的分割单元数(VAL1 NO.
Elem Divs)设置为 18,即在挤压出的每个体上将沿周向被分成 18 份。
拉伸比例为 0,保持
等截面拉伸。
将 Clear area(s) after ext (删除原始面)设置为 Yes,在挤压的单元完成之后将删除
原来的面以及其上的单元。
单击 ok 按钮,完成对单元选项的设置。
③ 绕轴旋转截面
Main Menu >Preprocessor >Modeling >Operate >Extrude >Areas >About Axis,将弹出 Sweep
Areas about Axis ( 绕轴扫描面)的拾取对话框。
从图形输出窗口中选择创建好的平面网格,单
击拾取对话框中的按钮。
然后从图形窗口中选取定义旋转轴的关键点 20,21,单击 ok 按钮,
将弹出 Sweep Areas about Axis 对话框,在对话框中输入旋转角度为(Arc length in degrees)360,
No. of volume segments (一周创建体的数目)为 4,单击 ok 按钮。
创建如图所示的整个盘的
有限元模型。
8
1
VOLUMES
TYPE NUMMAR 28 2013
20:
51:
16
Y
Z
X
1006024223chenli
④ 观察创建的网格形式。
Utility Menu >Plot >Element,图形窗口中将会显示出由平面网格扫掠而成的实体单元网格情
况。
存盘,SAVE_DB。
至此,完成了创建轮盘有限元模型的所有工作。
9
1
ELEMENTS
MAR 28 2013
20:
52:
24
Y
Z
X
1006024223chenli
7. 节点的坐标变换
根据轮盘的工作情况其约束条件为盘心轴向和周向约束,这种约束条件在直角坐标系下
无法定义,而柱坐标下可以非常方便地定义。
根据 ANSYS 程序中坐标系的定义规则,需要
将柱坐标系的 Z 轴和旋转轴重合,Y 轴表示转角,X 轴表示径向。
ANSYS 程序提供的全局柱
坐标系不满足要求。
通常可以有两种办法来解决这个问题:
a.将所建有限元模型进行旋转使其轴向和柱坐标 Z 轴方向一致。
b.重新建立一个柱坐标系使其的 Z 向和旋转轴一致。
本实例采用第二种方法。
具体操作过程如下:
① Utility Menu >WorkPlane >Offset WP by Increments ,弹出 Offset WP ( 工作平面偏移)菜单,
拖动 Degrees 滑动条,将 Degrees (旋转角度)值设置为 90。
单击按钮,使工作平面绕 Y
轴正向旋转 90 度,单击 ok 按钮,将工作平面的 WZ 轴和总体坐标系的 X 轴方向重合,。
② 在工作平面原点创建柱坐标系。
Utility Menu >WorkPlane >Local Coordinate Systems >Create Local CS >At WP Origin ,将弹出
Create Local CS at WP Origin (在工作平面原点创建本地坐标系)对话框,Ref number of new
10
coord sys (新坐标系的参考序号)缺省值为 11,一般就使用缺省值,也可自己设定。
在 Type of
coordinate system(坐标系类型)下拉框中选取 Cylindrical 1 (柱坐标系),其它设置为缺省值。
单
击 ok 按钮。
将完成要求的柱坐标系的创建,并且将新建坐标系定义为当前激活坐标系。
③ 将所有节点移到当前柱坐标系中。
Main Menu >Preprocessor >Modeling >Create >Nodes >Rotate Node CS >To Active CS将弹出
Rotate Nodes into CS 菜单。
在菜单中单击
按钮,将所有的节点都移到当前激活柱坐
标系下。
8. 进行静力分析
由于对轮盘模态的分析需要考虑离心力引起的应力对模态的影响,所以需要先对其进行
静力分析,求解出离心力产生的应力,及其对刚度阵的影响,将结果写入数据库文件。
有预
应力的轮盘静力分析具体过程如下。
① 指定分析类型及分析选项
a.Main Menu >Solution >New Analysis→“Static”
b.Main Menu >Solution > Sol’n Controls ,将弹出求解控制(Solution Controls) 对话框。
单击标
签“Basic”,在 Calculate Prestress effects 选项前打“√”。
打开预应力选项。
单击 OK
11
② 定义边界条件和转速
对于本实例分析的轮盘,由安装条件知道其边界条件应该是,在轮盘盘心的节点轴向和
周向固定,而径向自由。
其离心载荷是由于高速旋转产生的,因此需要在分析时指定轮盘的
旋转速度,具体操作过程如下:
a.将要加载模型放大。
在图形窗口中选取盘心部分,对其进行放大,以便在加载时能够准确
地选择所要加载的模型元素(本实例为节点)。
b.定义约束 Main Menu >Solution >Define Loads >Apply >Structural >Displacement >On Nodes,
将弹出 Apply U,ROT on Nodes ( 给节点施加约束)对话框,单击对话框中的 Circle 单选按钮,
选用圆形区域选择方式(因为需要约束的为盘心节点)。
将光标移至盘心,按住鼠标左键并拖
动光标,在图形窗口中将会出现一个选择用的圆形选择框随光标移动,当圆形框刚好将盘心
所有的第一排节点选中,而没有选择第二排节点时,放开按钮。
将会将所有盘心的节点选中,
后面的约束也将施加在这些节点上。
单击 ok 按钮,将弹出 Apply U,ROT on Nodes (给节点施
加约束)的对话框,在对话框的 DOFs to be constrained (被约束自由度)列表框中选择 UY,UZ 两
个自由度,分别表示对周向和轴向施加约束。
在 Displacement value ( 位移约束值)文本框中输
入 0(缺省值为:
0)。
12
1
ELEMENTS
UMAR 28 2013
21:
05:
53
Y
XZ
1006024223chenli
c.定义转速
Main Menu >Solution >Define Loads >Apply >Structural >Inertia >Angular Velocity >Global ,将
弹出 Apply angular velocity (角速度定义)对话框,在对话框中的 Global Cartesian X-comp (绕总
体坐标系中 X 轴旋转角速度)文本框中输入 1256.64,指定轮盘的旋转速度为 1256.64 弧度/秒。
单击 ok 按钮,完成对旋转角速度的定义。
存盘 SAVE_DB。
③ 进行求解
Main Menu | Solution | Current LS
9.进行模态分析设定
由于模态分析时位移约束条件和进行静力分析时一样,而静力分析时所定义的转速在模
态分析时将被忽略,因此进行模态分析时再不需要重新定义边条,可以直接进入求解菜单进
行求解。
① 指定分析选项
Main Menu >Solution >Analysis Type >New Analysis→“Modal”(模态分析)
13
② 分析选型设定
Main Menu >Solution >Analysis Type >Analysis Options
将弹出模态分析选项对话框,在对话框中,指定 Mode extraction method (模态提取方法)为 Block
Lanczos (分块兰索斯法),并指定 No. of modes extract (提取模态的阶数)为 10 ,将 Expand mode
shapes (模态扩展)设置为“YES”,在 No. of modes to expand (模态扩展阶数)文本框中输入 10。
这
样 ANSYS 程序在进行模态求解的同时将完成模态的扩展,不需要再单独进行模态扩展。
将
Incl prestress effects (预应力效应)设置为“YES”,这样在进行模态分析时 ANSYS 程序将会把
前面静力分析中求解得到的离心力产生的应力对刚度的影响考虑进去。
(重要) 单击 ok 按钮,
将会弹出 Block Lanczos Method (兰索斯法模态分析选项)对话框,在对话框中对模态提取的范
围进行定义,一般定义的范围要求将所关心的所有频率都包含在指定的范围之内,这里指定
的值为:
Start Freq (开始频率)是 0,End Frequency (结束频率)是 99999999。
单击 ok 按钮,完
成对分析选项的设置。
③ 进行求解
Main Menu | Solution | Current LS
10. 后处理
① 列出固有频率
Main Menu >General Postproc >Results Summary,将列出轮盘的所有求解的固有频率,在文本
框里列出了轮盘的前 10 阶固有频率,可以看出有些频率值相同,这是由于轮盘结构和边条都
是对称的,会出现振型和频率相同但相位不同的情况。
14
② 观察解得的模态
本实例中由于设置了对模态进行扩展,所以对于求得的每一阶固有频率,程序同时都求
解了其对应的模态振型,来反映在该固有频率时,轮盘的各节点的位移情况。
可以利用通用
后处理器方便地对其进行观察和分析,并且可以对各阶模态振型进行动画显示。
具体操作过
程如下:
a.选取菜单路径 Main Menu >General Postproc >Read Results >First Set,选择轮盘第一阶模态。
b.选取菜单路径 Main Menu >General Postproc >Plot Results >Nodal Solu
选择 DOF→
从结果可以看出,轮盘第一阶模态的固有频率为 544.004Hz,对应的振型为一节径。
c.观察第一阶模态的动画显示。
Utility Menu >PlotCtrls >Animate >Deformed Results
选择 DOF
图形显示窗口中将会显示本节模态的动画。
同时在工作目录下会将生成的动画文件保存下来,
以后可以观察。
观察完动画显示之后,单击 close 按钮,关闭动画控制器。
d. 选取菜单路径 Main Menu >General Postproc >Read Results >Next Set,选择轮盘第二阶模态。
15
e.选取菜单路径 Main Menu >General Postproc >Plot Results >Nodal Solu
选择 stress→von Mises→OK
图形窗口中将显示出第二阶模态振型。
f.选取菜单路径 Utility Menu >PlotCtrls >Animate >Deformed Results
选择 stress→von Mises→OK
图形显示窗口中将会显示二节模态的动画。
g.对余下的各阶模态重复步骤 d~f,可观察到所求解的各阶模态的振型和动画显示。
下面为本实例求解得到的其余不同频率时的振型图:
1
NODAL SOLUTION
STEP=1
SUB =1
FREQ=546.288
USUM(AVG)
RSYS=0
DMX =31.334
SMN =.003841
SMX =31.334
MAR 28 2013
21:
15:
56
MX
Y
Z
X
MN
.003841
6.966 13.928 20.891 27.853
3.48510.44717.40924.37231.334
1006024223chenli
16
1
NODAL SOLUTION
STEP=1
SUB =3
FREQ=616.279
MAR 28 2013
21:
24:
16
UX(AVG)
RSYS=0
DMX =21.948
SMN =-.497E-16
SMX =21.845
MX
Y
ZMN X
-.497E-16
4.854
9.709
14.563
19.418
2.427
7.282
12.136
16.991
21.845
1006024223 chenli
17
1
NODAL SOLUTION
STEP=1
SUB =5
FREQ=744.069
UX(AVG)
RSYS=0
DMX =32.818
SMN =-32.527
SMX =32.527
MAR 28 2013
21:
25:
03
MX
Y
X
Z
-32.527
MN
-18.071 -3.614 10.842 25.299
-25.299-10.8423.61418.07132.527
1006024223chenli
18