Télécharger la plaquette Cast3M


fr en

Télécharger hyshpt.eso

Retour à la liste

Numérotation des lignes : 

hyshpt

  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.  
  46. -INC PPARAM
  47. -INC CCOPTIO
  48. C
  49.       SEGMENT HYBSTO
  50.          REAL*8 HYBASE(NDIM,NBDDL,NBPP)
  51.       ENDSEGMENT
  52. C
  53.       SEGMENT SHXX
  54.          REAL*8 SHY(NDIM,NBDDL)
  55.       ENDSEGMENT
  56. C
  57. C= Initialisations
  58. C
  59.       MINTE = IPINTE
  60.       SEGACT MINTE
  61.       NBPP  = SHPTOT(/3)
  62.       NDIM  = IDIM * (IDIM+1)
  63.       SEGINI HYBSTO , SHXX
  64. *
  65. *= Calcul des fonctions de base en chaque point du segment d'integration
  66. *
  67.       DO 20 II=1,NBPP
  68.          XX = QSIGAU(II)
  69.          YY = ETAGAU(II)
  70.          ZZ = DZEGAU(II)
  71.          CALL HYSHP(XX,YY,ZZ,NEF,NBDDL,NDIM,SHY,IRET)
  72.          IF (IRET.NE.0) GOTO 100
  73.          DO 10 ID=1,NDIM
  74.             DO 10 NO=1,NBDDL
  75.                HYBASE(ID,NO,II) = SHY(ID,NO)
  76.   10     CONTINUE
  77.   20  CONTINUE
  78. C
  79. C= Desactivation et menage
  80. C
  81.       IPTHYB = HYBSTO
  82.       SEGDES MINTE , HYBSTO
  83.       SEGSUP SHXX
  84.       RETURN
  85. *
  86. *= ERREUR : Fonctions de forme non encore implementees
  87. *
  88.  100  CONTINUE
  89.       MOTERR(1:4) = NOMTP(NEF)
  90.       CALL ERREUR(68)
  91.       SEGSUP SHXX , HYBSTO
  92.       SEGDES MINTE
  93.       IPTHYB = 0
  94.       RETURN
  95.       END
  96.  
  97.  
  98.  
  99.  
  100.  

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