FEM Geometry Preparation and Meshing/fr: Difference between revisions

From FreeCAD Documentation
No edit summary
No edit summary
 
(36 intermediate revisions by 2 users not shown)
Line 7: Line 7:
|Author=[http://www.freecadweb.org/wiki/index.php?title=User:NewJoker NewJoker]
|Author=[http://www.freecadweb.org/wiki/index.php?title=User:NewJoker NewJoker]
|FCVersion=0.21 ou plus récente
|FCVersion=0.21 ou plus récente
|SeeAlso=[[FEM_Workbench/fr|Atelier FEM]]
}}
}}


Line 42: Line 43:
== Validité de la géométrie ==
== Validité de la géométrie ==


La géométrie utilisée pour la simulation par éléments finis doit être valide. Plus important encore, il ne doit pas y avoir d'intersections. Il s'agit d'un problème courant qui survient souvent lorsque des assemblages sont modélisés sans contraintes appropriées entre les pièces. L'outil [[Part_SectionCut/fr|Part Coupe persistante]] peut aider à trouver de telles interférences entre les pièces. Bien entendu, l'outil [[Part_Fuse/fr|Part Union]] peut aider à les résoudre si elles sont intentionnelles. D'autres problèmes liés à la géométrie (tels que les géométries non pliantes, les arêtes ou faces redondantes, etc.) doivent également être résolus avant de procéder au maillage. L'outil [[Part_CheckGeometry/fr|Part Vérifier la géométrie]] peut être utile, mais les contrôles visuels sont également importants. Lors de la préparation d'une simulation utilisant des éléments solides et en cas de doute quant à savoir si la pièce est réellement solide ou s'il s'agit simplement d'une coque fermée, les outils susmentionnés (Part Coupe persistante et l'onglet Contenu de la forme de Part Vérifier la géométrie) peuvent permettre de clarifier la situation.
La géométrie utilisée pour la simulation par éléments finis doit être valide. Plus important encore, il ne doit pas y avoir d'intersections. Il s'agit d'un problème courant qui survient souvent lorsque des assemblages sont modélisés sans contraintes appropriées entre les pièces. L'outil [[Part_SectionCut/fr|Part Coupe persistante]] peut aider à trouver de telles interférences entre les pièces. Bien entendu, l'outil [[Part_Fuse/fr|Part Union]] peut aider à les résoudre si elles sont intentionnelles. D'autres problèmes liés à la géométrie (tels que les géométries non pliantes, les arêtes ou faces redondantes, etc.) doivent également être résolus avant de procéder au maillage. L'outil [[Part_CheckGeometry/fr|Part Vérifier la géométrie]] peut être utile, mais les contrôles visuels sont également importants. Lors de la préparation d'une simulation utilisant des éléments solides et en cas de doute quant à savoir si la pièce est réellement solide ou s'il s'agit simplement d'une coque fermée, les outils susmentionnés (Part Coupe persistante et l'onglet ''Contenu de la forme'' de Part Vérifier la géométrie) peuvent permettre de clarifier la situation.


