Scribd
Importer
288 vues7 pages

Mini Projet MEF

Ce document décrit un projet de modélisation de vibrations d'une poutre en utilisant la méthode des éléments finis. Le document explique le développement des matrices de rigidité et de masse pour chaque élément, puis leur assemblage pour obtenir les matrices globales et calculer les valeurs propres.

Transféré par

ANDAHMOU
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 PDF, TXT ou lisez en ligne sur Scribd
288 vues7 pages

Mini Projet MEF

Ce document décrit un projet de modélisation de vibrations d'une poutre en utilisant la méthode des éléments finis. Le document explique le développement des matrices de rigidité et de masse pour chaque élément, puis leur assemblage pour obtenir les matrices globales et calculer les valeurs propres.

Transféré par

ANDAHMOU
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 PDF, TXT ou lisez en ligne sur Scribd
Vous êtes sur la page 1/ 7

MS-IM

Mini Projet : La modélisation de vibrations


d’une poutre en utilisant la méthode des
éléments finis

Encadré par : Mme.Abdoune

Réalisé par : EL OMAYRY Aymane


ANDAHMOU Abdellah
BOUTAFDA Mohamed

Université Mohammed V – Souissi MS-IM 2017/2018


ENSET
1
Objectif :

Ecrire un programme MATLAB sur la modélisation de vibrations


d’une poutre en utilisant la méthode des éléments finis.

On va traiter les axes suivants :

 Matrice de rigidité de chaque élément + Assemblage.


 Matrice de masse de chaque élément + Assemblage.
 Calcul des valeurs propres.

On utilisera un élément poutre en flexion simple comme exemple pour


exécuter le programme.

Données :

F = 18000N ; q = 104 N/m ; L = 2 m ;

S = 1 m2 ; E = 210 GPa ; I = 4.10-4 m4 ;

ρ = 8000 Kg/m3 ;

Université Mohammed V – Souissi MS-IM 2017/2018


ENSET
2
 On va discrétiser la poutre en trois éléments chaque élément a
les caractéristiques (E, I, L) comme montre la figure ci-dessus.

Etude de l’élément poutre :


L’équation obtenue par le principe fondamentale de la dynamique :

𝑑 2𝑣 𝑑 4𝑣
𝜌𝑆 2 + 𝐸𝐼 4 = 𝑓 (∗)
𝑑𝑡 𝑑𝑥

On cherchera les fonctions d’interpolations N1, N2, N3 et N4 à partir de


ces deux équations :

𝑣 (𝑥 ) = 𝑎 + 𝑏. 𝑥 + 𝑐. 𝑥 2 + 𝑑. 𝑥 3
𝑑𝑣
𝜃(𝑥 ) =
𝑑𝑥

On obtient finalement :
𝑣 (𝑥 ) = 𝑁(𝑥)1 . 𝑣1 + 𝑁(𝑥)2. 𝜃1 + 𝑁(𝑥)3. 𝑣2 + 𝑁(𝑥)4 . 𝜃2

Université Mohammed V – Souissi MS-IM 2017/2018


ENSET
3
La forme intégrale de l’équation (*) :

𝑑 2𝑢
𝑀𝑒 . { 2 } + 𝐾𝑒 . {𝑢} = {𝐹} (∗∗)
𝑑𝑡

𝑙
Or 𝑀𝑒 = ∫0 𝜌 𝑆〈𝑁(𝑥)〉𝑡 〈𝑁(𝑥)〉𝑑𝑥
𝑙 𝑑2 𝑡〈
𝑑2
𝐾𝑒 = ∫ 𝐸 𝐼 〈 𝑁(𝑥)〉 𝑁(𝑥)〉 𝑑𝑥
0 𝑑𝑥 2 𝑑𝑥 2
𝑁1 (𝑥)
𝑁 (𝑥)
𝑁(𝑥 ) = { 2 }
𝑁3 (𝑥)
𝑁4 (𝑥)

Etude des cas :


 Cas d’analyse statique :

Dans le cas stationnaire, u ne dépend pas de temps c’est-à-dire


𝑑2 𝑢
{ }=0
𝑑𝑡 2

Alors : l’équation (**) devient

𝐾𝑒 . {𝑢} = {𝐹 }
D’après cette équation on peut déterminer les déplacements des nœuds
de chaque élément dans le cas statique après l’assemblage de tous les
éléments. ( K : matrice de rigidité globale)

𝐾. {𝑢} = {𝐹 }

Université Mohammed V – Souissi MS-IM 2017/2018


ENSET
4
 Cas d’analyse modale :
Dans ce cas, u dépend de temps 𝑢 = 𝑢0 . cos(𝜔𝑡). L’équation (**)
devient :

