Thank you very much, it works perfectly !!<br><br>Have a nice day, <br><br>Best regards,<br><br>Alex<br><br><div class="gmail_quote">Le 8 mars 2010 16:38, Ruth V. Sabariego <span dir="ltr"><<a href="mailto:r.sabariego@ulg.ac.be">r.sabariego@ulg.ac.be</a>></span> a écrit :<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">You have to define the holes at the cube faces. See corrected geometry attached.<br>
Regards,<br>
Ruth<div><div></div><div class="h5"><br>
<br>
On 08/03/10 16:16, Alexis Blasselle wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
// maillage d'un cube et du reseau de fibres interieur<br>
lc1 = 0.1 ;<br>
R = 0.05;<br>
H = 0.25 ;<br>
<br>
Point(100) = {0., 0., 0., lc1};<br>
// premier cercle<br>
Point(1) = { 0. , 0. , -H/2 , lc1 } ;<br>
Point(2) = { R , 0. , -H/2 , lc1 } ;<br>
Point(3) = { 0. , R , -H/2. , lc1 } ;<br>
Point(4) = { -R , 0. , -H/2 , lc1 } ;<br>
Point(5) = { 0. , -R , -H/2 , lc1 } ;<br>
<br>
Circle(23) = { 2 , 1 , 3 } ;<br>
Circle(34) = { 3 , 1 , 4 } ;<br>
Circle(45) = { 4 , 1 , 5 } ;<br>
Circle(52) = { 5 , 1 , 2 } ;<br>
<br>
<br>
// deuxieme cercle<br>
Point(6) = { 0. , 0. , H/2 , lc1 } ;<br>
Point(7) = { R , 0. , H/2 , lc1 } ;<br>
Point(8) = { 0. , R , H/2. , lc1 } ;<br>
Point(9) = { -R , 0. , H/2 , lc1 } ;<br>
Point(10) = { 0. , -R , H/2 , lc1 } ;<br>
<br>
Circle(78) = { 7 , 6 , 8 } ;<br>
Circle(89) = { 8 , 6 , 9 } ;<br>
Circle(910) = { 9 , 6 , 10 } ;<br>
Circle(107) = { 10 , 6 , 7 } ;<br>
<br>
<br>
<br>
// troisieme cercle<br>
Point(11) = { -H/2 , 0. , 0. , lc1 } ;<br>
Point(12) = { -H/2 , 0. , R , lc1 } ;<br>
Point(13) = { -H/2 , -R , 0 , lc1 } ;<br>
Point(14) = { -H/2 , 0. , -R , lc1 } ;<br>
Point(15) = { -H/2 , R , 0 , lc1 } ;<br>
<br>
Circle(1213) = { 12 , 11 , 13 } ;<br>
Circle(1314) = { 13 , 11 , 14 } ;<br>
Circle(1415) = { 14 , 11 , 15 } ;<br>
Circle(1512) = { 15 , 11 , 12 } ;<br>
<br>
<br>
<br>
// quatrieme cercle<br>
Point(16) = { H/2 , 0. , 0. , lc1 } ;<br>
Point(17) = { H/2 , 0. , R , lc1 } ;<br>
Point(18) = { H/2 , -R , 0 , lc1 } ;<br>
Point(19) = { H/2 , 0. , -R , lc1 } ;<br>
Point(20) = { H/2 , R , 0 , lc1 } ;<br>
<br>
Circle(1718) = { 17 , 16 , 18 } ;<br>
Circle(1819) = { 18 , 16 , 19 } ;<br>
Circle(1920) = { 19 , 16 , 20 } ;<br>
Circle(2017) = { 20 , 16 , 17 } ;<br>
<br>
<br>
<br>
// cinquieme cercle<br>
Point(21) = { 0. , -H/2 , 0. , lc1 } ;<br>
Point(22) = { 0 , -H/2 , R , lc1 } ;<br>
Point(23) = { -R , -H/2 , 0 , lc1 } ;<br>
Point(24) = { 0 , -H/2 , -R , lc1 } ;<br>
Point(25) = { R , -H/2 , 0 , lc1 } ;<br>
<br>
<br>
Circle(2223) = { 22 , 21 , 23 } ;<br>
Circle(2324) = { 23 , 21 , 24 } ;<br>
Circle(2425) = { 24 , 21 , 25 } ;<br>
Circle(2522) = { 25 , 21 , 22 } ;<br>
<br>
<br>
<br>
// sixieme cercle<br>
Point(26) = { 0. , H/2 , 0. , lc1 } ;<br>
Point(27) = { 0 , H/2 , R , lc1 } ;<br>
Point(28) = { -R , H/2 , 0 , lc1 } ;<br>
Point(29) = { 0 , H/2 , -R , lc1 } ;<br>
Point(30) = { R , H/2 , 0 , lc1 } ;<br>
<br>
Circle(2728) = { 27 , 26 , 28 } ;<br>
Circle(2829) = { 28 , 26 , 29 } ;<br>
Circle(2930) = { 29 , 26 , 30 } ;<br>
Circle(3027) = { 30 , 26 , 27 } ;<br>
<br>
// les cercles de la boule , plan x = 0<br>
Circle(192) = {19, 100, 2};<br>
Circle(414) = {4, 100, 14};<br>
Circle(129) = {12, 100, 9};<br>
Circle(717) = {7, 100, 17};<br>
<br>
// plan z = 0<br>
Circle(2030) = {20, 100, 30};<br>
Circle(2815) = {28, 100, 15};<br>
Circle(1323) = {13, 100, 23};<br>
Circle(2518) = {25, 100, 18};<br>
<br>
// plan y = 0<br>
Circle(329) = {3, 100, 29};<br>
Circle(278) = {27, 100, 8};<br>
Circle(1022) = {10, 100, 22};<br>
Circle(245) = {24, 100, 5};<br>
<br>
M = 2*H + H/2;<br>
<br>
// les points de la boiboite<br>
Point(101) = {M, M, -M, lc1};<br>
Point(102) = {-M, M, -M, lc1};<br>
Point(103) = {-M, -M, -M, lc1};<br>
Point(104) = {M, -M, -M, lc1};<br>
<br>
Point(105) = {M, M, M, lc1};<br>
Point(106) = {-M, M, M, lc1};<br>
Point(107) = {-M, -M, M, lc1};<br>
Point(108) = {M, -M, M, lc1};<br>
<br>
// carre du bas<br>
<br>
Line(101102) = {101, 102};<br>
Line(102103) = {102, 103};<br>
Line(103104) = {103, 104};<br>
Line(104101) = {104, 101};<br>
<br>
// carre du haut<br>
<br>
Line(105106) = {105, 106};<br>
Line(106107) = {106, 107};<br>
Line(107108) = {107, 108};<br>
Line(108105) = {108, 105};<br>
<br>
// les aretes<br>
<br>
Line(101105) = {101, 105};<br>
Line(102106) = {102, 106};<br>
Line(103107) = {103, 107};<br>
Line(104108) = {104, 108};<br>
<br>
// les line loop des faces du cube<br>
Line Loop(1234) = {101102, 102103, 103104, 104101};<br>
Plane Surface(1234) = {-1234};<br>
<br>
Line Loop(1265) = {101102, 102106, -105106, -101105};<br>
Plane Surface(1265) = {1265};<br>
<br>
Line Loop(2376) = {102103, 103107, -106107, -102106};<br>
Plane Surface(2376) = {2376};<br>
<br>
Line Loop(5678) = {105106, 106107, 107108, 108105};<br>
Plane Surface(5678) = {5678};<br>
<br>
Line Loop(3487) = {103104, 104108, -107108, -103107};<br>
Plane Surface(3487) = {3487};<br>
<br>
Line Loop(4158) = {104101, 101105, -108105, -104108};<br>
Plane Surface(4158) = {4158};<br>
<br>
<br>
// extrusions<br>
out23[] = Extrude {0,0,-2*H} { Line{23}; };<br>
out34[] = Extrude {0,0,-2*H} { Line{34}; };<br>
out45[] = Extrude {0,0,-2*H} { Line{45}; };<br>
out52[] = Extrude {0,0,-2*H} { Line{52}; };<br>
<br>
Line Loop (1) = {out23[0], out34[0], out45[0], out52[0] };<br>
<br>
out78[] = Extrude {0,0,2*H} { Line{78}; };<br>
out89[] = Extrude {0,0,2*H} { Line{89}; };<br>
out910[] = Extrude {0,0,2*H} { Line{910}; };<br>
out107[] = Extrude {0,0,2*H} { Line{107}; };<br>
<br>
Line Loop (2) = {out78[0], out89[0], out910[0], out107[0] };<br>
<br>
out1213[] = Extrude {-2*H,0,0} { Line{1213}; };<br>
out1314[] = Extrude {-2*H,0,0} { Line{1314}; };<br>
out1415[] = Extrude {-2*H,0,0} { Line{1415}; };<br>
out1512[] = Extrude {-2*H,0,0} { Line{1512}; };<br>
<br>
Line Loop (3) = {out1213[0], out1314[0], out1415[0], out1512[0] };<br>
<br>
out1718[] = Extrude {2*H,0,0} { Line{1718}; };<br>
out1819[] = Extrude {2*H,0,0} { Line{1819}; };<br>
out1920[] = Extrude {2*H,0,0} { Line{1920}; };<br>
out2017[] = Extrude {2*H,0,0} { Line{2017}; };<br>
<br>
Line Loop (4) = {out1718[0], out1819[0], out1920[0], out2017[0] };<br>
<br>
out2223[] = Extrude {0,-2*H,0} { Line{2223}; };<br>
out2324[] = Extrude {0,-2*H,0} { Line{2324}; };<br>
out2425[] = Extrude {0,-2*H,0} { Line{2425}; };<br>
out2522[] = Extrude {0,-2*H,0} { Line{2522}; };<br>
<br>
Line Loop (5) = {out2223[0], out2324[0], out2425[0], out2522[0] };<br>
<br>
out2728[] = Extrude {0,2*H,0} { Line{2728}; };<br>
out2829[] = Extrude {0,2*H,0} { Line{2829}; };<br>
out2930[] = Extrude {0,2*H,0} { Line{2930}; };<br>
out3027[] = Extrude {0,2*H,0} { Line{3027}; };<br>
<br>
Line Loop (6) = {out2728[0], out2829[0], out2930[0], out3027[0] };<br>
<br>
// les surfaces<br>
Plane Surface(1) = {1};<br>
Plane Surface(2) = {2};<br>
Plane Surface(3) = {3};<br>
Plane Surface(4) = {4};<br>
Plane Surface(5) = {5};<br>
Plane Surface(6) = {6};<br>
<br>
Circle(208) = {20, 100, 8};<br>
<br>
Line Loop(7) = {717, -2017, 208, -78};<br>
Ruled Surface(7) = {7};<br>
<br>
Line Loop(8) = {-208, 2030, 3027, 278};<br>
Ruled Surface(8) = {8};<br>
<br>
Circle(257) = {25, 100, 7};<br>
<br>
Line Loop(9) = {-107, 1022, -2522, 257};<br>
Ruled Surface(9) = {9};<br>
<br>
Line Loop(10) = {-257, 2518, -1718, -717};<br>
Ruled Surface(10) = {10};<br>
<br>
Circle(302) = {30, 100, 2};<br>
<br>
Line Loop(11) = {-302, -2030, -1920, 192};<br>
Ruled Surface(11) = {11};<br>
<br>
Line Loop(12) = {329, 2930, 302, 23};<br>
Ruled Surface(12) = {12};<br>
<br>
Circle(185) = {18, 100, 5};<br>
<br>
Line Loop(13) = {-192, -1819, 185, 52};<br>
Ruled Surface(13) = {13};<br>
<br>
Line Loop(14) = {-2518, -2425, 245, -185};<br>
Ruled Surface(14) = {14};<br>
<br>
Circle(239) = {23, 100, 9};<br>
<br>
Line Loop(15) = {-129, 1213, 1323, 239};<br>
Ruled Surface(15) = {15};<br>
<br>
Line Loop(16) = {2223, 239, 910, 1022};<br>
Ruled Surface(16) = {16};<br>
<br>
Circle(2914) = {29, 100, 14};<br>
<br>
Line Loop(17) = {-2914, -329, 34, 414};<br>
Ruled Surface(17) = {17};<br>
<br>
Line Loop(18) = {-2815, 2829, 2914, 1415};<br>
Ruled Surface(18) = {18};<br>
<br>
Circle(234) = {23, 100, 4};<br>
<br>
Line Loop(19) = {234, 45, -245, -2324};<br>
Ruled Surface(19) = {19};<br>
<br>
Line Loop(20) = {-1323, 1314, -414, -234};<br>
Ruled Surface(20) = {20};<br>
<br>
Circle(289) = {28, 100, 9};<br>
<br>
Line Loop(21) = {-89, -278, 2728, 289};<br>
Ruled Surface(21) = {21};<br>
<br>
Line Loop(22) = {2815, 1512, 129, -289};<br>
Ruled Surface(22) = {22};<br>
<br>
Surface Loop(23) = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22};<br>
//Volume(23) = {23};<br>
<br>
// les faces du cube<br>
Surface Loop(24) = {1234, 1265, 2376, 5678, 3487, 4158};<br>
Volume(24) = {24, -23};<br>
</blockquote>
<br>
-- <br></div></div>
Dr. Ir. Ruth V. Sabariego<br>
University of Liege, Dept. of Electrical Engineering& Computer Science,<br>
Applied& Computational Electromagnetics (ACE),<br>
phone: +32-4-3663737 - fax: +32-4-3662910 - <a href="http://ace.montefiore.ulg.ac.be/" target="_blank">http://ace.montefiore.ulg.ac.be/</a><br>
<br>
<br>_______________________________________________<br>
gmsh mailing list<br>
<a href="mailto:gmsh@geuz.org">gmsh@geuz.org</a><br>
<a href="http://www.geuz.org/mailman/listinfo/gmsh" target="_blank">http://www.geuz.org/mailman/listinfo/gmsh</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><a href="http://www.ann.jussieu.fr/~blasselle/">http://www.ann.jussieu.fr/~blasselle/</a><br>