Esope source files

   1 : C INTFIB    SOURCE    PV        11/03/07    21:17:04     6885           
   2 : C sub 
   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 INTFIB(XMAT,XSECT,DEPST,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),XSECT,DEPST 
  19 :         REAL*8 SIG0(3),SIGF,EINT 
  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    = VAR0(17) 
  87 :  
  88 :         XLAM    = VAR0(21) 
  89 :         EPTO    = VAR0(22) 
  90 :         INDITER = nint(VAR0(23))
  91 : C 
  92 : C-----AVANCE EN DEFORMATION 
  93 : C 
  94 :         EPTO = (EPTO + DEPST)*LCCOEF*(1.0D0-XLAM) 
  95 : C 
  96 : C-----PREDICTION ELASTIQUE 
  97 : C 
  98 :       SIGMA = 2.0D0*GCOEFF*EPTO*(1.0D0-ENDI) + 2.0D0*GCOEFF*(EPTO-EPSF)* 
  99 :      &      ENDI 
 100 :         XENER = GCOEFF*EPTO*EPTO 
 101 :         XENE0 = GCOEFF*EPS0CO*EPS0CO 
 102 :         SIGMP = 2.0D0*GCOEFF*ENDI*(EPTO-EPSF) 
 103 : C 
 104 : C-----CALCUL DE YR POUR L ACCROISSEMENT DE LA ROUILLE 
 105 : C 
 105 :         CALL APPFIB(TCI,XENRR,15) 
 106 : C 
 107 : C-----CALCUL DU CRITERE D ENDOMMAGEMENT 
 108 : C 
 109 :         CRIT1 = XENER - (XENE0 + ZVAR + XENRR) 
 110 :  
 111 :         IF (CRIT1.GT.0.0D0) THEN 
 112 :                 ENDI = 1.0D0-1.0D0/(1.0D0+AD*(XENER-(XENE0+XENRR))) 
 113 :                 ZVAR = XENER-XENE0-XENRR 
 114 :         END IF 
 115 : C 
 116 : C-----CALCUL DU CRITERE POUR LE FROTTEMENT 
 117 : C 
 118 :         CRIT2 = DABS(SIGMP-XVAR) 
 119 :  
 120 :         IF (CRIT2.GT.0.0D0) THEN 
 121 :  
 122 :                 CRIT20 = DABS(CRIT2) 
 123 :  
 124 :                 DO I1=1,ITERMAX 
 125 :  
 126 :                  CRIT2 = DABS(SIGMP-XVAR) 
 127 :  
 128 :                   IF (DABS(CRIT2/CRIT20).LE.1.0D-5) THEN 
 129 :                              GOTO 10                          
 130 :                   END IF 
 131 :  
 132 :                   DFSIGF =  DABS(SIGMP-XVAR)/(SIGMP-XVAR) 
 133 :                   DFXVAR = -DABS(SIGMP-XVAR)/(SIGMP-XVAR) 
 134 :                   DGSIGF =  DABS(SIGMP-XVAR)/(SIGMP-XVAR) 
 135 :                   DGXVAR = -DABS(SIGMP-XVAR)/(SIGMP-XVAR) + ACOEFF*XVAR 
 136 :  
 137 :                   DLAMPI = CRIT2/(ENDI*2.0D0*GCOEFF*DFSIGF*DGSIGF+ 
 138 :      &                     GAMCOE*DFXVAR*DGXVAR) 
 139 :          
 140 :                   XVAR   = XVAR  - GAMCOE*DLAMPI*DGXVAR 
 141 :                   SIGMP  = SIGMP - 2.0D0*GCOEFF*ENDI*DGSIGF*DLAMPI 
 142 :  
 143 :                 END DO 
 144 :  
 145 :  10   CONTINUE           
 146 :  
 147 :         ALPH = XVAR/GAMCOE 
 148 :         EPSF = EPTO - SIGMP/(2.0D0*ENDI*GCOEFF) 
 149 :  
 150 :         END IF 
 151 : C 
 152 : C-----ON REAJUSTE LES CONTRAINTES 
 153 : C 
 154 :       SIGMA = 2.0D0*GCOEFF*EPTO*(1.0D0-ENDI) + 2.0D0*GCOEFF*(EPTO-EPSF)* 
 155 :      &      ENDI 
 156 : C 
 157 : C-----LES SORTIES 
 158 : C 
 159 :       VARF(11) = ENDI 
 160 :       VARF(12) = ALPH 
 161 :       VARF(13) = EPSF 
 162 :       VARF(14) = XVAR 
 163 :       VARF(15) = ZVAR 
 164 :       VARF(16) = EINT 
 165 :       VARF(17) = SIGMA 
 166 :  
 167 :       VARF( 7) = EPTO 
 168 :       VARF( 8) = SIGMA 
 169 :  
 170 : C 
 171 : C-----LES CONTRAINTES EN SORTIE 
 172 : C 
 173 :         SIGF     = SIGMA 
 174 :  
 175 :         RETURN 
 176 :         END 
 177 :  
 178 :  
 179 :  
 180 :