Scribd
Importer
27 vues3 pages

Analyse de la Vibration d'un Système

Ce document traite de la modélisation d'un système mécanique soumis à des vibrations. Il présente les équations régissant le mouvement d'un tel système en régime libre ou forcé, amorti ou non, et donne les expressions analytiques des déplacements, vitesses et accélérations en fonction du temps.

Transféré par

Hachetag
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 TXT, PDF, TXT ou lisez en ligne sur Scribd
27 vues3 pages

Analyse de la Vibration d'un Système

Ce document traite de la modélisation d'un système mécanique soumis à des vibrations. Il présente les équations régissant le mouvement d'un tel système en régime libre ou forcé, amorti ou non, et donne les expressions analytiques des déplacements, vitesses et accélérations en fonction du temps.

Transféré par

Hachetag
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 TXT, PDF, TXT ou lisez en ligne sur Scribd

# ---------- TP de Vibration ------------

#Param�tres :
m = 39007e3; #masse
k = 385e6; #raideur
c = 9804; #Amortissement

# 1.2.1 - R�gime libre - Syst�me concervatif


#Question 1 : Determiner la pulsation propre et la fr�quence propre du syst�me

wo = sqrt(k/m); #(3.1417)
fo = wo/(2*pi); #(4.0001e-05)

#Question 2 : Determiner l'expression de la r�ponse en d�placement, vitesse et


acc�l�ration

#m*upp + k*u = 0
#upp+wo�*u = 0

#La solution de u :
#u(t) = Acos(wo*t) + B*sin(wo*t)
#up(t) est la d�riv�e premi�re de u(t) =>
#up(t) = -wo * A * sin(wo*t) + wo*B*cos(wo*t)

#Condition initiale :
#u(t=0) = -10e-3
#up(t=0) = 0

#Ainsi, grace aux conditions initales, on calcul A et B :


A = -10e-3;
B = 0;
t = [0:0.02:5];
u = A*cos(wo*t) + B*sin(wo*t);
up = -wo * A * sin(wo*t) + wo*B*cos(wo*t) ;
upp = (-wo^2 *A * cos(wo*t)) + (-wo^2 * B * sin(wo*t));
figure(1)
plot(t,u,'r')
hold on
plot(t,up,'b')
hold on
plot (t,upp,'g')

#A expliquer dans le rapport

#1.2.2 Regime forc� - Syst�me conservatif

#Question 2 : Determiner l'expression de la r�ponse en d�placement, vitesse et


acc�l�ration

#m*upp + k*u = f
#upp+wo�*u = Fo*sin(Om*t)

#La solution de l'�quation homog�ne :


#uho(t) = Acos(wo*t) + B*sin(wo*t)

#La solution particuli�re :


#uh1(t) = Uh * sin(Om*t)
#u(t) = uho(t) + uh1(t) = (Acos(wo*t) + B*sin(wo*t)) + (Uh * sin(Om*t))
#up(t) = -wo * A * sin(wo*t) + wo*B*cos(wo*t) + (Om*Uh*cos(Om*t))
#upp(t) = (-wo^2*A*cos(wo*t)) +(-wo^2*B*sin(wo*t)) + (-Om^2*Uh*sin(Om*t))

#upp(t) = (-wo^2*A*cos(wo*t)) +(-wo^2*B*sin(wo*t)) + (-Om^2*Uh*sin(Om*t)) = (-


Om^2*F0*sin(Om*t))

#Le d�placement tends vers infinit lorsque Om = wo puisque nous n'avons pas
d'amortissement

# 1.2.3 - R�gime libre - Syst�me dissipatif

#Question 4 : Determiner la pulsation propre et la fr�quence propre du syst�me

wo = sqrt(k/m); #(3.1417)
Ta = c/(2*sqrt(k*m)); #On est sur un mouvement sous-amorti : Ta< 1

#Question 5 : Determiner l'expression de la r�ponse en d�placement

#upp+2Ta*wo*up+wo�*u = 0 avec wo = sqrt(k/m) et ta = c/(2/*sqrt(k*m))

Ta = c/(2*sqrt(k*m)); #On est sur un mouvement sous-amorti : Ta< 1

wd = wo*sqrt(1-Ta^2);

#u = exp(-Ta*wo*t)*(A*cos(wd*t)+B*sin(wd*t));
#up = -ta*wo*exp(-Ta*wo*t)*(A*cos(wd*t)+B*sin(wd*t)) + exp(-Ta*wo*t)*(-
wd*A*sin(wd*t)+wd*B*cos(wd*t));

#Condition initiale :
#u(t=0) = -10e-3
#up(t=0) = 0

A3 = -10e-3;
%B3 = (Ta*wo*u + up)/wd

#u = exp(-Ta*wo*t)*(A3*cos(wd*t))+((Ta*wo*u + up)/wd)*sin(wd*t);
function up = f(t,u,c,k)
up(1,1) = u(2);
up(2,1) = -c*u(2)-k*u(1);
end

% Conditions initiales
u01 = -10e-3;
u02 = 0;

#R�solution
[t,u] = ode45( @(t, u) f(t, u, c, k), [0:0.1:10*pi], [u01 u02]);

#Extraction de y1 et y2
u = u(:,1);
up = u(:,2);

figure(2)
#u fonction du temps
plot(t,u)

#up fonction du temps


hold on
plot(t,up)

#Plan de phase
hold on
plot(u,up)

Vous aimerez peut-être aussi