Télécharger mhybr3.eso

Retour à la liste

Numérotation des lignes :

  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