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