Re: [Gmsh] bad orientation

2020-05-15 Thread Ruth Vazquez Sabariego 
Dear all,

I am not familiar with FreeFEM++...
The mesh generation happens without any issue.

If the problem is the normal orientation, maybe you could invert the normal at 
the end of the geo file, for whatever surface you want.
e.g. with the file you sent, for having outward normals in the mesh you should 
add
Reverse Surface{1,13,15,14};

Alternative to changing the loops.
Best regards,
Ruth




—
Prof. Ruth V. Sabariego
KU Leuven, Dept. Electrical Engineering ESAT/Electa, EnergyVille
http://www.esat.kuleuven.be/electa
http://www.energyville.be








On 15 May 2020, at 08:55, Jose Juan Alonso del Rosario 
mailto:josejuan.alo...@gm.uca.es>> wrote:

Dear Ugis and colleages of GMSH,

Thank you for your suggestions. I have changed the loop orientation of the 
upper and lower faces of the pentagon and the whole volumen was meshed (no hole 
at the end). In addition the same error still happens.

"read mesh ok 0Mesh3, num Tetra:= 24011, num Vertice:= 5665 num boundary 
Triangles:= 6078
Mesh3::meshS, num Triangles:= 6078, num Vertice:= 3039 num boundary Edges:= 354
read mesh ok 0Mesh3, num Tetra:= 24011, num Vertice:= 5665 num boundary 
Triangles:= 6078
Mesh3::meshS, num Triangles:= 6078, num Vertice:= 3039 num boundary Edges:= 354
 Bad orientation: The adj border element  defined by [  3 210 ]  is oriented in 
the same direction in element 3519 and in the element 367 ** bug in mesh 
construction? orientation parameter?
S, num Triangles:= 6078, num Vertice:= 3039 num boundary Edges:= 354
 Bad orientation: The adj border element  defined by [  3 210 ]  is oriented in 
the same direction in element 3519 and in the element 367 ** bug in mesh 
construction? orientation parameter? "

This is after several Optimize3D and Optimize 3d with NetGen and reading it 
with freefem++

I include the .geo code. It is a thick slice with a pentagonal hole. The gmsh 
is now 4.5.6 and the meshing options are MeshAdapt and Frontal. I have played 
with the signs of plane surfaces 14 and 15.

Thanks in advance!!!

Jose





El jue., 14 may. 2020 a las 13:54, Ugis Lacis 
(mailto:ugis.la...@gmail.com>>) escribió:
Dear Jose,

I think I have had similar error, when importing GMSH mesh into FreeFEM++, few 
years ago (although it seems that the error message might have changed). It 
turned out that the direction of loop definition in GMSH was the issue for me. 
See below my exchange with FreeFEM++ mailing list.

Maybe it can be helpful for you.

Best regards,
Ugis Lacis
Researcher
KTH Mechanics

--
 Forwarded Message 
Subject:Re: [Freefempp] Issue with loading mesh from GMSH.
Date:   Fri, 3 Aug 2018 09:41:29 +0200
From:   Uģis Lācis 
To: FreeFem mailing list 



Hi everyone,

I managed to solve the issue myself. It turned out that the direction of GMSH 
loop was surviving till FreeFEM++ and while GMSH did not complain about 
anything, FreeFEM++ couldn't use the mesh. mesh2 could be fixed by changing 
signs in all lines within the line loop, i.e.,
change
Line Loop(5) = {18, 19, 20, 12, 7, 8};
to
Line Loop(5) = {-18, -19, -20, -12, -7, -8};

Maybe this helps to someone at some point.

Best,
Ugis

On 02/08/18 15:30, Uģis Lācis wrote:
Dear FreeFEM++ users and developers,

I have run into an issue which I am unable to solve by myself. I have been 
using FreeFEM++ together with GMSH without any problems for quite some time. 
Now, however, trying to load a 2D mesh generated by GMSH, gives me following 
error:

Assertion fail : (area>0)
line :281, in file ./include/fem.hpp

I attach the corresponding geo and msh files of two meshes (mesh1 - working 
fine, mesh2 - does not load), as well as edp file for reading and plotting the 
mesh.

