[Gmsh] Problème de pivot nul sur le logiciel de FEM

Christophe Geuzaine cag32 at case.edu
Thu Nov 24 19:16:13 CET 2005


Gaétan GREGOIRE wrote:
> Bonjour à tous.
> 
>  
> 
> J’ai retravaillé mon fichier pour extruder le maillage, cette fois ci 
> plus de problèmes de points coplanaires (merci pour vos précieux conseils !)
> 
>  
> 
> Cependant, lorsque mon maillage est fini on peut voir 5 couleurs 
> différentes, correspondants aux 5 volumes élémentaires. Sous le logiciel 
> de FEM que j’utilise (Code_Aster), un problème de pivot nul est détecté.
> 
> Est-ce que cela veut dire que mes volumes ne sont pas reliés ? Où est le 
> problème SVP ?
> 

You still have two overlapping surfaces (324 and 231), with 
non-compatible meshes: cf. attached image. This is probably the root of 
your problem.



>  
> 
> Merci d’avance pour vos conseils…
> 
>  
> 
> Voici le fichier .Geo :
> 
>  
> 
> //////////////////////////////////////////////////////////////
> 
> // Maillage du panneau - GMSH version 1.2
> 
> //////////////////////////////////////////////////////////////
> 
> // Variables
> 
> //////////////////////////////////////////////////////////////
> 
> // Espace avant commencer
> 
> Rint= 0;
> 
> // Hauteur du panneau
> 
> H= 1060;
> 
> // Largeur du panneau béton
> 
> L= 1000;
> 
> // Epaisseur du panneau béton
> 
> e= 2;
> 
> // Débord
> 
> D= 78;
> 
> // hauteur fixation - jambe de force
> 
> Lf1= 760;
> 
> //Largeur de fixation - jambe de force
> 
> Lf2= 450;
> 
> // diamètre JdF
> 
> D2= 50;
> 
> // Epaisseur JdF
> 
> e2= 1.5;
> 
> // Distance fond du panneau/fixation JdF
> 
> H2= 30;
> 
> // Taille de la maille
> 
> LgCarra = 50;
> 
> // Points
> 
> //////////////////////////////////////////////////////////////
> 
> Point(1) = {Rint,0,0,1};
> 
> Point(2) = {Rint+D,0,0,1};
> 
> Point(3) = {Rint+D,e,0,1};
> 
> Point(4) = {Rint+e,e,0,1};
> 
> Point(5) = {Rint+e,H-e,0,1};
> 
> Point(6) = {Rint+D,H-e,0,1};
> 
> Point(7) = {Rint+D,H,0,1};
> 
> Point(8) = {Rint,H,0,1};
> 
> Point(9) = {Rint+e,e,0,1};
> 
> Point(10) = {Rint+D,e,0,1};
> 
> Point(11) = {Rint+D,H-e,0,1};
> 
> Point(12) = {Rint+e,H-e,0,1};
> 
> Point(13) = {Rint+Lf2+H2,0,e,1};
> 
> Point(14) = {Rint+Lf2+H2,D2,e,1};
> 
> Point(15) = {Rint+Lf2+H2,D2,e+e2,1};
> 
> Point(16) = {Rint+Lf2+H2,e2,e+e2,1};
> 
> Point(17) = {Rint+Lf2+H2,e2,e+D2,1};
> 
> Point(18) = {Rint+Lf2+H2,0,e+D2,1};
> 
> Line(1) = {1,2};
> 
> Line(2) = {2,3};
> 
> Line(3) = {3,4};
> 
> Line(4) = {4,5};
> 
> Line(5) = {5,6};
> 
> Line(6) = {6,7};
> 
> Line(7) = {7,8};
> 
> Line(8) = {8,1};
> 
> Line(9) = {9,10};
> 
> Line(10) = {10,11};
> 
> Line(11) = {11,12};
> 
> Line(12) = {12,9};
> 
> Line(13) = {13,14};
> 
> Line(14) = {14,15};
> 
> Line(15) = {15,16};
> 
> Line(17) = {16,17};
> 
> Line(18) = {17,18};
> 
> Line(19) = {18,13};
> 
> Line Loop(100) = {7,8,1,2,3,4,5,6};
> 
> Plane Surface(101) = {100};
> 
> Line Loop(102) = {10,11,12,9};
> 
> Plane Surface(103) = {102};
> 
> Extrude Surface {101, 
> {0.0,0.0,-L/2}}{Layers{Ceil(L/2/LgCarra),1};Recombine;};
> 
> Extrude Surface {101, 
> {0.0,0.0,L/2}}{Layers{Ceil(L/2/LgCarra),1};Recombine;};
> 
> Extrude Surface {103, {0.0,0.0,-e}}{Layers{Ceil(e/2/LgCarra),1};Recombine;};
> 
> Extrude Surface {103, {0.0,0.0,e}}{Layers{Ceil(e/2/LgCarra),1};Recombine;};
> 
> Line Loop(300) = {17,18,19,13,14,15};
> 
> Plane Surface(301) = {300};
> 
> Angle = Atan(Lf1/Lf2);
> 
> Rotate { {0.0,0,-1},{Rint+Lf2+H2,0,e,1},Angle } {Surface{301};}
> 
> Extrude Surface {301, 
> {-Lf2,Lf1,0.0}}{Layers{Ceil(Sqrt(Lf1^2+Lf2^2)/LgCarra),1};Recombine;};
> 
> Physical Surface(1000) = {124,166};
> 
> Physical Surface(1100) = {120,162};
> 
> Physical Volume(2000) = {1,2,3,4,5};
> 
> Mesh.CharacteristicLengthFactor= LgCarra;
> 
> Coherence;
> 
>  
> 
>  
> 
> **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
-------------- next part --------------
A non-text attachment was scrubbed...
Name: overlap.png
Type: image/png
Size: 7978 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20051124/d0a76802/attachment.png>