1 : C TPRI6A SOURCE CHAT 05/01/13 03:43:32 5004 2 : SUBROUTINE TPRI6A(POINTS,QSIGAU,ETAGAU,DZEGAU,POIGAU,SHPTOT) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : ************************************************************************ 6 : * 7 : * T P R I 6 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 PRI6. 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(6),ETAGAU(6),DZEGAU(6),POIGAU(6) 33 : REAL*8 SHPTOT(6,6,6) 34 : CHARACTER*(*) POINTS 35 : * 36 : * CONSTANTES: 37 : * ----------- 38 : * 39 : * 40 : PARAMETER ( O0=0.D0 ) 41 : PARAMETER ( O1=1.D0 ) 42 : PARAMETER ( O5=0.5D0 ) 43 : PARAMETER ( W=1.D0/6.D0 ) 44 : PARAMETER ( Z=0.577350269189626D0 ) 45 : * 46 : * 47 : * AUTEUR, DATE DE CREATION: 48 : * ------------------------- 49 : * 50 : * J-P.SAUSSAIS OCT.77,REPRIS PAR DENIS ROBERT,LE 26 OCTOBRE 1987. 51 : * 52 : * LANGAGE: 53 : * -------- 54 : * 55 : * FORTRAN77 56 : * 57 : ************************************************************************ 58 : * 59 : IF (POINTS(1:5).EQ.'NOEUD') GOTO 66 60 : IF (POINTS(1:5).NE.'GAUSS') GOTO 200 61 : * 62 : * COORDONNEES DES 6 POINTS DE GAUSS 63 : * 64 : NBPGAU=6 65 : NBNO=6 66 : QSIGAU(1)= O5 67 : ETAGAU(1)= O0 68 : DZEGAU(1)=-Z 69 : QSIGAU(2)= O5 70 : ETAGAU(2)= O5 71 : DZEGAU(2)=-Z 72 : QSIGAU(3)= O0 73 : ETAGAU(3)= O5 74 : DZEGAU(3)=-Z 75 : QSIGAU(4)= O5 76 : ETAGAU(4)= O0 77 : DZEGAU(4)= Z 78 : QSIGAU(5)= O5 79 : ETAGAU(5)= O5 80 : DZEGAU(5)= Z 81 : QSIGAU(6)= O0 82 : ETAGAU(6)= O5 83 : DZEGAU(6)= Z 84 : * 85 : * POIDS 86 : * 87 : DO 10 I=1,NBPGAU 88 : POIGAU(I)=W 89 : 10 CONTINUE 90 : GOTO 2 91 : * 92 : * COORDONNEES INTRINSEQUES DES NOEUDS 93 : * 94 : 66 CONTINUE 95 : NBNO=6 96 : NBPGAU=NBNO 97 : QSIGAU(1)=O0 98 : ETAGAU(1)=O0 99 : DZEGAU(1)=-O1 100 : QSIGAU(2)=O1 101 : ETAGAU(2)=O0 102 : DZEGAU(2)=-O1 103 : QSIGAU(3)=O0 104 : ETAGAU(3)=O1 105 : DZEGAU(3)=-O1 106 : QSIGAU(4)=O0 107 : ETAGAU(4)=O0 108 : DZEGAU(4)=O1 109 : QSIGAU(5)=O1 110 : ETAGAU(5)=O0 111 : DZEGAU(5)=O1 112 : QSIGAU(6)=O0 113 : ETAGAU(6)=O1 114 : DZEGAU(6)=O1 115 : * 116 : * FONCTIONS DE FORME 117 : * 118 : 2 DO 20 I=1,NBPGAU 119 : XI=QSIGAU(I) 120 : YI=ETAGAU(I) 121 : ZI=DZEGAU(I) 122 : * 123 : SHPTOT(1,1,I)=O5*(O1-XI-YI)*(O1-ZI) 124 : SHPTOT(1,2,I)=O5* XI *(O1-ZI) 125 : SHPTOT(1,3,I)=O5* YI *(O1-ZI) 126 : SHPTOT(1,4,I)=O5*(O1-XI-YI)*(O1+ZI) 127 : SHPTOT(1,5,I)=O5* XI *(O1+ZI) 128 : SHPTOT(1,6,I)=O5* YI *(O1+ZI) 129 : * 130 : * DERIVEES EN XI DES FONCTIONS DE FORME 131 : * 132 : SHPTOT(2,1,I)=-O5*(O1-ZI) 133 : SHPTOT(2,2,I)= O5*(O1-ZI) 134 : SHPTOT(2,3,I)= O0 135 : SHPTOT(2,4,I)=-O5*(O1+ZI) 136 : SHPTOT(2,5,I)= O5*(O1+ZI) 137 : SHPTOT(2,6,I)= O0 138 : * 139 : * DERIVEES EN YI DES FONCTIONS DE FORME 140 : * 141 : SHPTOT(3,1,I)=-O5*(O1-ZI) 142 : SHPTOT(3,2,I)= O0 143 : SHPTOT(3,3,I)= O5*(O1-ZI) 144 : SHPTOT(3,4,I)=-O5*(O1+ZI) 145 : SHPTOT(3,5,I)= O0 146 : SHPTOT(3,6,I)= O5*(O1+ZI) 147 : * 148 : * DERIVEES EN ZI DES FONCTIONS DE FORME 149 : * 150 : SHPTOT(4,1,I)=-O5*(O1-XI-YI) 151 : SHPTOT(4,2,I)=-O5*XI 152 : SHPTOT(4,3,I)=-O5*YI 153 : SHPTOT(4,4,I)= O5*(O1-XI-YI) 154 : SHPTOT(4,5,I)= O5*XI 155 : SHPTOT(4,6,I)= O5*YI 156 : 20 CONTINUE 157 : RETURN 158 : 200 CONTINUE 158 : CALL ERREUR(19) 159 : * 160 : END 161 :
© Cast3M 2003 - All rights reserved.
Disclaimer