Télécharger la plaquette Cast3M


fr en

Télécharger avanx1.eso

Retour à la liste

Numérotation des lignes : 

avanx1

  1. C AVANX1    SOURCE    CB215821  17/11/30    21:15:12     9639           
  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 =  MAX(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.  
  85.  

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