C TTET4F    SOURCE    OF166741  23/12/04    21:15:23     11800          
 
C=======================================================================
C=                            T T E T 4 F                              =
C=                            -----------                              =
C=                                                                     =
C=  Fonction :                                                         =
C=  ----------                                                         =
C=   Acquisition des caracteristiques d'integration aux noeuds, aux    =
C=   points d'integration (Gauss) ou au centre de gravite pour         =
C=   l'element thermique support TET4.                                 =
C=                                                                     =
C=  Parametres :  (E)=Entree  (S)=Sortie                               =
C=  ------------                                                       =
C=   ISUPTH  (E)   Indique le support d'integration demande            =
C=   IPINTE  (S)   Pointeur sur le segment MINTE associe (ACTIF en S)  =
C=======================================================================

      SUBROUTINE TTET4F (ISUPTH,IPINTE)

      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)

-INC SMINTE

      PARAMETER ( UN=1.D0 )
      PARAMETER ( XZER=0.D0 )
      PARAMETER ( UNDEMI=.5D0 )
      PARAMETER ( UNQUA=.25D0 )
      PARAMETER ( UNSIX=1.D0/6.D0 )
      PARAMETER ( OHUIT=-.80D0 )
      PARAMETER ( OQUCI= .45D0 )
      PARAMETER ( UNDOU=1.D0/12.D0 )

      NBNO = 4
      IF      (ISUPTH.EQ.1) THEN
        NBPGAU = 4
      ELSE IF (ISUPTH.EQ.2) THEN
        NBPGAU = 5
      ELSE IF (ISUPTH.EQ.3) THEN
        NBPGAU = 1
      ELSE
        RETURN
      ENDIF
      SEGINI,MINTE

C COORDONNEES ET POIDS DES 4 NOEUDS
      IF (ISUPTH.EQ.1) THEN
        QSIGAU(1) = XZER
        QSIGAU(2) = UN
        QSIGAU(3) = XZER
        QSIGAU(4) = XZER
        ETAGAU(1) = XZER
        ETAGAU(2) = XZER
        ETAGAU(3) = UN
        ETAGAU(4) = XZER
        DZEGAU(1) = XZER
        DZEGAU(2) = XZER
        DZEGAU(3) = XZER
        DZEGAU(4) = UN
        POIGAU(1) = UNDOU
        POIGAU(2) = UNDOU
        POIGAU(3) = UNDOU
        POIGAU(4) = UNDOU

C COORDONNEES ET POIDS DES 5 POINTS DE GAUSS
      ELSE IF (ISUPTH.EQ.2) THEN

        QSIGAU(1) = UNQUA
        QSIGAU(2) = UNSIX
        QSIGAU(3) = UNSIX
        QSIGAU(4) = UNSIX
        QSIGAU(5) = UNDEMI

        ETAGAU(1) = UNQUA
        ETAGAU(2) = UNSIX
        ETAGAU(3) = UNSIX
        ETAGAU(4) = UNDEMI
        ETAGAU(5) = UNSIX

        DZEGAU(1) = UNQUA
        DZEGAU(2) = UNSIX
        DZEGAU(3) = UNDEMI
        DZEGAU(4) = UNSIX
        DZEGAU(5) = UNSIX

        POIGAU(1) = OHUIT*UNSIX
        POIGAU(2) = OQUCI*UNSIX
        POIGAU(3) = OQUCI*UNSIX
        POIGAU(4) = OQUCI*UNSIX
        POIGAU(5) = OQUCI*UNSIX

C COORDONNEES ET POIDS DU CENTRE DE GRAVITE
      ELSE IF (ISUPTH.EQ.3) THEN
        QSIGAU(1) = XZer
        ETAGAU(1) = XZer
        DZEGAU(1) = XZer
        POIGAU(1) = UN/(UN+UN+UN)

      ENDIF

C FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S)
      DO I = 1, NBPGAU
        XI = QSIGAU(I)
        YI = ETAGAU(I)
        ZI = DZEGAU(I)

        AUX = UN-XI-YI-ZI
        SHPTOT(1,1,I) = AUX
        SHPTOT(1,2,I) = XI
        SHPTOT(1,3,I) = YI
        SHPTOT(1,4,I) = ZI

        SHPTOT(2,1,I) = -UN
        SHPTOT(2,2,I) =  UN
        SHPTOT(2,3,I) =  XZER
        SHPTOT(2,4,I) =  XZER

        SHPTOT(3,1,I) = -UN
        SHPTOT(3,2,I) =  XZER
        SHPTOT(3,3,I) =  UN
        SHPTOT(3,4,I) =  XZER

        SHPTOT(4,1,I) = -UN
        SHPTOT(4,2,I) =  XZER
        SHPTOT(4,3,I) =  XZER
        SHPTOT(4,4,I) =  UN

      ENDDO

      IPINTE = MINTE

c      RETURN
      END