I have carefully checked that physical tags are applied correctly, I have tried 
exporting both in MSH and MESH (medit) formats, using GMSH versions 3.06 and 
2.16.0, as well as using FreeFEM++ versions 3.56 and 3.610002, all try has 
led me to the same behaviour. I have identified that the issue comes from the 
surface mesh, if I export only line mesh (comment out "Physical Surface(11) = 
{3,4};" from mesh2.geo and run "gmsh -2 mesh2.geo"), then it loads fine. But 
checking the mesh within GMSH, everything seems to be perfectly fine...

Do you have any suggestions on how to solve this issue?

Best regards,
Ugis Lacis
Researcher
KTH Mechanics
--


On 14/05/2020 12:03, Jose Juan Alonso del Rosario wrote:
Dear colleages

I am trying to build the mesh for a slice with a pentagonal hole.
GMSH runs properly with good quality tets but when trying to read it from 
freefem++ the error is:

" Bad

Re: [Gmsh] R: R: R: ellipsoid issue

2019-09-06 Thread Ruth Vazquez Sabariego 
No issues with my up to date and freshly compiled version (4.5.0-git-3284474ec) 
under MacOS 10.14.6.
My local version of OCC seems to be older 7.1.0.

Ruth

—
Prof. Ruth V. Sabariego
KU Leuven, Dept. Electrical Engineering ESAT/Electa, EnergyVille
http://www.esat.kuleuven.be/electa
http://www.energyville.be








On 6 Sep 2019, at 09:05, Alessandro Vicini 
mailto:alessandro.vic...@sitael.com>> wrote:

Well, it works under Linux though (with no warning messages). I guess the 
problem might be different OCC versions used by Gmsh under different OS (I see 
version 7.3.0 under Linux, and 7.3.1 under Win).

A.



Da: Peter Johnston 
mailto:p.johns...@griffith.edu.au>>
Inviato: venerdì 6 settembre 2019 08:54
A: Jeremy Theler mailto:jer...@seamplex.com>>; 
g...@geuz.org; Alessandro Vicini 
mailto:alessandro.vic...@sitael.com>>
Oggetto: Re: R: [Gmsh] R: ellipsoid issue

Thanks Alessandro,

I just ran mine again and I get the same behaviour as previously, with the same 
error messages.

I’m not sure what these means, but perhaps someone else can contribute?

Thanks again,

Peter.

-

Associate Professor Peter Johnston (FAustMS, FIMA)
School of Environment and Science
Griffith University | Nathan | QLD 4111 | Technology (N44) Room 3.19
T +61 7 373 57748| F +61 7 373 57656 Email 
p.johns...@griffith.edu.au
On 6 Sep 2019, 4:47 PM +1000, Alessandro Vicini 
mailto:alessandro.vic...@sitael.com>>, wrote:

Hello Peter,

in my case Gmsh just gets stuck when trying to read the geo file. No messages 
or whatsoever, I need to kill the program.
I just downloaded 4.4.1, but I get the same behavior.

A.


Da: Peter Johnston 
mailto:p.johns...@griffith.edu.au>>
Inviato: venerdì 6 settembre 2019 04:08
A: Jeremy Theler mailto:jer...@seamplex.com>>; Alessandro 
Vicini mailto:alessandro.vic...@sitael.com>>; 
g...@geuz.org
Oggetto: Re: [Gmsh] R: ellipsoid issue

Dear Jeremy and Alessandro,

Thanks for your input. Well, the simple ellipsoid problem seemed to work for me 
on my mac this morning. It guess it wouldn’t work yesterday after a long 
afternoon of frustration.

Anyway, I have a more complex example (see attached) that is causing some 
strange behaviour. When I run the script I get the following warning (see the 
first log file as well):

Warning : More than two degenerated edges in one model face of an OCC model

I am not sure what is causing this or how to get rid of it. Also, I am not sure 
if it is a problem or not. I can create the 1D mesh without any problem. 
However, when creating the 2D mesh, the following lines appear (see full log 
file for complete print out):

Info: :-( There are 2 intersections in the 1D mesh (curves 4 12)
Info: 8-| Splitting those edges and trying again
Info: :-) All edges recovered after 1 iteration

I am not sure if this is a problem or not. Any thoughts? At this point, I have 
not completed the final geometry that I want and am concerned these warnings 
might lead to other problems further down the track.

Any advice would be welcome.

Thanks again,

Peter.

-

Associate Professor Peter Johnston (FAustMS, FIMA)
School of Environment and Science
Griffith University | Nathan | QLD 4111 | Technology (N44) Room 3.19
T +61 7 373 57748| F +61 7 373 57656 Email 
p.johns...@griffith.edu.au
On 5 Sep 2019, 11:38 PM +1000, Peter Johnston 
mailto:p.johns...@griffith.edu.au>>, wrote:
Strangely, it works under 4.4.1 on my linux machine. I'll check the mac version 
again tomorrow.

Thanks,

Peter.

Associate Professor Peter Johnston (FAustMS, FIMA)
School of Environment and Science
Griffith University | Nathan | QLD 4111 | Technology (N44) Room 3.19
T +61 7 373 57748| F +61 7 373 57656 Email 
p.johns...@griffith.edu.au

From: Peter Johnston
Sent: Thursday, 5 September 2019 8:10 PM
To: Jeremy Theler mailto:jer...@seamplex.com>>; Alessandro 
Vicini mailto:alessandro.vic...@sitael.com>>; 
g...@geuz.org mailto:g...@geuz.org>>
Subject: Re: [Gmsh] R: ellipsoid issue

Thanks guys,

Yes it worked for me last week under 4.3.0. Today I downloaded 4.4.1 and it no 
longer works. I am using the Mac version.

Regards,

Peter.

Peter Johnston, Griffith University.

From: Jeremy Theler mailto:jer...@seamplex.com>>
Sent: Thursday, September 5, 2019 8:02:01 PM
To: Alessandro Vicini 
mailto:alessandro.vic...@sitael.com>>; Peter 
Johnston mailto:p.johns...@griffith.edu.au>>; 
g...@geuz.org mailto:g...@geuz.org>>
Subject: Re: [Gmsh] R: ellipsoid issue

Indeed, see attached png.
I am using latest Git linked against OCE 7.3.0 in GNU/Linux.

On Thu, 2019-09-05 at 08:13 +, Alessandro Vicini wrote:

Hello Peter,



I’m using 4.4.0

Re: [Gmsh] Gmsh with petsc/slepc 3.10?

2019-02-27 Thread Ruth Vazquez Sabariego 
Hi Christophe,

My Gmsh version (from git, up to date) is compiled (today!) with
petsc-3.10.1 and slepc-3.10.0

without any problem
petsc and slepc (both compiled in October 2018)

Best,
Ruth



—
Prof. Ruth V. Sabariego
KU Leuven
Dept. Electrical Engineering ESAT/Electa, EnergyVille
http://www.esat.kuleuven.be/electa
http://www.energyville.be

Free software: http://gmsh.info | http://getdp.info | http://onelab.info







On 27 Feb 2019, at 13:42, Christophe Trophime 
mailto:christophe.troph...@lncmi.cnrs.fr>> 
wrote:

Hi,
did someone tried to rebuild Gmsh with petsc/slepc 3.10?
Is it working?

Best
C



Christophe TROPHIME
Research Engineer

CNRS - LNCMI
25, rue des Martyrs
BP 166
38042 GRENOBLE Cedex 9
FRANCE

Tel : +33 (0)4 76 88 90 02
Fax : +33 (0) 4 76 88 10 01
Office U 19
M@il : 
christophe.troph...@lncmi.cnrs.fr
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

Re: [Gmsh] Question about Progression and Bump

2018-10-24 Thread Ruth Vazquez Sabariego 
Dear Nathan,

As you have noticed, the Bump does not do a double geometrical progression.
The formula used is a bit more complicated, you can see what it actually does 
in the code (meshGEdge.cpp).

A double progression is not (yet?) available.
Work around by dividing your second region in two?

Regards,
Ruth


—
Prof. Ruth V. Sabariego
KU Leuven
Dept. Electrical Engineering ESAT/Electa, EnergyVille
http://www.esat.kuleuven.be/electa
http://www.energyville.be

Free software: http://gmsh.info | http://getdp.info | http://onelab.info







On 23 Oct 2018, at 22:07, Nathan J. Neeteson 
mailto:nneete...@rglinc.com>> wrote:

Hello,

I have a question about using Progression and Bump in neighboring blocks in a 
block-structured mesh. What I want is for the two blocks to have cells of the 
same height where they meet. I know the thickness of the first cell from the 
wall (dx0), the length of each line (L), and the progression I want (r), so 
when I’m using Progression I can calculate the number of nodes to use as N = 
log(1+(L*(r-1))/(dx0)) / log(r) – 1. Then I can set the appropriate lines to 
Transfinite and assign them N Using Progression r and it works as expected.

However, in one of my blocks I want to refine towards both the top and bottom, 
so I need to use the Bump option. The problem is that I have no idea what “r” 
value to use after Bump (using the same r value as is used for Progression does 
not give me the results I want) and I also have no idea how to calculate the 
number of points to use along this line. MY first instinct was to calculate the 
number of points needed as double the number of points for a single progression 
with half of the length (two progressions end to end). I think this is right, 
but when I use r after Bump I get no refinement. What value am I supposed to 
use for number of points and the rate of growth when using the Bump function so 
that I get the equivalent of the product of two Progressions end to end?

Here is a zoomed in look at where I want the first cells to match heights:

The bottom block uses single progression and the upper block you are seeing the 
lower end of the double progression.


Here is my .geo file:
---
// orifice properties
orificeLength = 0.03;
orificeRadius = 0.002;

// pipe properties
pipeRadius = 0.05;

// distance from orifice to inlet
inletDist = 5*pipeRadius;

// distance from orifice to outlet
outletDist = 5*pipeRadius;

// grid size parameters
r = 1.1; // growth parameter
dx0 = 1*10^(-4); // first wall dist for y+=1

upstreamNx = (Log(1 + ((inletDist*(r-1))/(dx0))) / Log(r)) - 1; // number of 
elements in x direction upstream
downstreamNx = (Log(1 + ((outletDist*(r-1))/(dx0))) / Log(r)) - 1; // number of 
elements in x direction upstream
orificeNx = orificeLength/dx0; // number of elements in x direction inside 
orifice

upstreamOrificeNy = orificeRadius/dx0;
upstreamPipeNy = 2*((Log(1 + pipeRadius-orificeRadius)/2)*(r-1))/(dx0))) / 
Log(r)) - 1);


