Télécharger bpel.eso

Retour à la liste

Numérotation des lignes :

  1. C BPEL SOURCE PV 07/11/23 21:15:32 5978
  2. SUBROUTINE BPEL(RM,RH,FC28,T0,T,FLU)
  3. *
  4. * Calcul de la fonction de fluage selon l'eurocode 2
  5. *
  6. *
  7. * entree
  8. * RM rayon moyen de la piece [cm]
  9. * RH pourcentage d'humidite (0.<roh<100.)
  10. * FC28 resistance moyenne en compression a 28 jours [MPa] type de ciment
  11. * T0 temps de chargement [j]
  12. * T temps actuel [j]
  13. *
  14. * sortie
  15. * FLU J(to,t) [Pa-1]
  16. *
  17. *
  18. IMPLICIT INTEGER(I-N)
  19. IMPLICIT REAL*8(A-H,O-Z)
  20. real*8 j
  21. *
  22. *
  23. J=T0
  24. FCJ=J/(4.76D0+0.83D0*J)*FC28
  25. IF (FCJ.GT.(1.1D0*FC28)) THEN
  26. FCJ=1.1D0*FC28
  27. ENDIF
  28. IF (FCJ.GT.40.D0) THEN
  29. FCJ=J/(1.4D0+0.95D0*J)*FC28
  30. ENDIF
  31. IF (FCJ.GT.(1.1D0*FC28)) THEN
  32. FCJ=1.1D0*FC28
  33. ENDIF
  34. ECT0=11000.D0*FCJ**(1.D0/3.D0)
  35. EC28=11000.D0*FC28**(1.D0/3.D0)
  36. AKT0=100.D0/(100.D0+T0)
  37. AKC=(120.D0-RH)/30.D0+2.D0*(100.D0-RH)/
  38. & (3.D0*(20.D0+RM))
  39. AKFL=0.4D0+AKT0*AKC
  40. FTT0=SQRT(T-T0)/(SQRT(T-T0)+5.D0*SQRT(RM))
  41. PHI=AKFL*FTT0
  42. FLU=(1.D0/ECT0+PHI/EC28)/1.D6
  43. * write(6,*) ' rm rh fc28 t t0 ', rm,th,fc28,t ,t0
  44. * write(6,*) ' flu ' , flu
  45. END
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  

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