Bonjour,
J'essaie de générer le maillage d'une géométrie comportant un cercle que
je devrais raffiner impérativement du centre vers le bord pour créer une
sorte de couche limite à ce niveau. Le cercle ne constitue pas un trou
et le nombre de points sur son bord est déjà imposé par l'anneau autour.
même en utilisant un raffinement sur les diagonales sécantes, je
n'arrive pas à obtenir quelque chose de propre. Auriez vous une
suggestion (.geo joint)?
Bien cordialement,
Ouafa
//---------------------------------------------------------------------------
// centre
Point(1)={0, 0, 0, lc};
// cercle intérieur
Point(6)={-d2, 0, d2, lc};
Point(7)={d2, 0, d2, lc};
Point(8)={d2, 0, -d2, lc};
Point(9)={-d2, 0, -d2, lc};
// cercle extérieur
Point(10)={-d3, 0, d3, lc};
Point(11)={d3, 0, d3, lc};
Point(12)={d3, 0, -d3, lc};
Point(13)={-d3, 0, -d3, lc};
// hublot
Point(14)={-3*d3, 0, d4, lc};
Point(15)={-d3, 0, d4, lc};
Point(16)={d3, 0, d4, lc};
Point(17)={3*d3, 0, d4, lc};
Point(18)={3*d3, 0, d3, lc};
Point(19)={3*d3, 0, -d3, lc};
Point(20)={3*d3, 0, -d4, lc};
Point(21)={d3, 0, -d4, lc};
Point(22)={-d3, 0, -d4, lc};
Point(23)={-3*d3, 0, -d4, lc};
Point(24)={-3*d3, 0, -d3, lc};
Point(25)={-3*d3, 0, d3, lc};
Point(26)={L2, 0, d4, lc};
Point(27)={L2, 0, d3, lc};
Point(28)={L2, 0, -d3, lc};
Point(29)={L2, 0, -d4, lc};
Point(30)={-L1, 0, d4, lc};
Point(31)={-L1, 0, d3, lc};
Point(32)={-L1, 0, -d3, lc};
Point(33)={-L1, 0, -d4, lc};
// aval
Point(34)={L3, 0, d4, lc};
Point(35)={L3, 0, d3, lc};
Point(36)={L3, 0, -d3, lc};
Point(37)={L3, 0, -d4, lc};
//---------------------------------------------------------------------------
// cercle intérieur
Circle(5) = {6,1,7};
Circle(6) = {7,1,8};
Circle(7) = {8,1,9};
Circle(8) = {9,1,6};
// carré extérieur
Line(9 ) = {14,15};
Line(10) = {15,16};
Line(11) = {16,17};
Line(12) = {17,18};
Line(13) = {18,19};
Line(14) = {19,20};
Line(15) = {20,21};
Line(16) = {21,22};
Line(17) = {22,23};
Line(18) = {23,24};
Line(19) = {24,25};
Line(20) = {25,14};
Line(32) = {21,12};
Line(33) = {24,13};
Line(34) = {13,22};
Line(35) = {15,10};
Line(36) = {10,25};
Line(37) = {18,11};
Line(38) = {11,16};
// aval hublot
Line(21) = {17,26};
Line(22) = {26,27};
Line(23) = {27,28};
Line(24) = {28,29};
Line(25) = {29,20};
Line(41) = {27,18};
Line(42) = {28,19};
Line(31) = {12,19};
// amont
Line(26) = {23,33};
Line(27) = {33,32};
Line(28) = {32,31};
Line(29) = {31,30};
Line(30) = {30,14};
Line(39) = {25,31};
Line(40) = {24,32};
// aval
Line(43) = {26,34};
Line(44) = {27,35};
Line(45) = {28,36};
Line(46) = {29,37};
Line(47) = {34,35};
Line(48) = {35,36};
Line(49) = {36,37};
// diagonales du cercle intérieur
//Line(50) = {1,6};
//Line(51) = {1,7};
//Line(52) = {1,8};
//Line(53) = {1,9};
// diagonales du cercle extérieur
Line(54) = {6,10};
Line(55) = {7,11};
Line(56) = {8,12};
Line(57) = {9,13};
// cercle extérieur
Circle(58) = {10,1,11};
Circle(59) = {10,1,13};
Circle(60) = {13,1,12};
Circle(61) = {12,1,11};
// boucles de lignes
Line Loop(1) = {5,6,7,8};
//Line Loop(2) = {-51,50,5};
//Line Loop(3) = {6,-52,51};
//Line Loop(4) = {52,7,-53};
//Line Loop(5) = {53,8,-50};
Line Loop(6) = {-55,-5,54,58};
Line Loop(7) = {-61,-56,-6,55};
Line Loop(8) = {56,-60,-57,-7};
Line Loop(9) = {-8,57,-59,-54};
Line Loop(10) = {13,-31,61,-37};
Line Loop(11) = {-32,16,-34,60};
Line Loop(12) = {59,-33,19,-36};
Line Loop(13) = {-38,-58,-35,10};
Line Loop(14) = {12,37,38,11};
Line Loop(15) = {14,15,32,31};
Line Loop(16) = {34,17,18,33};
Line Loop(17) = {35,36,20,9};
Line Loop(18) = {22,41,-12,21};
Line Loop(19) = {23,42,-13,-41};
Line Loop(20) = {24,25,-14,-42};
Line Loop(21) = {-20,39,29,30};
Line Loop(22) = {-19,40,28,-39};
Line Loop(23) = {-18,26,27,-40};
Line Loop(24) = {47,-44,-22,43};
Line Loop(25) = {48,-45,-23,44};
Line Loop(26) = {49,-46,-24,45};
// surfaces
s1 = newreg; Ruled Surface(s1) = {1};
//s2 = newreg; Plane Surface(s2) = {2};
//s3 = newreg; Plane Surface(s3) = {3};
//s4 = newreg; Plane Surface(s4) = {4};
//s5 = newreg; Plane Surface(s5) = {5};
s6 = newreg; Plane Surface(s6) = {6};
s7 = newreg; Plane Surface(s7) = {7};
s8 = newreg; Plane Surface(s8) = {8};
s9 = newreg; Plane Surface(s9) = {9};
s10 = newreg; Plane Surface(s10) = {10};
s11 = newreg; Plane Surface(s11) = {11};
s12 = newreg; Plane Surface(s12) = {12};
s13 = newreg; Plane Surface(s13) = {13};
s14 = newreg; Plane Surface(s14) = {14};
s15 = newreg; Plane Surface(s15) = {15};
s16 = newreg; Plane Surface(s16) = {16};
s17 = newreg; Plane Surface(s17) = {17};
s18 = newreg; Plane Surface(s18) = {18};
s19 = newreg; Plane Surface(s19) = {19};
s20 = newreg; Plane Surface(s20) = {20};
s21 = newreg; Plane Surface(s21) = {21};
s22 = newreg; Plane Surface(s22) = {22};
s23 = newreg; Plane Surface(s23) = {23};
s24 = newreg; Plane Surface(s24) = {24};
s25 = newreg; Plane Surface(s25) = {25};
s26 = newreg; Plane Surface(s26) = {26};
//---------------------------------------------------------------------------
// étirement diagonales
Transfinite Line{54,55,56,57}=arctheta1 Using Progression f;
//Transfinite Line{50,51,52,53}=arctheta2 Using Progression 1/f;
//---------------------------------------------------------------------------
// maillage structuré
Transfinite Line {19,59,61,13,58,10,16,5,7,60,61,6,8} = N1 Using
Progression 1.0;
Transfinite Line {23,13,19,48,28} = N1 Using Progression 1.0;
Transfinite Line {-37,31,36,-33,11,-15,-9,17} = N1 Using Progression 1.06;
Transfinite Line {20,18,35,34,38,32,12,14,22,24,47,49,27,29} = N2 Using
Progression 1.0;
Transfinite Line {21,41,42,25} = N2 Using Progression 1.0;
Transfinite Line {43,44,45,46,30,39,40,26} = N3 Using Progression 1.0;
// raffinement parois latérales
//Transfinite Line {-20,-29,35,-38,12,22,47,27,18,-34,32,-14,-24,-49} =
N2 Using Progression 1.05; //resserement vers la bord
Transfinite Line {20,29,-35,38,-12,-22,-47,-27,-18,34,-32,14,24,49} = N2
Using Progression 1.04; // étirement vers le bord
Transfinite Surface {s6} = {11,7,6,10};
Transfinite Surface {s7} = {11,7,8,12};
Transfinite Surface {s8} = {8,12,13,9};
Transfinite Surface {s9} = {9,13,10,6};
Transfinite Surface {s10} = {11,18,19,12};
Transfinite Surface {s11} = {12,21,22,13};
Transfinite Surface {s12} = {13,24,25,10};
Transfinite Surface {s13} = {11,10,15,16};
Transfinite Surface {s14} = {18,11,16,17};
Transfinite Surface {s15} = {12,19,20,21};
Transfinite Surface {s16} = {13,22,23,24};
Transfinite Surface {s17} = {10,25,14,15};
Transfinite Surface {s18} = {17,26,27,18};
Transfinite Surface {s19} = {18,27,28,19};
Transfinite Surface {s20} = {19,28,29,20};
Transfinite Surface {s21} = {14,25,31,30};
Transfinite Surface {s22} = {25,24,32,31};
Transfinite Surface {s23} = {23,24,32,33};
Transfinite Surface {s24} = {34,26,27,35};
Transfinite Surface {s25} = {36,35,27,28};
Transfinite Surface {s26} = {37,36,28,29};
Recombine Surface {s6:s26};
------------------------------------------------------------------------
Windows 7 : Trouvez le PC qui vous convient! Découvrez notre offre !
<http://clk.atdmt.com/FRM/go/181574577/direct/01/>
_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh
