计算转子——轴承系统不平衡响应

这是用于计算转子——轴承系统不平衡响应的程序，转子用集总参数法简化，轴承支撑为各项同性支撑

% ================= 初始条件 ===================
function [NKm_kJp_kJd_klEIrr]=initial_parameters
N=35;              % 轴段数
E=2.06e11;         % 弹性模量
u=0.3;             % 泊松比
a=0.886;               % 截面影响系数
G=E/（2*（1+u））;
titlestr=‘转子参数‘;
pathstr=‘E:\传递矩阵法\高中压转子计算‘;
m_k=xlsread（[pathstr‘\高中压转子初始参数表.xls‘]titlestr‘B2:B37‘）;   % 从exel中导入转子的参数
Jp_k=xlsread（[pathstr‘\高中压转子初始参数表.xls‘]titlestr‘C2:C37‘）;
Jd_k=xlsread（[pathstr‘\高中压转子初始参数表.xls‘]titlestr‘D2:D37‘）;
D=xlsread（[pathstr‘\高中压转子初始参数表.xls‘]titlestr‘E2:E36‘）;
EI=xlsread（[pathstr‘\高中压转子初始参数表.xls‘]titlestr‘F2:F36‘）;
l=xlsread（[pathstr‘\高中压转子初始参数表.xls‘]titlestr‘G2:G36‘）;
d=xlsread（[pathstr‘\高中压转子初始参数表.xls‘]titlestr‘H2:H36‘）;
K=xlsread（[pathstr‘\高中压转子初始参数表.xls‘]titlestr‘I2:I37‘）;   
for i=1:1:N
    rr（i）=6*EI（i）/（a*G*（pi*（D（i）^2-d（i）^2）/4）*l（i）^2）;   % 剪切影响系数
end
% ============================================

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

     文件       1044  2011-05-13 10:33  riccati传递矩阵法计算转子不平衡响应\initial_parameters.m

     文件       4859  2011-05-11 18:59  riccati传递矩阵法计算转子不平衡响应\transfermatrix_calculate.m

     文件       2136  2011-06-29 14:55  riccati传递矩阵法计算转子不平衡响应\unbalance.m

     目录          0  2011-08-26 17:16  riccati传递矩阵法计算转子不平衡响应

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

                 8039                    4