// element index
i = 1;

// ~~-~~ POINTS ~~-~~
// order of point definition doesn't matter, just give descriptive names

inletCenterPoint = i;
Point(i) = {-orificeLength-inletDist,0,0,1};
i=i+1;

inletPipeWallPoint = i;
Point(i) = {-orificeLength-inletDist,pipeRadius,0,1};
i=i+1;

inletOrificeWallPoint = i;
Point(i) = {-orificeLength-inletDist,orificeRadius,0,1};
i=i+1;

upOrificeCenterPoint = i;
Point(i) = {-orificeLength,0,0,1};
i=i+1;

upOrificePipeWallPoint = i;
Point(i) = {-orificeLength,pipeRadius,0,1};
i=i+1;

upOrificeOrificeWallPoint = i;
Point(i) = {-orificeLength,orificeRadius,0,1};
i=i+1;


// ~~-~~ LINES ~~-~~
// all lines should be defined going up and/or to the right for consistency

inletOrificeLine = i;
Line(i) = {inletCenterPoint,inletOrificeWallPoint};
i=i+1;

inletPipeLine = i;
Line(i) = {inletOrificeWallPoint,inletPipeWallPoint};
i=i+1;

upstreamPipeWallLine = i;
Line(i) = {inletPipeWallPoint,upOrificePipeWallPoint};
i=i+1;

