[Gmsh] (no subject)

David Colignon David.Colignon at ulg.ac.be
Thu Mar 4 07:59:45 CET 2010


Bonjour Mohammed,

tu dois définir des "Physical Surface"

http://geuz.org/gmsh/doc/texinfo/gmsh.html#Elementary-vs-physical-entities

http://geuz.org/gmsh/doc/texinfo/gmsh.html#Surfaces

et

http://geuz.org/gmsh/doc/texinfo/gmsh.html#t1_002egeo

Cordialement,

Dave

--
David Colignon, Ph.D.
Collaborateur Logistique du F.R.S.-FNRS
CÉCI - Consortium des Équipements de Calcul Intensif
ACE - Applied & Computational Electromagnetics
Institut Montefiore B28
Université de Liège
4000 Liège - BELGIQUE
Tél: +32 (0)4 366 37 32
Fax: +32 (0)4 366 29 10
WWW:    http://hpc.montefiore.ulg.ac.be/
Agenda: http://www.google.com/calendar/embed?src=david.colignon%40gmail.com

On 03/03/10 20:04, Mohammed Mobtil wrote:
> 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 <mailto:mobtil at ensm-douai.fr>
>
>
>
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh