版塊導(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í)熱門版塊排行榜

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

返回列表

hustzcw

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 33.5
帖子: 6
在線: 2.5小時(shí)
蟲號(hào): 2519149
注冊(cè): 2013-06-24
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

[求助] 直接傳遞矩陣法求雙轉(zhuǎn)子系統(tǒng)臨界轉(zhuǎn)速報(bào)錯(cuò)Conversion to logical from sym

CODE:

rou=7754;

E=2.06e11;

%L1=1.20;內(nèi)轉(zhuǎn)軸的總長(zhǎng)

%L2=0.448;外轉(zhuǎn)軸的總長(zhǎng)

K1=1.236e6; K2=2.472e6;K3=4.025e6;K4=1.342e6;%彈性支撐剛度

L1=[0.11 0.04 0.11 0.11 0.11 0.11 0.11 0.130 0.068 0.106 0.146 0.05]; %內(nèi)轉(zhuǎn)軸的分段

L2=[0.09 0.05 0.11 0.06 0.06 0.06 0.018];%外轉(zhuǎn)軸的分段

r1=0.01;r2=0.018;R2=0.0225;%內(nèi)外軸半徑

rp1=0.09;h1=0.03;rp2=0.115;h2=0.018;%兩輪盤半徑及厚度

J1=pi*(r1^4)/4;J2=pi*(R2^4-r2^4)/4;%內(nèi)外軸截面慣性矩

Ip1=rou*pi*rp1^4/2;Id1=Ip1/2;Ip2=rou*pi*rp2^4/2;Id2=Ip2/2;%內(nèi)外軸輪盤極轉(zhuǎn)動(dòng)慣量及直徑轉(zhuǎn)動(dòng)慣量

M1=zeros(1,length(L1)+1);M2=zeros(1,length(L2)+1);

for jk1=1:length(L1)

m1(jk1)=rou*pi*r1^2*L1(jk1)/2 ;%內(nèi)轉(zhuǎn)軸的質(zhì)量分布

end

for jk2=1:length(L2)

m2(jk2)=rou*pi*(R2^2-r1^2)*L2(jk2)/2 ;%外轉(zhuǎn)軸的質(zhì)量分布

end

M1(:,2:length(L1))=m1(:,1:length(L1)-1)+m1(:,2:length(L1));

M1(1)=m1(1);

M1(length(L1)+1)=m1(length(L1));

M2(:,2:length(L2))=m2(:,1:length(L2)-1)+m2(:,2:length(L2));

M2(1)=m2(1);

M2(length(L2)+1)=m2(length(L2));

M1(10)=M1(10)+rou*pi*(rp1^2-r1^2)*h1;%低壓軸輪盤的等效質(zhì)量

M2(4)=M1(4)+rou*pi*(rp2^2-R2^2)*h2;%高壓軸輪盤的等效質(zhì)量

k=0;

 Tit=['第一階頻率的振型和彎矩圖';'第二階頻率的振型和彎矩圖';'第三階頻率的振型和彎矩圖'];

 syms w;%for w=0:0.01:2000;

for i1=1:length(L1)

T1(:,:,i1)=[1 L1(i1) L1(i1)^2/(2*E*J1) L1(i1)^3/(6*E*J1);0 1 L1(i1)/(E*J1) L1(i1)^2/(2*E*J1);0 0 1 L1(i1);0 0 0 1]*[1 0 0 0;0 1 0 0;0 0 1 0;M1(i1)*w^2 0 0 1];

T2(:,:,i1)=[1 0 0 0;0 1 0 0;0 0 1 0;M1(i1+1)*w^2 0 0 1]*[1 L1(i1) L1(i1)^2/(2*E*J1) L1(i1)^3/(6*E*J1);0 1 L1(i1)/(E*J1) L1(i1)^2/(2*E*J1);0 0 1 L1(i1);0 0 0 1];

T3=[1 L1(3) L1(3)^2/(2*E*J1) L1(3)^3/(6*E*J1);0 1 L1(3)/(E*J1) L1(3)^2/(2*E*J1);0 0 1 L1(3);0 0 0 1]*[1 0 0 0;0 1 0 0;0 0 1 0;(M1(3)*w^2-K1) 0 0 1];

T4=[1 0 0 0;0 1 0 0;0 0 1 0;(M1(12)*w^2-K2) 0 0 1]*[1 L1(11) L1(11)^2/(2*E*J1) L1(11)^3/(6*E*J1);0 1 L1(11)/(E*J1) L1(11)^2/(2*E*J1);0 0 1 L1(11);0 0 0 1];

T5=[1 0 0 0;0 1 0 0;0 (Ip1*w^2/1.2-Id1*w^2) 1 0;(M1(10)*w^2-K4) 0 0 1];%Hn1R(n1L)

