Esope source files

   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------------------------------------------------------------------
  71 :           CALL DERTRA(NYUNG,YUNG,T,ET,EPRIM,TINF,TSUP)
  72 :           XMATT(1)=ET
  73 : C
  73 :           CALL DERTRA(NXNU,XNU,T,XNUT,XNURIM,TINF,TSUP)
  74 :           XMATT(2)=XNUT
  75 : C
  75 :           CALL DERTRA(NRHO,RHO,T,RHOT,RHORIM,TINF,TSUP)
  76 :           XMATT(3)=RHOT
  77 : C
  77 :           CALL DERTRA(NALFA,ALFA,T,ALFAT,ALFRIM,TINF,TSUP)
  78 :           XMATT(4)=ALFAT
  79 : C
  79 :           CALL DERTRA(NEPSD,EPSD,T,EPSDT,EPSDIM,TINF,TSUP)
  80 :           XMATT(6)=EPSDT
  81 : C
  81 :           CALL DERTRA(NDC,DC,T,DCT,DCPRIM,TINF,TSUP)
  82 :           XMATT(7)=DCT
  83 : C
  83 :           CALL DERTRA(NEPSR,EPSR,T,EPSRT,EPSRIM,TINF,TSUP)
  84 :           XMATT(8)=EPSRT
  85 : C
  85 :           CALL VALPAR(TRAC,NPNTER,NTRAC,NC,X,Y,YPRIM,T,TINF,TSUP,
  86 :      1                                               TRUC,NCOURB)
  87 :           XMATT(5)=Y
  88 : C---------------------------------------------------------------
  89 : C    Y= ordonnee du point d'abscisse X de la nouvelle courbe de
  90 : C       traction a T
  91 : C---------------------------------------------------------------
  92 : C
  93 :       RETURN
  94 :       END
  95 : **********************************************************************