Télécharger incree.eso

Retour à la liste

Numérotation des lignes :

incree
  1. C INCREE SOURCE CHAT 05/01/13 00:35:37 5004
  2. SUBROUTINE INCREE (TAU,SIG,P,QTLD,DECRA,DESWA,EPSVPT,EC0PT,ESW0PT,
  3. & NSTRS1,IVIEX,MFR1,IFOURB)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. C-----------------------------------------------------------------------
  7. C
  8. C DESCRIPTION FONCTIONNELLE :
  9. C -------------------------
  10. C
  11. C Lois 'VISCO_EXTERNE' : exploitation des sorties de CREEP pour le
  12. C calcul du tenseur des vitesses de deformation inelastique
  13. C
  14. C Entrees : TAU,SIG,P,QTLD,DECRA,DESWA
  15. C Sorties : EPSVPT,EC0PT,ESW0PT
  16. C
  17. C Parametres : NSTRS1,IVIEX,MFR1,IFOURB
  18. C Les indicateurs IVIEX,MFR1,IFOURB ne servent pas
  19. C pour l'instant mais seront utiles si extensions :
  20. C autres modeles que 'GENERAL', autres formulations EF,
  21. C autres options de calcul
  22. C
  23. C-----------------------------------------------------------------------
  24. C Arguments de l'interface
  25. C
  26. REAL*8 SIG(*), DECRA(*), DESWA(*), EPSVPT(*)
  27. C
  28. C Variables locales
  29. C
  30. C Tableaux de travail
  31. C
  32. REAL*8 SS(8)
  33. C
  34. C------------------- Debut du code executable --------------------------
  35. C
  36. C 1 - Initialisation du tenseur des vitesses de deformation inelastique
  37. C
  38. DO 10 I=1,NSTRS1
  39. EPSVPT(I)=0.0D0
  40. 10 CONTINUE
  41. C
  42. C 2 - Calcul des vitesses de deformation inelastique
  43. C (grandeurs scalaires)
  44. C
  45. EC0PT = DECRA(1)/TAU
  46. ESW0PT = DESWA(1)/TAU
  47. C
  48. C 3 - Premiere contribution au tenseur des vitesses de deformation
  49. C inelastique : 1/3 * ESW0PT * [Id]
  50. C
  51. CONS=0.333333333333333D0*ESW0PT
  52. DO 30 I=1,3
  53. EPSVPT(I)=EPSVPT(I)+CONS
  54. 30 CONTINUE
  55. C
  56. C 4 - Seconde contribution au tenseur des vitesses de deformation
  57. C inelastique : (3/2)/QTLD * [SS]
  58. C
  59. IF ( QTLD.LE.1.D-10 ) GOTO 1000
  60. C
  61. C Calcul du deviateur des contraintes
  62. C
  63. DO 40 I=1,NSTRS1
  64. A=0.0D0
  65. IF (I.LE.3) A=1.0D0
  66. SS(I)=SIG(I)-A*P
  67. 40 CONTINUE
  68. C
  69. CONS=(1.5D0*EC0PT)/QTLD
  70. DO 41 I=1,NSTRS1
  71. EPSVPT(I)=EPSVPT(I)+CONS*SS(I)
  72. 41 CONTINUE
  73. C
  74. 1000 CONTINUE
  75. RETURN
  76. END
  77.  
  78.  
  79.  

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