版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖

chunshui2003

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 16312.9
散金: 1992
紅花: 14
帖子: 1397
在線(xiàn): 168小時(shí)
蟲(chóng)號(hào): 1000957
注冊(cè): 2010-04-19
專(zhuān)業(yè): 結(jié)構(gòu)工程

[交流] 【求助】Riccati傳遞矩陣法求解臨界轉(zhuǎn)速（有部分程序）已有3人參與

我的專(zhuān)業(yè)是水輪發(fā)電機(jī)組振動(dòng)，最近在研究Riccati傳遞矩陣法。但沒(méi)有人交流，自己的思路總是有限，可能到了某一個(gè)節(jié)點(diǎn)就繞不過(guò)去了。下面的一段程序是根據(jù)聞邦椿院士編寫(xiě)的《高等轉(zhuǎn)子動(dòng)力學(xué)》P11的一個(gè)算例而來(lái)，涉及到Riccati傳遞矩陣法，希望在這里能有朋友幫忙看一下是哪里出了問(wèn)題。畢竟我的這個(gè)專(zhuān)業(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)生曲線(xiàn)的橫坐標(biāo)
grid on
plot(n,x)

也許我發(fā)貼的版塊可能不太對(duì)，但我覺(jué)得有一分希望永遠(yuǎn)比沒(méi)有好。

[ Last edited by nono2009 on 2010-5-14 at 13:01 ]

回復(fù)此樓

» 猜你喜歡

一志愿大連理工大學(xué)材料求調(diào)劑已經(jīng)有5人回復(fù)
材料與化工調(diào)劑一志愿大連海事085600，349 已經(jīng)有6人回復(fù)
一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
材料科學(xué)與工程調(diào)劑已經(jīng)有11人回復(fù)
311求調(diào)劑一志愿合肥工業(yè)大學(xué) 已經(jīng)有4人回復(fù)
一志愿食品科學(xué)與工程083200求調(diào)劑已經(jīng)有4人回復(fù)
材料科學(xué)與工程求調(diào)劑已經(jīng)有10人回復(fù)
一志愿北化085600材料專(zhuān)碩275|有文章專(zhuān)利｜求調(diào)劑已經(jīng)有11人回復(fù)
一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩已經(jīng)有4人回復(fù)
085600材料與化工調(diào)劑已經(jīng)有15人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求解行列式。。。已經(jīng)有3人回復(fù)
病態(tài)矩陣的精確解求解已經(jīng)有17人回復(fù)
求助轉(zhuǎn)動(dòng)慣量的問(wèn)題已經(jīng)有15人回復(fù)
下面的程序哪里有問(wèn)題已經(jīng)有7人回復(fù)
求解ＡＢＣＤ矩陣傳輸理論已經(jīng)有4人回復(fù)
需要利用Maple求解的程序已經(jīng)有6人回復(fù)
關(guān)于轉(zhuǎn)動(dòng)慣量的問(wèn)題已經(jīng)有11人回復(fù)
求助一個(gè)矩陣求解問(wèn)題已經(jīng)有11人回復(fù)
關(guān)于多體系統(tǒng)傳遞矩陣法的問(wèn)題已經(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ù)

1樓 2010-05-12 19:51:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bluesine

鐵桿木蟲(chóng) (職業(yè)作家)

科苑小木蟲(chóng)

應(yīng)助: 132 (高中生)
貴賓: 1.991
金幣: 9852.4
散金: 89
紅花: 19
帖子: 3623
在線(xiàn): 377.5小時(shí)
蟲(chóng)號(hào): 869544
注冊(cè): 2009-10-12
性別: GG
專(zhuān)業(yè): 數(shù)學(xué)物理

nono2009:下次讓樓主重新發(fā)帖比較好。轉(zhuǎn)貼容易引起版面管理上的問(wèn)題。謝謝。 2010-05-14 13:02:17
chunshui2003(金幣+1):謝謝bluesine的建議，剛來(lái)本論壇不久，希望以后能和大家交流，我的專(zhuān)業(yè)水輪發(fā)電機(jī)組振動(dòng)，也涉及到大量的非線(xiàn)性方面內(nèi)容。 2010-05-15 12:26:24

編程的話(huà)請(qǐng)轉(zhuǎn)到程序板塊。

贊一下(1人)

回復(fù)此樓

板凳要做十年冷文章不發(fā)一個(gè)字

4樓2010-05-13 17:08:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 6 個(gè)回答