end

H1=T1(:,:,4);

for i2=5:9;

H1=T1(:,:,i2)*H1;

end

H10L=H1*T3*T1(:,:,2)*T1(:,:,1);

H10R=T2(:,:,12)*T4*T2(:,:,10);

Hd1=H10R*T5*H10L;%Hn2R(1L)

Hd2=H10R*T5;%Hn2R(n1L)

for i3=1:length(L2)

T6(:,:,i3)=[1 L2(i3) L2(i3)^2/(2*E*J2) L2(i3)^3/(6*E*J2);0 1 L2(i3)/(E*J2) L2(i3)^2/(2*E*J2);0 0 1 L2(i3);0 0 0 1]*[1 0 0 0;0 1 0 0;0 0 1 0;M2(i3)*w^2 0 0 1];

T7=[1 L2(2) L2(2)^2/(2*E*J2) L2(2)^3/(6*E*J2);0 1 L2(2)/(E*J2) L2(2)^2/(2*E*J2);0 0 1 L2(2);0 0 0 1]*[1 0 0 0;0 1 0 0;0 0 1 0;(M2(2)*w^2-K3) 0 0 1];

T8=[1 L2(4) L2(4)^2/(2*E*J2) L2(4)^3/(6*E*J2);0 1 L2(4)/(E*J2) L2(4)^2/(2*E*J2);0 0 1 L2(4);0 0 0 1]*[1 0 0 0;0 1 0 0;0 (Ip2*w^2-Id2*w^2) 1 0;M2(4)*w^2 0 0 1];

T9(:,:,i3)=[1 0 0 0;0 1 0 0;0 0 1 0;M2(i3+1)*w^2 0 0 1]*[1 L2(i3) L2(i3)^2/(2*E*J2) L2(i3)^3/(6*E*J2);0 1 L2(i3)/(E*J2) L2(i3)^2/(2*E*J2);0 0 1 L2(i3);0 0 0 1];

T10=[1 0 0 0;0 1 0 0;0 0 1 0;(M2(7)*w^2-K4) 0 0 1];

end

H20L=T6(:,:,6)*T6(:,:,5)*T8*T6(:,:,3)*T7*T6(:,:,1);

syms x;syms y;

Zn3L=[x,y,0,0]';

Zn4L=H20L*Zn3L;

syms x10;

Zn5R=T9(:,:,7)*T10*(Zn4L+[0 0 0 K4*x10]');

%%

X1=Zn5R(3,:);

X2=Zn5R(4,:);

S=solve(X1,X2,x,y);

U=S.x/x10;

V=S.y/x10;

a=H20L(1,1)*U+H20L(1,2)*V;

G1=Hd1(3,1)+a*Hd2(3,4)*K4*H10L(1,1);

G2=Hd1(3,2)+a*Hd2(3,4)*K4*H10L(1,2);

G3=Hd1(4,1)+a*Hd2(4,4)*K4*H10L(1,1);

G4=Hd1(4,2)+a*Hd2(4,4)*K4*H10L(1,2);

F=G1*G4-G2*G3 ; %剩余量

for w=0:0.01:2000;

if F*(-1)^k<0 %求解臨界轉(zhuǎn)速

k=k+1;

        wi(k)=w ;    %固有圓頻率

        w=wi(k);

         ni(k)=wi(k)*30/pi ; %臨界轉(zhuǎn)速

end

end

 ni=ni'  

 wi=wi' 

f=ni/60

參考南航胡絢的博士論文《反向旋轉(zhuǎn)雙轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)特性研究》第11頁(yè)的直接傳遞矩陣法編寫的計(jì)算臨界轉(zhuǎn)速的程序。報(bào)錯(cuò)Conversion to logical from sym is not possible. 該如何修改？

[ Last edited by jjdg on 2013-8-28 at 10:07 ]

回復(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材料專碩275|有文章專利｜求調(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ì)您同樣有幫助:

【求助】Riccati傳遞矩陣法求解臨界轉(zhuǎn)速（有部分程序）已經(jīng)有5人回復(fù)

1樓 2013-08-27 22:11:24

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

wkl1990

鐵桿木蟲 (正式寫手)

應(yīng)助: 17 (小學(xué)生)
金幣: 6473.4
散金: 618
紅花: 20
帖子: 505
在線: 266小時(shí)
蟲號(hào): 3203530
注冊(cè): 2014-05-13
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)

樓主你好，極轉(zhuǎn)動(dòng)慣量和直徑轉(zhuǎn)動(dòng)慣量是什么區(qū)別，我看到你上面有公式，但是我想找到它們的物理意義并理解一下，能告訴我這些公式的出處嗎？謝謝！

贊一下

回復(fù)此樓

2樓2014-07-01 16:36:48

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 hustzcw 的主題更新