| 查看: 2177 | 回复: 5 | |||
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 ] |
回复此楼» 猜你喜欢
能否退出参与的面上项目解除限项
已经有18人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有3人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
-
我面上完蛋了
已经有9人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求解行列式。。。
已经有3人回复
- 病态矩阵的精确解求解
已经有17人回复
- 求助转动惯量的问题
已经有15人回复
- 下面的程序哪里有问题
已经有7人回复
- 求解ＡＢＣＤ矩阵传输理论
已经有4人回复
- 需要利用Maple求解的程序
已经有6人回复
- 关于转动惯量的问题
已经有11人回复
- 求助一个矩阵求解问题
已经有11人回复
- 关于多体系统传递矩阵法的问题
已经有9人回复
- 一个矩阵证明题，求助！
已经有4人回复
- matlab 求解riccati 方程
已经有4人回复
- 矩阵相乘，逆向求解，这种方法可行吗？
已经有8人回复
- 求能得到跃迁矩阵元的程序
已经有4人回复
- 【求助/交流】500g 是多大的转速
已经有3人回复
- 【求助】临界相对湿度程序
已经有16人回复
- 【求助】求精馏塔理论板求解程序，简捷法和逐板法的应用程序
已经有4人回复
- 【求助】关于点估计和矩估计法
已经有4人回复
luomingqi
木虫 (正式写手)
- 应助: 47 (小学生)
- 金币: 4356.8
- 散金: 20
- 红花: 1
- 帖子: 870
- 在线: 66.5小时
- 虫号: 896596
- 注册: 2009-11-07
- 性别: GG
- 专业: 控制理论与方法
2楼2010-05-12 22:13:35
chunshui2003
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 17587.7
- 散金: 1992
- 红花: 14
- 帖子: 1546
- 在线: 168小时
- 虫号: 1000957
- 注册: 2010-04-19
- 专业: 结构工程
3楼2010-05-13 12:11:38
bluesine
至尊木虫 (职业作家)
科苑小木虫
- 应助: 132 (高中生)
- 贵宾: 1.991
- 金币: 10425.7
- 散金: 89
- 红花: 19
- 帖子: 3803
- 在线: 382.1小时
- 虫号: 869544
- 注册: 2009-10-12
- 性别: GG
- 专业: 数学物理
4楼2010-05-13 17:08:54
5楼2012-05-02 13:48:16
WGY348
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2051.6
- 散金: 86
- 红花: 1
- 帖子: 261
- 在线: 37.2小时
- 虫号: 3503349
- 注册: 2014-10-28
- 性别: GG
- 专业: 机械动力学
6楼2014-10-28 16:44:59
20