[Gmsh] Cubic Mesh

Ruth V. Sabariego r.sabariego at ulg.ac.be
Sat Jun 9 08:51:20 CEST 2012


Hi Paolo,
Your surface is in the plane yz, so you should extrude it along x and not along z:

> out[] = Extrude{L, 0, 0}{Surface{700}; Layers{nb_layers}; }; 
> IndexOutlet = out[0];


HTH,
Ruth

--
Dr. Ir. Ruth V. Sabariego
University of Liege, Electrical Engineering & Computer Science, 
Applied & Computational Electromagnetics (ACE),
phone: +32-4-3663737 - fax: +32-4-3662910 - http://ace.montefiore.ulg.ac.be/



On 08 Jun 2012, at 14:24, Paolo Tricerri wrote:

> Dear Sirs,
> 
> I am trying to create a mesh of a quarter of a cube. What I need to do is to create a 2D mesh of a surface and then extrude it along the length of the cube edge.
> 
> Tha script I have written is the following.
> 
> lc = 0.1;
> R= 0.5;
> L = 2.0;
> 
> // ############# POINTS ##############
> //Base
> Point(0) = {0, R, -R};
> Point(1) = {0, -R, -R};
> Point(2) = {0, -R, R};
> Point(3) = {0, R, R};
> 
> //############### LINES ##############
> //Base`
> //Constrained face
> Line (10) = {0,1};
> Line (20) = {1,2};
> Line (30) = {2,3};
> Line (40) = {3,0};
> 
> Physical Line(500) = {10};
> Physical Line(300) = {20};// 
> 
> ############# LINE LOOPS ############
> Line Loop(1000) = {10,20,30,40}; //creation face
> 
> // ############ SURFACES ############
> 
> Plane Surface(700) = {1000}; //on line loops
> Physical Surface(400) = {700};
> 
> Transfinite Line{10,20,30,40}=5;
> Transfinite Surface { 700 } = {0,1,2,3} ; //on points
> 
> nb_layers=8;
> out[] = Extrude{0, 0, L}{Surface{700}; Layers{nb_layers}; }; 
> IndexOutlet = out[0];
> 
> ////////////////////////////////////////////////////////////////////
> 
> I get the following error:
> 
> Meshing surface 1013 (extruded)
> Error   : Uncoherent extruded quadrangle in surface 1013
> Error   : Uncoherent extruded quadrangle in surface 1013
> Error   : Uncoherent extruded quadrangle in surface 1013
> Error   : Uncoherent extruded quadrangle in surface 1013
> Error   : Uncoherent extruded quadrangle in surface 1013
> 
> ...and so on for other extruded surface.
> 
> Any suggestion about what is going wrong? The code is basically copied and pasted from the manual of gmsh so at the moment, I do not know what is going wrong.
> 
> I appreciate any suggestion I'll receive.
> 
> Thanks a lot,
> 
> Kind regards,
> 
> Paolo Tricerri
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20120609/ffd31492/attachment.html>