CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Open Source Meshers: Gmsh, Netgen, CGNS, ... (http://www.cfd-online.com/Forums/openfoam-meshing-open/)
-   -   Error : Something wrong in edge loop of size=2, no sign ! (http://www.cfd-online.com/Forums/openfoam-meshing-open/87966-error-something-wrong-edge-loop-size-2-no-sign.html)

Intern May 4, 2011 06:09

Error : Something wrong in edge loop of size=2, no sign !
 
Hi ;
I try to modelize different layers with a hole inside 2 of them but I have this error message and I don't know how to deal with...

My code is :

// Gmsh project created on Mon Mar 28 16:28:26 2011
// all dimensions in mm
lc = 2; //0.005;

// Definition dimensions //

// PVC1
wPVC1=20;
hPVC1=20;
tPVC1=0.035;
zPVC1 =0;
lcPVC1=0.2;

//PVC2
wPVC2=20;
hPVC2=20;
tPVC2=0.1;
zPVC2=zPVC1+tPVC1;
lcPVC2=0.2;

//Alu
wAlu = 20;
hAlu = 20;
tAlu = 0.2;
lcAlu= 0.2;
zAlu=zPVC2+tPVC2;

// Definition structure //

//PVC1 (layer 3)
p1 = newp; Point(p1) = {0, 0, zPVC1, lcPVC1};
p2 = newp; Point(p2) = {wPVC1, 0, zPVC1, lcPVC1};
p3 = newp; Point(p3) = {wPVC1, hPVC1, zPVC1, lcPVC1};
p4 = newp; Point(p4) = {0, hPVC1, zPVC1, lcPVC1};
l1 = newc; Line(l1)={p1, p2};
l2 = newc; Line(l2)={p2, p3};
l3 = newc; Line(l3)={p3, p4};
l4 = newc; Line(l4)={p4, p1};
m=18;
Transfinite Line{l1,l2,l3,l4} = m Using Bump 10.0;
ll1 = newreg; Line Loop(ll1) = {l1, l2, l3, l4};

//Volume 3
S1 = newreg; Plane Surface(S1) = {ll1} ;
VolPVC1[]=Extrude {0, 0, tPVC1} {
Surface{S1}; Layers{ {5,5}, {0.5,1} }; //{5,5}=number of layers in the fist X part where X=first value in second brackets
};
Physical Volume("PVC1") = {VolPVC1[1]};


//PVC2 (layer 2)
p5 = newp; Point(p5) = {0, 0, zPVC2, lcPVC2};
p6 = newp; Point(p6) = {wPVC2, 0, zPVC2, lcPVC2};
p7 = newp; Point(p7) = {wPVC2, hPVC2, zPVC2, lcPVC2};
p8 = newp; Point(p8) = {0, hPVC2, zPVC2, lcPVC2};
l5 = newc; Line(l5)={p5, p6};
l6 = newc; Line(l6)={p6, p7};
l7 = newc; Line(l7)={p7, p8};
l8 = newc; Line(l8)={p8, p5};
Transfinite Line{l5,l6,l7,l8} = m Using Bump 10.0;
ll2 = newreg; Line Loop(ll2) = {l5, l6, l7, l8};

//hole
p37 = newp; Point(p37) = {20.15, 20, zPVC2, lcPVC2};
p38 = newp; Point(p38) = {19, 19.85, zPVC2, lcPVC2};
p39 = newp; Point(p39) = {19.85, 20, zPVC2, lcPVC2};
p40 = newp; Point(p40) = {20, 20.15, zPVC2, lcPVC2};
l40 = newc ; Spline(l40) ={p37,p38,p39,p40,p37};
Transfinite Line{l40} = 34;
ll15 = newreg ; Line Loop(ll15)={l40};

//Volume 2
S2 = newreg; Plane Surface(S2) = {ll2, ll15} ;
VolPVC2[]=Extrude {0, 0, tPVC2} {
Surface{S2}; Layers{ {5,5}, {0.5,1} }; //{5,5}=number of layers in the fist X part where X=first value in second brackets
};
Physical Volume("PVC2") = {VolPVC2[1]};


//Alu (layer 1)
p9 = newp; Point(p9) = {0, 0, zAlu, lcAlu};
p10 = newp; Point(p10) = {wAlu, 0, zAlu, lcAlu};
p11 = newp; Point(p11) = {wAlu, hAlu, zAlu, lcAlu};
p12 = newp; Point(p12) = {0, hAlu, zAlu, lcAlu};
l10 = newc; Line(l10)={p9, p10};
l11 = newc; Line(l11)={p10, p11};
l12 = newc; Line(l12)={p11, p12};
l13 = newc; Line(l13)={p12, p9};
Transfinite Line{l10,l11,l12,l13} = m Using Bump 10.0;

//hole
p41 = newp; Point(p41) = {20.15, 20, zAlu, lcAlu};
p42 = newp; Point(p42) = {20, 19.85, zAlu, lcAlu};
p43 = newp; Point(p43) = {19.85, 20, zAlu, lcAlu};
p44 = newp; Point(p44) = {20, 20.15, zAlu, lcAlu};
l41 = newc ; Spline(l41) ={p41,p42,p43,p44,p41};
Transfinite Line{l41} = 34;
ll16 = newreg ; Line Loop(ll16)={l41};

//Volume
ll3 = newreg; Line Loop(ll3) = {l10, l11, l12, l13};
S3 = newreg; Plane Surface(S3) = {ll3, ll16} ;
VolAlu[]=Extrude {0, 0, tAlu} {
Surface{S3}; Layers{ {5,5}, {0.5,1} };
};
Physical Volume("Alu") = {VolAlu[1]};


//Volume 4 (hole)

S12 = newreg; Plane Surface(S12) = {ll15} ;
VolSi2[]=Extrude {0, 0, tPVC2+tAlu} {
Surface{S12}; Layers{ {5,5}, {0.5,1} };
};
Physical Volume("Si2") = {VolSi2[1]} ;


If anyone can help me...Because I looked on forum but didn't find anything.

Thank you very much in advance.

Antoine


All times are GMT -4. The time now is 18:28.