Télécharger grain.eso

Retour à la liste

Numérotation des lignes :

  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 CCOPTIO
  33. DIMENSION EVP1(*),EVP2(*),SIG1(*),VAR1(*),SIG2(*),VAR2(*)
  34. DIMENSION XMAT(*),SIGS2(6)
  35. C
  36. C--------------RAPPEL DES DONNES MATERIAU CALCULEE EN DEBUT DE PAS DE TEMPS
  37. C--------------ATTENTION, DANS L'ETAT ACTUEL, CE SONT LES SEULES DONNEES
  38. C--------------DISPONIBLES. DANS L'AVENIR, IL FAUDRA REACTUALISER LES K0, K1
  39. C--------------K2, K3 EN COURS DE PAS DE TEMPS.
  40. C
  41. GR=VAR1(NVARI-2)
  42. DE=VAR1(NVARI-1)
  43. G0 = XMAT(20)
  44. D0 = XMAT(21)
  45. RM = XMAT(22)
  46. RK0 = XMAT(23)
  47. RK1 = XMAT(24)
  48. RK2 = XMAT(25)
  49. RK3 = XMAT(26)
  50. C
  51. C--------------ADIMENSIONNALISATION DES VARIABLES INTERNES
  52. C
  53. DE1=DE/D0
  54. GR1=GR/G0
  55. C
  56. C--------------CALCUL DE L'ENERGIE DE DEFORMATION AU TEMPS i
  57. C
  58. ENERI=PROCON(EVP1,SIG1,NSTRS)
  59. C
  60. C--------------CALCUL DE L'ENERGIE DE DEFORMATION AU TEMPS i+TAU/2
  61. C
  62. DO 10 I=1,NSTRS
  63. SIGS2(I)= 0.5*(SIG1(I)+SIG2(I))
  64. 10 CONTINUE
  65. ENERIS2=PROCON(EVP2,SIGS2,NSTRS)
  66. C
  67. C--------------CALCUL DE D AU TEMPS i+TAU/2
  68. C
  69. D1S2=DE1+((TAU/2.D0)*(((RK1/D0)*ENERI)-(RK2*DE1)))
  70. C
  71. C--------------CALCUL G AU TEMPS i+TAU/2
  72. C
  73. G1S2=GR1+(0.5D0*TAU*(((RK0/((G0*GR1)**RM))+(RK3*D0*DE1))/G0))
  74. C
  75. C--------------CALCUL DE G AU TEMPS i+1
  76. C
  77. GR2=GR1+TAU*(((RK0/((G0*G1S2)**RM))+RK3*D0*D1S2)/G0)
  78. C
  79. C--------------CALCUL DE D AU TEMPS i+1
  80. C
  81. DE2=DE1+(TAU*(((RK1/D0)*ENERIS2)-(RK2*D1S2)))
  82. C
  83. C REMPLISSAGE DU TABLEAU DE VARIABLES INTERNES
  84. C
  85. VAR2(NVARI-2)=GR2*G0
  86. VAR2(NVARI-1)=DE2*D0
  87. C
  88. RETURN
  89. END
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  

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