Esope source files

   1 : C ENPLAS    SOURCE    LJASON    07/09/24    21:15:08     5866           
   2 :        SUBROUTINE ENPLAS(XMAT,NMATT,VAR0,VARF,NVARI,SIG0,
   3 :      & SIGF,DEPST,NSTRS,KERRE, ISTEP)
   4 : 
   5 : *******************************************************
   6 : *     Variables                                       *
   7 : *                                                     *
   8 : *     XMAT : Variables matériau                       *
   9 : *                                                     *
  10 : *     VAR0 : Variables internes à l'état initial      *
  11 : *     VARF : Variables internes après incrément       *
  12 : *                                                     *
  13 : *     SIG0 : Contrainte à l'incrément précédent       *
  14 : *     SIGF : Contrainte à l'incrément suivant         *
  15 : *                                                     *
  16 : *     DEPST : Incrément de déplacement total imposé   *
  17 : *******************************************************
  18 : 
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8(A-H,O-Z)
  21 : 
  22 : *****************************
  23 : * Déclaration des variables *
  24 : *****************************
  25 : 
  26 :       INTEGER NMATT, NVARI, NSTRS, KERRE, ISTEP
  27 : 
  28 : * Entrees
  29 :       REAL*8 XMAT(NMATT), VAR0(NVARI), SIG0(NSTRS)
  30 :       REAL*8 DEPST(NSTRS) , SIG02(6)
  31 : 
  32 : * Sorties
  33 :       REAL*8 SIGF(NSTRS), VARF(NVARI)
  34 : 
  35 : * Matrice tangente
  36 :       REAL*8 DSDE(6,6)
  37 : 
  38 : * Variables materiau
  39 :       REAL*8 E,NU
  40 :       REAL*8 AC, AT, BC, BT, EPSD0, RC, RT
  41 :       REAL*8 P, AH, BH, CH, GAMMA, ALPHA
  42 :       REAL*8 A, K0, BETA
  43 : 
  44 : * Transport des deformations
  45 :       DO I = 4,NSTRS
  46 :         DEPST(I) = DEPST(I)/2.d0
  47 :       END DO
  48 : 
  49 :       DO I=1,NSTRS
  50 :          SIG02(I) = SIG0(I)
  51 :       END DO
  52 : 
  53 : ********************
  54 : * Partie plastique *
  55 : ********************
  56 : 
  56 :          CALL PPLAST(XMAT,NMATT,VAR0,VARF,NVARI,SIG0,SIGF,DEPST,
  57 :      &   NSTRS,KERRE)
  58 : 
  59 :          DO I = 1,NSTRS
  60 :            SIG0(I) = SIGF(I)
  61 :          END DO
  62 : 
  63 : ************************
  64 : * Partie endommagement *
  65 : ************************
  66 : 
  67 : * Appel a la routine MAZARS modifiee
  68 : 
  68 :          CALL NMAZAR(XMAT,NMATT,VAR0,VARF,NVARI,SIG0,SIGF,
  69 :      &   DEPST,NSTRS,ISTEP)
  70 : 
  71 :        DO I = 4,NSTRS
  72 :          DEPST(I)=2.d0*DEPST(I)
  73 :        END DO        
  74 : 
  75 :        IF (ISTEP.EQ.1) THEN
  76 :          DO I=1,NSTRS
  77 :             SIGF(I)=SIG02(I)
  78 :          END DO
  79 :        END IF
  80 : 
  81 :        IF (KERRE.NE.0.d0) THEN
  82 :          WRITE(*,*) KERRE
  83 :        END IF
  84 : 
  85 :        RETURN
  86 :        END
  87 :