wangli20060913 发表于 2006-9-23 15:12

请教：扭转振动的圆频率问题

在扭转振动分析中，遇到如下问题：
      转动惯量J的数量级为e+(-2),刚度K的数量级为e+10,运用MATLAB运行eig(inv(J)*K)后得到的圆频率数量级为e+5，这个结果是不可能的阿。我检查了转动惯量和刚度计算，没有问题，问题出在哪里，查不出，请高手执教。谢谢！

xinyuxf 发表于 2006-9-23 15:27

最好能把部分代码也贴出来，这样大家也好帮你找找原因，你这样说的有点简单了，大家可能会不明白

wangli20060913 发表于 2006-9-23 19:12

请教：扭转振动的圆频率问题

谢谢！
       对于扭转自由振动方程： JX''+KX=0
      
   J=
    K=
除了第一行和最后一行只有两个非零值，其他都是只有连续三个非零值，转动惯量很小，刚度很大，我参考别人的解法eig(inv(J)*K)结果其圆频率是几百，而我的结果是几万多，不知为什么，请高手指教，谢谢!!!

puyen 发表于 2006-9-26 17:19

对于求解 MX''+KX=0的固有频率matlab的命令是
E=eig(K,M);

w=sqrt(E)
其中w为固有频率

wangli20060913 发表于 2006-9-27 10:49

谢谢puyen!!!!

p123456q 发表于 2006-9-28 15:32

V=eig(inv(J)*K)与=eig(K,J)是一样的，求解出的V再开平方，就得出圆频率了。
页: [1]
查看完整版本: 请教：扭转振动的圆频率问题