[Gmsh] Sphere problem
Markus Hengsberger
mhengs at gmx.at
Fri Oct 29 09:53:52 CEST 2004
Yes, I tried to optimize mesh, but nothing really happens.
It looks the same after the opt. process.
Markus
> mhengs at gmx.at wrote:
>
> >Hi,
> >
> >i am working with 2D structured meshed spheres. Everything works fine,
> but
> >i have problem with the small elements on the top of the spheres. My
> solver
> >cannot work with this elements.
> >
> >Maybe you have a solution how to cope with this problem.
> >
> >Markus
> >
> >Here is an exmaple of 2x2x2 sphereblock:
> >
> >/*********************************************************************
> > *
> > *
> > * source file: example.igs
> > * translater version: jIGStoGEO 0.2
> > *
> > * output files:
> > * example.geo
> > * example.msh
> > *
> > *********************************************************************/
> >
> >
> >Function Sphere
> >
> > p1 = newp; Point(p1) = {x, y, z, 1.0} ;
> > p2 = newp; Point(p2) = {x+r,y, z, 1.0} ;
> > p3 = newp; Point(p3) = {x, y+r,z, 1.0} ;
> > p4 = newp; Point(p4) = {x, y, z+r, 1.0} ;
> > p5 = newp; Point(p5) = {x-r,y, z, 1.0} ;
> > p6 = newp; Point(p6) = {x, y-r,z, 1.0} ;
> > p7 = newp; Point(p7) = {x, y, z-r, 1.0} ;
> >
> > c1 = newreg; Circle(c1) = {p2,p1,p7};
> > c2 = newreg; Circle(c2) = {p7,p1,p5};
> > c3 = newreg; Circle(c3) = {p5,p1,p4};
> > c4 = newreg; Circle(c4) = {p4,p1,p2};
> > c5 = newreg; Circle(c5) = {p2,p1,p3};
> > c6 = newreg; Circle(c6) = {p3,p1,p5};
> > c7 = newreg; Circle(c7) = {p5,p1,p6};
> > c8 = newreg; Circle(c8) = {p6,p1,p2};
> > c9 = newreg; Circle(c9) = {p7,p1,p3};
> > c10 = newreg; Circle(c10) = {p3,p1,p4};
> > c11 = newreg; Circle(c11) = {p4,p1,p6};
> > c12 = newreg; Circle(c12) = {p6,p1,p7};
> >
> >
> > l1 = newreg; Line Loop(l1) = {c5,c10,c4};
> > l2 = newreg; Line Loop(l2) = {c9,-c5,c1};
> > l3 = newreg; Line Loop(l3) = {c12,-c8,-c1};
> > l4 = newreg; Line Loop(l4) = {c8,-c4,c11};
> > l5 = newreg; Line Loop(l5) = {-c10,c6,c3};
> > l6 = newreg; Line Loop(l6) = {-c11,-c3,c7};
> > l7 = newreg; Line Loop(l7) = {-c2,-c7,-c12};
> > l8 = newreg; Line Loop(l8) = {-c6,-c9,c2};
> >
> >
> > s1 = newreg; Ruled Surface(newreg) = {l1};
> > s2 = newreg; Ruled Surface(newreg) = {l2};
> > s3 = newreg; Ruled Surface(newreg) = {l3};
> > s4 = newreg; Ruled Surface(newreg) = {l4};
> > s5 = newreg; Ruled Surface(newreg) = {l5};
> > s6 = newreg; Ruled Surface(newreg) = {l6};
> > s7 = newreg; Ruled Surface(newreg) = {l7};
> > s8 = newreg; Ruled Surface(newreg) = {l8};
> >
> > Transfinite Line {c5,c10,c4} = 15 Using Progression 1.0;
> > Transfinite Line {c9,-c5,c1} = 15 Using Progression 1.0;
> > Transfinite Line {c12,-c8,-c1} = 15 Using Progression 1.0;
> > Transfinite Line {c8,-c4,c11} = 15 Using Progression 1.0;
> > Transfinite Line {-c10,c6,c3} = 15 Using Progression 1.0;
> > Transfinite Line {-c11,-c3,c7} = 15 Using Progression 1.0;
> > Transfinite Line {-c2,-c7,-c12} = 15 Using Progression 1.0;
> > Transfinite Line {-c6,-c9,c2} = 15 Using Progression 1.0;
> >
> > Transfinite Surface {s1} = {p4,p3,p2};
> > Transfinite Surface {s2} = {p7,p3,p2};
> > Transfinite Surface {s3} = {p7,p6,p2};
> > Transfinite Surface {s4} = {p4,p6,p2};
> > Transfinite Surface {s5} = {p4,p5,p3};
> > Transfinite Surface {s6} = {p4,p6,p5};
> > Transfinite Surface {s7} = {p7,p6,p5};
> > Transfinite Surface {s8} = {p7,p3,p5};
> >
> > Recombine Surface {s1};
> > Recombine Surface {s2};
> > Recombine Surface {s3};
> > Recombine Surface {s4};
> > Recombine Surface {s5};
> > Recombine Surface {s6};
> > Recombine Surface {s7};
> > Recombine Surface {s8};
> >
> >Return
> >
> >x = 0;
> >y = 0;
> >z = 0;
> >r = 0.5;
> >
> >
> >For w In {1:2}
> >For u In {1:2}
> >For t In {1:2}
> >
> >
> > x += 1.0;
> >
> > Call Sphere ;
> >
> >
> >EndFor
> >x = 0;
> >y += 1.0;
> >
> >EndFor
> >z += 1.0;
> >x = 0;
> >y = 0;
> >
> >
> >EndFor
> >
> >
> >
>
>
>
> Have you tryed to optimize the mesh ? It takes time but
> gives very good enhancement.
>
> JFR
>
>
>
>
> --
> Prof. Jean-François Remacle
> Department of Civil Engineering
> Center for Systems Engineering and Applied Mechanics
> Université catholique de Louvain, Belgium
> Tel : +32 10 472082
> Fax : +32 10 472210
>
--
Geschenkt: 3 Monate GMX ProMail + 3 Ausgaben der TV Movie mit DVD
++++ Jetzt anmelden und testen http://www.gmx.net/de/go/mail ++++