upstreamCenterLine = i;
Line(i) = {inletCenterPoint,upOrificeCenterPoint};
i=i+1;

upstreamOrificeLine = i;
Line(i) = {inletOrificeWallPoint,upOrificeOrificeWallPoint};
i=i+1;

upstreamOrificePlateLine = i;
Line(i) = {upOrificeOrificeWallPoint,upOrificePipeWallPoint};
i=i+1;

upstreamOrificeEntryLine = i;
Line(i) = {upOrificeCenterPoint,upOrificeOrificeWallPoint};
i=i+1;

// ~~-~~ LINE LOOPS ~~-~~
// all line loops should be oriented clockwise

upstreamPipeLoop = i;
Line Loop(i) = 
{inletPipeLine,upstreamPipeWallLine,-upstreamOrificePlateLine,-upstreamOrificeLine};
i=i+1;

upstreamOrificeLoop = i;
Line Loop(i) = 
{inletOrificeLine,upstreamOrificeLine,-upstreamOrificeEntryLine,-upstreamCenterLine};
i=i+1;



// ~~-~~ PLANE SURFACES ~~-~~

upstreamPipePlane = i;
Plane

Re: [Gmsh] command line options

2018-09-21 Thread Ruth Vazquez Sabariego 
You can get all the available options by typing in the command line:
gmsh -—help

