基于传递矩阵法和有限元法的转子动力学分析.pdfVIP

业信息化 所 以系统的传递矩阵可写为 7800kg／m ，端面直径 0．Olm，分段数 20段 ，弹性模量 {z}fl={ { I．．·{t{刁oR (4) 2．1xl0MPa，点击计算后 ，计算完成用户界面就会输 出各 阶固有频率 ，并且也会输 出各阶振型图，如 图2所示。 4 Y 或 0 4 0 M M g Q 葶．擅 ， m 啪， ～ m 厂 一 l 以Q0H=峨 =QnH=0为边界条件得到满足此边界条件的 t、 ’ol1m 频率方程为 一甍黔撩越时 骱霹吐后吝阶巅 §膏序．蓝，嚣，虹。鬻 羹1， ／： ＼ ， 骥 = ／ ＼ … ～0。 ，琵髓蔷- §0 盛 既 。 (5) ， ／ ； {， 嬲霪 j囊毫囊 ll 其中各向量均为含有60 的多项式 。这一齐次式有解的 l _ - - 条件是剩余量为零 。即：△ ((I)：)=0。按一定步长在感兴 ／／ 一 ■ 、0 1 i ：＼ ： i i j ；、 j i—寸 趣 的范围内选定 60，60+A60．~o+2Ato，…为试算频率 。通 ) o 。 o30·n50 n 0．9 一_ 。0 援谤謦 薯 《 过矩阵连乘及剩余量 的计算公式 (5)式，便可求得对应 囊嚣 l|l | 于上述各试算频率 的剩余量 A ((1)z)。如果对应于相邻两试 图2 传递矩阵法计算结果 算频率 的A (60z)异号 ，则说 明在这两试算频率之间，必 有一频率能使 △ (60：)=0，即为一个满足边界条件 的 0【 值 ，也是转子的一个 同步正 向涡动角速度或临界转速。用 3．2运用有限元法对转子固有频率和振型计算 二分法仔细搜索就可 以逐步逼近求 出这一临界转速 2¨]。 为 了实现更好 的比较 ．为此运用有 限元法计算来进 3转子固有频率与振型研究 行 比较 ．首先建立有 限元模 型 ，转子长度为 1m，半径 为 0．005m．有限元模型为 beam单元 ．模型划分网格分 3．1MATJ【AB创建图形用户界面 为 20段 ，边界条件是转子 的两端为铰支 。转子材料 的 MATLAB可 以创建 图形用户界面 GUI(GraphicalUser 弹 性