Suiv.: 1.14 La SD TAE Sup.: 1 Description des SD par type Préc.: 1.12 La SD NDL1 Index Table des matières

1.13 La SD NOPO

1.13.1 Contenu

Cette SD est l'organisation choisie pour stocker les maillages bidimensionnels ou tridimensionnels. Elle correspond à une description du maillage, élément par élément. Chaque élément est fourni via les informations suivantes :

• son type géométrique,
• son numéro de sous-domaine,
• la liste de ses sommets,
• la liste de ses noeuds,
• les attributs (numéro de référence) de ses sommets, arêtes et faces.

La SD NOPO est composée de six tableaux dont l'organisation est prédéfinie.

Tableau NOP0 : Informations générales.

De type entier, ce tableau possède 32 variables. Il contient une description générale du travail (titre, date, nom), de la SD NOPO (type, niveau, ...) et indique la présence ou non du tableau NOP1.

1:20 TITRE

le titre du travail en 20 mots de 4 caractères,

21:22 DATE

la date de création en 2 mots de 4 caractères,

23:28 NOMCRE

le nom du créateur en 6 mots de 4 caractères,

29 'NOPO'

le type de la SD,

30 NIVEAU

le paramètre de niveau de la SD,

31 ETAT

un paramètre réservé,

32 NTACM

le nombre de tableaux supplémentaires associés à la SD
(ils seront décrits dans le tableau NOP1).

Tableau NOP1 : Descripteur des éventuels tableaux supplémentaires.

Ce tableau est analogue au tableau B1 de la SD B (voir cette SD).

Tableau NOP2 : Description générale du maillage.

Ce tableau de type entier contient 27 valeurs.

1 NDIM

la dimension de l'espace (2 ou 3),

2 NDSR

le maximum des numéros de référence,

3 NDSD

le maximum des numéros de sous-domaine,

4 NCOPNP

le code de coïncidence des noeuds et des sommets : 1 si il y a coïncidence partout, 0 sinon,

5 NE

le nombre d'éléments du maillage,

6 NEPO

le nombre d'éléments réduits à un point,

7 NSEG

le nombre de segments,

8 NTRI

le nombre de triangles,

9 NQUA

le nombre de quadrangles,

10 NTET

le nombre de tétraèdres,

11 NPEN

le nombre de pentaèdres,

12 NHEX

le nombre d'hexaèdres,

13 NSUP

le nombre de super-éléments,

14 NEF

le nombre d'éléments frontaliers,

15 NOE

le nombre de noeuds,

16 N1

le nombre de noeuds sur un segment ou une arête (hors extrémités),

17 ISET

le nombre de noeuds internes à un triangle ou une face triangulaire,

18 ISEQ

le nombre de noeuds internes à un quadrangle ou une face quadrangulaire,

19 ISETE

le nombre de noeuds internes à un tétraèdre,

20 ISEPE

le nombre de noeuds internes à un pentaèdre,

21 ISEHE

le nombre de noeuds internes à un hexaèdre,

22 NP

le nombre de points,

23 NTYCOO

le type des valeurs des coordonnées (2 ici, cf. NOP1),

24 LPGDN

la plus grande différence entre les numéros des noeuds d'un élément +1,

25 NBEGM

le nombre de super-éléments ou de descriptions dans le tableau NOP3,

26 LNOP5

le nombre de mots du tableau NOP5,

27 NTACOO

le type des axes de coordonnées : 1 x,y,z ; 2 ; 3 .

Tableau NOP3 : Pointeur éventuel.

