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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 【求助】Riccati传递矩阵法求解临界转速（有部分程序）
北京石油化工学院2026年研究生招生接收调剂公告
61/1返回列表
查看: 1909  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288资源与环境专硕求调剂，不限专业，有学上就行 +11 lllllos 2026-03-30 11/550 2026-03-30 21:19 by 研究僧导导
[考研] 材料工程专硕求调剂 +8 hyl3153942 2026-03-29 8/400 2026-03-30 20:45 by dophin1985
[考研] 085600 材料与化工 329分求调剂 +18 Mr. Z 2026-03-25 19/950 2026-03-30 19:01 by Wang200018
[考研] 303求调剂 +7 DLkz1314. 2026-03-30 7/350 2026-03-30 16:05 by shuang5186
[考研] 085701求调剂初试286分 +5 secret0328 2026-03-28 5/250 2026-03-30 12:54 by fangnagu
[考研] 求化学调剂 +11 wulanna 2026-03-28 11/550 2026-03-30 10:59 by 探123
[考研] 296求调剂 +10 彼岸t 2026-03-29 10/500 2026-03-30 10:50 by 探123
[考研] 一志愿武汉理工，总分321，英一数二，求老师收留。 +11 nnnnnnn5 2026-03-25 11/550 2026-03-29 20:42 by 无际的草原
[考研] 321求调剂 +7 璞玉~~ 2026-03-25 8/400 2026-03-29 06:41 by 544594351
[考研] 求调剂 +6 芦lty 2026-03-25 7/350 2026-03-28 13:13 by 唐沐儿
[考研] 286求调剂 +12 PolarBear11 2026-03-26 12/600 2026-03-28 12:14 by zllcz
[考研] 0856，材料与化工321分求调剂 +12 大馋小子 2026-03-27 13/650 2026-03-28 10:56 by self2008
[考研] 308求调剂 +7 墨墨漠 2026-03-27 7/350 2026-03-28 07:43 by 热情沙漠
[考研] 265求调剂11408 +3 刘小鹿lu 2026-03-27 3/150 2026-03-27 20:53 by nihaoar
[考研] 085600，材料与化工321分调剂 +4 大馋小子 2026-03-27 6/300 2026-03-27 14:11 by 松花缸1201
[考研] 305求调剂 +5 哇卢卡库 2026-03-26 5/250 2026-03-27 14:01 by laoshidan
[硕博家园] 北京林业大学硕导招生广告 +6 kongweilin 2026-03-26 8/400 2026-03-27 10:18 by FF_16
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-26 3/150 2026-03-26 19:57 by nihaoar
[考研] 293求调剂 +7 加一一九 2026-03-24 7/350 2026-03-25 12:02 by userper
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
信息提示
关闭
请填处理意见
关闭