 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1895 | 回復(fù): 5 | |||
chunshui2003木蟲 (著名寫手)
|
[交流]
【求助】Riccati傳遞矩陣法求解臨界轉(zhuǎn)速（有部分程序） 已有3人參與
|
|
我的專業(yè)是水輪發(fā)電機(jī)組振動，最近在研究Riccati傳遞矩陣法。但沒有人交流，自己的思路總是有限，可能到了某一個節(jié)點就繞不過去了。下面的一段程序是根據(jù)聞邦椿院士編寫的《高等轉(zhuǎn)子動力學(xué)》P11的一個算例而來，涉及到Riccati傳遞矩陣法，希望在這里能有朋友幫忙看一下是哪里出了問題。畢竟我的這個專業(yè)有些偏，在算法和編程上力學(xué)的朋友更強(qiáng)。 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個結(jié)點 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]; %不計轉(zhuǎn)動慣量和陀螺力矩 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]; %參振質(zhì)量 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)生曲線的橫坐標(biāo) grid on plot(n,x) 也許我發(fā)貼的版塊可能不太對，但我覺得有一分希望永遠(yuǎn)比沒有好。 [ Last edited by nono2009 on 2010-5-14 at 13:01 ] |
回復(fù)此樓» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有23人回復(fù)
-
265求調(diào)劑
已經(jīng)有9人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有10人回復(fù)
-
招收調(diào)劑碩士
已經(jīng)有12人回復(fù)
-
296求調(diào)劑
已經(jīng)有8人回復(fù)
-
招收博士1-2人
已經(jīng)有3人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有11人回復(fù)
-
274求調(diào)劑
已經(jīng)有8人回復(fù)
-
290求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求解行列式。。。
已經(jīng)有3人回復(fù)
- 病態(tài)矩陣的精確解求解
已經(jīng)有17人回復(fù)
- 求助轉(zhuǎn)動慣量的問題
已經(jīng)有15人回復(fù)
- 下面的程序哪里有問題
已經(jīng)有7人回復(fù)
- 求解ＡＢＣＤ矩陣傳輸理論
已經(jīng)有4人回復(fù)
- 需要利用Maple求解的程序
已經(jīng)有6人回復(fù)
- 關(guān)于轉(zhuǎn)動慣量的問題
已經(jīng)有11人回復(fù)
- 求助一個矩陣求解問題
已經(jīng)有11人回復(fù)
- 關(guān)于多體系統(tǒng)傳遞矩陣法的問題
已經(jīng)有9人回復(fù)
- 一個矩陣證明題，求助！
已經(jīng)有4人回復(fù)
- matlab 求解riccati 方程
已經(jīng)有4人回復(fù)
- 矩陣相乘，逆向求解，這種方法可行嗎？
已經(jīng)有8人回復(fù)
- 求能得到躍遷矩陣元的程序
已經(jīng)有4人回復(fù)
- 【求助/交流】500g 是多大的轉(zhuǎn)速
已經(jīng)有3人回復(fù)
- 【求助】臨界相對濕度程序
已經(jīng)有16人回復(fù)
- 【求助】求精餾塔理論板求解程序，簡捷法和逐板法的應(yīng)用程序
已經(jīng)有4人回復(fù)
- 【求助】關(guān)于點估計和矩估計法
已經(jīng)有4人回復(fù)
luomingqi
木蟲 (正式寫手)
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 4356.8
- 散金: 20
- 紅花: 1
- 帖子: 870
- 在線: 66.5小時
- 蟲號: 896596
- 注冊: 2009-11-07
- 性別: GG
- 專業(yè): 控制理論與方法
2樓2010-05-12 22:13:35
chunshui2003
木蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 16242.9
- 散金: 1992
- 紅花: 14
- 帖子: 1387
- 在線: 168小時
- 蟲號: 1000957
- 注冊: 2010-04-19
- 專業(yè): 結(jié)構(gòu)工程
3樓2010-05-13 12:11:38
bluesine
鐵桿木蟲 (職業(yè)作家)
科苑小木蟲
- 應(yīng)助: 132 (高中生)
- 貴賓: 1.991
- 金幣: 9844.6
- 散金: 89
- 紅花: 19
- 帖子: 3618
- 在線: 377.4小時
- 蟲號: 869544
- 注冊: 2009-10-12
- 性別: GG
- 專業(yè): 數(shù)學(xué)物理
4樓2010-05-13 17:08:54
jackie_qqq
銀蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 317
- 帖子: 19
- 在線: 8.9小時
- 蟲號: 1788854
- 注冊: 2012-05-01
- 專業(yè): 金屬材料的凝固與結(jié)晶學(xué)
5樓2012-05-02 13:48:16
WGY348
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2051.6
- 散金: 86
- 紅花: 1
- 帖子: 261
- 在線: 37.2小時
- 蟲號: 3503349
- 注冊: 2014-10-28
- 性別: GG
- 專業(yè): 機(jī)械動力學(xué)
6樓2014-10-28 16:44:59
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 chunshui2003 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 279分求調(diào)劑 一志愿211 +9 | chaojifeixia 2026-03-19 | 10/500 |
|
|---|---|---|---|---|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +15 | 于海海海海 2026-03-19 | 15/750 |
|
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +8 | zlingli 2026-03-13 | 8/400 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|