Télécharger avanxx.eso

Retour à la liste

Numérotation des lignes :

  1. C AVANXX SOURCE CHAT 05/01/12 21:30:06 5004
  2. SUBROUTINE AVANXX(TAU,SIG,EPSV,VAR,SIG1,EPSV1,VAR1,
  3. &SIGD,EPSVD,VARD,NSTRS,NVARI)
  4. C
  5. C----------------------------------------------------------------
  6. C Objet: Mise a jour des variables internes p,r,D
  7. C des deformations plastiques
  8. C des contraintes
  9. C----------------------------------------------------------------
  10. C
  11. C----------------------------------------------------------------
  12. C Entree: TAU pas de temps
  13. C SIG(NSTRS,NSTRS) contraintes au debut du pas
  14. C EPSV(NSTRS,NSTRS) deformations plastiques
  15. C au debut du pas
  16. C VAR(NVARI) variables internes au debut du pas
  17. C EPSVD(NSTRS,NSTRS) vitesse de deformations
  18. C visco-plastiques calculees au debut du pas
  19. C VARD(NVARI) vitesse des variables internes au debut du pas
  20. C SIGD(NSTRS) vitesse des contraintes
  21. C-------------------------------------------------------------
  22. C
  23. C-------------------------------------------------------------
  24. C Sortie: SIG1(NSTRS,NSTRS) etat de contraintes a la fin du pas
  25. C EPSV1(NSTRS,NSTRS) deformations plastiques la fin
  26. C du pas
  27. C VAR1(NVARI) variables internes a la fin du pas
  28. C--------------------------------------------------------------
  29. C
  30. IMPLICIT INTEGER(I-N)
  31. IMPLICIT REAL*8(A-H,O-Z)
  32. logical dtlibr
  33. DIMENSION SIG(*),EPSV(*),VAR(*)
  34. DIMENSION SIG1(*),EPSV1(*),VAR1(*)
  35. DIMENSION EPSVD(*),VARD(*),SIGD(*)
  36. C
  37. C------------------------------------
  38. C CALCULATE VARIABIES : P , R , D
  39. C------------------------------------
  40. 20 DO 50 I = 1,NVARI
  41. 50 VAR1(I) = VAR(I)+TAU*VARD(I)
  42. IF (VAR1(3).LE.0.99) GOTO 60
  43. c modif antoine
  44. c TAU = (0.989D0-VAR(3))/VARD(3)
  45. c GOTO 20
  46. C
  47. C-------------------------------
  48. C CALCULATE PLASTIC STRAINS
  49. C-------------------------------
  50. 60 DO 70 I = 1,NSTRS
  51. 70 EPSV1(I) = EPSV(I)+TAU*EPSVD(I)
  52. C
  53. C -------------------------------
  54. C CALCULATE STRESS
  55. C--------------------------------
  56. DO 80 I=1,NSTRS
  57. 80 SIG1(I)=SIG(I)+TAU*SIGD(I)
  58. RETURN
  59. END
  60.  
  61.  
  62.  
  63.  
  64.  

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