1 : C SHAPE SOURCE PV 08/09/11 21:17:04 6150 2 : SUBROUTINE SHAPE(QSI,ETA,DZE,IELE,SHP,IRET) 3 : C======================================================================= 4 : C 5 : C FONCTIONS DE FORME ET LEUR DERIVEES 6 : C 7 : C======================================================================= 8 : IMPLICIT INTEGER(I-N) 9 : IMPLICIT REAL*8(A-H,O-Z) 10 : PARAMETER(XZER=0.D0,UNDEMI=.5D0,UN=1.D0,DEUX=2.D0) 11 : PARAMETER(TROIS=3.D0,QUATRE=4.D0) 12 : DIMENSION SHP(6,*) 13 : C======================================================================= 14 : C QSI ETA DZE COORDONNEES REDUITES DU POINT DE GAUSS 15 : C 16 : C SHP FONCTIONS DE FORME SHP(1,INO) FONCTIONS DE FORME DU ND INO 17 : C SHP(2,INO) DERIVEES PR RPPRT A QSI 18 : C SHP(3,INO) DERIVEES PR RPPRT A ETA 19 : C SHP(4,INO) DERIVEES PR RPPRT A DZE 20 : C SHP(5,INO) DERIVEES PR RPPRT A 4 21 : C COORDONNEES BARYCENTRIQUES 22 : C CAS DE DIMENSION 1 23 : C 2 OU 3 NOEUDS LIGNE 24 : C 2 NOEUDS POLYNOME DE HERMITE ET DERIVEES 25 : C CAS DE DIMENSION 2 26 : C 3 OU 6 NOEUDS TRIANGLE 27 : C 7 NOEUDS TRIANGLE 28 : C 4 OU 8 NOEUDS QUADRILATERE 29 : C 9 NOEUDS QUADRILATERE 30 : C CAS DE DIMENSION 3 FAIT DANS SHAPE3 ETIQUETTE 50 31 : C 4 OU 10 NOEUDS TETRAEDRE 32 : C 5 OU 13 NOEUDS PYRAMIDE 33 : C 6 OU 15 NOEUDS PRISME 34 : C 8 OU 20 NOEUDS CUBE 35 : C CAS DE DIMENSION 3 QUAF FAIT DANS SHAPEQ ETIQUETTE 51 36 : C 15 NOEUDS TETRAEDRE 37 : C 19 NOEUDS PYRAMIDE 38 : C 21 NOEUDS PRISME 39 : C 27 NOEUDS CUBE 40 : C======================================================================= 41 : IRET=1 42 : C 43 : C 66 INDIQUE QUE L ELEMENT PAS ENCORE IMPLEMENTE 44 : C 45 : C P1 S2 S3 T3 T4 T6 T7 Q4 Q5 Q8 Q9 R2 R3 C8 C20 P6 P15 46 : GOTO ( 1, 2, 3, 4,66, 6,7, 8,66,10,10,12,13,50, 50,50, 50, 47 : C L3 L4 L6 L8 MU T4 T10 P5 P13 ATT SUP RP3 LP6 LP8 POLY 48 : 1 18,19,66,66,66,50, 50,50, 50, 66, 66, 29, 30, 31, 66, 49 : C CU27 PR21 TE15 PY19 50 : 1 51,51,51,51 ) IELE 51 : if(iele.eq.260) then 52 : write(6,*) ' on passe dans shape bizarre' 53 : endif 54 : C 55 : 66 IRET=0 56 : GOTO 666 57 : C 58 : C FONCTIONS DE FORME MONODIMENSIONELLES 59 : C 60 : C 61 : C FONCTION FORME 1 NOEUD 62 : C 63 : 1 CONTINUE 64 : SHP(1,1)=UN 65 : GOTO 666 66 : C FONCTIONS DE FORME 2 NOEUDS 67 : C 68 : 2 CONTINUE 69 : SHP(1,1)=(UN-QSI)/DEUX 70 : SHP(1,2)=(UN+QSI)/DEUX 71 : SHP(2,1)=-UNDEMI 72 : SHP(2,2)= UNDEMI 73 : GOTO 666 74 : C 75 : C FONCTIONS DE FORME 3 NOEUDS 76 : C 77 : 3 CONTINUE 78 : SHP(1,1)=(UN-QSI)/DEUX 79 : SHP(1,3)=(UN+QSI)/DEUX 80 : SHP(2,1)=-UNDEMI 81 : SHP(2,3)= UNDEMI 82 : SHP(1,2)=UN-QSI*QSI 83 : SHP(2,2)=-DEUX*QSI 84 : SHP(1,1)=SHP(1,1)-SHP(1,2)/DEUX 85 : SHP(1,3)=SHP(1,3)-SHP(1,2)/DEUX 86 : SHP(2,1)=SHP(2,1)-SHP(2,2)/DEUX 87 : SHP(2,3)=SHP(2,3)-SHP(2,2)/DEUX 88 : GOTO 666 89 : C 90 : C POLYNOME DE HERMITE 2 NOEUDS 1 FCTN DE FORME ET SA DERIVEE 91 : C 92 : C104 SHP(1,1)=(UN-QSI)*(UN-QSI)*(DEUX+QSI)/QUATRE 93 : C SHP(2,1)= TROIS*(QSI*QSI-UN)/QUATRE 94 : C SHP(1,2)=(UN-QSI*QSI)*(UN-QSI)/QUATRE 95 : C SHP(2,2)=(QSI-UN)*(UN+TROIS*QSI)/QUATRE 96 : C SHP(1,3)=(UN+QSI)*(UN+QSI)*(DEUX-QSI)/QUATRE 97 : C SHP(2,3)=TROIS*(UN-QSI*QSI)/QUATRE 98 : C SHP(1,4)=(QSI*QSI-UN)*(UN+QSI)/QUATRE 99 : C SHP(2,4)=(UN+QSI)*(TROIS*QSI-UN)/QUATRE 100 : C GOTO 666 101 : C 102 : C FONCTIONS DE FORME BIDIMENSIONNELLES 103 : C 104 : 4 CONTINUE 105 : C 106 : C FONCTIONS DE FORME TRIANGLE A 3 NOEUDS 107 : C 108 : AUX=UN-QSI-ETA 109 : SHP(1,1)=AUX 110 : SHP(1,2)=QSI 111 : SHP(1,3)=ETA 112 : SHP(2,1)=-UN 113 : SHP(2,2)=UN 114 : SHP(2,3)=XZER 115 : SHP(3,1)=-UN 116 : SHP(3,2)=XZER 117 : SHP(3,3)=UN 118 : GOTO 666 119 : 6 CONTINUE 120 : C 121 : C FONCTIONS DE FORME TRIANGLE A 6 NOEUDS 122 : C 123 : AUX=UN-QSI-ETA 124 : SHP(1,1)=AUX 125 : SHP(1,3)=QSI 126 : SHP(1,5)=ETA 127 : SHP(2,1)=-UN 128 : SHP(2,3)=UN 129 : SHP(2,5)=XZER 130 : SHP(3,1)=-UN 131 : SHP(3,3)=XZER 132 : SHP(3,5)=UN 133 : C AJOUTER NOEUD 2 134 : SHP(1,2)=QUATRE*QSI*AUX 135 : SHP(2,2)=QUATRE*(AUX-QSI) 136 : SHP(3,2)=-QUATRE*QSI 137 : DO 1010 I=1,3 138 : SHP(I,1)=SHP(I,1)-SHP(I,2)/DEUX 139 : 1010 SHP(I,3)=SHP(I,3)-SHP(I,2)/DEUX 140 : C AJOUTER NOEUD 4 141 : SHP(1,4)=QUATRE*QSI*ETA 142 : SHP(2,4)=QUATRE*ETA 143 : SHP(3,4)=QUATRE*QSI 144 : DO 1020 I=1,3 145 : SHP(I,3)=SHP(I,3)-SHP(I,4)/DEUX 146 : 1020 SHP(I,5)=SHP(I,5)-SHP(I,4)/DEUX 147 : C AJOUTER NOEUD 6 148 : SHP(1,6)=QUATRE*ETA*AUX 149 : SHP(2,6)=-QUATRE*ETA 150 : SHP(3,6)=QUATRE*(AUX-ETA) 151 : DO 1030 I=1,3 152 : SHP(I,5)=SHP(I,5)-SHP(I,6)/DEUX 153 : 1030 SHP(I,1)=SHP(I,1)-SHP(I,6)/DEUX 154 : GOTO 666 155 : 7 CONTINUE 156 : C 157 : C FONCTIONS DE FORME TRIANGLE A 7 NOEUDS 158 : C 159 : AUX=UN-QSI-ETA 160 : SHP(1,1)=AUX*(DEUX*AUX-UN)+TROIS*AUX*QSI*ETA 161 : SHP(1,3)=QSI*(DEUX*QSI-UN)+TROIS*AUX*QSI*ETA 162 : SHP(1,5)=ETA*(DEUX*ETA-UN)+TROIS*AUX*QSI*ETA 163 : 164 : SHP(1,2)=QUATRE*(AUX*QSI-TROIS*AUX*QSI*ETA) 165 : SHP(1,4)=QUATRE*(QSI*ETA-TROIS*AUX*QSI*ETA) 166 : SHP(1,6)=QUATRE*(AUX*ETA-TROIS*AUX*QSI*ETA) 167 : 168 : SHP(1,7)=TROIS*TROIS*TROIS*AUX*QSI*ETA 169 : 170 : C GRADIENTS DES FONCTIONS DE FORME A 7 NOEUDS 171 : C NOEUD 1 172 : SHP(2,1)=-3.D0*ETA*ETA+4.D0*QSI+7.D0*ETA- 173 : * 6.D0*QSI*ETA-3.D0 174 : SHP(3,1)=-3.D0*QSI*QSI+4.D0*ETA+7.D0*QSI- 175 : * 6.D0*QSI*ETA-3.D0 176 : C NOEUD 3 177 : SHP(2,3)=-3.D0*ETA*ETA+4.D0*QSI+3.D0*ETA- 178 : * 6.D0*QSI*ETA-1.D0 179 : SHP(3,3)=-3.D0*QSI*QSI+3.D0*QSI-6.D0*QSI*ETA 180 : C NOEUD 5 181 : SHP(2,5)=-3.D0*ETA*ETA+3.D0*ETA-6.D0*QSI*ETA 182 : SHP(3,5)=-3.D0*QSI*QSI+4.D0*ETA+3.D0*QSI- 183 : * 6.D0*QSI*ETA-1.D0 184 : C NOEUD 2 185 : SHP(2,2)=QUATRE*TROIS*ETA*ETA-QUATRE*QUATRE*ETA-DEUX*QUATRE*QSI+ 186 : * DEUX*TROIS*QUATRE*QSI*ETA+QUATRE 187 : SHP(3,2)=QUATRE*TROIS*QSI*QSI-QUATRE*QUATRE*QSI+ 188 : * DEUX*TROIS*QUATRE*QSI*ETA 189 : C NOEUD 4 190 : SHP(2,4)=QUATRE*TROIS*ETA*ETA-DEUX*QUATRE*ETA+ 191 : * DEUX*TROIS*QUATRE*QSI*ETA 192 : SHP(3,4)=QUATRE*TROIS*QSI*QSI-DEUX*QUATRE*QSI+ 193 : * DEUX*TROIS*QUATRE*QSI*ETA 194 : C NOEUD 6 195 : SHP(2,6)=QUATRE*TROIS*ETA*ETA-QUATRE*QUATRE*ETA+ 196 : * DEUX*TROIS*QUATRE*QSI*ETA 197 : SHP(3,6)=QUATRE*TROIS*QSI*QSI-QUATRE*QUATRE*QSI-DEUX*QUATRE*ETA+ 198 : * DEUX*TROIS*QUATRE*QSI*ETA+QUATRE 199 : C NOEUD 7 200 : SHP(2,7)=TROIS*TROIS*TROIS*ETA*(UN-DEUX*QSI-ETA) 201 : SHP(3,7)=TROIS*TROIS*TROIS*QSI*(UN-QSI-DEUX*ETA) 202 : GOTO 666 203 : 8 CONTINUE 204 : C 205 : C FONCTIONS DE FORME QUADRANGLE A 4 NOEUDS 206 : C 207 : SHP(1,1)=(UN-QSI)*(UN-ETA)/QUATRE 208 : SHP(1,2)=(UN+QSI)*(UN-ETA)/QUATRE 209 : SHP(1,3)=(UN+QSI)*(UN+ETA)/QUATRE 210 : SHP(1,4)=(UN-QSI)*(UN+ETA)/QUATRE 211 : SHP(2,1)=(ETA-UN)/QUATRE 212 : SHP(2,2)=-SHP(2,1) 213 : SHP(2,3)=(ETA+UN)/QUATRE 214 : SHP(2,4)=-SHP(2,3) 215 : SHP(3,1)=(QSI-UN)/QUATRE 216 : SHP(3,2)=-(UN+QSI)/QUATRE 217 : SHP(3,3)=-SHP(3,2) 218 : SHP(3,4)=-SHP(3,1) 219 : GOTO 666 220 : C 221 : C FONCTIONS DE FORME QUADRANGLE A 8 NOEUDS 222 : C 223 : 10 CONTINUE 224 : SHP(1,1)=(UN-QSI)*(UN-ETA)/QUATRE 225 : SHP(1,3)=(UN+QSI)*(UN-ETA)/QUATRE 226 : SHP(1,5)=(UN+QSI)*(UN+ETA)/QUATRE 227 : SHP(1,7)=(UN-QSI)*(UN+ETA)/QUATRE 228 : SHP(2,1)=(ETA-UN)/QUATRE 229 : SHP(2,3)=-SHP(2,1) 230 : SHP(2,5)=(ETA+UN)/QUATRE 231 : SHP(2,7)=-SHP(2,5) 232 : SHP(3,1)=(QSI-UN)/QUATRE 233 : SHP(3,3)=-(UN+QSI)/QUATRE 234 : SHP(3,5)=-SHP(3,3) 235 : SHP(3,7)=-SHP(3,1) 236 : C AJOUTER NOEUD 2 237 : SHP(1,2)=(UN-QSI*QSI)*(UN-ETA)/DEUX 238 : SHP(2,2)=(ETA-UN)*QSI 239 : SHP(3,2)=(QSI*QSI-1)/DEUX 240 : DO 2010 I=1,3 241 : SHP(I,1)=SHP(I,1)-SHP(I,2)/DEUX 242 : 2010 SHP(I,3)=SHP(I,3)-SHP(I,2)/DEUX 243 : C AJOUTER NOEUD 4 244 : SHP(1,4)=(UN-ETA*ETA)*(UN+QSI)/DEUX 245 : SHP(2,4)=(UN-ETA*ETA)/DEUX 246 : SHP(3,4)=-ETA*(UN+QSI) 247 : DO 2020 I=1,3 248 : SHP(I,3)=SHP(I,3)-SHP(I,4)/DEUX 249 : 2020 SHP(I,5)=SHP(I,5)-SHP(I,4)/DEUX 250 : C AJOUTER NOEUD 6 251 : SHP(1,6)=(UN-QSI*QSI)*(UN+ETA)/DEUX 252 : SHP(2,6)=-QSI*(UN+ETA) 253 : SHP(3,6)=(UN-QSI*QSI)/DEUX 254 : DO 2030 I=1,3 255 : SHP(I,5)=SHP(I,5)-SHP(I,6)/DEUX 256 : 2030 SHP(I,7)=SHP(I,7)-SHP(I,6)/DEUX 257 : C AJOUTER NOEUD 8 258 : SHP(1,8)=(UN-ETA*ETA)*(UN-QSI)/DEUX 259 : SHP(2,8)=(ETA*ETA-UN)/DEUX 260 : SHP(3,8)=ETA*(QSI-UN) 261 : DO 2040 I=1,3 262 : SHP(I,7)=SHP(I,7)-SHP(I,8)/DEUX 263 : 2040 SHP(I,1)=SHP(I,1)-SHP(I,8)/DEUX 264 : C 265 : C FONCTIONS DE FORME QUADRANGLE A 9 NOEUDS 266 : C 267 : IF (IELE.NE.11) GOTO 666 268 : SHP(1,9)=(UN-QSI*QSI)*(UN-ETA*ETA) 269 : SHP(2,9)=DEUX*ETA*(QSI*QSI-UN) 270 : SHP(3,9)=DEUX*QSI*(ETA*ETA-UN) 271 : DO 2051 I=1,3 272 : DO 2051 J=1,7,2 273 : SHP(I,J)=SHP(I,J)+SHP(I,9)/QUATRE 274 : 2051 SHP(I,J+1)=SHP(I,J+1)-SHP(I,9)/DEUX 275 : GOTO 666 276 : C 277 : C FONCTIONS DE FORME ELEMENTS DE RACCORD 278 : C LIQUIDE - MASSIF LINEAIRE 279 : C 280 : C FONCTIONS DE FORME RAC2 281 : C 282 : 12 CONTINUE 283 : SHP(1,1)=(1.0-QSI)/DEUX 284 : SHP(1,2)=(1.0+QSI)/DEUX 285 : SHP(2,1)=-1.0/DEUX 286 : SHP(2,2)= 1.0/DEUX 287 : SHP(1,3)=SHP(1,2) 288 : SHP(1,4)=SHP(1,1) 289 : SHP(2,3)=SHP(2,2) 290 : SHP(2,4)=SHP(2,1) 291 : GOTO 666 292 : C 293 : C FONCTIONS DE FORME RAC3 294 : C 295 : 13 CONTINUE 296 : SHP(1,1)=(-QSI+QSI*QSI)*(UN-ETA)/QUATRE 297 : SHP(1,2)=( UN -QSI*QSI)*(UN-ETA)/DEUX 298 : SHP(1,3)=( QSI+QSI*QSI)*(UN-ETA)/QUATRE 299 : SHP(1,4)=( QSI+QSI*QSI)*(UN+ETA)/QUATRE 300 : SHP(1,5)=( UN -QSI*QSI)*(UN+ETA)/DEUX 301 : SHP(1,6)=(-QSI+QSI*QSI)*(UN+ETA)/QUATRE 302 : C 303 : SHP(2,1)=(-UN+DEUX*QSI)*(UN-ETA)/QUATRE 304 : SHP(2,2)=( -DEUX*QSI)*(UN-ETA)/DEUX 305 : SHP(2,3)=( UN+DEUX*QSI)*(UN-ETA)/QUATRE 306 : SHP(2,4)=( UN+DEUX*QSI)*(UN+ETA)/QUATRE 307 : SHP(2,5)=( -DEUX*QSI)*(UN+ETA)/DEUX 308 : SHP(2,6)=(-UN+DEUX*QSI)*(UN+ETA)/QUATRE 309 : C 310 : SHP(3,1)=-(-QSI+QSI*QSI)/QUATRE 311 : SHP(3,2)=-( UN -QSI*QSI)/DEUX 312 : SHP(3,3)=-( QSI+QSI*QSI)/QUATRE 313 : SHP(3,4)= ( QSI+QSI*QSI)/QUATRE 314 : SHP(3,5)= ( UN -QSI*QSI)/DEUX 315 : SHP(3,6)= (-QSI+QSI*QSI)/QUATRE 316 : GOTO 666 317 : 18 CONTINUE 318 : C 319 : C FONCTIONS DE FORME LIA3 320 : C 321 : AUX=UN-QSI-ETA 322 : SHP(1,1)=AUX 323 : SHP(1,2)=QSI 324 : SHP(1,3)=ETA 325 : SHP(2,1)=-UN 326 : SHP(2,2)=UN 327 : SHP(2,3)=XZER 328 : SHP(3,1)=-UN 329 : SHP(3,2)=XZER 330 : SHP(3,3)=UN 331 : SHP(1,4)=SHP(1,1) 332 : SHP(1,5)=SHP(1,2) 333 : SHP(1,6)=SHP(1,3) 334 : SHP(2,4)=SHP(2,1) 335 : SHP(2,5)=SHP(2,2) 336 : SHP(2,6)=SHP(2,3) 337 : SHP(3,4)=SHP(3,1) 338 : SHP(3,5)=SHP(3,2) 339 : SHP(3,6)=SHP(3,3) 340 : GOTO 666 341 : 19 CONTINUE 342 : C 343 : C FONCTIONS DE FORME LIA4 344 : C 345 : SHP(1,1)=(UN-QSI)*(UN-ETA)/QUATRE 346 : SHP(1,2)=(UN+QSI)*(UN-ETA)/QUATRE 347 : SHP(1,3)=(UN+QSI)*(UN+ETA)/QUATRE 348 : SHP(1,4)=(UN-QSI)*(UN+ETA)/QUATRE 349 : SHP(2,1)=(ETA-UN)/QUATRE 350 : SHP(2,2)=-SHP(2,1) 351 : SHP(2,3)=(ETA+UN)/QUATRE 352 : SHP(2,4)=-SHP(2,3) 353 : SHP(3,1)=(QSI-UN)/QUATRE 354 : SHP(3,2)=-(UN+QSI)/QUATRE 355 : SHP(3,3)=-SHP(3,2) 356 : SHP(3,4)=-SHP(3,1) 357 : SHP(1,5)=SHP(1,1) 358 : SHP(1,6)=SHP(1,2) 359 : SHP(1,7)=SHP(1,3) 360 : SHP(1,8)=SHP(1,4) 361 : SHP(2,5)=SHP(2,1) 362 : SHP(2,6)=SHP(2,2) 363 : SHP(2,7)=SHP(2,3) 364 : SHP(2,8)=SHP(2,4) 365 : SHP(3,5)=SHP(3,1) 366 : SHP(3,6)=SHP(3,2) 367 : SHP(3,7)=SHP(3,3) 368 : SHP(3,8)=SHP(3,4) 369 : GOTO 666 370 : 29 CONTINUE 371 : C 372 : C FONCTIONS DE FORME ELEMENTS DE RACCORD POREUX RAP3 373 : C 374 : SHP(1,1)=(-QSI+QSI*QSI)/DEUX 375 : SHP(1,2)=( UN -QSI*QSI) 376 : SHP(1,3)=( QSI+QSI*QSI)/DEUX 377 : SHP(1,4)=( QSI+QSI*QSI)/DEUX 378 : SHP(1,5)=( UN -QSI*QSI) 379 : SHP(1,6)=(-QSI+QSI*QSI)/DEUX 380 : SHP(2,1)=(-UN+DEUX*QSI)/DEUX 381 : SHP(2,2)=( -DEUX*QSI) 382 : SHP(2,3)=( UN+DEUX*QSI)/DEUX 383 : SHP(2,4)=( UN+DEUX*QSI)/DEUX 384 : SHP(2,5)=( -DEUX*QSI) 385 : SHP(2,6)=(-UN+DEUX*QSI)/DEUX 386 : c 387 : SHP(1,7)= ( UN-QSI)/DEUX 388 : SHP(1,8)= ( UN+QSI)/DEUX 389 : SHP(1,9)= ( UN+QSI)/DEUX 390 : SHP(1,10)=( UN-QSI)/DEUX 391 : SHP(1,11)=( UN-QSI)/DEUX 392 : SHP(1,12)=( UN+QSI)/DEUX 393 : SHP(2,7)= -UN/DEUX 394 : SHP(2,8)= UN/DEUX 395 : SHP(2,9)= UN/DEUX 396 : SHP(2,10)=-UN/DEUX 397 : SHP(2,11)=-UN/DEUX 398 : SHP(2,12)= UN/DEUX 399 : GOTO 666 400 : 30 CONTINUE 401 : C 402 : C FONCTIONS DE FORME ELEMENTS DE LIAISON POREUX RAP6 403 : C 404 : SHP(1,1)= (UN-QSI-ETA)*(UN-DEUX*QSI-DEUX*ETA) 405 : SHP(1,2)= QUATRE*QSI*(UN-QSI-ETA) 406 : SHP(1,3)= QSI*(DEUX*QSI-UN) 407 : SHP(1,4)= QUATRE*QSI*ETA 408 : SHP(1,5)= ETA*(DEUX*ETA-UN) 409 : SHP(1,6)= QUATRE*ETA*(UN-QSI-ETA) 410 : SHP(2,1)=-(UN-DEUX*QSI-DEUX*ETA)-DEUX*(UN-QSI-ETA) 411 : SHP(2,2)= QUATRE*(UN-DEUX*QSI-ETA) 412 : SHP(2,3)= QUATRE*QSI-UN 413 : SHP(2,4)= QUATRE*ETA 414 : SHP(2,5)= XZER 415 : SHP(2,6)=-QUATRE*ETA 416 : SHP(3,1)=-(UN-DEUX*QSI-DEUX*ETA)-DEUX*(UN-QSI-ETA) 417 : SHP(3,2)=-QUATRE*QSI 418 : SHP(3,3)= XZER 419 : SHP(3,4)= QUATRE*QSI 420 : SHP(3,5)= QUATRE*ETA-UN 421 : SHP(3,6)= QUATRE*(UN-QSI-DEUX*ETA) 422 : DO 111 I=7,12 423 : SHP(1,I)=SHP(1,I-6) 424 : SHP(2,I)=SHP(2,I-6) 425 : SHP(3,I)=SHP(3,I-6) 426 : 111 CONTINUE 427 : c 428 : SHP(1,13)= UN-QSI-ETA 429 : SHP(1,14)= QSI 430 : SHP(1,15)= ETA 431 : SHP(2,13)=-UN 432 : SHP(2,14)= UN 433 : SHP(2,15)= XZER 434 : SHP(3,13)=-UN 435 : SHP(3,14)= XZER 436 : SHP(3,15)= UN 437 : DO 222 I=16,18 438 : SHP(1,I)=SHP(1,I-3) 439 : SHP(2,I)=SHP(2,I-3) 440 : SHP(3,I)=SHP(3,I-3) 441 : SHP(1,I+3)=SHP(1,I-3) 442 : SHP(2,I+3)=SHP(2,I-3) 443 : SHP(3,I+3)=SHP(3,I-3) 444 : 222 CONTINUE 445 : GOTO 666 446 : 31 CONTINUE 447 : C 448 : C FONCTIONS DE FORME ELEMENTS DE LIAISON POREUX RAP8 449 : C 450 : SHP(1,2)= (UN-QSI*QSI)*(UN-ETA)/DEUX 451 : SHP(1,4)= (UN+QSI)*(UN-ETA*ETA)/DEUX 452 : SHP(1,6)= (UN-QSI*QSI)*(UN+ETA)/DEUX 453 : SHP(1,8)= (UN-QSI)*(UN-ETA*ETA)/DEUX 454 : SHP(1,1)= (UN-QSI)*(UN-ETA)/QUATRE-UN/DEUX*(SHP(1,8)+SHP(1,2)) 455 : SHP(1,3)= (UN+QSI)*(UN-ETA)/QUATRE-UN/DEUX*(SHP(1,2)+SHP(1,4)) 456 : SHP(1,5)= (UN+QSI)*(UN+ETA)/QUATRE-UN/DEUX*(SHP(1,4)+SHP(1,6)) 457 : SHP(1,7)= (UN-QSI)*(UN+ETA)/QUATRE-UN/DEUX*(SHP(1,6)+SHP(1,8)) 458 : SHP(2,2)=-QSI*(UN-ETA) 459 : SHP(2,4)= (UN-ETA*ETA)/DEUX 460 : SHP(2,6)=-QSI*(UN+ETA) 461 : SHP(2,8)=-(UN-ETA*ETA)/DEUX 462 : SHP(2,1)=-(UN-ETA)/QUATRE-UN/DEUX*(SHP(2,8)+SHP(2,2)) 463 : SHP(2,3)= (UN-ETA)/QUATRE-UN/DEUX*(SHP(2,2)+SHP(2,4)) 464 : SHP(2,5)= (UN+ETA)/QUATRE-UN/DEUX*(SHP(2,4)+SHP(2,6)) 465 : SHP(2,7)=-(UN+ETA)/QUATRE-UN/DEUX*(SHP(2,6)+SHP(2,8)) 466 : SHP(3,2)=-(UN-QSI*QSI)/DEUX 467 : SHP(3,4)=-ETA*(UN+QSI) 468 : SHP(3,6)= (UN-QSI*QSI)/DEUX 469 : SHP(3,8)=-ETA*(UN-QSI) 470 : SHP(3,1)=-(UN-QSI)/QUATRE-UN/DEUX*(SHP(3,8)+SHP(3,2)) 471 : SHP(3,3)=-(UN+QSI)/QUATRE-UN/DEUX*(SHP(3,2)+SHP(3,4)) 472 : SHP(3,5)= (UN+QSI)/QUATRE-UN/DEUX*(SHP(3,4)+SHP(3,6)) 473 : SHP(3,7)= (UN-QSI)/QUATRE-UN/DEUX*(SHP(3,6)+SHP(3,8)) 474 : DO 333 I=9,16 475 : SHP(1,I)=SHP(1,I-8) 476 : SHP(2,I)=SHP(2,I-8) 477 : SHP(3,I)=SHP(3,I-8) 478 : 333 CONTINUE 479 : c 480 : SHP(1,17)= (UN-QSI)*(UN-ETA)/QUATRE 481 : SHP(1,18)= (UN+QSI)*(UN-ETA)/QUATRE 482 : SHP(1,19)= (UN+QSI)*(UN+ETA)/QUATRE 483 : SHP(1,20)= (UN-QSI)*(UN+ETA)/QUATRE 484 : SHP(2,17)=-(UN-ETA)/QUATRE 485 : SHP(2,18)= (UN-ETA)/QUATRE 486 : SHP(2,19)= (UN+ETA)/QUATRE 487 : SHP(2,20)=-(UN+ETA)/QUATRE 488 : SHP(3,17)=-(UN-QSI)/QUATRE 489 : SHP(3,18)=-(UN+QSI)/QUATRE 490 : SHP(3,19)= (UN+QSI)/QUATRE 491 : SHP(3,20)= (UN-QSI)/QUATRE 492 : DO 444 I=21,24 493 : SHP(1,I)=SHP(1,I-4) 494 : SHP(2,I)=SHP(2,I-4) 495 : SHP(3,I)=SHP(3,I-4) 496 : SHP(1,I+4)=SHP(1,I-4) 497 : SHP(2,I+4)=SHP(2,I-4) 498 : SHP(3,I+4)=SHP(3,I-4) 499 : 444 CONTINUE 500 : GOTO 666 501 : C 502 : C FONCTIONS DE FORME TRIDIMENSIONNELLES 503 : C 503 : 50 CALL SHAPE3(QSI,ETA,DZE,IELE,SHP,IRT1) 504 : IRET=IRT1 505 : GOTO 666 506 : C 507 : C FONCTIONS DE FORME TRIDIMENSIONNELLES QUAF 508 : C 508 : 51 CALL SHAPEQ(QSI,ETA,DZE,IELE,SHP,IRT1) 509 : IRET=IRT1 510 : C 511 : 666 CONTINUE 512 : RETURN 513 : END 514 : 515 : 516 : 517 : 518 : 519 : 520 : 521 : 522 :
© Cast3M 2003 - All rights reserved.
Disclaimer