Télécharger la plaquette Cast3M


fr en

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