[Gmsh] Problens in a second order 3D mesh.

David Colignon David.Colignon at ulg.ac.be
Tue Mar 3 16:52:43 CET 2009


Hi Antonio,

it seems there are some problems in your geometry:

a) Line(2) = {2, 2};

b)
Surface Loop(60) = {51, 12, 47, 56, 55, 21};
Volume(61) = {60};
Surface Loop(62) = {21, 33, 10, 25, 29, 34};
Volume(63) = {62};
Physical Volume(64)={61,63};

The extrude command "Extrude {0, 0, -1} {Surface{10, 12};}" automagically creates the volumes n°1 and 2 (To know the
volume numbers, check Tools->Visibility

see aaaa.geo

Cheers,

Dave


-- 
David Colignon, Ph.D.
Collaborateur Logistique du F.R.S.-FNRS
CÉCI - Consortium des Équipements de Calcul Intensif
ACE - Applied & Computational Electromagnetics
Sart-Tilman B28
Université de Liège
4000 Liège - BELGIQUE
Tél: +32 (0)4 366 37 32
Fax: +32 (0)4 366 29 10
WWW:    http://hpc.montefiore.ulg.ac.be/
Agenda: http://www.google.com/calendar/embed?src=david.colignon%40gmail.com



Antonio Frasson wrote:
> Firstly I would like to thank you for the excellent software and for 
> yours attention with all that still have doubts.
> 
> When I ask a mesh of second order for the geometry described below, the 
> gmsh freezes. I have a 2.2.6 source  version compiled here andI test a 
> 3.0 binary version. I am using a UBUNTU 8.10 LINUX distribution.
> 
> Thanks,
> Antonio Frasson
> 
> 
> Point(1) = {-0.5, -0.25, 1};
> Point(2) = {0.5, -0.25, 1};
> Point(3) = {0.5, 0.25, 1};
> Point(4) = {-0.5, 0.25, 1};
> Point(5) = {0.5, 0.0, 1};
> Point(6) = {-0.5, 0.0, 1};
> 
> Line(1) = {1, 2};
> Line(2) = {2, 2};
> Line(3) = {2, 5};
> Line(4) = {5, 3};
> Line(5) = {3, 4};
> Line(6) = {4, 6};
> Line(7) = {6, 1};
> Line(8) = {6, 5};
> Line Loop(9) = {8, -3, -1, -7};
> Plane Surface(10) = {9};
> Line Loop(11) = {8, 4, 5, 6};
> Plane Surface(12) = {11};
> Extrude {0, 0, -1} {
>   Surface{10, 12};
> }
> Physical Surface(57) = {56, 34};
> Physical Surface(58) = {21};
> Physical Surface(59) = {55, 12, 51, 47, 25, 29, 10, 33};
> Surface Loop(60) = {51, 12, 47, 56, 55, 21};
> Volume(61) = {60};
> Surface Loop(62) = {21, 33, 10, 25, 29, 34};
> Volume(63) = {62};
> Physical Volume(64)={61,63};
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 2009/3/3 Antonio Frasson <frassoneletrica at gmail.com 
> <mailto:frassoneletrica at gmail.com>>
> 
>     Firstly I would like to thank you for the excellent software and for
>     yours attention with all that still have doubts.
> 
>     When I ask a mesh of second order for the geometry described below,
>     the gmsh freezes. I have a 2.2.6 source  version compiled here andI
>     test a 3.0 binary version. I am using a UBUNTU 8.10 LINUX distribution.
> 
>     Thanks,
>     Antonio Frasson
> 
> 
>     2009/3/2 Jean-Francois Remacle <jean-francois.remacle at uclouvain.be
>     <mailto:jean-francois.remacle at uclouvain.be>>
> 
> 
>         Le 02-mars-09 à 20:36, antonio at ele.ufes.br
>         <mailto:antonio at ele.ufes.br> a écrit :
> 
> 
>             Why the gmsh stop to work if I ask to a second order 3D
>             tetrahedron mesh?
> 
> 
>         could you be more precise ?
> 
> 
>             Thanks,
>             Antonio Frasson
> 
>             _______________________________________________
>             gmsh mailing list
>             gmsh at geuz.org <mailto:gmsh at geuz.org>
>             http://www.geuz.org/mailman/listinfo/gmsh
> 
> 
>         ----------------
>         Prof. Jean-Francois Remacle
>         Universite catholique de Louvain (UCL)
>         Tel : +32-10-472082 -- Mobile : +32-473-909930
> 
> 
> 
> 
> 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: aaaa.geo
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20090303/cb99b322/attachment.geo>