1 : C INCVAR SOURCE ROUX 07/10/25 21:15:48 5933 2 : C INCREP SOURCE ROUX G-M 10/11/04 3 : SUBROUTINE INCVAR(SIG,VAR,EPSVPT,VARPT,XMAT,NSTRS0, 4 : & MFR,NVARI,NCOMAT,INPLAS) 5 : IMPLICIT REAL*8(A-H,O-Z) 6 : DIMENSION SIG(*),VAR(*),EPSVPT(*),VARPT(*),XMAT(*) 7 : DIMENSION XX(6),ST(3,3),D(3),XP(3,3) 8 : 9 : C*********************************************** 10 : C CALCUL MECANIQUE 11 : C*********************************************** 12 : 13 : C------------------------------ 14 : C Initialisation des variables 15 : C------------------------------ 16 : IF (INPLAS.EQ.137.OR.INPLAS.EQ.139) THEN 17 : IF (INPLAS.EQ.137) THEN 18 : C1=XMAT(11) 19 : gam1=XMAT(12) 20 : C2=XMAT(13) 21 : ELSE 22 : C1=XMAT(10) 23 : gam1=XMAT(11) 24 : C2=XMAT(12) 25 : ENDIF 26 : ENDIF 27 : 28 : C----------------- 29 : C Calcul de J2(ss) 30 : C----------------- 31 : C calcul de ss, déviateur des contraintes 32 : TRACE =(SIG(1)+SIG(2)+SIG(3))/3.0D0 33 : DO 70 I=1,NSTRS0 34 : A = 0.0D0 35 : IF (I.LE.3) A=1.0D0 36 : IF (INPLAS.EQ.136.OR.INPLAS.EQ.138) THEN 37 : XX(I) = SIG(I)-A*TRACE 38 : ENDIF 39 : IF (INPLAS.EQ.137.OR.INPLAS.EQ.139) THEN 40 : XXX = SIG(I)-A*TRACE-VAR(I)*C1 41 : XX(I) = XXX-VAR(I+NSTRS0)*C2 42 : ENDIF 43 : 70 CONTINUE 44 : AJ2 = PROCON (XX,XX,NSTRS0) 45 : AJ2 = SQRT(1.5D0*AJ2) 46 : 47 : C------------------ 48 : C test d'élasticité 49 : C------------------ 50 : IF (INPLAS.EQ.136.OR.INPLAS.EQ.138) THEN 51 : RR=VAR(1)*XMAT(6)*XMAT(7) 52 : ENDIF 53 : IF (INPLAS.EQ.137.OR.INPLAS.EQ.139) THEN 54 : RR=VAR((2*NSTRS0)+1)*XMAT(6)*XMAT(7) 55 : ENDIF 56 : SY=XMAT(5) 57 : EE=AJ2-SY-RR 58 : IF (EE.GT.0D0) THEN 59 : TESTE=0.0D0 60 : ELSE 61 : TESTE=1.0D0 62 : ENDIF 63 : 64 : C------------------------------------------------- 65 : C vitesse de déformation plastique cumulée ppt=rpt 66 : C------------------------------------------------- 67 : IF (TESTE.EQ.1.0D0) THEN 68 : RPT=0.0D0 69 : ELSE 70 : IF (INPLAS.EQ.136.OR.INPLAS.EQ.137) THEN 71 : RPT1=EE/XMAT(8) 72 : RPT2=(exp(RPT1)-exp(-1.D0*RPT1))/2.d0 73 : RPT=XMAT(10)*(RPT2)**XMAT(9) 74 : ENDIF 75 : IF (INPLAS.EQ.138.OR.INPLAS.EQ.139) THEN 76 : RPT=(EE/XMAT(8))**XMAT(9) 77 : ENDIF 78 : ENDIF 79 : IF (INPLAS.EQ.136.OR.INPLAS.EQ.138) THEN 80 : VARPT(2)=RPT 81 : ENDIF 82 : IF (INPLAS.EQ.137.OR.INPLAS.EQ.139) THEN 83 : VARPT((2*NSTRS0)+2)=RPT 84 : ENDIF 85 : 86 : C---------------------------------- 87 : C vitesses de déformation plastique 88 : C---------------------------------- 89 : IF (TESTE.EQ.1.0D0) THEN 90 : DO 80 I=1,NSTRS0 91 : EPSVPT(I)=0.0D0 92 : 80 CONTINUE 93 : ELSE 94 : DO 90 I=1,NSTRS0 95 : EPSVPT(I)=(1.5D0*RPT*XX(I))/AJ2 96 : 90 CONTINUE 97 : ENDIF 98 : 99 : C---------------------------------------------- 100 : C vitesse d'évolution de l'écrouissage isotrope 101 : C---------------------------------------------- 102 : IF (INPLAS.EQ.136.OR.INPLAS.EQ.138) THEN 103 : Dr=(1-XMAT(6)*VAR(1))*RPT 104 : VARPT(1)=Dr 105 : ENDIF 106 : IF (INPLAS.EQ.137.OR.INPLAS.EQ.139) THEN 107 : Dr=(1-XMAT(6)*VAR((2*NSTRS0)+1))*RPT 108 : VARPT((2*NSTRS0)+1)=Dr 109 : ENDIF 110 : 111 : C------------------------------------------------- 112 : C vitesse d'évolution de l'écrouissage cinematique 113 : C------------------------------------------------- 114 : IF (INPLAS.EQ.137.OR.INPLAS.EQ.139) THEN 115 : C======NON-LINEAIRE 116 : DO 96 I=1,NSTRS0 117 : VARPT(I)=EPSVPT(I)-gam1*RPT*VAR(I) 118 : 96 CONTINUE 119 : C======NON-LINEAIRE 120 : DO 97 I=1,NSTRS0 121 : VARPT(I+NSTRS0)=EPSVPT(I) 122 : 97 CONTINUE 123 : ENDIF 124 : 125 : 126 : 127 : 128 : RETURN 129 : END 130 :
© Cast3M 2003 - All rights reserved.
Disclaimer