【论坛出品】转子系统临界转速Matlab计算程序（1.1.0版）

jackyyjy · 发表于 2013-1-21 14:42

本帖最后由 jackyyjy 于 2013-1-21 14:44 编辑

为什么我用ANSYS算出来的算例一，结果不一样啊？？？？
算例一：如图所示两端简支的光轴系统，轴长为2m，轴直径为0.1m，转子弹性模量为2.095e11Pa，转子材料密度为7.85e3kg/m3。
用论坛程序算出该转子可以求得临界转速的理论解，理论解为：
一阶临界转速：3043 转/分，换算成固有频率除以60在50HZ左右
但是我用ANSYS计算，两边固定支撑，结果是111HZ为什么？？？为什么？？？

111HZ啊。。。

茹果刘儿 · 发表于 2013-3-14 10:27

怎么研究？也没有源代码。。。如何用riccati 传递矩阵法呢？编程总是出错误，求指点。。。

z308725721 · 发表于 2013-4-16 17:52

目前正在学习这方面的东西，收藏着先，感觉很有用啊

c39850477 · 发表于 2013-4-17 23:36

多谢高手

茹果刘儿 · 发表于 2013-5-14 10:31

本帖最后由茹果刘儿于 2013-5-14 10:46 编辑

谢谢谢谢谢谢！！！！楼主可以让我们看看源文件吗？

mashimaro001 · 发表于 2013-5-15 01:41

请问楼主，是使用的传递矩阵法，还是使用的有限元法求的临界转速呢？

风驰涧 · 发表于 2013-5-23 15:54

灰常喜欢呀！！1

ME! · 发表于 2013-6-3 22:24

论坛出品发表于 2012-11-21 10:51

登录/注册后可看大图

程序暂时没有考虑陀螺力矩的功能

请问是通过单元矩阵，然后迭代运行的吧，我想问的是你用的是求特征值的matlab函数eig吗？我想考虑陀螺效应，然后该怎么做？求解钟一谔 p186面的特征值方程

ME! · 发表于 2013-6-3 22:24

论坛出品发表于 2012-11-21 10:51

登录/注册后可看大图

程序暂时没有考虑陀螺力矩的功能

请问是通过单元矩阵，然后迭代运行的吧，我想问的是你用的是求特征值的matlab函数eig吗？我想考虑陀螺效应，然后该怎么做？求解钟一谔 p186面的特征值方程

ME! · 发表于 2013-6-3 22:37

论坛出品发表于 2008-10-13 14:33

登录/注册后可看大图

算例二：如图所示转子，转子具体尺寸见程序中的输入参数，转子弹性模量为2.095e11Pa，转子材料密度为7.85e3 ...

请问我没有附加轮盘，但是我想将装在轴上的零件作为附加质量添加进去，应该怎么做？

ME! · 发表于 2013-6-4 21:45

论坛出品发表于 2008-10-13 14:32

登录/注册后可看大图

算例一：如图所示两端简支的光轴系统，轴长为2m，轴直径为0.1m，转子弹性模量为2.095e11Pa，转子材料密度为 ...

我很疑惑理论解不是简支梁的理论公式吗？
w=(pi)^2*sqrt(2.095e11*pi*(0.2)^4/64/(7850*pi*0.1^2*0.25^4)) %理论解
n=30*w/pi
结果和楼主的不同

ME! · 发表于 2013-6-4 21:49

maosuizijian 发表于 2009-10-31 11:36

登录/注册后可看大图

{很显然，该转子可以求得临界转速的理论解，理论解为：
一阶临界转速：3043 转/分

同问，我求出来的理论解不同？

ME! · 发表于 2013-6-4 22:17

jackyyjy 发表于 2012-12-12 16:16

登录/注册后可看大图

是啊,这个P文件应该怎么调用啊，调用不了。。。

补充内容 (2013-1-5 09:42):

我用ansys计算了一下第一个简支梁的模型，得出的固有频率也和楼主不一样

finish
/clear
/prep7

RO = 7850
PEX = 2.095e11
MP,EX,1,PEX !ex: 弹性模量
MP,DENS,1,RO  !密度
MP,GXY, 1,EX/2.6
MP,PRXY,1,0.3
ET,1,beam3,,,
R,1,0.007854,4.9087e-6

N,1,0
N,2,0.25
N,3,0.25*2
N,4,0.25*3
N,5,0.25*4
N,6,0.25*5
N,7,0.25*6
N,8,0.25*7
N,9,0.25*8

TYPE,1
E,1,2
E,2,3
E,3,4
E,4,5
E,5,6
E,6,7
E,7,8
E,8,9

finish

/SOLU
D,1,ALL
D,9,ALL

/SOLU
ANTYPE,MODAL
MODOPT,lanb,10,,,ON
MXPAND,10
SOLVE
finish

ME! · 发表于 2013-6-4 22:25

本帖最后由 ME! 于 2013-6-4 22:32 编辑

jackyyjy 发表于 2013-1-21 14:42

登录/注册后可看大图

为什么我用ANSYS算出来的算例一，结果不一样啊？？？？
算例一：如图所示两端简支的光轴系统，轴长 ...

我用ansys算了下，结果前三阶固有频率为       1  114.87HZ          1       1       1
   2  315.98HZ          1       2       2
   3  618.00HZ          1       3       3
第一阶临界转速=114.87*60用matlab计算了下，前三阶固有频率和临界转速为
f =  1.0e+002 *
1.149803836884732
3.171213383661561
6.227643357869332

n =  1.0e+004 *
0.689882302130839
1.902728030196937
3.736586014721600
所以我现在质疑楼主的程序了

茹果刘儿 · 发表于 2013-7-24 15:04

好东西~~谢谢~~~

账号		自动登录	找回密码
密码			我要加入

[转子动力学] 【论坛出品】转子系统临界转速Matlab计算程序（1.1.0版）