Ajout d'un critère de rupture sur les déformations principales adapté à la modélisation de la fracture osseuse.
Ajout de deux nouveaux descripteurs de la qualité du maillage, qui viennent remplacer le précédant portant sur les angles que forment les différentes faces/côtés des éléments. Les nouveaux descripteurs fonctionnent et possèdent la même définition sur tous les types d'éléments en 2D et 3D.
- Le premier concerne un aspect ratio
de l'élément, c'est-à-dire le ratio entre le plus grand et le plus petit segment d'un élément (code IF_MESH_QUALITY_AR
). L'aspect ratio n'est pas normé: un élémént “idéal” possède un aspect ratio égal à 1, et il peut aller théoriquement jusque l'infini.
- Le second appelé squish index
est tiré de la doc Ansys/Fluent et est défini à partir du produit scalaire entre le vecteur (normalisé) allant du centre de l'élément au centre de chaque côté (2D) ou face (3D) et le vecteur normal (normalisé) au côté / à la face correspondante (code IF_MESH_QUALITY
).
$ squish index= max_i \| 1-\frac{\overrightarrow{n_i}}{\|\overrightarrow{n_i}\|}.\frac{\overrightarrow{C-c_i}}{\| \overrightarrow{C-c_i} \|} \|$
avec $ C $ le centre de l'élément, $ c_i $ le centre de la face ou du côté $ i $, et $ \overrightarrow{n_i} $ la normale à la face ou au côté $ i $. Les “pires” éléments possèdent un squish_index égal à 1, tandis qu'un élément “idéal” possède un squish index nul.
Note: le calcul du squish index ne fonctionne pas pour les quadratic_tetra, car la surface associée à une face triangulaire quadratique n'est pas (encore) définie proprement.
Ajout de deux extracteurs sur les interactions permettant de compter le nombres d'éléments actifs et inactifs, i.e. rompus ou non.
A oo_meta\mtElements\rupture\MaximumPrincipalStrainRuptureCriterion.cpp A oo_meta\mtElements\rupture\MaximumPrincipalStrainRuptureCriterion.h A oo_meta\mtFEM\extractors\NumberOfActiveElementsExtractor.cpp A oo_meta\mtFEM\extractors\NumberOfActiveElementsExtractor.h A oo_meta\mtFEM\extractors\NumberOfInactiveElementsExtractor.cpp A oo_meta\mtFEM\extractors\NumberOfInactiveElementsExtractor.h
A oo_meta\apps\biomec\bending_rupture.py A oo_meta\apps\biomec\bending_rupture_sample.py A oo_meta\apps\monosMeca\quality_hexa.py A oo_meta\apps\monosMeca\quality_penta.py A oo_meta\apps\monosMeca\quality_quad.py A oo_meta\apps\monosMeca\quality_quadratic_hexa.py A oo_meta\apps\monosMeca\quality_quadratic_quad.py A oo_meta\apps\monosMeca\quality_quadratic_tetra.py A oo_meta\apps\monosMeca\quality_quadratic_triangle.py A oo_meta\apps\monosMeca\quality_tetra.py A oo_meta\apps\monosMeca\quality_triangle.py
— Cédric Laurent 2013/02/15