C ADVA22 SOURCE CHAT 05/01/12 21:18:38 5004 SUBROUTINE ADVA22(TAU,SIG,EPSV,VAR,SIG1,EPSV1,VAR1,SIGPT,EPSVPT, & VARPT,XMAT,NSTRS,NVARI) C ----------------------------------------------------------------- IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) -INC CCOPTIO DIMENSION SIG(*),EPSV(*),VAR(*),SIG1(*),EPSV1(*),VAR1(*) DIMENSION SIGPT(*),EPSVPT(*),VARPT(*),XMAT(*) C C ICI : le modele est obligatoirement CERAMIQUE C---------------------------------------------------------------- C THIS ROUTINE CALCULATES THE LATEST VALUES OF SIG,EPSV,VAR, C FROM THE INITIAL VALUES OF SIG,EPSV,VAR AND DERIVATIVES : C SIGPT,EPSVPT,VARPT C---------------------------------------------------------------- C T1 = XMAT(1)/(1.D0-XMAT(2)*XMAT(2)) T2 = T1 * XMAT(2) T3 = XMAT(1)/(1.D0+XMAT(2)) NS = NSTRS/2 C C________________________________________________________________ C STRAIN COMPONENTS C________________________________________________________________ DO 10 I =1,NSTRS EPSV1(I) = EPSV(I) + TAU*EPSVPT(I) 10 CONTINUE C C________________________________________________________________ C STRESS COMPONENTS C________________________________________________________________ IF(NS.EQ.2) THEN SIG1(1) = SIG(1) + TAU*(SIGPT(1)- T1*EPSVPT(1)-T2*EPSVPT(2)) SIG1(2) = SIG(2) + TAU*(SIGPT(2)- T1*EPSVPT(2)-T2*EPSVPT(1)) SIG1(3) = SIG(3) + TAU*(SIGPT(3)- T1*EPSVPT(3)-T2*EPSVPT(4)) SIG1(4) = SIG(4) + TAU*(SIGPT(4)- T1*EPSVPT(4)-T2*EPSVPT(3)) IF(IFOUR.EQ.-2) THEN SIG1(2)=0.D0 SIG1(4)=0.D0 ENDIF ELSE SIG1(1) = SIG(1) + TAU*(SIGPT(1)- T1*EPSVPT(1)-T2*EPSVPT(2)) SIG1(2) = SIG(2) + TAU*(SIGPT(2)- T1*EPSVPT(2)-T2*EPSVPT(1)) SIG1(3) = SIG(3) + TAU*(SIGPT(3)- T3*EPSVPT(3)) SIG1(4) = SIG(4) + TAU*(SIGPT(4)- T1*EPSVPT(4)-T2*EPSVPT(5)) SIG1(5) = SIG(5) + TAU*(SIGPT(5)- T1*EPSVPT(5)-T2*EPSVPT(4)) SIG1(6) = SIG(6) + TAU*(SIGPT(6)- T3*EPSVPT(6)) IF(IFOUR.EQ.-2) THEN SIG1(2)=0.D0 SIG1(5)=0.D0 ENDIF ENDIF C________________________________________________________________ C INTERNAL VARIABLES C________________________________________________________________ C La première variable interne qui était la derniere dans advan2 correspond à EPSE VAR1(1) = VAR(1)+TAU*VARPT(1) C La dernière variable interne correspond à EPSE fluage VAR1(NVARI) = VAR1(NVARI)+TAU*VARPT(NVARI) C l'avant dernière variable interne correspond à l'endommagement C on n'y touche pas ici VAR1(NVARI-1) = VAR(NVARI-1) CCC IF ((IFOUR.EQ.2.AND.MFR.EQ.3).OR.IFOUR.EQ.-2) THEN CCC IBID = 11 CCC ELSE IF(IFOUR.EQ.-3.OR.IFOUR.EQ.-1 CCC 1 .OR.IFOUR.EQ.0) THEN CCC IBID = 14 CCC ELSE IF(IFOUR.EQ.2.AND.MFR.EQ.1) THEN CCC IBID = 19 CCC ENDIF CCC Eloi : correction de IBID CCC de plus : ici on a obligatoirement MFR = 3 IBID=NVARI-(2*NSTRS)-5 C de 2 à IBID on a les variables d'ottosen qu'on ne touche pas C On ne change pas les variables de OTTOSEN DO 20 I=2,IBID VAR1(I)=VAR(I) 20 CONTINUE C de IBID + 1 à IBID + 2*NSTRS+3 DO 30 I=(IBID+1),(IBID+2*NSTRS+3) VAR1(I)=VARPT(I) 30 CONTINUE RETURN END C pour mémoire voici ce qu'il y avait dans advan2 ************************************************** C C IF((INPLAS.GE.19.AND.INPLAS.LE.23).OR.INPLAS.EQ.61) THEN C DO 30 I=1,2*NSTRS+3 C VAR1(I)=VARPT(I) C 30 CONTINUE C en fait dans IDVAR1 on a 2*NSTRS+4 = NVARI C il s'agit donc de la variable EPSE qu'on met à jour C DO 40 I=2*NSTRS+4,NVARI C VAR1(I)=VAR(I)+TAU*VARPT(I) C 40 CONTINUE ***********************************************