Si NBEGM (tableau NOP2) est non nul, ce tableau de type entier contient sur une variable une information relative à chaque super-élément ou description (ce tableau n'est pas actuellement utilisé).

Tableau NOP4 : Coordonnées des sommets.

De type NTYCOO ce tableau contient les NDIM.NP coordonnées des sommets :

X1, Y1, X2, Y2, ... ou X1, Y1, Z1, X2, ...

Tableau NOP5 : Description séquentielle des éléments.

De type entier, ce tableau décrit séquentiellement chaque élément du maillage, il contient :

Boucle i de 1 à NE

• NCGE, le code de la géométrie de l'élément (1 = noeud, 2 = segment, 3 = triangle, 4 = quadrangle, 5 = tétraèdre, 6 = pentaèdre, 7 = hexaèdre, 8 = super-élément, ...)
• NMAE, le nombre de mots nécessaires au stockage des numéros de référence des faces, des arêtes et des sommets,
• NDSDE, le numéro de sous-domaine de l'élément,
• NNO, le nombre de noeuds de l'élément,

  Boucle j de 1 à NNO

  • NONO(j), le numéro du noeud j de l'élément,
  Fin de la boucle j.
• si NCOPNP (cf. tableau NOP2) est nul :
  • NPO, le nombre de points de l'élément,

    Boucle j de 1 à NPO

    1. NOPO(j), le numéro du point j de l'élément,
    Fin de la boucle j.
• si NMAE est non nul :
  • INING, l'indicateur de niveau des informations sur les numéros de référence. Si INING = 3 seuls les numéros de référence des sommets sont donnés, pour INING = 2 les numéros de référence des arêtes puis des sommets sont donnés et pour INING = 1 seront présents les numéros des faces, des arêtes puis des sommets.

    Boucle de 2 à NMAE

    1. si INING = 3, les références des sommets,
    2. si INING = 2, les références des arêtes puis des sommets,
    3. si INING = 1, les références des faces, des arêtes puis des sommets.
    Fin de la boucle.

Définitions relatives à un élément

Chaque élément est défini, le premier sommet étant choisi, par :

 
Figure 1.4: Segment, triangle et quadrangle 


1. Le segment  :
   • (1) (2)

2. Le triangle  : numérotation dans le sens direct (en dimension 2 ; dans le cas de la dimension 3, la surface n'est pas signée).
   • arête [1] : (1) (2) (noeuds 4, 5, ... s'il y a des noeuds supplémentaires)
   • arête [2] : (2) (3) (...)
   • arête [3] : (3) (1)

3. Le quadrangle  : numérotation dans le sens direct (même remarque que pour le triangle)
   • arête [1] : (1) (2) (noeuds 5, 6, ...)
   • arête [2] : (2) (3) (...)
   • arête [3] : (3) (4)
   • arête [4] : (4) (1)

    
   Figure: Tétraèdre, pentaèdre et hexaèdre 
   

4. Le tétraèdre  : trièdre direct. ( désigne le vecteur d'origine le point i et d'extrémité le point j).
   • arête [1] : (1) (2) (noeuds 5, ...)
   • arête [2] : (2) (3) (...)
   • arête [3] : (3) (1)
   • arête [4] : (1) (4)
   • arête [5] : (2) (4)
   • arête [6] : (3) (4) (noeuds ..., p)

   Toute face, vue de l'extérieur, est dans le sens direct :

   • face [1] : (1) (3) (2) (noeuds p+1, ...)
   • face [2] : (1) (4) (3) (...)
   • face [3] : (1) (2) (4)
   • face [4] : (2) (3) (4)

5. Le pentaèdre  : trièdre direct
   • arête [1] : (1) (2) (noeuds 7, ...)
   • arête [2] : (2) (3) (...)
   • arête [3] : (3) (1)
   • arête [4] : (1) (4)
   • arête [5] : (2) (5)
   • arête [6] : (3) (6)
   • arête [7] : (4) (5)
   • arête [8] : (5) (6)
   • arête [9] : (6) (4) (noeuds ..., p)

   Toute face, vue de l'extérieur, est dans le sens direct.

   • face [1] : (1) (3) (2) (noeuds p+1, ...)
   • face [2] : (1) (4) (6) (3)
   • face [3] : (1) (2) (5) (4)
   • face [4] : (4) (5) (6)
   • face [5] : (2) (3) (5) (6)

6. L'hexaèdre  : trièdre direct
   • arête [1] : (1) (2) (noeuds 9, ...)
   • arête [2] : (2) (3) (...)
   • arête [3] : (3) (4)
   • arête [4] : (4) (1)
   • arête [5] : (1) (5)
   • arête [6] : (2) (6)
   • arête [7] : (3) (7)
   • arête [8] : (4) (8)
   • arête [9] : (5) (6)
   • arête [10] : (6) (7)
   • arête [11] : (7) (8)
   • arête [12] : (8) (5) (noeuds ..., p)

   Toute face, vue de l'extérieur, est dans le sens direct.

   • face [1] : (1) (4) (3) (2) (noeuds p+1, ...)
   • face [2] : (1) (5) (8) (4)
   • face [3] : (1) (2) (6) (5)
   • face [4] : (5) (6) (7) (8)
   • face [5] : (2) (3) (7) (6)
   • face [6] : (3) (4) (8) (7)

Remarque : A ces éléments correspondant aux formes géométriques habituelles s'ajoutent d'autres éléments développés pour traiter des problèmes spécifiques ; on peut citer, par exemple, les éléments de raccord  permettant la simulation d'une charnière  (raccord entre deux poutres , soit deux segments ; entre deux plans, soit deux quadrangles dans l'espace).

Remarque : Au fur et à mesure des besoins, de nouveaux éléments pourront apparaître. Il faudra alors les définir de manière consistante.

Outils associés

Ecriture et lecture d'une SD NOPO

Une SD, résidant en mémoire centrale, est écrite en totalité (elle est de catégorie 1) sur fichier séquentiel via le module SDECRI.

De même, une SD, résidant sur un fichier, est lue et reportée en totalité en mémoire centrale via le module SDLECT.

Impression d'une SD NOPO et dessin d'un maillage

Le contenu (total ou partiel) d'une SD NOPO peut être imprimé en utilisant les modules suivants :

• IMNOPO : la SD en entier (selon le taux d'impression IMPRE),

      SUBROUTINE  IMNOPO(M,NFNOPO,NINOPO,IMPRE)
C BUT : IMRESSION DES TABLEAUX DE LA SD NOPO, CROISSANTE AVEC IMPRE
C ---   FICHIER NFNOPO (OU MC SI 0), NIVEAU NINOPO

• IMNOPS : le tableau NOP2 et certains éléments du maillage,

      SUBROUTINE IMNOPS(M,NFNOPO,NINOPO,ISTRES,INDIC,LINDIC)
C BUT : IMRESSION DES TABLEAUX 2 ET 5 DE LA SD NOPO
C ---   FICHIER NFNOPO ( OU MC SI 0 ), NIVEAU NINOPO
C       ISTRES = 0 PAS DE NOP2, SINON IMPRESSION DE NOP2
C       INDIC(*) LES NUMEROS DES LINDIC ELEMENTS DE NOP5 A IMPRIMER

• INFONO : telle ou telle valeur contenue dans la SD

      SUBROUTINE INFONO(M,IOPT,TYPESD,NIVEAU,TAB2,TAB4,TAB5,NUMREC,
     +                  XREC,YREC,ZREC)
C BUT : DONNER DES INFORMATIONS SUR UNE SD DE MAILLAGE ( NOPO )

Les préprocesseurs IMAGXX et INFOXX permettent l'appel conversationnel de IMNOPO et INFONO respectivement.

Le maillage contenu dans un fichier peut être dessiné via le préprocesseur TRNOXX [RM 96].

Modules de création et de manipulation d'une SD NOPO

On se reportera à [Guide Modulef - 3] pour les méthodes et les modules de création et de modification de maillages.

Utilitaires de manipulation d'une SD NOPO

Il existe de nombreux utilitaires permettant de manipuler le contenu d'une SD NOPO, parmi ceux-ci, on trouve des utilitaires facilitant l'exploration du tableau NOP5.

 
      SUBROUTINE CLNOPO (M,XM,IANOP2,IANOP4,IA,NCGE,NDSDE,NNO,NONO,
     +                   NPO,NOPOI,COORP,NREF)
C BUT : LIRE LE TABLEAU NOP5  ET LE METTRE EN CLAIR (!)
C ---
C IN :
C --
C M,XM    : LE SUPER TABLEAU
C IA      : ADRESSE DU DEBUT DU TABLEAU NOP5  ( IANOP5 )
C OUT :
C ---
C MEMES MNEMONIQUES QUE DANS LA BROCHURE < DESCRIPTION DES SD >
C IA      : ACTUALISE ( DEBUT DE L'ELEMENT SUIVANT )
C COORP   : COORDONNEES DES POINTS
C NREF    : TABLEAU DES NUMEROS DE REFERENCE
C              3-D : POINTS,ARETES,FACES
C              2-D : POINTS,ARETES


      SUBROUTINE CLNOIA(NCOPNP,NOP5,IA,NCGE,IAS,IAR,IAF,NPO,NOPO)
C BUT : CLEAR LE TABLEAU NOP5 (CF. CLNOPO)
C ---
C OUT : IA, IAS, IAR, IAF, NOPO ET NPO
C ---   ADRESSES DES REFERENCES OU 0, TABLEAU DES POINTS, LEUR NOMBRE 


      SUBROUTINE CLNOIN(NCOPNP,NOP5,IA,NCGE,IAS,IAR,IAF,NPO,NNO,NONO)
C BUT : CLEAR LE TABLEAU NOP5 ( CF. CLNOPO )
C --- 
C OUT :  IA, IAS, IAR, IAF, NPO, NNO, NONO  


      SUBROUTINE CLNO5C(NDIM,NCOPNP,NOP5,IA,NCGE,NDSDE,NNO,NONO,
     +                  NPO,NOPOI,COOR4,COORP)
C BUT : LIRE NOP5 ET COOR4 ( ALIAS NOP4 )  POUR REMPLIR NONO,NOPOI,...

Les utilitaires NBSOME, NBARET et NBFACE donnent le nombre de sommets, arêtes et faces de chaque type d'élément.