奔腾6 发表于 2007-7-22 20:22

Riccati传递矩阵

   各位师兄师姐：大家好，很抱歉打扰你们，我最近开始做课题，是转子动力学方面的，用matlab编程求临界转速和振型，是按照闻邦椿《高等转子动力学》上面的Riccati传递矩阵公式编的，算例也是用上面的数据，但是得出的剩余量曲线和书上的相差很大，反复检查好长时间，就是不知哪里的问题，很头疼啊，只能向前辈请教了，请帮忙看一看，或者有正确的程序给我发一个作为参考，多谢，祝身体健康。
clear;
l=;      %13个结点
m=;
Jp=;
Jd=;       %不计转动惯量和陀螺力矩
I=;
E=4393;
v=;   %不计剪切影响
k=;    %油膜刚度
kb=;       %轴承座刚度
mb=;   %参振质量
S=;    %Riccati第一矩阵
s=1;
x=[];   %用于记录剩余量
for n=1864:1864       %试算频率
for i=1:13
K=k(i)*(kb(i)-mb(i)*n^2)/(k(i)+kb(i)-mb(i)*n^2)      %总刚度
u11=;
u12=;
u21=(l(i)/(E*I(i))).*;
u22=;
S=*inv()
end
D=det(S);
x=;      %记录剩余量
end
n=1864:1:1864;   %产生曲线的横坐标
grid on
plot(n,x)

sssssxxxxx921 发表于 2007-7-26 19:34

什么年代了怎么还用这种方法做啊找九十年代的研究生问吧

奔腾6 发表于 2007-7-28 10:09

回复 #2 sssssxxxxx921 的帖子

那用什么方法？？？

hector1982 发表于 2007-7-30 18:19

看＜转子动力学＞吧，
用传递矩阵法吧
参考
http://forum.vibunion.com/forum/viewthread.php?tid=44831&extra=page%3D1%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D157

yangyongwhu@163 发表于 2010-4-16 13:48

Riccati传递矩阵
比普通的传递矩阵法更为优化啊，你们不知道的吗？Riccati传递矩阵
法就是在普通传递矩阵法的基础上修改而来的

lhy 发表于 2010-12-20 22:15

本帖最后由 lhy 于 2010-12-20 22:16 编辑

Riccati传递矩阵法及其改进算法在进行振动分析时，只需对一些阶次很低的传递矩阵进行连续的矩阵乘法运算，并且和系统自由度没有直接关系，大大节省计算工作量，特别是在只要求计算若干阶低阶的固有频率和振型时更加方便。Riccati传递矩阵法的基本思想是：把系统分割成一系列具有简单动力学特性的单元（两端，三端或者多端），振动时系统的状态可以用该单元端点的状态矢量来表示，相邻单元之间状态矢量的关系，用传递矩阵表示。

yejet 发表于 2011-1-11 11:04

sssssxxxxx921 发表于 2007-7-26 19:34 static/image/common/back.gif
什么年代了怎么还用这种方法做啊找九十年代的研究生问吧

看来你已经很高级了，不知你用什么方法呢

yejet 发表于 2011-1-11 11:17

本帖最后由 yejet 于 2011-1-11 11:18 编辑

一个例子，仅作参考，比保证程序的正确性clc
clear
format short e
m1=2940;m13=2940;m2=5880;
M=;
l=1.3;
L=;
a=2.9592e-9;
A=;
k1=1.960000000e9;
k2=2.7048e9;
mb=3577;
k=0;
S=;
for w=0:0.01:5000
    for i=1:13
       k3=1.0e9;
       K=;
    U11=;
    U12=;
    U21=;
    U22=;
    S=(U11*S+U12)/(U21*S+U22);
   end
   F=S(1,1)*S(2,2)-S(1,2)*S(2,1);
if F*(-1)^k< 0 %求解临界转速
       k=k+1;
      wi(k)=w;
       w=wi(k)
       ni(k)=wi(k)*30/pi
    end

end

huangziyuan212 发表于 2014-6-5 16:22

{:{05}:}

huangziyuan212 发表于 2014-6-17 10:17

单元长度为0是什么意思啊？？？

璀璨星辰 发表于 2014-6-20 17:03

{:{39}:}

huangziyuan212 发表于 2014-6-30 08:59

{:{39}:}

dollfish000 发表于 2014-11-1 00:34

二楼挺搞笑，我们要相互学习，取长补短！
页: [1]
查看完整版本: Riccati传递矩阵