[Gmsh] overlapping lines between two surfaces

[Jose Javier Muñoz Criollo]

Hi

I think so. At the beginning I wrote a function in deal.ii to "manually"
delete the repeated points. Later I found that I can use a command in Gmsh
to do the same.
The procedure I followed was:

1 - write the geometry script into .geo file
2 - import .geo file into Gmsh. Generate mesh (2D or 3D). Save mesh into
.mesh file.
3 - write another .geo file with the commands

> Merge "filename.msh";
> Coherence Mesh;
>
4 - import last .geo file into gmsh. In the message console you should read
the amount of repeated points deleted.
5 - Save last mesh into another .mesh file.

I tried this with a much more complex 3D geometry and it worked quite well.
Although I didn't tried with the original example I posted.

By the way, if you're using Gmsh with Deal.ii don't forget to check the
last tutorial example posted "step-49". It deals with meshes in general and
specifically it offers some hints about the interaction between Gmsh and
Deal.ii.

Hope this helps.

Regards
Javier


2013/7/10 <ragab.fayez at gmail.com>

> Hi Javier
> Have you got a solution for this problem because i faced the same problem
>
> Thanks
>
> On Saturday, 29 October 2011 00:55:33 UTC+2, Javier Muñoz wrote:
>
>> Hi
>>
>> I'm trying to extrude two lines:
>>
>> lc = 0.01;
>> p1 = newp; Point (p1) = {0,0,0,lc};
>> p2 = newp; Point (p2) = {0.5,0,0,lc};
>> p3 = newp; Point (p3) = {1.0,0,0,lc};
>>
>> l1 = newl; Line (l1) = {p1,p2};
>> l2 = newl; Line (l2) = {p2,p3};
>>
>> p1 -------------- p2 ---------------p3
>>
>>
>> Extrusion height is 1 for line 1 (with 1 layer). Line two has a height
>> of z = 0.5 (with 1 layer) and then I want to extrude the 'top' face
>> generated to a height of z = 1 (also with 1 layer).
>>
>>
>> out1[] = Extrude {0,1,0} {
>>                 Line {1};
>>                 Layers {{1},{1}};
>>                 Recombine;
>>                 };
>> out2[] = Extrude {0,0.5,0} {
>>                 Line {2};
>>                 Layers {{1},{1}};
>>                 Recombine;
>>                 };
>> out3[] = Extrude {0,0.5,0} {
>>                 Line {out2[0]};
>>                 Layers {{1},{1}};
>>                 Recombine;
>>                 };
>> Coherence;
>>
>> |---------------||------------**--|
>> |                ||           **     |
>> |                ||           **     |
>> |                ||           **     |
>> |                ||           **     |
>> |                ||           **     |
>> |                ||           **     |
>> |                ||-----------**---|
>> |                ||           **     |
>> |                ||           **     |
>> |                ||           **     |
>> |                ||           **     |
>> |                ||           **     |
>> |                ||           **     |
>> p1 ------------ p2 ------------ p3
>>
>> The problem is: I have lines overlapping in the interface of surface
>> out1. The line on the right side of out1 overlaps with the lines on the
>> left side of out2 and out3. Is there anyway to avoid this? As you can
>> see I already try to use Coherence (maybe in the wrong way?).
>>
>> I want to import the mesh generated into deal.ii but as long as I have
>> overlapping lines, deal.ii appears to recognize the interfaces as
>> boundaries.
>>
>> I appreciate your help.
>>
>> Regards
>> Javier Munoz
>>
>>
>> ______________________________**_________________
>> gmsh mailing list
>> gm... at geuz.org
>> http://www.geuz.org/mailman/**listinfo/gmsh<http://www.geuz.org/mailman/listinfo/gmsh>
>>
>>


-- 
José Javier Muñoz Criollo
PhD Student | Cardiff University School of Engineering
E-mail: josejmcriollo at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20130710/971542ce/attachment.html>