[Gmsh] (no subject)

[Mohammed Mobtil]

Bonjour,

 

J'utilise sur Gmsh une géométrie 3D comportant des cercle (demi cercle) dont les points centre des ces cercle n'appartient pas à mon domaine de calcul.

Le but est de récupérer la matrice des coordonnées ainsi que la matrice de connectivité (les éléments triangulaire) pour résoudre mon problème physique par éléments finis.

 

Or le fichier que je récupère contient les coordonnées des centres des cercles qui ne font pas partis de mon domaine de calcul mais qui sont un passage obligé pour dessiner des arcs.

 

Auriez-vous une suggestion pour générer des cercles sans utiliser les centres?? Ou comment je peux récupérer les données du maillage (coordonnées + connectivité) de mon domaine de calcul sans les informations des points centre?

 

Exemple: voir le scripte au dessous.

 

Ci joint le fichier (.msh) que je récupère: 

Dans la matrice des coordonnées, on trouve aussi les 2 centres (points 3 et 8) qui ne font pas partis de mon domaine.

Même si je peux les supprimer du fichier, il y aura un problème de décalage dans la matrice de connectivité (les sommets des éléments ne seront plus représentatif de ma géométrie)


Bien cordialement,



%%%%%%%%%%%%%%%
c = 0.1;
 
 //points suivant x (pour y=z=0)
 Point(1) = {0,0,0,lc};
 Point(2) = {4,0,0,lc};
 Point(3) = {5,0,0,lc};
 Point(4) = {6,0,0,lc};
 Point(5) = {10,0,0,lc};
 
 //Point suivant x (pour y=3 et z=0)
 Point(6) = {0,3,0,lc};
 Point(7) = {4,3,0,lc};
 Point(8) = {5,3,0,lc};
 Point(9) = {6,3,0,lc};
 Point(10) = {10,3,0,lc};
 
 //Point suivant z (pour x=5 et y=0&y=3)
 Point(11) = {5,0,1,lc};
 Point(12) = {5,3,1,lc};
 
  //LIGNES
 Line(1)   = {1,2};
 Line(2)   = {4,5};
 Line(3)   = {6,7};
 Line(4)   = {9,10};
 Line(5)   = {1,6};
 Line(6)   = {2,7};
 Line(7)   = {11,12};
 Line(8)   = {4,9};
 Line(9)   = {5,10};
 
 //ARCS
 Circle(10) = {2,3,11};
 Circle(11) = {11,3,4};
 Circle(12) = {7,8,12};
 Circle(13) = {12,8,9};
 
 //SURFACE
 Line Loop(1) = {1,6,-3,-5};
 Ruled Surface(1) = {1};
 
 Line Loop(2) = {10,7,-12,-6};
 Ruled Surface(2) = {2};
 
 Line Loop(3) = {11,8,-13,-7};
 Ruled Surface(3) = {3};
 
 Line Loop(4) = {2,9,-4,-8};
 Ruled Surface(4) = {4};
 
  //MAILLAGE
 Transfinite Line {1,3,2,4} = 4 Using Progression 1.0;
 Transfinite Line {10,12,11,13} = 5 Using Progression 1.0;
 Transfinite Line {5,6,7,8,9} = 6 Using Progression 1.0;
 
 Transfinite Surface {1} = {1,2,7,6};  
 Transfinite Surface {2} = {2,11,12,7}; 
 Transfinite Surface {3} = {11,4,9,12};  
 Transfinite Surface {4} = {4,5,10,9}; 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


----------------------------------------------------------------
Mohammed MOBTIL
Doctorant
Dépt  Energétique Industrielle -  Ecole des Mines
941, rue Charles BOURSEUL
BP 10838 - F59508 DOUAI CEDEX
Tel : 33 (0)3 27 71 23 96  -  Fax : 33 (0)3 27 71 29 15
mobtil at ensm-douai.fr
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20100303/f495ac12/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Text123.geo
Type: application/octet-stream
Size: 1313 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20100303/f495ac12/attachment.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Text123.msh
Type: application/octet-stream
Size: 7813 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20100303/f495ac12/attachment.msh>