Metafor

ULiege - Aerospace & Mechanical Engineering

User Tools

Site Tools


commit:2013:02_14

Commit 2013-02-14

Modifs


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.

Fichiers ajoutés/supprimés

 
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

Tests ajoutés/supprimés

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

commit/2013/02_14.txt · Last modified: 2016/03/30 15:23 (external edit)