−𝜔2 𝑀. {𝑢} + 𝐾. {𝑢} = {𝐹 }

(M : Matrice de masse globale)


Il faut faire la résolution du système suivant :

det(𝐾−𝜔2 𝑀) = 0

 Les w obtenus sont les fréquences et les modes propres.

Programme MATLAB :

Les sous-programmes :
 Matrice de masse pour un élément :

function ME=matrice_masse(R,S,l)

syms x L
N1=[1-3*(x/L)^2+2*(x/L)^3];
N2=[x-2*(x^2/L)+(x^3/L^2)];
N3=[3*(x/L)^2-2*(x/L)^3];
N4=[-(x^2/L)+(x^3/L^2)];
N=[N1,N2,N3,N4];
% Matrice M en fonction de x
M=transpose(N)*N;

% M(i,j) sont le composantes des matrices M=M(x)


% l'intégration sur [0,l]
% Pour M(1,1)
f11=@(x) ((2.*x.^3)./l.^3 - (3.*x.^2)./l.^2 + 1).^2;
M11=integral(f11,0,l);

Université Mohammed V – Souissi MS-IM 2017/2018


ENSET
5
% Pour M(2,2)
f22=@(x)(x - (2.*x.^2)./l + x.^3./l.^2).^2;
M22=integral(f22,0,l);

% Pour M(3,3)
f33=@(x)((3.*x.^2)./l.^2 - (2.*x.^3)./l.^3).^2;
M33=integral(f33,0,l);

% Pour M(4,4)
f44=@(x)(x.^2./l - x.^3./l.^2).^2;
M44=integral(f44,0,l);

% Pour M(1,2)
f12=@(x)((2.*x.^3)./l.^3 - (3.*x.^2)./l.^2 + 1).*(x -
(2.*x.^2)/l + x.^3./l.^2);
M12=integral(f12,0,l);

% Pour M(1,3)
f13=@(x)((3.*x.^2)./l.^2 -
(2.*x.^3)./l.^3).*((2.*x.^3)./l.^3 - (3.*x.^2)./l.^2 +
1);
M13=integral(f13,0,l);

% Pour M(1,4)
f14=@(x)-(x.^2./l - x.^3./l.^2).*((2.*x.^3)./l.^3 -
(3.*x.^2)./l.^2 + 1);
M14=integral(f14,0,l);

% Pour M(2,3)
f23=@(x)((3.*x.^2)./l.^2 - (2.*x.^3)./l.^3).*(x -
(2.*x.^2)./l + x.^3./l.^2);
M23=integral(f23,0,l);

% Pour M(2,4)
f24=@(x)-(x.^2./l - x.^3./l.^2).*(x - (2.*x.^2)./l +
x.^3./l.^2);
M24=integral(f24,0,l);

% Pour M(3,4)
f34=@(x)-(x.^2./l - x.^3./l.^2).*((3.*x.^2)./l.^2 -
(2.*x.^3)./l.^3);
M34=integral(f34,0,l);

% Matrice de masse ME=R*S*Me est symétrique


Me=[M11 M12 M13 M14;
M12 M22 M23 M24;
M13 M23 M33 M34;
M14 M24 M34 M44;];
ME=R*S*Me; % Matrice de masse pour un élément de (R,S,l)

Université Mohammed V – Souissi MS-IM 2017/2018


ENSET
6
 Matrice de rigidité pour un élément :
function [k, qe]=matrice_rigidite(I,E,qi,qj,xe)
% Calcul de [k] et [qe]=[fq+fT] d'une poutre 1D "flexion
plane"
% Syntaxe: [nom1, nom2]=poutre1d(I,E,qi,qj,xe)
% I = moment quadratique de section;
% E = module d'élasticité,
% qi, qj = charges/unité de longeur (N/m ou N/mm)(vers y
positif)
% xe = coord. [xi xj] des 2noeuds i et j
% Exemple : [k, qe]=poutre1d(1, 1e7, 0, 0, [0 100]);

xi=xe(1,1); xj=xe(1,2);
L=(xj-xi); % Longueur de l'élément
% matrice de rigidité pour un élément
k=(E*I/(L^3))*[ 12 6*L -12 6*L;
6*L 4*L*L -6*L 2*L*L;
-12 -6*L 12 -6*L;
6*L 2*L*L -6*L 4*L*L];
% Forces et moments de charge répartie
qe=[L*(7*qi+3*qj)/20; L*L*(3*qi+2*qj)/60;
L*(3*qi+7*qj)/20; -L*L*(2*qi+3*qj)/60];

 Le programme de l’exemple vous la trouvez ci-joint avec ce


document.

Fin

Université Mohammed V – Souissi MS-IM 2017/2018


ENSET
7

Vous aimerez peut-être aussi