“Riccati传递矩阵法计算汽轮机转子模型临界转速_不平衡响应_转子传递矩阵_转子临界转速计算_转子_

function output=calculate3(in) syms p k; p=in; kz=2.0e6; %轴承支撑节点刚度 %定义转轴参数，共分14段 H1=H(0.15,1.02e-7,8.85,p); % in1=l in2=J in3=m in4=p H2=H(0.25,1.84e-7,11.86,p); H3=H(0.25,2.61e-7,14.11,p); H4=H(0.25,3.59e-7,16.56,p); H5=H(0.25,2.20e-7,12.96,p); H6=H(0.35,6.56e-7,7.08,p); H7=H(0.13,6.56e-7,7.08,p); H8=H(0.22,1.84e-7,11.86,p); H9=H(0.22,2.61e-7,14.11,p); H10=H(0.22,3.59e-7,16.56,p); H11=H(0.22,2.61e-7,14.11,p); H12=H(0.22,1.84e-7,11.86,p); %转盘 Ip1=0.186;Id1=0.093;m1=16.15; Ip2=0.345;Id2=0.172;m2=22.05; Ip3=0.588;Id3=0.294;m3=28.98; Ip4=0.06;Id4=0.03;m4=5; Ip5=0.588;Id5=0.294;m5=29.02; Ip6=0.186;Id6=0.093;m6=16.08; Ip7=0.186;Id7=0.093;m7=16.08; Ip8=0.588;Id8=0.294;m8=29.02; %求解过程 S1=[0 0;0 0]; T22=[H1(1,1:2);H1(2,1:2)]; T21=[H1(1,3:4);H1(2,3:4)]; T12=[H1(3,1:2);H1(4,1:2)]; T11=[H1(3,3:4);H1(4,3:4)]; S2=(T11*S1+T12)*inv(T21*S1+T22); iS2=inv(S2); KK=[0 0;-kz 0]; I=eye(2); iS3=iS2*inv(I+KK*iS2); T11=[H2(1,1:2);H2(2,1:2)]; T12=[H2(1,3:4);H2(2,3:4)]; T21=[H2(3,1:2);H2(4,1:2)]; T22=[H2(3,3:4);H2(4,3:4)]; S4=(T21*iS3+T22)*inv(T11*iS3+T12); K1=[1 0 0 0;0 1 0 0;0 -(1-Ip1/Id1)*Id1*p^2 1 0;m1*p^2 0 0 1]; T22=[K1(1,1:2);K1(2,1:2)]; T21=[K1(1,3:4);K1(2,3:4)]; T12=[K1(3,1:2);K1(4,1:2)]; T11=[K1(3,3:4);K1(4,3:4)]; S5=(T11*S4+T12)*inv(T21*S4+T22); T22=[H3(1,1:2);H3(2,1:2)]; T21=[H3(1,3:4);H3(2,3:4)]; T12=[H3(3,1:2);H3(4,1:2)]; T11=[H3(3,3:4);H3(4,3:4)]; S6=(T11*S5+T12)*inv(T21*S5+T22); K2=[1 0 0 0;0 1 0 0;0 -(1-Ip2/Id2)*Id2*p^2 1 0;m2*p^2 0 0 1]; T22=[K2(1,1:2);K2(2,1:2)]; T21=[K2(1,3:4);K2(2,3:4)]; T12=[K2(3,1:2);K2(4,1:2)]; T11=[K2(3,3:4);K2(4,3:4)]; S7=(T11*S6+T12)*inv(T21*S6+T22); T22=[H4(1,1:2);H4(2,1:2)]; T21=[H4(1,3:4);H4(2,3:4)]; T12=[H4(3,1:2);H4(4,1:2)]; T11=[H4(3,3:4);H4(4,3:4)]; S8=(T11*S7+T12)*inv(T21*S7+T22); K3=[1 0 0 0;0 1 0 0;0 -(1-Ip3/Id3)*Id3*p^2 1 0;m3*p^2 0 0 1]; T22=[K3(1,1:2);K3(2,1:2)]; T21=[K3(1,3:4);K3(2,3:4)]; T12=[K3(3,1:2);K3(4,1:2)]; T11=[K3(3,3:4);K3(4,3:4)]; S9=(T11*S8+T12)*inv(T21*S8+T22); T22=[H5(1,1:2);H5(2,1:2)]; T21=[H5(1,3:4);H5(2,3:4)]; T12=[H5(3,1:2);H5(4,1:2)]; T11=[H5(3,3:4);H5(4,3:4)]; S10=(T11*S9+T12)*inv(T21*S9+T22); iS10=inv(S10); iS11=iS10*inv(I+KK*iS10); T11=[H6(1,1:2);H6(2,1:2)]; T12=[H6(1,3:4);H6(2,3:4)]; T21=[H6(3,1:2);H6(4,1:2)]; T22=[H6(3,3:4);H6(4,3:4)]; S12=(T21*iS11+T22)*inv(T11*iS11+T12); K4=[1 0 0 0;0 1 0 0;0 -(1-Ip4/Id4)*Id4*p^2 1 0;m4*p^2 0 0 1]; T22=[K4(1,1:2);K4(2,1:2)]; T21=[K4(1,3:4);K4(2,3:4)]; T12=[K4(3,1:2);K4(4,1:2)]; T11=[K4(3,3:4);K4(4,3:4)]; S13=(T11*S12+T12)*inv(T21*S12+T22); T22=[H7(1,1:2);H7(2,1:2)]; T21=[H7(1,3:4);H7(2,3:4)]; T12=[H7(3,1:2);H7(4,1:2)]; T11=[H7(3,3:4);H7(4,3:4)]; S14=(T11*S13+T12)*inv(T21*S13+T22); iS14=inv(S14); iS15=iS14*inv(I+KK*iS14); T11=[H8(1,1:2);H8(2,1:2)]; T12=[H8(1,3:4);H8(2,3:4)]; T21=[H8(3,1:2);H8(4,1:2)]; T22=[H8(3,3:4);H8(4,3:4)]; S16=(T21*iS15+T22)*inv(T11*iS15+T12); K5=[1 0 0 0;0 1 0 0;0 -(1-Ip5/Id5)*Id5*p^2 1 0;m5*p^2 0 0 1]; T22=[K5(1,1:2);K5(2,1:2)]; T21=[K5(1,3:4);K5(2,3:4)]; T12=[K5(3,1:2);K5(4,1:2)]; T11=[K5(3,3:4);K5(4,3:4)]; S17=(T11*S16+T12)*inv(T21*S16+T22); T22=[H9(1,1:2);H9(2,1:2)]; T21=[H9(1,3:4);H9(2,3:4)]; T12=[H9(3,1:2);H9(4,1:2)]; T11=[H9(3,3:4);H9(4,3:4)]; S18=(T11*S17+T12)*inv(T21*S17+T22); K6=[1 0 0 0;0 1 0 0;0 -(1-Ip6/Id6)*Id6*p^2 1 0;m6*p^2 0 0 1]; T22=[K6(1,1:2);K6(2,1:2)]; T21=[K6(1,3:4);K6(2,3:4)]; T12=[K6(3,1:2);K6(4,1:2)]; T11=[K6(3,3:4);K6(4,3:4)]; S19=(T11*S18+T12)*inv(T21*S18+T22); T22=[H10(1,1:2);H10(2,1:2)]; T21=[H10(1,3:4);H10(2,3:4)]; T12=[H10(3,1:2);H10(4,1:2)]; T11=[H10(3,3:4);H10(4,3:4)]; S20=(T11*S19+T12)*inv(T21*S19+T22); K7=[1 0 0 0;0 1 0 0;0 -(1-Ip7/Id7)*Id7*p^2 1 0;m7*p^2 0 0 1]; T22=[K7(1,1:2);K7(2,1:2)]; T21=[K7(1,3:4);K7(2,3:4)]; T12=[K7(3,1:2);K7(4,1:2)]; T11=[K7(3,3:4);K7(4,3:4)]; S21=(T11*S20+T12)*inv(T21*S20+T22); T22=[H11(1,1:2);H11(2,1:2)]; T21=[H11(1,3:4);H11(2,3:4)]; T12=[H11(3,1:2);H11(4,1:2)]; T11=[H11(3,3:4);H11(4,3:4)]; S22=(T11*S21+T12)*inv(T21*S21+T22); K8=[1 0 0 0;0 1 0 0;0 -(1-Ip8/Id8)*Id8*p^2 1 0;m8*p^2 0 0 1]; T22=[K8(1,1:2);K8(2,1:2)]; T21=[K8(1,3:4);K8(2,3:4)]; T12=[K8(3,1:2);K8(4,1:2)]; T11=[K8(3,3:4);K8(4,3:4)]; S23=(T11*S22+T12)*inv(T21*S22+T22); T22=[H12(1,1:2);H12(2,1:2)]; T21=[H12(1,3:4);H12(2,3:4)]; T12=[H12(3,1:2);H12(4,1:2)]; T11=[H12(3,3:4);H12(4,3:4)]; S24=(T11*S23+T12)*inv(T21*S23+T22); % iS24=inv(S24); % iS25=iS24*inv(I+KK*iS24); ee=det(24); output=double(ee);