Télécharger hyshpt.eso

Retour à la liste

Numérotation des lignes :

  1. C HYSHPT SOURCE BP208322 15/06/22 21:19:21 8543
  2. SUBROUTINE HYSHPT(NEF,NBDDL,IPINTE,IPTHYB)
  3. C-----------------------------------------------------------------------
  4. C Calcul les fonctions de bases de la formulation hybride pour l'ef NEF
  5. C aux points contenus dans le segment d'integration IPINTE.
  6. C-----------------------------------------------------------------------
  7. C
  8. C---------------------------
  9. C Parametres Entree/Sortie :
  10. C---------------------------
  11. C
  12. C E/ NEF : Numero de l'element fini dans NOMTP (cf common CCHAMP)
  13. C E/ NBDDL : Nombre de noeuds portant des degres de liberte
  14. C E/ IPINTE : Pointeur vers le segment d'integration des fonctions
  15. C de forme de l'element fini NEF
  16. C /S IPTHYB : Pointeur du segment HYBSTO contenant les fonctions de
  17. C bases en formulation hybride
  18. C
  19. C----------------------
  20. C Variables en COMMON :
  21. C----------------------
  22. C
  23. C E/ IDIM : Dimension de l'espace (cf CCOPTIO)
  24. C E/ NOMTP : Noms des elements finis (cf CCHAMP)
  25. C /S MOTERR : Chaine de char pour erreur (cf CCOPTIO)
  26. C
  27. C----------------------
  28. C Tableaux de travail :
  29. C----------------------
  30. C
  31. C SHP : Contient les fonctions de base en un point
  32. C
  33. C-----------------------------------------------------------------------
  34. C
  35. C Langage : ESOPE + FORTRAN77
  36. C
  37. C Auteurs : F.DABBENE 08/93
  38. C
  39. C-----------------------------------------------------------------------
  40. IMPLICIT INTEGER(I-N)
  41. IMPLICIT REAL*8(A-H,O-Z)
  42. C
  43. -INC CCHAMP
  44. -INC SMINTE
  45. -INC CCOPTIO
  46. C
  47. SEGMENT HYBSTO
  48. REAL*8 HYBASE(NDIM,NBDDL,NBPP)
  49. ENDSEGMENT
  50. C
  51. SEGMENT SHXX
  52. REAL*8 SHY(NDIM,NBDDL)
  53. ENDSEGMENT
  54. C
  55. C= Initialisations
  56. C
  57. MINTE = IPINTE
  58. SEGACT MINTE
  59. NBPP = SHPTOT(/3)
  60. NDIM = IDIM * (IDIM+1)
  61. SEGINI HYBSTO , SHXX
  62. *
  63. *= Calcul des fonctions de base en chaque point du segment d'integration
  64. *
  65. DO 20 II=1,NBPP
  66. XX = QSIGAU(II)
  67. YY = ETAGAU(II)
  68. ZZ = DZEGAU(II)
  69. CALL HYSHP(XX,YY,ZZ,NEF,NBDDL,NDIM,SHY,IRET)
  70. IF (IRET.NE.0) GOTO 100
  71. DO 10 ID=1,NDIM
  72. DO 10 NO=1,NBDDL
  73. HYBASE(ID,NO,II) = SHY(ID,NO)
  74. 10 CONTINUE
  75. 20 CONTINUE
  76. C
  77. C= Desactivation et menage
  78. C
  79. IPTHYB = HYBSTO
  80. SEGDES MINTE , HYBSTO
  81. SEGSUP SHXX
  82. RETURN
  83. *
  84. *= ERREUR : Fonctions de forme non encore implementees
  85. *
  86. 100 CONTINUE
  87. MOTERR(1:4) = NOMTP(NEF)
  88. CALL ERREUR(68)
  89. SEGSUP SHXX , HYBSTO
  90. SEGDES MINTE
  91. IPTHYB = 0
  92. RETURN
  93. END
  94.  
  95.  
  96.  
  97.  
  98.  

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