Analyse structurelleProduction et robotique

FEA structurelle et modale

L’outil de simulation le plus utilisé en ingénierie — discrétiser une structure en éléments, résoudre pour le déplacement et la contrainte sous charge, et extraire ses fréquences propres et modes de vibration d’un problème aux valeurs propres. Contrainte de poutre encastrée, convergence de maillage, analyse modale, concentration de contrainte, flambement et un support réel.

DjiniousLabModélisation et simulation

Le champ de contrainte de von Mises d’une poutre encastrée chargée à son extrémité, maximal à l’encastrement et décroissant vers l’extrémité libre

déflexion de la poutre encastrée contre la théorie des poutres: 0.58%déflexion de la poutre encastrée contre la théorie des poutres
fréquences modales contre analytique: 0.0001%fréquences modales contre analytique
facteur de concentration de contrainte à un trou: 3.3facteur de concentration de contrainte à un trou
charge de flambement contre Euler: 0.00%charge de flambement contre Euler
notebooks de conception: 10notebooks de conception
exigences RÉUSSIES: 5 / 5exigences RÉUSSIES

Va-t-elle casser, et comment vibre-t-elle ?

L’analyse par éléments finis répond aux deux questions que doit passer chaque structure : de combien fléchit-elle et où casse-t-elle sous charge, et comment vibre-t-elle. Vous discrétisez la pièce en milliers de petits éléments, assemblez une matrice de raideur géante, et résolvez — pour le déplacement et la contrainte, ou pour les fréquences propres et les modes de vibration issus d’un problème aux valeurs propres. C’est l’outil de simulation le plus largement utilisé en ingénierie (ANSYS Mechanical, Nastran), et ce programme le construit depuis l’élément de base — validé, à chaque étape, face aux formules que vous pouvez vérifier à la main.

Là où le matériau travaille le plus dur.

Chargez une poutre encastrée à son extrémité et la contrainte n’est pas uniforme — elle se concentre à l’encastrement, sur les surfaces supérieure et inférieure où la flexion est maximale, et tombe à zéro à l’axe neutre et à l’extrémité libre. Le solveur par éléments finis calcule le champ complet de contrainte de von Mises qui montre exactement cela, et la déflexion à l’extrémité qu’il prédit se situe à 0,58 % près de la formule de poutre d’Euler-Bernoulli. Le champ vous dit où ajouter du matériau ; la concordance vous dit de lui faire confiance.

Le champ de contrainte de von Mises d’une poutre encastrée chargée à son extrémité, concentré à l’encastrement — Le héros : le champ de contrainte de von Mises d’une poutre encastrée chargée à son extrémité. La contrainte se concentre sur les fibres supérieure et inférieure de l’encastrement — là où une poutre casse — et décroît vers l’extrémité libre. La déflexion prédite correspond à la théorie des poutres d’Euler-Bernoulli à 0,58 % près. C’est à cela que sert la FEA : voir où une pièce travaille le plus dur.

Chaque structure a des notes qu’elle veut faire résonner.

Au-delà de la résistance statique, une structure a des fréquences propres — les régimes auxquels elle veut vibrer — et un mode de vibration pour chacune. Trouvez-les en résolvant un problème aux valeurs propres généralisé sur les matrices de raideur et de masse, et vous obtenez les résonances qu’il faut éviter pour le trafic d’un pont ou les vibrations d’un moteur. Sur une poutre encastrée, les fréquences calculées correspondent aux valeurs analytiques à quatre décimales près, et les modes de vibration sont les formes de flexion classiques des manuels.

Les quatre premiers modes de flexion d’une poutre encastrée, avec un nombre croissant de nœuds — Analyse modale : les quatre premiers modes de vibration d’une poutre encastrée, chacun avec un nœud de plus que le précédent, issus d’un problème aux valeurs propres sur les matrices de raideur et de masse. Les fréquences propres correspondent à la théorie analytique des poutres à 0,0001 % près. Connaître ces résonances, c’est ce qui empêche une structure de se secouer jusqu’à se détruire.

Trous, coins et colonnes élancées.

