luckylin 发表于 2012-3-26 08:33

哪位热心人士能帮我分析下我用matlab所建立转子模型

为什么画出来的转子后半部分悬臂轴显示不出来？
clear
format short e
close all

set(0,'defaultaxesfontsize',24)
set(0,'defaultaxesfontname','Times New Roman')
set(0,'defaulttextfontsize',24)
set(0,'defaulttextfontname','Times New Roman')

E = 2.11e11;
Poisson = 0.3;
G = E/(2*(1+Poisson));
rho = 7810;
damping_factor = 0;
Length = 630;

Rotor.node = [1 0;2 17;3 34;4 51;5 68;
    6 85;   % bearing 1
    7 107.2;8 129.4;9 151.6;10 173.8;
    11 196;   %left side of disc1
    12 204;   %middle of disc1
    13 212;   %right side of disc1
    14 243;15 274;16 305;
    17 336;   
    18 367;   
    19 398;   %left side of disc2
    20 406;   %middle of disc1
    21 414;   %right side of disc2
    22 446;23 478;24 510;
    25 542;   % bearing 2
    26 564;27 586;28 608;29 630];
% 各种节点定义：
Node_bearing = ;
Node_disc = ;
Rotor.shaft = [ 212 13 0 rho E G damping_factor;
    223 13 0 rho E G damping_factor;
    234 13 0 rho E G damping_factor;
    245 13 0 rho E G damping_factor;
    256 13 0 rho E G damping_factor;
    267 13 0 rho E G damping_factor; %bearing 1
    278 13 0 rho E G damping_factor;
    289 13 0 rho E G damping_factor;
    2910 13 0 rho E G damping_factor;
    21011 13 0 rho E G damping_factor;
    21112 13 0 rho E G damping_factor; %disc 1
    21213 13 0 rho E G damping_factor; %disc 1
    21314 13 0 rho E G damping_factor;
    21415 13 0 rho E G damping_factor;
    215 16 13 0 rho E G damping_factor;
    216 17 13 0 rho E G damping_factor;
    217 18 13 0 rho E G damping_factor;
    218 19 13 0 rho E G damping_factor;
    21920 13 0 rho E G damping_factor; %disc 2
    22021 13 0 rho E G damping_factor; %disc 2
    22122 13 0 rho E G damping_factor;
    22223 13 0 rho E G damping_factor;
    223 24 13 0 rho E G damping_factor;
    224 25 13 0 rho E G damping_factor;
    225 26 0.04 0 rho E G damping_factor;%bearing 2
    226 27 0.04 0 rho E G damping_factor;
    2 27 28 0.04 0 rho E G damping_factor;
    2 28 29 0.04 0 rho E G damping_factor;
    ];
Rotor.disc = [
                1 Node_disc(1) rho 16 152 13; %%%%% Disc 1
                1 Node_disc(2) rho 16 152 13; %%%%% Disc 2
                ];      
% Define bearing characteristics - and put bearings at nodes 6 and 25
Bear_Char_1 = ;
Bear_Char_2 = ;
Bear_Char_3 = ;
Bear_Char_4 = ;

Bear_Char_5 = ;
Bear_Char_6 = ;
Bear_Char_7 = ;
Bear_Char_8 = ;

Bear_Char_9 = ;
Bear_Char_10 = ;

Bear_Char_11 = ;

Rotor.bearing = [3 Node_bearing(1) Bear_Char_9; ...
                     3 Node_bearing(2) Bear_Char_9];
% draw picture of the rotor         
figure(1)
clf
picrotor(Rotor)

ChaChing 发表于 2012-3-27 01:09

picrotor ???
??? Undefined function or method 'picrotor' for input arguments of type 'struct'.
页: [1]
查看完整版本: 哪位热心人士能帮我分析下我用matlab所建立转子模型