Metafor

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