非线性振动 matlab,故障转子系统的非线性振动分析与诊断方法附录A-matlab程序

A.1 传递距阵法分析程序

%main_critical.m

%该程序使用Riccati传递距阵法计算转子系统的临界转速及振型

%本函数中均采用国际单位制

% 第一步：设置初始条件(调用函数shaft_parameters)

%初始值设置包括：轴段数N，搜索次数M

%输入轴段参数：内径d,外径D，轴段长度l，支撑刚度K，单元质量mm，极转动惯量Jpp[N,M,d,D,l,K,mm,Jpp]=shaft_parameters;

% 第二步：计算单元的5个特征值(调用函数shaft_pra_cal)

%单元的5个特征值：

%m_k：:质量

%Jp_k：极转动惯量

%Jd_k：直径转动惯量

%EI：弹性模量与截面对中性轴的惯性矩的乘积

%rr：剪切影响系数

[m_k,Jp_k,EI,rr]=shaft_pra_cal(N,D,d,l,Jpp,mm);

% 第三步：计算剩余量(调用函数surplus_calculate),并绘制剩余量图

%剩余量：D1

for i=1:1:M

ptx(i)=0;

pty(i)=0;

end