Scribd
Importer
34 vues2 pages

Partie 2 Script:: ALL All

Le script MATLAB calcule les fréquences et modes propres d'un rotor en utilisant deux cas de raideur de palier. Dans le premier cas, les raideurs des paliers sont égales, tandis que dans le second, elles diffèrent. Les résultats sont affichés pour chaque cas, incluant les fréquences propres en Hz et les modes propres associés.

Transféré par

zakiben305
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats RTF, PDF, TXT ou lisez en ligne sur Scribd
34 vues2 pages

Partie 2 Script:: ALL All

Le script MATLAB calcule les fréquences et modes propres d'un rotor en utilisant deux cas de raideur de palier. Dans le premier cas, les raideurs des paliers sont égales, tandis que dans le second, elles diffèrent. Les résultats sont affichés pour chaque cas, incluant les fréquences propres en Hz et les modes propres associés.

Transféré par

zakiben305
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats RTF, PDF, TXT ou lisez en ligne sur Scribd

PARTIE 2 SCRIPT :

clc;
clear ALL;
close all
% Données physiques
L = 0.5; % Longueur du rotor (m)
d = 0.2; % Diamètre du rotor (m)
rho = 7810; % Masse volumique (kg/m^3)
k1 = 1e6; % Raideur palier 1 (N/m)
k2 = 1e6; % Raideur palier 2 (N/m) pour cas 1

% Volume et masse
V = pi * (d/2)^2 * L;
m = rho * V;

% Moment d'inertie autour du centre (basculement)


J = (1/12) * m * L^2;

% Matrices pour système 2 ddl : [ y_G ; theta ]


M = [m, 0;
0, J];

% Cas 1 : k1 = k2 = 1 MN/m
K1 = [k1 + k2, (k2 - k1)*L/2;
(k2 - k1)*L/2, (k1 + k2)*(L^2)/4];

% Résolution du problème aux valeurs propres


[Vecs1, D1] = eig(K1, M);
freqs1 = sqrt(diag(D1)) / (2*pi); % en Hz

disp('--- Cas 1 : k1 = k2 = 1e6 N/m ---');


disp('Fréquences propres (Hz) :');
disp(freqs1);
disp('Modes propres :');
disp(Vecs1);

% Cas 2 : k1 = 1 MN/m, k2 = 1.3 MN/m


k1 = 1e6;
k2 = 1.3e6;

K2 = [k1 + k2, (k2 - k1)*L/2;


(k2 - k1)*L/2, (k1 + k2)*(L^2)/4];

[Vecs2, D2] = eig(K2, M);


freqs2 = sqrt(diag(D2)) / (2*pi); % en Hz

disp('--- Cas 2 : k1 = 1e6, k2 = 1.3e6 N/m ---');


disp('Fréquences propres (Hz) :');
disp(freqs2);
disp('Modes propres :');
disp(Vecs2);

Vous aimerez peut-être aussi