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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 【求助】Riccati传递矩阵法求解临界转速（有部分程序）
61/1返回列表
查看: 2132  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 jRl3mE6ddZGq 2026-08-03 4/200 2026-08-04 04:54 by k7OM8YghWbkC
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 04:50 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:46 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 5/250 2026-08-04 04:44 by k7OM8YghWbkC
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:41 by k7OM8YghWbkC
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 by k7OM8YghWbkC
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[基金申请] 微信指数没变化，科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
关闭
请填处理意见
关闭