[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>