Skip to content
Snippets Groups Projects
Commit ed9c60ca authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

some benchmarks contributed by other users

parent ae09b447
No related branches found
No related tags found
No related merge requests found
unit = 1.e-2;
size = 0.1 * unit;
Point(1) = {4.8*unit,0.8*unit,0, size};
Point(2) = {4.8*unit,1.2*unit,0, size/2};
Point(3) = {4.8*unit,1.4*unit,0, size/2};
Point(4) = {4*unit,1.4*unit,0, size/2};
Point(5) = {4*unit,1.2*unit,0, size/2};
Point(6) = {3.5*unit,1.2*unit,0, size};
Point(7) = {2.1*unit,0.4*unit,0, size*2};
Point(8) = {0,0.4*unit,0, size*2};
Point(9) = {0,0,0, size*2};
Point(10) = {2.2*unit,0,0, size*2};
Point(11) = {3.6*unit,0.8*unit,0, size};
Point(13) = {0.0445,0.008,0,size};
Point(15) = {0.0435,0.008,0,size};
Point(16) = {0.0435,0.012,0,size/4};
Point(18) = {0.0445,0.012,0,size/4};
Point(19) = {0.058,0.008,0,size};
Point(20) = {0.058,0.010,0,size};
Point(21) = {0.048,0.010,0,size};
Point(22) = {0.0435,0.01,0,size/4};
Point(23) = {0.0445,0.01,0,size/4};
Line(1) = {9,10};
Line(2) = {10,11};
Line(3) = {11,15};
Line(4) = {15,22};
Line(5) = {22,23};
Line(6) = {23,13};
Line(7) = {13,1};
Line(8) = {1,19};
Line(9) = {19,20};
Line(10) = {20,21};
Line(11) = {21,2};
Line(12) = {2,18};
Line(13) = {18,16};
Line(14) = {16,5};
Line(15) = {5,6};
Line(16) = {6,7};
Line(17) = {7,8};
Line(18) = {8,9};
Line(19) = {5,4};
Line(20) = {4,3};
Line(21) = {3,2};
Line(22) = {22,16};
Line(23) = {23,18};
Line Loop(24) = {15,16,17,18,1,2,3,4,22,14};
Plane Surface(25) = {24};
Line Loop(26) = {23,13,-22,5};
Plane Surface(27) = {26};
Line Loop(28) = {12,-23,6,7,8,9,10,11};
Plane Surface(29) = {28};
Line Loop(30) = {14,19,20,21,12,13};
Plane Surface(31) = {30};
/* entree flux */
/* pastille */
/* contact fixe */
/* contact fixe ou air */
Physical Line(1000) = {20};
Physical Surface(2000) = {31};
Physical Surface(3000) = {25,29};
Physical Surface(3001) = {27};
This diff is collapsed.
size = 0.002;
z = 0.006 ;
Point(1) = {0.048, 0.008, 0, size};
Point(2) = {0.048, 0.012, 0, size/2};
Point(3) = {0.048, 0.014, 0, size/2};
Point(4) = {0.04, 0.014, 0, size/2};
Point(5) = {0.04, 0.012, 0, size/2};
Point(6) = {0.035, 0.012, 0, size};
Point(7) = {0.021, 0.004, 0, size*2};
Point(8) = {0, 0.004, 0, size*2};
Point(9) = {0, 0, 0, size*2};
Point(10) = {0.022, 0, 0, size*2};
Point(11) = {0.036, 0.008, 0, size};
Point(12) = {0.058, 0.008, 0, size};
Point(13) = {0.058, 0.010, 0, size};
Point(14) = {0.048, 0.010, 0, size};
Point(21) = {0.048, 0.008, z, size};
Point(22) = {0.048, 0.012, z, size/2};
Point(23) = {0.048, 0.014, z, size/2};
Point(24) = {0.04, 0.014, z, size/2};
Point(25) = {0.04, 0.012, z, size/2};
Point(26) = {0.035, 0.012, z, size};
Point(27) = {0.021, 0.004, z, size*2};
Point(28) = {0, 0.004, z, size*2};
Point(29) = {0, 0, z, size*2};
Point(30) = {0.022, 0, z, size*2};
Point(31) = {0.036, 0.008, z, size};
Point(32) = {0.058, 0.008, z, size};
Point(33) = {0.058, 0.010, z, size};
Point(34) = {0.048, 0.010, z, size};
/* pts suppl dans la pastille */
Point(35) = {0.048, 0.014, z/3, size/2};
Point(36) = {0.04, 0.014, z/3, size/2};
Point(37) = {0.048, 0.014, 2*z/3, size/2};
Point(38) = {0.04, 0.014, 2*z/3, size/2};
/* trou de mesure */
Point(40) = {0.044, 0.008, z/2, size/2};
Point(41) = {0.044, 0.011, z/2, size/2};
Point(42) = {0.044, 0.012, z/2, size/2};
Point(43) = {0.0435, 0.008, z/2, size/4};
Point(44) = {0.0445, 0.008, z/2, size/4};
Point(45) = {0.0435, 0.011, z/2, size/4};
Point(46) = {0.0445, 0.011, z/2, size/4};
Point(47) = {0.0435, 0.012, z/2, size/4};
Point(48) = {0.0445, 0.012, z/2, size/4};
Point(49) = {0.044, 0.008, z/2-0.0005, size/4};
Point(50) = {0.044, 0.008, z/2+0.0005, size/4};
Point(51) = {0.044, 0.011, z/2-0.0005, size/4};
Point(52) = {0.044, 0.011, z/2+0.0005, size/4};
Point(53) = {0.044, 0.012, z/2-0.0005, size/4};
Point(54) = {0.044, 0.012, z/2+0.0005, size/4};
/* trou trou de droite */
Point(60) = {0.053, 0.008, z/2, size/2};
Point(61) = {0.053, 0.01, z/2, size/2};
Point(62) = {0.052, 0.008, z/2, size/2};
Point(63) = {0.054, 0.008, z/2, size/2};
Point(64) = {0.052, 0.01, z/2, size/2};
Point(65) = {0.054, 0.01, z/2, size/2};
Point(66) = {0.053, 0.008, z/2-0.001, size/2};
Point(67) = {0.053, 0.008, z/2+0.001, size/2};
Point(68) = {0.053, 0.01, z/2-0.001, size/2};
Point(69) = {0.053, 0.01, z/2+0.001, size/2};
Line(101) = {9,10};
Line(102) = {10,11};
Line(103) = {11,1};
Line(104) = {1,12};
Line(105) = {12,32};
Line(106) = {32,21};
Line(107) = {21,31};
Line(108) = {31,30};
Line(109) = {30,29};
Line(110) = {29,9};
Line(111) = {30,10};
Line(112) = {31,11};
Line(113) = {29,28};
Line(114) = {8,28};
Line(115) = {8,9};
Line(116) = {8,7};
Line(117) = {7,10};
Line(118) = {27,7};
Line(119) = {27,30};
Line(120) = {27,28};
Line(121) = {27,26};
Line(122) = {7,6};
Line(123) = {6,26};
Line(124) = {26,31};
Line(125) = {11,6};
Line(126) = {26,25};
Line(127) = {25,5};
Line(128) = {5,6};
Line(129) = {25,24};
Line(130) = {24,38};
Line(131) = {38,36};
Line(132) = {36,4};
Line(133) = {4,5};
Line(134) = {22,2};
Line(135) = {2,3};
Line(136) = {3,35};
Line(137) = {35,37};
Line(138) = {37,23};
Line(139) = {23,22};
Line(140) = {24,23};
Line(141) = {37,38};
Line(142) = {36,35};
Line(143) = {3,4};
Line(144) = {25,22};
Line(145) = {5,2};
Line(146) = {14,2};
Line(147) = {22,34};
Line(148) = {34,14};
Line(149) = {34,33};
Line(150) = {33,13};
Line(151) = {13,14};
Line(152) = {33,32};
Line(153) = {12,13};
Line(154) = {21,34};
Line(155) = {1,14};
Line(156) = {21,1};
Circle(157) = {48,42,53};
Circle(158) = {53,42,47};
Circle(159) = {47,42,54};
Circle(160) = {54,42,48};
Circle(161) = {51,41,45};
Circle(162) = {45,41,52};
Circle(163) = {52,41,46};
Circle(164) = {46,41,51};
Circle(165) = {44,40,49};
Circle(166) = {49,40,43};
Circle(167) = {43,40,50};
Circle(168) = {50,40,44};
Circle(169) = {65,61,68};
Circle(170) = {68,61,64};
Circle(171) = {64,61,69};
Circle(172) = {69,61,65};
Circle(173) = {66,60,62};
Circle(174) = {62,60,67};
Circle(175) = {67,60,63};
Circle(176) = {63,60,66};
Line(177) = {64,62};
Line(178) = {69,67};
Line(179) = {63,65};
Line(180) = {66,68};
Line(181) = {43,45};
Line(182) = {49,51};
Line(183) = {44,46};
Line(184) = {50,52};
Line(185) = {52,54};
Line(186) = {46,48};
Line(187) = {51,53};
Line(188) = {45,47};
Line Loop(189) = {-114,115,-110,113};
Plane Surface(190) = {189};
Line Loop(191) = {-114,116,-118,120};
Plane Surface(192) = {191};
Line Loop(193) = {-101,-115,116,117};
Plane Surface(194) = {193};
Line Loop(195) = {109,113,-120,119};
Plane Surface(196) = {195};
Line Loop(197) = {-101,-110,-109,111};
Plane Surface(198) = {197};
Line Loop(199) = {-123,-122,-118,121};
Plane Surface(200) = {199};
Line Loop(201) = {-122,117,102,125};
Plane Surface(202) = {201};
Line Loop(203) = {108,-119,121,124};
Plane Surface(204) = {203};
Line Loop(205) = {-112,108,111,102};
Plane Surface(206) = {205};
Line Loop(207) = {-128,145,-146,-155,-103,125};
Plane Surface(208) = {207};
Line Loop(209) = {-103,-112,-107,156};
Line Loop(210) = {168,165,166,167};
Plane Surface(211) = {209,210};
Plane Surface(212) = {210};
Line Loop(213) = {162,163,164,161};
Plane Surface(214) = {213};
Line Loop(215) = {159,160,157,158};
Plane Surface(216) = {215};
Line Loop(217) = {123,126,127,128};
Plane Surface(218) = {217};
Line Loop(219) = {-134,-144,127,145};
Plane Surface(220) = {219,215};
Line Loop(221) = {147,-154,107,-124,126,144};
Plane Surface(222) = {221};
Line Loop(223) = {-148,-154,156,155};
Plane Surface(224) = {223};
Line Loop(225) = {-134,147,148,146};
Plane Surface(226) = {225};
Line Loop(227) = {-127,129,130,131,132,133};
Plane Surface(228) = {227};
Line Loop(229) = {139,-144,129,140};
Plane Surface(230) = {229};
Line Loop(231) = {136,137,138,139,134,135};
Plane Surface(232) = {231};
Line Loop(233) = {145,135,143,133};
Plane Surface(234) = {233};
Line Loop(235) = {-132,142,-136,143};
Plane Surface(236) = {235};
Line Loop(237) = {137,141,131,142};
Plane Surface(238) = {237};
Line Loop(239) = {-130,140,-138,141};
Plane Surface(240) = {239};
Line Loop(241) = {152,106,154,149};
Plane Surface(242) = {241};
Line Loop(243) = {104,105,106,156};
Line Loop(244) = {174,175,176,173};
Plane Surface(245) = {243,244};
Plane Surface(246) = {244};
Line Loop(247) = {-151,-153,-104,155};
Plane Surface(248) = {247};
Line Loop(249) = {-150,152,-105,153};
Plane Surface(250) = {249};
Line Loop(251) = {-151,-150,-149,148};
Line Loop(252) = {171,172,169,170};
Plane Surface(253) = {251,252};
Plane Surface(254) = {252};
Line Loop(255) = {174,-178,-171,177};
Ruled Surface(256) = {255};
Line Loop(257) = {-179,-175,-178,172};
Ruled Surface(258) = {257};
Line Loop(259) = {-169,-179,176,180};
Ruled Surface(260) = {259};
Line Loop(261) = {-173,180,170,177};
Ruled Surface(262) = {261};
Line Loop(263) = {162,-184,-167,181};
Ruled Surface(264) = {263};
Line Loop(265) = {-183,-168,184,163};
Ruled Surface(266) = {265};
Line Loop(267) = {182,-164,-183,165};
Ruled Surface(268) = {267};
Line Loop(269) = {-181,-166,182,161};
Ruled Surface(270) = {269};
Line Loop(271) = {-159,-188,162,185};
Ruled Surface(272) = {271};
Line Loop(273) = {-160,-185,163,186};
Ruled Surface(274) = {273};
Line Loop(275) = {157,-187,-164,186};
Ruled Surface(276) = {275};
Line Loop(277) = {188,-158,-187,161};
Ruled Surface(278) = {277};
Surface Loop(279) = {266,-268,270,264,-272,-220,226,-222,224,211,208,218,200,-202,194,-198,190,-192,-196,-204,-206,-274,276,-278};
Complex Volume(280) = {279};
Surface Loop(281) = {258,-260,-253,248,250,-242,245,-224,-256,262};
Complex Volume(282) = {281};
Surface Loop(283) = {214,-272,-216,-274,276,-278};
Complex Volume(284) = {283};
Surface Loop(285) = {216,-220,-232,-236,-228,230,-240,238,234};
Complex Volume(286) = {285};
/* entree flux */
Physical Surface(1000) = {236,238 /*,240 */ };
/* pastille */
Physical Volume(2000) = {286};
/* cuivre */
Physical Volume(3000) = {280,282 /* ,284 */};
//Fichier test
mm = 1e-3;
s = 1.;
prec4 = 15 * mm * s;
prec1 = 7 * mm * s;
prec2 = 10 * mm * s;
prec3 = 4 * mm * s;
//Donnees geometriques
d = 20 * mm;
largt = 50 * mm;
longt = 60 * mm;
h = 40 * mm;
c = 25 * mm;
k = 1.5;
bl = 10 * mm;
bk = k * bl;
long_ei = largt + 2 * c;
long_ii = long_ei - 2 * bk;
larg_ei = 2 * bk;
larg_ii = larg_ei - 2 * bl;
ep_i = 10 * mm;
htot = 200 * mm;
mult = 1.5;
// pts ind 1
Point(1) = { larg_ii/2. , 0 , h/2. , prec1 };
Point(2) = { larg_ii/2. , long_ii/2. - larg_ii/2. , h/2. , prec1 };
Point(3) = { 0 , long_ii/2. , h/2. , prec1 };
Point(4) = { -larg_ii/2. , long_ii/2. - larg_ii/2. , h/2. , prec1 };
Point(5) = { -larg_ii/2. , 0 , h/2. , prec1 };
Point(6) = { larg_ii/2. , 0 , h/2. + ep_i , prec1 };
Point(7) = { larg_ii/2. , long_ii/2. - larg_ii/2. , h/2. + ep_i , prec1 };
Point(8) = { 0 , long_ii/2. , h/2. + ep_i , prec1 };
Point(9) = { -larg_ii/2. , long_ii/2. - larg_ii/2. , h/2. + ep_i , prec1 };
Point(10) = { -larg_ii/2. , 0 , h/2. + ep_i , prec1 };
Point(11) = { 0 , long_ii/2. - larg_ii/2 , h/2. , prec1 };
Point(12) = { 0 , long_ii/2. - larg_ii/2 , h/2.+ ep_i , prec1 };
Point(13) = { larg_ei/2. , 0 , h/2. , prec1 };
Point(14) = { larg_ei/2. , long_ei/2. - larg_ei/2. , h/2. , prec1 };
Point(15) = { 0 , long_ei/2. , h/2. , prec1 };
Point(16) = { -larg_ei/2. , long_ei/2. - larg_ei/2. , h/2. , prec1 };
Point(17) = { -larg_ei/2. , 0 , h/2. , prec1 };
Point(18) = { larg_ei/2. , 0 , h/2. + ep_i , prec1 };
Point(19) = { larg_ei/2. , long_ei/2. - larg_ei/2. , h/2. + ep_i , prec1 };
Point(20) = { 0 , long_ei/2. , h/2. + ep_i , prec1 };
Point(21) = { -larg_ei/2. , long_ei/2. - larg_ei/2. , h/2. + ep_i , prec1 };
Point(22) = { -larg_ei/2. , 0 , h/2. + ep_i , prec1 };
h2 =h+3*ep_i;
//Pts ind 2
Point(401) = { larg_ii/2. , 0 , h2/2. , prec1 };
Point(402) = { larg_ii/2. , long_ii/2. - larg_ii/2. , h2/2. , prec1 };
Point(403) = { 0 , long_ii/2. , h2/2. , prec1 };
Point(404) = { -larg_ii/2. , long_ii/2. - larg_ii/2. , h2/2. , prec1 };
Point(405) = { -larg_ii/2. , 0 , h2/2. , prec1 };
Point(406) = { larg_ii/2. , 0 , h2/2. + ep_i , prec1 };
Point(407) = { larg_ii/2. , long_ii/2. - larg_ii/2. , h2/2. + ep_i , prec1 };
Point(408) = { 0 , long_ii/2. , h2/2. + ep_i , prec1 };
Point(409) = { -larg_ii/2. , long_ii/2. - larg_ii/2. , h2/2. + ep_i , prec1 };
Point(410) = { -larg_ii/2. , 0 , h2/2. + ep_i , prec1 };
Point(411) = { 0 , long_ii/2. - larg_ii/2 , h2/2. , prec1 };
Point(412) = { 0 , long_ii/2. - larg_ii/2 , h2/2.+ ep_i , prec1 };
Point(413) = { larg_ei/2. , 0 , h2/2. , prec1 };
Point(414) = { larg_ei/2. , long_ei/2. - larg_ei/2. , h2/2. , prec1 };
Point(415) = { 0 , long_ei/2. , h2/2. , prec1 };
Point(416) = { -larg_ei/2. , long_ei/2. - larg_ei/2. , h2/2. , prec1 };
Point(417) = { -larg_ei/2. , 0 , h2/2. , prec1 };
Point(418) = { larg_ei/2. , 0 , h2/2. + ep_i , prec1 };
Point(419) = { larg_ei/2. , long_ei/2. - larg_ei/2. , h2/2. + ep_i , prec1 };
Point(420) = { 0 , long_ei/2. , h2/2. + ep_i , prec1 };
Point(421) = { -larg_ei/2. , long_ei/2. - larg_ei/2. , h2/2. + ep_i , prec1 };
Point(422) = { -larg_ei/2. , 0 , h2/2. + ep_i , prec1 };
// points volume air et points tole
largeur = 100 * mm;
longueur = 100 * mm;
hauteur = 2./3.*100 * mm;
ept = 1 * mm;
Point(23) = { -largeur/2. , 0 , 0 , prec3 };
Point(25) = { largeur/2. , 0 , 0 , prec3 };
Point(26) = { largeur/2. , longueur , 0 , prec3 };
Point(27) = { -largeur/2. , longueur , 0 , prec3 };
Point(28) = { -largeur/2. , 0 , ept , prec3 };
Point(29) = { largeur/2. , 0 , ept , prec3 };
Point(30) = { largeur/2. , longueur , ept , prec3 };
Point(31) = { -largeur/2. , longueur , ept , prec3 };
Point(32) = { -largeur/2. , 0 , hauteur , prec4 };
Point(33) = { largeur/2. , 0 , hauteur , prec4 };
Point(34) = { largeur/2. , longueur , hauteur , prec4 };
Point(35) = { -largeur/2. , longueur , hauteur , prec4 };
Line(1) = {1,13};
Line(2) = {13,18};
Line(3) = {18,6};
Line(4) = {6,1};
Line(5) = {5,10};
Line(6) = {10,22};
Line(7) = {22,17};
Line(8) = {17,5};
Line(9) = {6,7};
Circle(10) = {7,12,8};
Circle(11) = {8,12,9};
Line(12) = {9,9,10};
Line(13) = {1,2};
Circle(14) = {2,11,3};
Circle(15) = {3,11,4};
Line(16) = {4,5};
Line(17) = {13,14};
Circle(18) = {14,3,15};
Circle(19) = {15,3,16};
Line(20) = {16,17};
Line(21) = {18,19};
Circle(22) = {19,8,20};
Circle(23) = {20,8,21};
Line(24) = {21,22};
Line(25) = {23,25};
Line(26) = {25,26};
Line(27) = {26,27};
Line(28) = {27,23};
Line(29) = {28,29};
Line(30) = {29,30};
Line(31) = {30,31};
Line(32) = {31,28};
Line(33) = {33,34};
Line(34) = {34,35};
Line(35) = {35,32};
Line(36) = {32,33};
Line(37) = {25,29};
Line(38) = {29,33};
Line(39) = {26,30};
Line(40) = {30,34};
Line(41) = {27,31};
Line(42) = {31,35};
Line(43) = {23,28};
Line(44) = {28,32};
Line(45) = {14,19};
Line(46) = {2,7};
Line(47) = {4,9};
Line(48) = {16,21};
Line(49) = {15,20};
Line(50) = {3,8};
Line(51) = {5,1};
Line(52) = {10,6};
Line(1002) = {401,413};
Line(1003) = {413,418};
Line(1004) = {418,406};
Line(1005) = {406,401};
Line(1006) = {405,410};
Line(1007) = {410,422};
Line(1008) = {422,417};
Line(1009) = {417,405};
Line(1010) = {418,419};
Line(1012) = {414,413};
Line(1013) = {422,421};
Line(1014) = {417,416};
Line(1015) = {406,407};
Line(1016) = {401,402};
Line(1017) = {405,404};
Line(1018) = {410,409};
Circle(1019) = {404,411,403};
Circle(1020) = {403,411,402};
Circle(1023) = {416,403,415};
Circle(1024) = {415,403,414};
Circle(1025) = {409,412,408};
Circle(1026) = {408,412,407};
Circle(1027) = {421,408,420};
Circle(1028) = {420,408,419};
Line(1029) = {415,415,420};
Line(1030) = {408,403};
Line(1037) = {410,406};
Line(1038) = {405,401};
Line(1041) = {407,402};
Line(1042) = {414,419};
Line(1043) = {404,409};
Line(1044) = {416,421};
Line Loop(53) = {25,26,27,28};
Plane Surface(54) = {53};
Line Loop(55) = {30,-39,-26,37};
Plane Surface(56) = {55};
Line Loop(57) = {-31,-39,27,41};
Plane Surface(58) = {57};
Line Loop(59) = {43,-32,-41,28};
Plane Surface(60) = {59};
Line Loop(61) = {-37,-25,43,29};
Plane Surface(62) = {61};
Line Loop(63) = {33,-40,-30,38};
Plane Surface(64) = {63};
Line Loop(65) = {-34,-40,31,42};
Plane Surface(66) = {65};
Line Loop(67) = {44,-35,-42,32};
Plane Surface(68) = {67};
Line Loop(69) = {36,-38,-29,44};
//Line Loop(70) = {1,2,3,-52,6,7,8,51};
//Plane Surface(71) = {69,70};
Line Loop(72) = {29,30,31,32};
Plane Surface(73) = {72};
Line Loop(74) = {1,2,3,4};
Plane Surface(75) = {74};
Line Loop(76) = {52,4,-51,5};
Plane Surface(77) = {76};
Line Loop(78) = {7,8,5,6};
Plane Surface(79) = {78};
Line Loop(80) = {21,-45,-17,2};
Plane Surface(81) = {80};
Line Loop(82) = {-49,-18,45,22};
Ruled Surface(83) = {82};
Line Loop(84) = {-23,-49,19,48};
Ruled Surface(85) = {84};
Line Loop(86) = {-20,48,24,7};
Ruled Surface(87) = {86};
Line Loop(88) = {-10,-46,14,50};
Ruled Surface(89) = {88};
Line Loop(90) = {-11,-50,15,47};
Ruled Surface(91) = {90};
Line Loop(92) = {13,46,-9,4};
Plane Surface(93) = {92};
Line Loop(94) = {5,-12,-47,16};
Plane Surface(95) = {94};
Line Loop(96) = {8,-16,-15,-14,-13,1,17,18,19,20};
Plane Surface(97) = {96};
Line Loop(98) = {51,13,14,15,16};
Plane Surface(99) = {98};
Line Loop(100) = {52,9,10,11,12};
Plane Surface(101) = {100};
Line Loop(102) = {10,11,12,6,-24,-23,-22,-21,3,9};
Plane Surface(103) = {102};
Line Loop(104) = {36,33,34,35};
Plane Surface(105) = {104};
Line Loop(1031) = {1006,1007,1008,1009};
Plane Surface(1032) = {1031};
Line Loop(1033) = {1002,1003,1004,1005};
Plane Surface(1034) = {1033};
Line Loop(1035) = {-1004,1010,-1028,-1027,-1013,-1007,1018,1025,1026,-1015};
Plane Surface(1036) = {1035};
Line Loop(1039) = {1020,-1016,1002,-1012,-1024,-1023,-1014,1009,1017,1019};
Plane Surface(1040) = {1039};
Line Loop(1045) = {1028,-1042,-1024,1029};
Ruled Surface(1046) = {1045};
Line Loop(1047) = {1029,-1027,-1044,1023};
Ruled Surface(1048) = {1047};
Line Loop(1049) = {1030,-1019,1043,1025};
Ruled Surface(1050) = {1049};
Line Loop(1051) = {-1020,-1030,1026,1041};
Ruled Surface(1052) = {1051};
Line Loop(1053) = {1010,-1042,1012,1003};
Plane Surface(1054) = {1053};
Line Loop(1055) = {1016,-1041,-1015,1005};
Plane Surface(1056) = {1055};
Line Loop(1057) = {-1013,1008,1014,1044};
Plane Surface(1058) = {1057};
Line Loop(1059) = {-1018,-1006,1017,1043};
Plane Surface(1060) = {1059};
Line Loop(1061) = {1037,1005,-1038,1006};
Plane Surface(1062) = {1061};
Line Loop(1063) = {1019,1020,-1016,-1038,1017};
Plane Surface(1064) = {1063};
Line Loop(1065) = {1026,-1015,-1037,1018,1025};
Plane Surface(1066) = {1065};
Line Loop(1079) = {8,51,1,2,3,-52,6,7};
Line Loop(1080) = {1038,1002,1003,1004,-1037,1007,1008,1009};
Plane Surface(1081) = {69,1079,1080};
Surface Loop(1067) = {1060,1066,-1052,-1064,-1050,-1056,1062};
Complex Volume(1068) = {1067};
Surface Loop(1070) = {1060,1036,1034,-1040,-1052,-1050,-1056,-1054,1046,-1048,-1058,1032};
Complex Volume(1071) = {1070};
Surface Loop(1072) = {89,101,-77,93,-99,91,95};
Complex Volume(1073) = {1072};
Surface Loop(1074) = {95,-79,87,97,91,103,89,93,-75,81,83,-85};
Complex Volume(1075) = {1074};
Surface Loop(1076) = {62,56,-58,54,-60,-73};
Complex Volume(1077) = {1076};
Surface Loop(1082) = {1066,-1036,-1081,105,-64,66,-73,68,97,99,81,103,-101,87,-85,83,-1064,1040,1054,-1046,1048,1058};
Complex Volume(1083) = {1082};
//Physical datas
Physical Volume(502) = { 1075 , 1071 }; //ind inf et sup
Physical Volume(501) = { 1077 }; //tole
Physical Volume(500) = { 1083 }; // Volume air
Physical Volume(550) = { 1073 }; //boite air inf
Physical Volume(551) = { 1068 }; //boite air sup
Physical Surface(600) = {1054,1046,1048,1058,1060,1050,1052,1056,1036,1040,97,87,85,83,81,103,93,95,91,89}; //skinind
Physical Surface(601) = { 99,1064 }; //cutind
Physical Surface(602) = {93,89,91,95,1056,1052,1050,1060}; //skindomainc2_Mag
Physical Surface(603) = { 73 }; //skintole
Physical Surface(1001) = {64,66,68,105,54,58}; //gh0
Physical Surface(700) = {75,1034}; //elec0
Physical Surface(701) = {79,1032}; //elec1
mm = 1e-3;
s = 0.6;
prec1 = 20 * mm * s;
prec2 = 35 * mm * s;
prec3 = 60 * mm * s;
/* Donnees geometriques */
d = 10 * mm;
largt = 500 * mm;
longt = 600 * mm;
h = 50 * mm;
c = 30 * mm;
k = 1.5;
bl = 100 * mm;
bk = k * bl;
long_ei = largt + 2 * c;
long_ii = long_ei - 2 * bk;
larg_ei = 2 * bk;
larg_ii = larg_ei - 2 * bl;
ep_i = 10 * mm;
htot = 200 * mm;
mult = 1.5;
/* pts ind */
Point(1) = { larg_ii/2. , 0 , h/2. , prec1 };
Point(2) = { larg_ii/2. , long_ii/2. - larg_ii/2. , h/2. , prec1 };
Point(3) = { 0 , long_ii/2. , h/2. , prec1 };
Point(4) = { -larg_ii/2. , long_ii/2. - larg_ii/2. , h/2. , prec1 };
Point(5) = { -larg_ii/2. , 0 , h/2. , prec1 };
Point(6) = { larg_ii/2. , 0 , h/2. + ep_i , prec1 };
Point(7) = { larg_ii/2. , long_ii/2. - larg_ii/2. , h/2. + ep_i , prec1 };
Point(8) = { 0 , long_ii/2. , h/2. + ep_i , prec1 };
Point(9) = { -larg_ii/2. , long_ii/2. - larg_ii/2. , h/2. + ep_i , prec1 };
Point(10) = { -larg_ii/2. , 0 , h/2. + ep_i , prec1 };
Point(11) = { 0 , long_ii/2. - larg_ii/2 , h/2. , prec1 };
Point(12) = { 0 , long_ii/2. - larg_ii/2 , h/2.+ ep_i , prec1 };
Point(13) = { larg_ei/2. , 0 , h/2. , prec1 };
Point(14) = { larg_ei/2. , long_ei/2. - larg_ei/2. , h/2. , prec1 };
Point(15) = { 0 , long_ei/2. , h/2. , prec1 };
Point(16) = { -larg_ei/2. , long_ei/2. - larg_ei/2. , h/2. , prec1 };
Point(17) = { -larg_ei/2. , 0 , h/2. , prec1 };
Point(18) = { larg_ei/2. , 0 , h/2. + ep_i , prec1 };
Point(19) = { larg_ei/2. , long_ei/2. - larg_ei/2. , h/2. + ep_i , prec1 };
Point(20) = { 0 , long_ei/2. , h/2. + ep_i , prec1 };
Point(21) = { -larg_ei/2. , long_ei/2. - larg_ei/2. , h/2. + ep_i , prec1 };
Point(22) = { -larg_ei/2. , 0 , h/2. + ep_i , prec1 };
/* pts tole */
Point(23) = { - longt/2. , 0 , 0 , prec1 };
Point(24) = { longt/2. , 0 , 0 , prec1 };
Point(25) = { longt/2. , largt/2. , 0 , prec1 };
Point(26) = { - longt/2. , largt/2. , 0 , prec1 };
Point(27) = { - longt/2. , 0 , d/2. , prec1 };
Point(28) = { longt/2. , 0 , d/2. , prec1 };
Point(29) = { longt/2. , largt/2. , d/2. , prec1 };
Point(30) = { - longt/2. , largt/2. , d/2. , prec1 };
/* pts supp contour air */
Point(31) = { longt/2. , mult * largt/2. , 0 , prec2 };
Point(32) = { - longt/2. , mult * largt/2. , 0 , prec2 };
Point(33) = { longt/2. , mult * largt/2. , htot , prec3 };
Point(34) = { - longt/2. , mult * largt/2. , htot , prec3 };
Point(35) = { longt/2. , 0 , htot , prec3 };
Point(36) = { - longt/2. , 0 , htot , prec3 };
Line(1) = {1,13};
Line(2) = {13,18};
Line(3) = {18,6};
Line(4) = {6,1};
Line(5) = {5,10};
Line(6) = {10,22};
Line(7) = {22,17};
Line(8) = {17,5};
Line(9) = {6,7};
Circle(10) = {7,12,8};
Circle(11) = {8,12,9};
Line(12) = {9,9,10};
Line(13) = {1,2};
Circle(14) = {2,11,3};
Circle(15) = {3,11,4};
Line(16) = {4,5};
Line(17) = {13,14};
Circle(18) = {14,3,15};
Circle(19) = {15,3,16};
Line(20) = {16,17};
Line(21) = {18,19};
Circle(22) = {19,8,20};
Circle(23) = {20,8,21};
Line(24) = {21,22};
Line(25) = {23,27};
Line(26) = {27,36};
Line(27) = {26,30};
Line(28) = {23,26};
Line(29) = {26,32};
Line(30) = {32,34};
Line(31) = {34,36};
Line(32) = {30,27};
Line(33) = {32,31};
Line(34) = {31,33};
Line(35) = {33,34};
Line(36) = {33,35};
Line(37) = {35,36};
Line(38) = {31,25};
Line(39) = {25,24};
Line(40) = {24,28};
Line(41) = {28,35};
Line(42) = {25,29};
Line(43) = {29,28};
Line(44) = {28,27};
Line(45) = {24,23};
Line(46) = {29,30};
Line(47) = {25,26};
Line(62) = {9,4};
Line(63) = {8,3};
Line(64) = {7,2};
Line(65) = {21,16};
Line(66) = {20,15};
Line(67) = {19,14};
Line(68) = {5,1};
Line(69) = {6,10};
/* tole */
Line Loop(70) = {-25,28,27,32};
Plane Surface(71) = {70};
Line Loop(72) = {-25,-45,40,44};
Plane Surface(73) = {72};
Line Loop(74) = {40,-43,-42,39};
Plane Surface(75) = {74};
Line Loop(76) = {27,-46,-42,47};
Plane Surface(77) = {76};
Line Loop(78) = {-32,-46,43,44};
Plane Surface(79) = {78};
Line Loop(80) = {-47,39,45,28};
Plane Surface(81) = {80};
Surface Loop(82) = {79,71,-73,-81,-77,75};
Complex Volume(83) = {82};
/* boite air */
Line Loop(84) = {-13,-4,9,64};
Plane Surface(85) = {84};
Line Loop(86) = {-69,4,-68,5};
Plane Surface(87) = {86};
Line Loop(88) = {5,-12,62,16};
Plane Surface(89) = {88};
Line Loop(90) = {-63,-10,64,14};
Ruled Surface(91) = {90};
Line Loop(92) = {-62,-11,63,15};
Ruled Surface(93) = {92};
Line Loop(94) = {-69,9,10,11,12};
Plane Surface(95) = {94};
Line Loop(96) = {16,68,13,14,15};
Plane Surface(97) = {96};
Surface Loop(98) = {95,-87,-85,-97,89,93,91};
Complex Volume(99) = {98};
/* inducteur */
Line Loop(100) = {-67,-21,-2,17};
Plane Surface(101) = {100};
Line Loop(102) = {-66,-22,67,18};
Ruled Surface(103) = {102};
Line Loop(104) = {-65,-23,66,19};
Ruled Surface(105) = {104};
Line Loop(106) = {-20,-65,24,7};
Plane Surface(107) = {106};
Line Loop(108) = {7,8,5,6};
Plane Surface(109) = {108};
Line Loop(110) = {1,2,3,4};
Plane Surface(111) = {110};
Line Loop(112) = {14,15,16,-8,-20,-19,-18,-17,-1,13};
Plane Surface(113) = {112};
Line Loop(114) = {11,12,6,-24,-23,-22,-21,3,9,10};
Plane Surface(115) = {114};
Surface Loop(116) = {91,93,89,-109,107,-113,-105,115,-103,-101,-111,-85};
Complex Volume(117) = {116};
/* air */
Line Loop(118) = {-35,-34,-33,30};
Plane Surface(119) = {118};
Line Loop(120) = {-41,-43,-42,-38,34,36};
Plane Surface(121) = {120};
Line Loop(122) = {33,38,47,29};
Plane Surface(123) = {122};
Line Loop(124) = {32,26,-31,-30,-29,27};
Plane Surface(125) = {124};
Line Loop(126) = {31,-37,-36,35};
Plane Surface(127) = {126};
Line Loop(128) = {-37,-41,44,26};
Line Loop(129) = {8,68,1,2,3,69,6,7};
Plane Surface(130) = {128,129};
Surface Loop(131) = {97,-113,-130,127,125,79,-77,121,123,119,-101,-103,-105,107,115,-95};
Complex Volume(132) = {131};
Physical Volume(500) = { 99 , 131 }; /* boite air + air */
Physical Volume(501) = { 83 }; /* tole */
Physical Volume(501) = { 117 }; /* ind */
Physical Surface(1000) = { 87 , 109 , 111 , 119 , 121 , 123 , 125 , 127 , 130 };
Physical Surface(700) = { 109 };
Physical Surface(701) = { 111 };
This diff is collapsed.
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment