Esope source files

   1 : C LCPC      SOURCE    CHAT      05/01/13    01:13:35     5004           
   2 :         SUBROUTINE LCPC(H0,RHOS,E1AN,T0,T,FLU)
   3 : *
   4 : *  Calcul de la fonction de fluage selon la 
   5 : *  methode LCPC-Granger adaptee a FLAMANVILLE 2
   6 : *
   7 : *
   8 : * entree
   9 : *        H0    epaisseur du mur [mm], permettant de calculer FACK2
  10 : *              facteur d'echelle traduisant le passage d'une eprouvette
  11 : *              de diametre phi a un mur d'epaisseur H0.
  12 : *              FACK2 = ((2 * phi)/(3 * H0))**2
  13 : *              Valeurs donnees pour phi = 160 mm.
  14 : *        RHOS  taux d'armatures passives permettant de calculer KSIGP
  15 : *              KSIGP = (1+20*RhoS)**(-1)
  16 : *        E1AN  Module d'Young mesure en laboratoire a 1 an d'age.
  17 : *        T0   temps de chargement [j]
  18 : *        T    temps actuel [j]
  19 : *
  20 : * sortie
  21 : *        FLU  J(to,t) [Pa-1]
  22 : *
  23 : *        
  24 :       IMPLICIT INTEGER(I-N)
  25 :         IMPLICIT REAL*8(A-H,O-Z)
  26 : 
  27 :         REAL*8 KSA,KSIGP
  28 : *
  29 :         EPSFP = 65.63D-6
  30 :         AFP = 0.4969242D0
  31 :         BFPJ = 21.97750113D0
  32 :         EPSFT = 262.6D-6
  33 :         AFT = 0.30071D0
  34 :         BFTJ = 17.39335116D0
  35 :         PHI = 160.D0
  36 : *
  37 : *        WRITE(6,*) 'H0 = ',H0
  38 : *        WRITE(6,*) 'RHOS = ',RHOS
  39 : *        WRITE(6,*) 'T0 = ',T0
  40 : *        WRITE(6,*) 'T = ',T
  41 : *
  42 :         KSA = 1.D0/(1.D0 + (16.D0*1.071D-2))
  43 :         KSIGP = 1.D0/(1.D0 + (20.D0 * RHOS))
  44 :         FACK2 = ((2.D0 * PHI)/(3.D0 * H0))**(2.)       
  45 :         GCT028 = (0.1D0 + (28.D0**0.2))/(0.1D0 + (T0**0.2))
  46 :         FPTMT0 = EPSFP*((T-T0)**AFP)/(((T-T0)**AFP) + BFPJ)
  47 :         FPK2 = EPSFP*((FACK2*(T-T0))**AFP)/(((FACK2*(T-T0))**AFP)
  48 :      &                                                    + BFPJ)
  49 :         FTK2 = EPSFT*((FACK2*(T-T0))**AFT)/(((FACK2*(T-T0))**AFT)
  50 :      &                                                    + BFTJ)
  51 : *
  52 :         FLU = (1.D0/E1AN + (KSA * KSIGP * GCT028 * 
  53 :      &                      (FPTMT0 + FTK2 - FPK2)))/1.D6
  54 : *        WRITE(6,*) 'FLU = ',FLU
  54 : *        CALL ERREUR(7)
  55 : *        IF (FLU.GT.1D-15) THEN
  56 : *          STOP
  57 : *       ENDIF
  58 :         END
  59 :  
  60 :  
  61 :  
  62 :