C TTRI3F SOURCE OF166741 23/12/04 21:15:24 11800 ======================================================================= C= T T R I 3 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 TRI3. = 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 TTRI3F (ISUPTH,IPINTE) IMPLICIT INTEGER(I-N) IMPLICIT REAL*8 (A-H,O-Z) -INC SMINTE PARAMETER ( O0=0.D0 ) PARAMETER ( O1=1.D0 ) PARAMETER ( O1s2=0.5D0 ) PARAMETER ( O1s3=0.333333333333333333333333333333333333D0 ) PARAMETER ( O66=0.666666666666666666666666666666666667D0 ) PARAMETER ( O16=0.166666666666666666666666666666666667D0 ) NBNO = 3 IF (ISUPTH.EQ.1 .OR. ISUPTH.EQ.2) THEN NBPGAU = 3 ELSE IF (ISUPTH.EQ.3) THEN NBPGAU = 1 ELSE CALL ERREUR(5) RETURN ENDIF SEGINI,MINTE C COORDONNEES DES 3 NOEUDS IF (ISUPTH.EQ.1) THEN QSIGAU(1) = O0 QSIGAU(2) = O1 QSIGAU(3) = O0 ETAGAU(1) = O0 ETAGAU(2) = O0 ETAGAU(3) = O1 POIGAU(1) = O16 POIGAU(2) = O16 POIGAU(3) = O16 C COORDONNEES DES 3 POINTS DE GAUSS ELSE IF (ISUPTH.EQ.2) THEN QSIGAU(1) = O16 QSIGAU(2) = O16 QSIGAU(3) = O66 ETAGAU(1) = O16 ETAGAU(2) = O66 ETAGAU(3) = O16 POIGAU(1) = O16 POIGAU(2) = O16 POIGAU(3) = O16 C COORDONNEES DU CENTRE DE GRAVITE ELSE IF (ISUPTH.EQ.3) THEN QSIGAU(1) = O1s3 ETAGAU(1) = O1s3 POIGAU(1) = O1s2 ENDIF * FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S) DO I = 1, NBPGAU XI = QSIGAU(I) YI = ETAGAU(I) SHPTOT(1,1,I) = O1-XI-YI SHPTOT(1,2,I) = XI SHPTOT(1,3,I) = YI SHPTOT(2,1,I) = -O1 SHPTOT(2,2,I) = O1 SHPTOT(2,3,I) = O0 SHPTOT(3,1,I) = -O1 SHPTOT(3,2,I) = O0 SHPTOT(3,3,I) = O1 ENDDO IPINTE = MINTE C RETURN END