C BBGRDK    SOURCE    CHAT      05/01/12    21:36:32     5004
      SUBROUTINE BBGRDK(QSI,ETA,GEOM,BB)
C |====================================================================|
C |                                                                    |
C |   EVALUATION DE LA MATRICE BB(2,9) AU POINT QSI ETA                |
C |   POUR LES TERMS DE GRADIAN DE L'ELEMENT DKT                       |
C |     DUXDZ  DUYDZ   DUZDX  DUZDY        MAIS                        |
C |     DUXDZ=-DUZDX=+BETA(X)=BB(1,I)   (I=1,2...9)                    |
C |     DUYDZ=-DUZDY=+BETA(Y)=BB(2,I)   (I=1,2...9)                    |
C |                                                                    |
C |   CODE SUO X.Z.  23/01/1987                                        |
C |====================================================================|
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION GEOM(*),BB(2,*)
C
      DATA TRDE/1.5D0/,DEUX/2.D0/,QUATRE/4.D0/,UN/1.D0/
C
      A4=GEOM(1)
      A5=GEOM(2)
      A6=GEOM(3)
      B4=GEOM(4)
      B5=GEOM(5)
      B6=GEOM(6)
      C4= GEOM(7)
      C5= GEOM(8)
      C6= GEOM(9)
      D4=GEOM(10)
      D5=GEOM(11)
      D6=GEOM(12)
      E4=GEOM(13)
      E5=GEOM(14)
      E6=GEOM(15)
      AUX=UN-QSI-ETA
C
C       FONCTIONS DE FORME DE TRIANGLE A 6 NOEUDS
C
        SN1=-AUX*(1-DEUX*AUX)
        SN2=-QSI*(1-DEUX*QSI)
        SN3=-ETA*(1-DEUX*ETA)
        SN4=QUATRE*QSI*ETA
        SN5=QUATRE*ETA*AUX
        SN6=QUATRE*QSI*AUX
C
C       FONCTIONS DE FORME DE TRIANGLE DKT
C
        BB(1,1)= TRDE*(A6*SN6-A5*SN5)
        BB(1,2)= B5*SN5+B6*SN6
        BB(1,3)= SN1-C5*SN5-C6*SN6
        BB(1,4)= TRDE*(A4*SN4-A6*SN6)
        BB(1,5)= B6*SN6+B4*SN4
        BB(1,6)= SN2-C6*SN6-C4*SN4
        BB(1,7)= TRDE*(A5*SN5-A4*SN4)
        BB(1,8)= B4*SN4+B5*SN5
        BB(1,9)= SN3-C4*SN4-C5*SN5
        BB(2,1)= TRDE*(D6*SN6-D5*SN5)
        BB(2,2)= -SN1+E5*SN5+E6*SN6
        BB(2,3)= -B5*SN5-B6*SN6
        BB(2,4)= TRDE*(D4*SN4-D6*SN6)
        BB(2,5)= -SN2+E6*SN6+E4*SN4
        BB(2,6)= -B6*SN6-B4*SN4
        BB(2,7)= TRDE*(D5*SN5-D4*SN4)
        BB(2,8)= -SN3+E4*SN4+E5*SN5
        BB(2,9)= -B4*SN4-B5*SN5
C
C
      RETURN
      END