[Gmsh] Décomposition d'un solide en deux entités - problème de maillage (French + english message)

Christophe Geuzaine christophe.geuzaine at case.edu
Wed Nov 23 17:14:21 CET 2005


Gaétan GREGOIRE wrote:
> Bonjour à tous / hi all.
> 
>  
> 
> J’ai une question à propos de la décomposition de solides en différentes 
> entités. Chacune des deux entités que j’ai défini peut être maillée 
> séparémment, mais ensemble beaucoup de messages d’erreurs apparaissent 
> et il n’y a pas de maillage.
> 
>  
> 
> Quelle commande me manque-t’il ?
> 

The problem is due to the fact that you have overlapping surfaces in the 
geometry (Gmsh does not compute intersections in the CAD).



>  
> 
> Merci d’avance pour votre patience.
> 
>  
> 
> Voici le contenu du fichier .geo :
> 
>  
> 
> /I have a question about decomposition of several physical entities. 
> Each physical entity in my example can be easily meshed without the 
> other part, but I can’t mesh these two entities together./
> 
> / /
> 
> /What command line is missing, please ?/
> 
> / /
> 
> /Thanks by advance for your help. I included the content of the .geo file:/
> 
>  
> 
>  
> 
> //////////////////////////////////////////////////////////////
> 
> // Maillage du panneau - GMSH version 1.2
> 
> //////////////////////////////////////////////////////////////
> 
>  
> 
> // Variables
> 
> //////////////////////////////////////////////////////////////
> 
>  
> 
> // Espace avant commencer
> 
> Rint= 0;
> 
>  
> 
> // Hauteur du panneau
> 
> H= 1.06;
> 
>  
> 
> // Largeur du panneau béton
> 
> L= 1;
> 
>  
> 
> // Epaisseur du panneau béton
> 
> e= 0.02;
> 
>  
> 
> // Débord
> 
> D= 0.078;
> 
>  
> 
> // hauteur fixation - jambe de force
> 
> Lf1= 0.760;
> 
>  
> 
> //Largeur de fixation - jambe de force
> 
> Lf2= 0.450;
> 
>  
> 
> // diamètre JdF
> 
> D2= 0.05;
> 
>  
> 
> // Epaisseur JdF
> 
> e2= 0.015;
> 
>  
> 
> // Distance fond du panneau/fixation JdF
> 
> H2=0.03;
> 
>  
> 
> // Modélisation des 2 structures
> 
> //////////////////////////////////////////////////////////////
> 
>  
> 
> // Points, panneau
> 
> Point(1) = {Rint,0,-L/2,1};
> 
> Point(2) = {Rint,0,L/2,1};
> 
> Point(3) = {Rint,H-e,L/2,1};
> 
> Point(4) = {Rint,H-e,-L/2,1};
> 
> Point(5) = {Rint+e,e,-L/2,1};
> 
> Point(6) = {Rint+e,e,L/2,1};
> 
> Point(7) = {Rint+e,H-e,L/2,1};
> 
> Point(8) = {Rint+e,H-e,-L/2,1};
> 
> Point(9) = {Rint+e,e,-e,1};
> 
> Point(10) = {Rint+e,e,e,1};
> 
> Point(11) = {Rint+e,H-e,e,1};
> 
> Point(12) = {Rint+e,H-e,-e,1};
> 
> Point(13) = {Rint+D,e,-e,1};
> 
> Point(14) = {Rint+D,e,e,1};
> 
> Point(15) = {Rint+D,H-e,e,1};
> 
> Point(16) = {Rint+D,H-e,-e,1};
> 
> Point(17) = {Rint+D,H-e,-L/2,1};
> 
> Point(18) = {Rint+D,H-e,L/2,1};
> 
> Point(19) = {Rint+D,H,-L/2,1};
> 
> Point(20) = {Rint+D,H,L/2,1};
> 
> Point(21) = {Rint,H,-L/2,1};
> 
> Point(22) = {Rint,H,L/2,1};
> 
> Point(23) = {Rint+D,0,-L/2,1};
> 
> Point(24) = {Rint+D,0,L/2,1};
> 
> Point(25) = {Rint+D,e,-L/2,1};
> 
> Point(26) = {Rint+D,e,L/2,1};
> 
>  
> 
> // Points, JdF
> 
> Point(27) = {Rint+H2+Lf2-D2/2,0,e,1};
> 
> Point(28) = {Rint+H2+Lf2-D2/2,0,e+e2,1};
> 
> Point(29) = {Rint+H2+Lf2+D2/2,0,e,1};
> 
> Point(30) = {Rint+H2+Lf2+D2/2-e2,0,e+e2,1};
> 
> Point(31) = {Rint+H2,Lf1-D2/2,e,1};
> 
> Point(32) = {Rint+H2,Lf1-D2/2,e+e2,1};
> 
> Point(33) = {Rint+H2,Lf1+D2/2,e,1};
> 
> Point(34) = {Rint+H2,Lf1+D2/2-e2,e+e2,1};
> 
> Point(35) = {D,Lf1+D2/2-Lf1*(D-H2)/Lf2,e,1};
> 
> Point(36) = {D,Lf1-D2/2-Lf1*(D-H2)/Lf2,e,1};
> 
> Point(37) = {Rint+H2,Lf1+D2/2-e2,e+D2,1};
> 
> Point(38) = {Rint+H2,Lf1+D2/2,e+D2,1};
> 
> Point(39) = {Rint+H2+Lf2+D2/2,0,e+D2,1};
> 
> Point(40) = {Rint+H2+Lf2+D2/2-e2,0,e+D2,1};
> 
>  
> 
> // Lignes, panneau
> 
> Line(1) = {1,23};
> 
> Line(2) = {5,9};
> 
> Line(3) = {24,2};
> 
> Line(4) = {2,1};
> 
> Line(5) = {1,21};
> 
> Line(6) = {21,22};
> 
> Line(7) = {22,20};
> 
> Line(8) = {20,18};
> 
> Line(9) = {18,7};
> 
> Line(10) = {7,6};
> 
> Line(11) = {19,17};
> 
> Line(12) = {17,8};
> 
> Line(13) = {8,5};
> 
> Line(14) = {6,10};
> 
> Line(15) = {10,14};
> 
> Line(17) = {13,9};
> 
> Line(18) = {13,16};
> 
> Line(19) = {16,12};
> 
> Line(20) = {12,9};
> 
> Line(21) = {14,15};
> 
> Line(22) = {15,11};
> 
> Line(23) = {2,22};
> 
> Line(24) = {10,11};
> 
> Line(25) = {21,19};
> 
> Line(26) = {8,12};
> 
> Line(27) = {7,11};
> 
> Line(28) = {17,16};
> 
> Line(29) = {19,20};
> 
> Line(30) = {23,24};
> 
> Line(31) = {23,25};
> 
> Line(32) = {25,13};
> 
> Line(33) = {24,26};
> 
> Line(34) = {6,26};
> 
> Line(35) = {5,25};
> 
> Line(37) = {15,18};
> 
> Line(39) = {14,26};
> 
>  
> 
> // Lignes, JdF
> 
> Line(40) = {27,29};
> 
> Line(42) = {30,28};
> 
> Line(43) = {28,27};
> 
> Line(44) = {33,31};
> 
> Line(46) = {32,34};
> 
> Line(50) = {30,34};
> 
> Line(51) = {28,32};
> 
> Line(52) = {33,29};
> 
> Line(53) = {31,27};
> 
> Line(56) = {31,32};
> 
> Line(57) = {34,37};
> 
> Line(58) = {37,38};
> 
> Line(59) = {38,33};
> 
> Line(60) = {29,39};
> 
> Line(61) = {39,40};
> 
> Line(62) = {40,30};
> 
> Line(64) = {38,39};
> 
> Line(65) = {37,40};
> 
>  
> 
> // Surfaces, panneau
> 
> Line Loop(78) = {24,-22,-21,-15};
> 
> Plane Surface(79) = {78};
> 
> Line Loop(80) = {17,-20,-19,-18};
> 
> Plane Surface(81) = {80};
> 
> Line Loop(82) = {22,-27,-9,-37};
> 
> Plane Surface(83) = {82};
> 
> Line Loop(84) = {7,8,9,10,34,-33,3,23};
> 
> Plane Surface(85) = {84};
> 
> Line Loop(86) = {14,15,39,-34};
> 
> Plane Surface(87) = {86};
> 
> Line Loop(88) = {17,-2,35,32};
> 
> Plane Surface(89) = {88};
> 
> Line Loop(90) = {1,31,-35,-13,-12,-11,-25,-5};
> 
> Plane Surface(91) = {90};
> 
> Line Loop(92) = {28,19,-26,-12};
> 
> Plane Surface(93) = {92};
> 
> Line Loop(94) = {29,-7,-6,25};
> 
> Plane Surface(95) = {94};
> 
> Line Loop(96) = {6,-23,4,5};
> 
> Plane Surface(97) = {96};
> 
> Line Loop(98) = {4,1,30,3};
> 
> Plane Surface(99) = {98};
> 
> Line Loop(100) = {21,37,-8,-29,11,28,-18,-32,-31,30,33,-39};
> 
> Plane Surface(101) = {100};
> 
> Line Loop(102) = {26,20,-2,-13};
> 
> Plane Surface(103) = {102};
> 
> Line Loop(104) = {14,24,-27,10};
> 
> Plane Surface(105) = {104};
> 
>  
> 
> // Surfaces, JdF
> 
> Line Loop(106) = {59,44,56,46,57,58};
> 
> Plane Surface(107) = {106};
> 
> Line Loop(108) = {52,60,-64,59};
> 
> Plane Surface(109) = {108};
> 
> Line Loop(110) = {62,50,57,65};
> 
> Plane Surface(111) = {110};
> 
> Line Loop(112) = {56,-51,43,-53};
> 
> Plane Surface(113) = {112};
> 
> Line Loop(114) = {65,-61,-64,-58};
> 
> Plane Surface(115) = {114};
> 
> Line Loop(116) = {46,-50,42,51};
> 
> Plane Surface(117) = {116};
> 
> Line Loop(118) = {40,-52,44,53};
> 
> Plane Surface(119) = {118};
> 
> Line Loop(120) = {62,42,43,40,60,61};
> 
> Plane Surface(121) = {120};
> 
>  
> 
> // Volume, panneau
> 
> Surface Loop(200) = {101,79,-105,87,85,95,97,-99,91,89,-81,-103,-93,83};
> 
> Volume(201) = {200};
> 
>  
> 
> // Volume, JdF
> 
> Surface Loop(202) = {111,-121,117,-107,109,119,113,-115};
> 
> Volume(203) = {202};
> 
>  
> 
> // Sauvegarde des entités des 2 structures
> 
> //////////////////////////////////////////////////////////////
> 
>  
> 
> // Panneau
> 
> Physical Surface(1000) = {79};
> 
> Physical Surface(1001) = {81};
> 
> Physical Surface(1002) = {83};
> 
> Physical Surface(1003) = {85};
> 
> Physical Surface(1004) = {87};
> 
> Physical Surface(1005) = {89};
> 
> Physical Surface(1006) = {91};
> 
> Physical Surface(1007) = {93};
> 
> Physical Surface(1008) = {95};
> 
> Physical Surface(1009) = {97};
> 
> Physical Surface(1010) = {99};
> 
> Physical Surface(1011) = {101};
> 
> Physical Surface(1012) = {103};
> 
> Physical Surface(1013) = {105};
> 
> Physical Volume(2000) = {201};
> 
>  
> 
> // JdF
> 
> Physical Surface(1100) = {107};
> 
> Physical Surface(1101) = {109};
> 
> Physical Surface(1102) = {111};
> 
> Physical Surface(1103) = {113};
> 
> Physical Surface(1104) = {115};
> 
> Physical Surface(1105) = {117};
> 
> Physical Surface(1106) = {119};
> 
> Physical Surface(1107) = {121};
> 
> Physical Volume(2100) = {203};
> 
>  
> 
> // Ensemble
> 
> Physical Volume(2200) = {201,203};
> 
>  
> 
> // Caractéristiques du mesh
> 
> //////////////////////////////////////////////////////////////
> 
> Mesh.CharacteristicLengthFactor= 0.05;
> 
>  
> 
>  
> 
>  
> 
> **Gaétan GREGOIRE**
> 
> Responsable Recherche & Développement
> 
> -------------------------------------------------------------------------------------------
> 
> *Groupe ATLANTIS*
> 
> /Parc d'activité du Baconnet/
> 
> /69700 Montagny - France/
> 
> Tél/fax: +33 (0) 4 72 49 17 42
> 
> Mél: g.gregoire.atlantis at groupeatlantis.com 
> <mailto:g.gregoire.atlantis at groupeatlantis.com>
> 
> -------------------------------------------------------------------------------------------
> 
>  
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh


-- 
Christophe Geuzaine
Assistant Professor, Case Western Reserve University, Mathematics
http://www.case.edu/artsci/math/geuzaine