GENERATION_DE_MAILLAGE ---------------------- 0) lancer emc2 : taper la commande emc2 1) création du rectangle - menu haut/POINT - menu bas/XY PT/x (ou PT XY/x) - pour créer les 4 sommets du rectangle (-3,-1.2), (3,-1.2), (3,1.2), (-3,1.2), taper au clavier : -3=-1.2=3=-1.2=3=1.2=-3=1.2= - menu droite/SHOW_ALL (pour un zoom adapté à la figure) ou VOIR_TOUT - pour créer les extrémités (-3,-0.16) et (-3,0.64) du segment à travers lequel va entrer la condition aux limites, taper au clavier : -3=-0.16=-3=0.64= - menu haut/SEGMENT - menu bas/POINT/p : cliquer à la souris près d'un sommet défini précédemment, puis près du sommet suivant. Cela crée un segment joignant les deux points. Construire de même les autres côtés du rectangle. 2) création du disque - menu haut/ARC - menu haut/CENTER (ou CENTRE) - menu bas/XY PT/x (ou PT XY/x) - pour créer le centre du cercle de coordonnées (0,0), taper au clavier 0=0= - menu haut/RADIUS (ou RAYON) - pour définir le rayon, taper au clavier 0,39894228= - le cercle est alors tracé 3) préparation du maillage : découpage des frontières en intervalles - menu gauche/PREP_MESH - donnée du nombre d'intervalles : menu haut/NB_INTERVAL clavier : 30= cliquer près du cercle, et près du bord du bord droit du rectangle : ces 3 frontières sont découpées en 30 intervalles clavier : 70= cliquer près des bords haut puis bas du rectangle : ces 2 frontières sont découpées en 70 intervalles clavier : 12= cliquer près du segment inférieur du côté gauche clavier : 10= cliquer près du segment du milieu du côté gauche clavier : 8= cliquer près du segment supérieur du côté gauche - sauvegarde du travail effectué jusqu'ici (géométrie de la frontière et nb d'intervalles sur les bords) : cliquer menu gauche/SAVE (ou SAUVER) clavier : geometrie_cylindre, puis taper sur ENTREE le fichier geometrie_cylindre.emc2_bd est créé. 3bis) attribution d'un numéro de référence au bord d'entrée : pour repérer le bord sur lequel on imposera une CL d'entrée : - menu haut/REF_LIGNE - clavier : 10= - cliquer ensuite près du segment du milieu du bord gauche - ainsi, les noeuds de ce bord seront affectés d'une référence égale à 10 4) génération du maillage - menu gauche EDIT_MESH - taper ENTREE pour toutes les questions posées (on choisit ainsi les options par défaut) - la maillage s'affiche, et l'on voit que l'intérieur du disque a été maillé - pour supprimer l'intérieur du disque : menu bas/S_DOM/r menu haut/SUPPRESS (ou SUPPRIMER) cliquer dans le disque : la maillage est supprimé - sauvegarde du maillage : menu gauche/SAVE (ou SAUVER) taper au clavier le type de fichier de sortie (choisir am_fmt), taper ENTREE choisir et taper le nom donné au fichier (exemple maillage_cylindre) taper ENTREE - sortie du EMC2 : menu gauche/QUIT IMPORTANT : ----------- tout ceci est tiré de la doc, sections 2.4 et 7.1. À vous de vous familiariser avec ce logiciel. FORMAT DE SORTIE : ------------------ 1) données en sortie : nbt est le nombre de triangles nbs est le nombre de sommets nu(1:3,1:nbt) est un tableau d'entiers donnant les numéros des trois sommets de chaque triangle, dans le sens trigonométrique c(1:2,nbs) est un tableau de réels donnant les coordonnées des sommets refs(nbs) donnant les numéros de références des sommets référencés (pour repérer les noeuds de la frontière et les différentes CL à utiliser) reft(nbs) donnant les numéros de références des triangles référencés (pour repérer les triangles de la frontière et les différentes CL à utiliser) 2) lecture au format am_fmt : open(1,file=xxx.am_fmt,form=formatted,status=old) read (1,*) nbs,nbt read (1,*) ((nu(i,j),i=1,3),j=1,nbt) read (1,*) ((c(i,j),i=1,2),j=1,nbs) read (1,*) (reft(i),i=1,nbt) read (1,*) (refs(i),i=1,nbs) close(1)