1 : C TTE10A SOURCE CHAT 05/01/13 03:54:13 5004 2 : SUBROUTINE TTE10A(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 1 0 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 TE10 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(*),ETAGAU(*),DZEGAU(*),POIGAU(*) 33 : REAL*8 SHPTOT(6,10,*) 34 : CHARACTER*(*) POINTS 35 : * 36 : * CONSTANTES: 37 : * ----------- 38 : * 39 : * 40 : PARAMETER ( DEUX=2.D0 ) 41 : PARAMETER ( UN=1.D0 ) 42 : PARAMETER ( QUATRE=4.D0 ) 43 : PARAMETER ( XZER=0.D0 ) 44 : PARAMETER ( UNDEMI=.5D0 ) 45 : PARAMETER ( UNQUA=.25D0 ) 46 : PARAMETER ( UNSIX=.166666666666666666D0 ) 47 : PARAMETER ( OHUIT=-.80D0 ) 48 : PARAMETER ( OQUCI= .45D0 ) 49 : PARAMETER ( UN30=1.D0/30.D0 ) 50 : * 51 : * 52 : * AUTEUR, DATE DE CREATION: 53 : * ------------------------- 54 : * 55 : * DENIS ROBERT,LE 21 JANVIER 1988. 56 : * 57 : * LANGAGE: 58 : * -------- 59 : * 60 : * FORTRAN77 61 : * 62 : ***' tte10a ' 63 : ********************************************************************* 64 : * 65 : IF (POINTS(1:5).EQ.'NOEUD') GOTO 66 66 : IF (POINTS(1:5).NE.'GAUSS') GOTO 200 67 : * 68 : * COORDONNEES DES 5 POINTS DE GAUSS 69 : * 70 : NBPGAU=5 71 : NBNO=10 72 : QSIGAU(1)=UNQUA 73 : QSIGAU(2)=UNSIX 74 : QSIGAU(3)=UNSIX 75 : QSIGAU(4)=UNSIX 76 : QSIGAU(5)=UNDEMI 77 : ETAGAU(1)=UNQUA 78 : ETAGAU(2)=UNSIX 79 : ETAGAU(3)=UNSIX 80 : ETAGAU(4)=UNDEMI 81 : ETAGAU(5)=UNSIX 82 : DZEGAU(1)=UNQUA 83 : DZEGAU(2)=UNSIX 84 : DZEGAU(3)=UNDEMI 85 : DZEGAU(4)=UNSIX 86 : DZEGAU(5)=UNSIX 87 : POIGAU(1)=OHUIT*UNSIX 88 : POIGAU(2)=OQUCI*UNSIX 89 : POIGAU(3)=OQUCI*UNSIX 90 : POIGAU(4)=OQUCI*UNSIX 91 : POIGAU(5)=OQUCI*UNSIX 92 : GOTO 2 93 : * 94 : * COORDONNEES DES NOEUDS 95 : * 96 : 66 CONTINUE 97 : C OPTION NON IMPLEMENTEE 97 : C CALL ERREUR (251) 98 : C RETURN 99 : NBNO=10 100 : NBPGAU=NBNO 101 : QSIGAU(2)=UNDEMI 102 : QSIGAU(4)=UNDEMI 103 : QSIGAU(8)=UNDEMI 104 : ETAGAU(4)=UNDEMI 105 : ETAGAU(6)=UNDEMI 106 : ETAGAU(9)=UNDEMI 107 : DZEGAU(7)=UNDEMI 108 : DZEGAU(8)=UNDEMI 109 : DZEGAU(9)=UNDEMI 110 : QSIGAU(3)=UN 111 : ETAGAU(5)=UN 112 : DZEGAU(10)=UN 113 : C DO 241 IA=1,10 114 : C POIGAU(IA)=UN30 115 : C 241 CONTINUE 116 : * 117 : * FONCTIONS DE FORME ET LEURS DERIVEES 118 : * 119 : 2 CONTINUE 120 : DO 1 I=1,NBPGAU 121 : AUX=UN-QSIGAU(I)-ETAGAU(I)-DZEGAU(I) 122 : SHPTOT(1,1,I)=AUX*(DEUX*AUX-UN) 123 : SHPTOT(1,2,I)=QUATRE*QSIGAU(I)*AUX 124 : SHPTOT(1,3,I)=QSIGAU(I)*(DEUX*QSIGAU(I)-UN) 125 : SHPTOT(1,4,I)=QUATRE*QSIGAU(I)*ETAGAU(I) 126 : SHPTOT(1,5,I)=ETAGAU(I)*(DEUX*ETAGAU(I)-UN) 127 : SHPTOT(1,6,I)=QUATRE*ETAGAU(I)*AUX 128 : SHPTOT(1,7,I)=QUATRE*DZEGAU(I)*AUX 129 : SHPTOT(1,8,I)=QUATRE*QSIGAU(I)*DZEGAU(I) 130 : SHPTOT(1,9,I)=QUATRE*ETAGAU(I)*DZEGAU(I) 131 : SHPTOT(1,10,I)=DZEGAU(I)*(DEUX*DZEGAU(I)-UN) 132 : SHPTOT(2,1,I)=UN-QUATRE*AUX 133 : SHPTOT(2,2,I)=QUATRE*(AUX-QSIGAU(I)) 134 : SHPTOT(2,3,I)=QUATRE*QSIGAU(I)-UN 135 : SHPTOT(2,4,I)=QUATRE*ETAGAU(I) 136 : SHPTOT(2,5,I)=XZER 137 : SHPTOT(2,6,I)=-QUATRE*ETAGAU(I) 138 : SHPTOT(2,7,I)=-QUATRE*DZEGAU(I) 139 : SHPTOT(2,8,I)= QUATRE*DZEGAU(I) 140 : SHPTOT(2,9,I)= XZER 141 : SHPTOT(2,10,I)=XZER 142 : SHPTOT(3,1,I)= UN-QUATRE*AUX 143 : SHPTOT(3,2,I)= -QUATRE*QSIGAU(I) 144 : SHPTOT(3,3,I)= XZER 145 : SHPTOT(3,4,I)= QUATRE*QSIGAU(I) 146 : SHPTOT(3,5,I)= QUATRE*ETAGAU(I)-UN 147 : SHPTOT(3,6,I)= QUATRE*(AUX-ETAGAU(I)) 148 : SHPTOT(3,7,I)=-QUATRE*DZEGAU(I) 149 : SHPTOT(3,8,I)= XZER 150 : SHPTOT(3,9,I)= QUATRE*DZEGAU(I) 151 : SHPTOT(3,10,I)=XZER 152 : SHPTOT(4,1,I)= UN-QUATRE*AUX 153 : SHPTOT(4,2,I)= -QUATRE*QSIGAU(I) 154 : SHPTOT(4,3,I)= XZER 155 : SHPTOT(4,4,I)= XZER 156 : SHPTOT(4,5,I)= XZER 157 : SHPTOT(4,6,I)= -QUATRE*ETAGAU(I) 158 : SHPTOT(4,7,I)=QUATRE*(AUX-DZEGAU(I)) 159 : SHPTOT(4,8,I)=QUATRE*QSIGAU(I) 160 : SHPTOT(4,9,I)=QUATRE*ETAGAU(I) 161 : SHPTOT(4,10,I)=QUATRE*DZEGAU(I)-UN 162 : 1 CONTINUE 163 : RETURN 164 : 200 CONTINUE 164 : CALL ERREUR(19) 165 : END 166 : 167 :
© Cast3M 2003 - All rights reserved.
Disclaimer