Suiv.: Remontée aux flux ou aux Sup.: 2 Manipulation de solutions Préc.: Calcul des forces d'interface et Index Table des matières

Remontée aux contraintes (S.D. TAE)

Le préprocesseur STREXX permet de construire la S.D. TAE des contraintes obtenue en remontant les tableaux des contraintes élémentaires de la S.D. TAE, établie au moment où l'on calcule matrices et seconds membres élémentaires, par la S.D. B contenant les déplacements, S.D. obtenue lors du calcul de la solution du problème d'élasticité considéré.

Tous les éléments finis en élasticité sont possibles.

Le module STRESS, appelé par STREXX calcule les contraintes dans chaque élément (au barycentre ou aux points d'intégration numériques (cf. fiches techniques)) et, de plus, peut imprimer les contraintes de tout ou partie des éléments du maillage. Dans le cas d'un problème bidimensionnel (la visualisation du cas 3D n'étant pas implémentée), la structure TAE issue de STRESS est utilisée pour visualiser les contraintes (cf. TRSTXX dans la partie 2 de ce guide).

Le programme STRESS se présente sous la forme:

      SUBROUTINE STRESS(M,XM,ISTRES,LINDIC,INDIC,NFBE,NIBE,NFBT,NIBT,
     +                  NFTAE,NITAE,NFNDL1,NINDL1,NFNDL2,NINDL2,
     +                  NFMAIL,NIMAIL,NFTAES,NITAES,NTTAES)
C   ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C   MODULE  STRESS : CALCUL DES CONTRAINTES LORS DE LA RESOLUTION D UN
C   --------------   PROBLEME  D'ELASTICITE A PARTIR DES CONTRAINTES
C                    ELEMENTAIRES ET DE LA SOLUTION
C   ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

où:

• M et XM représentent le tableau de travail,
• ISTRES est une option permettant l'impression des contraintes. ISTRES vaut IJK avec:

  .

  I = 1, impression de tous les éléments,

  .

  I = 2, impression des éléments compris entre les numéros N1 et N2 donnés dans INDIC,

  .

  I = 3, impression des LINDIC éléments dont les numéros sont donnés dans INDIC,

  .

  I = 4, impression des éléments dont les numéros (LINDIC numéros) de sous-domaine sont donnés dans INDIC,

  .

  J = 1, si l'impression des contraintes principales est demandée, 0 sinon,

  .

  K = 1, si l'impression des directions principales est demandée (I forcé à 1), 0 sinon,

• LINDIC, fonction de I, est la longueur du tableau INDIC,
• INDIC, fonction de I, est le tableau suivant:

  .

  I = 1, le tableau est inutilisé,

  .

  I = 2, de longueur 2, le tableau contient le numéro du premier et du dernier élément à imprimer,

  .

  I = 3 , le numéro des LINDIC éléments à imprimer,

  .

  I = 4, le numéro des LINDIC sous-domaines à traiter. Seront imprimées les contraintes de tous les éléments ayant ces numéros de sous-domaine.

• NFBE, NIBE sont le numéro de fichier et le niveau de la S.D. B contenant les déplacements (NFBE > 0: le fichier est séquentiel, NFBE = 0, la S.D. est en mémoire centrale, NFBE < 0, le fichier est en accès direct),
• NFBT, NIBT, idem pour la S.D. B contenant les températures (cas thermoélastique),
• NFTAE, NITAE sont le numéro de fichier et le niveau de la S.D. TAE contenant les tableaux élémentaires (dont celui des contraintes élémentaires). NFTAE > 0.
• NFNDL1, NINDL1, utiles si ND = 0, sont le numéro de fichier et le niveau de la S.D. NDL1 élastique (ND = 0 indique que le nombre de degrés de liberté par noeud n'est pas constant),
• NFNDL2, NINDL2, utiles si ND = 0, sont le numéro de fichier et le niveau de la S.D. NDL1 thermique,
• NFMAIL, NIMAIL sont le numéro de fichier et le niveau de la S.D. MAIL,
• NFTAES, NITAES, NTTAES sont le numéro de fichier, le niveau et le nombre de tableaux à associer de la S.D.S. TAE qui contiendra les contraintes élément par élément (NFTAES > 0).

Les fichiers d'entrée nécessaires sont NFBE, NFTAE, NFMAIL et éventuellement NFBT, NFNDL1 et NFNDL2. Le fichier de sortie est NFTAES.

Il n'y a pas de données à fournir spécifiquement pour ce module sauf dans le cas où NTTAES est non nul. Alors, il faudra fournir:

Boucle I de 1 à NTTAES

• nom du tableau I sur 4 caractères,
• type de ce tableau (1 = entier, 2 = réel simple précision, 4 = character, 5 = réel double précision, ...)
• nombre de mots de ce tableau
• les variables de ce tableau
• commentaire sur le contenu du tableau I sur 18 mots de 4 caractères.

La structure de sortie TAE est décrite dans le [Guide Modulef - 2].

Notons que selon ISTRES, le calcul des valeurs et vecteurs propres est effectué ou non. Le sous-programme VPM2DD calcule les valeurs et vecteurs propres dans le cas bidimensionnel. VAM3DD calcule les valeurs propres tandis que VPM3DD calcule les vecteurs propres en 3-dimensions.

La structure TAE issue du module STRESS contient les informations relatives aux contraintes de tous les éléments du maillage. Le module COMTAE permet de sélectionner certains éléments seulement et de construire la S.D. TAE correspondante.

      SUBROUTINE COMTAE(M,ISTRES,INDIC,LINDIC,NFTAE,NITAE,NFTAES,NITAES)
C  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C  MODULE COMTAE : COMPRESSER LA S.D. TAES ( ISSUE DE STRESS )
C  -------------   ON NE GARDE QUE LES ELEMENTS DE INDIC ET ISTRES
C  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

où:

• M est le tableau de travail,
• ISTRES est une option permettant le choix des éléments. ISTRES vaut I00 avec:

  .

  I = 1, le choix de tous les éléments,

  .

  I = 2, le choix des éléments compris entre les numéros N1 et N2 donnés dans INDIC,

  .

  I = 3, le choix des LINDIC éléments dont les numéros sont donnés dans INDIC,

  .

  I = 4, le choix des éléments dont les numéros (LINDIC numéros) de sous-domaine sont donnés dans INDIC,

• INDIC, fonction de I, est le tableau suivant:

  .

  I = 1, le tableau est inutilisé,

  .

  I = 2, de longueur 2, le tableau contient le numéro du premier et du dernier élément à choisir,

  .

  I = 3 , le numéro des LINDIC éléments à choisir,

  .

  I = 4, le numéro des LINDIC sous-domaines à traiter. Seront sélectionnés tous les éléments ayant ces numéros de sous-domaine.

• LINDIC, fonction de I, est la longueur du tableau INDIC,
• NFTAE, NITAE sont le numéro de fichier et le niveau de la S.D. TAE issue de STRESS,
• NFTAES, NITAES sont le numéro de fichier et le niveau de la S.D.S. TAE qui contiendra les contraintes des élément sélectionnés (NFTAES > 0 et NFTAES NFTAE).

Le fichier d'entrée nécessaire est NFTAE. Le fichier de sortie est NFTAES.

Il n'y a pas de données à fournir spécifiquement pour ce module.