Esope source files

   1 : C DTCHPL    SOURCE    CHAT      05/01/12    22:57:15     5004           
   2 : c     Sous-Programme du module MISTRAL1
   3 : C     --------------------------------------------------------------------------   
   4 :       SUBROUTINE DTCHPL (KPLAST, SEQA,SEQ,RRA,RR,VEPEA,VEPSE,
   5 :      &                   DTA, DT, ICHPL)
   6 : C     --------------------------------------------------------------------------
   7 : C       Determination d'un eventuel changement d'etat :
   8 : C                   elastique (KPLAST=0) / plastique (KPLAST=1)
   9 : C     entre l'instant T et T+DT en fonction :
  10 : C     - des contraintes equivalentes SEQ a T et SEQA a T-DTA,
  11 : C     - des contraintes seuil RR a T et RRA a T-DTA,
  12 : C     ou :
  13 : C     - des vitesses de deformation equivalente VEPSE a T et VEPEA a T-DTA.
  14 : C     Le cas echeant : reduction du pas de temps DT prevu et ICHPL = 1,
  15 : C     sinon : ICHPL = 0.
  16 : C     --------------------------------------------------------------------------
  17 :       IMPLICIT INTEGER(I-N)
  18 :       IMPLICIT REAL*8 (A-H, O-Z)
  19 :         DTCHP = DT
  20 :         IF (KPLAST.EQ.0) THEN
  21 :          SMRA = SEQA-RRA
  22 :          SMR = SEQ-RR
  23 :          IF (SMR.GT.SMRA) THEN
  24 :           DTCHP = DTA*SMR/(SMRA-SMR)
  25 :          END IF
  26 :         ELSE
  27 :          IF (VEPSE.LT.VEPEA) THEN
  28 :           DTCHP = DTA*VEPSE/(VEPEA-VEPSE)         
  29 :          END IF
  30 :         END IF
  31 :         IF (DTCHP.LT.DT) THEN
  32 :          ICHPL = 1
  33 :          DT = DTCHP
  34 :         ELSE
  35 :          ICHPL = 0
  36 :         END IF
  37 :         RETURN
  38 :         END
  39 :       
  40 :