Esope source files

   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
  56 :       CALL MATJAC(XEL,SHP,IDIM,NBNN,ZJAC,DJAC)
  57 : C
  58 : C= Recuperation des fonctions de base hybride
  59 : C
  60 :       DO 40 NP=1,NBDDL
  61 :          DO 30 I=1,NDIM2
  62 :             SHY(I,NP) = HYBASE(I,NP,IGAU)
  63 :  30      CONTINUE
  64 :  40   CONTINUE
  65 : *                                                                               
  66 :       RETURN
  67 :       END                                                                       
  68 :