[Gmsh] Transfinite Volume command - Example of use -

Paolo Tricerri paolotricerri at gmail.com
Tue Oct 23 23:48:15 CEST 2012


Dear Ruth,

I really thank you for your help. It worked!
The physical quantities are created because actually I have to use the
mesh in FE code.

Thanks again,

Paolo

2012/10/23 Ruth V. Sabariego <r.sabariego at ulg.ac.be>:
> 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
>
>