Bonjour à tous,
Je sollicite votre aide pour trouver un code Matlab permettant de convertir des éléments triangulaires de type T3 (3 noeuds) vers T6 (6 noeuds). Malgré 2 heures de recherches sous Google en Français et en Anglais, je ne trouve pas mon bonheur.
J'utilise actuellement Distmesh mais celui-ci ne génere que des maillages avec des triangles T3.
Merci d'avance !
Bonjour,
Par coïncidence, je suis moi-même en train de regarder le code de "distmesh" pour une autre application (remaillage d'une surface à la géométrie complexe). Je suis donc en mesure de vous répondre.
"distmesh" retourne:
-) Une liste de n vertex sous forme d'une matrice p de taille (n x 3) ou (n x 2) contenant les coordonnées des vertex le long de chaque ligne.
-) Une liste de nt triangles sous forme d'une matrice t de taille (nt x 3) listant les index des vertex pour chaque triangle de type t3; chaque index pointant vers un vertex formant un sommet d'un triangle.
Si vous souhaitez des éléments t6 à partir de ces triangles t3, une manière de faire est:
-) Pour chaque triangle t3 calculer trois vertex v_sup1, v_sup2, v_sup3 supplémentaires qui seront placés aux milieux des arêtes du triangle courant. Par exemple, pour le ième triangle, ces trois vertex supplémentaires seront:
-) Ajouter (concaténer) ces nouveaux vertex à la fin de la liste des vertex p et garder trace de leurs index dans la liste. Si p possède m lignes, alors ces index valent, respectivement pour v_sup1, v_sup2 et v_sup3, m+1, m+2 et m+3.Code:v_sup1 = (p(t(i, 1),:) + p(t(i, 2), :)) / 2; v_sup2 = (p(t(i, 2),:) + p(t(i, 3), :)) / 2; v_sup3 = (p(t(i, 3),:) + p(t(i, 1), :)) / 2;
-) Ajouter ces index (m+1), (m+2) et (m+3) à la suite de la i ème ligne de la matrice t dans les colonnes respectives 4, 5 et 6. (Matrice t qui devra préalablement être redimensionnée à la taille (nt x 6) tout en conservant les données initiales dans les colonnes 1, 2 et 3).
-) Recommencer pour le triangle suivant (i+1 ème) jusqu'à avoir parcouru tous les triangles.
En espérant avoir été clair.
Petite remarque: la procédure que j'ai décrite ci-dessous dupliquera inutilement des vertex. Pour éviter cela, il faudrait maintenir une liste des arêtes pour lesquelles un vertex a été ajouté et de l'index qui lui est associé.
Bonjour et merci pour votre retour.
C'est limpide, j'espère pouvoir tester ça rapidement.
Je vous tiendrai au courant .
Merci encore !
Cette méthode fonctionne sans pb.
Plus qu'à écrire le reste du code.
Merci encore !
Re : Matlab : Conversion éléments T3 vers T6
De rien.......
