Suiv.: Appel des modules dérivés Sup.: Méthode de blocs Préc.: Description rapide de la méthode Index Table des matières

3.2.2 Appel du module COLIB2

Le préprocesseur COLIXX  (cf. plus bas et voir également les versions simplifiées et le module COLIBR) permet l'appel conversationnel du module dans le cas où le nombre  de données à fournir est peu important. Dans le cas général, il parait plus simple d'appeler directement le module en écrivant un programme principal comme suit:

• dimensionner le tableau M (dans le common blanc)
• déclarer le tableau entier NARET(5,NA) avec NA le nombre d'arêtes grossières
• déclarer le tableau entier NTGRO(11,NBGRO) avec NBGRO le nombre d'éléments grossiers
• déclarer le tableau entier NTRS(NBS) avec NBS le nombre de sommets grossiers
• déclarer le tableau entier NTFR(6,NBFR) avec NBFR le nombre de faces grossières ayant un numéro de référence non nul (si NBFR = 0, on déclarera NTFR comme NTFR(6,1))
• déclarer le tableau réel XYZ(3,NBS) avec NBS le nombre de sommets grossiers
• déclarer le tableau réel XYZINT(3,LINT) avec LINT le nombre total de points de contrôle donnés par l'utilisateur sur les arêtes grossières (si LINT = 0, on déclarera XYZINT comme XYZINT(3,1))
• appeler INITI

appel du module

CALL COLIB2 (M,NFNOPO,NINOPO,NDIMS,NBS,NA,NBGRO,
. NBFR,XYZ,NTRS,NARET,XYZINT,NTFR,NTGRO)

avec :

• M : le super-tableau
• NF(NI)NOPO : le numéro de support de la S.D.S. NOPO et son niveau
• NDIMS : la dimension de l'espace de sortie, si NDIMS = 2, la cote Z est ignorée (sans contrôle), ceci permet de créer des segments, triangles ou quadrangles dans l'espace R²
• NBS : le nombre de sommets du maillage grossier
• NA : le nombre d'arêtes du maillage grossier
• NBGRO : le nombre de blocs (ou éléments) du maillage grossier
• NBFR : le nombre de faces à décrire en vue de leurs affecter un numéro de référence non nul
• XYZ(3,.) : les coordonnées des sommets grossiers (XYZ(J,I) est la coordonnée J du sommet I)
• NARET(5,.) : le tableau des arêtes grossières avec:
  • NARET(1,I) : le numéro du sommet origine de l'arête I
  • NARET(2,I) : le numéro de son extrémité
  • NARET(3,I) : le code de découpage de l'arête: 0 pour une arête droite découpée automatiquement en segments de même longueur; 1 si on donne les points intermédiaires (cf. XYZINT)
  • NARET(4,I) : le nombre de points de contrôle à créer sur l'arête (extrémités exclues)
  • NARET(5,I) : la référence de l'arête
• XYZINT(3,.) : les coordonnées des points de contrôle des arêtes donnés effectivement par l'utilisateur. Pour chaque arête I dont NARET(3,I) = 1, on donne ces valeurs de la façon suivante: XYZINT(J,I) est la coordonnée J du point I, où I est le numéro global du point; on compte d'abord les points (extrémités exclues) de la première arête ainsi décrite, puis ceux de la deuxième, etc.
• NTFR(6,.) : le tableau des faces à spécifier avec:
  • NTFR(1,I) : la référence de la face I
  • NTFR(2,I) : son code géométrique (3: triangle, 4: quadrangle)
  • NTFR(3,I) à NTFR(6,I) : la liste de ses sommets (0 pour le sommet '4' d'un triangle)
• NTGRO(11,.) : le tableau des éléments grossiers avec:
  • NTGRO(1,I) : le type géométrique du bloc I (2: segment, 3: triangle,
    4: quadrangle, 5: tétraèdre, 6: pentaèdre, 7: hexaèdre)
  • NTGRO(2,I) à NTGRO(9,I) : la liste de ses sommets (0 pour les sommets n'existant pas)
  • NTGRO(10,I) : le type géométrique des éléments finis issus du découpage du bloc (cf. remarque)
  • NTGRO(11,I) : le numéro de sous-domaine du bloc I (donc celui de tous les sous-éléments créés)

Remarques :

• le module vérifiera la positivité des volumes des éléments créés, dans le cas où un ou plusieurs de ceux-ci sont négatifs, cela signifie que soit les points de contrôle des arêtes sont mal donnés (à l'envers en particulier), soit que la géométrie du bloc considéré est trop 'tordue'
• dans la présente version, on a NTGRO(10,I) = NTGRO(1,I)