[Gmsh] question about volume with holes

David Colignon David.Colignon at ulg.ac.be
Tue Apr 15 16:37:08 CEST 2008


Hi,

t[] cannot be a hole into t1 because they have surfaces in common ...

Cheers,

Dave


-- 
David Colignon, Ph.D.
Collaborateur Logistique F.R.S.-FNRS (Equipements de Calcul Intensif)
ACE - Applied & Computational Electromagnetics
Institut Montefiore B28
Université de Liège
4000 Liège - BELGIQUE
Tél: +32 (0)4 366 37 32
Fax: +32 (0)4 366 29 10
WWW:    http://www.montefiore.ulg.ac.be/personnel.php?op=detail&id=898
Agenda: http://www.google.com/calendar/embed?src=david.colignon%40gmail.com



hubrys at virgilio.it wrote:
> hi!
> I'm not able to obtain a volume with holes using t[] and t1[] 
> extrusions. I mean, t[] has to makes a hole into t1[].
> Here follows my 
> project:
> 
> 
> 
> Function structur
> 
>   a1=newp; Point(a1) = {a,0,0,lc1};
>   
> a2=newp; Point(a2) = {a-r,0,0,lc1};
>   a3=newp; Point(a3) = {a-h,0,0,
> lc1};
>   a4=newp; Point(a4) = {a-r,r,0,lc1};
>   a5=newp; Point(a5) = {a-r,
> -r,0,lc1};
> 
>   b1=newp; Point(b1) = {0,b ,0,lc1};
>   b2=newp; Point(b2) = 
> {0,b-r ,0,lc1};
>   b3=newp; Point(b3) = {0,b-h ,0,lc1};
>   b4=newp; Point
> (b4) = {r,b-r ,0,lc1};
>   b5=newp; Point(b5) = {-r,b-r ,0,lc1};
> 
>   
> c1=newreg; Circle(c1) = {a1, a2, a4}; arc[i] = c1; i++;
>   c2=newreg; 
> Circle(c2) = {a4, a2, a3}; arc[i] = c2; i++;
>   c3=newreg; Circle(c3) = 
> {a3, a2, a5}; arc[i] = c3; i++;
>   c4=newreg; Circle(c4) = {a5, a2, a1}; 
> arc[i] = c4; i++;
>   c5=newreg; Circle(c5) = {b1, b2, b4}; arc[i] = c5; 
> i++;
>   c6=newreg; Circle(c6) = {b4, b2, b3}; arc[i] = c6; i++;
>   
> c7=newreg; Circle(c7) = {b3, b2, b5}; arc[i] = c7; i++;
>   c8=newreg; 
> Circle(c8) = {b5, b2, b1}; arc[i] = c8; i++;
> 
> Return
> 
> 
> 
> 
> 
>  lc=8; lc1 = 
> 1; a=27; b=a; R=50; 
>  r=1; h=2*r; sx=a; dx=-a; z=50; z1=80; 
> 
> 
> Point
> (11) = {0,0,0,lc};
> Point(21) = {-R,0,0,lc};
> Point(31) = {R,0,0,lc};
> Point(41) = {0,-R,0,lc};
> Point(51) = {0,R,0,lc};
> 
> Circle(11) = 
> {21,11,51};
> Circle(21) = {51,11,31};
> Circle(31) = {31,11,41};
> Circle
> (41) = {41,11,21};
> 
> 
> 
> i = 0;
> 
> For a In {sx:1:-3}
> 
>   Call structur;
>   b-
> =3;
> 
> EndFor
> 
> b=-1;
> 
> For a In {-1:dx:-3}
> 
>   Call structur;
>   b-=3;
> 
> EndFor
> 
> 
>  
> 
> s1 = newreg; Line Loop(s1) = {arc[]};
> s2 = newreg; Plane 
> Surface(s2) = {s1};           
> 
> s3 = newreg; Line Loop(s3) = 
> {11,21,31,41};
> s4 = newreg; Plane Surface(s4) = {s3};
> 
> 
> t[] = Extrude 
> {0,0,z}{Surface{s2};Layers{20}; Recombine;};
> t1[] = Extrude {0,0,z1}
> {Surface{s4};Layers{20}; Recombine;};
> 
> 
> 
> Thanks for the help!
> 
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh