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