Télécharger la plaquette Cast3M


fr en

Télécharger grain.eso

Retour à la liste

Numérotation des lignes : 

grain

  1. C GRAIN     SOURCE    CHAT      05/01/13    00:20:57     5004
  2. C GRAIN     SOURCE    LLB
  3.       SUBROUTINE GRAIN(TAU,EVP1,EVP2,SIG1,SIG2,VAR1,VAR2,XMAT,
  4.      &     NSTRS,NVARI,KERRE)
  5. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  6. C
  7. C       CALCUL DE LA TAILLE DE GRAIN A L'AIDE D'UN SHEMA DE TYPE
  8. C    RUNGE KUTTA 2.1. LE MODELE UTILISE EST CELUI DE BESSON.
  9. C ________________________________________________________________________
  10. C
  11. C   EN ENTREE :
  12. C
  13. C    DE1     : LE TAUX DE DEFAUTS AU TEMPS i
  14. C    GR1     : LA TAILLE DE GRAIN AU TEMPS i
  15. C    TAU     : LE PAS DE TEMPS
  16. C    EVP1    : LA VITESSE DE DEFORMATION VISCOPLASTIQUE AU TEMPS i
  17. C    EVP2    : LA VITESSE DE DEFORAMTIOM VISCOPLATIQUE AU TEMPS i+TAU/2
  18. C    SIG1    : LA CONTRAINTE AU TEMPS i
  19. C    SIG2    : LA CONTRAINTE AU TEMPS i+TAU
  20. C    XMAT    : LE TABLEAU DES CARACTERISTIQUES MATERIELLES AU TEMPS i+TAU
  21. C    NSTRS   : LA DIMENSION DES VECTEURS EVP1...SIG1...
  22. C
  23. C   EN SORTIE :
  24. C
  25. C    DE2     : LE TAUX DE DEFAUT AU TEMPS i+TAU
  26. C    GR2     : LA TAILLE DE GRAIN AU TEMPS i+TAU
  27. C
  28. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  29. C
  30.       IMPLICIT INTEGER(I-N)
  31.       IMPLICIT REAL*8(A-H,O-Z)
  32. -INC PPARAM
  33. -INC CCOPTIO
  34.       DIMENSION EVP1(*),EVP2(*),SIG1(*),VAR1(*),SIG2(*),VAR2(*)
  35.       DIMENSION XMAT(*),SIGS2(6)
  36. C
  37. C--------------RAPPEL DES DONNES MATERIAU CALCULEE EN DEBUT DE PAS DE TEMPS
  38. C--------------ATTENTION, DANS L'ETAT ACTUEL, CE SONT LES SEULES DONNEES
  39. C--------------DISPONIBLES. DANS L'AVENIR, IL FAUDRA REACTUALISER LES K0, K1
  40. C--------------K2, K3 EN COURS DE PAS DE TEMPS.
  41. C
  42.       GR=VAR1(NVARI-2)
  43.       DE=VAR1(NVARI-1)
  44.       G0 = XMAT(20)
  45.       D0 = XMAT(21)
  46.       RM = XMAT(22)
  47.       RK0 = XMAT(23)
  48.       RK1 = XMAT(24)
  49.       RK2 = XMAT(25)
  50.       RK3 = XMAT(26)
  51. C
  52. C--------------ADIMENSIONNALISATION DES VARIABLES INTERNES
  53. C
  54.       DE1=DE/D0
  55.       GR1=GR/G0
  56. C
  57. C--------------CALCUL DE L'ENERGIE DE DEFORMATION AU TEMPS i
  58. C
  59.       ENERI=PROCON(EVP1,SIG1,NSTRS)
  60. C
  61. C--------------CALCUL DE L'ENERGIE DE DEFORMATION AU TEMPS i+TAU/2
  62. C
  63.       DO 10 I=1,NSTRS
  64.        SIGS2(I)= 0.5*(SIG1(I)+SIG2(I))
  65. 10    CONTINUE
  66.       ENERIS2=PROCON(EVP2,SIGS2,NSTRS)
  67. C
  68. C--------------CALCUL DE D AU TEMPS i+TAU/2
  69. C
  70.       D1S2=DE1+((TAU/2.D0)*(((RK1/D0)*ENERI)-(RK2*DE1)))
  71. C
  72. C--------------CALCUL G AU TEMPS i+TAU/2
  73. C
  74.       G1S2=GR1+(0.5D0*TAU*(((RK0/((G0*GR1)**RM))+(RK3*D0*DE1))/G0))
  75. C
  76. C--------------CALCUL DE G AU TEMPS i+1
  77. C
  78.       GR2=GR1+TAU*(((RK0/((G0*G1S2)**RM))+RK3*D0*D1S2)/G0)
  79. C
  80. C--------------CALCUL DE D AU TEMPS i+1
  81. C
  82.       DE2=DE1+(TAU*(((RK1/D0)*ENERIS2)-(RK2*D1S2)))
  83. C
  84. C   REMPLISSAGE DU TABLEAU DE VARIABLES INTERNES
  85. C
  86.       VAR2(NVARI-2)=GR2*G0
  87.       VAR2(NVARI-1)=DE2*D0
  88. C
  89.       RETURN
  90.       END
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  

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