1 : C INTFIC SOURCE PV 11/03/07 21:17:05 6885 2 : C sub intfibc 3 : C====&===1=========2=========3=========4=========5=========6=========7== 4 : C Commentaires : Subroutine permettant de mettre en oeuvre le 5 : C modele RICINT pour une formulation à fibre 6 : C Auteurs : B. Richard (doctorant) 7 : C Date : Février 2008 8 : C====&===1=========2=========3=========4=========5=========6=========7== 9 : SUBROUTINE INTFIC(XMAT,DEPS,SIG0,VAR0,SIGF,VARF) 10 : C 11 : C-----DECLARATION GENERALE 12 : C 13 : IMPLICIT INTEGER(I-N) 14 : IMPLICIT REAL*8(A-H,O-Z) 15 : C 16 : C-----DECLARATION DES VARIABLES 17 : C 18 : REAL*8 XMAT(18),VAR0(25),VARF(25),DEPST 19 : REAL*8 SIG0(3),SIGF(3),EINT,DEPS(3) 20 : REAL*8 GCOEFF,AD,ACOEFF,GAMCOE,LCCOEF,EPS0CO,TCI 21 : REAL*8 ENDI,ALPH,EPSF,XVAR,ZVAR,XLAM,EPTO,SIGMA 22 : INTEGER INDITER,I1,ITERMAX 23 : REAL*8 XENER,SIGMP,CRIT1,XENRR,CRIT2,DLAMPI,CRIT20 24 : REAL*8 DFSIGF,DFXVAR,DGSIGF,DGXVAR,XENE0 25 : C 26 : C-----LES PARAMETRES 27 : C 28 : ITERMAX = 100000 29 : C 30 : C-----NOMENCLATURE DES PARAMETRES MATERIAUX (XMAT) ET VARIABLES INTERNES 31 : C (VAR0 ET VARF) 32 : C 33 : 34 : C --------------------------- 35 : C | INTERFACE A/B CORRODEE | 36 : C --------------------------- 37 : C XMAT(11) GCOEFF 38 : C XMAT(12) AD 39 : C XMAT(13) ACOEFF 40 : C XMAT(14) GAMCOE 41 : C XMAT(15) LCCOEF 42 : C XMAT(16) EPS0CO 43 : C XMAT(17) TCI 44 : C XMAT(18) CALA 45 : C 46 : C --------------------------- 47 : C | INTERFACE A/B CORRODEE | 48 : C --------------------------- 49 : C VAR0(11) VARF(11) ENDI 50 : C VAR0(12) VARF(12) ALPH 51 : C VAR0(13) VARF(13) EPSF 52 : C VAR0(14) VARF(14) XVAR 53 : C VAR0(15) VARF(15) ZVAR 54 : C VAR0(16) VARF(16) EINT 55 : C VAR0(17) VARF(17) SIGF 56 : C VAR0(18) VARF(18) 57 : C VAR0(19) VARF(19) 58 : C VAR0(20) VARF(20) 59 : C 60 : C --------------------------- 61 : C | AUTRES | 62 : C --------------------------- 63 : C VAR0(21) VARF(21) XLAM 64 : C VAR0(22) VARF(22) EPTO 65 : C VAR0(23) VARF(23) INDITER 66 : 67 : C 68 : C-----LES PARAMETRES MATERIAUX 69 : C 70 : GCOEFF = XMAT(11) 71 : AD = XMAT(12) 72 : ACOEFF = XMAT(13) 73 : GAMCOE = XMAT(14) 74 : LCCOEF = XMAT(15) 75 : EPS0CO = XMAT(16) 76 : TCI = XMAT(17) 77 : C 78 : C-----LES VARIABLES INTERNES 79 : C 80 : ENDI = VAR0(11) 81 : ALPH = VAR0(12) 82 : EPSF = VAR0(13) 83 : XVAR = VAR0(14) 84 : ZVAR = VAR0(15) 85 : EINT = VAR0(16) 86 : SIGF(1) = VAR0(17) 87 : 88 : XLAM = 0.0D0 89 : EPTO = VAR0(22) 90 : INDITER = nint(VAR0(23)) 91 : C 92 : C-----AVANCE EN DEFORMATION 93 : C 94 : DEPST= DEPS(1) 95 : EPTO = (EPTO + DEPST)*LCCOEF*(1.0D0-XLAM) 96 : C 97 : C-----PREDICTION ELASTIQUE 98 : C 99 : SIGMA = 2.0D0*GCOEFF*EPTO*(1.0D0-ENDI) + 2.0D0*GCOEFF*(EPTO-EPSF)* 100 : & ENDI 101 : XENER = GCOEFF*EPTO*EPTO 102 : XENE0 = GCOEFF*EPS0CO*EPS0CO 103 : SIGMP = 2.0D0*GCOEFF*ENDI*(EPTO-EPSF) 104 : C 105 : C-----CALCUL DE YR POUR L ACCROISSEMENT DE LA ROUILLE 106 : C 106 : CALL APPFIB(TCI,XENRR,15) 107 : C 108 : C-----CALCUL DU CRITERE D ENDOMMAGEMENT 109 : C 110 : CRIT1 = XENER - (XENE0 + ZVAR + XENRR) 111 : 112 : IF (CRIT1.GT.0.0D0) THEN 113 : ENDI = 1.0D0-1.0D0/(1.0D0+AD*(XENER-(XENE0+XENRR))) 114 : ZVAR = XENER-XENE0-XENRR 115 : END IF 116 : C 117 : C-----CALCUL DU CRITERE POUR LE FROTTEMENT 118 : C 119 : CRIT2 = DABS(SIGMP-XVAR) 120 : 121 : IF (CRIT2.GT.0.0D0) THEN 122 : 123 : CRIT20 = DABS(CRIT2) 124 : 125 : DO I1=1,ITERMAX 126 : 127 : CRIT2 = DABS(SIGMP-XVAR) 128 : 129 : IF (DABS(CRIT2/CRIT20).LE.1.0D-5) THEN 130 : GOTO 10 131 : END IF 132 : 133 : DFSIGF = DABS(SIGMP-XVAR)/(SIGMP-XVAR) 134 : DFXVAR = -DABS(SIGMP-XVAR)/(SIGMP-XVAR) 135 : DGSIGF = DABS(SIGMP-XVAR)/(SIGMP-XVAR) 136 : DGXVAR = -DABS(SIGMP-XVAR)/(SIGMP-XVAR) + ACOEFF*XVAR 137 : 138 : DLAMPI = CRIT2/(ENDI*2.0D0*GCOEFF*DFSIGF*DGSIGF+ 139 : & GAMCOE*DFXVAR*DGXVAR) 140 : 141 : XVAR = XVAR - GAMCOE*DLAMPI*DGXVAR 142 : SIGMP = SIGMP - 2.0D0*GCOEFF*ENDI*DGSIGF*DLAMPI 143 : 144 : END DO 145 : 146 : 10 CONTINUE 147 : 148 : ALPH = XVAR/GAMCOE 149 : EPSF = EPTO - SIGMP/(2.0D0*ENDI*GCOEFF) 150 : 151 : END IF 152 : C 153 : C-----ON REAJUSTE LES CONTRAINTES 154 : C 155 : SIGMA = 2.0D0*GCOEFF*EPTO*(1.0D0-ENDI) + 2.0D0*GCOEFF*(EPTO-EPSF)* 156 : & ENDI 157 : C 158 : C-----LES SORTIES 159 : C 160 : VARF(11) = ENDI 161 : VARF(12) = ALPH 162 : VARF(13) = EPSF 163 : VARF(14) = XVAR 164 : VARF(15) = ZVAR 165 : VARF(16) = EINT 166 : VARF(17) = SIGMA 167 : 168 : VARF( 7) = EPTO 169 : VARF( 8) = SIGMA 170 : VARF(22) = EPTO 171 : VARF(23) = INDITER 172 : 173 : C 174 : C-----LES CONTRAINTES EN SORTIE 175 : C 176 : SIGF(1) = SIGMA 177 : 178 : RETURN 179 : END 180 : 181 : 182 : 183 :
© Cast3M 2003 - All rights reserved.
Disclaimer