Télécharger typfac.eso

Retour à la liste

Numérotation des lignes :

typfac
  1. C TYPFAC SOURCE KICH 18/01/11 21:16:00 9690
  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. IF (MELM.EQ.4.OR.MELM.EQ.5.OR.MELM.EQ.8.OR.MELM.EQ.9.OR.
  16. . MELM.EQ.69.OR.MELM.EQ.70.OR.MELM.EQ.183.or.
  17. . MELM.eq.263) THEN
  18. C
  19. C ELEMENTS MASSIFS BIDIMENSIONNELS LINEAIRES
  20. C
  21. IRET=2
  22. ELSE IF (MELM.GE.111.AND.MELM.LE.122) THEN
  23. C
  24. C ELEMENTS MASSIFS BIDIMENSIONNELS POLYGONAUX
  25. C
  26. IRET=2
  27. ELSE IF (MELM.EQ.6.OR.MELM.EQ.7.OR.MELM.EQ.10.OR.MELM.EQ.11
  28. . .OR.MELM.EQ.79.OR.MELM.EQ.80.OR.MELM.EQ.71.OR.MELM.EQ.72
  29. . .OR.MELM.EQ.184) THEN
  30. C
  31. C ELEMENTS MASSIFS BIDIMENSIONNELS QUADRATIQUES
  32. C
  33. IRET=3
  34. ELSE IF (MELM.EQ.14.OR.MELM.EQ.16.OR.MELM.EQ.23.OR.MELM.EQ.25.OR.
  35. . MELM.EQ.73.OR.MELM.EQ.74.OR.MELM.EQ.75.OR.MELM.EQ.264.OR.
  36. . MELM.EQ.273)THEN
  37. C
  38. C ELEMENTS MASSIFS TRIDIMENSIONNELS LINEAIRES
  39. C
  40. IF (NBNN.EQ.3) THEN
  41. IRET=31
  42. ELSE
  43. IRET=32
  44. ENDIF
  45. ELSE IF (MELM.EQ.15.OR.MELM.EQ.17.OR.MELM.EQ.24.OR.MELM.EQ.26
  46. . .OR.MELM.EQ.76.OR.MELM.EQ.77.OR.MELM.EQ.78.OR.MELM.EQ.274
  47. . .OR.MELM.EQ.81.OR.MELM.EQ.82.OR.MELM.EQ.83) THEN
  48. C
  49. C ELEMENTS MASSIFS TRIDIMENSIONNELS QUADRATIQUES
  50. C
  51. IF (NBNN.EQ.6) THEN
  52. IRET=33
  53. ELSE
  54. IRET=34
  55. ENDIF
  56. C
  57. C elements fluides 2d et 3d
  58. C
  59. else if (melm.eq.35 .or. melm.eq.36) then
  60. iret = 97
  61. else if (melm.eq.38. or. melm.eq.40) then
  62. if (nbnn.eq.3) then
  63. iret = 35
  64. else
  65. iret = 36
  66. endif
  67. else if (melm.eq.37) then
  68. iret = 36
  69. else if (melm.eq.39) then
  70. iret = 35
  71.  
  72. C= Elements MASSIFs UNIDIMENSIONNELs (1D)
  73. ELSE IF (MELM.EQ.193.OR.MELM.EQ.194) THEN
  74. IRET=45
  75. NBNN=1
  76.  
  77. ELSE
  78. IRET=0
  79. ENDIF
  80. RETURN
  81. END
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  

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