Télécharger avanx1.eso

Retour à la liste

Numérotation des lignes :

  1. C AVANX1 SOURCE AM 08/12/19 21:15:03 6237
  2. SUBROUTINE AVANX1(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
  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. DIMENSION D1(3,3),Di(3),ROT(3,3)
  37. C write(6,*) 'Entrée dans AVANX1'
  38. C
  39. C----------------------
  40. C CALCULATE VARIABLES
  41. C----------------------
  42. 20 DO 50 I = 1,7
  43. 50 VAR1(I) = VAR(I)+TAU*VARD(I)
  44. XD11=VAR1(2)
  45. XD22=VAR1(3)
  46. XD33=VAR1(4)
  47. XD12=VAR1(5)
  48. XD13=VAR1(6)
  49. XD23=VAR1(7)
  50. CALL ZERO(D1,3,3)
  51. D1(1,1) = XD11
  52. D1(2,2) = XD22
  53. D1(3,3) = XD33
  54. D1(1,2) = XD12
  55. D1(1,3) = XD13
  56. D1(2,3) = XD23
  57. D1(2,1) = D1(1,2)
  58. D1(3,1) = D1(1,3)
  59. D1(3,2) = D1(2,3)
  60. C
  61. CALL JACOB3(D1,3,Di,ROT)
  62. C
  63. ADMAX = DMAX1(Di(1),Di(2),Di(3))
  64. VAR1(8)=ADMAX
  65. C
  66. C-------------------------------
  67. C CALCULATE PLASTIC STRAINS
  68. C-------------------------------
  69. DO 60 I = 1,NSTRS
  70. 60 EPSV1(I) = EPSV(I)+TAU*EPSVD(I)
  71. C
  72. C -------------------------------
  73. C CALCULATE STRESS
  74. C--------------------------------
  75. DO 70 I=1,NSTRS
  76. 70 SIG1(I)=SIG(I)+TAU*SIGD(I)
  77. RETURN
  78. END
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  

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