转子临界转速-Matlab

转子系统的临界转速计算，Matlab程序

clc
clear all;
Nshaft=8;                            %轴段数量；
%%
RotorE = 2.095e11;                                                                        %转子弹性模量；
RotorM = 7.85e3;                                                                          %转子材料密度
ShaftL = [0.250.250.250.250.250.250.250.25];                                       %各轴段长度
ShaftDI = ones（1Nshaft）*0.1;                                                             %各轴段外径
ShaftDO = ones（1Nshaft）*0.0;                                                             %各轴段内径；
SPtype= 1;                                                                                %支承类型选择， 1 -- 刚性支承， 2 -- 弹性支承；
LocationF=[19];                                                                          %支承所在节点编号；
SPstiffness=[];                                                                           %支承刚度，SPtype = 1时，无需进行任何设置；
%%
addtionN = [];                                 %附加轮盘编号
addtionM = [];                                 %附加轮盘质量
addtionJ = [];                                 %附加轮盘转动惯量
%%
CSN = 5;                                     %输出的临界转速数量
CriticalSpeeds=Chinavib_CriticalSpeeds（NshaftRotorERotorMShaftLShaftDIShaftDOSPtypeLocationFSPstiffnessaddtionNaddtionMaddtionJCSN）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2198  2014-10-09 10:35  Critical Speed\Chinavib_CriticalSpeeds.p

     文件       1386  2014-10-09 10:35  Critical Speed\shaft.m

     目录          0  2014-10-09 10:37  Critical Speed

----------- ---------  ---------- -----  ----

                 3584                    3