2D mesh => -2
3D mesh => -3
gmsh -smooth number_steps
gmsh -optimize


Best,
Ruth


On 21 Sep 2018, at 14:47, Alessandro Vicini 
mailto:alessandro.vic...@sitael.com>> wrote:


I usually work interactively, and I go through the following steps to generate 
a 3D mesh:

1)  2D mesh;
2)  Smoothing of the 2D mesh (repeated a few times);
3)  3D mesh;
4)  Optimization of 3D mesh.

How would I carry out the same process from command line? Thank you.

A.


___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

Re: [Gmsh] Create a 3D cube with each face having a physical surface

2018-02-09 Thread Ruth Vazquez Sabariego 
Using your second option, you may use BooleanFragments.
Try at the end of your file:

// Surface{7,8,9,10,11,12} In Volume{1}; // Comment out this line
Physical Volume(100)  = {1};

aux() = BooleanFragments{ Volume{1}; Delete; }{ Surface{7:12}; Delete;};
Printf("", aux()); // For checking the numbers of the generated geometrical 
entities

// For visibility
Recursive Color Red{Volume{1};}
Recursive Color Yellow{Surface{7:12};}

With regard to your first option, you can only create rectangles in the XY 
plane, but you can translate and rotate to place them where you want.

HTH
Ruth

—
Prof. Ruth V. Sabariego
KU Leuven
Dept. Electrical Engineering ESAT/Electa, EnergyVille
http://www.esat.kuleuven.be/electa
http://www.energyville.be

Free software: http://gmsh.info | http://getdp.info | http://onelab.info







On 8 Feb 2018, at 22:13, Sathyanarayan Rao 
> wrote:

Dear GMSH team,

I am trying to generate a cube with each face having a rectangular surface that 
will be later identified as a sensor.

1) In case I want to use rectangle command of opencascade is it possible to 
create rectangle other than in XY plane ?



```
SetFactory("OpenCASCADE");
Mesh.Algorithm3D = 4; // 3D mesh algorithm (1=Delaunay, 4=Frontal, 5=Frontal 
Delaunay, 6=Frontal Hex, 7=MMG3D, 9=R-tree)
R = 0.05;
w = R;
Box(1) = {-R,-R,-R, 2*R,2*R,2*R};
Rectangle(10) = {-R/2,-R/2,-R, R,R};
Rectangle(11) = {-R/2,-R/2,R, R,R};
```

 2)  When I manually create surfaces in each face of the cube, they come as 
separate regions in 3D mesh.

