基于ansys经典界面的偏心圆盘启动的转子动力学分析-2.docxVIP

基于ANSYS经典界面的偏心圆盘启动的转子动力学分析－２ ?这是＜基于ANSYS经典界面的偏心圆盘启动的转子动力学分析－１＞的续篇，为方便阅读，重抄题目如下． 【问题描述】 一个转子如下图所示。 该转子两端简支，在距离下端三分之一处有一个轴承支撑；而在距离上端三分之一处有一个刚性圆盘，该圆盘存在一个偏心质量。 现在该转轴从静止开始转动，经过４秒转速达到５０００ＲＰＭ。要求对该转子的启动过程进行仿真，考察偏心圆盘所在处的轴位移和弯曲应力随时间变化的过程。 其它已知条件如下： （１）几何参数 轴的长度: 0.4 m 轴的半径: 0.01 m （２）圆盘的惯性参数: 质量= 16.47 kg 惯性矩 (XX,YY) = 9.47e-2 kg.m2 惯性矩 (ZZ) = 0.1861 kg.m2 偏心质量 (0.1g)位于圆盘上，距离轴心的距离为0.15 m . （３）材料参数： 弹性模量(E) = 2.0e+11 N/m2 泊松比 (υ) = 0.3 密度 = 7800 kg/m3 《注》该算例来自于ANSYS APDL转子动力学部分的帮助实例。 【问题分析】 １.这是一个瞬态动力学问题。 ２.该瞬态动力学问题要考虑科里奥利效应，因此在分析中要打开coriolis开关。 ３.使用ＢＥＡＭ１８８建模转子，使用ＭＡＳＳ２１建模圆盘，使用两个一维的弹簧单元ＣＯＭＢＩＮ１４建模轴承，对这两个弹簧设置不同的刚度，以模拟非轴对称的轴承。 ４.轴上所有节点限制轴向位移以及围绕轴转动的位移。轴的两个端点进一步限制住所有平移，而代表轴承与机架连接的ＣＯＭＢＩＮ１４的机架端节点位移全固定。 ５.使用数组存储不平衡力的两个分量。这里要使用一系列中间变量。具体算法请参见＜基于ANSYS经典界面的偏心圆盘启动的转子动力学分析－１＞。 ６.做瞬态动力学分析，在加载时，对节点施加集中力，并把不平衡力数组设置为其参数，从而施加变化的力。 ７.使用时间历程后处理器，首先提取所需节点的位移和弯曲应力，然后绘制其图形。 ８.本文使用命令流的方式进行讲解。 【求解步骤】 一 建模 １.定义参数 在命令窗口中输入 length = 0.4ro_shaft = 0.01ro_disk = 0.15md = 16.47id = 9.427e-2ip = 0.1861kxx = 2.0e+5kyy = 5.0e+5beta = 2.e-4这里定义了轴的长度等变量，这些变量在下面会用到。至于其具体含义，读者在下面的命令中发现有不知道含义的变量出现时，可以返回来查找，就可以 确定其确切含义。命令运行后，在参数对话框中显示如下２.定义材料模型继续输入/prep7mp,ex,1,2.0e+11mp,nuxy,1,.3mp,dens,1,7800其中，第１行：进入到前处理器第２－４行：定义材料的弹性模量，泊松比和密度。 ３.定义单元类型 （１）定义建模轴的单元 继续输入 et,1,188sect,1,beam,csolidsecdata,ro_shaft,20 其中，第１行：定义单元类型是ＢＡＥＭ１８８第２行：设置该单元类型是实心圆截面第３行：定义该圆截面的半径是ro_shaft。查阅前面的参数，知道 ro_shaft = 0.01 ，即其半径是１０ｍｍ（２）定义建模圆盘的质量单元继续输入et,2,21r,2,md,md,md,id,id,ip其中，第１行：定义单元类型是 ＭＡＳＳ２１第２行：设置该质量单元的惯性量，其中参数md,id,ip 的具体数据请参见前面的参数定义。（３）定义建模轴的两中ＣＯＭＢＩＮ１４单元继续输入et,3,14,,1r,3,kxx,beta*kxx其中，第１ 行：定义单元类型是ＣＯＭＢＩＮ１４.第２行：定义该单元的刚度和阻尼。刚度是kxx；阻尼是beta*kxx 。其大小请参见前面的参数定义。继续输入et,4,14,,2r,4,kyy,beta*kyy其中，第１行：定义单元类型是ＣＯＭＢＩＮ１４.第２行： 定义该单元的刚度和阻尼。刚度是kyy；阻尼是beta*kyy。其大小请参见前面的参数定义。之所以要定义两个ＣＯＭＢＩＮ１４单元，是因为该轴承两个 方向的刚度和阻尼不相同。此时如果打开实常数对话框，可以看到可见，已经输入了３个实常数。 ４.创建轴 继续输入 type,1secn,1mat,1k,1k,2,,,lengthl,1,2lesize,1,,,9lmesh,all 其中，第１－３行：设定下面要用到的单元类型，截面以及材料。第４－５行：创建轴的两个端点第６行：连接这两个端点，得到一条直线第７行：指定将该直线分 为９个单元第８行：划分网格。为了能够看清楚结果继续输入 /PNUM,KP,1?/PNUM,LINE,1/PNUM,NODE,1/PNUM,E