Télécharger fibeta.eso

Retour à la liste

Numérotation des lignes :

  1. C FIBETA SOURCE CHAT 11/03/16 21:22:27 6902
  2. C FIBETA SOURCE AM 96/12/24 21:19:32 2448
  3. C======================================================================
  4. C FIBETA - D. COMBESCURE et P. PEGON - ELSA- 1996
  5. C======================================================================
  6. C
  7. C MODELE GLOBAL D'ETAGE
  8. C (Sur des elements de poutre TIMO - Fibre
  9. C (Contrainte de cisaillement/Cisail.)
  10. C
  11. SUBROUTINE FIBETA(XMAT,XCAR,SIG0,VAR0,SIGF,VARF,DEPS,
  12. . WRK2,NCURVP,NCURVN,KERRE)
  13. IMPLICIT INTEGER(I-N)
  14. IMPLICIT REAL*8(A-H,O-Z)
  15. C
  16. C=======================================================================
  17. C CETTE ROUTINE EST APPELE DANS ECO100
  18. C
  19. C
  20. C WRK0 = Segment materiaux
  21. C WRK1 = Segment contraintes
  22. C WRK2 = Segment courbe
  23. C NCURVP = Longueur courbe sens positif
  24. C NCURVN = Longueur courbe sens negatif
  25. C
  26. C======================================================================
  27. C XCAR = Caracteristique de la section
  28. C DEPS = Increment de deformation axiale
  29. C SIG0 = Contrainte initiale
  30. C VAR0 = Variables internes initiales
  31. C SIGF = Contrainte finale
  32. C VARF = Variables internes finales
  33. C DEFP = Deformation plastique
  34. C
  35. C=======================================================================
  36. C
  37. PARAMETER (XZER=0.D0,UN=1.D0,EPSILO=1.D-16)
  38. DIMENSION XMAT(*),XCAR(*),SIG0(*),VAR0(*),SIGF(*),VARF(*),DEPS(*)
  39. C
  40. C=======================================================================
  41. C VARIABLES ET SEGMENTS NECESSAIRES
  42. C=========================================================================
  43. -INC CCOPTIO
  44. C Segment de la courbe
  45. SEGMENT WRK2
  46. REAL*8 TRAC(LTRAC)
  47. ENDSEGMENT
  48. C
  49. C
  50. C Lecture materiau
  51. C
  52. YOUN = XMAT(1)
  53. XNU = XMAT(2)
  54. XDELAP = XMAT(5)
  55. XDELAN = (-1.D0)*XMAT(6)
  56. XDMAXP = XMAT(7)
  57. XDMAXN = XMAT(8)
  58. XBETA = XMAT(9)
  59. XALPH = XMAT(10)
  60. XTETA = XMAT(11)
  61. C
  62. C
  63. IF (IDIM.EQ.3) THEN
  64. SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(2)
  65. C
  66. XDDEP = DEPS(3)
  67. XFOR0 = SIG0(3)
  68. C
  69. ELSE
  70. C
  71. SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(1)
  72. C
  73. XDDEP = DEPS(2)
  74. XFOR0 = SIG0(2)
  75. ENDIF
  76. C
  77. C Appel de la loi
  78. C
  79. CALL ETAGEL(XDDEP,XFOR0,XFORF,VAR0,VARF,SECZ,
  80. & XDELAP,XDELAN,XDMAXP,XDMAXN,XBETA,XALPH,XTETA,
  81. & WRK2,NCURVP,NCURVN,KERRE)
  82. C
  83. C==================================================================
  84. C
  85. C On remplit les tableaux avant de sortir
  86. C
  87. C==================================================================
  88. IF (IDIM.EQ.3) THEN
  89. SIGF(1) = XZER
  90. SIGF(2) = XZER
  91. SIGF(3) = XFORF
  92. SIGF(4) = XZER
  93. SIGF(5) = XZER
  94. SIGF(6) = XZER
  95. ELSE
  96. SIGF(1) = XZER
  97. SIGF(2) = XFORF
  98. SIGF(3) = XZER
  99. ENDIF
  100. C
  101. C=====================================================
  102. C FIN DE LA ROUTINE DU MODELE
  103. C======================================================
  104. RETURN
  105.  
  106. END
  107.  
  108.  
  109.  
  110.  
  111.  

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