Télécharger flure2.eso

Retour à la liste

Numérotation des lignes :

  1. C FLURE2 SOURCE CHAT 05/01/13 00:04:43 5004
  2. SUBROUTINE FLURE2(H0,RHOS,E1AN,T0,TABLE,IDIM,REL)
  3. *
  4. * inversion de la fonction de fluage de LCPC - Granger
  5. *
  6. * entree
  7. * H0,RHOS,E1AN,T0 : parametres de la fonction
  8. * de fluage de LCPC-Granger
  9. * TABLE : tableau des temps de calcul
  10. * IDIM : dimension de TABLE
  11. *
  12. * sortie
  13. * REL : tableau des valeurs de la fonction de relaxation
  14. * calculee aux temps de TABLE
  15. *
  16. IMPLICIT INTEGER(I-N)
  17. IMPLICIT REAL*8(A-H,O-Z)
  18. *
  19. REAL*8 TABLE(0:1000),REL(0:1000)
  20. *
  21. * initialisation de REL a 0
  22. *
  23. DO 5 I=0,1000
  24. REL(I)=0.
  25. 5 CONTINUE
  26. *
  27. * calcul de REL(0)=rel(t0,t0)
  28. *
  29. CALL LCPC(H0,RHOS,E1AN,T0,T0,FLU)
  30. REL(0)=1./FLU
  31. T1=TABLE(1)
  32. *
  33. * calcul de REL(K)=rel(t0,table(k))
  34. *
  35. DO 10 K=1,IDIM
  36. TK=TABLE(K)
  37. TKM=TABLE(K-1)
  38. CALL LCPC(H0,RHOS,E1AN,T0,TK,FLU0)
  39. CALL LCPC(H0,RHOS,E1AN,T1,TK,FLU1)
  40. SOM=2.+(FLU1-FLU0)*REL(0)
  41. IF (K.GE.2) THEN
  42. DO 20 I=1,K-1
  43. TIM=TABLE(I-1)
  44. TIP=TABLE(I+1)
  45. CALL LCPC(H0,RHOS,E1AN,TIP,TK,FLUP)
  46. CALL LCPC(H0,RHOS,E1AN,TIM,TK,FLUM)
  47. SOM=SOM+(FLUP-FLUM)*REL(I)
  48. 20 CONTINUE
  49. ENDIF
  50. CALL LCPC(H0,RHOS,E1AN,TK,TK,FLUK)
  51. CALL LCPC(H0,RHOS,E1AN,TKM,TK,FLUKM)
  52. REL(K)=SOM/(FLUK+FLUKM)
  53. 10 CONTINUE
  54. END
  55.  
  56.  
  57.  

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