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

Gaétan GREGOIRE g.gregoire.atlantis at groupeatlantis.com
Wed Nov 23 14:26:53 CET 2005


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 ?

 

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

-------------------------------------------------------------------------------------------

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20051123/172d07d9/attachment.html>