[Gmsh] (no subject)

Mihai Victor Pricop vpricop at aero.incas.ro
Sat Aug 4 21:35:13 CEST 2001


 Hello!
 I'm a beginer in CFD domain and several weeks ago, searching for free 
software for mesh generation, I was finding GMSH.
 I developed some modest examples, which are attached to this message.
 I ask you for the permission to use those examples for a presentation which 
I will send for studies applications.
 I mention that GMSH is the only software which I use for 3D mesh generation.

 Best regards
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wing1.geo
Type: text/x-c++
Size: 7361 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20010804/7dd8171c/attachment.geo>
-------------- next part --------------
r=1;
R=8*r;
l=30*r;
lr=r/15;
lR=R/15;
p1=newp; Point(p1)={0,0,0,lr};
p2=newp; Point(p2)={0,0,r,lr};
p3=newp; Point(p3)={-r,0,0,lr};
p4=newp; Point(p4)={0,0,-r,lr};
p5=newp; Point(p5)={0,-r,0,lr};
p6=newp; Point(p6)={0,r,0,lr};
p17=newp; Point(p17)={r,0,0,lr};

p7=newp; Point(p7)={0,0,R,lR};
p8=newp; Point(p8)={-R,0,0,lR};
p9=newp; Point(p9)={0,0,-R,lR};
p10=newp; Point(p10)={0,-R,0,lR};
p11=newp; Point(p11)={0,R,0,lR};
p12=newp; Point(p12)={l,0,0,lR};

lR=2.5*lR;
p13=newp; Point(p13)={l,0,R,lR};
p14=newp; Point(p14)={l,-R,0,lR};
p15=newp; Point(p15)={l,0,-R,lR};
p16=newp; Point(p16)={l,R,0,lR};

c1=newreg; Circle(c1)={p7,p1,p8};
c2=newreg; Circle(c2)={p8,p1,p9};
c3=newreg; Circle(c3)={p7,p1,p10};
c4=newreg; Circle(c4)={p10,p1,p9};
c5=newreg; Circle(c5)={p9,p1,p11};
c6=newreg; Circle(c6)={p11,p1,p7};
c7=newreg; Circle(c7)={p13,p12,p14};
c8=newreg; Circle(c8)={p14,p12,p15};
c9=newreg; Circle(c9)={p15,p12,p16};
c10=newreg; Circle(c10)={p16,p12,p13};
c11=newreg; Circle(c11)={p8,p1,p10};
c12=newreg; Circle(c12)={p8,p1,p11};

// sfera mica
c13=newreg; Circle(c13)={p2,p1,p5};
c14=newreg; Circle(c14)={p5,p1,p4};
c15=newreg; Circle(c15)={p4,p1,p6};
c16=newreg; Circle(c16)={p6,p1,p2};
c17=newreg; Circle(c17)={p2,p1,p3};
c18=newreg; Circle(c18)={p3,p1,p4};
c19=newreg; Circle(c19)={p4,p1,p17};
c20=newreg; Circle(c20)={p17,p1,p2};
c21=newreg; Circle(c21)={p3,p1,p5};
c22=newreg; Circle(c22)={p5,p1,p17};
c23=newreg; Circle(c23)={p17,p1,p6};
c24=newreg; Circle(c24)={p6,p1,p3};

l1=newreg; Line(l1)={p7,p13};
l2=newreg; Line(l2)={p10,p14};
l3=newreg; Line(l3)={p9,p15};
l4=newreg; Line(l4)={p11,p16};
//l5=newreg; Line(l5)={p8,p1};
//l6=newreg; Line(l6)={p8,p1};
//l7=newreg; Line(l7)={p10,p1};
//l8=newreg; Line(l8)={p9,p1};
//l9=newreg; Line(l9)={p11,p1};
//l10=newreg; Line(l10)={p14,p12};
//l11=newreg; Line(l11)={p13,p12};
//l12=newreg; Line(l12)={p16,p12};
//l13=newreg; Line(l13)={p15,p12};

ls1=newreg; Line Loop(ls1)={c1,c11,-c3};
s1=newreg; Ruled Surface(s1)={ls1};
Physical Surface(5)={s1};

ls2=newreg; Line Loop(ls2)={c6,c1,c12};
s2=newreg; Ruled Surface(s2)={ls2};
Physical Surface(5)={s2};

ls3=newreg; Line Loop(ls3)={c2,-c4,-c11};
s3=newreg; Ruled Surface(s3)={ls3};
Physical Surface(5)={s3};

ls4=newreg; Line Loop(ls4)={c2,c5,-c12};
s4=newreg; Ruled Surface(s4)={ls4};
Physical Surface(5)={s4};

ls5=newreg; Line Loop(ls5)={c3,l2,-c7,-l1};
s5=newreg; Ruled Surface(s5)={ls5};
Physical Surface(4)={s5};

ls6=newreg; Line Loop(ls6)={c4,l3,-c8,-l2};
s6=newreg; Ruled Surface(s6)={ls6};
Physical Surface(4)={s6};

ls7=newreg; Line Loop(ls7)={c5,l4,-c9,-l3};
s7=newreg; Ruled Surface(s7)={ls7};
Physical Surface(4)={s7};

ls8=newreg; Line Loop(ls8)={c6,l1,-c10,-l4};
s8=newreg; Ruled Surface(s8)={ls8};
Physical Surface(4)={s8};

ls9=newreg; Line Loop(ls9)={c7,c8,c9,c10};
s9=newreg; Plane Surface(s9)={ls9};
Physical Surface(4)={s9};

// Sfera mica
//Partea frontala
ls10=newreg; Line Loop(ls10)={c17,c21,-c13};
s10=newreg; Ruled Surface(s10)={ls10};
Physical Surface(3)={s10};

ls11=newreg; Line Loop(ls11)={c21,c14,-c18};
s11=newreg; Ruled Surface(s11)={ls11};
Physical Surface(3)={s11};

ls12=newreg; Line Loop(ls12)={c15,c24,c18};
s12=newreg; Ruled Surface(s12)={ls12};
Physical Surface(3)={s12};

ls13=newreg; Line Loop(ls13)={c16,c17,-c24};
s13=newreg; Ruled Surface(s13)={ls13};
Physical Surface(3)={s13};

// Partea posterioara
ls14=newreg; Line Loop(ls14)={c13,c22,c20};
s14=newreg; Ruled Surface(s14)={ls14};
Physical Surface(3)={s14};

ls15=newreg; Line Loop(ls15)={c14,c19,-c22};
s15=newreg; Ruled Surface(s15)={ls15};
Physical Surface(3)={s15};

ls16=newreg; Line Loop(ls16)={c19,c23,-c15};
s16=newreg; Ruled Surface(s16)={ls16};
Physical Surface(3)={s16};

ls17=newreg; Line Loop(ls17)={c23,c16,-c20};
s17=newreg; Ruled Surface(s17)={ls17};
Physical Surface(3)={s17};

lv1=newreg; Surface Loop(lv1)={s1,s2,s3,s4,s5,s6,s7,s8,s9,s10,s11,s12,s13,s14,s15,s16,s17};
v1=newreg; Volume(v1)={lv1};
Physical Volume(0)={v1};