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

luckylin

为什么画出来的转子后半部分悬臂轴显示不出来？
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 = [6 25];
Node_disc = [12 20];
Rotor.shaft = [ 2  1  2 13 0 rho E G damping_factor;
    2  2  3 13 0 rho E G damping_factor;
    2  3  4 13 0 rho E G damping_factor;
    2  4  5 13 0 rho E G damping_factor;
    2  5  6 13 0 rho E G damping_factor;
    2  6  7 13 0 rho E G damping_factor; %bearing 1
    2  7  8 13 0 rho E G damping_factor;
    2  8  9 13 0 rho E G damping_factor;
    2  9  10 13 0 rho E G damping_factor;
    2  10  11 13 0 rho E G damping_factor;
    2  11  12 13 0 rho E G damping_factor; %disc 1
    2  12  13 13 0 rho E G damping_factor; %disc 1
    2  13  14 13 0 rho E G damping_factor;
    2  14  15 13 0 rho E G damping_factor;
    2  15 16 13 0 rho E G damping_factor;
    2  16 17 13 0 rho E G damping_factor;
    2  17 18 13 0 rho E G damping_factor;
    2  18 19 13 0 rho E G damping_factor;  
    2  19  20 13 0 rho E G damping_factor; %disc 2
    2  20  21 13 0 rho E G damping_factor; %disc 2
    2  21  22 13 0 rho E G damping_factor;
    2  22  23 13 0 rho E G damping_factor;
    2  23 24 13 0 rho E G damping_factor;  
    2  24 25 13 0 rho E G damping_factor;
    2  25 26 0.04 0 rho E G damping_factor;  %bearing 2
    2  26 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 = [5e4 5e4 0.0 0.0];
Bear_Char_2 = [5e4 1e5 0.0 0.0];
Bear_Char_3 = [5e4 1e5 2e3 2e3];
Bear_Char_4 = [5e4 1e5 1e4 1e4];

Bear_Char_5 = [10e6 10e6 0.0 0.0];
Bear_Char_6 = [10e6 20e6 0.0 0.0];
Bear_Char_7 = [10e6 20e6 60e3 60e3];
Bear_Char_8 = [10e6 20e6 400e3 400e3];

Bear_Char_9 = [2e5 4e5 0 0];
Bear_Char_10 = [10e6 10e6 30e3 30e3];

Bear_Char_11 = [3.5e5 4e5 0 0];

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

picrotor ???

??? Undefined function or method 'picrotor' for input arguments of type 'struct'.