luomingqi

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 47 (小學(xué)生)
金幣: 4356.8
散金: 20
紅花: 1
帖子: 870
在線(xiàn): 66.5小時(shí)
蟲(chóng)號(hào): 896596
注冊(cè): 2009-11-07
性別: GG
專(zhuān)業(yè): 控制理論與方法

★
小雨萌萌(金幣+1):辛苦獎(jiǎng)，鼓勵(lì)一下 2010-05-12 22:22:56

不是很理解的，我搞的是線(xiàn)性系統(tǒng)的，常用MATLAB的，但是你的這個(gè)東西，我怎么的都沒(méi)有看懂的啊

贊一下(1人)

回復(fù)此樓

跟蹤

2樓2010-05-12 22:13:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

chunshui2003

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 16312.9
散金: 1992
紅花: 14
帖子: 1397
在線(xiàn): 168小時(shí)
蟲(chóng)號(hào): 1000957
注冊(cè): 2010-04-19
專(zhuān)業(yè): 結(jié)構(gòu)工程

這是求解一個(gè)軸系統(tǒng)臨界轉(zhuǎn)速的程序，用的是傳遞矩陣法，Riccati的傳遞矩陣法。如果你要是搞轉(zhuǎn)子應(yīng)該知道。上面的程序是我在別人那里獲得的，但他說(shuō)有問(wèn)題。我覺(jué)得也是，主要是建立“S”之后，我不太清楚如果求轉(zhuǎn)速下一步的程序如何。我表達(dá)的可能不清楚，如果你有《高等轉(zhuǎn)子動(dòng)力學(xué)》看看就明白了。

贊一下(1人)

回復(fù)此樓

3樓2010-05-13 12:11:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jackie_qqq

銀蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 317
帖子: 19
在線(xiàn): 8.9小時(shí)
蟲(chóng)號(hào): 1788854
注冊(cè): 2012-05-01
專(zhuān)業(yè): 金屬材料的凝固與結(jié)晶學(xué)

★
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖

我現(xiàn)在做的是雙轉(zhuǎn)子的子結(jié)構(gòu)傳遞矩陣法！參考下！謝謝啦

贊一下

回復(fù)此樓

5樓2012-05-02 13:48:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 6 個(gè)回答

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿北化085600材料專(zhuān)碩275\|有文章專(zhuān)利｜求調(diào)劑 +11	Micky11223 2026-03-25	11/550	2026-03-31 04:05 by fmesaito
[考研] 307求調(diào)劑 +18	超級(jí)伊昂大王 2026-03-24	20/1000	2026-03-31 01:16 by BruceLiu320
[考研] 085600 286分材料求調(diào)劑 +11	麻辣魷魚(yú) 2026-03-27	12/600	2026-03-30 19:33 by Wang200018
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +7	nhj_ 2026-03-25	7/350	2026-03-30 12:44 by fangnagu
[考研] 求化學(xué)調(diào)劑 +11	wulanna 2026-03-28	11/550	2026-03-30 10:59 by 探123
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +7	傾____蕭 2026-03-27	8/400	2026-03-30 09:37 by longlotian
[考研] 290求調(diào)劑 +3	dfffsar 2026-03-29	3/150	2026-03-29 22:38 by 毛毛毛阿莫2
[考研] 086000生物與醫(yī)藥調(diào)劑 +5	Feisty。 2026-03-28	9/450	2026-03-29 12:02 by longlotian
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8	邱gl 2026-03-27	8/400	2026-03-28 12:42 by 唐沐兒
[考研] 0856，材料與化工321分求調(diào)劑 +12	大饞小子 2026-03-27	13/650	2026-03-28 10:56 by self2008
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3	@taotao 2026-03-28	3/150	2026-03-28 10:26 by JourneyLucky
[考研] 265求調(diào)劑11408 +3	劉小鹿lu 2026-03-27	3/150	2026-03-27 20:53 by nihaoar
[考研] 305求調(diào)劑 +5	哇盧卡庫(kù) 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 085601 材料工程 313分求調(diào)劑 +5	Ong3 2026-03-27	5/250	2026-03-27 12:24 by goldfish51
[考研] 0703化學(xué)338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 336材料求調(diào)劑 +7	陳瀅瑩 2026-03-26	9/450	2026-03-27 00:20 by wxiongid
[考研] 321求調(diào)劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚(yú)的貓
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4	w320357296 2026-03-25	6/300	2026-03-25 17:59 by 255671
[考研] 【2026考研調(diào)劑】制藥工程 284分求相關(guān)專(zhuān)業(yè)調(diào)劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt
[考研] 材料考研調(diào)劑生 +3	黃粱一夢(mèng)千年 2026-03-24	3/150	2026-03-24 17:00 by barlinike