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

Antonio Frasson frassoneletrica at gmail.com
Tue Mar 3 19:35:58 CET 2009


OOPS!
Sorry, I did not see my error,
thanks a lot!

Frasson

2009/3/3 David Colignon <David.Colignon at ulg.ac.be>

> 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
>>
>
> 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)={1,2};
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20090303/85ab4124/attachment.html>