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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 【求助】Riccati传递矩阵法求解临界转速（有部分程序）
51/1返回列表
查看: 2223  |  回复: 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的回帖

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的回帖
查看全部 6 个回答

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的回帖

jackie_qqq

银虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我现在做的 是双转子的子结构传递矩阵法！参考下！谢谢啦
赞一下 回复此楼
5楼2012-05-02 13:48:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 18:01 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:25 by Aj1rhIDL5ixY
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 16:25 by Aj1rhIDL5ixY
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:12 by Aj1rhIDL5ixY
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 16:00 by Aj1rhIDL5ixY
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:57 by Aj1rhIDL5ixY
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:45 by Aj1rhIDL5ixY
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:37 by Aj1rhIDL5ixY
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:13 by Aj1rhIDL5ixY
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:01 by Aj1rhIDL5ixY
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:00 by Aj1rhIDL5ixY
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 10:11 by 23jxep3nCNZb
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 10:01 by 23jxep3nCNZb
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
[基金申请] 面上有专家说实验设备不是我们单位的 +6 zj_muchong 2026-09-08 8/400 2026-09-11 06:56 by suifengdao
[有机交流] 科研求助 10+4 我吃一只鱼 2026-09-07 4/200 2026-09-10 08:32 by 0940628
信息提示
关闭
请填处理意见
关闭