1 : C IDVAR5 SOURCE BR232186 12/08/01 22:28:53 7456 2 : 3 : *_______________________________________________________________________ 4 : * RECHERCHE DES NOMS DE COMPOSANTES DES VARIABLES INTERNES 5 : * -------------------------------------------------------- 6 : * CAS DE L'ENDOMMAGEMENT 7 : * --------------------- 8 : * 9 : * ENTREES : 10 : * --------- 11 : * MATEPL NUMERO DU MATERIAU ENDOMAGEABLE 12 : * MFR NUMERO DE FORMULATION 13 : * IFOUR VALEUR DE IFOU DANS CCOPTIO 14 : * IFLAG 0=TOUTES LES VARIABLES 15 : * 1=FONCTIONS SEUIL 16 : * 17 : * SORTIES : 18 : * --------- 19 : * IPCOMP POINTEUR SUR LA LISTE DES NOMS DE COMPOSANTES 20 : * OBLIGATOIRES ET FACULTATIVES 21 : * NBROBL NOMBRE DE COMPOSANTES OBLIGATOIRES 22 : * NBRFAC NOMBRE DE COMPOSANTES FACTULTATIVES 23 : * 24 : * PP 16/8/91 ajout var.int pour l'endommagement 25 : *_______________________________________________________________________ 26 : 27 : SUBROUTINE IDVAR5 (MATEPL,MFR,IFOUR,IFLAG,IPCOMP,NBROBL,NBRFAC) 28 : 29 : IMPLICIT INTEGER(I-N) 30 : IMPLICIT REAL*8 (A-H,O-Z) 31 : 32 : SEGMENT NOMID 33 : CHARACTER*8 LESOBL(NBROBL),LESFAC(NBRFAC) 34 : ENDSEGMENT 35 : 36 : NBROBL=0 37 : NBRFAC=0 38 : NOMID=0 39 : C PRINt *,'ENTREE IDVAR5 MATEPL, MFR =',MATEPL,MFR 40 : 41 : C 1 - Cas du modele SCALAIRE 42 : C ============================ 43 : IF (MATEPL.EQ.30.OR.MATEPL.EQ.118) THEN 44 : IF (IFLAG.EQ.0) THEN 45 : NBROBL=2 46 : SEGINI,NOMID 47 : LESOBL(1)='EPTI' 48 : LESOBL(2)='D' 49 : ELSE 50 : NBROBL=1 51 : SEGINI,NOMID 52 : LESOBL(1)='EPTI' 53 : ENDIF 54 : 55 : C 2 - Cas du modele UNILATERAL simple ou avec viscosite 56 : C ======================================================= 57 : ELSE IF (MATEPL.EQ.31.OR.MATEPL.EQ. 97) THEN 58 : IF (IFLAG.EQ.0) THEN 59 : NBROBL=4 60 : SEGINI,NOMID 61 : LESOBL(1)='YLI1' 62 : LESOBL(2)='YLI2' 63 : LESOBL(3)='DOM1' 64 : LESOBL(4)='DOM2' 65 : ELSE 66 : NBROBL=2 67 : SEGINI,NOMID 68 : LESOBL(1)='YLI1' 69 : LESOBL(2)='YLI2' 70 : ENDIF 71 : 72 : C 3 - Cas du modele ROTATING_CRACK 73 : C ================================== 74 : ELSE IF (MATEPL.EQ.37) THEN 75 : NBROBL=2 76 : SEGINI,NOMID 77 : LESOBL(1)='EPSM' 78 : LESOBL(2)='EPSN' 79 : 80 : C 4 - Cas du modele SIC_SIC 81 : C Modele disponible seulement en MASSIF 3D 82 : C ============================================== 83 : ELSE IF (MATEPL.EQ.88) THEN 84 : IF (MFR.EQ.1) THEN 85 : NBROBL=10 86 : SEGINI,NOMID 87 : LESOBL( 1)='EPSE' 88 : LESOBL( 2)='DOM1' 89 : LESOBL( 3)='DOM2' 90 : LESOBL( 4)='DOM3' 91 : LESOBL( 5)='DF11' 92 : LESOBL( 6)='DF22' 93 : LESOBL( 7)='DF33' 94 : LESOBL( 8)='DF23' 95 : LESOBL( 9)='DF13' 96 : LESOBL(10)='DF12' 97 : ELSE 97 : CALL ERREUR (19) 98 : ENDIF 99 : 100 : C 5 - Cas du modele MICROPLANS ORTHOTROPE 101 : C ========================================= 102 : ELSE IF (MATEPL.EQ.96) THEN 103 : NBROBL=14 104 : SEGINI,NOMID 105 : LESOBL( 1)='EPSE' 106 : LESOBL( 2)='R ' 107 : C= Tenseur d'endommagement 108 : LESOBL( 3)='D11 ' 109 : LESOBL( 4)='D22 ' 110 : LESOBL( 5)='D33 ' 111 : LESOBL( 6)='D12 ' 112 : LESOBL( 7)='D13 ' 113 : LESOBL( 8)='D23 ' 114 : C= Tenseur des contraintes anelastiques dues a l'endommagement 115 : LESOBL( 9)='SA11' 116 : LESOBL(10)='SA22' 117 : LESOBL(11)='SA33' 118 : LESOBL(12)='SA12' 119 : LESOBL(13)='SA13' 120 : LESOBL(14)='SA23' 121 : 122 : C 5 - Cas du modele MICROPLANS ISOTROPE 123 : C ======================================= 124 : ELSE IF (MATEPL.EQ.98) THEN 125 : NBROBL=9 126 : SEGINI,NOMID 127 : LESOBL(1)='EPSE' 128 : LESOBL(2)='R ' 129 : C= Scalaire d'endommagement 130 : LESOBL(3)='DOM' 131 : C= Tenseur des contraintes anelastiques dues a l'endommagement 132 : LESOBL(4)='SA11' 133 : LESOBL(5)='SA22' 134 : LESOBL(6)='SA33' 135 : LESOBL(7)='SA12' 136 : LESOBL(8)='SA13' 137 : LESOBL(9)='SA23' 138 : 139 : * RDj 140 : * 141 : ELSEIF(MATEPL.EQ.124)THEN 142 : * 143 : * CAS DU MODELE SCALAIRE AVEC DEFORMATION ANELASTIQUE 144 : * 145 : NBROBL=24 146 : SEGINI,NOMID 147 : LESOBL(1) ='ECM' 148 : LESOBL(2) ='D' 149 : LESOBL(3) ='DELT' 150 : LESOBL(4) ='EIC1' 151 : LESOBL(5) ='EIC2' 152 : LESOBL(6) ='EIC3' 153 : LESOBL(7) ='EIC4' 154 : LESOBL(8) ='EIC5' 155 : LESOBL(9) ='EIC6' 156 : LESOBL(10) ='EIT1' 157 : LESOBL(11) ='EIT2' 158 : LESOBL(12) ='EIT3' 159 : LESOBL(13) ='EIT4' 160 : LESOBL(14) ='EIT5' 161 : LESOBL(15) ='EIT6' 162 : LESOBL(16) ='EPTI' 163 : LESOBL(17) ='VER1' 164 : LESOBL(18) ='VER2' 165 : LESOBL(19) ='VER3' 166 : LESOBL(20) ='VER4' 167 : LESOBL(21) ='VER5' 168 : LESOBL(22) ='VER6' 169 : LESOBL(23) ='ENDP' 170 : LESOBL(24) ='ENDM' 171 : * RDj 172 : * 173 : * CCd 174 : * 175 : ELSEIF(MATEPL.EQ.125)THEN 176 : * 177 : * MODEL DAMAGE_TC - R. Faria & J. Oliver 178 : * 179 : C NBROBL=12 180 : NBROBL=6 181 : SEGINI NOMID 182 : LESOBL(1)='THRT' 183 : LESOBL(2)='EQUT' 184 : LESOBL(3)='DAMT' 185 : LESOBL(4)='THRC' 186 : LESOBL(5)='EQUC' 187 : LESOBL(6)='DAMC' 188 : * 189 : C LESOBL(7) ='SE11' 190 : C LESOBL(8) ='SE22' 191 : C LESOBL(9) ='SE33' 192 : C LESOBL(10)='SE12' 193 : C LESOBL(11)='SE13' 194 : C LESOBL(12)='SE23' 195 : * 196 : 197 : C 11 - Cas du modele SICf_SIC Scalaire ONERA 198 : C SICSCAL 199 : C Modele disponible seulement en MASSIF 3D 200 : C ============================================== 201 : ELSE IF (MATEPL.EQ.131) THEN 202 : IF (MFR.EQ.1) THEN 203 : NBROBL=19 204 : SEGINI,NOMID 205 : * Somme des defo inel 206 : LESOBL( 1)='EPSE' 207 : * Dommage dans les différentes directions 208 : LESOBL( 2)='DOM1' 209 : LESOBL( 3)='DOM2' 210 : LESOBL( 4)='DOM3' 211 : * Déformation repère orthotrope ordre ONERA 212 : LESOBL( 5)='DF11' 213 : LESOBL( 6)='DF22' 214 : LESOBL( 7)='DF33' 215 : LESOBL( 8)='DF23' 216 : LESOBL( 9)='DF13' 217 : LESOBL(10)='DF12' 218 : * Indice d'activation des fissures 219 : LESOBL(11)='NUA1' 220 : LESOBL(12)='NUA2' 221 : LESOBL(13)='NUA3' 222 : * Déformations inélastiques repère orthotrope ordre ONERA 223 : LESOBL(14)='EPI1' 224 : LESOBL(15)='EPI2' 225 : LESOBL(16)='EPI3' 226 : LESOBL(17)='EPI4' 227 : LESOBL(18)='EPI5' 228 : LESOBL(19)='EPI6' 229 : ELSE 229 : CALL ERREUR (19) 230 : ENDIF 231 : C 12 - Cas du modele SICf_SIC pseudo Tensoriel ONERA 232 : C SICTENS 233 : C Modele disponible seulement en MASSIF 3D 234 : C ============================================== 235 : ELSE IF (MATEPL.EQ.132) THEN 236 : IF (MFR.EQ.1) THEN 237 : NBROBL=23 238 : SEGINI,NOMID 239 : * somme defo inel 240 : LESOBL( 1)='EPSE' 241 : * 5 variables de dommage 242 : LESOBL( 2)='DOM1' 243 : LESOBL( 3)='DOM2' 244 : LESOBL( 4)='DOM3' 245 : LESOBL( 5)='DOM4' 246 : LESOBL( 6)='DOM5' 247 : * Défo dans le repère ortho ordre ONERA 248 : LESOBL( 7)='DF11' 249 : LESOBL( 8)='DF22' 250 : LESOBL( 9)='DF33' 251 : LESOBL( 10)='DF23' 252 : LESOBL( 11)='DF13' 253 : LESOBL( 12)='DF12' 254 : * Indices d'activation 255 : LESOBL(13)='NUA1' 256 : LESOBL(14)='NUA2' 257 : LESOBL(15)='NUA3' 258 : LESOBL(16)='NUA4' 259 : LESOBL(17)='NUA5' 260 : * Déformations inélastique repère ortho ordre ONERA 261 : LESOBL(18)='EPI1' 262 : LESOBL(19)='EPI2' 263 : LESOBL(20)='EPI3' 264 : LESOBL(21)='EPI4' 265 : LESOBL(22)='EPI5' 266 : LESOBL(23)='EPI6' 267 : ELSE 267 : CALL ERREUR (19) 268 : ENDIF 269 : * 270 : ELSEIF(MATEPL.EQ.134)THEN 271 : * 272 : * 13. CAS DU MODELE DESMORAT 273 : * 274 : NBROBL=23 275 : SEGINI,NOMID 276 : * 277 : * flag pour le non local 278 : LESOBL(1) ='EPTI' 279 : * 280 : * tenseur d'endommagement 281 : * 282 : LESOBL(2) ='D11 ' 283 : LESOBL(3) ='D22 ' 284 : LESOBL(4) ='D33 ' 285 : LESOBL(5) ='D12 ' 286 : LESOBL(6) ='D13 ' 287 : LESOBL(7) ='D23 ' 288 : * 289 : * tenseur de la base d'endommagement 290 : * 291 : LESOBL(8) ='RT11' 292 : LESOBL(9) ='RT12' 293 : LESOBL(10) ='RT13' 294 : LESOBL(11) ='RT21' 295 : LESOBL(12) ='RT22' 296 : LESOBL(13) ='RT23' 297 : LESOBL(14) ='RT31' 298 : LESOBL(15) ='RT32' 299 : LESOBL(16) ='RT33' 300 : * 301 : * flag pour l'endommagement maximum 302 : LESOBL(17) ='FAIL' 303 : * 304 : * tenseur des deformations 305 : * 306 : LESOBL(18) ='EP11' 307 : LESOBL(19) ='EP22' 308 : LESOBL(20) ='EP33' 309 : LESOBL(21) ='EP12' 310 : LESOBL(22) ='EP13' 311 : LESOBL(23) ='EP23' 312 : ELSEIF(MATEPL.EQ.141)THEN 313 : * 314 : C 14 - Cas du modele 'FATSIN' 315 : C ============================================== 316 : IF (IFLAG.EQ.0) THEN 317 : NBROBL=3 318 : SEGINI,NOMID 319 : LESOBL(1)='EPTI' 320 : LESOBL(2) ='D' 321 : LESOBL(3) ='EPT0' 322 : ELSE 323 : NBROBL=1 324 : SEGINI,NOMID 325 : LESOBL(1)='EPTI' 326 : ENDIF 327 : ELSEIF(MATEPL.EQ.144)THEN 328 : * 329 : C 15 - Cas du modele 'RICRAG' 330 : C ============================================== 331 : NBROBL=26 332 : SEGINI,NOMID 333 : * 334 : * POUR LE NON LOCAL 335 : * 336 : LESOBL(1) ='EPTI' 337 : * 338 : * TENSEUR DE GLISSEMENT 339 : * 340 : LESOBL(2) ='EPI1' 341 : LESOBL(3) ='EPI2' 342 : LESOBL(4) ='EPI3' 343 : LESOBL(5) ='EPI4' 344 : LESOBL(6) ='EPI5' 345 : LESOBL(7) ='EPI6' 346 : * 347 : * TENSEUR D ECROUISSAGE CINEMATIQUE 348 : * 349 : LESOBL(8) ='ABC1' 350 : LESOBL(9) ='ABC2' 351 : LESOBL(10) ='ABC3' 352 : LESOBL(11) ='ABC4' 353 : LESOBL(12) ='ABC5' 354 : LESOBL(13) ='ABC6' 355 : * 356 : * QUELQUES VARIABLES INTERNES SCALAIRES 357 : * 358 : LESOBL(14) ='D ' 359 : LESOBL(15) ='DOND' 360 : LESOBL(16) ='DONI' 361 : LESOBL(17) ='FFFZ' 362 : LESOBL(18) ='TYPC' 363 : * 364 : * TENSEUR DE DEFORMATION TOTALES 365 : * 366 : LESOBL(19) ='EPA1' 367 : LESOBL(20) ='EPA2' 368 : LESOBL(21) ='EPA3' 369 : LESOBL(22) ='EPA4' 370 : LESOBL(23) ='EPA5' 371 : LESOBL(24) ='EPA6' 372 : * 373 : * FLAGS POUR ENDOMMAGEMENT ECCESSIF ET ECROUISSAGE 374 : * 375 : LESOBL(25) ='FZF1' 376 : LESOBL(26) ='FZF2' 377 : C +BR 378 : ELSEIF(MATEPL.EQ.157)THEN 379 : * 380 : C 16 - Cas du modele 'GLRC_DM' 381 : C ============================================== 382 : NBROBL=12 383 : SEGINI,NOMID 384 : LESOBL( 1) ='VA01' 385 : LESOBL( 2) ='VA02' 386 : LESOBL( 3) ='VA03' 387 : LESOBL( 4) ='VA04' 388 : LESOBL( 5) ='VA05' 389 : LESOBL( 6) ='VA06' 390 : LESOBL( 7) ='VA07' 391 : LESOBL( 8) ='VA08' 392 : LESOBL( 9) ='XD1 ' 393 : LESOBL(10) ='XD2 ' 394 : LESOBL(11) ='XDM ' 395 : LESOBL(12) ='XDF ' 396 : 397 : ELSEIF(MATEPL.EQ.158)THEN 398 : * 399 : C 17 - Cas du modele 'RICBET' 400 : C ============================================== 401 : NBROBL=42 402 : SEGINI,NOMID 403 : 404 : LESOBL( 1)='EPTI' 405 : 406 : LESOBL( 2)='EPI1' 407 : LESOBL( 3)='EPI2' 408 : LESOBL( 4)='EPI3' 409 : LESOBL( 5)='EPI4' 410 : LESOBL( 6)='EPI5' 411 : LESOBL( 7)='EPI6' 412 : 413 : LESOBL( 8)='ABC1' 414 : LESOBL( 9)='ABC2' 415 : LESOBL(10)='ABC3' 416 : LESOBL(11)='ABC4' 417 : LESOBL(12)='ABC5' 418 : LESOBL(13)='ABC6' 419 : 420 : LESOBL(14)='D ' 421 : 422 : LESOBL(15)='FZFT' 423 : 424 : LESOBL(16)='EPA1' 425 : LESOBL(17)='EPA2' 426 : LESOBL(18)='EPA3' 427 : LESOBL(19)='EPA4' 428 : LESOBL(20)='EPA5' 429 : LESOBL(21)='EPA6' 430 : 431 : LESOBL(22)='SIC1' 432 : LESOBL(23)='SIC2' 433 : LESOBL(24)='SIC3' 434 : LESOBL(25)='SIC4' 435 : LESOBL(26)='SIC5' 436 : LESOBL(27)='SIC6' 437 : 438 : LESOBL(28)='DEFF' 439 : 440 : LESOBL(29)='SIP1' 441 : LESOBL(30)='SIP2' 442 : LESOBL(31)='SIP3' 443 : LESOBL(32)='SIP4' 444 : LESOBL(33)='SIP5' 445 : LESOBL(34)='SIP6' 446 : 447 : LESOBL(35)='EPP1' 448 : LESOBL(36)='EPP2' 449 : LESOBL(37)='EPP3' 450 : LESOBL(38)='EPP4' 451 : LESOBL(39)='EPP5' 452 : LESOBL(40)='EPP6' 453 : 454 : LESOBL(41)='RECO' 455 : 456 : LESOBL(42)='DEPP' 457 : 458 : ELSEIF(MATEPL.EQ.159)THEN 459 : * 460 : C 18 - Cas du modele 'RICCOQ' 461 : C ============================================== 462 : NBROBL=18 463 : SEGINI,NOMID 464 : 465 : LESOBL(1) ='EPM1' 466 : LESOBL(2) ='EPM2' 467 : LESOBL(3) ='EPM3' 468 : 469 : LESOBL(4) ='EPF1' 470 : LESOBL(5) ='EPF2' 471 : LESOBL(6) ='EPF3' 472 : 473 : LESOBL(7) ='DT ' 474 : LESOBL(8) ='ZT ' 475 : LESOBL(9) ='DC ' 476 : LESOBL(10)='ZC ' 477 : 478 : LESOBL(11)='EPP1' 479 : LESOBL(12)='EPP2' 480 : LESOBL(13)='EPP3' 481 : 482 : LESOBL(14)='XEC1' 483 : LESOBL(15)='XEC2' 484 : LESOBL(16)='XEC3' 485 : 486 : LESOBL(17)='XYDT' 487 : LESOBL(18)='XYDC' 488 : 489 : C -BR 490 : ENDIF 491 : IPCOMP=NOMID 492 : SEGDES,NOMID 493 : 494 : 495 : RETURN 496 : END 497 : 498 : 499 : 500 : 501 : 502 : 503 : 504 :
© Cast3M 2003 - All rights reserved.
Disclaimer