Calcul matrices de masse et de flexibilite
Outil interactif premium pour estimer la matrice de masse, la matrice de rigidite, la matrice de flexibilite et des indicateurs modaux elementaires d un systeme mecanique a 2 degres de liberte de type batiment cisaille ou modele concentre.
Resultats
Visualisation
Le graphique compare les grandeurs principales du modele choisi. Les frequences sont affichees en hertz, les deplacements statiques dans les unites derivees du systeme selectionne.
Guide expert sur le calcul des matrices de masse et de flexibilite
Le calcul des matrices de masse et de flexibilite constitue une base incontournable en dynamique des structures, en vibration des systemes mecaniques, en calcul par elements finis et en ingenierie sismique. Des que l on passe d un modele a un seul degre de liberte a un systeme a plusieurs degres de liberte, les grandeurs scalaires deviennent des objets matriciels. La masse n est plus seulement une valeur, elle est organisee sous forme de matrice, souvent notee M. De la meme maniere, la flexibilite, qui exprime la capacite d un systeme a se deformer sous charge, devient une matrice notee F, generalement egale a l inverse de la matrice de rigidite K lorsque celle ci est inversible. Comprendre ces relations permet de passer du raisonnement intuitif a une modelisation predictive solide.
Dans un cadre simple, la matrice de masse rassemble l inertie associee a chaque degre de liberte. Dans un modele concentre de batiment a deux niveaux, on place une masse sur chaque plancher. On obtient alors une matrice diagonale avec les termes m1 et m2 sur la diagonale principale. La matrice de rigidite, elle, traduit les liens elastiques entre les degres de liberte. Pour un batiment cisaille a deux etages, une expression classique est :
- M = [[m1, 0], [0, m2]]
- K = [[k1 + k2, -k2], [-k2, k2]]
- F = K-1
Ces matrices servent ensuite a plusieurs objectifs : calcul des deplacements statiques sous une charge laterale, determination des frequences propres, identification des modes de vibration, verification du confort vibratoire, estimation de la reponse sismique, et calibration des modeles numeriques. En pratique, les ingenieurs utilisent la relation dynamique matricielle M x¨ + C x· + K x = p(t), ou C represente la matrice d amortissement. Meme si l amortissement est un sujet important, l ossature du probleme repose d abord sur les matrices de masse et de rigidite, puis sur la flexibilite qui s en deduit.
Pourquoi la matrice de masse est essentielle
La masse controle l inertie du systeme. Plus la masse est elevee, plus il faut d effort pour imposer une acceleration. En dynamique, cela se traduit par des frequences naturelles plus basses a rigidite constante. Une matrice de masse mal evaluee provoque donc des erreurs directes sur les frequences propres et sur les reponses temporelles. En elements finis, on distingue souvent la matrice de masse concentree et la matrice de masse consistante. La premiere est simple, diagonale et robuste pour de nombreux calculs pratiques. La seconde est plus riche et plus fidele a la repartition continue de la masse, particulierement utile dans les modeles de poutres, plaques et coques.
Dans les structures civiles, la masse n est pas seulement le poids propre des materiaux. Il faut aussi integrer une partie des charges d exploitation, les equipements fixes, les facades, les cloisons et parfois des masses mobiles selon le cas d etude. En mecanique des machines, la matrice de masse inclut les masses concentrees, l inertie rotationnelle, les arbres, disques et organes tournants. Plus le modele est complexe, plus la definition precise des degres de liberte devient importante.
Que represente vraiment la matrice de flexibilite
La flexibilite est l inverse conceptuel de la rigidite. Si la rigidite mesure la resistance a la deformation, la flexibilite mesure la deformation obtenue pour une charge donnee. C est un outil tres utile pour interpreter les deplacements, les couplages entre degres de liberte et l influence de chaque ressort ou element structural. Dans un systeme lineaire elastique, si u = F p, alors la composante Fij represente le deplacement au degre de liberte i provoque par une charge unitaire appliquee au degre de liberte j. Cette lecture physique est extremement precieuse.
La matrice de flexibilite est particulierement importante dans les methodes basees sur les deplacements et dans les procedures d identification structurale. Elle aide aussi a detecter les zones faibles d un assemblage. Une augmentation locale de flexibilite peut signaler une perte de rigidite, une deterioration materielle, une fissuration ou une modification de condition d appui. Dans la surveillance de sante structurale, comparer les flexibilites mesurees et calculees est une approche classique.
Exemple de calcul pour un systeme a 2 degres de liberte
Considerons un modele a deux masses et deux ressorts. On choisit les degres de liberte lateraux aux niveaux 1 et 2. La matrice de masse est diagonale si l on concentre la masse sur les planchers. La matrice de rigidite prend en compte la rigidite du premier niveau k1 et celle du second niveau k2. A partir de la, on calcule la flexibilite par inversion de K. Ensuite, pour une charge laterale p, le vecteur de deplacement statique est u = F p. Enfin, les frequences naturelles se deduisent du probleme aux valeurs propres :
- Construire M
- Construire K
- Verifier que K est inversible
- Calculer F = K^-1
- Calculer les deplacements statiques u = F p
- Resoudre det(K – omega² M) = 0
- Convertir omega en frequence avec f = omega / 2pi
Cette logique simple est la colonne vertebrale de modeles beaucoup plus grands. Pour un systeme a n degres de liberte, les matrices deviennent n x n, mais les principes restent identiques. Ce qui change est surtout la taille numerique, la qualite des donnees d entree et la methode de resolution.
Ordres de grandeur utiles en ingenierie
Pour concevoir un modele credible, il faut s appuyer sur des valeurs physiques realistes. Le tableau suivant rassemble des ordres de grandeur largement utilises pour la densite et le module d elasticite de materiaux structuraux courants. Ces donnees sont des valeurs typiques employees en calcul preliminaire, avec des variations selon les formulations, les alliages, le taux d humidite, l age du beton et les conditions de mise en oeuvre.
| Materiau | Densite typique | Module d elasticite typique | Impact sur M et F |
|---|---|---|---|
| Acier structural | 7850 kg/m³ | 200 GPa | Masse elevee, rigidite elevee, flexibilite faible |
| Aluminium | 2700 kg/m³ | 69 GPa | Masse reduite, rigidite moderee, frequences souvent favorables |
| Beton ordinaire | 2300 a 2400 kg/m³ | 25 a 35 GPa | Masse importante, rigidite dependante de la fissuration |
| Bois de structure | 350 a 700 kg/m³ | 8 a 14 GPa | Masse faible, rigidite anisotrope, sensibilite aux liaisons |
Dans une perspective de vibration, les plages de frequences naturelles de certaines structures et systemes sont egalement instructives. Ces ordres de grandeur aident a verifier si un resultat numerique a du sens. Une frequence propre completement hors plage est souvent le signe d une erreur d unite, de rigidite, de masse, ou de conditions aux limites.
| Systeme | Frequence fondamentale typique | Commentaire |
|---|---|---|
| Plancher de bureau | 4 a 10 Hz | Zone sensible au confort pieton |
| Pont pietonnier leger | 1.5 a 5 Hz | Risque de resonance avec la marche |
| Batiment bas a moyen | 0.5 a 5 Hz | Depend de la hauteur, du systeme lateral et de la masse |
| Machine industrielle rigide | 10 a 60 Hz | Verification importante pour l isolation vibratoire |
Erreurs frequentes dans le calcul des matrices
Les erreurs les plus communes ne viennent pas de l algebra lineaire, mais des hypotheses de modelisation. Voici les plus courantes :
- Confondre masse et poids. En SI, la masse est en kg, le poids est en N.
- Melanger t, kg, kN/m et N/m sans conversion coherente.
- Oublier une contribution de masse, par exemple les equipements ou cloisons.
- Utiliser une rigidite excessive en ignorant la fissuration, le jeu ou la souplesse des appuis.
- Construire une matrice K non symetrique dans un probleme qui devrait l etre.
- Inverser une matrice proche de la singularite sans verifier les conditions aux limites.
- Employer un trop petit nombre de degres de liberte pour un phenomene distribue.
Une bonne pratique consiste a realiser plusieurs tests de coherence : verifier les dimensions physiques, comparer les ordres de grandeur a des references, controler la symetrie des matrices et examiner les vecteurs propres obtenus. Si les frequences augmentent quand les rigidites augmentent, et diminuent quand les masses augmentent, le modele suit deja une logique physique saine.
Difference entre rigidite, compliance et flexibilite
Dans la litterature, les termes compliance et flexibilite sont souvent utilises comme synonymes. En mecanique lineaire, ils designent la reponse deformable du systeme a une charge. La rigidite est son inverse dans un sens matriciel. En revanche, dans certains contextes de materiaux, on parle aussi de matrice de compliance constitutive, qui relie contraintes et deformations a l echelle du materiau. Il faut donc distinguer la flexibilite structurelle globale de la compliance locale de la loi de comportement. Les deux notions sont liees, mais elles n agissent pas au meme niveau du modele.
Application en dynamique sismique et vibration
En ingenierie parasismique, le couple matrice de masse et matrice de rigidite gouverne directement les modes propres et la participation modale. La masse modale effective permet de savoir quels modes contribuent vraiment a la reponse. Dans les structures regulieres, quelques modes dominent souvent la reponse laterale. Dans des structures irregulieres, la situation est plus complexe et exige un maillage plus fin ou des modeles plus riches. En vibration des machines, ces matrices servent aussi a eviter les resonances de service. Le but est de separer les frequences d excitation des frequences propres du systeme.
Quand l ingenieur calcule la flexibilite, il obtient une lecture immediate de la deformabilite sous force statique. Cette lecture est utile pour les verifications de derive, de deplacement en tete, de tassement relatif, ou de precision d un systeme mecanique. En couplant cette approche aux frequences naturelles, on obtient une vision a la fois statique et dynamique du comportement.
Comment exploiter l outil de calcul ci dessus
Le calculateur de cette page est volontairement centre sur un systeme a 2 DDL, ce qui en fait un excellent support pedagogique et un bon outil de pre dimensionnement. Entrez m1 et m2, puis les rigidites k1 et k2. La matrice de masse est construite automatiquement. La matrice de rigidite est ensuite assemblee, puis inversee pour obtenir la matrice de flexibilite. Les deplacements statiques sont calcules a partir des forces appliquees sur les deux degres de liberte. Enfin, les frequences naturelles sont determinees par resolution analytique du probleme quadratique pour ce cas particulier.
Le graphique permet ensuite de choisir une lecture simple :
- Comparer les masses et les flexibilites diagonales
- Visualiser les deplacements statiques des deux niveaux
- Observer les deux frequences naturelles du systeme
Cette approche n a pas vocation a remplacer un solveur complet d elements finis, mais elle permet de valider rapidement un ordre de grandeur, de tester la sensibilite a une variation de masse ou de rigidite, et de renforcer la comprehension des relations fondamentales entre inertie, rigidite et deformation.