```
SetFactory("OpenCASCADE");
Mesh.Algorithm3D = 4; // 3D mesh algorithm (1=Delaunay, 4=Frontal, 5=Frontal 
Delaunay, 6=Frontal Hex, 7=MMG3D, 9=R-tree)
R = 0.05;
R1 = R;
w = R;
Box(1) = {-R,-R,-R, 2*R,2*R,2*R};
q = 0.005;

/* Right electrode */
RE = 300;

Point(RE+1) = { R1,  R1/2,  R1/2, q};
Point(RE+2) = { R1, -R1/2,  R1/2, q};
Point(RE+3) = { R1, -R1/2, -R1/2, q};
Point(RE+4) = { R1,  R1/2, -R1/2, q};
Line(RE+13) = {RE+4, RE+1};
Line(RE+14) = {RE+1, RE+2};
Line(RE+15) = {RE+2, RE+3};
Line(RE+16) = {RE+3, RE+4};
Line Loop(RE+7) = {RE+16, RE+13, RE+14, RE+15};
Plane Surface(7) = {RE+7};
Physical Surface(1) = {7};

/* Left electrode */

RE = 400;

Point(RE+1) = { -R1,  R1/2,  R1/2, q};
Point(RE+2) = { -R1, -R1/2,  R1/2, q};
Point(RE+3) = { -R1, -R1/2, -R1/2, q};
Point(RE+4) = { -R1,  R1/2, -R1/2, q};
Line(RE+13) = {RE+4, RE+1};
Line(RE+14) = {RE+1, RE+2};
Line(RE+15) = {RE+2, RE+3};
Line(RE+16) = {RE+3, RE+4};
Line Loop(RE+7) = {RE+16, RE+13, RE+14, RE+15};
Plane Surface(8) = {RE+7};
Physical Surface(10001) = {8};


// create plate electrodes
/* back electrode */

RE = 500;

Point(RE+1) = {  R1/2,  -R1, R1/2, q};
Point(RE+2) = { -R1/2,   -R1, R1/2, q};
Point(RE+3) = {  -R1/2,   -R1,-R1/2, q};
Point(RE+4) = { R1/2,  -R1, -R1/2, q};
Line(RE+13) = {RE+4, RE+1};
Line(RE+14) = {RE+1, RE+2};
Line(RE+15) = {RE+2, RE+3};
Line(RE+16) = {RE+3, RE+4};
Line Loop(RE+7) = {RE+16, RE+13, RE+14, RE+15};
Plane Surface(9) = {RE+7};
Physical Surface(10002) = {9};


/* front electrode */

RE = 600;

Point(RE+1) = {  R1/2,  R1, R1/2, q};
Point(RE+2) = { -R1/2,  R1, R1/2, q};
Point(RE+3) = {  -R1/2,  R1,-R1/2, q};
Point(RE+4) = { R1/2, R1, -R1/2, q};
Line(RE+13) = {RE+4, RE+1};
Line(RE+14) = {RE+1, RE+2};
Line(RE+15) = {RE+2, RE+3};
Line(RE+16) = {RE+3, RE+4};
Line Loop(RE+7) = {RE+16, RE+13, RE+14, RE+15};
Plane Surface(10) = {RE+7};
Physical Surface(10003) = {10};

/* top electrode */

RE = 700;

Point(RE+1) = {  R1/2,R1/2,R1, q};
Point(RE+2) = { -R1/2,R1/2, R1,q};
Point(RE+3) = {  -R1/2, -R1/2,R1, q};
Point(RE+4) = { R1/2,  -R1/2,R1, q};
Line(RE+13) = {RE+4, RE+1};
Line(RE+14) = {RE+1, RE+2};
Line(RE+15) = {RE+2, RE+3};
Line(RE+16) = {RE+3, RE+4};
Line Loop(RE+7) = {RE+16, RE+13, RE+14, RE+15};
Plane Surface(11) = {RE+7};
Physical Surface(10004) = {11};


/* bottom electrode */

RE = 800;

Point(RE+1) = {  R1/2,R1/2,-R1, q};
Point(RE+2) = { -R1/2,R1/2, -R1,q};
Point(RE+3) = {  -R1/2, -R1/2,-R1, q};
Point(RE+4) = { R1/2,  -R1/2,-R1, q};
Line(RE+13) = {RE+4, RE+1};
Line(RE+14) = {RE+1, RE+2};
Line(RE+15) = {RE+2, RE+3};
Line(RE+16) = {RE+3, RE+4};
Line Loop(RE+7) = {RE+16, RE+13, RE+14, RE+15};
Plane Surface(12) = {RE+7};
Physical Surface(10005) = {12};

Surface{7,8,9,10,11,12} In Volume{1};
Physical Volume(100)  = {1};
```

Somehow the command "Surface{} In Volume" seems to not care for what I want.  
Please let me know
how can I get a cube with each face having a rectangular surface with its own 
physical number.






Best Regards,

