1 : C AVANP SOURCE PV 09/12/21 21:15:00 6589 2 : C 3 : C PLG93Debut 4 : C 5 : SUBROUTINE AVANP(TAU,SIG,EPSV,VAR,SIG1,EPSV1,VAR1,DSPT, 6 : & EPSVPT,VARPT,XMAT,NSTRS,NVARI,IFOUR,NCOMAT) 7 : IMPLICIT INTEGER(I-N) 8 : IMPLICIT REAL*8(A-H,O-Z) 9 : DIMENSION SIG(*),EPSV(*),VAR(*),SIG1(*),EPSV1(*),VAR1(*) 10 : DIMENSION DSPT(*),EPSVPT(*),VARPT(*),XMAT(*) 11 : C ____________________________________________________________________ 12 : C 13 : C THIS ROUTINE UPDATES THE STRAINS ,STRESSES & INTERNAL VARIABLES 14 : C Version modifiee de la procedure advan1 pour le modele 15 : C viscoplastique poudre. 16 : C ____________________________________________________________________ 17 : C 18 : C connaissant l'etat c'est a dire SIG, EPSV et VAR ainsi que les 19 : C derivees EPSVPT et DSPT on calcule l'etat SIG1, EPSV1 et VAR1 20 : C au temps T + TAU 21 : C ____________________________________________________________________ 22 : C 23 : C STEP-1 : UPDATING STRAINS 24 : C ---------------------------- 25 : YOU = XMAT(1) 26 : XNU = XMAT(2) 27 : DO 10 I=1,NSTRS 28 : EPSV1(I)=EPSV(I)+TAU*EPSVPT(I) 29 : 10 CONTINUE 30 : C 31 : C STEP-2 : UPDATING STESSES 32 : C ---------------------------- 33 : IF(IFOUR.EQ.-2) THEN 34 : SIG(3)=0.0D0 35 : X2MU= YOU/(1.0D0+XNU) 36 : XCO = X2MU/(1.0D0-XNU) 37 : SIG1(1)=SIG(1)+TAU*(DSPT(1)-XCO*(EPSVPT(1)+XNU*EPSVPT(2))) 38 : SIG1(2)=SIG(2)+TAU*(DSPT(2)-XCO*(EPSVPT(2)+XNU*EPSVPT(1))) 39 : SIG1(4)=SIG(4)+TAU*(DSPT(4)-X2MU*EPSVPT(4)) 40 : SIG1(3)=0.0D0 41 : ELSEIF (IFOUR.EQ.2.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-1.OR. 42 : & IFOUR.EQ.1.OR.IFOUR.EQ.-3) THEN 43 : TRACE=EPSVPT(1)+EPSVPT(2)+EPSVPT(3) 44 : X2MU = YOU / (1.0D0+XNU) 45 : ALMB = X2MU*XNU/(1.0D0-2.0D0*XNU) 46 : DO 20 I=1,NSTRS,1 47 : A = 0.0D0 48 : IF (I.LE.3) A=1.0D0 49 : SIG1(I) = SIG(I) + TAU*(DSPT(I)-ALMB*A*TRACE-X2MU*EPSVPT(I)) 50 : 20 CONTINUE 51 : C= Modes de calcul 1D contraintes planes suivant z (DYCZ et GYCZ) 52 : ELSE IF (IFOUR.EQ.4.OR.IFOUR.EQ.8) THEN 53 : SIG(3) = 0.D0 54 : XCO = YOU/(1.D0-XNU*XNU) 55 : SIG1(1) = SIG(1) + TAU*(DSPT(1)-XCO*(EPSVPT(1)+XNU*EPSVPT(2))) 56 : SIG1(2) = SIG(2) + TAU*(DSPT(2)-XCO*(EPSVPT(2)+XNU*EPSVPT(1))) 57 : SIG1(3) = 0.D0 58 : C= Modes de calcul 1D contraintes planes suivant y (CYDZ et CYGZ) et 59 : C= 1D axisymetrique contraintes axiales nulles (AXCZ) 60 : ELSE IF (IFOUR.EQ.5.OR.IFOUR.EQ.10.OR.IFOUR.EQ.13) THEN 61 : SIG(2) = 0.D0 62 : XCO = YOU/(1.D0-XNU*XNU) 63 : SIG1(1) = SIG(1) + TAU*(DSPT(1)-XCO*(EPSVPT(1)+XNU*EPSVPT(3))) 64 : SIG1(2) = 0.D0 65 : SIG1(3) = SIG(3) + TAU*(DSPT(3)-XCO*(EPSVPT(3)+XNU*EPSVPT(1))) 66 : C= Modes de calcul 1D contraintes planes suivant y et z (CYCZ) 67 : ELSE IF (IFOUR.EQ.6) THEN 68 : SIG(2) = 0.D0 69 : SIG(3) = 0.D0 70 : SIG1(1) = SIG(1) + TAU*(DSPT(1)-YOU*EPSVPT(1)) 71 : SIG1(2) = 0.D0 72 : SIG1(3) = 0.D0 73 : C= Autres modes de calcul 1D deformations planes (DYDZ GYDZ DYGZ GYGZ) 74 : C= 1D axisymetrique (AXDZ AXGZ) et 1D spherique 75 : ELSE IF (IFOUR.EQ. 3.OR.IFOUR.EQ. 7.OR.IFOUR.EQ. 9.OR.IFOUR.EQ.11 76 : . .OR.IFOUR.EQ.12.OR.IFOUR.EQ.14.OR.IFOUR.EQ.15) THEN 77 : X2MU = YOU/(1.D0+XNU) 78 : TRACE = EPSVPT(1)+EPSVPT(2)+EPSVPT(3) 79 : XCO = X2MU*XNU*TRACE/(1.D0-XNU) 80 : SIG1(1) = SIG(1) + TAU*(DSPT(1)-X2MU*EPSVPT(1)-XCO) 81 : SIG1(2) = SIG(2) + TAU*(DSPT(2)-X2MU*EPSVPT(2)-XCO) 82 : SIG1(3) = SIG(3) + TAU*(DSPT(3)-X2MU*EPSVPT(3)-XCO) 83 : ENDIF 84 : C 85 : C STEP-3 : UPDATING THE INTERNAL VARIABLES 86 : C------------------------------------------ 87 : DO 30 I=1,NSTRS 88 : VAR1(1+I)=EPSV1(I) 89 : 30 CONTINUE 90 : C 91 : C majoration de la densite relative : RHOR=RHOR-TraceEPSV*RHOR 92 : C 93 : TRE = EPSVPT(1)+EPSVPT(2)+EPSVPT(3) 94 : VAR1(NVARI) = VAR(NVARI) * EXP(-TRE*TAU) 95 : IF (VAR1(NVARI).GT.1D0) VAR1(NVARI) = 1.D0 96 : IF (VAR1(NVARI).ne.var1(nvari)) VAR1(NVARI) = 1.D0 97 : C 98 : C majoration de EPSE : 99 : C 100 : VAR1(1) = VAR(1)+TAU*VARPT(1) 101 : C 102 : C _____________________________________________________________________ 103 : RETURN 104 : END 105 : C 106 : 107 : 108 : 109 : 110 : 111 : 112 : 113 :
© Cast3M 2003 - All rights reserved.
Disclaimer