1 : C INTDEL SOURCE CHAT 05/01/13 00:39:49 5004 2 : SUBROUTINE INTDEL(XNTH,XNTB,XNTT,LRN,MELE) 3 : C 4 : C------------------------------------------------------------------------ 5 : C 6 : C CETTE ROUTINE CONTIENT LES INTEGRALS QUI PRESENTENT DES 7 : C FONCTIONS DELTA DE DIRAC DANS LE CAS DES ELEMENTS JOINTS 8 : C POREUX (BALD) 9 : C CETTE ROUTINE EST APPELEE PAR: PERMAB.ESO, RIGI2.ESO, 10 : C EPSI2.ESO, SIGPOR.ESO. 11 : C 12 : C------------------------------------------------------------------------ 13 : C 14 : C ENTREES 15 : C MELE = NUMERO DE L ELEMENT FINI DANS NOMTP 16 : C LRN = DIMENSIONS DES TABLEAUX 17 : C 18 : C SORTIES 19 : C XNTH = TABLEAU QUI SERT POUR LA MATRICE DE PERMEABILITE 20 : C XNTB = TABLEAU QUI SERT POUR LA MATRICE DE PERMEABILITE 21 : C XNTT = TABLEAU QUI SERT POUR LES MATRICES DE PERMEABILITE 22 : C ET DE RIGIDITE 23 : C 24 : C AUTEUR: BALD 20/01/96 25 : C 26 : C------------------------------------------------------------------------ 27 : C 28 : IMPLICIT INTEGER(I-N) 29 : IMPLICIT REAL*8(A-H,O-Z) 30 : C 31 : -INC CCOPTIO 32 : C 33 : DIMENSION XNTH(LRN,LRN),XNTB(LRN,LRN),XNTT(LRN) 34 : PARAMETER(XZER=0.D0,UN=1.D0) 35 : C 35 : CALL ZERO(XNTH,LRN,LRN) 35 : CALL ZERO(XNTB,LRN,LRN) 35 : CALL ZERO(XNTT,1,LRN) 36 : C 37 : IF(MELE.EQ.108)THEN 38 : * 39 : * CAS DE L ELEMENT JOP3 40 : * 41 : XNTH(1,1)=UN 42 : XNTH(1,2)=UN 43 : XNTH(1,3)=XZER 44 : XNTH(1,4)=XZER 45 : XNTH(1,5)=-UN 46 : XNTH(1,6)=-UN 47 : * 48 : XNTH(2,1)=UN 49 : XNTH(2,2)=UN 50 : XNTH(2,3)=XZER 51 : XNTH(2,4)=XZER 52 : XNTH(2,5)=-UN 53 : XNTH(2,6)=-UN 54 : * 55 : XNTH(5,1)=-UN 56 : XNTH(5,2)=-UN 57 : XNTH(5,3)=XZER 58 : XNTH(5,4)=XZER 59 : XNTH(5,5)=UN 60 : XNTH(5,6)=UN 61 : * 62 : XNTH(6,1)=-UN 63 : XNTH(6,2)=-UN 64 : XNTH(6,3)=XZER 65 : XNTH(6,4)=XZER 66 : XNTH(6,5)=UN 67 : XNTH(6,6)=UN 68 : *** 69 : XNTB(3,1)=XZER 70 : XNTB(3,2)=XZER 71 : XNTB(3,3)=UN 72 : XNTB(3,4)=UN 73 : XNTB(3,5)=-UN 74 : XNTB(3,6)=-UN 75 : * 76 : XNTB(4,1)=XZER 77 : XNTB(4,2)=XZER 78 : XNTB(4,3)=UN 79 : XNTB(4,4)=UN 80 : XNTB(4,5)=-UN 81 : XNTB(4,6)=-UN 82 : * 83 : XNTB(5,1)=XZER 84 : XNTB(5,2)=XZER 85 : XNTB(5,3)=-UN 86 : XNTB(5,4)=-UN 87 : XNTB(5,5)=UN 88 : XNTB(5,6)=UN 89 : * 90 : XNTB(6,1)=XZER 91 : XNTB(6,2)=XZER 92 : XNTB(6,3)=-UN 93 : XNTB(6,4)=-UN 94 : XNTB(6,5)=UN 95 : XNTB(6,6)=UN 96 : *** 97 : XNTT(1)=XZER 98 : XNTT(2)=XZER 99 : XNTT(3)=XZER 100 : XNTT(4)=XZER 101 : XNTT(5)=UN 102 : XNTT(6)=UN 103 : ELSE IF(MELE.EQ.109)THEN 104 : * 105 : * CAS DE L ELEMENT JOP6 106 : * 107 : XNTH(1,1)=UN 108 : XNTH(1,2)=UN 109 : XNTH(1,3)=UN 110 : XNTH(1,4)=XZER 111 : XNTH(1,5)=XZER 112 : XNTH(1,6)=XZER 113 : XNTH(1,7)=-UN 114 : XNTH(1,8)=-UN 115 : XNTH(1,9)=-UN 116 : * 117 : XNTH(7,1)=-UN 118 : XNTH(7,2)=-UN 119 : XNTH(7,3)=-UN 120 : XNTH(7,4)=XZER 121 : XNTH(7,5)=XZER 122 : XNTH(7,6)=XZER 123 : XNTH(7,7)=UN 124 : XNTH(7,8)=UN 125 : XNTH(7,9)=UN 126 : * 127 : DO 1 M=1,7,3 128 : DO 2 K=1,2 129 : DO 3 L=1,9 130 : XNTH(K+M,L)=XNTH(M,L) 131 : 3 CONTINUE 132 : 2 CONTINUE 133 : 1 CONTINUE 134 : *** 135 : XNTB(4,1)=XZER 136 : XNTB(4,2)=XZER 137 : XNTB(4,3)=XZER 138 : XNTB(4,4)=UN 139 : XNTB(4,5)=UN 140 : XNTB(4,6)=UN 141 : XNTB(4,7)=-UN 142 : XNTB(4,8)=-UN 143 : XNTB(4,9)=-UN 144 : * 145 : XNTB(7,1)=XZER 146 : XNTB(7,2)=XZER 147 : XNTB(7,3)=XZER 148 : XNTB(7,4)=-UN 149 : XNTB(7,5)=-UN 150 : XNTB(7,6)=-UN 151 : XNTB(7,7)=UN 152 : XNTB(7,8)=UN 153 : XNTB(7,9)=UN 154 : * 155 : DO 7 M=1,7,3 156 : DO 8 K=1,2 157 : DO 9 L=1,9 158 : XNTB(K+M,L)=XNTB(M,L) 159 : 9 CONTINUE 160 : 8 CONTINUE 161 : 7 CONTINUE 162 : *** 163 : XNTT(1)=XZER 164 : XNTT(2)=XZER 165 : XNTT(3)=XZER 166 : XNTT(4)=XZER 167 : XNTT(5)=XZER 168 : XNTT(6)=XZER 169 : XNTT(7)=UN 170 : XNTT(8)=UN 171 : XNTT(9)=UN 172 : ELSE IF(MELE.EQ.110)THEN 173 : * 174 : * CAS DE L ELEMENT JOP8 175 : * 176 : XNTH(1,1) =UN 177 : XNTH(1,2) =UN 178 : XNTH(1,3) =UN 179 : XNTH(1,4) =UN 180 : XNTH(1,5) =XZER 181 : XNTH(1,6) =XZER 182 : XNTH(1,7) =XZER 183 : XNTH(1,8) =XZER 184 : XNTH(1,9) =-UN 185 : XNTH(1,10)=-UN 186 : XNTH(1,11)=-UN 187 : XNTH(1,12)=-UN 188 : * 189 : XNTH(9,1) =-UN 190 : XNTH(9,2) =-UN 191 : XNTH(9,3) =-UN 192 : XNTH(9,4) =-UN 193 : XNTH(9,5) =XZER 194 : XNTH(9,6) =XZER 195 : XNTH(9,7) =XZER 196 : XNTH(9,8) =XZER 197 : XNTH(9,9) =UN 198 : XNTH(9,10)=UN 199 : XNTH(9,11)=UN 200 : XNTH(9,12)=UN 201 : * 202 : DO 4 M=1,9,4 203 : DO 5 K=1,3 204 : DO 6 L=1,12 205 : XNTH(K+M,L)=XNTH(M,L) 206 : 6 CONTINUE 207 : 5 CONTINUE 208 : 4 CONTINUE 209 : *** 210 : XNTB(5,1) =XZER 211 : XNTB(5,2) =XZER 212 : XNTB(5,3) =XZER 213 : XNTB(5,4) =XZER 214 : XNTB(5,5) =UN 215 : XNTB(5,6) =UN 216 : XNTB(5,7) =UN 217 : XNTB(5,8) =UN 218 : XNTB(5,9) =-UN 219 : XNTB(5,10)=-UN 220 : XNTB(5,11)=-UN 221 : XNTB(5,12)=-UN 222 : * 223 : XNTB(9,1) =XZER 224 : XNTB(9,2) =XZER 225 : XNTB(9,3) =XZER 226 : XNTB(9,4) =XZER 227 : XNTB(9,5) =-UN 228 : XNTB(9,6) =-UN 229 : XNTB(9,7) =-UN 230 : XNTB(9,8) =-UN 231 : XNTB(9,9) =UN 232 : XNTB(9,10)=UN 233 : XNTB(9,11)=UN 234 : XNTB(9,12)=UN 235 : * 236 : DO 10 M=1,9,4 237 : DO 11 K=1,3 238 : DO 12 L=1,12 239 : XNTB(K+M,L)=XNTB(M,L) 240 : 12 CONTINUE 241 : 11 CONTINUE 242 : 10 CONTINUE 243 : *** 244 : XNTT(1) =XZER 245 : XNTT(2) =XZER 246 : XNTT(3) =XZER 247 : XNTT(4) =XZER 248 : XNTT(5) =XZER 249 : XNTT(6) =XZER 250 : XNTT(7) =XZER 251 : XNTT(8) =XZER 252 : XNTT(9) =UN 253 : XNTT(10)=UN 254 : XNTT(11)=UN 255 : XNTT(12)=UN 256 : ELSE 256 : CALL ERREUR(5) 257 : END IF 258 : C 259 : RETURN 260 : END 261 : 262 : 263 :
© Cast3M 2003 - All rights reserved.
Disclaimer