Télécharger barfra.eso

Retour à la liste

Numérotation des lignes :

barfra
  1. C BARFRA SOURCE CHAT 05/01/12 21:33:57 5004
  2. C BARFRA SOURCE DC 99/04/28
  3. SUBROUTINE BARFRA(XMAT,XCAR,DEPST,VAR0,SIGF,VARF,DEFP)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. DIMENSION XMAT(10),VAR0(4),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 ELEMENT DE BARRE
  12. C----------------------------------------------------------------------
  13. C Didier COMBESCURE - EMSI - Avril 1999
  14. C----------------------------------------------------------------------
  15. C
  16. C XMAT( 1) YOUNG
  17. C XMAT( 2) NU
  18. C XMAT( 3) RHO
  19. C XMAT( 4) ALPHA
  20. C XMAT( 5) FC
  21. C XMAT( 6) FC_R
  22. C XMAT( 7) STC
  23. C XMAT( 8) FT
  24. C XMAT( 9) FT_R
  25. C XMAT(10) STT
  26. C
  27. C SIGF STRESS
  28. C DEPST,DEFP
  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. DEPS = DEPST
  36. C
  37. GAMMA=XMAT(1)/2/(XMAT(2)+1)
  38. EYOUN = XMAT( 1)
  39. FC = (-1.)*XMAT(5)
  40. FC_R = (-1.)*XMAT(6)
  41. STC = (-1.)*XMAT(7)
  42. XEPSOC = FC/EYOUN
  43. FT = XMAT(8)
  44. FT_R = XMAT(9)
  45. STT = XMAT(10)
  46. XEPSOT = FT/EYOUN
  47. C
  48. XDAMC = VAR0(2)
  49. XDAMT = VAR0(3)
  50. C
  51. XEPSO = VAR0(1) + DEPST
  52. C
  53. C CALCUL DE LA VARIABLE D'ENDOMMAGEMENT
  54. C
  55. IF (XEPSO.GT. 0.D0) THEN
  56. C
  57. C CAS DE LA TRACTION
  58. C
  59. XDAMT0 = 1.0D0 - ( (XEPSOT/XEPSO)* ( (FT_R / FT) +
  60. & ((1.D0 - (FT_R / FT))*EXP((XEPSOT - XEPSO)/ STT) )) )
  61. C
  62. IF (XDAMT0. GT . XDAMT) THEN
  63. XDAMT = XDAMT0
  64. ENDIF
  65. C
  66. XDAM = XDAMT
  67. C
  68. ELSEIF (XEPSO.LT. 0.D0) THEN
  69. C
  70. C CAS DE LA COMPRESSION
  71. C
  72. XDAMC0 = 1.0D0 - ( (XEPSOC/XEPSO)* ( (FC_R /FC) +
  73. & ((1.D0 - (FC_R /FC))*EXP((XEPSOC - XEPSO)/STC) )) )
  74. C
  75. IF (XDAMC0. GT . XDAMC) THEN
  76. XDAMC = XDAMC0
  77. ENDIF
  78. C
  79. XDAM = XDAMC
  80. C
  81. ELSE
  82. XDAM= 0.D0
  83. ENDIF
  84. C
  85. C
  86. SIGF= EYOUN* (1.D0 - XDAM) * XEPSO * XCAR
  87. C
  88. VARF(1) = XEPSO
  89. VARF(2) = XDAMC
  90. VARF(3) = XDAMT
  91. VARF(4) = (1.D0 - XDAM) * EYOUN
  92. C
  93. DEFP = 0.D0
  94. C
  95. RETURN
  96. END
  97.  
  98.  
  99.  
  100.  

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