应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 【求助】Riccati传递矩阵法求解临界转速（有部分程序）
北京石油化工学院2026年研究生招生接收调剂公告
61/1返回列表
查看: 1914  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

chunshui2003

木虫 (著名写手)

[交流] 【求助】Riccati传递矩阵法求解临界转速（有部分程序） 已有3人参与

我的专业是水轮发电机组振动，最近在研究Riccati传递矩阵法。但没有人交流，自己的思路总是有限，可能到了某一个节点就绕不过去了。下面的一段程序是根据闻邦椿院士编写的《高等转子动力学》P11的一个算例而来，涉及到Riccati传递矩阵法，希望在这里能有朋友帮忙看一下是哪里出了问题。毕竟我的这个专业有些偏，在算法和编程上力学的朋友更强。


l=[1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,1.3,0];      %13个结点
m=[2940,5880,5880,5880,5880,5880,5880,5880,5880,5880,5880,5880,2940];
Jp=[0,0,0,0,0,0,0,0,0,0,0,0,0];
Jd=[0,0,0,0,0,0,0,0,0,0,0,0,0];       %不计转动惯量和陀螺力矩
I=[100,100,100,100,100,100,100,100,100,100,100,100,100];
E=4393;
v=[0,0,0,0,0,0,0,0,0,0,0,0,0];     %不计剪切影响
k=[1.96*10^9,0,0,1.96*10^9,0,0,1.96*10^9,0,0,1.96*10^9,0,0,1.96*10^9];    %油膜刚度
kb=[2.7048*10^9,0,0,2.7048*10^9,0,0,2.7048*10^9,0,0,2.7048*10^9,0,0,2.7048*10^9];       %轴承座刚度
mb=[3577,3577,3577,3577,3577,3577,3577,3577,3577,3577,3577,3577,3577];     %参振质量
S=[0,0;0,0];    %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=[1,l(i);0,1];
u12=[l(i)*(m(i)*n^2-K),(Jp(i)-Jd(i))*n^2;m(i)*n^2-K,0];
u21=(l(i)/(E*I(i))).*[l(i)/2,l(i)^2*(1-v(i))/6;1,l(i)/2];
u22=[1+l(i)^3*(1-v(i))*(m(i)*n^2-K)/(6*E*I(i)),l(i)+l(i)^2*(Jp(i)-Jd(i))*n^2/(2*E*I(i));l(i)^2*(m(i)*n^2-K)/(2*E*I(i)),1+l(i)*(Jp(i)-Jd(i))*n^2/(E*I(i))];
S=[u11*S+u12]*inv([u21*S+u22])
end
D=det(S);
x=[x,D];        %记录剩余量
end
n=1864:1:1864;     %产生曲线的横坐标
grid on
plot(n,x)

也许我发贴的版块可能不太对，但我觉得有一分希望永远比没有好。

[ Last edited by nono2009 on 2010-5-14 at 13:01 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-05-12 19:51:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luomingqi

木虫 (正式写手)
★
小雨萌萌(金币+1):辛苦奖，鼓励一下 2010-05-12 22:22:56
不是很理解的，我搞的是线性系统的，常用MATLAB的，但是你 的这个东西，我怎么的都没有看懂的啊
赞一下(1人) 回复此楼
跟踪
2楼2010-05-12 22:13:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chunshui2003

木虫 (著名写手)
这是求解一个轴系统临界转速的程序，用的是传递矩阵法，Riccati的传递矩阵法。如果你要是搞转子应该知道。上面的程序是我在别人那里获得的，但他说有问题。我觉得也是，主要是建立“S”之后，我不太清楚如果求转速下一步的程序如何。我表达的可能不清楚，如果你有《高等转子动力学》看看就明白了。
赞一下(1人) 回复此楼
3楼2010-05-13 12:11:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫
nono2009:下次让楼主重新发帖比较好。转贴容易引起版面管理上的问题。谢谢。 2010-05-14 13:02:17
chunshui2003(金币+1):谢谢bluesine的建议，刚来本论坛不久，希望以后能和大家交流，我的专业水轮发电机组振动，也涉及到大量的非线性方面内容。 2010-05-15 12:26:24
编程的话请转到程序板块。
赞一下(1人) 回复此楼
板凳要做十年冷文章不发一个字
4楼2010-05-13 17:08:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jackie_qqq

银虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我现在做的 是双转子的子结构传递矩阵法！参考下！谢谢啦
赞一下 回复此楼
5楼2012-05-02 13:48:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WGY348

铜虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
您好，我也是研究转子动力学的，您的这个程序还在吗？我能看一下吗？
赞一下 回复此楼
只做第一个我
6楼2014-10-28 16:44:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chunshui2003 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本科211总分289，08工学真心求调剂 +3 utopiaE 2026-03-30 3/150 2026-03-30 23:42 by ms629
[考研] 11408总分309，一志愿东南大学求调剂，不挑专业 +5 天赋带到THU 2026-03-29 6/300 2026-03-30 20:49 by dick_runner
[考研] 286分调剂 +10 Faune 2026-03-30 11/550 2026-03-30 20:35 by 啊李999
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
[考研] 322求调剂：一志愿湖南大学 材料与化工（085600），已过六级。 +9 XX小邓 2026-03-29 9/450 2026-03-30 17:18 by limeifeng
[考研] 284求调剂 +14 junqihahaha 2026-03-26 15/750 2026-03-30 14:12 by 探123
[考研] 070300化学354求调剂 +7 101次希望 2026-03-28 7/350 2026-03-30 13:22 by Iveryant
[考研] 085701求调剂初试286分 +5 secret0328 2026-03-28 5/250 2026-03-30 12:54 by fangnagu
[考研] 283求调剂（080500） +14 A child 2026-03-27 14/700 2026-03-30 12:06 by 探123
[考研] 318求调剂 +7 陈晨79 2026-03-30 7/350 2026-03-30 10:49 by 探123
[考研] 329求调剂 +10 钮恩雪 2026-03-25 10/500 2026-03-29 13:32 by peike
[考研] 352分-085602-一志愿985 +5 海纳百川Ly 2026-03-29 5/250 2026-03-29 09:57 by Sjndkwm
[考研] 一志愿太原理工安全工程300分，求调剂 +5 0857求调剂. 2026-03-24 6/300 2026-03-28 22:04 by zhq0425
[考研] 材料求调剂一志愿哈工大324 +7 闫旭东 2026-03-28 9/450 2026-03-28 08:51 by Xu de nuo
[考研] 0703化学求调剂，各位老师看看我！！！ +5 祁祺祺 2026-03-25 5/250 2026-03-27 21:44 by 东方猪猪
[考研] 351求调剂 +4 麦克阿磊 2026-03-24 4/200 2026-03-27 00:32 by wxiongid
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
[考研] 302求调剂 +4 锦衣卫藤椒 2026-03-25 4/200 2026-03-25 16:29 by 功夫疯狂
[考研] 各位老师您好：本人初试372分 +5 jj涌77 2026-03-25 6/300 2026-03-25 14:15 by mapenggao
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
信息提示
关闭
请填处理意见
关闭