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] <mailto:[email protected]>
_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh
