1 : C TTRI3A SOURCE CHAT 05/01/13 03:54:37 5004 2 : SUBROUTINE TTRI3A(POINTS,QSIGAU,ETAGAU,POIGAU,SHPTOT) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : ************************************************************************ 6 : * 7 : * T T R I 3 A 8 : * ----------- 9 : * 10 : * FONCTION: ON ACQUIERT: 11 : * --------- - LES COORDONNEES INTRINSEQUES 12 : * - LES FONCTIONS DE FORME ET LEURS DERIVEES 13 : * POUR L'ELEMENT DE THERMIQUE TRI3. 14 : * 15 : * MODULES UTILISES: 16 : * ----------------- 17 : * 18 : -INC CCOPTIO 19 : * 20 : * PARAMETRES: (E)=ENTREE (S)=SORTIE (+ = CONTENU DANS UN COMMUN) 21 : * ----------- 22 : * 23 : * POINTS (E) CHAINE DE CARACTERES: 24 : * 'NOEUDS' SI ON VEUT LES VALEURS AUX NOEUDS 25 : * 'GAUSS' POUR LES VALEURS AUX POINTS DE GAUSS. 26 : * QSIGAU (S) | 27 : * ETAGAU (S) | COORDONNEES INTRINSEQUES. 28 : * POIGAU (S) POIDS D'INTEGRATION,SI ON A 'GAUSS'. 29 : * SHPTOT (S) FONCTIONS DE FORME ET LEURS DERIVEES. 30 : * 31 : REAL*8 QSIGAU(3),ETAGAU(3),POIGAU(3) 32 : REAL*8 SHPTOT(6,3,3) 33 : CHARACTER*(*) POINTS 34 : * 35 : * CONSTANTES: 36 : * ----------- 37 : * 38 : * 39 : PARAMETER ( O0=0.D0 ) 40 : PARAMETER ( O1=1.D0 ) 41 : PARAMETER ( O5=0.5D0 ) 42 : PARAMETER ( O66=0.666666666666666666D0 ) 43 : PARAMETER ( O16=0.166666666666666666D0 ) 44 : * 45 : * 46 : * AUTEUR, DATE DE CREATION: 47 : * ------------------------- 48 : * 49 : * DENIS ROBERT,LE 03 FEVRIER 1988. 50 : * 51 : * LANGAGE: 52 : * -------- 53 : * 54 : * FORTRAN77 55 : * 56 : ************************************************************************ 57 : * 58 : IF (POINTS(1:5).EQ.'NOEUD') GOTO 66 59 : IF (POINTS(1:5).NE.'GAUSS') GOTO 100 60 : * 61 : * COORDONNEES DES 3 POINTS DE GAUSS 62 : * 63 : NBPGAU=3 64 : NBNO=3 65 : QSIGAU(1)=O16 66 : ETAGAU(1)=O16 67 : QSIGAU(2)=O16 68 : ETAGAU(2)=O66 69 : QSIGAU(3)=O66 70 : ETAGAU(3)=O16 71 : POIGAU(1)=O16 72 : POIGAU(2)=O16 73 : POIGAU(3)=O16 74 : GOTO 2 75 : * 76 : * COORDONNEES DES NOEUDS 77 : * 78 : 66 CONTINUE 79 : NBNO=3 80 : NBPGAU=NBNO 81 : QSIGAU(2)=O1 82 : ETAGAU(3)=O1 83 : POIGAU(1)=O16 84 : POIGAU(2)=O16 85 : POIGAU(3)=O16 86 : * 87 : * FONCTIONS DE FORME ET LEURS DERIVEES 88 : * 89 : 2 DO 1 I=1,NBPGAU 90 : XI=QSIGAU(I) 91 : YI=ETAGAU(I) 92 : AUX=O1-XI-YI 93 : SHPTOT(1,1,I)=AUX 94 : SHPTOT(1,2,I)=XI 95 : SHPTOT(1,3,I)=YI 96 : SHPTOT(2,1,I)=-O1 97 : SHPTOT(2,2,I)=O1 98 : SHPTOT(2,3,I)=O0 99 : SHPTOT(3,1,I)=-O1 100 : SHPTOT(3,2,I)=O0 101 : SHPTOT(3,3,I)=O1 102 : * 103 : 1 CONTINUE 104 : RETURN 105 : 100 CONTINUE 105 : CALL ERREUR(19) 106 : * 107 : END 108 :
© Cast3M 2003 - All rights reserved.
Disclaimer