 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1896 | 回復(fù): 5 | |||
chunshui2003木蟲 (著名寫手)
|
[交流]
【求助】Riccati傳遞矩陣法求解臨界轉(zhuǎn)速（有部分程序） 已有3人參與
|
|
我的專業(yè)是水輪發(fā)電機(jī)組振動(dòng)，最近在研究Riccati傳遞矩陣法。但沒有人交流，自己的思路總是有限，可能到了某一個(gè)節(jié)點(diǎn)就繞不過去了。下面的一段程序是根據(jù)聞邦椿院士編寫的《高等轉(zhuǎn)子動(dòng)力學(xué)》P11的一個(gè)算例而來，涉及到Riccati傳遞矩陣法，希望在這里能有朋友幫忙看一下是哪里出了問題。畢竟我的這個(gè)專業(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個(gè)結(jié)點(diǎn) 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]; %不計(jì)轉(zhuǎn)動(dòng)慣量和陀螺力矩 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]; %不計(jì)剪切影響 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ā)貼的版塊可能不太對(duì)，但我覺得有一分希望永遠(yuǎn)比沒有好。 [ Last edited by nono2009 on 2010-5-14 at 13:01 ] |
回復(fù)此樓» 猜你喜歡
317求調(diào)劑
已經(jīng)有12人回復(fù)
-
一志愿武漢理工材料工程專碩調(diào)劑
已經(jīng)有9人回復(fù)
-
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑
已經(jīng)有7人回復(fù)
-
324求調(diào)劑
已經(jīng)有5人回復(fù)
-
288求調(diào)劑
已經(jīng)有16人回復(fù)
-
317求調(diào)劑
已經(jīng)有9人回復(fù)
-
287求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。
已經(jīng)有4人回復(fù)
-
求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有4人回復(fù)
-
085600材料與化工
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求解行列式。。。
已經(jīng)有3人回復(fù)
- 病態(tài)矩陣的精確解求解
已經(jīng)有17人回復(fù)
- 求助轉(zhuǎn)動(dòng)慣量的問題
已經(jīng)有15人回復(fù)
- 下面的程序哪里有問題
已經(jīng)有7人回復(fù)
- 求解ＡＢＣＤ矩陣傳輸理論
已經(jīng)有4人回復(fù)
- 需要利用Maple求解的程序
已經(jīng)有6人回復(fù)
- 關(guān)于轉(zhuǎn)動(dòng)慣量的問題
已經(jīng)有11人回復(fù)
- 求助一個(gè)矩陣求解問題
已經(jīng)有11人回復(fù)
- 關(guān)于多體系統(tǒng)傳遞矩陣法的問題
已經(jīng)有9人回復(fù)
- 一個(gè)矩陣證明題，求助！
已經(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ù)
- 【求助】臨界相對(duì)濕度程序
已經(jīng)有16人回復(fù)
- 【求助】求精餾塔理論板求解程序，簡(jiǎn)捷法和逐板法的應(yīng)用程序
已經(jīng)有4人回復(fù)
- 【求助】關(guān)于點(diǎn)估計(jì)和矩估計(jì)法
已經(jīng)有4人回復(fù)
luomingqi
木蟲 (正式寫手)
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 4356.8
- 散金: 20
- 紅花: 1
- 帖子: 870
- 在線: 66.5小時(shí)
- 蟲號(hào): 896596
- 注冊(cè): 2009-11-07
- 性別: GG
- 專業(yè): 控制理論與方法
2樓2010-05-12 22:13:35
chunshui2003
木蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 16242.9
- 散金: 1992
- 紅花: 14
- 帖子: 1387
- 在線: 168小時(shí)
- 蟲號(hào): 1000957
- 注冊(cè): 2010-04-19
- 專業(yè): 結(jié)構(gòu)工程
3樓2010-05-13 12:11:38
bluesine
鐵桿木蟲 (職業(yè)作家)
科苑小木蟲
- 應(yīng)助: 132 (高中生)
- 貴賓: 1.991
- 金幣: 9834.6
- 散金: 89
- 紅花: 19
- 帖子: 3619
- 在線: 377.5小時(shí)
- 蟲號(hào): 869544
- 注冊(cè): 2009-10-12
- 性別: GG
- 專業(yè): 數(shù)學(xué)物理
4樓2010-05-13 17:08:54
jackie_qqq
銀蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 317
- 帖子: 19
- 在線: 8.9小時(shí)
- 蟲號(hào): 1788854
- 注冊(cè): 2012-05-01
- 專業(yè): 金屬材料的凝固與結(jié)晶學(xué)
5樓2012-05-02 13:48:16
WGY348
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2051.6
- 散金: 86
- 紅花: 1
- 帖子: 261
- 在線: 37.2小時(shí)
- 蟲號(hào): 3503349
- 注冊(cè): 2014-10-28
- 性別: GG
- 專業(yè): 機(jī)械動(dòng)力學(xué)
6樓2014-10-28 16:44:59
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 chunshui2003 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！ +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 14/700 |
|
|
|
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 5/250 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|