Télécharger mobflu.eso

Retour à la liste

Numérotation des lignes :

  1. C MOBFLU SOURCE CB215821 16/04/21 21:17:46 8920
  2. C MOBFLU SOURCE INSL 24/10/96
  3. SUBROUTINE MOBFLU(SIGR,STFL,D1,NSTRS,IFLU,RMOY,PHIF,TEMP0,
  4. & DTEMP,TEMP1,POAR)
  5. C
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. DIMENSION SIGR(NSTRS),STFL(NSTRS),D1(NSTRS,NSTRS)
  9. C----------------------------------------------------------------------
  10. *
  11. IF(IFLU.EQ.0) GOTO 10
  12. IF(IFLU.EQ.10.OR.IFLU.EQ.11.OR.IFLU.EQ.20.OR.IFLU.EQ.21) GOTO 20
  13. IF(IFLU.EQ.30.OR.IFLU.EQ.31) GOTO 30
  14. *
  15. WRITE(*,*) '!! ATTENTION CE MODELE DE FLUAGE N EST PAS DISPONIBLE'
  16. WRITE(*,*) ' (DANS MOBFLU) IFLU= ',IFLU
  17. STOP
  18. c+mdj
  19. C GOTO 100
  20. c+mdj
  21. *
  22. 10 CONTINUE
  23. *
  24. * CE MATERIAUX NE FLUE PAS
  25. *
  26. GOTO 100
  27. *
  28. 20 CONTINUE
  29. *
  30. * MODELE DE FLUAGE BETON : LOI BPEL91
  31. *
  32. * IFLU = 10 : FLUAGE (BPEL91_K(t)..) + BETON ELASTIQUE
  33. * IFLU = 11 : FLUAGE (BPEL91_K(t)..) + BETON ELASTOPLASTIQUE
  34. * IFLU = 20 : FLUAGE (BPEL91_K(t)=2) + BETON ELASTIQUE
  35. * IFLU = 21 : FLUAGE (BPEL91_K(t)=2) + BETON ELASTOPLASTIQUE
  36. *
  37. TEMPF= TEMP0 + DTEMP
  38. TEMPR= TEMPF - TEMP1
  39. IF(IFLU.EQ.10.OR.IFLU.EQ.11) THEN
  40. DKS = 1.D0/(1.D0+(20.D0*POAR))
  41. DKE = 0.4D0
  42. DKC = (120.D0-PHIF)/30.D0+.6667D0*(100.D0-PHIF)/(20.D0+RMOY)
  43. DKT1 = 100.D0/(100.D0+TEMP1)
  44. FACT = DKS*(DKE+DKC*DKT1)
  45. ELSE
  46. FACT = PHIF
  47. ENDIF
  48. *
  49. COEF0 = SQRT(TEMPR)
  50. COEF1 = 5.D0*SQRT(RMOY)
  51. *
  52. * CALCUL EN DEFORMATION DE FLUAGE TOTALE
  53. *
  54. COEFT = FACT*COEF0/(COEF0+COEF1)
  55. *
  56. CALL BST(D1,SIGR,NSTRS,NSTRS,STFL)
  57. DO 2 I=1,NSTRS
  58. STFL(I)=COEFT*STFL(I)
  59. 2 CONTINUE
  60. TEMP0 = TEMPF
  61. GOTO 100
  62. *
  63. 30 CONTINUE
  64. *
  65. * MODELE DE FLUAGE BETON :LOI1 ----> EPSfl = A1 * T**A2 * Sig/E
  66. *
  67. * IFLU = 31 : FLUAGE (LOI1) + BETON ELASTOPLASTIQUE
  68. * IFLU = 30 : FLUAGE (LOI1) + BETON ELASTIQUE
  69. *
  70. TEMPF = TEMP0 + DTEMP
  71. COEFT = RMOY*(TEMPF**PHIF)
  72. CALL BST(D1,SIGR,NSTRS,NSTRS,STFL)
  73. DO 3 I=1,NSTRS
  74. STFL(I)=COEFT*STFL(I)
  75. 3 CONTINUE
  76. TEMP0 = TEMPF
  77. GOTO 100
  78. C
  79. 100 CONTINUE
  80. C---------------------------------------------------------------------
  81. 1991 FORMAT(18(1X,E12.5))
  82. RETURN
  83. END
  84.  
  85.  
  86.  
  87.  
  88.  

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