Differences

This shows you the differences between two versions of the page.

--- commit:2013:02_14 [2013/07/02 12:31] – laurent
+++ commit:2013:02_14 [2016/03/30 15:23] (current) – external edit 127.0.0.1
@@ Line 1: / Line 1: @@
+===== 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 //