enfr

Calcul de volumogramme

31 mai 2011

par N. Tardieu, EDF R&D / AMA

Un volumogramme donne une représentation de la composante d’un champ en terme de pourcentage du volume de la pièce contenant une plage de valeurs donnée.

Ainsi, si l’on étudie la pièce présentée Figure 1 et que l’on s’intéresse à la répartition de la contrainte de Von Mises, le volumogramme est le suivant :

Figure 1 : Maillage du pallier étudié
   INTERVALLE         BORNE_INF      BORNE_SUP       DISTRIBUTION
------------------------------------------------------------------------------------
       1              3.47944E-02    4.00516E+01     2.21569E+01
       2              4.00516E+01    8.00685E+01     2.26231E+01
       3              8.00685E+01    1.20085E+02     1.98340E+01
       4              1.20085E+02    1.60102E+02     1.24788E+01
       5              1.60102E+02    2.00119E+02     8.81198E+00
       6              2.00119E+02    2.40136E+02     5.74077E+00
       7              2.40136E+02    2.80153E+02     3.48664E+00
       8              2.80153E+02    3.20170E+02     2.12184E+00
       9              3.20170E+02    3.60186E+02     1.47352E+00
      10              3.60186E+02    4.00203E+02     1.27246E+00

La plage de variation de la contrainte de Von Mises s’étend de 3.47944E-02 à 4.00203E+02 MPa et la découpe en 10 intervalles :

  • la partie de la structure qui voit cette contrainte varier de 3.47944E-02 à 4.00516E+01 MPa représente 22% du volume de structure ;
  • la partie de la structure qui voit cette contrainte varier de 4.00516E+01 à 8.00685E+01 MPa représente 22% du volume de structure ;
  • etc…
Figure 2 : Répartition de la contrainte de Von Mises

Le volumogramme peut être aisément visualisé sous forme d’histogramme (Figure 3) où il fournit une vision très synthétique de la répartition d’un champ.

Figure 3 : Représentation graphique du volumogramme

La mise en œuvre dans Code_Aster est très simple grâce à la commande POST_ELEM, où l’utilisateur précise le nom du champ, la composante à traiter.

Table=POST_ELEM(RESULTAT=RESU,
     VOLUMOGRAMME=_F(NOM_CHAM='SIEQ_ELGA',
     NOM_CMP='VMIS',),);