[Gmsh] Problème de pivot nul sur le logiciel de FEM
Gaétan GREGOIRE
g.gregoire.atlantis at groupeatlantis.com
Thu Nov 24 18:58:35 CET 2005
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 ?
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
-------------------------------------------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20051124/28e8c292/attachment.html>