[Gmsh] overlapping lines between two surfaces

ragab.fayez at gmail.com ragab.fayez at gmail.com
Wed Jul 10 12:08:01 CEST 2013


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 <javascript:>
> http://www.geuz.org/mailman/listinfo/gmsh
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20130710/5022513c/attachment.html>