Sathyanarayan Rao, PhD student
Earth and Life Institute/Environmental Sciences (ELI-e)
Université catholique de Louvain
c.037, Croix du Sud 2, Louvain-la-Neuve, Belgium
Phone: 010473827 ( intercom 73827)



___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

Re: [Gmsh] Missing surfaces in mesh when using Boolean operators

2018-02-07 Thread Ruth Vazquez Sabariego 
Hi Isabela,

You may use BooleanFragments instead of BooleanDifference. That seems to work:
Try the following at the end of your file:

/*
voxcyl() = BooleanDifference{Volume{v};}{Volume{union19()}; Delete;};
//Physical Volume(0) = {voxcyl()};
*/

// --- Definition of Final Volume 
voxcyl() += BooleanFragments{ Volume{v}; Delete;}{ Volume{union19()}; Delete; };
Physical Volume(0) = {voxcyl({#voxcyl()-1})}; // around the vessels
Physical Volume(1) = {voxcyl({0:#voxcyl()-2})}; // the vessels

Printf("",voxcyl()); // Just for checking

#voxcyl() gives you the number of elements in the array.

Best regards,
Ruth

—
Prof. Ruth V. Sabariego
KU Leuven
Dept. Electrical Engineering ESAT/Electa, EnergyVille
http://www.esat.kuleuven.be/electa
http://www.energyville.be

Free software: http://gmsh.info | http://getdp.info | http://onelab.info







On 7 Feb 2018, at 11:50, ISABELA PAREDES CISNEROS 
> wrote:

Dear all,

I am trying to mesh a cube with cylindrical holes in it. Some of the holes 
intersect. I am using a BooleanUnion of all the cylinders and then a 
BooleanDifference of the cube and the cylinders. In the meshing, I get no error 
but the mesh on the cylinders boundary is wrong. A lot of surfaces are missing.

Do you have any idea of what am I doing wrong?

Best regards,

Isabela

P.D: I have commented some lines (10 and 113) in the GEO file to visualize the 
mesh over the cylinders.

___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

Re: [Gmsh] GMSH Linux32bit updates

2017-11-21 Thread Ruth Vazquez Sabariego 
Hi!
You can get official gmsh binary files for linux at:
http://gmsh.info/bin/Linux/

Regards,
Ruth Sabariego
—
Prof. Ruth V. Sabariego
KU Leuven
Dept. Electrical Engineering ESAT/Electa, EnergyVille
http://www.esat.kuleuven.be/electa
http://www.energyville.be

Free software: http://gmsh.info | http://getdp.info | http://onelab.info







On 21 Nov 2017, at 15:16, La alll 
> wrote:

Hello!
I had troubles to compile GMSH under linux32 bit (opencascade is not linking)
Is it possible to make official binary file?
Linux32 is very popular under virtualbox, especially for education
Regards
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

Re: [Gmsh] Structured 2D mesh for structure more than 8 corners !

2017-02-16 Thread Ruth Vazquez Sabariego 
Transfinite surfaces have to be defined with only 4 points.
The number of divisions have to be adapted accordingly.

Try:
Transfinite Line {2}=5;
Transfinite Line {6}=6;
Transfinite Line {8}= (5+6+6*2+2)-5+1;
Transfinite Line {7}= 11;
Transfinite Line {1}= 11;
Transfinite Line {3,5} =6;
Transfinite Line {4} = 2;

//Transfinite Surface {10}={5,6,7,8,1,2,3,4};
Transfinite Surface {10}={1,2,7,8};
Recombine Surface {10};


Mesh.Smoothing = 5; // For a nicer mesh ;-)

HTH
Ruth






—
Prof. Ruth V. Sabariego
KU Leuven
Dept. Electrical Engineering ESAT/Electa, EnergyVille
http://www.esat.kuleuven.be/electa
http://www.energyville.be

Free software: http://gmsh.info | http://getdp.info | http://onelab.info







On 15 Feb 2017, at 14:41, Gladiator Gladiator 
> wrote:

Dear all,

I am trying to make a structured mesh for following algorithm, but I could not 
make it.
And my mesh is still unstructured !
Would you please let me know where I made a mistake ?
Many thanks,

this is my algorithm:
//**
Point(1) = {0, 0, 0, 1.0};
Point(2) = {10, 0, 0, 1.0};
Point(3) = {10, 4, 0, 1.0};
Point(4) = {5, 4, 0, 1.0};
Point(5) = {5, 5, 0, 1.0};
Point(6) = {10, 5, 0, 1.0};
Point(7) = {10, 10, 0, 1.0};
Point(8) = {0, 10, 0, 1.0};


Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 5};
Line(5) = {5, 6};
Line(6) = {6, 7};
Line(7) = {7, 8};
Line(8) = {8, 1};

Line Loop(9) = {7, 8, 1, 2, 3, 4, 5, 6};

Plane Surface(10) = {9};
Physical Surface(11) = {10};

Transfinite Line {2}=5;
Transfinite Line {6}=6;
Transfinite Line {8}= 11;
Transfinite Line {7}= 11;
Transfinite Line {1}= 11;
Transfinite Line {3,5} =6;
Transfinite Line {4} = 2;

Transfinite Surface {10}={5,6,7,8,1,2,3,4};
Recombine Surface {10};


Again, thank you all,
Best,

Gladi, PhD student
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

Re: [Gmsh] Newbie question

2016-05-26 Thread Ruth Vazquez Sabariego 
Gmsh is the mesher and a possible interface to other FEA software.

It is fully interfaced with the Fe solftware Getdp in Onelab
Maybe you should start having a look at the models in:
http://onelab.info/wiki/ONELAB


Regards,
Ruth

—
Prof. Ruth V. Sabariego
KU Leuven
Dept. Electrical Engineering ESAT/Electa, EnergyVille
http://www.esat.kuleuven.be/electa
http://www.energyville.be

Free software: http://gmsh.info | http://getdp.info | http://onelab.info







On 26 May 2016, at 02:12, Karen Pease > 
wrote:

Hi all.  I've spent all evening trying out different meshing and FEA
software and gotten precisely nowhere... so I decided it's time to
ask. Where is a good place to start?  All I'm wanting to do is:

1) Take an existing model (say, exported as STL)
2) Import it
3) Assign material properties, anchor points, and assign forces
4) Run a FEA simulation
5) See the stresses and deflections

