Télécharger breva1.eso

Retour à la liste

Numérotation des lignes :

  1. C BREVA1 SOURCE AM 11/12/19 21:15:00 7227
  2. SUBROUTINE BREVA1(SIG0,DEPST,VAR0,XMAT,XCAR,SIGF,VARF,DEFP)
  3. C-----------------------------------------------------------------------
  4. C
  5. C DRIVER OF THE BILINEAR PLASTIC BENDING MODEL
  6. C
  7. C INPUT:
  8. C
  9. C SIG0 INITIAL GENERALIZED STRESS
  10. C DEPST INCREMENT OF GENERALIZED STRAIN
  11. C XMAT MATERIAL CHARACTERISTICS
  12. C XCAR GEOMETRIC CHARACTERISTICS
  13. C VAR0 INITIAL INTERNAL VARIABLE
  14. C
  15. C OUTPUT:
  16. C
  17. C SIGF FINAL GENERALIZED STRESS
  18. C VARF FINAL INTERNAL VARIABLE
  19. C DEFP INCREMENT OF PLASTIC DEFORMATION
  20. C
  21. C GENERALIZED STRESS (.=0 OR F):
  22. C
  23. C MOMY = SIG.(5)
  24. C
  25. C MATERIAL CHARACTERISTICS:
  26. C
  27. C E = XMAT(1) ELASTIC YOUNG MODULUS
  28. C NU = XMAT(2) POISSON RATIO
  29. C ETAN = XMAT(5) PLASTIC MODULUS
  30. C YIEF = XMAT(6) YIELDING FORCE
  31. C
  32. C GEOMETRIC CHARACTERISTICS:
  33. C
  34. C INRX = XCAR(1) X-INERTIAL TORSION
  35. C INRY = XCAR(2) Y-INERTIAL BENDING <-----
  36. C INRZ = XCAR(3) Z-INERTIAL BENDING
  37. C SECT = XCAR(4) SECTION
  38. C SRDY = XCAR(5) TRANSVERSE SHEAR'S REDUCED SECTION
  39. C SRDZ = XCAR(6) TRANSVERSE SHEAR'S REDUCED SECTION
  40. C
  41. C INTERNAL VARIABLE (.=0 OR F):
  42. C
  43. C EPSE = VAR.(1) CUMULATED PLASTIC AXIAL STRAIN
  44. C EITA = VAR.(2) TANGENT X-AXIAL STIFFNESS
  45. C MOEY = VAR.(3) ELASTIC AXIAL FORCE
  46. C EPSX = VAR.(4) AXIAL STRAIN
  47. C
  48. IMPLICIT INTEGER(I-N)
  49. IMPLICIT REAL*8(A-H,O-Z)
  50. -INC CCOPTIO
  51. C
  52. DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*),
  53. . SIGF(*),VARF(*) ,DEFP(*)
  54. C
  55. C-----------------------------------------------------------------------
  56. C
  57. C ELASTIC STRESS INCREMENT FOR ALL STRESS COMPONENT EXCEPT EFFX
  58. C
  59. CISA=0.5D0*XMAT(1)/(1.D0+XMAT(2))
  60. SIGF(2)=SIG0(2) + CISA*XCAR(5)*DEPST(2)
  61. SIGF(3)=SIG0(3) + CISA*XCAR(6)*DEPST(3)
  62. SIGF(4)=SIG0(4) + CISA*XCAR(1)*DEPST(4)
  63. SIGF(5)=SIG0(5) + XMAT(1)*XCAR(2)*DEPST(5)
  64. SIGF(6)=SIG0(6) + XMAT(1)*XCAR(3)*DEPST(6)
  65. C
  66. C ELASTIC AND PLASTIC AXIAL STIFFNESS
  67. C
  68. EA = XMAT(1)*XCAR(4)
  69. ETA = XMAT(5)*XCAR(4)
  70. C
  71. C
  72. C AXIAL STRAIN INCREMENT
  73. C
  74. DEPSX = DEPST(1)
  75. C
  76. C CALL TO THE SCALAR BILINEAR PLASTIC MODEL COMPUTING MOMY
  77. C
  78. FEP0 = SIG0(1) - VAR0(3)
  79. FE0 = VAR0(3)
  80. CALL BILINE(DEPSX,XMAT(6),EA, ETA, VARF(2),
  81. > FE0,FEP0,SIGF(1),DEPSE)
  82. C
  83. C FINAL INTERNAL VARIABLE
  84. C
  85. VARF(1)=VAR0(1)+DEPSE
  86. VARF(3)=FE0
  87. VARF(4)=VAR0(4)+DEPSX
  88. C
  89. C PLASTIC DEFORMATION INCREMENT
  90. C
  91. DEFP(1) = DEPSE
  92. DEFP(2) = 0.D0
  93. DEFP(3) = 0.D0
  94. DEFP(4) = 0.D0
  95. DEFP(5) = 0.D0
  96. DEFP(6) = 0.D0
  97. RETURN
  98. END
  99.  
  100.  
  101.  
  102.  
  103.  

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