<span id="Geometry_simplification"></span>
<span id="Geometry_simplification"></span>
Line 65: Line 66:
{{Caption|Géométrie du support simplifiée en utilisant uniquement l'outil Part Supprimer la fonctionnalité}}
{{Caption|Géométrie du support simplifiée en utilisant uniquement l'outil Part Supprimer la fonctionnalité}}


Dans le cas des assemblages (plus d'informations à ce sujet dans l'une des sections suivantes), certaines pièces peuvent souvent être exclues des simulations et remplacées par des conditions aux limites comme si elles étaient attachées aux pièces à analyser. Cette approche est valable si les pièces exclues sont significativement plus rigides (en termes de rigidité structurelle, en tenant compte non seulement de l'élasticité du matériau mais aussi de la géométrie de la pièce) que les pièces à analyser auxquelles elles étaient connectées. En effet, les conditions aux limites fixes introduisent de la rigidité (comme si la pièce analysée était attachée à un composant infiniment rigide) et les supports flexibles (comme les ressorts) ne sont pas disponibles dans l'atelier FEM de FreeCAD.
Dans le cas des assemblages (plus d'informations à ce sujet dans l'une des sections suivantes), certaines pièces peuvent souvent être exclues des simulations et remplacées par des conditions aux limites comme si elles étaient attachées aux pièces à analyser. Cette approche est valable si les pièces exclues sont significativement plus rigides (en termes de rigidité structurelle, en tenant compte non seulement de l'élasticité du matériau mais aussi de la géométrie de la pièce) que les pièces à analyser auxquelles elles étaient connectées. En effet, les conditions aux limites fixes introduisent de la rigidité (comme si la pièce analysée était attachée à un composant infiniment rigide) et les supports flexibles comme les éléments ressorts ne sont pas disponibles dans l'atelier FEM de FreeCAD lorsque l'on utilise CalculiX (Elmer a un [[FEM_ConstraintSpring/fr|FEM Ressort]]).


La simplification de la géométrie pour la simulation par éléments finis peut également impliquer de la couper dans l'un des plans de symétrie afin d'utiliser l'hypothèse de symétrie plane dans l'analyse. Cette hypothèse n'est valable que lorsque tous les aspects suivants du modèle présentent une symétrie dans un plan donné :
La simplification de la géométrie pour la simulation par éléments finis peut également impliquer de la couper dans l'un des plans de symétrie afin d'utiliser l'hypothèse de symétrie plane dans l'analyse. Cette hypothèse n'est valable que lorsque tous les aspects suivants du modèle présentent une symétrie dans un plan donné :
Line 84: Line 85:
== Partitionnement de la géométrie ==
== Partitionnement de la géométrie ==


Le partitionnement est une division de la géométrie en segments plus petits. Dans d'autres logiciels, il est généralement utilisé pour permettre le maillage hexagonal, mais dans FreeCAD, il peut être utile pour d'autres raisons. La principale application du partitionnement est lorsqu'une charge (ou une condition limite) doit être appliquée uniquement à une région sélectionnée de la surface de la pièce. La façon la plus simple d'y parvenir est de créer une esquisse avec un contour approprié sur cette face et d'utiliser l'outil [[Part_BooleanFragments/fr|Part Fragments booléens]] pour diviser la face avec l'esquisse. Une autre raison d'utiliser le partitionnement est lorsque plusieurs matériaux doivent être appliqués à une seule pièce (sans avoir à utiliser plusieurs pièces connectées entre elles). Dans ce cas, le partitionnement peut être effectué à l'aide d'un [[PartDesign_Plane/fr|plan]] et de l'outil Fragments booléens en mode ''Compsolid''. Le partitionnement peut également être utilisé pour créer des régions pour [[FEM_MeshRegion/fr|affiner le maillage]].
<div lang="en" dir="ltr" class="mw-content-ltr">
So-called partitioning is a division of the geometry into smaller segments. In other software, it's commonly used to allow hex meshing but in FreeCAD it can be useful for other reasons too. The main application of partitioning is when a load (or a boundary condition) has to be applied only to a selected region of the part's surface. The easiest way to achieve it is to create a sketch with a proper contour on that face and use the [[Part_BooleanFragments|Part Boolean Fragments]] tool to split the face with the sketch. Another reason for partitioning is when multiple materials have to be applied to a single part (without having to use multiple parts connected with each other). Then partitioning can be done using a [[PartDesign_Plane|datum plane]] and Boolean Fragments tool with the ''Compsolid'' mode.
</div>


<div lang="en" dir="ltr" class="mw-content-ltr">
[[Image:FEM_partition.JPG|400px]]
[[Image:FEM_partition.JPG|400px]]
</div>


{{Caption|Pièce avec une cloison de séparation pour l'application d'une charge ou d'une condition limite}}
<div lang="en" dir="ltr" class="mw-content-ltr">
{{Caption|Part with a face partition for load or boundary condition application}}
</div>


<span id="Assembly_geometries"></span>
<div lang="en" dir="ltr" class="mw-content-ltr">
== Assembly geometries ==
== Assemblage de géométries ==
</div>


L'une des principales limitations actuelles de l'atelier FEM est que les maillages multiples ne sont pas pris en charge. En pratique, cela signifie qu'il n'est pas possible de mailler chaque pièce de l'assemblage individuellement, puis de connecter les pièces avec les contraintes appropriées pour l'analyse. Il est donc nécessaire de créer un objet unique contenant toutes les pièces de l'assemblage et de le mailler. Il existe plusieurs options différentes, toutes basées sur les [[Part_Module/fr#Booléen|Part Outils booléens]]. Le choix dépend de l'effet désiré : si chaque pièce/volume et leurs limites doivent être sélectionnables (par exemple pour l'affectation de matériaux ou les définitions des conditions limites agissant sur les faces internes) ou non :
<div lang="en" dir="ltr" class="mw-content-ltr">
* [[Part_Fuse/fr|Part Union]] : fusionne les pièces, rendant impossible leur sélection individuelle, par exemple pour les définitions de matériaux,
One of the current major limitations of the FEM workbench is that multiple meshes are not supported. In practice, this means that one cannot mesh each part of the assembly individually and then connect the parts with proper constraints for the analysis. Instead, it's necessary to create a single object containing all the parts of the assembly and mesh it. There are several different options here, all relying on [[Part_Module#Boolean|Part boolean tools]]. The choice depends on the desired effect - whether the individual parts/volumes and their boundaries should be selectable (e.g. for material assignments or definitions of boundary conditions acting on internal faces) or not:
* [[Part_Compound/fr|Part Composé]] : crée un objet composé, ce qui permet de sélectionner chaque pièce,
* [[Part_Fuse|Part Fuse]] - merges the parts, making it impossible to select them individually e.g. for material definitions,
* [[Part_JoinConnect/fr|Part Connecter]] : fonctionne comme Part Union, mais fusionne les parties, ce qui rend impossible leur sélection individuelle,
* [[Part_Compound|Part Compound]] - creates a compound object, making it possible to select individual parts,
* [[Part_BooleanFragments/fr|Part Fragments booléens]] : fonctionne comme Part Composé, ce qui permet de sélectionner chaque partie.
* [[Part_JoinConnect|Part JoinConnect]] - works like Part Fuse, merges the parts, making it impossible to select them individually,
Il est important de mentionner que si les pièces se touchent, un maillage continu sera créé sur l'objet booléen et aucune contrainte ne sera nécessaire pour la simulation. S'il y a un petit espace entre les pièces (ou une intersection à l'intérieur d'un Part Composé), le maillage ne sera pas continu et des contraintes comme la [[FEM_ConstraintTie/fr|Contrainte de liaison]] ou la [[FEM_ConstraintContact/fr|Contrainte de contact]] seront nécessaires. L'exécution d'une analyse en fréquence est un bon moyen de révéler si le maillage est continu ou non. Si les pièces ne sont pas connectées, les premières formes des modes avec déformation visualisées à l'aide du [[FEM_PostFilterWarp/fr|filtre des déformations]] montreront une séparation, les pièces "s'envoleront".
* [[Part_BooleanFragments|Part BooleanFragments]] - works like Part Compound, making it possible to select individual parts.
It's important to mention that if the parts are touching, a continuous mesh will be created on the boolean object and no constraints will be needed for the simulation. If there's even a small gap between the parts, the mesh won't be continuous and constraints like [[FEM_ConstraintTie|tie]] or [[FEM_ConstraintContact|contact]] will be needed. Running a frequency analysis is a good way to reveal if the mesh is continuous or not - if the parts are not connected, the first mode shapes with deformation visualized using [[FEM_PostFilterWarp|Warp filter]] will show separation - the parts will "fly away".
</div>


<div lang="en" dir="ltr" class="mw-content-ltr">
[[Image:FEM_modal_separation.JPG|400px]]
[[Image:FEM_modal_separation.JPG|400px]]
</div>


{{Caption|La forme du premier mode d'une analyse de fréquence visualisée avec le filtre des déformations - deux cubes avec un petit écart initial ont été analysés.}}
<div lang="en" dir="ltr" class="mw-content-ltr">
{{Caption|The first mode shape of a frequency analysis visualized with the Warp filter - two cubes with a small initial gap were analyzed}}
</div>


La sélection des régions internes (faces/volumes) peut s'avérer délicate. Elle peut être nécessaire pour l'application de différents matériaux, charges ou conditions aux limites (en particulier dans les analyses thermiques et électromagnétiques). Plusieurs méthodes sont possibles :
<div lang="en" dir="ltr" class="mw-content-ltr">
* activer [[Std_ToggleClipPlane/fr|plan de coupe]] pour la sélection et le choix des faces internes,
Selection of internal regions (faces/volumes) can be tricky. It might be needed for the application of different materials, body loads or boundary conditions (especially in thermal and electromagnetic analyses). Several ways are possible:
* masquer l'objet booléen, ne montrer qu'une des parties auxquelles il a été appliqué et la sélectionner,
* enabling [[Std_ToggleClipPlane|clipping plane]] for the time of selection and picking internal faces,
* sélectionner un autre objet externe et éditer la propriété ''References'' dans l'onglet ''Données'' d'un élément d'analyse donné (nécessite la spécification manuelle du numéro de l'objet géométrique).
* hiding the boolean object, showing only one of the parts it was applied to and selecting it,
* selecting another, external object and editing the References property in the Data tab of a given analysis feature (requires manual specification of the geometric object's number).
</div>


<span id="Meshing_basics"></span>
<div lang="en" dir="ltr" class="mw-content-ltr">
== Meshing basics ==
== Les bases du maillage ==
</div>


Un maillage trop grossier est l'une des sources les plus courantes d'inexactitudes et d'autres problèmes dans les simulations par éléments finis. Il s'agit souvent d'une erreur partielle des paramètres du maillage automatique. Ils génèrent généralement des maillages très grossiers et inadaptés lorsque la taille des éléments n'est pas spécifiée manuellement, mais laissée avec une valeur par défaut. Il faut toujours connaître les dimensions approximatives de la pièce, en particulier la taille de la plus petite caractéristique pertinente (l'outil [[Part_Measure_Linear/fr|Part Mesure linéaire]] peut aider à cet égard) et spécifier la taille maximale appropriée de l'élément sur cette base. Lors du maillage avec Gmsh, il existe également un paramètre de taille minimale des éléments qui peut empêcher la création d'éléments trop petits autour de petites caractéristiques géométriques, ce qui peut conduire à des maillages inutilement denses (et parfois même à un plantage ou à un gel de FreeCAD lorsqu'on essaie de générer de tels maillages). D'une manière générale, il est préférable de commencer par un maillage plus grossier (moins long à générer), de voir à quoi il ressemble (une certaine expérience est nécessaire) et de l'affiner si nécessaire. Il est souvent judicieux d'utiliser un maillage dense uniquement autour des zones d'intérêt (endroits présentant des gradients/concentrations de contraintes importants - encoches) et un maillage relativement grossier à l'écart de ces zones. De cette manière, le nombre d'éléments peut être considérablement réduit, ce qui permet de diminuer les temps de résolution. Le raffinement local du maillage est défini à l'aide de [[FEM_MeshRegion/fr|FEM Région de maillage FEM]].
<div lang="en" dir="ltr" class="mw-content-ltr">
Too coarse mesh is one of the most common sources of inaccuracies and other issues in FEM. It's often a partial fault of automatic mesher settings - they typically generate very coarse, unsuitable meshes when the element size is not manually specified but left with a default value. One should always know the approximate dimensions of the part, especially the size of the smallest relevant feature ([[Part_Measure_Linear|Part Measure Linear]] tool can help with that) and specify the proper maximum element size based on that. It's better to start with a coarser mesh (taking less time to generate), see what it looks like (some experience is necessary) and refine it if necessary. It often makes sense to use dense mesh only around the areas of interest (locations with large stress gradients/concentrations - notches) and relatively coarse mesh away from them. This way, the number of elements can be significantly reduced, leading to shorter solving times. Local mesh refinement is defined using [[FEM_MeshRegion|FEM MeshRegion]].
</div>


<div lang="en" dir="ltr" class="mw-content-ltr">
[[Image:FEM_default_mesh.PNG|400px]]
[[Image:FEM_default_mesh.PNG|400px]]
</div>


{{Caption|Par défaut, maillage trop grossier}}
<div lang="en" dir="ltr" class="mw-content-ltr">
{{Caption|Default, too coarse mesh}}
</div>


<div lang="en" dir="ltr" class="mw-content-ltr">
[[Image:FEM_globally_refined_mesh.PNG|400px]]
[[Image:FEM_globally_refined_mesh.PNG|400px]]
</div>


{{Caption|Maillage amélioré globalement}}
<div lang="en" dir="ltr" class="mw-content-ltr">
{{Caption|Globally refined mesh}}
</div>


<div lang="en" dir="ltr" class="mw-content-ltr">
[[Image:FEM_locally_refined_mesh.PNG|400px]]
[[Image:FEM_locally_refined_mesh.PNG|400px]]
</div>


{{Caption|Maillage amélioré localement}}
<div lang="en" dir="ltr" class="mw-content-ltr">
{{Caption|Locally refined mesh}}
</div>


Le choix du type d'élément n'est pas facile et dépend de nombreux facteurs, mais la règle générale est que les éléments hexaédriques et quadruples sont préférables aux éléments tétraédriques et triangulaires. Cependant, les géométries complexes ne peuvent pas être maillées avec des éléments hexaédriques et FreeCAD ne peut pas les générer du tout (seuls les maillages quadruples peuvent être générés sur les surfaces, voir [https://forum.freecad.org/viewtopic.php?t=20351 ce fil de discussion]). Les éléments hexaédriques peuvent être importés à partir de maillages externes comme [https://gmsh.info Gmsh] et utilisés dans l'atelier FEM comme le montre [https://www.youtube.com/watch?v=vylt24G7qj4&t=932s cette vidéo].
<div lang="en" dir="ltr" class="mw-content-ltr">
The choice of element type is not easy and depends on many factors but the general rule is that hexahedral and quad elements are preferable over tetrahedral and triangular ones. However, complex geometries can't be meshed with hexahedral elements and FreeCAD can't generate them at all (only quad meshes can be generated on surfaces - see [https://forum.freecad.org/viewtopic.php?t=20351 this forum thread]). Hexahedral elements can be imported from external meshers like [https://gmsh.info Gmsh] and used in the FEM workbench as shown in [https://www.youtube.com/watch?v=vylt24G7qj4&t=932s this video].
</div>


Le choix de l'ordre des éléments (premier ou second) dépend des conditions d'analyse, mais dans la plupart des cas, les éléments de second ordre sont préférables. Ceci est particulièrement vrai pour les éléments triangulaires et tétraédriques, leurs versions de premier ordre (linéaires) ne sont normalement pas recommandées pour un usage régulier et ne devraient être utilisées que comme éléments de remplissage dans des régions de faible importance. Cependant, comme FreeCAD ne peut pas générer d'éléments hexaédriques, les tétraèdres linéaires peuvent être utilisés dans certains cas, si les maillages sont suffisamment denses. En particulier lors d'analyses avec une [[FEM_ConstraintContact/fr|contrainte de contact]].
<div lang="en" dir="ltr" class="mw-content-ltr">
The choice of element order (first or second) depends on the analysis conditions but in most cases, second-order elements are preferred. This is particularly the case with triangular and tetrahedral elements - their first-order (linear) versions are normally not recommended for regular usage and they should be used only as filler elements in regions of low importance. However, since FreeCAD can't generate hexahedral elements, linear tetrahedrons can be used in some cases, if the meshes are dense enough. Especially when performing analyses with [[FEM_ConstraintContact|contact]].
</div>


<span id="Negative_Jacobians"></span>
<div lang="en" dir="ltr" class="mw-content-ltr">
== Jacobiens négatifs ==
== Mesh convergence studies ==
</div>


Si les règles ci-dessus sont respectées (notamment en ce qui concerne la validité de la géométrie, le defeaturing et la sélection de la taille des éléments), le maillage devrait être généré correctement. Cependant, dans certains cas, la géométrie ne peut pas être trop simplifiée, ou la procédure de modélisation est appropriée mais conduit quand même à des arêtes et des faces de petite taille. Dans ce cas, le maillage avec des éléments de second ordre peut échouer en raison de jacobiens négatifs. La raison en est que les maillages doivent suivre le modèle CAO et placer les nœuds médians des éléments de second ordre sur la géométrie. Dans le cas de formes plus complexes, les éléments peuvent être tellement étirés qu'ils deviennent inversés. Le jacobien est l'une des mesures de qualité de maillage les plus courantes. Il représente l'écart de l'élément par rapport à la forme idéale. Elle devient négative lorsque l'élément se retourne (devient inversé), soit en raison d'une déformation importante au cours de l'analyse (non prise en compte ici), soit en raison des problèmes de maillage susmentionnés. Les jacobiens négatifs dans l'atelier FEM de FreeCAD peuvent être signalés par Gmsh ou par CalculiX. Leurs emplacements dans le maillage sont mis en évidence lorsque les analyses CalculiX sont soumises en utilisant le bouton [[FEM_SolverRun/fr|Lancer les calculs du solveur]]. Les conseils suivants peuvent aider à les éliminer :
<div lang="en" dir="ltr" class="mw-content-ltr">
* Définir la propriété ''Second Order Linear'' de l'objet [[FEM_MeshGmshFromShape/fr|Maillage FEM à partir d'une forme avec Gmsh]] à ''true'', ceci a pour résultat que les nœuds du milieu des éléments de second ordre sont simplement ajoutés au milieu des bords droits (initialement) des éléments de premier ordre, sans les attacher à la géométrie et résout le problème dans la plupart des cas,
Mesh convergence studies are recommended in all serious projects requiring accurate results. The reason is that the results can change a lot, approaching correct values when the mesh is refined. The following approach should be used:
* Utiliser Netgen au lieu de Gmsh. Netgen est connu pour être moins sujet aux problèmes de jacobiens négatifs mais ne les signale pas non plus, de sorte que l'utilisateur ne peut s'en rendre compte qu'au moment de soumettre l'analyse,
# After obtaining the first results and noting them (usually maximum von Mises stress, von Mises stress at a given location and maximum displacement), refine the mesh (globally or better locally - with [[FEM_MeshRegion|FEM MeshRegion]]) and rerun the simulation.
* Réduire davantage la taille des éléments,
# Check the results and note their new values. If they differ significantly from the initial results, refine the mesh further and rerun the analysis.
* Exporter la géométrie, essayer de la mailler dans l'interface graphique de Gmsh ou Netgen (NGSolve) ou dans un autre maillage autonome (comme Salome_Meca). Ces outils ont des fonctions supplémentaires qui peuvent aider à se débarrasser des jacobiens négatifs (par exemple, Gmsh a ce que l'on appelle des "outils d'ordre élevé"),
# Repeat the process if the results still change (usually grow) significantly with mesh refinement.
* Utiliser des éléments de premier ordre. Cela ne doit être fait qu'en dernier recours car les tétraèdres de premier ordre sont connus pour leur imprécision.
It usually helps to create a plot with a given result vs mesh density. This way it's easier to notice when the results start to converge. The acceptable difference in results between two runs is usually around a few percent (e.g. below 5%).
</div>


Indépendamment de ces conseils, il est important de souligner une fois de plus que les jacobiens négatifs sont généralement dus à des approches de modélisation désordonnées et à un manque de préparation de la géométrie pour l'analyse (ce qui est particulièrement courant avec les modèles STEP téléchargés à partir de divers sites web). Même si le maillage est finalement généré dans de tels cas, les résultats risquent d'être de mauvaise qualité (rappelons la règle GIGO mentionnée dans le premier paragraphe). Par conséquent, le nettoyage de la géométrie et la préparation pour la FEM devraient toujours être la priorité.
<div lang="en" dir="ltr" class="mw-content-ltr">

In some cases, it may happen that the maximum stress will be growing indefinitely regardless of how dense the mesh will be. Such a non-physical effect is known as a stress singularity. It may occur due to the following reasons:
<span id="Mesh_convergence_studies"></span>
* concentrated forces applied to solid and shell models,
== Études sur la convergence des maillages ==
* boundary conditions applied to points (individual nodes),

* sharp corners,
Les études de convergence de maillage sont recommandées dans tous les projets sérieux nécessitant des résultats précis. En effet, les résultats peuvent changer considérablement et s'approcher des valeurs correctes lorsque le maillage est affiné. L'approche suivante doit être utilisée :
* contact occurring at a corner.
# Après avoir obtenu les premiers résultats et les avoir notés (généralement la contrainte de von Mises maximale, la contrainte de von Mises à un endroit donné et le déplacement maximal), affiner le maillage (globalement ou mieux localement avec [[FEM_MeshRegion/fr|FEM Région de maillage FEM]]) et réexécuter la simulation.
Typical ways of dealing with stress singularities are:
# Vérifiez les résultats et notez leurs nouvelles valeurs. S'ils diffèrent significativement des résultats initiaux, affinez encore le maillage et relancez l'analyse.
* applying loads and boundary conditions to small areas instead of points - see the section about partitioning above,
# Répétez le processus si les résultats changent encore (généralement en augmentant) de manière significative avec l'affinement du maillage.
* adding small fillets to sharp corners (an exception to the rule of omitting small fillets when simplifying the geometry for FEM),
Il est généralement utile de créer un graphique avec un résultat donné en fonction de la densité du maillage. De cette façon, il est plus facile de remarquer quand les résultats commencent à converger. La différence acceptable dans les résultats entre deux exécutions est généralement de l'ordre de quelques pour cent (par exemple, moins de 5 %).
* including plasticity in material behavior to enable stress redistribution and limit the stresses to values allowed by the plasticity definition while observing the proper level of yielding (plastic strain),

* ignoring singularities and reading stresses away from them if possible (based on St. Venant’s Principle).
Dans certains cas, il peut arriver que la contrainte maximale augmente indéfiniment, quelle que soit la densité du maillage. Un tel effet non physique est connu sous le nom de singularité de contrainte. Il peut se produire pour les raisons suivantes
</div>
* les forces concentrées appliquées aux modèles de solides et de coques,
* conditions limites appliquées à des points (nœuds individuels),
* angles vifs,
* contact survenant dans un coin.
Les moyens typiques de traiter les singularités de contrainte sont les suivants :
* appliquer des charges et des conditions aux limites à de petites zones plutôt qu'à des points - voir la section sur le partitionnement ci-dessus,
* ajouter de petits filets aux angles vifs (une exception à la règle qui consiste à omettre les petits filets lors de la simplification de la géométrie pour les éléments finis),
* inclure la plasticité dans le comportement du matériau pour permettre la redistribution des contraintes et limiter les contraintes aux valeurs autorisées par la définition de la plasticité tout en observant le niveau approprié de déformation (déformation plastique),
* ignorer les singularités et lire les contraintes loin d'elles si possible (sur la base du principe de St. Venant).


<div lang="en" dir="ltr" class="mw-content-ltr">
[[Image:FEM_mesh_convergence.PNG|400px]]
[[Image:FEM_mesh_convergence.PNG|400px]]
</div>


{{Caption|Tracés de convergence de maillage typiques :<br>
<div lang="en" dir="ltr" class="mw-content-ltr">
- le déplacement ('''courbe verte''') converge rapidement,<br>
{{Caption|Typical mesh convergence plots - displacement ('''green curve''') converges quickly, maximum stress at a notch like a hole ('''blue curve''') needs more iterations of mesh refinement to converge while maximum stress at a sharp corner with a fixed boundary condition ('''red curve''') doesn't converge at all - stress singularity occurs and a small fillet would have to be added and the connection would have to be modeled in more realistic, flexible way to avoid this behavior.}}
- la contrainte maximale au niveau d'une entaille comme un trou ('''courbe bleue''') nécessite davantage d'itérations de raffinement du maillage pour converger, <br>
</div>
- tandis que la contrainte maximale au niveau d'un angle aigu avec une condition limite fixe ('''courbe rouge''') ne converge pas du tout. Une singularité de contrainte se produit et un petit congé devrait être ajouté et la connexion devrait être modélisée de manière plus réaliste et plus flexible pour éviter ce comportement.<br>}}


{{FEM Tools navi{{#translation:}}}}
{{FEM Tools navi{{#translation:}}}}

Latest revision as of 21:21, 26 March 2024

Other languages:
Tutoriel
Thème
Analyse par éléments finis
Niveau
Débutant
Temps d'exécution estimé
Non applicable
Auteurs
NewJoker
Version de FreeCAD
0.21 ou plus récente
Fichiers exemples
None
Voir aussi
Atelier FEM

Contexte

La préparation de la géométrie et le maillage sont des éléments cruciaux du prétraitement des simulations utilisant la méthode des éléments finis (FEM). Bien que les logiciels de simulation facilement accessibles et connectés à un environnement CAO (comme l'atelier FEM dans FreeCAD) rendent tentante la réalisation immédiate d'analyses sur de nouvelles conceptions, il est très important de se rappeler que la FEM est une méthode avancée et qu'elle nécessite une géométrie et un maillage correctement préparés afin de fournir des résultats raisonnables et précis. La vieille règle garbage in, garbage out est particulièrement importante ici. Il existe également d'autres paramètres cruciaux dont la précision de la FEM dépend fortement (tels que les propriétés des matériaux et les conditions aux limites), mais les premières étapes et certaines des sources de problèmes les plus courantes sont la préparation de la géométrie et le maillage, dont il est question sur cette page.

Types de géométrie utilisés pour la méthode des éléments finis dans FreeCAD

  • Lignes (polylignes) - utilisées pour les analyses avec des éléments de poutre
  • Surfaces - utilisées pour les analyses avec des éléments de coque
  • Solides - utilisés pour les analyses avec des éléments solides

Choix du type de géométrie

Bien que la plupart des conceptions soient constituées de solides, il est souvent fortement recommandé d'utiliser des polylignes ou des surfaces pour les éléments finis si la structure le permet :

  • Si une pièce est mince (longue et fine) et en forme de poutre et a une section transversale régulière de l'un des types de section de poutre actuellement supportés (rectangulaire, circulaire ou tube), elle doit être analysée à l'aide d'éléments de poutre (à moins qu'il n'y ait des formes spécifiques de charge, de réponse ou des détails géométriques inévitables qui invalident cette supposition). Fondamentalement, il faut dessiner une ligne centrale (des conseils sur la manière de l'extraire d'une géométrie solide existante peuvent être trouvés dans ce fil de discussion) et appliquer une section de poutre appropriée avec une éventuelle rotation. Il n'y a pas de règle unique dictant quand les éléments de poutre peuvent être utilisés, mais il est souvent conseillé que les dimensions de la section transversale soient < 1/10 de la longueur de la pièce pour que l'hypothèse de la poutre soit valide.

Pièce mince convenant à l'analyse avec des éléments de poutre - axe central mis en évidence

  • Si une pièce a des parois minces (comme les pièces en tôle), elle doit être analysée à l'aide d'éléments de coque (sauf si des résultats de contact précis sont nécessaires ou si certaines limites des éléments de coque sont rencontrées). Ce point est très important et souvent négligé. Pour obtenir des résultats précis (en particulier en cas de flexion), il faut quelques éléments (au moins 3 à 5) dans la direction de l'épaisseur. Dans le cas de pièces à parois minces, cela se traduit généralement par de grandes mailles (en particulier lorsque des tétraèdres sont utilisés, car les éléments hexaédriques ne peuvent pas être générés dans FreeCAD) et un coût de calcul élevé - puissance informatique élevée et long temps de résolution. Pour obtenir une géométrie adaptée à l'analyse avec des éléments coques, il faut dessiner une surface médiane de la pièce (des conseils sur la manière de l'extraire d'une géométrie solide existante peuvent être trouvés dans ce fil de discussion, celui-ci et celui-ci) et appliquer une épaisseur appropriée. Là encore, il n'y a pas de règle unique, mais il est généralement recommandé que l'épaisseur soit < 1/10 d'une dimension globale typique (longueur/largeur) pour que l'hypothèse de la coque soit valide.

Pièce à paroi mince adaptée à l'analyse avec des éléments de coque - mise en évidence de la surface médiane

Il convient de rappeler que les éléments de poutre et de coque utilisés dans CalculiX ne sont pas de véritables éléments de poutre/coque (ils n'utilisent pas les formulations d'éléments de poutre/coque connues dans la littérature et dans d'autres logiciels) - ils sont étendus de manière interne aux solides. Néanmoins, leur utilisation est recommandée dans les cas susmentionnés.

Validité de la géométrie

La géométrie utilisée pour la simulation par éléments finis doit être valide. Plus important encore, il ne doit pas y avoir d'intersections. Il s'agit d'un problème courant qui survient souvent lorsque des assemblages sont modélisés sans contraintes appropriées entre les pièces. L'outil Part Coupe persistante peut aider à trouver de telles interférences entre les pièces. Bien entendu, l'outil Part Union peut aider à les résoudre si elles sont intentionnelles. D'autres problèmes liés à la géométrie (tels que les géométries non pliantes, les arêtes ou faces redondantes, etc.) doivent également être résolus avant de procéder au maillage. L'outil Part Vérifier la géométrie peut être utile, mais les contrôles visuels sont également importants. Lors de la préparation d'une simulation utilisant des éléments solides et en cas de doute quant à savoir si la pièce est réellement solide ou s'il s'agit simplement d'une coque fermée, les outils susmentionnés (Part Coupe persistante et l'onglet Contenu de la forme de Part Vérifier la géométrie) peuvent permettre de clarifier la situation.

Simplification de la géométrie

Les designs préparés dans les logiciels de CAO sont généralement trop détaillés pour être adaptées aux simulations par éléments finis. Dans de nombreux cas, il est nécessaire de les simplifier au préalable. Cette étape est souvent négligée, mais elle est très importante car il peut être difficile d'obtenir un bon maillage lorsque la pièce est trop détaillée et même si un tel maillage est finalement obtenu, il peut être très dense, ce qui entraîne des temps de résolution déraisonnables. Il faut donc toujours examiner la conception et essayer de la simplifier autant que possible, en ne laissant que les caractéristiques géométriques qui peuvent avoir un impact significatif sur les résultats (résistance/rigidité) et qui ne peuvent donc pas être ignorées. Les caractéristiques suivantes sont généralement omises

  • les petits congés et chanfreins,
  • petits trous,
  • autres petits détails,
  • soudures,
  • boulons, filetages,
  • les éléments décoratifs (logos, gravures).

L'outil Part Supprimer la fonctionnalité et l'extension Defeaturing peuvent être utiles pour simplifier les pièces pour les simulations.

Géométrie originale du support

Géométrie du support simplifiée en utilisant uniquement l'outil Part Supprimer la fonctionnalité

Dans le cas des assemblages (plus d'informations à ce sujet dans l'une des sections suivantes), certaines pièces peuvent souvent être exclues des simulations et remplacées par des conditions aux limites comme si elles étaient attachées aux pièces à analyser. Cette approche est valable si les pièces exclues sont significativement plus rigides (en termes de rigidité structurelle, en tenant compte non seulement de l'élasticité du matériau mais aussi de la géométrie de la pièce) que les pièces à analyser auxquelles elles étaient connectées. En effet, les conditions aux limites fixes introduisent de la rigidité (comme si la pièce analysée était attachée à un composant infiniment rigide) et les supports flexibles comme les éléments ressorts ne sont pas disponibles dans l'atelier FEM de FreeCAD lorsque l'on utilise CalculiX (Elmer a un FEM Ressort).

La simplification de la géométrie pour la simulation par éléments finis peut également impliquer de la couper dans l'un des plans de symétrie afin d'utiliser l'hypothèse de symétrie plane dans l'analyse. Cette hypothèse n'est valable que lorsque tous les aspects suivants du modèle présentent une symétrie dans un plan donné :

  • géométrie,
  • charges,
  • conditions aux limites,
  • réponse (il faut être prudent avec les analyses de fréquence/flexion utilisant la symétrie - on n'obtiendra pas de formes de mode antisymétriques).

L'utilisation de la symétrie (1/2, 1/4 ou 1/8 du modèle) est recommandée dans la mesure du possible, car elle permet de réduire considérablement le coût de calcul de l'analyse. Un autre avantage est qu'elle élimine certains mouvements des corps rigides, ce qui facilite la contrainte de la pièce. Une condition limite de symétrie doit être appliquée aux faces appartenant au plan de coupe :

  • la translation dans la direction normale au plan de symétrie doit être bloquée pour les pièces solides,
  • la translation dans la direction normale au plan de symétrie et les rotations autres qu'autour de l'axe normal au plan de symétrie doivent être bloquées pour les pièces de coque et de poutre.

La force appliquée doit être correctement réduite si le plan de symétrie coupe la région à laquelle la force est appliquée (sans importance lorsque la charge de pression est utilisée).

Modèle de 1/8 d'un récipient sous pression cylindrique avec des conditions limites de symétrie et une charge de pression interne

Partitionnement de la géométrie

Le partitionnement est une division de la géométrie en segments plus petits. Dans d'autres logiciels, il est généralement utilisé pour permettre le maillage hexagonal, mais dans FreeCAD, il peut être utile pour d'autres raisons. La principale application du partitionnement est lorsqu'une charge (ou une condition limite) doit être appliquée uniquement à une région sélectionnée de la surface de la pièce. La façon la plus simple d'y parvenir est de créer une esquisse avec un contour approprié sur cette face et d'utiliser l'outil Part Fragments booléens pour diviser la face avec l'esquisse. Une autre raison d'utiliser le partitionnement est lorsque plusieurs matériaux doivent être appliqués à une seule pièce (sans avoir à utiliser plusieurs pièces connectées entre elles). Dans ce cas, le partitionnement peut être effectué à l'aide d'un plan et de l'outil Fragments booléens en mode Compsolid. Le partitionnement peut également être utilisé pour créer des régions pour affiner le maillage.

Pièce avec une cloison de séparation pour l'application d'une charge ou d'une condition limite

Assemblage de géométries

L'une des principales limitations actuelles de l'atelier FEM est que les maillages multiples ne sont pas pris en charge. En pratique, cela signifie qu'il n'est pas possible de mailler chaque pièce de l'assemblage individuellement, puis de connecter les pièces avec les contraintes appropriées pour l'analyse. Il est donc nécessaire de créer un objet unique contenant toutes les pièces de l'assemblage et de le mailler. Il existe plusieurs options différentes, toutes basées sur les Part Outils booléens. Le choix dépend de l'effet désiré : si chaque pièce/volume et leurs limites doivent être sélectionnables (par exemple pour l'affectation de matériaux ou les définitions des conditions limites agissant sur les faces internes) ou non :

  • Part Union : fusionne les pièces, rendant impossible leur sélection individuelle, par exemple pour les définitions de matériaux,
  • Part Composé : crée un objet composé, ce qui permet de sélectionner chaque pièce,
  • Part Connecter : fonctionne comme Part Union, mais fusionne les parties, ce qui rend impossible leur sélection individuelle,
  • Part Fragments booléens : fonctionne comme Part Composé, ce qui permet de sélectionner chaque partie.

Il est important de mentionner que si les pièces se touchent, un maillage continu sera créé sur l'objet booléen et aucune contrainte ne sera nécessaire pour la simulation. S'il y a un petit espace entre les pièces (ou une intersection à l'intérieur d'un Part Composé), le maillage ne sera pas continu et des contraintes comme la Contrainte de liaison ou la Contrainte de contact seront nécessaires. L'exécution d'une analyse en fréquence est un bon moyen de révéler si le maillage est continu ou non. Si les pièces ne sont pas connectées, les premières formes des modes avec déformation visualisées à l'aide du filtre des déformations montreront une séparation, les pièces "s'envoleront".

La forme du premier mode d'une analyse de fréquence visualisée avec le filtre des déformations - deux cubes avec un petit écart initial ont été analysés.

La sélection des régions internes (faces/volumes) peut s'avérer délicate. Elle peut être nécessaire pour l'application de différents matériaux, charges ou conditions aux limites (en particulier dans les analyses thermiques et électromagnétiques). Plusieurs méthodes sont possibles :

  • activer plan de coupe pour la sélection et le choix des faces internes,
  • masquer l'objet booléen, ne montrer qu'une des parties auxquelles il a été appliqué et la sélectionner,
  • sélectionner un autre objet externe et éditer la propriété References dans l'onglet Données d'un élément d'analyse donné (nécessite la spécification manuelle du numéro de l'objet géométrique).

Les bases du maillage

Un maillage trop grossier est l'une des sources les plus courantes d'inexactitudes et d'autres problèmes dans les simulations par éléments finis. Il s'agit souvent d'une erreur partielle des paramètres du maillage automatique. Ils génèrent généralement des maillages très grossiers et inadaptés lorsque la taille des éléments n'est pas spécifiée manuellement, mais laissée avec une valeur par défaut. Il faut toujours connaître les dimensions approximatives de la pièce, en particulier la taille de la plus petite caractéristique pertinente (l'outil Part Mesure linéaire peut aider à cet égard) et spécifier la taille maximale appropriée de l'élément sur cette base. Lors du maillage avec Gmsh, il existe également un paramètre de taille minimale des éléments qui peut empêcher la création d'éléments trop petits autour de petites caractéristiques géométriques, ce qui peut conduire à des maillages inutilement denses (et parfois même à un plantage ou à un gel de FreeCAD lorsqu'on essaie de générer de tels maillages). D'une manière générale, il est préférable de commencer par un maillage plus grossier (moins long à générer), de voir à quoi il ressemble (une certaine expérience est nécessaire) et de l'affiner si nécessaire. Il est souvent judicieux d'utiliser un maillage dense uniquement autour des zones d'intérêt (endroits présentant des gradients/concentrations de contraintes importants - encoches) et un maillage relativement grossier à l'écart de ces zones. De cette manière, le nombre d'éléments peut être considérablement réduit, ce qui permet de diminuer les temps de résolution. Le raffinement local du maillage est défini à l'aide de FEM Région de maillage FEM.

Par défaut, maillage trop grossier

Maillage amélioré globalement

Maillage amélioré localement

Le choix du type d'élément n'est pas facile et dépend de nombreux facteurs, mais la règle générale est que les éléments hexaédriques et quadruples sont préférables aux éléments tétraédriques et triangulaires. Cependant, les géométries complexes ne peuvent pas être maillées avec des éléments hexaédriques et FreeCAD ne peut pas les générer du tout (seuls les maillages quadruples peuvent être générés sur les surfaces, voir ce fil de discussion). Les éléments hexaédriques peuvent être importés à partir de maillages externes comme Gmsh et utilisés dans l'atelier FEM comme le montre cette vidéo.

Le choix de l'ordre des éléments (premier ou second) dépend des conditions d'analyse, mais dans la plupart des cas, les éléments de second ordre sont préférables. Ceci est particulièrement vrai pour les éléments triangulaires et tétraédriques, leurs versions de premier ordre (linéaires) ne sont normalement pas recommandées pour un usage régulier et ne devraient être utilisées que comme éléments de remplissage dans des régions de faible importance. Cependant, comme FreeCAD ne peut pas générer d'éléments hexaédriques, les tétraèdres linéaires peuvent être utilisés dans certains cas, si les maillages sont suffisamment denses. En particulier lors d'analyses avec une contrainte de contact.

Jacobiens négatifs

Si les règles ci-dessus sont respectées (notamment en ce qui concerne la validité de la géométrie, le defeaturing et la sélection de la taille des éléments), le maillage devrait être généré correctement. Cependant, dans certains cas, la géométrie ne peut pas être trop simplifiée, ou la procédure de modélisation est appropriée mais conduit quand même à des arêtes et des faces de petite taille. Dans ce cas, le maillage avec des éléments de second ordre peut échouer en raison de jacobiens négatifs. La raison en est que les maillages doivent suivre le modèle CAO et placer les nœuds médians des éléments de second ordre sur la géométrie. Dans le cas de formes plus complexes, les éléments peuvent être tellement étirés qu'ils deviennent inversés. Le jacobien est l'une des mesures de qualité de maillage les plus courantes. Il représente l'écart de l'élément par rapport à la forme idéale. Elle devient négative lorsque l'élément se retourne (devient inversé), soit en raison d'une déformation importante au cours de l'analyse (non prise en compte ici), soit en raison des problèmes de maillage susmentionnés. Les jacobiens négatifs dans l'atelier FEM de FreeCAD peuvent être signalés par Gmsh ou par CalculiX. Leurs emplacements dans le maillage sont mis en évidence lorsque les analyses CalculiX sont soumises en utilisant le bouton Lancer les calculs du solveur. Les conseils suivants peuvent aider à les éliminer :

  • Définir la propriété Second Order Linear de l'objet Maillage FEM à partir d'une forme avec Gmsh à true, ceci a pour résultat que les nœuds du milieu des éléments de second ordre sont simplement ajoutés au milieu des bords droits (initialement) des éléments de premier ordre, sans les attacher à la géométrie et résout le problème dans la plupart des cas,
  • Utiliser Netgen au lieu de Gmsh. Netgen est connu pour être moins sujet aux problèmes de jacobiens négatifs mais ne les signale pas non plus, de sorte que l'utilisateur ne peut s'en rendre compte qu'au moment de soumettre l'analyse,
  • Réduire davantage la taille des éléments,
  • Exporter la géométrie, essayer de la mailler dans l'interface graphique de Gmsh ou Netgen (NGSolve) ou dans un autre maillage autonome (comme Salome_Meca). Ces outils ont des fonctions supplémentaires qui peuvent aider à se débarrasser des jacobiens négatifs (par exemple, Gmsh a ce que l'on appelle des "outils d'ordre élevé"),
  • Utiliser des éléments de premier ordre. Cela ne doit être fait qu'en dernier recours car les tétraèdres de premier ordre sont connus pour leur imprécision.

Indépendamment de ces conseils, il est important de souligner une fois de plus que les jacobiens négatifs sont généralement dus à des approches de modélisation désordonnées et à un manque de préparation de la géométrie pour l'analyse (ce qui est particulièrement courant avec les modèles STEP téléchargés à partir de divers sites web). Même si le maillage est finalement généré dans de tels cas, les résultats risquent d'être de mauvaise qualité (rappelons la règle GIGO mentionnée dans le premier paragraphe). Par conséquent, le nettoyage de la géométrie et la préparation pour la FEM devraient toujours être la priorité.

Études sur la convergence des maillages

Les études de convergence de maillage sont recommandées dans tous les projets sérieux nécessitant des résultats précis. En effet, les résultats peuvent changer considérablement et s'approcher des valeurs correctes lorsque le maillage est affiné. L'approche suivante doit être utilisée :

  1. Après avoir obtenu les premiers résultats et les avoir notés (généralement la contrainte de von Mises maximale, la contrainte de von Mises à un endroit donné et le déplacement maximal), affiner le maillage (globalement ou mieux localement avec FEM Région de maillage FEM) et réexécuter la simulation.
  2. Vérifiez les résultats et notez leurs nouvelles valeurs. S'ils diffèrent significativement des résultats initiaux, affinez encore le maillage et relancez l'analyse.
  3. Répétez le processus si les résultats changent encore (généralement en augmentant) de manière significative avec l'affinement du maillage.

Il est généralement utile de créer un graphique avec un résultat donné en fonction de la densité du maillage. De cette façon, il est plus facile de remarquer quand les résultats commencent à converger. La différence acceptable dans les résultats entre deux exécutions est généralement de l'ordre de quelques pour cent (par exemple, moins de 5 %).

Dans certains cas, il peut arriver que la contrainte maximale augmente indéfiniment, quelle que soit la densité du maillage. Un tel effet non physique est connu sous le nom de singularité de contrainte. Il peut se produire pour les raisons suivantes

  • les forces concentrées appliquées aux modèles de solides et de coques,
  • conditions limites appliquées à des points (nœuds individuels),
  • angles vifs,
  • contact survenant dans un coin.

Les moyens typiques de traiter les singularités de contrainte sont les suivants :

  • appliquer des charges et des conditions aux limites à de petites zones plutôt qu'à des points - voir la section sur le partitionnement ci-dessus,
  • ajouter de petits filets aux angles vifs (une exception à la règle qui consiste à omettre les petits filets lors de la simplification de la géométrie pour les éléments finis),
  • inclure la plasticité dans le comportement du matériau pour permettre la redistribution des contraintes et limiter les contraintes aux valeurs autorisées par la définition de la plasticité tout en observant le niveau approprié de déformation (déformation plastique),
  • ignorer les singularités et lire les contraintes loin d'elles si possible (sur la base du principe de St. Venant).

Tracés de convergence de maillage typiques :
- le déplacement (courbe verte) converge rapidement,
- la contrainte maximale au niveau d'une entaille comme un trou (courbe bleue) nécessite davantage d'itérations de raffinement du maillage pour converger,
- tandis que la contrainte maximale au niveau d'un angle aigu avec une condition limite fixe (courbe rouge) ne converge pas du tout. Une singularité de contrainte se produit et un petit congé devrait être ajouté et la connexion devrait être modélisée de manière plus réaliste et plus flexible pour éviter ce comportement.