1 : C TTET4A SOURCE CHAT 05/01/13 03:54:25 5004 2 : SUBROUTINE TTET4A(POINTS,QSIGAU,ETAGAU,DZEGAU,POIGAU,SHPTOT) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : ************************************************************************ 6 : * 7 : * T T E T 4 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 TET4 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 : * DZEGAU (S) | 29 : * POIGAU (S) POIDS D'INTEGRATION,SI ON A 'GAUSS'. 30 : * SHPTOT (S) FONCTIONS DE FORME ET LEURS DERIVEES. 31 : * 32 : REAL*8 QSIGAU(5),ETAGAU(5),DZEGAU(5),POIGAU(5) 33 : REAL*8 SHPTOT(6,4,5) 34 : CHARACTER*(*) POINTS 35 : * 36 : * CONSTANTES: 37 : * ----------- 38 : * 39 : * 40 : PARAMETER ( UN=1.D0 ) 41 : PARAMETER ( XZER=0.D0 ) 42 : PARAMETER ( UNDEMI=.5D0 ) 43 : PARAMETER ( UNQUA=.25D0 ) 44 : PARAMETER ( UNSIX=1.D0/6.D0 ) 45 : PARAMETER ( OHUIT=-.80D0 ) 46 : PARAMETER ( OQUCI= .45D0 ) 47 : PARAMETER ( UNDOU=1.D0/12.D0 ) 48 : * 49 : * 50 : * AUTEUR, DATE DE CREATION: 51 : * ------------------------- 52 : * 53 : * DENIS ROBERT,LE 16 FEVRIER 1988. 54 : * 55 : * LANGAGE: 56 : * -------- 57 : * 58 : * FORTRAN77 59 : * 60 : ************************************************************************ 61 : * 62 : IF (POINTS(1:5).EQ.'NOEUD') GOTO 66 63 : IF (POINTS(1:5).NE.'GAUSS') GOTO 200 64 : * 65 : * COORDONNEES DES 5 POINTS DE GAUSS 66 : * 67 : NBPGAU=5 68 : NBNO=10 69 : QSIGAU(1)=UNQUA 70 : QSIGAU(2)=UNSIX 71 : QSIGAU(3)=UNSIX 72 : QSIGAU(4)=UNSIX 73 : QSIGAU(5)=UNDEMI 74 : ETAGAU(1)=UNQUA 75 : ETAGAU(2)=UNSIX 76 : ETAGAU(3)=UNSIX 77 : ETAGAU(4)=UNDEMI 78 : ETAGAU(5)=UNSIX 79 : DZEGAU(1)=UNQUA 80 : DZEGAU(2)=UNSIX 81 : DZEGAU(3)=UNDEMI 82 : DZEGAU(4)=UNSIX 83 : DZEGAU(5)=UNSIX 84 : POIGAU(1)=OHUIT*UNSIX 85 : POIGAU(2)=OQUCI*UNSIX 86 : POIGAU(3)=OQUCI*UNSIX 87 : POIGAU(4)=OQUCI*UNSIX 88 : POIGAU(5)=OQUCI*UNSIX 89 : GOTO 2 90 : * 91 : * COORDONNEES DES NOEUDS 92 : * 93 : 66 CONTINUE 94 : NBNO=4 95 : NBPGAU=NBNO 96 : QSIGAU(2)=UN 97 : ETAGAU(3)=UN 98 : DZEGAU(4)=UN 99 : DO 241 IA=1,4 100 : POIGAU(IA)=UNDOU 101 : 241 CONTINUE 102 : * 103 : * FONCTIONS DE FORME ET LEURS DERIVEES 104 : * 105 : 2 CONTINUE 106 : DO 1 I=1,NBPGAU 107 : AUX=UN-QSIGAU(I)-ETAGAU(I)-DZEGAU(I) 108 : SHPTOT(1,1,I)=AUX 109 : SHPTOT(1,2,I)=QSIGAU(I) 110 : SHPTOT(1,3,I)=ETAGAU(I) 111 : SHPTOT(1,4,I)=DZEGAU(I) 112 : SHPTOT(2,1,I)=-UN 113 : SHPTOT(2,2,I)= UN 114 : SHPTOT(2,3,I)= XZER 115 : SHPTOT(2,4,I)= XZER 116 : SHPTOT(3,1,I)=-UN 117 : SHPTOT(3,2,I)= XZER 118 : SHPTOT(3,3,I)= UN 119 : SHPTOT(3,4,I)= XZER 120 : SHPTOT(4,1,I)=-UN 121 : SHPTOT(4,2,I)= XZER 122 : SHPTOT(4,3,I)= XZER 123 : SHPTOT(4,4,I)= UN 124 : 1 CONTINUE 125 : RETURN 126 : 200 CONTINUE 126 : CALL ERREUR(19) 127 : END 128 :
© Cast3M 2003 - All rights reserved.
Disclaimer