C AVANX1 SOURCE CB215821 17/11/30 21:15:12 9639 SUBROUTINE AVANX1(TAU,SIG,EPSV,VAR,SIG1,EPSV1,VAR1, &SIGD,EPSVD,VARD,NSTRS,NVARI) C C---------------------------------------------------------------- C Objet: Mise a jour des variables internes C des deformations plastiques C des contraintes C---------------------------------------------------------------- C C---------------------------------------------------------------- C Entree: TAU pas de temps C SIG(NSTRS,NSTRS) contraintes au debut du pas C EPSV(NSTRS,NSTRS) deformations plastiques C au debut du pas C VAR(NVARI) variables internes au debut du pas C EPSVD(NSTRS,NSTRS) vitesse de deformations C visco-plastiques calculees au debut du pas C VARD(NVARI) vitesse des variables internes au debut du pas C SIGD(NSTRS) vitesse des contraintes C------------------------------------------------------------- C C------------------------------------------------------------- C Sortie: SIG1(NSTRS,NSTRS) etat de contraintes a la fin du pas C EPSV1(NSTRS,NSTRS) deformations plastiques la fin C du pas C VAR1(NVARI) variables internes a la fin du pas C-------------------------------------------------------------- C IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) logical dtlibr DIMENSION SIG(*),EPSV(*),VAR(*) DIMENSION SIG1(*),EPSV1(*),VAR1(*) DIMENSION EPSVD(*),VARD(*),SIGD(*) DIMENSION D1(3,3),Di(3),ROT(3,3) C write(6,*) 'Entrée dans AVANX1' C C---------------------- C CALCULATE VARIABLES C---------------------- 20 DO 50 I = 1,7 50 VAR1(I) = VAR(I)+TAU*VARD(I) XD11=VAR1(2) XD22=VAR1(3) XD33=VAR1(4) XD12=VAR1(5) XD13=VAR1(6) XD23=VAR1(7) CALL ZERO(D1,3,3) D1(1,1) = XD11 D1(2,2) = XD22 D1(3,3) = XD33 D1(1,2) = XD12 D1(1,3) = XD13 D1(2,3) = XD23 D1(2,1) = D1(1,2) D1(3,1) = D1(1,3) D1(3,2) = D1(2,3) C CALL JACOB3(D1,3,Di,ROT) C ADMAX = MAX(Di(1),Di(2),Di(3)) VAR1(8)=ADMAX C C------------------------------- C CALCULATE PLASTIC STRAINS C------------------------------- DO 60 I = 1,NSTRS 60 EPSV1(I) = EPSV(I)+TAU*EPSVD(I) C C ------------------------------- C CALCULATE STRESS C-------------------------------- DO 70 I=1,NSTRS 70 SIG1(I)=SIG(I)+TAU*SIGD(I) RETURN END