Télécharger lcpc.eso

Retour à la liste

Numérotation des lignes :

lcpc
  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
  55. * CALL ERREUR(7)
  56. * IF (FLU.GT.1D-15) THEN
  57. * STOP
  58. * ENDIF
  59. END
  60.  
  61.  
  62.  
  63.  
  64.  

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