Télécharger la plaquette Cast3M


fr en

Télécharger mhybr3.eso

Retour à la liste

Numérotation des lignes : 

mhybr3

  1. C MHYBR3    SOURCE    CHAT      05/01/13    01:44:47     5004
  2.       SUBROUTINE MHYBR3(IGAU,NBNN,NBDDL,NDIM,IDIM,NDIM2,XEL,HYBASE,
  3.      S                  SHPTOT,SHY,SHP,ZJAC,DJAC)
  4. C-----------------------------------------------------------------------
  5. C Recuperation de la matrice Jacobienne et de son determinant ainsi
  6. C que des fonctions de base hybride au point de gauss IGAU
  7. C-----------------------------------------------------------------------
  8. C
  9. C---------------------------
  10. C Parametres Entree/Sortie :
  11. C---------------------------
  12. C
  13. C E/  IGAU    : Numero du point de gauss considere.
  14. C E/  NBNN    : Nombre de noeud du support geometrique
  15. C E/  NBDDL   : Nombre de noeud portant les degre de liberte
  16. C E/  NDIM    : Premiere dimension de HYBASE (=IDIM*(IDIM+1))
  17. C E/  IDIM    : Dimension de l'espace (2 en 2D et 3 en 3D).
  18. C E/  NDIM2   : Dimension de SHY. Vaut IDIM si on ne recupere que la
  19. C               fonction de base et NDIM si on recupere aussi les
  20. C               derivees
  21. C E/  XEL     : Coordonnees reelles des noeuds de l'element.
  22. C E/  HYBASE  : Fonctions de base hybride pour tout les points de gauss
  23. C E/  SHPTOT  : Fonctions de forme pour tout les points de gauss
  24. C  /S SHY     : Fonctions de base hybride au point de gauss IGAU
  25. C  /S SHP     : Fonctions de forme au point de gauss IGAU
  26. C  /S ZJAC    : Matrice jacobienne au point de gauss IGAU
  27. C  /S DJAC    : Determinant de la matrice ZJAC
  28. C
  29. C----------------------
  30. C Tableaux de travail :
  31. C----------------------
  32. C
  33. C     SHP     : Bien que remplit par necessite, cette variable ne
  34. C               sert a rien en amont.
  35. C
  36. C-----------------------------------------------------------------------
  37. C
  38. C Langage : FORTRAN77
  39. C
  40. C Auteurs : F.DABBENE 08/93
  41. C
  42. C-----------------------------------------------------------------------
  43.       IMPLICIT INTEGER(I-N)
  44.       IMPLICIT REAL*8 (A-H,O-Z)
  45. C
  46.       REAL*8 XEL(3,NBNN),SHPTOT(6,NBNN,*),SHP(6,NBNN),
  47.      S       ZJAC(IDIM,IDIM),HYBASE(NDIM,NBDDL,*),SHY(NDIM2,NBDDL)
  48. C
  49. C= Calcul de la matrice jacobienne et de son determinant
  50. C
  51.       IFIN = IDIM + 1
  52.       DO 20 NP=1,NBNN
  53.          DO 10 I=1,IFIN
  54.             SHP(I,NP) = SHPTOT(I,NP,IGAU)
  55.  10      CONTINUE
  56.  20   CONTINUE
  57.       CALL MATJAC(XEL,SHP,IDIM,NBNN,ZJAC,DJAC)
  58. C
  59. C= Recuperation des fonctions de base hybride
  60. C
  61.       DO 40 NP=1,NBDDL
  62.          DO 30 I=1,NDIM2
  63.             SHY(I,NP) = HYBASE(I,NP,IGAU)
  64.  30      CONTINUE
  65.  40   CONTINUE
  66. *
  67.       RETURN
  68.       END
  69.  
  70.  

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