Télécharger la plaquette Cast3M


fr en

Télécharger typfac.eso

Retour à la liste

Numérotation des lignes : 

typfac

  1. C TYPFAC    SOURCE    OF166741  25/02/06    21:15:07     12146          
  2.  
  3. C=======================================================================
  4. C
  5. C     DONNE LA FORMULATION D UNE FACE DE ELEMENT MASSIF MELM
  6. C  ENTREE
  7. C     MELM = NUMERO DE L ELEMENT MASSIF DANS NOMTP
  8. C     NBNN = NOMBRE DE NOEUDS DE LA FACE
  9. C  SORTIE
  10. C     IRET = NUMERO DE LA FACE DANS NOMTP
  11. C            0 SI OPERATION IMPOSSIBLE
  12. C
  13. C======================================================================
  14.       SUBROUTINE TYPFAC(MELM,NBNN,IRET)
  15.  
  16.       IMPLICIT INTEGER(I-N)
  17.       IMPLICIT REAL*8(A-H,O-Z)
  18.  
  19. C==DEB= FORMULATION HHO == INCLUDE =====================================
  20. -INC CCHHOPA
  21. C==FIN= FORMULATION HHO ================================================
  22.  
  23. C
  24. C     ELEMENTS MASSIFS BIDIMENSIONNELS LINEAIRES
  25. C
  26.       IF (MELM.EQ. 4 .OR. MELM.EQ. 5 .OR. MELM.EQ.8 .OR. MELM.EQ.9 .OR.
  27.      .    MELM.EQ.69 .OR. MELM.EQ.70 .OR. MELM.EQ.183 .OR.
  28.      .    MELM.EQ.263) THEN
  29.         IRET = 2
  30. C
  31. C     ELEMENTS MASSIFS BIDIMENSIONNELS POLYGONAUX
  32. C
  33.       ELSE IF (MELM.GE.111.AND.MELM.LE.122) THEN
  34.         IRET = 2
  35. C
  36. C     ELEMENTS MASSIFS BIDIMENSIONNELS QUADRATIQUES
  37. C
  38.       ELSE IF (MELM.EQ.6.OR.MELM.EQ.7.OR.MELM.EQ.10.OR.MELM.EQ.11
  39.      .      .OR.MELM.EQ.79.OR.MELM.EQ.80.OR.MELM.EQ.71.OR.MELM.EQ.72
  40.      .      .OR.MELM.EQ.184) THEN
  41.         IRET = 3
  42. C
  43. C     ELEMENTS MASSIFS TRIDIMENSIONNELS LINEAIRES
  44. C
  45.       ELSE IF (MELM.EQ.14.OR.MELM.EQ.16.OR.MELM.EQ.23.OR.MELM.EQ.25.OR.
  46.      .   MELM.EQ.73.OR.MELM.EQ.74.OR.MELM.EQ.75.OR.MELM.EQ.264.OR.
  47.      .   MELM.EQ.273)THEN
  48.         IF (NBNN.EQ.3) THEN
  49.           IRET = 31
  50.         ELSE
  51.           IRET = 32
  52.         ENDIF
  53. C
  54. C     ELEMENTS MASSIFS TRIDIMENSIONNELS QUADRATIQUES
  55. C
  56.       ELSE IF (MELM.EQ.15.OR.MELM.EQ.17.OR.MELM.EQ.24 .OR. MELM.EQ.26
  57.      .     .OR.MELM.EQ.76.OR.MELM.EQ.77.OR.MELM.EQ.78 .OR. MELM.EQ.274
  58.      .     .OR.MELM.EQ.81.OR.MELM.EQ.82.OR.MELM.EQ.83
  59.      &     .OR.MELM.EQ.279 .OR. MELM.EQ.280) THEN
  60.         IF (NBNN.EQ.6) THEN
  61.           IRET = 33
  62.         ELSE
  63.           IRET = 34
  64.         ENDIF
  65. C
  66. C elements fluides 2d et 3d
  67. C
  68.       else if (melm.eq.35 .or. melm.eq.36) then
  69.         iret = 97
  70.       else if (melm.eq.38. or. melm.eq.40) then
  71.         if (nbnn.eq.3) then
  72.           iret = 35
  73.         else
  74.           iret = 36
  75.         endif
  76.       else if (melm.eq.37) then
  77.         iret = 36
  78.       else if (melm.eq.39) then
  79.         iret = 35
  80.  
  81. C= Elements MASSIFs UNIDIMENSIONNELs (1D)
  82.       ELSE IF (MELM.EQ.193.OR.MELM.EQ.194) THEN
  83.         IRET = 45
  84.         NBNN = 1
  85.  
  86. C==DEB= FORMULATION HHO ================================================
  87.       ELSE IF (MELM .EQ. HHO_NUM_ELEMENT) THEN
  88.         IF (NBNN.EQ.2) THEN
  89.           IRET = 2
  90.         ELSE
  91.           write(ioimp,*) 'HHO TYPFAC: not implemented'
  92.           call erreur(5)
  93.           return
  94.         ENDIF
  95. C==FIN= FORMULATION HHO ================================================
  96.  
  97.       ELSE
  98.         IRET=0
  99.       ENDIF
  100.  
  101. c      RETURN
  102.       END
  103.  
  104.  
  105.  
  106.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales