Télécharger etage.eso

Retour à la liste

Numérotation des lignes :

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

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