1 : C UBIQUI SOURCE AM 08/12/19 21:17:10 6237 2 : SUBROUTINE UBIQUI(DDAUX,CMATE,VALMAT,VALCAR,N2EL,N2PTEL, 3 : 1 IB,IGAU,EPAIST,MELE,NPINT,SECT,LHOOK,TXR,XLOC, 4 : 2 XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI,SIG0,DEPST,VAR0,XMAT, 5 : 3 NBPGAU,NCOMAT,XCAR,DSIGT,SIGF,VARF,DEFP,KERRE 6 : 4, ecou,necou,iecou) 7 : C 8 : C---------------------------------------------------------------------- 9 : C PLASTICITE MODELE UBIQUITOUS 10 : C 11 : C ENTREES 12 : C SIG0(NSTRS) = CONTRAINTES INITIALES 13 : C NSTRS = NOMBRE DE CONTRAINTES 14 : C DEPST(NSTRS)= INCREMENT DE DEFORMATIONS TOTALES 15 : C VAR0(NVARI) = VARIABLES INTERNES DEBUT 16 : C ( 1 ) = EPSE 17 : C ( 2 ) = EPN1 18 : C ( 3 ) = GAP1 19 : C ( 4 ) = LAM1 20 : C ( 5 ) = EPN2 21 : C ( 6 ) = GAP2 22 : C ( 7 ) = LAM2 23 : C XMAT(NCOMAT)= COMPOSANTES DE MATERIAU 24 : C NCOMAT = NOMBRE DE COMPOSANTES DE MATERIAU 25 : C XCAR(ICARA) = CARACTERISTIQUES 26 : * CMATE = NOM DU MATERIAU 27 : * VALMAT= TABLEAU DE CARACTERISTIQUES DU MATERIAU 28 : * VALCAR= TABLEAU DE CARACTERISTIQUES GEOMETRIQUES 29 : * N2EL = NBRE D ELEMENTS DANS SEGMENT DE HOOKE 30 : * N2PTEL= NBRE DE POINTS DANS SEGMENT DE HOOKE 31 : * MFR = NUMERO DE LA FORMULATION 32 : * IFOU = OPTION DE CALCUL 33 : * IB = NUMERO DE L ELEMENT COURANT 34 : * IGAU = NUMERO DU POINT COURANT 35 : * EPAIST= EPAISSEUR 36 : * NBPGAU= NBRE DE POINTS DE GAUSS 37 : * MELE = NUMERO DE L ELEMENT FINI 38 : * NPINT = NBRE DE POINTS D INTEGRATION 39 : * NBGMAT= NBRE DE POINTS DANS SEGMENT DE CARACTERISTIQUES 40 : * NELMAT= NBRE D ELEMENTS DANS SEGMENT DE CARACTERISTIQUES 41 : * SECT = SECTION 42 : * LHOOK = TAILLE DE LA MATRICE DE HOOKE 43 : * TXR,XLOC,XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI = TABLEAUX UTILISES 44 : * UTILISES POUR LE CALCUL DE LA MATRICE DE HOOKE 45 : * 46 : C SORTIES 47 : C SIGF(NSTRS) = CONTRAINTES FINALES 48 : C VARF(NVARI) = VARIABLES INTERNES FINALES 49 : C DEFP(NSTRS) = DEFORMATIONS PLASTIQUES 50 : C KERRE = 0 TOUT OK 51 : C 1 SI DLAMBDA NEGATIF 52 : C 2 NOMBRE MAX D ITERATIONS INTERNES DEPASSE 53 : C 21 ON NE TROUVE PAS L INTERSECTION AVEC LA SURFACE DE CHARGE 54 : C 22 SIG0 A L EXTERIEUR DE LA SURFACE DE CHARGE 55 : C 56 : C----------------------------------------------------------------------- 57 : C VARIABLES PASSEES PAR LES COMMONS COPTIO , ECOU ET NECOU 58 : C 59 : C IFOUR INDICE DU TYPE DE PROBLEME 60 : C -2 CONTRAINTES PLANES 61 : C -1 DEFORMATIONS PLANES 62 : C 0 AXISYMETRIQUE 63 : C 1 SERIE DE FOURIER 64 : C 2 TRIDIMENSIONNEL 65 : C ITYP TYPE DE FORMULATION MECANIQUE 66 : C ITYP=1 CAS DES ELEMENTS MASSIFS 67 : C ITYP=2 CAS DES COQUES 68 : C ITYP=3 CAS DES MEMBRANES 69 : C ITYP=4 CAS DES CABLES ET DES BARRES 70 : C ITYP=5 CAS QUELCONQUE 71 : C ITYP=6 CAS DES CONTRAINTES PLANES 72 : C ITYP=7 CAS DES COQUES A NU=0. OU CONTRAINTES PLANES 73 : C ITYP=8 CAS DES MEMBRANES A NU=0. OU CONTRAINTES PLANES 74 : C ITYP=9 CAS DES COQUES EPAISSES 75 : C ITYP=10 CAS DES JOINTS 76 : C ITYP=11 CAS DES POUTRES 77 : C ITYP=12 CAS DES TUYAUX 78 : C ITYP=13 CAS DES COQUES AVEC CISAILLEMENT TRANSVERSE 79 : C 80 : C----------------------------------------------------------------------- 81 : C 82 : IMPLICIT INTEGER(I-N) 83 : IMPLICIT REAL*8(A-H,O-Z) 84 : -INC CCOPTIO 85 : SEGMENT IECOU 86 : * COMMON/IECOU/NYOG,NYNU,NYALFA,NYSMAX,NYN,NYM,NYKK, 87 : INTEGER icow1,icow2,icow3,icow4,icow5,icow6,icow7, 88 : C INTEGER NYOG, NYNU, NYALFA,NYSMAX,NYN, NYM, NYKK, 89 : 1 icow8,icow9,icow10,icow11,icow12,icow13,icow14,icow15,icow16, 90 : C . NYALF1,NYBET1,NYR, NYA, NYRHO,NSIGY, NNKX, NYKX, IND, 91 : 2 icow17,icow18,icow19,icow20,icow21,icow22,icow23,icow24, 92 : C . NSOM, NINV, NINCMA,NCOMP, JELEM, LEGAUS,INAT, NCXMAT, 93 : 3 icow25,icow26,icow27,icow28,icow29,icow30,icow31, 94 : C . LTRAC, MFR, IELE, NHRM, NBNN, NBELEM,ICARA, 95 : 4 icow32,icow33,NSTRS,MFR,NBGMAT,NELMAT,icow38, 96 : C . LW2, NDEF, NSTRSS,MFR1, NBGMAT,NELMAT,MSOUPA, 97 : 5 icow39,icow40,icow41,icow42,icow43,icow44 98 : C . NUMAT1,LENDO, NBBB, NNVARI,KERR1, MELEME 99 : INTEGER icow45,icow46,icow47,icow48,icow49,icow50, 100 : . icow51,icow52,icow53,icow54,icow55,icow56 101 : . icow57,icow58 102 : ENDSEGMENT 103 : SEGMENT ECOU 104 : *** COMMON/ECOU/TEST,ALFAH, 105 : C REAL*8 ecow00,ecow0, 106 : REAL*8 TEST, ALFAH, 107 : 1 ecow1,ecow2,ecow3(6),ecow4(9),ecow5(6), 108 : C 1 HPAS,TEMPS,COVNMS(6),VECPRO(9),VALPRO(6), 109 : 2 ecow6(12),ecow7(6),ecow8(6),ecow9(6),ecow10(6),ecow11(6), 110 : 2 ecow12(6), 111 : C 2 CVNMSD(12),STOT(6),SIGEL(6),DSIGP(6),SIGT(6),W1(6),W2(6), 112 : 1 ecow13(6),ecow14(6),ecow15(12),ecow16(3), 113 : C 1 DALPHA(6),EPSPLA(6),E(12),XINV(3), 114 : 2 ecow17(6),ecow18(6),ecow19,ecow20 115 : C 2 SIPLAD(6),DSIGP0(6),TET,TETI 116 : ENDSEGMENT 117 : C COMMON/ECOU/TEST,ALFAH,ZOZO(112) 118 : SEGMENT NECOU 119 : * COMMON/NECOU/NCOURB,IPLAST,IT,IMAPLA,ISOTRO, 120 : C INTEGER ncow1, ncow2,ncow3,ncow4, ncow5, 121 : INTEGER NCOURB,IPLAST,IT, IMAPLA,ISOTRO, 122 : 1 ITYP, IFOURB, ncow8, 123 : C . ITYP, IFOURB,IFLUAG, 124 : 2 ncow9,ncow10,ncow11,ncow12,ncow13, 125 : C . ICINE,ITHER, IFLUPL,ICYCL, IBI, 126 : 3 ncow14,ncow15,ncow16, 127 : C . JFLUAG,KFLUAG,LFLUAG, 128 : 4 ncow17,ncow18,ncow19,ncow20,ncow21 129 : C . IRELAX,JNTRIN,MFLUAG,JSOUFL,JGRDEF 130 : ENDSEGMENT 131 : C COMMON/NECOU/NCOURB,IPLAST,IT,IMAPLA,ISOTRO,ITYP,IFOURB,IZOZO(14) 132 : C 133 : DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*) 134 : DIMENSION SIGF(*),VARF(*),DEFP(*) 135 : DIMENSION SIGMA(6),DSIGMA(6),STOT(6),SIGFIN(6) 136 : DIMENSION XLAMBD(6),DEFPLA(6),DSIGT(*) 137 : DIMENSION SPHER(6),AUXIL(6),DSIGZE(6) 138 : DIMENSION VALMAT(*),VALCAR(*) 139 : DIMENSION TXR(IDIM,*),CRIGI(12) 140 : DIMENSION DDAUX(LHOOK,*),DDHOMU(LHOOK,*) 141 : DIMENSION XLOC(IDIM,*),XGLOB(IDIM,*) 142 : DIMENSION D1HOOK(LHOOK,*),ROTHOO(LHOOK,*) 143 : CHARACTER*8 CMATE 144 : C 145 : C TEST DE CONSISTANCE DES DONNEES 146 : C 147 : DO 100 I=1,NCOMAT 148 : IF(XMAT(I).LT.0.D0) KERRE=538 149 : 100 CONTINUE 150 : IF(KERRE.EQ.538) RETURN 151 : C 152 : TRA1=0.D0 153 : PHI1=0.D0 154 : PSI1=0.D0 155 : COHE1=0.D0 156 : TPHI1=0.D0 157 : TPSI1=0.D0 158 : HACHE1=0.D0 159 : TRA2=0.D0 160 : PHI2=0.D0 161 : PSI2=0.D0 162 : COHE2=0.D0 163 : TPHI2=0.D0 164 : TPSI2=0.D0 165 : HACHE2=0.D0 166 : CO21=0.D0 167 : CO22=0.D0 168 : SI21=0.D0 169 : SI22=0.D0 170 : C 171 : YOUN=XMAT(1) 172 : XNU =XMAT(2) 173 : NCRI= INT(XMAT(3)) 174 : ANG1=XMAT(4) 175 : TRA1=XMAT(5) 176 : PHI1=XMAT(6) 177 : PSI1=XMAT(7) 178 : RHO =XMAT(8) 179 : ALPH=XMAT(9) 180 : ANG2=XMAT(10) 181 : TRA2=XMAT(11) 182 : PHI2=XMAT(12) 183 : PSI2=XMAT(13) 184 : C 185 : IF(IIMPI.EQ.28) WRITE(IOIMP,3000) (XMAT(I),I=1,13) 186 : C 187 : C PETIT TEST SUR NU POUR CERTAINS CAS 188 : C 189 : IF(MFR.EQ.3.AND.IFOURB.EQ.-2.AND.XNU.NE.0.D0) THEN 190 : KERRE=38 191 : RETURN 192 : ENDIF 193 : C 194 : C 195 : C INITIALISATION 196 : C 197 : NCAS=2 198 : UNIT=0.0174532925199432957692D0 199 : C 200 : C 201 : C CALCUL ET TEST DES PARAMETRES DU MATERIAU UBIQUITOUS 202 : C 203 : C TEST D'ERREUR SUR LE NOMBRE DE CRITERE 204 : C----------------------------------------------- 205 : C ----------------------------------------------------------- 206 : C ON SE LIMITE AUX DEFORMATIONS PLANES , CONTRAINTES PLANES 207 : C ET AXISYMETRIQUE 208 : C------------------------------------------------------------ 209 : IF(IFOURB.LE.0) GOTO 101 210 : C 211 : KERRE=539 212 : RETURN 213 : C------------------------------------------------------------- 214 : 101 IF(IFOURB.EQ.-2) NCAS=1 215 : IF(NCRI.EQ.1.OR.NCRI.EQ.2) GO TO 102 216 : C 217 : KERRE=540 218 : RETURN 219 : C---------------------------------------------------------------- 220 : C TEST SUR L'EGALITE DE ANG1 ET ANG2 (CAS DE DEUX CRITERES) 221 : C---------------------------------------------------------------- 222 : 102 IF (NCRI.EQ.1) GO TO 104 223 : ANG=MIN (ABS(ANG1-ANG2),ABS(ABS(ANG2-ANG1)-180.D0)) 224 : IF (ANG.GT.1D-3) GO TO 103 225 : C 226 : KERRE=541 227 : RETURN 228 : C------------------------------------ 229 : C TEST SUR LA VALEUR DE ANG2 230 : C------------------------------------ 231 : 103 IF(ANG2.GE.0.D0.AND.ANG2.LE.180.D0) GO TO 104 232 : C 233 : KERRE=542 234 : RETURN 235 : C 236 : C------ CAS D'UN SEUL CRITERE 237 : C 238 : C------------------------------------ 239 : C TEST SUR LA VALEUR DE ANG1 240 : C------------------------------------ 241 : 104 IF(ANG1.GE.0.D0.AND.ANG1.LE.180.D0) GO TO 105 242 : C 243 : KERRE=543 244 : RETURN 245 : C=============================================== 246 : C CALCUL DES PARAMETRES GEOMETRIQUES 247 : C=============================================== 248 : 105 IF (NCRI.EQ.1) GO TO 106 249 : C 250 : C CAS DE DEUX CRITERES 251 : C 252 : ANG2=ANG2*UNIT 253 : CO22=COS(ANG2) 254 : SI22=SIN(ANG2) 255 : ANG1=ANG1*UNIT 256 : CO11=COS(ANG1) 257 : SI11=SIN(ANG1) 258 : CO21=COS(ANG2-ANG1) 259 : SI21=SIN(ANG2-ANG1) 260 : GO TO 300 261 : C 262 : C CAS D'UN SEUL CRITERE 263 : C 264 : 106 ANG1=ANG1*UNIT 265 : CO11=COS(ANG1) 266 : SI11=SIN(ANG1) 267 : C 268 : 300 CONTINUE 269 : C*************************************************************** 270 : IF (NCRI.EQ.1) GO TO 303 271 : TRA2=XMAT(11) 272 : PHI2=XMAT(12) 273 : PSI2=XMAT(13) 274 : C--------------------------------------------------------------- 275 : C TEST SUR PHI2 (PHI2 COMPRIS ENTRE 0 ET 90 STRICTEMENT) 276 : C--------------------------------------------------------------- 277 : IF(PHI2.GE.0.D0.AND.PHI2.LT.90.D0) GO TO 301 278 : KERRE=544 279 : RETURN 280 : C--------------------------------------------------------------- 281 : C TEST SUR PSI2 (PSI2 COMPRIS ENTRE 0 ET 90 STRICTEMENT) 282 : C--------------------------------------------------------------- 283 : 301 IF(PSI2.GE.0.D0.AND.PSI2.LT.90.D0) GO TO 302 284 : KERRE=545 285 : RETURN 286 : C-------------------------------------------- 287 : C TEST SUR TRA2 (TRA2 POSITIF OU NUL) 288 : C-------------------------------------------- 289 : 302 IF(TRA2.GE.0.D0) GO TO 303 290 : KERRE=546 291 : RETURN 292 : C 293 : C--------------------------------------------------------------- 294 : C TEST SUR PHI1 (PHI1 COMPRIS ENTRE 0 ET 90 STRICTEMENT) 295 : C--------------------------------------------------------------- 296 : 303 TRA1=XMAT(5) 297 : PHI1=XMAT(6) 298 : PSI1=XMAT(7) 299 : 300 : IF(PHI1.GE.0.D0.AND.PHI1.LT.90.D0) GO TO 304 301 : KERRE=547 302 : RETURN 303 : C--------------------------------------------------------------- 304 : C TEST SUR PSI1 (PSI1 COMPRIS ENTRE 0 ET 90 STRICTEMENT) 305 : C--------------------------------------------------------------- 306 : 304 IF(PSI1.GE.0.D0.AND.PSI1.LT.90.D0) GO TO 305 307 : KERRE=548 308 : RETURN 309 : C----------------------------------------------------- 310 : C TEST SUR TRA1 (TRA1 DOIT ETRE POSTIF OU NUL) 311 : C----------------------------------------------------- 312 : 305 IF(TRA1.GE.0.D0) GO TO 306 313 : KERRE=549 314 : RETURN 315 : C=========================================================== 316 : C CALCUL DES PARAMETRES PHISIQUES 317 : C=========================================================== 318 : 306 IF(NCRI.EQ.1) GO TO 307 319 : C--------------------------- 320 : C CALCUL DES PARAMETRES 2 321 : C--------------------------- 322 : TPSI2=TAN(PSI2*UNIT) 323 : TPHI2=TAN(PHI2*UNIT) 324 : COHE2=TPHI2*TRA2 325 : C------------------------------------------------------------ 326 : C CALCUL DE HACHE2(CAS=1 CONT PLANE;CAS=2 DEFO PLANE) 327 : C------------------------------------------------------------ 328 : IF (NCAS.EQ.1) HACHE2=TPSI2*TPHI2*YOUN/(1.D0-XNU* 329 : 1 XNU)+YOUN/(XNU+1.D0) 330 : IF (NCAS.EQ.2) HACHE2=TPSI2*TPHI2*YOUN*(1.D0-XNU)/(1.D0+XNU)/ 331 : 1 (1.D0-2.D0*XNU)+YOUN/(XNU+1.D0) 332 : C--------------------------------------- 333 : C------------------------------ 334 : C CAS D'UN SEUL CRITERE 335 : C------------------------------ 336 : C--------------------------- 337 : C CALCUL DES PARAMETRES 1 338 : C--------------------------- 339 : 307 TPSI1=TAN(PSI1*UNIT) 340 : TPHI1=TAN(PHI1*UNIT) 341 : COHE1=TPHI1*TRA1 342 : C------------------------------------------------------------ 343 : C CALCUL DE HACHE1(CAS=1 CONT PLANE;CAS=2 DEFO PLANE) 344 : C------------------------------------------------------------ 345 : IF (NCAS.EQ.1) HACHE1=TPSI1*TPHI1*YOUN/(1.D0-XNU* 346 : 1 XNU)+YOUN/(XNU+1.D0) 347 : IF (NCAS.EQ.2) HACHE1=TPSI1*TPHI1*YOUN*(1.D0-XNU)/(1.D0+XNU)/ 348 : 1 (1.D0-2.D0*XNU)+YOUN/(XNU+1.D0) 349 : C 350 : C IMPRESSION DES VARIABLES 351 : C 352 : IF(IIMPI.EQ.28) THEN 353 : WRITE(IOIMP,1000) 354 : WRITE(IOIMP,1001) NCRI,NCAS 355 : WRITE(IOIMP,1002) ANG1,TRA1,COHE1 356 : WRITE(IOIMP,1003) PHI1,PSI1,HACHE1 357 : WRITE(IOIMP,1004) ANG2,TRA2,COHE2 358 : WRITE(IOIMP,1005) PHI2,PSI2,HACHE2 359 : ENDIF 360 : C 361 : C CALCUL DE L INCREMENT DE CONTRAINTE 362 : C 363 : ifou=ifourb 364 : nstrbi=nstrs 365 : nbgmab=nbgmat 366 : nlmatb=nelmat 367 : mfrbi=mfr 367 : CALL CALSIG(DEPST,DDAUX,NSTRbi,CMATE,VALMAT,VALCAR, 368 : 1 N2EL,N2PTEL,MFR,IFOU,IB,IGAU,EPAIST,NBPGAU, 369 : 2 MELE,NPINT,NBGMAb,NEMATb,SECT,LHOOK,TXR,XLOC, 370 : 3 XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI,DSIGT,IRTD) 371 : 372 : * 373 : mfr=mfrbi 374 : ifourb=ifou 375 : nstrs=nstrbi 376 : nbgmat=nbgmab 377 : nelmat=nlmatb 378 : IF(IRTD.NE.1) THEN 379 : KERRE=69 380 : GOTO 1900 381 : ENDIF 382 : * 383 : C 384 : C ECOULEMENT PLASTIQUE 385 : C 386 : MCOD=1 387 : nstrbi=nstrs 388 : mfrbi=mfr 388 : CALL VISAVI(SIG0,DSIGT,VAR0,SIGMA,DSIGMA,SPHER,AUXIL, 389 : . SIGF,DEFP,VARF,SIGFIN,DEFPLA, 390 : . DSIGZE,ICENT2,MCOD,IBOU,MFRbi,NSTRbi,XCAR,CMATE,ecou,necou) 391 : mfr=mfrbi 392 : nstrs=nstrbi 393 : IF(ITYP.EQ.0) THEN 394 : KERRE=269 395 : RETURN 396 : ENDIF 397 : C 398 : C CAS DES CONTRAINTES PLANES : 399 : C IFOURB = -2 400 : C AVEC DEUX TYPES DE FORMULATION MECANIQUE : 401 : C CAS DES COQUES : ITYP = 2 ==> ALFAH = 0 402 : C OU ITYP = 7 ==> ALFAH = 0 403 : C CAS DES CONTRAINTES PLANES : ITYP = 6 404 : C 405 : C DANS LE CAS DES COQUES ITYP = 2 , NOUS METTONS ALFAH = 0 406 : C 407 : IF(ITYP.EQ.2.OR.ITYP.EQ.7) ALFAH=0.D0 408 : C 409 : C INITIALISATIONS DES PARAMETRES 410 : C 411 : EPSE =VAR0(1) 412 : XLAMBD(1)=VAR0(2) 413 : XLAMBD(2)=VAR0(3) 414 : XLAMBD(5)=VAR0(4) 415 : XLAMBD(3)=VAR0(5) 416 : XLAMBD(4)=VAR0(6) 417 : XLAMBD(6)=VAR0(7) 418 : C 419 : DO 200 IBA=1,IBOU 420 : STOT(IBA)=SIGMA(IBA)+DSIGMA(IBA) 421 : 200 CONTINUE 422 : IF(IIMPI.EQ.28) THEN 423 : WRITE(IOIMP,2000) 424 : WRITE(IOIMP,2001) (SIGMA(I),I=1,IBOU) 425 : WRITE(IOIMP,2002) (DSIGMA(I),I=1,IBOU) 426 : WRITE(IOIMP,2003) (STOT(I),I=1,IBOU) 427 : ENDIF 428 : C 429 : C ECOULEMENT PLASTIQUE 430 : C 430 : CALL ECUBI(SIGMA,DSIGMA,XLAMBD,DEFPLA,YOUN,XNU, 431 : . NCRI,ANG1,TRA1,TPHI1,TPSI1,HACHE1,COHE1,CO11,SI11, 432 : . NCAS,ANG2,TRA2,TPHI2,TPSI2,HACHE2,COHE2,CO22,SI22, 433 : . CO21,SI21, 434 : . SIGFIN,IDAM,KERRE) 435 : C 436 : IF(IIMPI.EQ.28) THEN 437 : WRITE(IOIMP,2000) 438 : WRITE(IOIMP,2001) (SIGMA(I),I=1,IBOU) 439 : WRITE(IOIMP,2002) (DSIGMA(I),I=1,IBOU) 440 : WRITE(IOIMP,2003) (STOT(I),I=1,IBOU) 441 : WRITE(IOIMP,2004) (SIGFIN(I),I=1,IBOU) 442 : WRITE(IOIMP,2005) (DEFPLA(I),I=1,IBOU) 443 : WRITE(IOIMP,2006) (XLAMBD(I),I=1,6) 444 : ENDIF 445 : C 446 : IF(KERRE.NE.0) RETURN 447 : C 448 : C PAS DE PLASTICITE COMPORTEMENT ELASTIQUE 449 : C 450 : IF(IDAM.EQ.0) THEN 451 : DO 400 IBA=1,IBOU 452 : SIGFIN(IBA)=STOT(IBA) 453 : DEFPLA(IBA)=0.D0 454 : 400 CONTINUE 455 : ENDIF 456 : C 457 : C NOUS AVONS ENDOMMAGE LE MATERIAU 458 : C 459 : C MISE A JOUR DES VARIABLES INTERNES 460 : C 461 : VARF(1)=EPSE+(SQRT(XLAMBD(1)*XLAMBD(1)+XLAMBD(2)*XLAMBD(2)+ 462 : . XLAMBD(3)*XLAMBD(3)+XLAMBD(4)*XLAMBD(4))) 463 : VARF(2)=XLAMBD(1) 464 : VARF(3)=XLAMBD(2) 465 : VARF(4)=XLAMBD(5) 466 : VARF(5)=XLAMBD(3) 467 : VARF(6)=XLAMBD(4) 468 : VARF(7)=XLAMBD(6) 469 : C 470 : MCOD=2 471 : nstrbi=nstrs 472 : mfrbi=mfr 472 : CALL VISAVI(SIG0,DSIGT,VAR0,SIGMA,DSIGMA,SPHER,AUXIL, 473 : . SIGF,DEFP,VARF,SIGFIN,DEFPLA, 474 : . DSIGZE,ICENT2,MCOD,IBOU,MFRbi,NSTRbi,XCAR,CMATE,ecou,necou) 475 : mfr=mfrbi 476 : nstrs=nstrbi 477 : C 478 : C LES FORMATS D IMPRESSION 479 : C 480 : 1000 FORMAT(1X,'CONSTANTES DU MATERIAU UBIQUITOUS') 481 : 1001 FORMAT(1X,'NCRI =',I4,3X,'NCAS =',I4) 482 : 1002 FORMAT(1X,'ANG1 =',E12.5,2X,'TRA1 =',E12.5,2X,'COHE1 =',E12.5) 483 : 1003 FORMAT(1X,'PHI1 =',E12.5,2X,'PSI1 =',E12.5,2X,'HACH1 =',E12.5) 484 : 1004 FORMAT(1X,'ANG2 =',E12.5,2X,'TRA2 =',E12.5,2X,'COHE2 =',E12.5) 485 : 1005 FORMAT(1X,'PHI1 =',E12.5,2X,'PSI1 =',E12.5,2X,'HACH2 =',E12.5) 486 : 2000 FORMAT(1X,'RESULTATS DE L ECOULEMENT PLASTIQUE') 487 : 2001 FORMAT(1X,'SIGMA =',6(1X,1PE12.5)) 488 : 2002 FORMAT(1X,'DSIGMA =',6(1X,1PE12.5)) 489 : 2003 FORMAT(1X,'STOT =',6(1X,1PE12.5)) 490 : 2004 FORMAT(1X,'SIGFIN =',6(1X,1PE12.5)) 491 : 2005 FORMAT(1X,'DEFPLA =',6(1X,1PE12.5)) 492 : 2006 FORMAT(1X,'XLAMBD =',6(1X,1PE12.5)) 493 : 3000 FORMAT(1X,'XMAT =',/,3(6(1X,1PE12.5),/)) 494 : C 495 : 1900 RETURN 496 : END 497 : 498 : 499 : 500 : 501 : 502 : 503 :
© Cast3M 2003 - All rights reserved.
Disclaimer