alxhslm/ROLADYN
Fork: 6 Star: 36 (更新于 2024-11-16 03:20:05)
license: GPL-3.0
Language: MATLAB .
MATLAB toolbox for simulating the lateral dynamics of rotating machines.
ROLADYN
MATLAB suite for simulating the lateral dynamics of rotating machines.
Example rotor-bearing system
Define a rotating machine with a flexible rotor as follows:
P.Rotor = {};
P.Bearing = {};
P.Stator = {};
% Rotor
P.Rotor{1}.Name = 'Main';
P.Rotor{1}.Speed = 1;
P.Rotor{1}.Nodes = [0, 0.25, 0.5, 0.75, 1.0]; % Nodes of the FEM
%% Flexible shaft
P.Rotor{1}.Shaft = {};
P.Rotor{1}.Shaft{1}.Name = 'Shaft';
P.Rotor{1}.Shaft{1}.iNodes = 1:5;
P.Rotor{1}.Shaft{1}.Section.ro = 0.0;
P.Rotor{1}.Shaft{1}.Section.ri = 0.025;
P.Rotor{1}.Shaft{1}.Material.Name = 'steel';
P.Rotor{1}.Shaft{1}.Material.E = 210E9;
P.Rotor{1}.Shaft{1}.Material.rho = 7800;
%% Disc
P.Rotor{1}.Disc = {};
P.Rotor{1}.Disc{1}.Name = 'EndCap';
P.Rotor{1}.Disc{1}.iNode = 3;
P.Rotor{1}.Disc{1}.Type = 'Rigid';
P.Rotor{1}.Disc{1}.Options.bGyro = 1;
P.Rotor{1}.Disc{1}.Material.Name = 'steel';
P.Rotor{1}.Disc{1}.Material.E = Inf; %Rigid disc
P.Rotor{1}.Disc{1}.Material.rho = 7800;
P.Rotor{1}.Disc{1}.Ring.Geometry.R = [0 0.3];
P.Rotor{1}.Disc{1}.Ring.Geometry.t = 15E-3;
% Bearings
P.Bearing{1}.Name = 'LeftBearing';
P.Bearing{1}.kyy = 1E3;
P.Bearing{1}.kxx = 1E3;
P.Bearing{1}.cxx = 1E-3;
P.Bearing{1}.cyy = 1E-3;
P.Bearing{1}.Node{1}.Type = 'rotor';
P.Bearing{1}.Node{1}.iNode = 1;
P.Bearing{1}.Node{1}.iRotor = 1;
P.Bearing{1}.z = 0.0;
P.Bearing{2}.Name = 'RightBearing';
P.Bearing{2}.kyy = 1E3;
P.Bearing{2}.kxx = 1E3;
P.Bearing{2}.cxx = 1E-3;
P.Bearing{2}.cyy = 1E-3;
P.Bearing{2}.Node{1}.Type = 'rotor';
P.Bearing{2}.Node{1}.iNode = 5;
P.Bearing{2}.Node{1}.iRotor = 1;
P.Bearing{2}.z = 1.0;
% Excitation
P.Excite{1}.Name = 'Unbalance';
P.Excite{1}.Type = 'unbalance';
P.Excite{1}.Mode = 'sync';
P.Excite{1}.iRotor = 1;
P.Excite{1}.iDisc = 1;
P.Excite{1}.m = 0.01;
P.Excite{1}.r = 0.1;
P.Excite{1}.Angle = 0.0;
P = setupsystem(P);
You can visualise the system with:
plot_system(P);
which will generate a plot of the rototing machine like this:
Modal analysis
To compute the natural frequencies at a range of speeds:
O = linspace(0,50,100)*2*pi; % Speed in rad/s
[V,d,W] = rotor_eig(P.Model,O);
where:
-
V
andW
contain the left and right eigevectors -
d
contains the eigenvalues
at each speed.
To convert these eigenvalues and eigenvectors vectors into natural frequencies and mode shapes, run the following:
iNode = 2;
[omega,zeta, modes] = eig2modes(V,d,W);
% sort the modes
[~,ii] = sort(omega(:,1));
omega = omega(ii,:);
zeta = zeta(ii,:);
modes = modes(:,ii,:);
kappa = compute_whirl(P,modes,iNode);
where:
-
omega
contains the natural frequencies in rad/s -
zeta
contains the non-dimensionalised damping ratios -
kappa
contains the whirl direction, defined in terms of the motion at the central disc (node 2).
Campbell diagram
To plot the Campbell diagram with the first 2 modes, run the following
options.xunits = 'Hz';
iPlot = 1:2;
plot_campbell(P,O,omega(iPlot,:),zeta(iPlot,:),kappa(:,iPlot,:),iNode, options);
which will generate a plot of the natural frequencies as a function of rotation speed like this:
Whirl
To plot the shape of the response at each resonance (ie each mode shape), run the following:
plot_whirl(P,O,modes(:,iPlot,:),omega(iPlot,:),kappa(:,iPlot,:));
There is a slider which you can use to change which rotation speed the modes are plotted from.
Synchronous response
Compute the response at all frequencies with:
[U,Udot,Uddot] = excitation_frf(P,O);
F = P.Model.Excite.K*U + P.Model.Excite.C*Udot + P.Model.Excite.M*Uddot;
q = rotor_frf_modal(V,d,W,F,O);
And then plot the FRF of the disc response (DOF 9 and 10) using:
options.xunits = 'Hz';
options.iPlot = 9:10;
plot_frf(P,q,O,options);
最近版本更新:(数据更新于 2024-09-10 10:11:50)
主题(topics):
bearing, model, rotor, rotor-dynamic, structural-dynamics, vibration
alxhslm/ROLADYN同语言 MATLAB最近更新仓库
2024-11-05 14:39:56 BIMK/PlatEMO
2023-11-01 01:41:52 cvnlab/GLMsingle
2022-09-02 00:24:19 kmansouri/OPERA
2020-11-30 22:26:07 ultralytics/functions-matlab
1970-01-01 00:00:00 KeshengZhang/PlatEMOv2.6
1970-01-01 00:00:00 tomtkg/Test_Functions_for_Multi-objective_Optimization