Les pièces réelles cassent à leurs détails. Un trou dans une plaque chargée triple la contrainte locale — le facteur de concentration de contrainte d’environ 3 que prédit la solution analytique de Kirsch, et que le solveur retrouve à mesure que le trou rétrécit. Une colonne élancée ne cède pas, elle flambe, à la charge d’Euler que le solveur reproduit au chiffre près. Et le coin rentrant d’un support en L est un concentrateur de contrainte qui détermine son coefficient de sécurité. Ce sont les détails qui décident si une pièce survit.

Une plaque chargée avec un trou, montrant une concentration de contrainte au bord du trou d’environ trois fois la contrainte en champ lointain — Concentration de contrainte : un trou dans une plaque chargée élève la contrainte locale d’environ 3× le champ lointain — le facteur de concentration de contrainte que connaît chaque concepteur, retrouvé ici (Kt ≈ 3,3, convergeant vers l’analytique 3 à mesure que le trou rétrécit). La raison pour laquelle c’est un trou, pas la masse, qui fissure généralement une pièce en premier.

Un treillis de pont Warren avec des membrures colorées selon la tension et la compression sous une charge à mi-portée — Un treillis, l’autre classique : un pont Warren dont chaque membrure est colorée selon qu’elle est en tension ou en compression sous une charge à mi-portée. La même FEM, appliquée à une ossature de membrures élancées — le pain quotidien du génie des structures.

D’un champ de contrainte à un coefficient de sécurité.

Assemblez le tout sur un composant réel — un support en L encastré le long du haut avec une charge sur son bras. Le solveur trouve la contrainte de pointe au coin rentrant, et diviser la limite d’élasticité du matériau par cette contrainte donne le coefficient de sécurité, le chiffre unique sur lequel repose une revue de conception. Balayez la charge et vous lisez la charge admissible : jusqu’où vous pouvez pousser la pièce avant que la marge de sécurité ne s’épuise. C’est la FEA qui fait son véritable travail — transformer une géométrie et une charge en un feu vert ou rouge.

Le champ de contrainte de von Mises d’un support en L avec un pic au coin rentrant — Une pièce réelle : le champ de von Mises d’un support en L, avec la contrainte de pointe au coin rentrant — 78,7 MPa sous charge, un coefficient de sécurité de 3,2 contre une limite d’élasticité de l’acier. Balayez la charge et vous obtenez l’admissible. La transformation géométrie-plus-charge en feu vert/rouge que la FEA existe pour fournir.

Chaque valeur re-dérivée au moment de la validation finale.

Le notebook V&V reconstruit le solveur depuis zéro et re-dérive chaque exigence face à la théorie sous forme close, en affichant un tableau RÉUSSI/ÉCHEC.

Résultat

Flèche de la poutre encastrée-libre vs théorie des poutres : 0,58 %
Fréquences modales vs solution analytique : 0,0001 %
Convergence du maillage : monotone
Facteur de concentration de contrainte : 3,3 → 3
Charge de flambement vs Euler : 0,00 %

Exigence

Flèche de la poutre encastrée-libre vs théorie des poutres : R-01 < 3 %
Fréquences modales vs solution analytique : R-02 < 0,1 %
Convergence du maillage : R-03
Facteur de concentration de contrainte : R-04
Charge de flambement vs Euler : R-05 < 1 %

Un vrai solveur, à fidélité pédagogique.

Le livrable est dix notebooks et le dossier — un champ structurel discrétisé avec un problème aux valeurs propres, pas un réseau acausal à paramètres localisés, donc pas de bloc personnalisé ni de canevas. Il utilise des éléments Q4 en contraintes planes et des éléments de poutre d’Euler-Bernoulli, une élasticité linéaire en petites déformations, un solveur direct dense sur des maillages de taille modeste, et une récupération des contraintes au point central — sans contact, plasticité ni grandes déformations. Ce que le programme démontre, c’est la capacité qui manquait le plus au portefeuille : la mécanique des structures par éléments finis, à la fois en contrainte statique et en vibration modale, validée à chaque étape par rapport à la théorie sous forme close — et c’est le même solveur EF que fait tourner dans sa boucle le programme phare d’optimisation topologique, le pendant en champ structurel des solveurs de champ CFD et électromagnétique.