Télécharger dtchpl.eso

Retour à la liste

Numérotation des lignes :

  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.  
  41.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales