[Kratos] Nelson: some idea

Enrique Escolano escolano en cimne.upc.edu
Jue Jun 20 12:52:23 CEST 2013


Nelson, si quieres puedes hablar en español, llevas demasiado tiempo fuera...

Para detectar el caso de aristas que conectan dos mallas de tetraedros que solo se tocan en dichas aristas:
1- extrae los triangulos de contorno de los tetraedros 
(enumeras todas las caras de todos los tetrahedros y eliminas de la lista los que se repitan dos veces, considerando caras iguales si tiene los mismos tres nodos en cualquier orden)
2- enumeras todas las aristas de estos triángulos de contorno, y eliminas de la lista los que se repitan exactamente dos veces (considerando aristas iguales si tienen los mismos dos nodos en cualquier orden)
y lo que te quedarán son las aristas raras, que se repetirán cuatro veces (dos por cada 'dominio')

Para detectar el caso de nodos que conectan dos mallas de tetraedros que solo se tocan en dichos nodos se me ocurre una solución algo mas complicada:
1- extrae los triangulos de contorno de los tetraedros 
2- 'Colorea' (marca con un número) las distintas regiones que puedan existir, para ello escoges un triangulo que no tenga 'color' asignado, le das el siguente número de color libre, y propagas este color a todos los triángulos adyacentes por arista, y haces que se repita lo mismo para cada uno de los tres vecinos (que no tengan ya un color asignado).
3- para cada triangulosmarcas sus tres nodos con el color del elemento, y si el nodo ya tenía asignado un color distinto es que es uno de los nodos especiales que buscabas.

Enrique
  ----- Original Message ----- 
  From: Nelson Lafontaine 
  To: Kratos maining list ; Enrique Escolano 
  Sent: Thursday, June 20, 2013 12:12 PM
  Subject: Nelson: some idea


  Hi friends


  In this case i need an implementation idea related we geometry.

  Consider a two volume domain discretized in tetrahedra element.

  both domains is connected  with only one node or one side on the

  tettahedra (a line). How can i identify  this nodes?? there a numerical

  tricks to do do it??? The main objective of this is to force the 3D discrete

  fracture.


  Any ida is welcome. Thanks in advance.


  Nelson 



  -- 
  Ing. Nelson M. Lafontaine

------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: http://listas.cimne.upc.edu/pipermail/kratos/attachments/20130620/37eeb2ea/attachment.htm 


Más información sobre la lista de distribución Kratos