1 : C SHAPE3 SOURCE CHAT 05/01/13 03:15:24 5004 2 : SUBROUTINE SHAPE3(QSI,ETA,DZE,IELE,SHP,IRET) 3 : C======================================================================= 4 : C 5 : C FONCTIONS DE FORME TRIDIMENSIONNELLES 6 : C EBERSOLT 7 : C======================================================================= 8 : IMPLICIT INTEGER(I-N) 9 : IMPLICIT REAL*8(A-H,O-Z) 10 : DIMENSION SHP(6,*) 11 : DATA XZER/0.D0/,UNDEMI/.5D0/,UN/1.D0/,DEUX/2.D0/,QUATRE/4.D0/ 12 : DATA CINQ/5.D0/,HUIT/8.D0/,TROIS/3.D0/ 13 : DATA UNQUA/.25D0/ 14 : C======================================================================= 15 : C QSI ETA DZE COORDONNEES REDUITES DU POINT DE GAUSS 16 : C 17 : C SHP FONCTIONS DE FORME SHP(1,INO) FONCTIONS DE FORME DU ND INO 18 : C SHP(2,INO) DERIVEES PR RPPRT A QSI 19 : C SHP(3,INO) DERIVEES PR RPPRT A ETA 20 : C SHP(4,INO) DERIVEES PR RPPRT A DZE 21 : C SHP(5,INO) DERIVEES PR RPPRT A 4 22 : C DERIVEES BARYCENTRIQUE 23 : C CAS DE DIMENSION 3 24 : C 4 OU 10 NOEUDS TETRAEDRE 25 : C 5 OU 13 NOEUDS PYRAMIDE 26 : C 6 OU 15 NOEUDS PRISME 27 : C 8 OU 20 NOEUDS CUBE 28 : C======================================================================= 29 : C 30 : IRET=1 31 : C P1 S2 S3 T3 T4 T6 T7 Q4 Q5 Q8 Q9 R2 R3 C8 C20 P6 P15 32 : GOTO (66,66,66,66,66,66,66,66,66,66,66,66,66,14, 15,16, 17, 33 : C L3 L4 L6 L8 MU T4 T10 P5 P13 34 : 1 66,66,66,66,66,23, 24,25, 30) IELE 35 : C 36 : 66 IRET=0 37 : GOTO 666 38 : C 39 : 16 CONTINUE 40 : C 41 : C FONCTIONS DE FORME PRISMATIQUES A 6 NOEUDS 42 : C 43 : AUX=UN-QSI-ETA 44 : AUX1=(UN-DZE)/DEUX 45 : AUX2=(UN+DZE)/DEUX 46 : SHP(1,1)=AUX*AUX1 47 : SHP(1,2)=QSI*AUX1 48 : SHP(1,3)=ETA*AUX1 49 : SHP(2,1)=-AUX1 50 : SHP(2,2)=AUX1 51 : SHP(2,3)=XZER 52 : SHP(3,1)=-AUX1 53 : SHP(3,2)=XZER 54 : SHP(3,3)=AUX1 55 : SHP(4,1)=-AUX/DEUX 56 : SHP(4,2)=-QSI/DEUX 57 : SHP(4,3)=-ETA/DEUX 58 : SHP(1,4)=AUX*AUX2 59 : SHP(1,5)=QSI*AUX2 60 : SHP(1,6)=ETA*AUX2 61 : SHP(2,4)=-AUX2 62 : SHP(2,5)=AUX2 63 : SHP(2,6)=XZER 64 : SHP(3,4)=-AUX2 65 : SHP(3,5)=XZER 66 : SHP(3,6)=AUX2 67 : SHP(4,4)=AUX/DEUX 68 : SHP(4,5)=QSI/DEUX 69 : SHP(4,6)=ETA/DEUX 70 : GOTO 666 71 : 17 CONTINUE 72 : C 73 : C FONCTION DE FORME PRISMATIQUES A 15 NOEUDS 74 : C 75 : AUX=UN-QSI-ETA 76 : DZEM=(UN-DZE) 77 : DZEP=(UN+DZE) 78 : PAUX =QSI-UNDEMI*DZE-UN 79 : PAUX1=ETA-UNDEMI*DZE-UN 80 : PAUX2=QSI+ETA+UNDEMI*DZE 81 : PAUX3=QSI+UNDEMI*DZE-UN 82 : PAUX4=ETA+UNDEMI*DZE-UN 83 : PAUX5=QSI+ETA-UNDEMI*DZE 84 : SHP(1,1)=-AUX*DZEM*PAUX2 85 : SHP(1,2)=DEUX*QSI*AUX*DZEM 86 : SHP(1,3)=QSI*DZEM*PAUX 87 : SHP(2,1)=(PAUX2-AUX)*DZEM 88 : SHP(2,2)=DEUX*(AUX-QSI)*DZEM 89 : SHP(2,3)=(PAUX+QSI)*DZEM 90 : SHP(3,1)=(PAUX2-AUX)*DZEM 91 : SHP(3,2)=-DEUX*QSI*DZEM 92 : SHP(3,3)=XZER 93 : SHP(4,1)=AUX*(PAUX2-UNDEMI*DZEM) 94 : SHP(4,2)=-DEUX*QSI*AUX 95 : SHP(4,3)=-QSI*(UNDEMI*DZEM+PAUX) 96 : SHP(1,4)=DEUX*QSI*ETA*DZEM 97 : SHP(1,5)=ETA*DZEM*PAUX1 98 : SHP(1,6)=DEUX*ETA*AUX*DZEM 99 : SHP(2,4)=DEUX*ETA*DZEM 100 : SHP(2,5)=XZER 101 : SHP(2,6)=-DEUX*ETA*DZEM 102 : SHP(3,4)=DEUX*QSI*DZEM 103 : SHP(3,5)=(ETA+PAUX1)*DZEM 104 : SHP(3,6)=DEUX*(AUX-ETA)*DZEM 105 : SHP(4,4)=-DEUX*QSI*ETA 106 : SHP(4,5)=-ETA*(PAUX1+UNDEMI*DZEM) 107 : SHP(4,6)=-DEUX*ETA*AUX 108 : SHP(1,7)=AUX*DZEM*DZEP 109 : SHP(1,8)=QSI*DZEM*DZEP 110 : SHP(1,9)=ETA*DZEM*DZEP 111 : SHP(2,7)=-DZEM*DZEP 112 : SHP(2,8)=DZEM*DZEP 113 : SHP(2,9)=XZER 114 : SHP(3,7)=-DZEM*DZEP 115 : SHP(3,8)=XZER 116 : SHP(3,9)=DZEM*DZEP 117 : SHP(4,7)=-DEUX*DZE*AUX 118 : SHP(4,8)=-DEUX*DZE*QSI 119 : SHP(4,9)=-DEUX*DZE*ETA 120 : SHP(1,10)=-AUX*DZEP*PAUX5 121 : SHP(1,11)=DEUX*QSI*AUX*DZEP 122 : SHP(1,12)=QSI*DZEP*PAUX3 123 : SHP(2,10)=(PAUX5-AUX)*DZEP 124 : SHP(2,11)=DEUX*(AUX-QSI)*DZEP 125 : SHP(2,12)=(PAUX3+QSI)*DZEP 126 : SHP(3,10)=(PAUX5-AUX)*DZEP 127 : SHP(3,11)=-DEUX*QSI*DZEP 128 : SHP(3,12)=XZER 129 : SHP(4,10)=AUX*(UNDEMI*DZEP-PAUX5) 130 : SHP(4,11)=DEUX*QSI*AUX 131 : SHP(4,12)=QSI*(PAUX3+UNDEMI*DZEP) 132 : SHP(1,13)=DEUX*QSI*ETA*DZEP 133 : SHP(1,14)=ETA*DZEP*PAUX4 134 : SHP(1,15)=DEUX*ETA*AUX*DZEP 135 : SHP(2,13)=DEUX*ETA*DZEP 136 : SHP(2,14)=XZER 137 : SHP(2,15)=-DEUX*ETA*DZEP 138 : SHP(3,13)=DEUX*QSI*DZEP 139 : SHP(3,14)=(ETA+PAUX4)*DZEP 140 : SHP(3,15)=DEUX*(AUX-ETA)*DZEP 141 : SHP(4,13)=DEUX*QSI*ETA 142 : SHP(4,14)=ETA*(PAUX4+UNDEMI*DZEP) 143 : SHP(4,15)=DEUX*ETA*AUX 144 : GOTO 666 145 : 23 CONTINUE 146 : C 147 : C FONCTIONS DE FORME TETRAEDRE A 4 NOEUDS 148 : C 149 : AUX=UN-QSI-ETA-DZE 150 : SHP(1,1)=AUX 151 : SHP(1,2)=QSI 152 : SHP(1,3)=ETA 153 : SHP(1,4)=DZE 154 : SHP(2,1)=-UN 155 : SHP(2,2)= UN 156 : SHP(2,3)= XZER 157 : SHP(2,4)= XZER 158 : SHP(3,1)=-UN 159 : SHP(3,2)= XZER 160 : SHP(3,3)= UN 161 : SHP(3,4)= XZER 162 : SHP(4,1)=-UN 163 : SHP(4,2)= XZER 164 : SHP(4,3)= XZER 165 : SHP(4,4)= UN 166 : GOTO 666 167 : 24 CONTINUE 168 : C 169 : C FONCTIONS DE FORME TETRAEDRE A 10 NOEUDS 170 : C 171 : AUX=UN-QSI-ETA-DZE 172 : SHP(1,1)=AUX*(DEUX*AUX-UN) 173 : SHP(1,2)=QUATRE*QSI*AUX 174 : SHP(1,3)=QSI*(DEUX*QSI-UN) 175 : SHP(1,4)=QUATRE*QSI*ETA 176 : SHP(1,5)=ETA*(DEUX*ETA-UN) 177 : SHP(1,6)=QUATRE*ETA*AUX 178 : SHP(1,7)=QUATRE*DZE*AUX 179 : SHP(1,8)=QUATRE*QSI*DZE 180 : SHP(1,9)=QUATRE*ETA*DZE 181 : SHP(1,10)=DZE*(DEUX*DZE-UN) 182 : SHP(2,1)=UN-QUATRE*AUX 183 : SHP(2,2)=QUATRE*(AUX-QSI) 184 : SHP(2,3)=QUATRE*QSI-UN 185 : SHP(2,4)=QUATRE*ETA 186 : SHP(2,5)=XZER 187 : SHP(2,6)=-QUATRE*ETA 188 : SHP(2,7)=-QUATRE*DZE 189 : SHP(2,8)= QUATRE*DZE 190 : SHP(2,9)= XZER 191 : SHP(2,10)=XZER 192 : SHP(3,1)= UN-QUATRE*AUX 193 : SHP(3,2)= -QUATRE*QSI 194 : SHP(3,3)= XZER 195 : SHP(3,4)= QUATRE*QSI 196 : SHP(3,5)= QUATRE*ETA-UN 197 : SHP(3,6)= QUATRE*(AUX-ETA) 198 : SHP(3,7)=-QUATRE*DZE 199 : SHP(3,8)= XZER 200 : SHP(3,9)= QUATRE*DZE 201 : SHP(3,10)=XZER 202 : SHP(4,1)= UN-QUATRE*AUX 203 : SHP(4,2)= -QUATRE*QSI 204 : SHP(4,3)= XZER 205 : SHP(4,4)= XZER 206 : SHP(4,5)= XZER 207 : SHP(4,6)= -QUATRE*ETA 208 : SHP(4,7)=QUATRE*(AUX-DZE) 209 : SHP(4,8)=QUATRE*QSI 210 : SHP(4,9)=QUATRE*ETA 211 : SHP(4,10)=QUATRE*DZE-UN 212 : GOTO 666 213 : 14 CONTINUE 214 : C 215 : C FONCTIONS DE FORME CUBE A 8 NOEUDS 216 : C 217 : QSIM= UN-QSI 218 : QSIP= UN+QSI 219 : ETAM= UN-ETA 220 : ETAP= UN+ETA 221 : DZEM= UN-DZE 222 : DZEP= UN+DZE 223 : SHP(1,1)=QSIM*ETAM*DZEM/HUIT 224 : SHP(1,2)=QSIP*ETAM*DZEM/HUIT 225 : SHP(1,3)=QSIP*ETAP*DZEM/HUIT 226 : SHP(1,4)=QSIM*ETAP*DZEM/HUIT 227 : SHP(1,5)=QSIM*ETAM*DZEP/HUIT 228 : SHP(1,6)=QSIP*ETAM*DZEP/HUIT 229 : SHP(1,7)=QSIP*ETAP*DZEP/HUIT 230 : SHP(1,8)=QSIM*ETAP*DZEP/HUIT 231 : SHP(2,1)= -ETAM*DZEM/HUIT 232 : SHP(2,2)= ETAM*DZEM/HUIT 233 : SHP(2,3)= ETAP*DZEM/HUIT 234 : SHP(2,4)= -ETAP*DZEM/HUIT 235 : SHP(2,5)= -ETAM*DZEP/HUIT 236 : SHP(2,6)= ETAM*DZEP/HUIT 237 : SHP(2,7)= ETAP*DZEP/HUIT 238 : SHP(2,8)= -ETAP*DZEP/HUIT 239 : SHP(3,1)=-QSIM*DZEM/HUIT 240 : SHP(3,2)=-QSIP*DZEM/HUIT 241 : SHP(3,3)= QSIP*DZEM/HUIT 242 : SHP(3,4)= QSIM*DZEM/HUIT 243 : SHP(3,5)=-QSIM*DZEP/HUIT 244 : SHP(3,6)=-QSIP*DZEP/HUIT 245 : SHP(3,7)= QSIP*DZEP/HUIT 246 : SHP(3,8)= QSIM*DZEP/HUIT 247 : SHP(4,1)=-QSIM*ETAM/HUIT 248 : SHP(4,2)=-QSIP*ETAM/HUIT 249 : SHP(4,3)=-QSIP*ETAP/HUIT 250 : SHP(4,4)=-QSIM*ETAP/HUIT 251 : SHP(4,5)= QSIM*ETAM/HUIT 252 : SHP(4,6)= QSIP*ETAM/HUIT 253 : SHP(4,7)= QSIP*ETAP/HUIT 254 : SHP(4,8)= QSIM*ETAP/HUIT 255 : GOTO 666 256 : 15 CONTINUE 257 : C 258 : C FONCTIONS DE FORME CUBE A 20 NOEUDS 259 : C 260 : QSIM=(UN-QSI) 261 : QSIP=(UN+QSI) 262 : ETAM=(UN-ETA) 263 : ETAP=(UN+ETA) 264 : DZEM=(UN-DZE) 265 : DZEP=(UN+DZE) 266 : C NOEUDS AUX SOMMETS 1 3 5 7 13 15 17 19 267 : SHP(1,1 )=QSIM*ETAM*DZEM*(QSIM+ETAM+DZEM-CINQ)/HUIT 268 : SHP(1,3 )=QSIP*ETAM*DZEM*(QSIP+ETAM+DZEM-CINQ)/HUIT 269 : SHP(1,5 )=QSIP*ETAP*DZEM*(QSIP+ETAP+DZEM-CINQ)/HUIT 270 : SHP(1,7 )=QSIM*ETAP*DZEM*(QSIM+ETAP+DZEM-CINQ)/HUIT 271 : SHP(1,13)=QSIM*ETAM*DZEP*(QSIM+ETAM+DZEP-CINQ)/HUIT 272 : SHP(1,15)=QSIP*ETAM*DZEP*(QSIP+ETAM+DZEP-CINQ)/HUIT 273 : SHP(1,17)=QSIP*ETAP*DZEP*(QSIP+ETAP+DZEP-CINQ)/HUIT 274 : SHP(1,19)=QSIM*ETAP*DZEP*(QSIM+ETAP+DZEP-CINQ)/HUIT 275 : SHP(2,1 )= -ETAM*DZEM*(DEUX*QSIM+ETAM+DZEM-CINQ)/HUIT 276 : SHP(2,3 )= ETAM*DZEM*(DEUX*QSIP+ETAM+DZEM-CINQ)/HUIT 277 : SHP(2,5 )= ETAP*DZEM*(DEUX*QSIP+ETAP+DZEM-CINQ)/HUIT 278 : SHP(2,7 )= -ETAP*DZEM*(DEUX*QSIM+ETAP+DZEM-CINQ)/HUIT 279 : SHP(2,13)= -ETAM*DZEP*(DEUX*QSIM+ETAM+DZEP-CINQ)/HUIT 280 : SHP(2,15)= ETAM*DZEP*(DEUX*QSIP+ETAM+DZEP-CINQ)/HUIT 281 : SHP(2,17)= ETAP*DZEP*(DEUX*QSIP+ETAP+DZEP-CINQ)/HUIT 282 : SHP(2,19)= -ETAP*DZEP*(DEUX*QSIM+ETAP+DZEP-CINQ)/HUIT 283 : SHP(3,1 )= -QSIM*DZEM*(QSIM+DEUX*ETAM+DZEM-CINQ)/HUIT 284 : SHP(3,3 )= -QSIP*DZEM*(QSIP+DEUX*ETAM+DZEM-CINQ)/HUIT 285 : SHP(3,5 )= QSIP*DZEM*(QSIP+DEUX*ETAP+DZEM-CINQ)/HUIT 286 : SHP(3,7 )= QSIM*DZEM*(QSIM+DEUX*ETAP+DZEM-CINQ)/HUIT 287 : SHP(3,13)= -QSIM*DZEP*(QSIM+DEUX*ETAM+DZEP-CINQ)/HUIT 288 : SHP(3,15)= -QSIP*DZEP*(QSIP+DEUX*ETAM+DZEP-CINQ)/HUIT 289 : SHP(3,17)= QSIP*DZEP*(QSIP+DEUX*ETAP+DZEP-CINQ)/HUIT 290 : SHP(3,19)= QSIM*DZEP*(QSIM+DEUX*ETAP+DZEP-CINQ)/HUIT 291 : SHP(4,1 )= -QSIM*ETAM*(QSIM+ETAM+DEUX*DZEM-CINQ)/HUIT 292 : SHP(4,3 )= -QSIP*ETAM*(QSIP+ETAM+DEUX*DZEM-CINQ)/HUIT 293 : SHP(4,5 )= -QSIP*ETAP*(QSIP+ETAP+DEUX*DZEM-CINQ)/HUIT 294 : SHP(4,7 )= -QSIM*ETAP*(QSIM+ETAP+DEUX*DZEM-CINQ)/HUIT 295 : SHP(4,13)= QSIM*ETAM*(QSIM+ETAM+DEUX*DZEP-CINQ)/HUIT 296 : SHP(4,15)= QSIP*ETAM*(QSIP+ETAM+DEUX*DZEP-CINQ)/HUIT 297 : SHP(4,17)= QSIP*ETAP*(QSIP+ETAP+DEUX*DZEP-CINQ)/HUIT 298 : SHP(4,19)= QSIM*ETAP*(QSIM+ETAP+DEUX*DZEP-CINQ)/HUIT 299 : C NOEUDS SUR LES COTES PARALLELES A L AXE QSI 2 6 14 18 300 : SHP(1,2 )=UNQUA*(UN-QSI*QSI)*ETAM*DZEM 301 : SHP(1,6 )=UNQUA*(UN-QSI*QSI)*ETAP*DZEM 302 : SHP(1,14)=UNQUA*(UN-QSI*QSI)*ETAM*DZEP 303 : SHP(1,18)=UNQUA*(UN-QSI*QSI)*ETAP*DZEP 304 : SHP(2,2 )=-UNDEMI*QSI*ETAM*DZEM 305 : SHP(2,6 )=-UNDEMI*QSI*ETAP*DZEM 306 : SHP(2,14)=-UNDEMI*QSI*ETAM*DZEP 307 : SHP(2,18)=-UNDEMI*QSI*ETAP*DZEP 308 : SHP(3,2 )=-UNQUA*(UN-QSI*QSI)*DZEM 309 : SHP(3,6 )= UNQUA*(UN-QSI*QSI)*DZEM 310 : SHP(3,14)=-UNQUA*(UN-QSI*QSI)*DZEP 311 : SHP(3,18)= UNQUA*(UN-QSI*QSI)*DZEP 312 : SHP(4,2 )=-UNQUA*(UN-QSI*QSI)*ETAM 313 : SHP(4,6 )=-UNQUA*(UN-QSI*QSI)*ETAP 314 : SHP(4,14)= UNQUA*(UN-QSI*QSI)*ETAM 315 : SHP(4,18)= UNQUA*(UN-QSI*QSI)*ETAP 316 : C NOEUDS SUR LES COTES PARALELLES A L AXE ETA 4 8 16 20 317 : SHP(1,4 )=UNQUA*QSIP*(UN-ETA*ETA)*DZEM 318 : SHP(1,8 )=UNQUA*QSIM*(UN-ETA*ETA)*DZEM 319 : SHP(1,16)=UNQUA*QSIP*(UN-ETA*ETA)*DZEP 320 : SHP(1,20)=UNQUA*QSIM*(UN-ETA*ETA)*DZEP 321 : SHP(2,4 )= UNQUA*(UN-ETA*ETA)*DZEM 322 : SHP(2,8 )=-UNQUA*(UN-ETA*ETA)*DZEM 323 : SHP(2,16)= UNQUA*(UN-ETA*ETA)*DZEP 324 : SHP(2,20)=-UNQUA*(UN-ETA*ETA)*DZEP 325 : SHP(3,4 )=-UNDEMI*QSIP*ETA*DZEM 326 : SHP(3,8 )=-UNDEMI*QSIM*ETA*DZEM 327 : SHP(3,16)=-UNDEMI*QSIP*ETA*DZEP 328 : SHP(3,20)=-UNDEMI*QSIM*ETA*DZEP 329 : SHP(4,4 )=-UNQUA*QSIP*(UN-ETA*ETA) 330 : SHP(4,8 )=-UNQUA*QSIM*(UN-ETA*ETA) 331 : SHP(4,16)= UNQUA*QSIP*(UN-ETA*ETA) 332 : SHP(4,20)= UNQUA*QSIM*(UN-ETA*ETA) 333 : C NOEUDS SUR LES COTES PARALELLES A L AXE DZE 9 10 11 12 334 : SHP(1,9 )= UNQUA*QSIM*ETAM*(UN-DZE*DZE) 335 : SHP(1,10)= UNQUA*QSIP*ETAM*(UN-DZE*DZE) 336 : SHP(1,11)= UNQUA*QSIP*ETAP*(UN-DZE*DZE) 337 : SHP(1,12)= UNQUA*QSIM*ETAP*(UN-DZE*DZE) 338 : SHP(2,9 )=-UNQUA*ETAM*(UN-DZE*DZE) 339 : SHP(2,10)= UNQUA*ETAM*(UN-DZE*DZE) 340 : SHP(2,11)= UNQUA*ETAP*(UN-DZE*DZE) 341 : SHP(2,12)=-UNQUA*ETAP*(UN-DZE*DZE) 342 : SHP(3,9 )=-UNQUA*QSIM*(UN-DZE*DZE) 343 : SHP(3,10)=-UNQUA*QSIP*(UN-DZE*DZE) 344 : SHP(3,11)= UNQUA*QSIP*(UN-DZE*DZE) 345 : SHP(3,12)= UNQUA*QSIM*(UN-DZE*DZE) 346 : SHP(4,9 )=-UNDEMI*QSIM*ETAM*DZE 347 : SHP(4,10)=-UNDEMI*QSIP*ETAM*DZE 348 : SHP(4,11)=-UNDEMI*QSIP*ETAP*DZE 349 : SHP(4,12)=-UNDEMI*QSIM*ETAP*DZE 350 : GOTO 666 351 : 25 CONTINUE 352 : C 353 : C FONCTIONS DE FORME PYRAMIDE A 5 NOEUDS 354 : C 355 : DZEM =UN-DZE 356 : DZEM4=(UN-DZE)*QUATRE 357 : AUX = QSI+ETA+DZE-UN 358 : AUX1=-QSI+ETA+DZE-UN 359 : AUX2=-QSI-ETA+DZE-UN 360 : AUX3= QSI-ETA+DZE-UN 361 : IF(ABS(DZE-UN).LT.1.D-6) THEN 362 : DO 251 IB=1,5 363 : DO 251 IA=1,3 364 : SHP(IA,IB)= XZER 365 : 251 CONTINUE 366 : SHP(1,5)= UN 367 : SHP(2,1)= UNDEMI 368 : SHP(2,3)=-UNDEMI 369 : SHP(3,2)= UNDEMI 370 : SHP(3,4)=-UNDEMI 371 : DO 252 IA=1,4 372 : SHP(4,IA)=-UNDEMI 373 : 252 CONTINUE 374 : SHP(4,5)= UN 375 : GO TO 666 376 : ENDIF 377 : SHP(1,1)=AUX1*AUX2/DZEM4 378 : SHP(1,2)=AUX2*AUX3/DZEM4 379 : SHP(1,3)=AUX *AUX3/DZEM4 380 : SHP(1,4)=AUX *AUX1/DZEM4 381 : SHP(1,5)=DZE 382 : SHP(2,1)=(-AUX1-AUX2)/DZEM4 383 : SHP(2,2)=( AUX2-AUX3)/DZEM4 384 : SHP(2,3)=( AUX +AUX3)/DZEM4 385 : SHP(2,4)=( AUX1-AUX )/DZEM4 386 : SHP(2,5)=XZER 387 : SHP(3,1)=( AUX2-AUX1)/DZEM4 388 : SHP(3,2)=(-AUX2-AUX3)/DZEM4 389 : SHP(3,3)=( AUX3-AUX )/DZEM4 390 : SHP(3,4)=( AUX +AUX1)/DZEM4 391 : SHP(3,5)=XZER 392 : SHP(4,1)=(AUX1+AUX2+AUX1*AUX2/DZEM)/DZEM4 393 : SHP(4,2)=(AUX2+AUX3+AUX2*AUX3/DZEM)/DZEM4 394 : SHP(4,3)=(AUX3+AUX +AUX3*AUX /DZEM)/DZEM4 395 : SHP(4,4)=( AUX+AUX1+ AUX*AUX1/DZEM)/DZEM4 396 : SHP(4,5)=UN 397 : GOTO 666 398 : 30 CONTINUE 399 : C 400 : C FONCTIONS DE FORME PYRAMIDE A 13 NOEUDS 401 : C 402 : DZEM=UN-DZE 403 : DZEM2=(UN-DZE)*DEUX 404 : PAUX = QSI+ETA+DZE-UN 405 : PAUX1=-QSI+ETA+DZE-UN 406 : PAUX2=-QSI-ETA+DZE-UN 407 : PAUX3= QSI-ETA+DZE-UN 408 : QAUX = QSI-UNDEMI 409 : QAUX1= ETA-UNDEMI 410 : QAUX2=-QSI-UNDEMI 411 : QAUX3=-ETA-UNDEMI 412 : IF(ABS(DZE-UN).LT.1.D-6) THEN 413 : DO 301 IB=1,13 414 : DO 301 IA=1,3 415 : SHP(IA,IB)= XZER 416 : 301 CONTINUE 417 : SHP(1,13)= UN 418 : SHP(2,1 )=-UNDEMI 419 : SHP(2,5 )= UNDEMI 420 : SHP(2,9 )= DEUX 421 : SHP(2,11)=-DEUX 422 : SHP(3,3 )=-UNDEMI 423 : SHP(3,7 )= UNDEMI 424 : SHP(3,10)= DEUX 425 : SHP(3,12)=-DEUX 426 : DO 302 IA=1,7,2 427 : SHP(4,IA)= UN/QUATRE 428 : SHP(4,IA+1)= XZER 429 : 302 CONTINUE 430 : DO 303 IA=9,12 431 : SHP(4,IA)=-UN 432 : 303 CONTINUE 433 : SHP(4,13)= TROIS 434 : GO TO 666 435 : ENDIF 436 : SHP(1,1 )= PAUX1*PAUX2*QAUX /DZEM2 437 : SHP(1,2 )=-PAUX1*PAUX2*PAUX3/DZEM2 438 : SHP(1,3 )= PAUX2*PAUX3*QAUX1/DZEM2 439 : SHP(1,4 )=-PAUX2*PAUX3*PAUX /DZEM2 440 : SHP(1,5 )= PAUX3*PAUX *QAUX2/DZEM2 441 : SHP(1,6 )=-PAUX3*PAUX *PAUX1/DZEM2 442 : SHP(1,7 )= PAUX *PAUX1*QAUX3/DZEM2 443 : SHP(1,8 )=-PAUX *PAUX1*PAUX2/DZEM2 444 : SHP(1,9 )= DZE*PAUX1*PAUX2/DZEM 445 : SHP(1,10)= DZE*PAUX2*PAUX3/DZEM 446 : SHP(1,11)= DZE*PAUX3*PAUX /DZEM 447 : SHP(1,12)= DZE*PAUX *PAUX1/DZEM 448 : SHP(1,13)= DEUX*DZE*(DZE-UNDEMI) 449 : C 450 : SHP(2,1 )=(PAUX1*PAUX2-(PAUX1+PAUX2)*QAUX)/DZEM2 451 : SHP(2,2 )=(PAUX2*PAUX3+PAUX1*PAUX3-PAUX1*PAUX2)/DZEM2 452 : SHP(2,3 )=(PAUX2-PAUX3)*QAUX1/DZEM2 453 : SHP(2,4 )=(PAUX3*PAUX-PAUX2*PAUX-PAUX2*PAUX3)/DZEM2 454 : SHP(2,5 )=((PAUX+PAUX3)*QAUX2-PAUX3*PAUX)/DZEM2 455 : SHP(2,6 )=(PAUX3*PAUX-PAUX*PAUX1-PAUX3*PAUX1)/DZEM2 456 : SHP(2,7 )=(PAUX1-PAUX)*QAUX3/DZEM2 457 : SHP(2,8 )=(PAUX*PAUX2+PAUX*PAUX1-PAUX1*PAUX2)/DZEM2 458 : SHP(2,9 )=(-PAUX2-PAUX1)*DZE/DZEM 459 : SHP(2,10)=(PAUX2-PAUX3)*DZE/DZEM 460 : SHP(2,11)=(PAUX+PAUX3)*DZE/DZEM 461 : SHP(2,12)=(PAUX1-PAUX)*DZE/DZEM 462 : SHP(2,13)=XZER 463 : C 464 : SHP(3,1 )=(PAUX2-PAUX1)*QAUX/DZEM2 465 : SHP(3,2 )=(PAUX1*PAUX3+PAUX1*PAUX2-PAUX2*PAUX3)/DZEM2 466 : SHP(3,3 )=(PAUX2*PAUX3-(PAUX2+PAUX3)*QAUX1)/DZEM2 467 : SHP(3,4 )=(PAUX3*PAUX+PAUX2*PAUX-PAUX2*PAUX3)/DZEM2 468 : SHP(3,5 )=(PAUX3-PAUX)*QAUX2/DZEM2 469 : SHP(3,6 )=(PAUX*PAUX1-PAUX3*PAUX1-PAUX3*PAUX)/DZEM2 470 : SHP(3,7 )=((PAUX1+PAUX)*QAUX3-PAUX*PAUX1)/DZEM2 471 : SHP(3,8 )=(PAUX*PAUX1-PAUX1*PAUX2-PAUX*PAUX2)/DZEM2 472 : SHP(3,9 )=(PAUX2-PAUX1)*DZE/DZEM 473 : SHP(3,10)=(-PAUX3-PAUX2)*DZE/DZEM 474 : SHP(3,11)=(PAUX3-PAUX)*DZE/DZEM 475 : SHP(3,12)=(PAUX+PAUX1)*DZE/DZEM 476 : SHP(3,13)=XZER 477 : C 478 : SHP(4,1 )=(PAUX1+PAUX2+PAUX1*PAUX2/DZEM)*QAUX/DZEM2 479 : SHP(4,2 )=-UN*(PAUX2*PAUX3+PAUX1*PAUX3+PAUX1*PAUX2+ 480 : + PAUX1*PAUX2*PAUX3/DZEM)/DZEM2 481 : SHP(4,3 )=(PAUX2+PAUX3+PAUX2*PAUX3/DZEM)*QAUX1/DZEM2 482 : SHP(4,4 )=-UN*(PAUX3*PAUX +PAUX2*PAUX+PAUX2*PAUX3+ 483 : + PAUX2*PAUX3*PAUX/DZEM)/DZEM2 484 : SHP(4,5 )=(PAUX+PAUX3+PAUX*PAUX3/DZEM)*QAUX2/DZEM2 485 : SHP(4,6 )=-UN*(PAUX*PAUX1+PAUX3*PAUX1+PAUX3*PAUX+ 486 : + PAUX3*PAUX*PAUX1/DZEM)/DZEM2 487 : SHP(4,7 )=(PAUX1+PAUX+PAUX*PAUX1/DZEM)*QAUX3/DZEM2 488 : SHP(4,8 )=-UN*(PAUX1*PAUX2+PAUX*PAUX2+PAUX*PAUX1+ 489 : + PAUX*PAUX1*PAUX2/DZEM)/DZEM2 490 : SHP(4,9 )=PAUX1*PAUX2/DZEM/DZEM+(PAUX2+PAUX1)*DZE/DZEM 491 : SHP(4,10)=PAUX2*PAUX3/DZEM/DZEM+(PAUX3+PAUX2)*DZE/DZEM 492 : SHP(4,11)=PAUX3*PAUX /DZEM/DZEM+(PAUX +PAUX3)*DZE/DZEM 493 : SHP(4,12)=PAUX *PAUX1/DZEM/DZEM+(PAUX +PAUX1)*DZE/DZEM 494 : SHP(4,13)=QUATRE*DZE-UN 495 : GO TO 666 496 : C 497 : 666 CONTINUE 498 : RETURN 499 : END 500 :
© Cast3M 2003 - All rights reserved.
Disclaimer