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