Télécharger la plaquette Cast3M


fr en

Télécharger initte.eso

Retour à la liste

Numérotation des lignes : 

initte

  1. C INITTE    SOURCE    CHAT      05/01/13    00:38:20     5004
  2.  
  3.       SUBROUTINE INITTE(NYUNG,YUNG,NXNU,XNU,NRHO,RHO,NALFA,ALFA,NEPSD,
  4.      1EPSD,NDC,DC,NEPSR,EPSR,NTRAC,TRAC,NC,NPNTER,X,T,XMATT,NXMATT,
  5.      2TRUC,NCOURB)
  6. C
  7. C
  8. C-----------------------------------------------------------------
  9. C  OBJET: initialisation a la temperature T des caracteristiques
  10. C         mecaniques du materiau dans le cas de la plasticite
  11. C         avec ecrouissage et endommagement isotrope,en regime
  12. C         anisotherme
  13. C-----------------------------------------------------------------
  14. C
  15. C-----------------------------------------------------------------
  16. C ENTREE:
  17. C -------
  18. C          NYUNG = 2*nombre de points de la courbe T -->E(T)
  19. C          YUNG  = tableau contenant les valeurs connues du
  20. C                  module d'young E a des temperatures precisees Ti
  21. C          NXNU  = 2*nombre de points de la courbe T -->NU(T)
  22. C          XNU   = tableau contenant les valeurs connues du coeff.
  23. C                  de Poisson NU a des temperatures precisees Ti
  24. C          NRHO  = 2*nombre de points de la courbe T -->RHO(T)
  25. C          RHO   = tableau contenant les valeurs connues de la masse
  26. C                  volumique a des temperatures precisees Ti
  27. C          NALFA = 2*nombre de points de la courbe T -->ALFA(T)
  28. C          ALFA  = tableau contenant les valeurs connues du coeff.
  29. C                  de dilatation thermique ALFA a des temperatures
  30. C                  precisees Ti.(EPSth=ALFA*DT)
  31. C          NEPSD = 2*nombre de points de la courbe T -->EPSD(T)
  32. C          EPSD  = tableau contenant les valeurs connues du seuil
  33. C                  d'endomt EPSD a des temperatures precisees Ti
  34. C          NDC   = 2*nombre de points de la courbe T -->DC(T)
  35. C          DC    = tableau contenant les valeurs connues de l'endomt
  36. C                  critique DC a des temperatures precisees Ti
  37. C          NEPSR = 2*nombre de points de la courbe T -->EPSR(T)
  38. C          EPSR  = tableau contenant les valeurs connues de la def.
  39. C                  a rupture EPSR a des temperatures precisees Ti
  40. C          NTRAC = dimension de TRAC
  41. C          TRAC  = tableau contenant les courbes de traction connues
  42. C                  a des temperatures precisees Ti
  43. C          NC    = nombres de courbes de traction
  44. C          NPNTER = tableau contenant le nombre de points pour
  45. C                  chaque courbe de traction
  46. C          X     = abscisse pour laquelle on veut calculer l'ordonnee de
  47. C                  la nouvelle courbe de traction a la temperature T
  48. C          T     = nouvelle temperature pour laquelle on recalcule
  49. C                  les caracteristiques mecaniques du materiau
  50. C------------------------------------------------------------------
  51. C
  52. C------------------------------------------------------------------
  53. C SORTIES:
  54. C --------
  55. C  XMATT(NXMATT) = tableau contenant les nouvelles caracteristiques
  56. C                  mecaniques du materiau
  57. C------------------------------------------------------------------
  58.       IMPLICIT INTEGER(I-N)
  59.       IMPLICIT REAL*8 (A-H,O-Z)
  60. C
  61.       DIMENSION YUNG(*),XNU(*),RHO(*),ALFA(*),EPSD(*),DC(*)
  62.       DIMENSION EPSR(*),TRAC(*),NPNTER(*),XMATT(*),TRUC(*)
  63. C
  64. C------------------------------------------------------------------
  65. C  si les caracteristiques mecaniques du materiau sont constantes
  66. C  YUNG,XNU,ALFA,EPSD,DC,EPSR seront des tableaux de dimension 1
  67. C  TRAC ne change pas;on calcule l'ordonnee du point de la courbe
  68. C  de traction d'abscisse X;
  69. C  sinon on calcule les nouvelles caracteristiques du materiau a T
  70. C  et on les stocke dans XMATT.
  71. C------------------------------------------------------------------
  72.           CALL DERTRA(NYUNG,YUNG,T,ET,EPRIM,TINF,TSUP)
  73.           XMATT(1)=ET
  74. C
  75.           CALL DERTRA(NXNU,XNU,T,XNUT,XNURIM,TINF,TSUP)
  76.           XMATT(2)=XNUT
  77. C
  78.           CALL DERTRA(NRHO,RHO,T,RHOT,RHORIM,TINF,TSUP)
  79.           XMATT(3)=RHOT
  80. C
  81.           CALL DERTRA(NALFA,ALFA,T,ALFAT,ALFRIM,TINF,TSUP)
  82.           XMATT(4)=ALFAT
  83. C
  84.           CALL DERTRA(NEPSD,EPSD,T,EPSDT,EPSDIM,TINF,TSUP)
  85.           XMATT(6)=EPSDT
  86. C
  87.           CALL DERTRA(NDC,DC,T,DCT,DCPRIM,TINF,TSUP)
  88.           XMATT(7)=DCT
  89. C
  90.           CALL DERTRA(NEPSR,EPSR,T,EPSRT,EPSRIM,TINF,TSUP)
  91.           XMATT(8)=EPSRT
  92. C
  93.           CALL VALPAR(TRAC,NPNTER,NTRAC,NC,X,Y,YPRIM,T,TINF,TSUP,
  94.      1                                               TRUC,NCOURB)
  95.           XMATT(5)=Y
  96. C---------------------------------------------------------------
  97. C    Y= ordonnee du point d'abscisse X de la nouvelle courbe de
  98. C       traction a T
  99. C---------------------------------------------------------------
  100. C
  101.       RETURN
  102.       END
  103. **********************************************************************
  104.  

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