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  26/10/26    21:15:09     12682          
  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.       DO 50 I = 1,7
  43.         VAR1(I) = VAR(I)+TAU*VARD(I)
  44. 50    CONTINUE
  45.       XD11=VAR1(2)
  46.       XD22=VAR1(3)
  47.       XD33=VAR1(4)
  48.       XD12=VAR1(5)
  49.       XD13=VAR1(6)
  50.       XD23=VAR1(7)
  51.       CALL ZERO(D1,3,3)
  52.       D1(1,1) = XD11
  53.       D1(2,2) = XD22
  54.       D1(3,3) = XD33
  55.       D1(1,2) = XD12
  56.       D1(1,3) = XD13
  57.       D1(2,3) = XD23
  58.       D1(2,1) = D1(1,2)
  59.       D1(3,1) = D1(1,3)
  60.       D1(3,2) = D1(2,3)
  61. C
  62.       CALL JACOB3(D1,3,Di,ROT)
  63. C
  64.       ADMAX =  MAX(Di(1),Di(2),Di(3))
  65.       VAR1(8)=ADMAX
  66. C
  67. C-------------------------------
  68. C CALCULATE  PLASTIC  STRAINS
  69. C-------------------------------
  70.       DO  60 I = 1,NSTRS
  71.       EPSV1(I) = EPSV(I)+TAU*EPSVD(I)
  72. 60    CONTINUE
  73. C
  74. C -------------------------------
  75. C CALCULATE STRESS
  76. C--------------------------------
  77.       DO 70 I=1,NSTRS
  78.         SIG1(I)=SIG(I)+TAU*SIGD(I)
  79. 70    CONTINUE
  80.       RETURN
  81.       END
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  

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