[Gmsh] Transfinite Volume command - Example of use -

[Ruth V. Sabariego]

Hi Paolo,
First, you do not need to define physical quantities, points, surfaces, volumes, for generating the geometry. 
Those numbers are useful in e.g. a FE software. You can thus remove them.
For getting the volume, you need a surface loop, and the numbers in the loop are the surface numbers, not point numbers.
Change the end of your file to:

Surface Loop(100) = {2, 3, 1, 6, 4, 5};
Volume(1) = {100};
Transfinite Volume(1) ;
Physical Volume(1000) = {1};

Regards,
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 23 Oct 2012, at 13:56, Paolo Tricerri wrote:

> Dear all,
> 
> I am having problem at using the command Transfinite Volume. I am
> trying to give to the definition of the volume the list of the eight
> points at the corners of the geometry but the volume does not get
> meshed.
> 
> Hereafter you can see the code I have using. Btw, I would like to ask
> if any of you could give an example of use of the command Transfinite
> Volume {expression-list} | "*" < = {expression-list} since I do not
> understand if I have to pass to is Physical Point or Points.
> 
> Many thanks,
> 
> I wish you all a nice day,
> 
> Paolo Tricerri
> 
> Code:
> 
> lc = 0.1;
> H= 0.44;
> L = 0.48;
> pH = 0.16;
> thickness = pH/10;
> 
> // ############# POINTS ##############
> //Base
> Point(0) = {0, 0, 0, lc};
> Point(1) = {0, H, 0, lc};
> Point(2) = {L, H+pH, 0, lc};
> Point(3) = {L, H, 0, lc};
> 
> Point(4) = {0, 0, thickness, lc};
> Point(5) = {0, H, thickness, lc};
> Point(6) = {L, H+pH, thickness, lc};
> Point(7) = {L, H, thickness, lc};
> 
> Physical Point(1000) = {0};
> Physical Point(2000) = {1};
> Physical Point(3000) = {2};
> Physical Point(4000) = {3};
> Physical Point(5000) = {4};
> Physical Point(6000) = {5};
> Physical Point(7000) = {6};
> Physical Point(8000) = {7};
> 
> 
> //############### LINES ##############
> //Base`
> Line (1) = {0,1};
> Line (2) = {1,2};
> Line (3) = {2,3};
> Line (4) = {3,0};
> 
> Line (5) = {4,5};
> Line (6) = {5,6};
> Line (7) = {6,7};
> Line (8) = {7,4};
> 
> Line (9) = {0,4};
> Line (10) = {1,5};
> Line (11) = {3,7};
> Line (12) = {2,6};
> 
> // ############# LINE LOOPS ############
> Line Loop(1) = {1,2,3,4}; //on lines
> Line Loop(2) = {5,6,7,8}; //on lines
> Line Loop(3) = {1,10,-5,-9}; //on lines
> Line Loop(4) = {3,11,-7,-12}; //on lines
> Line Loop(5) = {9,-8,-11,4}; //on lines
> Line Loop(6) = {10,6,-12,-2}; //on lines
> 
> // ############ SURFACES ############
> 
> Plane Surface(1) = {1}; //on line loops
> Plane Surface(2) = {2}; //on line loops
> Plane Surface(3) = {3}; //on line loops
> Plane Surface(4) = {4}; //on line loops
> Plane Surface(5) = {5}; //on line loops
> Plane Surface(6) = {6}; //on line loops
> 
> // ############ TRANSFINITE LINES ############
> Transfinite Line{1,3,5,7}=23; //Number of points on a 1/4 of circle
> Transfinite Line{2,4,6,8}=15;       //Number of structure layers
> Transfinite Line{9,10,11,12}=4;
> // ############ TRANSFINITE SURFACES ############
> 
> Transfinite Surface { 1 } = { 0,1,2,3} ; //on points
> Transfinite Surface { 2 } = { 4,5,6,7} ; //on points
> Transfinite Surface { 3 } = { 0,4,5,1} ; //on points
> Transfinite Surface { 4 } = { 3,7,6,2} ; //on points
> Transfinite Surface { 5 } = { 0,4,7,3} ; //on points
> Transfinite Surface { 6 } = { 1,5,6,2} ; //on points
> Physical Surface(20)={1};
> Physical Surface(30)={2};
> Physical Surface(40)={3};
> Physical Surface(50)={4};
> Physical Surface(60)={5};
> Physical Surface(70)={6};
> 
> //Surface Loop(100) = {0,1,2,3,4,5,6,7};
> 
> Volume(1) = {0,1,2,3,4,5,6,7};
> //Transfinite Volume(1) = {0,1,2,3,4,5,6,7};
> Physical Volume(1000) = {1};
> //############# EXTRUSION TO 3D ###############
> 
> _______________________________________________
> 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/20121023/c60003b6/attachment.html>