Télécharger fibpar.eso

Retour à la liste

Numérotation des lignes :

  1. C FIBPAR SOURCE CHAT 12/04/06 21:15:12 7348
  2. C FIBFRA SOURCE DC 98/01/26
  3. SUBROUTINE FIBPAR(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(2),SIGF(3),VARF(2)
  7. C----------------------------------------------------------------------
  8. C
  9. C COMPORTEMENT DE FIBRE
  10. C MODELE BETON PLASTIQUE DU BAEL
  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
  18. C XMAT( 5) FY
  19. C XMAT( 6) HH
  20. C
  21. C DEPS(1) EPXX SIG0(1) SIGF(1) SMXX
  22. C DEPS(2) GAXY SIG0(2) SIGF(2) SMXY
  23. C DEPS(3) GAXZ SIG0(3) SIGF(3) SMXZ
  24. C
  25. C VAR0(1) VARF(1) EPSE DEFORMATION PLASTIQUE
  26. C
  27. C----------------------------------------------------------------------
  28. GAMMA=XMAT(1)/2.D0/(XMAT(2)+1.D0)
  29. EYOUN = XMAT( 1)
  30. FYP = XMAT(5)
  31. FYN = (-1.D0)*XMAT(5)
  32. HHH = XMAT(6)
  33. FYPE = FYP + ( HHH * VAR0(1))
  34. FYNE = FYN + ( HHH * VAR0(1))
  35. C
  36. XEPSE = VAR0(1)
  37. C
  38. XSIGT = SIG0(1) + EYOUN * DEPS(1)
  39. C
  40. IF (XSIGT . GT . FYPE) THEN
  41. C
  42. XDEP = ((EYOUN/(EYOUN+HHH))*DEPS(1))
  43. & + ((SIG0(1) - FYPE)/(EYOUN+HHH))
  44. XSIGF = SIG0(1) + (EYOUN*(DEPS(1) - XDEP))
  45. XTAN = (EYOUN*HHH)/(EYOUN+HHH)
  46. C
  47. ELSEIF (XSIGT . LT . FYNE) THEN
  48. C
  49. XDEP = ((EYOUN/(EYOUN+HHH))*DEPS(1))
  50. & + ((SIG0(1) - FYNE)/(EYOUN+HHH))
  51. XSIGF = SIG0(1) + (EYOUN*(DEPS(1) - XDEP))
  52. XTAN = (EYOUN*HHH)/(EYOUN+HHH)
  53. C
  54. ELSE
  55. C
  56. XSIGF = XSIGT
  57. XDEP = 0.D0
  58. XTAN = EYOUN
  59. C
  60. ENDIF
  61. C
  62. C
  63. SIGF(1) = XSIGF
  64. C
  65. SIGF(2) = SIG0(2) + GAMMA * DEPS(2)
  66. SIGF(3) = SIG0(3) + GAMMA * DEPS(3)
  67. C
  68. VARF(1) = XEPSE + XDEP
  69. VARF(2) = XTAN
  70. C
  71. RETURN
  72. END
  73.  
  74.  
  75.  
  76.  
  77.  
  78.  

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