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
[email protected]

Attachment: Text123.geo
Description: Binary data

Attachment: Text123.msh
Description: Binary data

_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh

Reply via email to