Télécharger fibfra.eso

Retour à la liste

Numérotation des lignes :

  1. C FIBFRA SOURCE CHAT 12/04/06 21:15:11 7348
  2. C FIBFRA SOURCE DC 98/01/26
  3. SUBROUTINE FIBFRA(XMAT,DEPS,SIG0,VAR0,SIGF,VARF)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. DIMENSION XMAT(*),DEPS(3),SIG0(3),VAR0(4),SIGF(3),VARF(4)
  7. C----------------------------------------------------------------------
  8. C
  9. C COMPORTEMENT DE FIBRE FRAGILE EN COMPRESSION ET TRACTION
  10. C MODELE PLASTIQUE POUR LE RENFORCEMENT TFC
  11. C----------------------------------------------------------------------
  12. C Didier COMBESCURE - EMSI - Janvier 1998
  13. C----------------------------------------------------------------------
  14. C
  15. C XMAT( 1) YOUNG
  16. C XMAT( 2) NU
  17. C XMAT( 3) RHO
  18. C XMAT( 4) ALPHA
  19. C XMAT( 5) FC
  20. C XMAT( 6) FC_R
  21. C XMAT( 7) STC
  22. C XMAT( 8) FT
  23. C XMAT( 9) FT_R
  24. C XMAT(10) STT
  25. C
  26. C DEPS(1) EPXX SIG0(1) SIGF(1) SMXX
  27. C DEPS(2) GAXY SIG0(2) SIGF(2) SMXY
  28. C DEPS(3) GAXZ SIG0(3) SIGF(3) SMXZ
  29. C
  30. C VAR0(1) VARF(1) EPSO ALLONGEMENT TOTAL
  31. C VAR0(2) VARF(2) DAMC ENDOMMAGEMENT EN COMP.
  32. C VAR0(3) VARF(3) DAMT ENDOMMAGEMENT EN TRACT.
  33. C
  34. C----------------------------------------------------------------------
  35. GAMMA=XMAT(1)/2.D0/(XMAT(2)+1.D0)
  36. EYOUN = XMAT( 1)
  37. FC = (-1.D0)*XMAT(5)
  38. FC_R = (-1.D0)*XMAT(6)
  39. STC = (-1.D0)*XMAT(7)
  40. XEPSOC = FC/EYOUN
  41. FT = XMAT(8)
  42. FT_R = XMAT(9)
  43. STT = XMAT(10)
  44. XEPSOT = FT/EYOUN
  45. C
  46. XDAMC = VAR0(2)
  47. XDAMT = VAR0(3)
  48. C
  49. XEPSO = VAR0(1) + DEPS(1)
  50. C
  51. C CALCUL DE LA VARIABLE D'ENDOMMAGEMENT
  52. C
  53. IF (XEPSO.GT. 0.D0) THEN
  54. C
  55. C CAS DE LA TRACTION
  56. C
  57. XDAMT0 = 1.0D0 - ( (XEPSOT/XEPSO)* ( (FT_R / FT) +
  58. & ((1.D0 - (FT_R / FT))*EXP((XEPSOT - XEPSO)/ STT) )) )
  59. C
  60. IF (XDAMT0. GT . XDAMT) THEN
  61. XDAMT = XDAMT0
  62. ENDIF
  63. C
  64. XDAM = XDAMT
  65. C
  66. ELSEIF (XEPSO.LT. 0.D0) THEN
  67. C
  68. C CAS DE LA COMPRESSION
  69. C
  70. XDAMC0 = 1.0D0 - ( (XEPSOC/XEPSO)* ( (FC_R /FC) +
  71. & ((1.D0 - (FC_R /FC))*EXP((XEPSOC - XEPSO)/STC) )) )
  72. C
  73. IF (XDAMC0. GT . XDAMC) THEN
  74. XDAMC = XDAMC0
  75. ENDIF
  76. C
  77. XDAM = XDAMC
  78. C
  79. ELSE
  80. XDAM= 0.D0
  81. ENDIF
  82. C
  83. C
  84. SIGF(1)= EYOUN* (1.D0 - XDAM) * XEPSO
  85. C
  86. SIGF(2) = SIG0(2) + GAMMA * DEPS(2)
  87. SIGF(3) = SIG0(3) + GAMMA * DEPS(3)
  88. C
  89. VARF(1) = XEPSO
  90. VARF(2) = XDAMC
  91. VARF(3) = XDAMT
  92. VARF(4) = (1.D0 - XDAM) * EYOUN
  93. C
  94. RETURN
  95. END
  96.  
  97.  
  98.  
  99.  
  100.  

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