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

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