 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1908 | 回復: 5 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
chunshui2003木蟲 (著名寫手)
|
[交流]
【求助】Riccati傳遞矩陣法求解臨界轉速（有部分程序） 已有3人參與
|
||
|
我的專業(yè)是水輪發(fā)電機組振動，最近在研究Riccati傳遞矩陣法。但沒有人交流，自己的思路總是有限，可能到了某一個節(jié)點就繞不過去了。下面的一段程序是根據(jù)聞邦椿院士編寫的《高等轉子動力學》P11的一個算例而來，涉及到Riccati傳遞矩陣法，希望在這里能有朋友幫忙看一下是哪里出了問題。畢竟我的這個專業(yè)有些偏，在算法和編程上力學的朋友更強。 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; %產(chǎn)生曲線的橫坐標 grid on plot(n,x) 也許我發(fā)貼的版塊可能不太對，但我覺得有一分希望永遠比沒有好。 [ Last edited by nono2009 on 2010-5-14 at 13:01 ] |
回復此樓» 猜你喜歡
0703化學
已經(jīng)有12人回復
-
材料與化工272求調劑
已經(jīng)有19人回復
-
一志愿河北工業(yè)大學0817化工278分求調劑
已經(jīng)有16人回復
-
317求調劑
已經(jīng)有8人回復
-
332求調劑
已經(jīng)有7人回復
-
材料與化工304求B區(qū)調劑
已經(jīng)有7人回復
-
085404求調劑，總分309，本科經(jīng)歷較為豐富
已經(jīng)有5人回復
-
332求調劑
已經(jīng)有3人回復
-
318求調劑
已經(jīng)有3人回復
-
295求調劑
已經(jīng)有5人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 求解行列式。。。
已經(jīng)有3人回復
- 病態(tài)矩陣的精確解求解
已經(jīng)有17人回復
- 求助轉動慣量的問題
已經(jīng)有15人回復
- 下面的程序哪里有問題
已經(jīng)有7人回復
- 求解ＡＢＣＤ矩陣傳輸理論
已經(jīng)有4人回復
- 需要利用Maple求解的程序
已經(jīng)有6人回復
- 關于轉動慣量的問題
已經(jīng)有11人回復
- 求助一個矩陣求解問題
已經(jīng)有11人回復
- 關于多體系統(tǒng)傳遞矩陣法的問題
已經(jīng)有9人回復
- 一個矩陣證明題，求助！
已經(jīng)有4人回復
- matlab 求解riccati 方程
已經(jīng)有4人回復
- 矩陣相乘，逆向求解，這種方法可行嗎？
已經(jīng)有8人回復
- 求能得到躍遷矩陣元的程序
已經(jīng)有4人回復
- 【求助/交流】500g 是多大的轉速
已經(jīng)有3人回復
- 【求助】臨界相對濕度程序
已經(jīng)有16人回復
- 【求助】求精餾塔理論板求解程序，簡捷法和逐板法的應用程序
已經(jīng)有4人回復
- 【求助】關于點估計和矩估計法
已經(jīng)有4人回復
chunshui2003
木蟲 (著名寫手)
- 應助: 0 (幼兒園)
- 金幣: 16312.9
- 散金: 1992
- 紅花: 14
- 帖子: 1397
- 在線: 168小時
- 蟲號: 1000957
- 注冊: 2010-04-19
- 專業(yè): 結構工程
3樓2010-05-13 12:11:38
luomingqi
木蟲 (正式寫手)
- 應助: 47 (小學生)
- 金幣: 4356.8
- 散金: 20
- 紅花: 1
- 帖子: 870
- 在線: 66.5小時
- 蟲號: 896596
- 注冊: 2009-11-07
- 性別: GG
- 專業(yè): 控制理論與方法
2樓2010-05-12 22:13:35
bluesine
鐵桿木蟲 (職業(yè)作家)
科苑小木蟲
- 應助: 132 (高中生)
- 貴賓: 1.991
- 金幣: 9850.4
- 散金: 89
- 紅花: 19
- 帖子: 3623
- 在線: 377.5小時
- 蟲號: 869544
- 注冊: 2009-10-12
- 性別: GG
- 專業(yè): 數(shù)學物理
4樓2010-05-13 17:08:54
jackie_qqq
銀蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 317
- 帖子: 19
- 在線: 8.9小時
- 蟲號: 1788854
- 注冊: 2012-05-01
- 專業(yè): 金屬材料的凝固與結晶學
5樓2012-05-02 13:48:16
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程專碩求調劑 +3 | hyl3153942 2026-03-29 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 281求調劑 +4 | 亞克西good 2026-03-26 | 6/300 |
|
|
|
[考研] 283求調劑（080500） +7 | A child 2026-03-27 | 7/350 |
|
|
|
[考研] 一志愿北京化工大學材料與化工（085600）296求調劑 +17 | 稻妻小編 2026-03-26 | 17/850 |
|
|
|
[考研] 329求調劑 +10 | 鈕恩雪 2026-03-25 | 10/500 |
|
|
|
[考研] 291求調劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
|
[考研] 一志愿211 初試270分 求調劑 +8 | 谷雨上岸 2026-03-23 | 9/450 |
|
|
|
[考研] 一志愿太原理工安全工程300分，求調劑 +5 | 0857求調劑. 2026-03-24 | 6/300 |
|
|
|
[考研] 343求調劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
|
[考研] 083000學碩274求調劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|
|
[考研] 299求調劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 材料277求調劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 291求調劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[考研] 291求調劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
|
[考研] 331環(huán)境科學與工程求調劑 +3 | 熠然好運氣 2026-03-27 | 3/150 |
|
|
|
[考研] 調劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[考研] 359求調劑 +4 | 王了個楠 2026-03-25 | 4/200 |
|
|
|
[考研] 342求調劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
|
[考研] 材料專碩331求調劑 +4 | 鮮當牛 2026-03-24 | 4/200 |
|