1 : C CREEP SOURCE FANDEUR 13/01/29 21:15:19 7683 2 : SUBROUTINE CREEP (DECRA,DESWA,STATEV,serd,EC0,ESW0,P,QTILD, 3 : & TEMP,DTEMP,PREDEF,DPRED,TIME,DTIME, 4 : & CMNAME,leximp,LEND,COORDS,NSTATV,NOEL,NPT, 5 : & layer,kspt,KSTEP,KINC) 6 : IMPLICIT INTEGER(I-N) 7 : IMPLICIT REAL*8(A-H,O-Z) 8 : C----------------------------------------------------------------------- 9 : C 10 : C DESCRIPTION FONCTIONNELLE : 11 : C ------------------------- 12 : C 13 : C Version du module externe CREEP integree au code pour eviter tout 14 : C probleme d'edition de lien et pour permettre le passage de tests 15 : C de non regression. 16 : C 17 : C En utilisation normale, ce module doit etre surcharge par celui 18 : C de meme nom fourni par l'utilisateur. 19 : C 20 : C Si cette version de CREEP est appelee, cela signifie qu'il n'y a 21 : C pas eu surcharge par le module utilisateur 22 : C ATTENTION : si l'utilisateur oublie de surcharger et affecte a sa 23 : C loi le meme numero que l'un de ceux pris en charge ci-apres, il y 24 : C a risque d'erreur non controlee a l'execution. 25 : C 26 : C CONVENTION : on autorise la programmation de l'utilisateur a se 27 : C servir de KSTEP comme code retour, de la meme maniere que l'on 28 : C s'en sert ci-dessous. Regles a respecter : 29 : C 1. Pas d'initialisation superflue de KSTEP en entrant dans CREEP. 30 : C KSTEP est initialise a 1 par WCREE0, avant l'appel a CREEP. 31 : C 2. En cas d'erreur, KSTEP est affecte d'une valeur differente de 1 32 : C 33 : C----------------------------------------------------------------------- 34 : C 35 : C DESCRIPTION DE L'INTERFACE : 36 : C -------------------------- 37 : C 38 : C OUT : DECRA : REAL*8(5) 39 : C DECRA(1) : increment scalaire de deformation de 40 : C fluage, au debut du pas si LEND=0, a la fin du 41 : C pas si LEND=1. 42 : C Les composantes 2 a 5 ne sont pas exploitees 43 : C dans le cas d'une adherence a CAST3M. 44 : C 45 : C OUT : DESWA : REAL*8(5) 46 : C DESWA(1) : increment scalaire de deformation de 47 : C gonflement, au debut du pas si LEND=0, a la fin 48 : C du pas si LEND=1. 49 : C Les composantes 2 a 5 ne sont pas exploitees 50 : C dans le cas d'une adherence a CAST3M. 51 : C 52 : C Remarque : la routine CREEP evalue des VITESSES 53 : C de deformation (fluage et gonflement) au debut 54 : C ou a la fin du pas, suivant la valeur de LEND. 55 : C Les increments DECRA(1) et DESWA(1) sont ensuite 56 : C determines par les produits de ces vitesses par 57 : C le pas de temps DTIME. 58 : C 59 : C IN/OUT : STATEV : REAL*8(*), variables internes supplementaires 60 : C Il s'agit des eventuelles variables internes 61 : C s'ajoutant aux 4 pre-definies 'EC0 ', 'ESW0', 62 : C 'P ' et 'QTLD'. 63 : C Valeurs au debut du pas si LEND=0. 64 : C Peuvent etre mises a jour si LEND=1 : valeurs a 65 : C la fin du pas. 66 : C 67 : C IN/OUT : serd : REAL*8, puissance volumique de deformation 68 : C inelastique (produit contracte du tenseur des 69 : C contraintes et du tenseur des vitesses de 70 : C deformation inelastique). 71 : C Entree/sortie facultative de la routine CREEP, 72 : C non exploitee par CAST3M pour l'instant. 73 : C 74 : C IN : EC0 : REAL*8, deformation de fluage cumulee. 75 : C 1ere variable interne parmi les 4 pre-definies. 76 : C Valeur au debut du pas. 77 : C 78 : C IN : ESW0 : REAL*8, deformation de gonflement cumulee. 79 : C 2eme variable interne parmi les 4 pre-definies. 80 : C Valeur au debut du pas. 81 : C 82 : C IN : P : REAL*8, 1er invariant du tenseur des contraintes 83 : C (1/3 trace(SIGMA)). 84 : C 3eme variable interne parmi les 4 pre-definies. 85 : C Valeur au debut du pas si LEND=0, a la fin du 86 : C pas si LEND=1. 87 : C 88 : C IN : QTILD : REAL*8, contrainte equivalente de Von Mises 89 : C (SQRT(3/2 S:S) ou S deviateur des contraintes). 90 : C 4eme variable interne parmi les 4 pre-definies. 91 : C Valeur au debut du pas si LEND=0, a la fin du 92 : C pas si LEND=1. 93 : C 94 : C IN : TEMP : REAL*8, temperature a la fin du pas. 95 : C DTEMP : REAL*8, increment de temperature sur le pas. 96 : C PREDEF : REAL*8(*), vecteur des parametres externes de la 97 : C loi de comportement, valeurs a la fin du pas. 98 : C DPRED : REAL*8(*), increments des parametres externes 99 : C sur le pas de temps. 100 : C 101 : C IN : TIME : REAL*8(2) 102 : C TIME(1) : duree cumulee des iterations internes 103 : C a la fin du pas en cours. 104 : C TIME(2) : instant intermediaire absolu a la fin 105 : C du pas en cours. 106 : C DTIME : REAL*8, valeur du pas de temps en cours. 107 : C 108 : C Remarque : l'instant absolu correspondant a la 109 : C precedente iteration interne convergee est 110 : C TIME(2)-DTIME. 111 : C 112 : C IN : CMNAME : CHARACTER*16, identifiant de la loi de 113 : C comportement. 114 : C On conserve le type 'chaine de caracteres' pour 115 : C l'identifiant de la loi, afin de preserver la 116 : C compatibilite avec ABAQUS. 117 : C Dans le cas d'une adherence a CAST3M, la loi est 118 : C identifiee par le numero qui lui a ete attribue: 119 : C argument 'NUME_LOI' de l'operateur MODE. 120 : C Par convention, ce numero est encode dans les 4 121 : C derniers caracteres de la chaine, et doit etre 122 : C recupere dans CREEP par une instruction du type 123 : C K4ILOI = CMNAME(13:16) 124 : C avec K4ILOI variable locale de type CHARACTER*4. 125 : C 126 : C IN : leximp : INTEGER, indicateur d'integration par un schema 127 : C Entree non active dans le cas d'une adherence 128 : C a CAST3M. 129 : C 130 : C IN : LEND : INTEGER, indicateur de debut/fin de pas. 131 : C LEND=0 : les entrees P, QTILD sont definies au 132 : C debut du pas ; les sorties DECRA(1), DESWA(1) 133 : C sont calculees au debut du pas. 134 : C LEND=1 : les entrees P, QTILD sont definies a la 135 : C fin du pas ; les sorties DECRA(1), DESWA(1) sont 136 : C calculees a la fin du pas. 137 : C 138 : C IN : COORDS : REAL*8(*), coordonnees cartesiennes du point 139 : C d'integration courant. 140 : C 141 : C IN NSTATV : INTEGER, nombre de variables internes 142 : C supplementaires (en plus des 4 pre-definies). 143 : C 144 : C IN : NOEL : INTEGER, numero de l'element courant. 145 : C NPT : INTEGER, numero du point d'integration courant. 146 : C 147 : C IN : layer : INTEGER, numero de couche pour des coques 148 : C composites ou des solides multi-couches. 149 : C kspt : INTEGER, numero de section dans la couche 150 : C courante. 151 : C Entrees non actives dans le cas d'une adherence 152 : C a CAST3M. 153 : C 154 : C IN/OUT : KSTEP : INTEGER 155 : C Entree n'ayant pas de sens dans le cas d'une 156 : C adherence a CAST3M. 157 : C On autorise la programmation de l'utilisateur a 158 : C se servir de KSTEP comme code retour de CREEP 159 : C (cf. description fonctionnelle ci-dessus). 160 : C 161 : C IN : KINC : INTEGER, compteur d'iterations internes. 162 : C Incrementation geree par l'appelant CCREEP. 163 : C 164 : C----------------------------------------------------------------------- 165 : C 166 : C MODELES PRIS EN CHARGE par le module CREEP 'bouchon' : 167 : C ---------------------------------------------------- 168 : C 169 : C ILOI = 21 : Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21 170 : C ILOI = 23 RESTRICTION A LA FORMULATION MASSIVE 171 : C Le modele ILOI=23 comprend deux variables internes 172 : C supplementaires artificielles : 173 : C 'PSUP' <=> 'P ' et 'QSUP' <=> 'QTLD' 174 : C 175 : C ILOI = 22 : Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21 176 : C ILOI = 24 RESTRICTION A LA FORMULATION MASSIVE 177 : C Instanciation des composantes de materiau au fur et a 178 : C mesure des iterations internes (uniquement les coef. 179 : C du polynome de la loi de fluage). 180 : C Le modele comprend cette fois des parametres externes. 181 : C Les composantes de materiau sont evaluees par les 182 : C memes fonctions externes que celles utilisees par 183 : C l'operateur VARI. 184 : C MODELE 'MAQUETTE' 185 : C Le modele ILOI=24 comprend deux variables internes 186 : C supplementaires artificielles : 187 : C 'PSUP' <=> 'P ' et 'QSUP' <=> 'QTLD' 188 : C 189 : C----------------------------------------------------------------------- 190 : C Arguments de l'interface 191 : C 192 : CHARACTER*16 CMNAME 193 : C 194 : INTEGER leximp, LEND, NSTATV, NOEL, NPT, layer, kspt, 195 : & KSTEP, KINC 196 : C 197 : REAL*8 DECRA(5), DESWA(5), STATEV(*), 198 : & serd, EC0, ESW0, P, QTILD, 199 : & TEMP, DTEMP, PREDEF(*), DPRED(*), TIME(2), DTIME, 200 : & COORDS(*) 201 : C 202 : C Variables locales 203 : C 204 : C Tableaux de travail 205 : C 206 : REAL*8 VALPAR(5) 207 : C 208 : C Variables scalaires 209 : C 210 : INTEGER NBPAR 211 : C 212 : REAL*8 TFIS, TUO2, FACF, DSIU, DGRA, 213 : & XAF0, XAF1, XAF2, XAF3, XAF4, XAF5, XAF6, 214 : & EC0PT 215 : C 216 : CHARACTER*4 K4ILOI 217 : C 218 : REAL*8 FVAF0, FVAF1, FVAF2, FVAF3, FVAF4, FVAF5, FVAF6 219 : EXTERNAL FVAF0, FVAF1, FVAF2, FVAF3, FVAF4, FVAF5, FVAF6 220 : C 221 : C 222 : C------------------- Debut du code executable -------------------------- 223 : C 224 : C 225 : K4ILOI = CMNAME(13:16) 226 : C 227 : C======================================================================= 228 : C ILOI = 21 : 229 : C Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21 230 : C RESTRICTION A LA FORMULATION MASSIVE 231 : C======================================================================= 232 : C 233 : IF (K4ILOI.EQ.' 21') THEN 234 : C 235 : C........Affectation des parametres pour evaluation des coefficients 236 : C du polynome de la loi de fluage 237 : C 238 : TFIS = 1.5D+19 239 : TUO2 = 1.5D+3 240 : FACF = 1.D0 241 : DSIU = 0.95D0 242 : DGRA = 8.0D-6 243 : C 244 : C........Evaluation des coefficients du polynome de la loi de fluage 245 : C 246 : XAF0 = 0.D0 * FACF 247 : XAF1 = ( ( DEXP(0.183D0*(100.D0-DSIU)) ) * 248 : & ( (0.2031D-13+(0.67861D-32*TFIS)) / (DGRA*DGRA) ) * 249 : & ( DEXP(-45294.4D0/TUO2) ) ) 250 : & * FACF 251 : XAF2 = 1.D0 252 : XAF3 = ( ( DEXP(0.241D0*(100.D0-DSIU)) ) * 253 : & ( 0.1524D-25+(0.50825D-44*TFIS) ) * 254 : & ( DEXP(-66431.8D0/TUO2) ) ) 255 : & * FACF 256 : XAF4 = 4.5D0 257 : XAF5 = ( (1.447D-34*TFIS) * (DEXP(-3624.D0/TUO2)) ) * FACF 258 : XAF6 = 1.D0 259 : C 260 : C........Calcul de la vitesse de deformation de fluage et de l'increment 261 : C........Pas de terme de gonflement 262 : C 263 : EC0PT = XAF0 + XAF1 * (QTILD**XAF2) + XAF3 * (QTILD**XAF4) 264 : & + XAF5 * (QTILD**XAF6) 265 : DECRA(1) = EC0PT * DTIME 266 : C 267 : DESWA(1) = 0.D0 268 : C 269 : C======================================================================= 270 : C ILOI = 22 : 271 : C Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21 272 : C RESTRICTION A LA FORMULATION MASSIVE 273 : C Instanciation des composantes de materiau au fur et a mesure des 274 : C iterations internes (uniquement les coef. du polynome de la loi 275 : C de fluage). 276 : C Le modele comprend cette fois des parametres externes. 277 : C Les composantes de materiau sont evaluees par les memes fonctions 278 : C externes que celles utilisees par l'operateur VARI. 279 : C MODELE 'MAQUETTE' 280 : C======================================================================= 281 : C 282 : ELSE IF (K4ILOI.EQ.' 22') THEN 283 : C 284 : IERUT = 0 285 : C 286 : C........Affectation des parametres et evaluation des coefficients 287 : C du polynome de la loi de fluage 288 : C 289 : IF ( LEND.EQ.0 ) THEN 290 : C 291 : NBPAR = 3 292 : VALPAR(1) = PREDEF(1)-DPRED(1) 293 : VALPAR(2) = PREDEF(2)-DPRED(2) 294 : VALPAR(3) = PREDEF(3)-DPRED(3) 295 : XAF0 = FVAF0(VALPAR,NBPAR,IERUT) 296 : IF (IERUT.GT.0) THEN 297 : KSTEP = -22000-IERUT 298 : RETURN 299 : ENDIF 300 : C 301 : XAF2 = FVAF2(VALPAR,NBPAR,IERUT) 302 : IF (IERUT.GT.0) THEN 303 : KSTEP = -22000-IERUT 304 : RETURN 305 : ENDIF 306 : C 307 : XAF4 = FVAF4(VALPAR,NBPAR,IERUT) 308 : IF (IERUT.GT.0) THEN 309 : KSTEP = -22000-IERUT 310 : RETURN 311 : ENDIF 312 : C 313 : XAF5 = FVAF5(VALPAR,NBPAR,IERUT) 314 : IF (IERUT.GT.0) THEN 315 : KSTEP = -22000-IERUT 316 : RETURN 317 : ENDIF 318 : C 319 : XAF6 = FVAF6(VALPAR,NBPAR,IERUT) 320 : IF (IERUT.GT.0) THEN 321 : KSTEP = -22000-IERUT 322 : RETURN 323 : ENDIF 324 : C 325 : NBPAR = 4 326 : VALPAR(4) = PREDEF(4)-DPRED(4) 327 : XAF3 = FVAF3(VALPAR,NBPAR,IERUT) 328 : IF (IERUT.GT.0) THEN 329 : KSTEP = -22000-IERUT 330 : RETURN 331 : ENDIF 332 : C 333 : NBPAR = 5 334 : VALPAR(5) = PREDEF(5)-DPRED(5) 335 : XAF1 = FVAF1(VALPAR,NBPAR,IERUT) 336 : IF (IERUT.GT.0) THEN 337 : KSTEP = -22000-IERUT 338 : RETURN 339 : ENDIF 340 : C 341 : ELSE 342 : C 343 : NBPAR = 3 344 : XAF0 = FVAF0(PREDEF,NBPAR,IERUT) 345 : IF (IERUT.GT.0) THEN 346 : KSTEP = -22000-IERUT 347 : RETURN 348 : ENDIF 349 : C 350 : XAF2 = FVAF2(PREDEF,NBPAR,IERUT) 351 : IF (IERUT.GT.0) THEN 352 : KSTEP = -22000-IERUT 353 : RETURN 354 : ENDIF 355 : C 356 : XAF4 = FVAF4(PREDEF,NBPAR,IERUT) 357 : IF (IERUT.GT.0) THEN 358 : KSTEP = -22000-IERUT 359 : RETURN 360 : ENDIF 361 : C 362 : XAF5 = FVAF5(PREDEF,NBPAR,IERUT) 363 : IF (IERUT.GT.0) THEN 364 : KSTEP = -22000-IERUT 365 : RETURN 366 : ENDIF 367 : C 368 : XAF6 = FVAF6(PREDEF,NBPAR,IERUT) 369 : IF (IERUT.GT.0) THEN 370 : KSTEP = -22000-IERUT 371 : RETURN 372 : ENDIF 373 : C 374 : NBPAR = 4 375 : XAF3 = FVAF3(PREDEF,NBPAR,IERUT) 376 : IF (IERUT.GT.0) THEN 377 : KSTEP = -22000-IERUT 378 : RETURN 379 : ENDIF 380 : C 381 : NBPAR = 5 382 : XAF1 = FVAF1(PREDEF,NBPAR,IERUT) 383 : IF (IERUT.GT.0) THEN 384 : KSTEP = -22000-IERUT 385 : RETURN 386 : ENDIF 387 : C 388 : ENDIF 389 : C 390 : C........Calcul de la vitesse de deformation de fluage et de l'increment 391 : C........Pas de terme de gonflement 392 : C 393 : EC0PT = XAF0 + XAF1 * (QTILD**XAF2) + XAF3 * (QTILD**XAF4) 394 : & + XAF5 * (QTILD**XAF6) 395 : DECRA(1) = EC0PT * DTIME 396 : C 397 : DESWA(1) = 0.D0 398 : C 399 : C======================================================================= 400 : C ILOI = 23 : 401 : C Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21 402 : C RESTRICTION A LA FORMULATION MASSIVE 403 : C Equivalent au modele ILOI=21 404 : C Variables internes supplementaires 'PSUP' et 'QSUP' 405 : C======================================================================= 406 : C 407 : ELSE IF (K4ILOI.EQ.' 23') THEN 408 : C 409 : C........Affectation des parametres pour evaluation des coefficients 410 : C du polynome de la loi de fluage 411 : C 412 : TFIS = 1.5D+19 413 : TUO2 = 1.5D+3 414 : FACF = 1.D0 415 : DSIU = 0.95D0 416 : DGRA = 8.0D-6 417 : C 418 : C........Evaluation des coefficients du polynome de la loi de fluage 419 : C 420 : XAF0 = 0.D0 * FACF 421 : XAF1 = ( ( DEXP(0.183D0*(100.D0-DSIU)) ) * 422 : & ( (0.2031D-13+(0.67861D-32*TFIS)) / (DGRA*DGRA) ) * 423 : & ( DEXP(-45294.4D0/TUO2) ) ) 424 : & * FACF 425 : XAF2 = 1.D0 426 : XAF3 = ( ( DEXP(0.241D0*(100.D0-DSIU)) ) * 427 : & ( 0.1524D-25+(0.50825D-44*TFIS) ) * 428 : & ( DEXP(-66431.8D0/TUO2) ) ) 429 : & * FACF 430 : XAF4 = 4.5D0 431 : XAF5 = ( (1.447D-34*TFIS) * (DEXP(-3624.D0/TUO2)) ) * FACF 432 : XAF6 = 1.D0 433 : C 434 : C........Calcul de la vitesse de deformation de fluage et de l'increment 435 : C........Pas de terme de gonflement 436 : C 437 : EC0PT = XAF0 + XAF1 * (QTILD**XAF2) + XAF3 * (QTILD**XAF4) 438 : & + XAF5 * (QTILD**XAF6) 439 : DECRA(1) = EC0PT * DTIME 440 : C 441 : DESWA(1) = 0.D0 442 : C 443 : C........MAJ des variables internes supplementaires le cas echeant 444 : C 445 : IF ( LEND.EQ.1 ) THEN 446 : STATEV(1) = P 447 : STATEV(2) = QTILD 448 : ENDIF 449 : C 450 : C======================================================================= 451 : C ILOI = 24 : 452 : C Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21 453 : C RESTRICTION A LA FORMULATION MASSIVE 454 : C Instanciation des composantes de materiau au fur et a mesure des 455 : C iterations internes (uniquement les coef. du polynome de la loi 456 : C de fluage). 457 : C Le modele comprend cette fois des parametres externes. 458 : C Les composantes de materiau sont evaluees par les memes fonctions 459 : C externes que celles utilisees par l'operateur VARI. 460 : C MODELE 'MAQUETTE' 461 : C Equivalent au modele ILOI=22 462 : C Variables internes supplementaires 'PSUP' et 'QSUP' 463 : C======================================================================= 464 : C 465 : ELSE IF (K4ILOI.EQ.' 24') THEN 466 : C 467 : C........Affectation des parametres et evaluation des coefficients 468 : C du polynome de la loi de fluage 469 : C MAJ des variables internes supplementaires le cas echeant 470 : C 471 : IF ( LEND.EQ.0 ) THEN 472 : C 473 : NBPAR = 3 474 : VALPAR(1) = PREDEF(1)-DPRED(1) 475 : VALPAR(2) = PREDEF(2)-DPRED(2) 476 : VALPAR(3) = PREDEF(3)-DPRED(3) 477 : XAF0 = FVAF0(VALPAR,NBPAR,IERUT) 478 : IF (IERUT.GT.0) THEN 479 : KSTEP = -24000-IERUT 480 : RETURN 481 : ENDIF 482 : C 483 : XAF2 = FVAF2(VALPAR,NBPAR,IERUT) 484 : IF (IERUT.GT.0) THEN 485 : KSTEP = -24000-IERUT 486 : RETURN 487 : ENDIF 488 : C 489 : XAF4 = FVAF4(VALPAR,NBPAR,IERUT) 490 : IF (IERUT.GT.0) THEN 491 : KSTEP = -24000-IERUT 492 : RETURN 493 : ENDIF 494 : C 495 : XAF5 = FVAF5(VALPAR,NBPAR,IERUT) 496 : IF (IERUT.GT.0) THEN 497 : KSTEP = -24000-IERUT 498 : RETURN 499 : ENDIF 500 : C 501 : XAF6 = FVAF6(VALPAR,NBPAR,IERUT) 502 : IF (IERUT.GT.0) THEN 503 : KSTEP = -24000-IERUT 504 : RETURN 505 : ENDIF 506 : C 507 : NBPAR = 4 508 : VALPAR(4) = PREDEF(4)-DPRED(4) 509 : XAF3 = FVAF3(VALPAR,NBPAR,IERUT) 510 : IF (IERUT.GT.0) THEN 511 : KSTEP = -24000-IERUT 512 : RETURN 513 : ENDIF 514 : C 515 : NBPAR = 5 516 : VALPAR(5) = PREDEF(5)-DPRED(5) 517 : XAF1 = FVAF1(VALPAR,NBPAR,IERUT) 518 : IF (IERUT.GT.0) THEN 519 : KSTEP = -24000-IERUT 520 : RETURN 521 : ENDIF 522 : C 523 : ELSE 524 : C 525 : NBPAR = 3 526 : XAF0 = FVAF0(PREDEF,NBPAR,IERUT) 527 : IF (IERUT.GT.0) THEN 528 : KSTEP = -24000-IERUT 529 : RETURN 530 : ENDIF 531 : C 532 : XAF2 = FVAF2(PREDEF,NBPAR,IERUT) 533 : IF (IERUT.GT.0) THEN 534 : KSTEP = -24000-IERUT 535 : RETURN 536 : ENDIF 537 : C 538 : XAF4 = FVAF4(PREDEF,NBPAR,IERUT) 539 : IF (IERUT.GT.0) THEN 540 : KSTEP = -24000-IERUT 541 : RETURN 542 : ENDIF 543 : C 544 : XAF5 = FVAF5(PREDEF,NBPAR,IERUT) 545 : IF (IERUT.GT.0) THEN 546 : KSTEP = -24000-IERUT 547 : RETURN 548 : ENDIF 549 : C 550 : XAF6 = FVAF6(PREDEF,NBPAR,IERUT) 551 : IF (IERUT.GT.0) THEN 552 : KSTEP = -24000-IERUT 553 : RETURN 554 : ENDIF 555 : C 556 : NBPAR = 4 557 : XAF3 = FVAF3(PREDEF,NBPAR,IERUT) 558 : IF (IERUT.GT.0) THEN 559 : KSTEP = -24000-IERUT 560 : RETURN 561 : ENDIF 562 : C 563 : NBPAR = 5 564 : XAF1 = FVAF1(PREDEF,NBPAR,IERUT) 565 : IF (IERUT.GT.0) THEN 566 : KSTEP = -24000-IERUT 567 : RETURN 568 : ENDIF 569 : C 570 : STATEV(1) = P 571 : STATEV(2) = QTILD 572 : C 573 : ENDIF 574 : C 575 : C........Calcul de la vitesse de deformation de fluage et de l'increment 576 : C........Pas de terme de gonflement 577 : C 578 : EC0PT = XAF0 + XAF1 * (QTILD**XAF2) + XAF3 * (QTILD**XAF4) 579 : & + XAF5 * (QTILD**XAF6) 580 : DECRA(1) = EC0PT * DTIME 581 : C 582 : DESWA(1) = 0.D0 583 : C 584 : C======================================================================= 585 : C Modele externe non reconnu 586 : C======================================================================= 587 : C 588 : ELSE 589 : C 590 : KSTEP = 0 591 : C 592 : ENDIF 593 : C 594 : RETURN 595 : END 596 : 597 : 598 :
© Cast3M 2003 - All rights reserved.
Disclaimer