Seems like a pretty basic workflow, but I can hardly get past the
first step.  GMSH got me the furthest - at least I could actually see
it.  But it's just a wireframe, I can't do anything more with it than
pan around it, there's no interactive points.  The model is just
window decor as far as I can tell.  I see nothing in the visibility
options that makes it any more interactive.

(At least this is further than I got with other tools... netgen might
as well have had the interface in Swahili).

I've tried reading the documentation, but what's called a "tutorial"
(http://gmsh.info/doc/texinfo/gmsh.html#Tutorial) is practically an
API manual.  I've tried watching videos on Youtube, but one, not only
do they rarely deal with importing geometry like STLs, but for most
programs the interface in the videos looks different from what I see
because I use Linux and they use Windows.  And the Windows versions
usually seem to have more user friendly versions (at least gmsh is
better than calculix in this regard, with calculix it's like half the
interface is missing)

What do I do?  How can I get my bearings here?

Any feedback at all would be greatly appreciated... thanks, all.

- kv, Karen

___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

Re: [Gmsh] Gmsh difficulty - Socket Listening issue

2016-05-26 Thread Ruth Vazquez Sabariego 
It looks as if the gmsh is not able to find the solver, by default getdp…
I would just remove the defined Solver (right click for getting the option).
When opening inductor.pro again, it will prompt for getting it back.

HTH, otherwise give us a bit more information e.g. operative system and 
versions you are using.

Regards,
Ruth


—
Prof. Ruth V. Sabariego
KU Leuven
Dept. Electrical Engineering ESAT/Electa, EnergyVille
http://www.esat.kuleuven.be/electa
http://www.energyville.be

Free software: http://gmsh.info | http://getdp.info | http://onelab.info







On 26 May 2016, at 05:40, Daniel West 
> wrote:

Hello,

Just while trying to view an example (http://onelab.info/wiki/Inductor), it 
keeps getting cut short by an error: "Abnormal server termination: http://127.0.0.1:0/>>"

I'm really not sure why this is happening; I have tried restoring default 
settings to no avail.

Thanks!
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh