[Gmsh] Epaisseur de paroi et maillage (message in english and french)
Christophe Geuzaine
christophe.geuzaine at case.edu
Wed Nov 23 17:17:35 CET 2005
Gaétan GREGOIRE wrote:
> Bonjour à tous.
>
>
>
> Lorsque j’utilise ce fichier, le maillage ne marche pas en dessous d’une
> certaine épaisseur de paroi :
>
>
>
> · Si e=0.02 et e2=0.015, le maillage marche très bien
>
> · Si e=0.002 et e2=0.0015, GMSH plante et se ferme.
>
>
>
> Quelqu’un peut-il m’expliquer ?
>
The automatic 3D algorithm is designed to generate isotropic elements.
For structures as thin as those in your example, it is probably better
to use extruded meshes.
>
>
> Merci d’avance. Cordialement,
>
>
>
> Hi all.
>
>
>
> When I use this geometry, I’m unable to mesh the volume with some
> parameter values:
>
>
>
> · With e=0.02 and e2=0.015, everything is ok
>
> · With e=0.002 and e2=0.0015, GMSH crash and close itself.
>
>
>
> Can anybody explain that to me?
>
>
>
> Thank you by advance. Best regards,
>
>
>
>
>
> //////////////////////////////////////////////////////////////
>
> // Maillage du panneau + jb force - 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.002;
>
>
>
> // 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.0015;
>
>
>
> // Distance fond du panneau/fixation JdF
>
> H2=0.03;
>
>
>
> // Points
>
> //////////////////////////////////////////////////////////////
>
>
>
> 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};
>
> 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
>
> //////////////////////////////////////////////////////////////
>
>
>
> 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,36};
>
> 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};
>
> Line(40) = {27,29};
>
> Line(42) = {30,28};
>
> Line(43) = {28,27};
>
> Line(44) = {33,31};
>
> Line(46) = {32,34};
>
> Line(48) = {29,35};
>
> Line(49) = {27,36};
>
> Line(50) = {30,34};
>
> Line(51) = {28,32};
>
> Line(52) = {33,35};
>
> Line(53) = {31,36};
>
> Line(54) = {36,35};
>
> Line(55) = {35,15};
>
> 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
>
> //////////////////////////////////////////////////////////////
>
>
>
> Line Loop(66) = {37,9,27,-22};
>
> Plane Surface(67) = {66};
>
> Line Loop(68) = {7,8,9,10,34,-33,3,23};
>
> Plane Surface(69) = {68};
>
> Line Loop(70) = {4,5,6,-23};
>
> Plane Surface(71) = {70};
>
> Line Loop(72) = {10,14,24,-27};
>
> Plane Surface(73) = {72};
>
> Line Loop(74) = {6,7,-29,-25};
>
> Plane Surface(75) = {74};
>
> Line Loop(76) = {26,-19,-28,12};
>
> Plane Surface(77) = {76};
>
> Line Loop(78) = {25,11,12,13,35,-31,-1,5};
>
> Plane Surface(79) = {78};
>
> Line Loop(80) = {2,-17,-32,-35};
>
> Plane Surface(81) = {80};
>
> Line Loop(82) = {15,39,-34,14};
>
> Plane Surface(83) = {82};
>
> Line Loop(84) = {18,19,20,-17};
>
> Plane Surface(85) = {84};
>
> Line Loop(86) = {8,-37,-55,-54,-21,39,-33,-30,31,32,18,-28,-11,29};
>
> Plane Surface(87) = {86};
>
> Line Loop(88) = {2,-20,-26,13};
>
> Plane Surface(89) = {88};
>
> Line Loop(90) = {22,-24,15,21,-53,-44,52,55};
>
> Plane Surface(91) = {90};
>
> Line Loop(92) = {4,1,30,3};
>
> Plane Surface(93) = {92};
>
> Line Loop(94) = {59,44,56,46,57,58};
>
> Plane Surface(95) = {94};
>
> Line Loop(96) = {42,43,40,60,61,62};
>
> Plane Surface(97) = {96};
>
> Line Loop(98) = {48,-54,-49,40};
>
> Plane Surface(99) = {98};
>
> Line Loop(100) = {52,-48,60,-64,59};
>
> Plane Surface(101) = {100};
>
> Line Loop(102) = {61,-65,58,64};
>
> Plane Surface(103) = {102};
>
> Line Loop(104) = {57,65,62,50};
>
> Plane Surface(105) = {104};
>
> Line Loop(106) = {42,51,46,-50};
>
> Plane Surface(107) = {106};
>
> Line Loop(108) = {56,-51,43,49,-53};
>
> Plane Surface(109) = {108};
>
>
>
> // Volume
>
> //////////////////////////////////////////////////////////////
>
>
>
> Surface Loop(110) =
> {93,-71,79,75,-69,87,67,73,-83,91,-109,95,-101,-99,97,-107,-105,-103,81,-89,-85,-77};
>
> Volume(111) = {110};
>
>
>
> // Svg des entités
>
> //////////////////////////////////////////////////////////////
>
>
>
> Physical Surface(1000) = {67};
>
> Physical Surface(1001) = {69};
>
> Physical Surface(1002) = {71};
>
> Physical Surface(1003) = {73};
>
> Physical Surface(1004) = {75};
>
> Physical Surface(1005) = {77};
>
> Physical Surface(1006) = {79};
>
> Physical Surface(1007) = {81};
>
> Physical Surface(1008) = {83};
>
> Physical Surface(1009) = {85};
>
> Physical Surface(1010) = {87};
>
> Physical Surface(1011) = {89};
>
> Physical Surface(1012) = {91};
>
> Physical Surface(1013) = {93};
>
> Physical Surface(1014) = {95};
>
> Physical Surface(1015) = {97};
>
> Physical Surface(1016) = {99};
>
> Physical Surface(1017) = {101};
>
> Physical Surface(1018) = {103};
>
> Physical Surface(1019) = {105};
>
> Physical Surface(1020) = {107};
>
> Physical Surface(1021) = {109};
>
> Physical Volume(2000) = {111};
>
>
>
> // Génération du maillage
>
> //////////////////////////////////////////////////////////////
>
> Mesh.CharacteristicLengthFactor= 0.02;
>
>
>
>
>
>
>
> **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