===== Commit 2013-02-14 =====

===== Modifs =====

------------
Ajout d'un [[doc:user:elements:volumes:ruptureCritere#MaximumPrincipalStrainRuptureCriterion|critère de rupture sur les déformations principales]] adapté à la modélisation de la fracture osseuse.

{{:commit:2013:rupt_strain0000.png?500|}}

------------
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} \|} \|${{:commit:2013:triangle.png?200|}}



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 [[doc:user:results:courbes_res#extracteur_du_nombre_d_elements_actifsinactifs_lies_a_une_interaction|extracteurs sur les interactions]] permettant de compter le nombres d'éléments actifs et inactifs, i.e. rompus ou non. 
  
===== Fichiers ajoutés/supprimés =====
<code> 
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
</code>
===== Tests ajoutés/supprimés =====
<code>
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
</code>

 --- //[[cedric.laurent@ulg.ac.be|Cédric Laurent]] 2013/02/15 //