Télécharger cetag.eso

Retour à la liste

Numérotation des lignes :

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

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