多情清秋 发表于 2005-5-20 16:24

关于无阻尼jeffcott转子的计算问题

本帖最后由 VibInfo 于 2016-4-15 14:41 编辑

　　

很简单的一个模型


　　

x"=-kx+lo*cos(t)


　　

y"=-kx+lo*sin(t)


　　

用Runge-kutta法计算，算了两天都没算出一个和理论解对的上的结果。


　　

包括自己编写的Runge-kutta程序以及调用V_Fortran程序库里的，结果都是一样，其轴心轨迹以一个环面的形式出现。


　　

后来在这个方程上加了一个微小的阻尼之后得出了正确的结果


　　

不知道原因何在?大家可以讨论一下

FSI 发表于 2005-5-20 19:58

[初学者必看] UG入门专辑

刚性方程

adminftp 发表于 2005-10-15 11:02

回复：（FSI） [初学者必看] UG入门专辑

本帖最后由 VibInfo 于 2016-4-15 14:41 编辑

　　

以下是引用FSI在2005-5-20 19:58:10的发言：
刚性方程

　　
什么是刚性方程?为什么刚性方程会造成这种现象?刚性方程应该如何求解?

阿H 发表于 2005-10-19 19:52

讨论：关于无阻尼jeffcott转子的计算问题

本帖最后由 VibInfo 于 2016-4-15 14:41 编辑

　　

我在做碰摩转子振动特性仿真时，也是以jeffcott转子为分析模型，以无阻尼计算时确实遇到这样的问题，当时没太注意，其实在不同的转速下，其轴心轨迹的图形不一样，更有趣的是还有一个“心”形。 有可能是在一定转速下，由于偏心的存在使得振动向某个方向发展，自然轴心轨迹也朝这个方位发展，无阻尼、两个振动方向振动向某方位增大的趋势，环面是向某个方向发展的吧! 欢迎继续讨论

frogfish 发表于 2005-10-19 21:06

回复：（FSI） [初学者必看] UG入门专辑

本帖最后由 VibInfo 于 2016-4-15 14:42 编辑

　　

以下是引用FSI在2005-5-20 19:58:10的发言：
刚性方程

　　
刚性方程的定义是

1. 特征值的实部为负数
2. 实部最大的的特征值的实部与实部最小的的特征值的实部的比值远大于1

这样的方程定义为刚性方程。

现在看楼主给出的两个方程，由于两个方程是解耦的
那么任意去一个方程，求其特征值，显见为 i*sqrt(k) -i*sqrt
显然两个条件都不满足，那怎么能称做刚性方程呢?

TNC 发表于 2005-10-20 09:44

回复：（阿H）讨论：关于无阻尼jeffcott转子的计算问题...

本帖最后由 VibInfo 于 2016-4-15 14:42 编辑

　　

以下是引用阿H在2005-10-19 19:52:00的发言：


　　

我在做碰摩转子振动特性仿真时，也是以jeffcott转子为分析模型，以无阻尼计算时确实遇到这样的问题，当时没太注意，其实在不同的转速下，其轴心轨迹的图形不一样，更有趣的是还有一个“心”形。 有可能是在一定转速下，由于偏心的存在使得振动向某个方向发展，自然轴心轨迹也朝这个方位发展，无阻尼、两个振动方向振动向某方位增大的趋势，环面是向某个方向发展的吧! 欢迎继续讨论


　　
这个应该和偏心无关，是数值分析上的问题，由于舍入误差造成的分析结果的错误
个人观点
页: [1]
查看完整版本: 关于无阻尼jeffcott转子的计算问题